Research Papers

Multi-Integral Method for Solving the Forward Dynamics of Stiff Multibody Systems

[+] Author and Article Information
Paul Milenkovic

Department of Electrical and Computer Engineering,
University of Wisconsin-Madison,
Madison, WI 53706
e-mail: phmilenk@wisc.edu

Contributed by the Dynamic Systems Division of ASME for publication in the JOURNAL OF DYNAMIC SYSTEMS, MEASUREMENT, AND CONTROL. Manuscript received August 23, 2012; final manuscript received April 19, 2013; published online June 27, 2013. Assoc. Editor: Sergey Nersesov.

J. Dyn. Sys., Meas., Control 135(5), 051014 (Jun 27, 2013) (9 pages) Paper No: DS-12-1274; doi: 10.1115/1.4024300 History: Received August 23, 2012; Revised April 19, 2013

The Hermite–Obreshkov–Padé (HOP) procedure is an implicit method for the numerical solution of a system of ordinary differential equations (ODEs) applicable to stiff dynamical systems. This procedure applies an Obreshkov condition to multiple derivatives of the system state vector, both at the start and end of a time step in the numerical solution. That condition is shown to be satisfied by the Hermite interpolating polynomial that matches the state vector and its derivatives, also at the start and end of a time step. The Hermite polynomial, in turn, can be specified in terms of the system state and its derivatives at the start of a step together with a collection of free parameters. Adjusting these free parameters to minimize magnitudes of the ODE residual and its derivatives at the end of a step serves as a proxy for matching the system state and its derivatives. A high-order Taylor expansion at the start of a time step interval models the residual and its derivatives over the entire interval. A variant of this procedure adjusts those parameters to match integrals of the system state over the duration of that interval. This is done by minimizing magnitudes of integrals of the ODE residual calculated from the extrapolating Taylor-series expansion, a process that avoids the need to determine integration constants for multiple integrals of the state. This alternative method eliminates the calculation of high-order derivatives of the system state and hence avoids loss in accuracy from floating-point round off. Numerical performance is evaluated on a dynamically unbalanced constant-velocity (CV) coupling having a high spring rate constraining shaft deflection.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Hairer, E., and Wanner, G., 1991, Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems, Springer, Berlin.
Gear, C., 1971, “Simultaneous Numerical Solution of Differential-Algebraic Equations,” IEEE Trans. Circuit Theory, 18(1), pp. 89–95. [CrossRef]
Negrut, D., Haug, E. J., and German, H. C., 2003, “An Implicit Runge–Kutta Method for Integration of Differential Algebraic Equations of Multibody Dynamics,” Multibody Syst. Dyn., 9(2), pp. 121–142. [CrossRef]
Eriksson, K., Johnson, C., and Logg, A., 2004, “Explicit Time-Stepping for Stiff ODEs,” SIAM J. Sci. Comput. (USA), 25(4), pp. 1142–1157. [CrossRef]
Gear, C. W., and Kevrekidis, I. G., 2003, “Projective Methods for Stiff Differential Equations: Problems With Gaps in Their Eigenvalue Spectrum,” SIAM J. Sci. Comput. (USA), 24(4), pp. 1091–1106. [CrossRef]
Milenkovic, P., “Numerical Solution of Stiff Multi-Body Dynamic Systems Based on Kinematic Derivatives,” ASME J. Dyn. Syst., Meas., Control (submitted).
Griewank, A., 1996, “ODE Solving via Automatic Differentiation and Rational Prediction,” Numerical Analysis 1995, Vol. 344 (Pitman Research Notes in Mathematics Series), D. F.Griffiths and G. A.Watson, eds., Longman, Harlow, UK, pp. 36–56.
Corliss, G., Griewank, A., and Henneberger, P., 1997, “High-Order Stiff ODE Solvers via Automatic Differentiation and Rational Prediction,” Numerical Analysis and Its Applications, Vol. 1196 (Lecture Notes in Computer Science), L.Vulkov, J.Wasniewski, and P.Yalamov, eds., Springer, Berlin, pp. 114–125.
Gad, E., Nakhla, M., and Achar, R., 2009, “A-Stable and L-Stable High-Order Integration Methods for Solving Stiff Differential Equations,” IEEE Trans. Comput.-Aided Des., 28(9), pp. 1359–1372. [CrossRef]
Zhou, Y., Gad, E., and Nakhla, M. S., 2012, “Structural Characterization and Efficient Implementation Techniques For-Stable High-Order Integration Methods,” IEEE Trans. Comput.-Aided Des., 31(1), pp. 101–108. [CrossRef]
Milenkovic, P., 2011, “Solution of the Forward Dynamics of a Single-Loop Linkage Using Power Series,” ASME J. Dyn. Syst., Meas., Control, 133(6), p. 061002. [CrossRef]
Griewank, A., and Walther, A., 2008, Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation, Society for Industrial and Applied Mathematics (SIAM), Philadelphia.
Griewank, A., and Walther, A., 2006, “On the Efficient Generation of Taylor Expansions for DAE Solutions by Automatic Differentiation,” Applied Parallel Computing, State of the Art in Scientific Computing, J.Dongarra, K.Madsen, and J.Wasniewski, eds., Springer, Berlin, pp. 1089–1098.
Milenkovic, P., 2012, “Series Solution for Finite Displacement of Single-Loop Spatial Linkages,” ASME J. Mech. Rob., 4(2), p. 021016. [CrossRef]
Scholz, H. E., 1998, “The Examination of Nonlinear Stability and Solvability of the Algebraic Equations for the Implicit Taylor Series Method,” Appl. Numer. Math., 28(2), pp. 439–458. [CrossRef]
Wall, H. S., 1948, Analytic Theory of Continued Fractions, D. Van Nostrand, New York.
Barrio, R., 2005, “Performance of the Taylor Series Method for ODEs/DAEs,” Appl. Math. Comput., 163(2), pp. 525–545. [CrossRef]
Wilf, H. S., 1994, Generatingfunctionology, Academic Press, San Diego, CA.
Milenkovic, P., 2013, “Projective Constraint Stabilization for a Power Series Forward Dynamics Solver,” ASME J. Dyn. Syst., Meas., Control, 135(3), p. 031004. [CrossRef]





Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In