BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David Martín de Diego

David Martín de Diego appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPmath.DGmath.NAmath.OCNumerical Analysismath.SGmath.DSeess.SYphysics.class-phSystems and Control

Catalog footprint

What is connected

26works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

A parallel iterative method for variational integration

Discrete variational methods show excellent performance in numerical simulations of different mechanical systems. In this paper, we introduce an iterative procedure for the solution of discrete variational equations for boundary value problems. More concretely, we explore a parallelization strategy that leverages the capabilities of multicore CPUs and GPUs (graphics cards). We study this parallel method for higher-order Lagrangian systems, which appear in fully-actuated problems and beyond. The most important part of the paper is devoted to a precise study of different convergence conditions for these methods. We illustrate their excellent behavior in some interesting examples, namely Zermelo's navigation problem, a fuel-optimal navigation problem, interpolation problems or in a fuel optimization problem for a controlled 4-body problem in astrodynamics showing the potential of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Presymplectic integrators for optimal control problems via retraction maps

Retractions maps are used to define a discretization of the tangent bundle of the configuration manifold as two copies of the configuration manifold where the dynamics take place. Such discretization maps can be conveniently lifted to the cotangent bundle so that symplectic integrators are constructed for Hamilton's equations. Optimal control problems are provided with a Hamiltonian framework by Pontryagin's Maximum Principle. That is why we use discretization maps and the integrability algorithm to obtain presymplectic integrators for optimal control problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Retraction maps: a seed of geometric integrators

The classical notion of retraction map used to approximate geodesics is extended and rigorously defined to become a powerful tool to construct geometric integrators and it is called discretization map. Using the geometry of the tangent and cotangent bundles, we are able to tangently and cotangent lift such a map so that these lifts inherit the same properties as the original one and they continue to be discretization maps. In particular, the cotangent lift of a discretization map is a natural symplectomorphism, what plays a key role for constructing geometric integrators and symplectic methods. As a result, a wide range of (higer-order) numerical methods are recovered and canonically constructed by using different discretization maps, as well as some operations with Lagrangian submanifolds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Variational integrators for non-autonomous systems with applications to stabilization of multi-agent formations

Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. Variational integrators are an important class of geometric integrators. The general idea for those variational integrators is to discretize Hamilton's principle rather than the equations of motion in a way that preserves some of the invariants of the original system. In this paper we construct variational integrators with fixed time step for time-dependent Lagrangian systems modelling an important class of autonomous dissipative systems. These integrators are derived via a family of discrete Lagrangian functions each one for a fixed time-step. This allows to recover at each step on the set of discrete sequences the preservation properties of variational integrators for autonomous Lagrangian systems, such as symplecticity or backward error analysis for these systems. We also present a discrete Noether theorem for this class of systems. Applications of the results are shown for the problem of formation stabilization of multi-agent systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Contact geometry for simple thermodynamical systems with friction

Using contact geometry we give a new characterization of a simple but important class of thermodynamical systems which naturally satisfy the first law of thermodynamics (total energy preservation) and the second law (increase of entropy). We completely clarify its qualitative dynamics, the underlying geometrical structures and we show how to use discrete gradient methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Geometric Optimal Trajectory Tracking of Nonholonomic Mechanical Systems

We study the tracking of a trajectory for a nonholonomic system by recasting the problem as a constrained optimal control problem. The cost function is chosen to minimize the error in positions and velocities between the trajectory of a nonholonomic system and the desired reference trajectory, both evolving on the distribution which defines the nonholonomic constraints. The problem is studied from a geometric framework. Optimality conditions are determined by the Pontryagin Maximum Principle and also from a variational point of view, which allows the construction of geometric integrators. Examples and numerical simulations are shown to validate the results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Variational order for forced Lagrangian systems II: Euler-Poincaré equations with forcing

In this paper we provide a variational derivation of the Euler-Poincaré equations for systems subjected to external forces using an adaptation of the techniques introduced by Galley and others. Moreover, we study in detail the underlying geometry which is related to the notion of Poisson groupoid. Finally, we apply the previous construction to the formal derivation of the variational error for numerical integrators of forced Euler-Poincaré equations and the application of this theory to the derivation of geometric integrators for forced systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy-Preserving Integrators Applied to Nonholonomic Systems

We introduce energy-preserving integrators for nonholonomic mechanical systems. We will see that the nonholonomic dynamics is completely determined by a triple $({\mathcal D}^*, Π, \mathcal{H})$, where ${\mathcal D}^*$ is the dual of the vector bundle determined by the nonholonomic constraints, $Π$ is an almost-Poisson bracket (the nonholonomic bracket) and $\mathcal{H}:{\mathcal D}^{*}\rightarrow \mathbb{R}$ is a Hamiltonian function. For this triple, we can apply energy-preserving integrators, in particular, we show that discrete gradients can be used in the numerical integration of nonholonomic dynamics. By construction, we achieve preservation of the constraints and of the energy of the nonholonomic system. Moreover, to facilitate their applicability to complex systems which cannot be easily transformed into the aforementioned almost-Poisson form, we rewrite our integrators using just the initial information of the nonholonomic system. The derived procedures are tested on several examples: A chaotic quartic nonholonomic mechanical system, the Chaplygin sleigh system, the Suslov problem and a continuous gearbox driven by an asymmetric pendulum. Their performace is compared with other standard methods in nonholonomic dynamics, and their merits verified in practice.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hamilton-Jacobi theory, Symmetries and Coisotropic Reduction

