0
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.

FIGURES IN THIS ARTICLE
<>
Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Fig. 1

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

Tables

Errata

Discussions

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