Symplectic integration of constrained Hamiltonian systems

A Hamiltonian system in potential form (H(q, p) = ptMp/2 + F(q)) subject to smooth constraints on q can be viewed as a Hamiltonian system on a manifold, but numerical computations must be performed in Rn . In this paper, methods which reduce "Hamiltonian differential-algebraic equations" to ODEs in Euclidean space are examined. The authors study the construction of canonical parametrizations or local charts as well as methods based on the construction of ODE systems in the space in which the constraint manifold is embedded which preserve the constraint manifold as an invariant manifold. In each case, a Hamiltonian system of ordinary differential equations is produced. The stability of the constraint-invariants and the behavior of the original Hamiltonian along solutions are investigated both numerically and analytically.

[1]  R. Courant,et al.  Methods of Mathematical Physics , 1962 .

[2]  T. E. Hull,et al.  Numerical solution of initial value problems , 1966 .

[3]  Ralph Abraham,et al.  Foundations Of Mechanics , 2019 .

[4]  J. Baumgarte.  Stabilization of constraints and integrals of motion in dynamical systems , 1972 .

[5]  V. Arnold.  Mathematical Methods of Classical Mechanics , 1974 .

[6]  P. Olver.  Applications of Lie Groups to Differential Equations , 1986 .

[7]  W. Rheinboldt.  Numerical analysis of parametrized nonlinear equations , 1986 .

[8]  Anthony M. Bloch,et al.  Control of constrained Hamiltonian systems and applications to control of constrained robots , 1986 .

[9]  G. J. Cooper.  Stability of Runge-Kutta Methods for Trajectory Problems , 1987 .

[10]  J. M. Sanz-Serna,et al.  Runge-kutta schemes for Hamiltonian systems , 1988 .

[11]  B. Leimkuhler,et al.  Stabilization and projection methods for multibody dynamics , 1990 .

[12]  F. Potra,et al.  On the Numerical Solution of Euler-Lagrange Equations∗ , 1991 .

[13]  J. M. Sanz-Serna,et al.  Symplectic integrators for Hamiltonian problems: an overview , 1992, Acta Numerica.

[14]  Reich Sebastian,et al.  Symplectic Integration of Constrained Hamiltonian Systems by Runge-Kutta Methods , 1993 .

[15]  Uri M. Ascher,et al.  Stability of Computational Methods for Constrained Dynamics Systems , 1991, SIAM J. Sci. Comput..

[16]  Benedict J. Leimkuhler,et al.  Discretization And Weak Invariants , 1994 .

[17]  W. Rheinboldt,et al.  On the numerical solution of the Euler-Lagrange equations , 1995 .

[18]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[19]  L. Jay.  Symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems , 1996 .