Application of Runge–Kutta–Rosenbrock Methods to the Analysis of Flexible Multibody Systems

Numerical integration methods are discussed for general equations of motion for multibody systems with flexible parts, which are fairly stiff, time-dependent and non-linear. A family of semi-implicit methods, which belong to the class of Runge–Kutta–Rosenbrock methods, with rather weak non-linear stability properties, are developed. These comprise methods of first, second and third order of accuracy that are A-stable and L-stable and hence introduce numerical damping and the filtering of high frequency components. It is shown, both from theory and examples, that it is generally preferable to use deformation mode coordinates to global nodal coordinates as independent variables in the formulation of the equations of motion. The methods are applied to a series of examples consisting of an elastic pendulum, a beam supported by springs, a four-bar mechanism, and a robotic manipulator with collocated control.

[1]  D. A. Calahan,et al.  A stable, accurate method of numerical integration for nonlinear systems , 1968 .

[2]  One-step methods with adaptive stability functions for the integration of differential equations , 1973 .

[3]  R. Alexander.  Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[4]  A. Prothero,et al.  On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations , 1974 .

[5]  W. A. Thornton,et al.  Optimality Criterion Techniques Applied to Mechanical Design , 1978 .

[6]  C. Lubich.  Integration of stiff mechanical systems by Runge-Kutta methods , 1993 .

[7]  Harold P. Frisch.  Industrial Practices and Examples Space Structures and Robots , 1990 .

[8]  Ernst Hairer,et al.  Error of Runge-Kutta methods for stiff problems studied via differential algebraic equations , 1988 .

[9]  J. P. Meijaard,et al.  A Comparison of Numerical Integration Methods With a View to Fast Simulation of Mechanical Dynamical Systems , 1990 .

[10]  R. Piché.  An L-stable Rosenbrock method for step-by-step time integration in structural dynamics , 1995 .

[11]  van Dh Dick Campen,et al.  Optimization of Multibody Systems Using Approximation Concepts , 1996 .

[12]  J. C. Simo,et al.  Non-linear dynamics of three-dimensional rods: Exact energy and momentum conserving algorithms , 1995 .

[13]  C. F. Haines.  Implicit integration processes with error estimate for the numerical solution of differential equations , 1969, Comput. J..

[14]  T. Steihaug,et al.  An attempt to avoid exact Jacobian and nonlinear equations in the numerical solution of stiff differential equations , 1979 .

[15]  J. Meijaard.  Direct determination of periodic solutions of the dynamical equations of flexible mechanisms and manipulators , 1991 .

[16]  L. Shampine,et al.  Computer solution of ordinary differential equations : the initial value problem , 1975 .

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

[18]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[19]  Adrian Sandu,et al.  Benchmarking stiff ode solvers for atmospheric chemistry problems II: Rosenbrock solvers , 1997 .

[20]  Ben Jonker,et al.  A Finite Element Dynamic Analysis of Flexible Manipulators , 1990, Int. J. Robotics Res..

[21]  J. B. Jonker.  A finite element dynamic analysis of spatial mechanisms with flexible links , 1989 .

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

[23]  J. P. Meijaard,et al.  Validation of flexible beam elements in dynamics programs , 1996 .

[24]  F. Kollmann,et al.  Application of modern time integrators to Hart's inelastic model , 1999 .

[25]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[26]  J. B. Jonker,et al.  A Finite Element Dynamic Analysis of Robotic Manipulators with Flexible Links , 1988 .

[27]  Arend L. Schwab,et al.  Small Vibrations Superimposed on a Prescribed Rigid Body Motion , 2002 .

[28]  J. C. Simo,et al.  Unconditionally stable algorithms for rigid body dynamics that exactly preserve energy and momentum , 1991 .

[29]  D. Arnold,et al.  Computer Solution of Ordinary Differential Equations. , 1981 .

[30]  I. W. Sandberg,et al.  Numerical integration of systems of stiff nonlinear differential equations , 1968 .

[31]  H. H. Rosenbrock,et al.  Some general implicit processes for the numerical solution of differential equations , 1963, Comput. J..

[32]  M. Crisfield,et al.  Energy‐conserving and decaying Algorithms in non‐linear structural dynamics , 1999 .

[33]  J. C. Simo,et al.  Exact energy-momentum conserving algorithms and symplectic schemes for nonlinear dynamics , 1992 .

[34]  Adrian Sandu,et al.  Benchmarking Stiff ODE Solvers for Atmospheric Chemistry Problems I: Implicit versus Explicit , 1996 .

[35]  P. Rentrop,et al.  Generalized Runge-Kutta methods of order four with stepsize control for stiff ordinary differential equations , 1979 .

[36]  J. C. Simo,et al.  How to render second order accurate time-stepping algorithms fourth order accurate while retaining the stability and conservation properties , 1994 .

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

[38]  G. N. Sandor,et al.  A general method of kineto-elastodynamic design of high speed mechanisms , 1973 .

[39]  Edward J. Haug,et al.  Design Sensitivity Analysis and Optimization of Kinematically Driven Systems , 1984 .

[40]  Ben Jonker,et al.  A finite element dynamic analysis of flexible spatial mechanisms and manipulators , 1988 .

[41]  M. Roche,et al.  Rosenbrock methods for Differential Algebraic Equations , 1987 .

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

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

[44]  Lesław K. Bieniasz,et al.  Finite-difference electrochemical kinetic simulations using the Rosenbrock time integration scheme , 1999 .