Implementation of linear multistep methods for solving constrained equations of motion

A class of multistep methods for numerical integration of the equations of motion of constrained mechanical systems is considered. They are based on different local parametrizations of the differentiable manifold defined by the kinematic constraints. During integration, the kinematic constraints are preserved at all levels: position, velocity, and acceleration. The numerical implementations described in the paper require only one symmetric matrix factorization at each integration step.

[1]  E. Haug,et al.  Generalized Coordinate Partitioning for Dimension Reduction in Analysis of Constrained Dynamic Systems , 1982 .

[2]  Gene H. Golub,et al.  Matrix computations , 1983 .

[3]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[4]  W. Rheinboldt.  Differential-algebraic systems as differential equations on manifolds , 1984 .

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

[6]  C. W. Gear,et al.  Maintianing solution invariants in the numerical solution of ODEs , 1986 .

[7]  Ronald L. Huston,et al.  Dynamics of Multibody Systems , 1988 .

[8]  Ernst Hairer,et al.  The numerical solution of differential-algebraic systems by Runge-Kutta methods , 1989 .

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

[10]  Ernst Hairer,et al.  Runge-Kutta methods for differential-algebraic equations , 1989 .

[11]  Florian A. Potra,et al.  Differential-Geometric Techniques for Solving Differential Algebraic Equations , 1990 .

[12]  U. Ascher,et al.  Projected implicit Runge-Kutta methods for differential-algebraic equations , 1990 .

[13]  Christian Lubich,et al.  Extrapolation integrators for constrained multibody systems , 1991, IMPACT Comput. Sci. Eng..

[14]  B. Leimkuhler,et al.  Numerical solution of differential-algebraic equations for constrained mechanical motion , 1991 .

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

[16]  J. Yen,et al.  Implicit Numerical Integration of Constrained Equations of Motion Via Generalized Coordinate Partitioning , 1992 .

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

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