An Implicit Runge–Kutta Method for Integration of Differential Algebraic Equations of Multibody Dynamics

When performing dynamic analysis of a constrained mechanical system, aset of index 3 Differential Algebraic Equations (DAE) describes the timeevolution of the model. This paper presents a state space DAE solutionframework that can embed an arbitrary implicit Ordinary DifferentialEquations (ODE) code for numerical integration of a reduced set of statespace ordinary differential equations. This solution framework isconstructed with the goal of leveraging with minimal effort establishedoff the shelf implicit ODE integrators for efficiently solving the DAEof multibody dynamics. This concept is demonstrated by embedding awell-known public domain singly diagonal implicit Runge–Kutta code inthe framework provided. The resulting L-stable, stiffly accurateimplicit algorithm is shown to be two orders of magnitude faster than astate of the art explicit algorithm when used to simulate a stiffvehicle model.

[1]  L. Shampine.  Conservation laws and the numerical solution of ODEs, II , 1999 .

[2]  U. Ascher,et al.  Stabilization of Constrained Mechanical Systems with DAEs and Invariant Manifolds , 1995 .

[3]  Edward J. Haug,et al.  Dynamic and sensitivity analysis of multibody systems , 1998 .

[4]  E. Eich.  Convergence results for a coordinate projection method applied to mechanical systems with algebraic constraints , 1993 .

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

[6]  N. K. Mani,et al.  Application of Singular Value Decomposition for Analysis of Mechanical System Dynamics , 1985 .

[7]  F. Potra,et al.  ON THE IMPLICIT INTEGRATION OF DIFFERENTIAL-ALGEBRAIC EQUATIONS OF MULTIBODY DYNAMICS by , 1998 .

[8]  Ch . Engstler,et al.  MEXX - Numerical Software for the Integration of Constrained Mechanical Multibody Systems , 1992 .

[9]  Lawrence F. Shampine,et al.  Art of writing a Runge--Kutta code , 1976 .

[10]  Arnold Neumaier,et al.  Introduction to Numerical Analysis , 2001 .

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

[12]  F. Potra,et al.  A Topology Based Approach for Exploiting Sparsity in Multibody Dynamics , 1995 .

[13]  F. Potra.  Implementation of linear multistep methods for solving constrained equations of motion , 1993 .

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

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

[16]  James Demmel,et al.  LAPACK Users' Guide, Third Edition , 1999, Software, Environments and Tools.

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

[18]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[19]  Roderic C. Deyo,et al.  Real-Time Integration Methods for Mechanical System Simulation , 1991 .

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

[21]  Lawrence F. Shampine,et al.  The art of writing a Runge-Kutta code. II. , 1979 .

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

[23]  U. Ascher,et al.  Stabilization of DAEs and invariant manifolds , 1994 .

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

[25]  Begnaud Francis Hildebrand,et al.  Introduction to numerical analysis: 2nd edition , 1987 .

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

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

[28]  Ö. Ólafsson,et al.  Modeling and velocity stabilization of constrained mechanical systems , 1994 .

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

[30]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[31]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[32]  George M. Lance,et al.  A Differentiable Null Space Method for Constrained Dynamic Analysis , 1987 .