A DISCUSSION OF LOW ORDER NUMERICAL INTEGRATION FORMULAS FOR RIGID AND FLEXIBLE MULTIBODY DYNAMICS

The premise of this work is that real-life mechanical systems limit the use of high order integration formulas due to the presence in the associated models of friction and contact/impact elements. In such cases producing a numerical solution necessarily relies on low order integration formulas. The resulting algorithms are generally robust and expeditious; their major drawback remains that they typically require small integration step-sizes in order to meet a user prescribed accuracy. This paper looks at three low order numerical integration formulas: Newmark, HHT, and BDF of order two. These formulas are used in two contexts. A first set of three methods is obtained by considering a direct index-3 discretization approach that solves for the equations of motion and imposes the position kinematic constraints. The second batch of three additional methods draws on the HHT and BDF integration formulas and considers in addition to the equations of motion both the position and velocity kinematic constraint equations. The first objective of this paper is to review the theoretical results available in the literature regarding the stability and convergence properties of these low order methods when applied in the context of multibody dynamics simulation. When no theoretical results are available, numerical experiments are carried out to gauge order behavior. The second objective is to perform a set of numerical experiments to compare these six methods in terms of several metrics: (a) efficiency, (b) velocity constraint drift, and (c) energy preservation. A set of simple mechanical systems is used for this purpose: a double pendulum, a slider crank with rigid bodies, and a slider crank with a flexible body represented in the floating frame of reference formulation.Copyright © 2007 by ASME

[1]  R. Y. Yakoub,et al.  Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Theory , 2001 .

[2]  N. Orlandea,et al.  A Sparsity-Oriented Approach to the Dynamic Analysis and Design of Mechanical Systems—Part 2 , 1977 .

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

[4]  J. Yen.  Constrained equations of motion in multibody dynamics as ODEs on manifolds , 1993 .

[5]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

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

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

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

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

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

[11]  C. W. Gear,et al.  Simultaneous Numerical Solution of Differential-Algebraic Equations , 1971 .

[12]  Dan Negrut,et al.  A MATLAB implementation of the seven-body mechanism for implicit integration of the constrained equations of motion , 2007 .

[13]  Dan Negrut,et al.  On the Use of the HHT Method in the Context of Index 3 Differential Algebraic Equations of Multibody Dynamics , 2005 .

[14]  Olivier A. Bauchau,et al.  Robust integration schemes for flexible multibody systems , 2003 .

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

[16]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[17]  Dan Negrut,et al.  An Implicit Runge–Kutta Method for Integration of Differential Algebraic Equations of Multibody Dynamics , 2003 .

[18]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[19]  M. Géradin,et al.  Mechanical Vibrations: Theory and Application to Structural Dynamics , 1994 .

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

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

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

[23]  R. März.  Differential algebraic systems anew , 2002 .

[24]  Dan Negrut,et al.  Implicit Integration in Molecular Dynamics Simulation , 2008 .

[25]  M. A. Chace,et al.  A Sparsity-Oriented Approach to the Dynamic Analysis and Design of Mechanical Systems—Part 1 , 1977 .

[26]  M. Iqbal,et al.  On photon correlation measurements of colloidal size distributions using Bayesian strategies , 2000 .

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

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

[29]  E. Haug,et al.  A State-Space-Based Implicit Integration Algorithm for Differential-Algebraic Equations of Multibody Dynamics* , 1997 .

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

[31]  Ahmed A. Shabana,et al.  Dynamics of Multibody Systems , 2020 .

[32]  Björn Engquist,et al.  Backward differentiation approximations of nonlinear differential/algebraic systems , 1988 .

[33]  Per Lötstedt,et al.  Numerical solution of nonlinear differential equations with algebraic constraints I: Convergence results for backward differentiation formulas , 1986 .

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

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

[36]  U. Nowak,et al.  Numerical Integration of Constrained Mechanical Systems Using MEXX , 1995 .

[37]  M. Arnold,et al.  Convergence of the generalized-α scheme for constrained mechanical systems , 2007 .

[38]  Christoph Lunk,et al.  Solving constrained mechanical systems by the family of Newmark and α‐methods , 2006 .

[39]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

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

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

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

[43]  Olivier A. Bauchau,et al.  Time-Step-Size-Independent Conditioning and Sensitivity to Perturbations in the Numerical Solution of Index Three Differential Algebraic Equations , 2007, SIAM J. Sci. Comput..

[44]  Linda,et al.  A Time Integration Algorithm forFlexible Mechanism Dynamics : TheDAE-Method , 1996 .

[45]  Dan Negrut,et al.  A Second Order Extension of the Generalized–α Method for Constrained Systems in Mechanics , 2009 .

[46]  Dan Negrut,et al.  On an Implementation of the Hilber-Hughes-Taylor Method in the Context of Index 3 Differential-Algebraic Equations of Multibody Dynamics (DETC2005-85096) , 2007 .