RATTLE Method for Dissipative Constrained Hamiltonian Systems

Firstly, the Hamiltonian formulation for the constrained dissipative systems was deduced, it is index-3 DAEs. The index-2 DAEs was obtained based on the GGL stabilized method. Then RATTLE method that proposed for motion equations of conservative constrained Hamiltonian systems was extended to solve dissipative constrained Hamiltonian systems. For the dissipative systems, the symplectic structure of RATTLE method is no longer preserved, but this method can capture the decay of the energy accurately because of no numerical dissipation. Numerical experiment results illustrate the effectiveness of the method.

[1]  L. Petzold Differential/Algebraic Equations are not ODE's , 1982 .

[2]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[3]  Werner C. Rheinboldt,et al.  Nonholonomic motion of rigid mechanical systems from a DAE viewpoint , 1987 .

[4]  採編典藏組 Society for Industrial and Applied Mathematics(SIAM) , 2008 .

[5]  Laurent O. Jay,et al.  Specialized Partitioned Additive Runge-Kutta Methods for Systems of Overdetermined DAEs with Holonomic Constraints , 2007, SIAM J. Numer. Anal..

[6]  Florian A. Potra,et al.  ODAE methods for the numerical solution of Euler-Lagrange equations , 1992 .

[7]  O. Bauchau,et al.  Review of Classical Approaches for Constraint Enforcement in Multibody Systems , 2008 .

[8]  J. Marsden,et al.  Discrete mechanics and variational integrators , 2001, Acta Numerica.

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

[10]  E. J. Haug,et al.  Computer aided kinematics and dynamics of mechanical systems. Vol. 1: basic methods , 1989 .

[11]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[12]  H. C. Andersen Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[13]  C. W. Gear,et al.  Automatic integration of Euler-Lagrange equations with constraints , 1985 .

[14]  Javier García de Jalón,et al.  Kinematic and Dynamic Simulation of Multibody Systems: The Real Time Challenge , 1994 .

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

[16]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[17]  Ernst Hairer,et al.  Simulating Hamiltonian dynamics , 2006, Math. Comput..

[18]  B. Leimkuhler,et al.  Symplectic Numerical Integrators in Constrained Hamiltonian Systems , 1994 .

[19]  Javier García de Jalón,et al.  Kinematic and Dynamic Simulation of Multibody Systems , 1994 .

[20]  O. Bauchau,et al.  Review of Contemporary Approaches for Constraint Enforcement in Multibody Systems , 2008 .