Reduction theory has played a major role in the study of Hamiltonian systems. On the other hand, the Hamilton-Jacobi theory is one of the main tools to integrate the dynamics of certain Hamiltonian problems and a topic of research on its own. Moreover, the construction of several symplectic integrators rely on approximations of a complete solution of the Hamilton-Jacobi equation. The natural question that we address in this paper is how these two topics (reduction and Hamilton-Jacobi theory) fit together. We obtain a reduction and reconstruction procedure for the Hamilton-Jacobi equation with symmetries, even in a generalized sense to be clarified below. Several applications and relations to other reductions of the Hamilton-Jacobi theory are shown in the last section of the paper. It is remarkable that as a by-product we obtain a generalization of the Ge-Marsden reduction procedure. Quite surprinsingly, the classical ansatzs available in the literature to solve the Hamilton-Jacobi equation are also particular instances of our framework.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverse problem for Lagrangian systems on Lie algebroids and applications to reduction by symmetries

The language of Lagrangian submanifolds is used to extend a geometric characterization of the inverse problem of the calculus of variations on tangent bundles to regular Lie algebroids. Since not all closed sections are locally exact on Lie algebroids, the Helmholtz conditions on Lie algebroids are necessary but not sufficient, so they give a weaker definition of the inverse problem. As an application the Helmholtz conditions on Atiyah algebroids are obtained so that the relationship between the inverse problem and the reduced inverse problem by symmetries can be described. Some examples and comparison with previous approaches in the literature are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Symplectic groupoids and discrete constrained Lagrangian mechanics

In this article, we generalize the theory of discrete Lagrangian mechanics and variational integrators in two principal directions. First, we show that Lagrangian submanifolds of symplectic groupoids give rise to discrete dynamical systems, and we study the properties of these systems, including their regularity and reversibility, from the perspective of symplectic and Poisson geometry. Next, we use this framework -- along with a generalized notion of generating function due to Sniatycki and Tulczyjew -- to develop a theory of discrete constrained Lagrangian mechanics. This allows for systems with arbitrary constraints, including those which are non-integrable (in an appropriate discrete, variational sense). In addition to characterizing the dynamics of these constrained systems, we also develop a theory of reduction and Noether symmetries, and study the relationship between the dynamics and variational principles. Finally, we apply this theory to discretize several concrete examples of constrained systems in mechanics and optimal control.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalized variational calculus for continuous and discrete mechanical systems

In this paper, we consider a generalization of variational calculus which allows us to consider in the same framework different cases of mechanical systems, for instance, Lagrangian mechanics, Hamiltonian mechanics, systems subjected to constraints, optimal control theory and so on. This generalized variational calculus is based on two main notions: the tangent lift of curves and the notion of complete lift of a vector field. Both concepts are also adapted for the case of skew-symmetric algebroids, therefore, our formalism easily extends to the case of Lie algebroids and nonholonomic systems. Hence, this framework automatically includes reduced mechanical systems subjected or not to constraints. Finally, we show that our formalism can be used to tackle the case of discrete mechanics, including reduced systems, systems subjected to constraints and discrete optimal control theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Geometric integrators for higher-order variational systems and their application to optimal control

Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. In this paper we present a method to construct symplectic-momentum integrators for higher-order Lagrangian systems. Given a regular higher-order Lagrangian $L\colon T^{(k)}Q\to\mathbb{R}$ with $k\geq 1$, the resulting discrete equations define a generally implicit numerical integrator algorithm on $T^{(k-1)}Q\times T^{(k-1)}Q$ that approximates the flow of the higher-order Euler--Lagrange equations for $L$. The algorithm equations are called higher-order discrete Euler--Lagrange equations and constitute a variational integrator for higher-order mechanical systems. The general idea for those variational integrators is to directly discretize Hamilton's principle rather than the equations of motion in a way that preserves the invariants of the original system, notably the symplectic form and, via a discrete version of Noether's theorem, the momentum map. We construct an exact discrete Lagrangian $L_d^e$ using the locally unique solution of the higher-order Euler--Lagrange equations for $L$ with boundary conditions. By taking the discrete Lagrangian as an approximation of $L_d^e$, we obtain variational integrators for higher-order mechanical systems. We apply our techniques to optimal control problems since, given a cost function, the optimal control problem is understood as a second-order variational problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Isotropic submanifolds and the inverse problem for mechanical constrained systems

