MATH-PHApr 16, 2011
On the geometry of higher-order variational problems on Lie groupsLeonardo Colombo, David Martin de Diego
In this paper, we describe a geometric setting for higher-order lagrangian problems on Lie groups. Using left-trivialization of the higher-order tangent bundle of a Lie group and an adaptation of the classical Skinner-Rusk formalism, we deduce an intrinsic framework for this type of dynamical systems. Interesting applications as, for instance, a geometric derivation of the higher-order Euler-Poincaré equations, optimal control of underactuated control systems whose configuration space is a Lie group are shown, among others, along the paper.
MATH-PHSep 22, 2011
Discrete second-order Euler-Poincaré equations. Applications to optimal controlLeonardo Colombo, Fernando Jimenez, David Martin de Diego
In this paper we will discuss some new developments in the design of numerical methods for optimal control problems of Lagrangian systems on Lie groups. We will construct these geometric integrators using discrete variational calculus on Lie groups, deriving a discrete version of the second-order Euler-Lagrange equations. Interesting applications as, for instance, a discrete derivation of the Euler-Poincaré equations for second-order Lagrangians and its application to optimal control of a rigid body, and of a Cosserat rod are shown at the end of the paper.
NAMar 4, 2018
Lie-Poisson integratorsDavid Martin de Diego
In this paper, we discuss the geometric integration of hamiltonian systems on Poisson manifolds, in particular, in the case, when the Poisson structure is induced by a Lie algebra, that is, it is a Lie-Poisson structure. A Hamiltonian system on a Poisson manifold $(P, Π)$ is a smooth manifold $P$ equipped with a bivector field $Π$ satisfying $[Π, Π]=0$ (Jacobi identity), inducing the Poisson bracket on $C^{\infty}(P)$, $\{f, g\}\equiv Π(df, dg)$ where $f, g\in C^{\infty}(P)$. For any $f\in C^{\infty}(P)$ the Hamiltonian vector field is defined by $X_f(g)=\{g, f\}$. The Hamiltonian vector fields $X_f$ generate an integrable generalized distribution on $P$ and the leaves of this foliation are symplectic. The flow of any hamiltonian vector field preserves the Poisson structure, it fixes each leaf and the hamiltonian itself is a first integral. It is important to characterize numerical methods preserving some of these fundamental properties of the hamiltonian flow on Poisson manifolds (geometric integrators). We discuss the difficulties of deriving these Poisson methods using standard techniques and we present some modern approaches to the problem.
MATH-PHJan 1, 2012
Hamiltonian dynamics and constrained variational calculus: continuous and discrete settingsManuel de Leon, Fernando Jimenez, David Martin de Diego
The aim of this paper is to study the relationship between Hamiltonian dynamics and constrained variational calculus. We describe both using the notion of Lagrangian submanifolds of convenient symplectic manifolds and using the so-called Tulczyjew's triples. The results are also extended to the case of discrete dynamics and nonholonomic mechanics. Interesting applications to geometrical integration of Hamiltonian systems are obtained.
47.0NAApr 3
Geometric Integrators for Nonholonomic Systems on Lie GroupsViyom Vivek, David Martin de Diego, Ravi N. Banavar
We present a general framework for constructing structure-preserving numerical integrators for nonholonomically constrained mechanical systems evolving on Lie groups using retraction maps. Retraction maps generalize the exponential map and provide a convenient tool for performing numerical integration on manifolds. In nonholonomic mechanics, the constraints restrict the dynamics to a nonintegrable distribution rather than the entire tangent bundle. Using the Hamel formulation, the equations of motion can be expressed in local coordinates adapted to this constraint distribution. We then specialize the framework to the case of Lie groups, where both the dynamics and the constraints exhibit symmetries, allowing a simplified formulation of the numerical scheme. The resulting integrator respects the constraint distribution and enforces the nonholonomic constraints at each discrete time step. The approach is illustrated using the Suslov problem.