The object of this study is to solve the stability problem for the numerical integration of constrained multibody mechanical systems. The dynamic equations of motion of the constrained multibody mechanical system are mixed differential-algebraic equations (DAE). In applying numerical integration methods to this equation, constrained equations and their first and second derivatives must be satisfied simultaneously. That is, the generalized coordinates and their derivatives are dependent. Direct integration methods do not consider this dependency and constraint violation occurs. To solve this problem, Baumgarte proposed a constraint stabilization method in which a position and velocity terms were added in the second derivative of the constraint equation. The disadvantage of this method is that there is no reliable method for selecting the coefficients of the position and velocity terms. Improper selection of these coefficients can lead to erroneous results. In this study, stability analysis methods in digital control theory are used to solve this problem. Correct choice of the coefficients for the Adams method are found for both fixed and variable integration step size.

1.
Haug, E. J., 1989, “Computer Aided Kinematics and Dynamics of Mechanical System,” Vol. I: Basic Methods, Allyn and Bacon, Massachusetts, USA.
2.
DADS User’ s Manual, 1987, Rev. 4.0, Computer Aided Design Software, Inc. Oakdale, Iowa.
3.
ADAMS User’ s Guide, 1981, Mechanical Dynamics, Inc., 555 South Fortest, Ann Arbor, Michigan.
4.
Bae, D. S., Hwang, R. S., and Haug, E. J., 1988, “A Recursive Formulation for Real-Time Dynamic Simulation,” Advances in Design Automation, ASME, New York, pp. 499–508.
5.
Hwang, R. S., Bae, D. S., Haug, E. J., and Kuhl, J. G., 1988, “Parallel Processing for Real-Time Dynamic System Simulation,” Advances in Design Automation, ASME, New York, pp. 509–518.
6.
Chung
S.
, and
Haug
E. J.
,
1993
, “
Real-Time Simulation of Multibody Dynamics on Shared Memory Multiprocessors
,”
ASME, JOURNAL OF DYNAMIC SYSTEMS, MEASUREMENT, AND CONTROL
, Vol.
115
, pp.
627
637
.
7.
Petzold
L.
,
1982
, “
Differential/Algebraic Equations are not ODEs
,”
SIAM J. Sci. Stat. Comput.
, Vol.
3
, No.
3
, pp.
367
384
.
8.
Wehage
R. A.
, and
Haug
E. J.
,
1982
, “
Generalized Coordinates Partitioning for Dimension Reduction in Analysis of Constrained Dynamic System
,”
ASME JOURNAL OF MECHANICAL DESIGN
, Vol.
104
, pp.
247
255
.
9.
Baumgarte
J.
,
1972
, “
Stabilization of Constraints and Integrals of Motion in Dynamical Systems
,”
Computer Methods in Applied Mechanics and Engineering
, Vol.
1
, pp.
1
16
.
10.
Chang
C. O.
, and
Nikravesh
P. E.
,
1985
, “
An Adaptive Constraint Violation Stabilization Method for Dynamic Analysis of Mechanical Systems
,”
ASME JOURNAL OF MECHANISMS, TRANSMISSIONS AND AUTOMATION IN DESIGN
, Vol.
17
, pp.
488
492
.
11.
Park
K. C.
, and
Chiou
J. C.
,
1988
, “
Stabilization of Computational Procedures for Constrained Dynamical Systems
,”
AIAA Journal of Guidance and Control
, Vol.
11
, pp.
365
370
.
12.
Bauchau
O. A.
,
Damilano
G.
, and
Theron
N. J.
,
1995
, “
Numerical Integration of Non-Linear Elastic Multi-Body Systems
,”
International Journal for Numerical Methods in Engineering
, Vol.
38
, pp.
2727
2751
.
13.
Rosen
A.
, and
Edelstein
E.
,
1997
, “
Investigation of a New Formulation of the Lagrange Method for Constrained Dynamic Systems
,”
ASME JOURNAL OF APPLIED MECHANICS
, Vol.
64
, pp.
116
122
.
14.
Yoon
S.
,
Howe
R. M.
, and
Greenwood
D. T.
,
1994
, “
Geometric Elimination of Constraint Violations in Numerical Simulation of Lagrange Equations
,”
ASME JOURNAL OF MECHANICAL DESIGN
, Vol.
116
, pp.
1058
1064
.
15.
Yoon
S.
,
Howe
R. M.
, and
Greenwood
D. T.
,
1995
, “
Stability and Accuracy Analysis of Baumgarte’s Constraint Stabilization Method
,”
ASME JOURNAL OF MECHANICAL DESIGN
, Vol.
117
, pp.
446
453
.
16.
Park
T.
,
1986
, “
A Hybrid Constraint Stabilization-Generalized Coordinate Partitioning Method for Machine Dynamics
,”
ASME JOURNAL OF MECHANISMS, TRANSMISSIONS, AND AUTOMATION IN DESIGN
, Vol.
108
, No.
2
, pp.
211
216
.
17.
Haug
E. J.
, and
Yen
J.
,
1992
, “
Implicit Numerical Integration of Constrained Equations of motion Via Generalized Coordinate Partitioning
,”
ASME JOURNAL OF MECHANICAL DESIGN
, Vol.
114
, pp.
296
304
.
18.
Ogata, K., 1995, Discrete-Time Control Systems, 2nd Edition, Prentice-Hall, Englewood Cliffs, New Jersey, USA.
19.
Franklin, G. F., Powell, A. E.-N., 1994, Feedback Control of Dynamic Systems, 3rd Edition, Addison-Wesley.
20.
Lambert, J. D., 1991, Numerical Methods for Ordinary Differential Systems, John Wiley & Sons, West Sussex, England.
This content is only available via PDF.
You do not currently have access to this content.