The inverse problem of the calculus of variations consists in determining if the solutions of a given system of second order differential equations correspond with the solutions of the Euler-Lagrange equations for some regular Lagrangian. This problem in the general version remains unsolved. Here, we contribute to it with a novel description in terms of Lagrangian submanifolds of a symplectic manifold, also valid under some adaptation for the non-autonomous version. One of the advantages of this new point of view is that we can easily extend our description to the study of the inverse problem of the calculus of variations for second order systems along submanifolds. In this case, instead of Lagrangian submanifolds we will use isotropic submanifolds, covering both the nonholonomic and holonomic constraints for autonomous and non-autonomous systems as particular examples. Moreover, we use symplectic techniques to extend these isotropic submanifolds to Lagrangian ones, allowing us to describe the constrained solutions as solutions of a variational problem now without constraints. Mechanical examples such as the rolling disk are provided to illustrate the main results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

New developments on the Geometric Nonholonomic Integrator

In this paper, we will discuss new developments regarding the Geometric Nonholonomic Integrator (GNI) [23, 24]. GNI is a discretization scheme adapted to nonholonomic mechanical systems through a discrete geometric approach. This method was designed to account for some of the special geometric structures associated to a nonholonomic motion, like preservation of energy, preservation of constraints or the nonholonomic momentum equation. First, we study the GNI versions of the symplectic-Euler methods, paying special attention to their convergence behavior. Then, we construct an extension of the GNI in the case of affine constraints. Finally, we generalize the proposed method to nonholonomic reduced systems, an important subclass of examples in nonholonomic dynamics. We illustrate the behavior of the proposed method with the example of the Chaplygin sphere, which accounts for the last two features, namely it is both a reduced and an affine system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variational integrators for underactuated mechanical control systems with symmetries

Optimal control problems for underactuated mechanical systems can be seen as a higher-order variational problem subject to higher-order constraints (that is, when the Lagrangian function and the constraints depend on higher-order derivatives such as the acceleration, jerk or jounces). In this paper we discuss the variational formalism for the class of underactuated mechanical control systems when the configuration space is a trivial principal bundle and the construction of variational integrators for such mechanical control systems. An interesting family of geometric integrators can be defined using discretizations of the Hamilton's principle of critical action. This family of geometric integrators is called variational integrators, being one of their main properties the preservation of geometric features as the symplecticity, momentum preservation and good behavior of the energy. We construct variational integrators for higher-order mechanical systems on trivial principal bundles and their extension for higher-order constrained systems and we devote special attention to the particular case of underactuated mechanical systems

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Higher-order discrete variational problems with constraints

An interesting family of geometric integrators for Lagrangian systems can be defined using discretizations of the Hamilton's principle of critical action. This family of geometric integrators is called variational integrators. In this paper, we derive new variational integrators for higher-order lagrangian mechanical system subjected to higher-order constraints. From the discretization of the variational principles, we show that our methods are automatically symplectic and, in consequence, with a very good energy behavior. Additionally, the symmetries of the discrete Lagrangian imply that momenta is conserved by the integrator. Moreover, we extend our construction to variational integrators where the lagrangian is explicitly time-dependent. Finally, some motivating applications of higher-order problems are considered; in particular, optimal control problems for explicitly time-dependent underactuated systems and an interpolation problem on Riemannian manifolds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The local description of discrete Mechanics

In this paper, we introduce local expressions for discrete Mechanics. To apply our results simultaneously to several interesting cases, we derive these local expressions in the framework of Lie groupoids, following the program proposed by Alan Weinstein in [19]. To do this, we will need some results on the geometry of Lie groupoids, as, for instance, the construction of symmetric neighborhoods or the existence of local bisections. These local descriptions will be particular useful for the explicit construction of geometric integrators for mechanical systems (reduced or not), in particular, discrete Euler-Lagrange equations, discrete Euler-Poincaré equations, discrete Lagrange-Poincaré equations... The results contained in this paper can be considered as a local version of the study that we have started in [13], on the geometry of discrete Mechanics on Lie groupoids.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Morse families in optimal control problems

We geometrically describe optimal control problems in terms of Morse families in the Hamiltonian framework. These geometric structures allow us to recover the classical first order necessary conditions for optimality and the starting point to run an integrability algorithm. Moreover the integrability algorithm is adapted to optimal control problems in such a way that the trajectories originated by discontinuous controls are also obtained. From the Hamiltonian viewpoint we obtain the equations of motion for optimal control problems in the Lagrangian formalism by means of a proper Lagrangian submanifold. Singular optimal control problems and overdetermined ones are also studied along the paper.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinematic reduction and the Hamilton-Jacobi equation

A close relationship between the classical Hamilton-Jacobi theory and the kinematic reduction of control systems by decoupling vector fields is shown in this paper. The geometric interpretation of this relationship relies on new mathematical techniques for mechanics defined on a skew-symmetric algebroid. This geometric structure allows us to describe in a simplified way the mechanics of nonholonomic systems with both control and external forces.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Unified formalism for non-autonomous mechanical systems

We present a unified geometric framework for describing both the Lagrangian and Hamiltonian formalisms of regular and non-regular time-dependent mechanical systems, which is based on the approach of Skinner and Rusk (1983). The dynamical equations of motion and their compatibility and consistency are carefully studied, making clear that all the characteristics of the Lagrangian and the Hamiltonian formalisms are recovered in this formulation. As an example, it is studied a semidiscretization of the nonlinear wave equation proving the applicability of the proposed formalism.

0 citations0 reviewsNo code linkedOpen access