Research Papers

Numerical Solution of Stiff Multibody Dynamic Systems Based on Kinematic Derivatives

[+] Author and Article Information
Paul Milenkovic

Department of Electrical
and Computer Engineering,
University of Wisconsin-Madison,
1415 Engineering Drive,
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 May 25, 2012; final manuscript received June 30, 2014; published online August 8, 2014. Assoc. Editor: YangQuan Chen.

J. Dyn. Sys., Meas., Control 136(6), 061001 (Aug 08, 2014) (9 pages) Paper No: DS-12-1158; doi: 10.1115/1.4028049 History: Received May 25, 2012; Revised June 30, 2014

The Hermite–Obreshkov–Padé (HOP) method of numerical integration is applicable to stiff systems of differential equations, where the linearization has large range of eigenvalues. A practical implementation of HOP requires the ability to determine high-order time derivatives of the system variables. In the case of a constrained multibody dynamical system, the power series solution for the kinematic differential equation is the foundation for an algorithmic differentiation (AD) procedure determining those derivatives. The AD procedure is extended in this paper to determine rates of change in the time derivatives with respect to variation in the position and velocity state variables of the multibody system. The coefficients of this variation form the Jacobian matrix required for Newton–Raphson iteration. That procedure solves the implicit relations for the state variables at the end of each integration time step. The resulting numerical method is applied to the rotation of a dynamically unbalanced constant-velocity (CV) shaft coupling, where the deflection angle of the output shaft is constrained to low levels by springs of high rate and damping.

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


Griewank, A., 1996, “ODE Solving via Automatic Differentiation and Rational Prediction,” Numerical Analysis 1995, Vol. 344, D. F.Griffiths, and G. A.Watson, eds., Pitman Research Notes in Mathematics Series, Longman, Harlow, United Kingdom, 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, L.Vulkov, J.Wasniewski, and P.Yalamov, eds., Lecture Notes in Computer Science, Springer, Berlin, Germany, pp. 114–125.
Barrio, R., 2005, “Performance of the Taylor Series Method for ODEs/DAEs,” Appl. Math. Comput., 163(2), pp. 525–545. [CrossRef]
Milenkovic, P., 2012, “Series Solution for Finite Displacement of Single-Loop Spatial Linkages,” ASME J. Mech. Rob., 4(2), p. 021016. [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, PA.
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, Germany, pp. 1089–1098.
Hairer, E., Nørsett, S. P., and Wanner, G., 1993, Solving Ordinary Differential Equations: Nonstiff Problems, Springer, Berlin, Germany.
Hairer, E., and Wanner, G., 1991, Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems, Springer, Berlin, Germany.
Iserles, A., and Nørsett, S. P., 1991, Order Stars, Chapman & Hall, London, United Kingdom.
Martín-Vaquero, J., 2010, “A 17th-Order Radau IIA Method for Package RADAU. Applications in Mechanical Systems,” Comput. Math. Appl., 59(8), pp. 2464–2472. [CrossRef]
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. Integr. Circuits Syst., 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. Integr. Circuits Syst., 31(1), pp. 101–108. [CrossRef]
Milenkovic, P., 2013, “Projective Constraint Stabilization for a Power Series Forward Dynamics Solver,” ASME J. Dyn. Syst., Meas. Control, 135(3), p. 031004. [CrossRef]
Scholz, H. E., 1998, “The Examination of Nonlinear Stability and Solvability of the Algebraic Equations for the Implicit Taylor Series Method,” Appl. Num. Math., 28(2), pp. 439–458. [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. Scientific Comput., 25(4), pp. 1142–1157. [CrossRef]
Featherstone, R., 2010, “Exploiting Sparsity in Operational-Space Dynamics,” Int. J. Rob. Res., 29(10), pp. 1353–1368. [CrossRef]
Featherstone, R., 2005, “Efficient Factorization of the Joint-Space Inertia Matrix for Branched Kinematic Trees,” Int. J. Rob. Res., 24(6), pp. 487–500. [CrossRef]
R.Featherstone, 2011, personal communication.


Grahic Jump Location
Fig. 1

Clemens CV coupling with torsion springs and dampers on the output-shaft pitch and yaw axes



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