Constraint violation stabilization using gradient feedback in constrained dynamics simulation

Conventional holonomic or nonholonomic constraints are defined as geometric constraints in this paper. If the total energy of a dynamic system can be computed from the initial energy plus the time integral of the energy input rate due to external or internal forces, then the total energy can be artificially treated as a constraint. The violation of the total energy constraint due to numerical errors during simulation can be used as information to control these errors. When geometric constraint control is combined with energy constraint control, numerical simulation of a constrained dynamic system becomes more accurate. An energy constraint control based on the gradient feedback of the energy constraint violation leads to a new method to control both geometric and energy constraint violations, so-called constraint violation stabilization using gradient feedback. A new convenient and effective method to implement energy constraint control in numerical simulation is developed based on the geometric interpretation of the relation between constraints in the phase space. Several combinations of energy constraint control with either Baumgarte's constraint violation stabilization method or the new constraint violation stabilization using gradient feedback are also addressed. Finally, a new method for implementing constraint controls is developed by using the Euler method for integrating constraint control terms, even when higher-order integration methods are used for all other integrations.

[1]  Donald A. Smith,et al.  DAMN - Digital Computer Program for the Dynamic Analysis of Generalized Mechanical Systems , 1971 .

[2]  L. Petzold.  Differential/Algebraic Equations are not ODE's , 1982 .

[3]  J. Baumgarte.  A New Method of Stabilization for Holonomic Constraints , 1983 .

[4]  Parviz E. Nikravesh,et al.  Some Methods for Dynamic Analysis of Constrained Mechanical Systems: A Survey , 1984 .

[5]  K. Park,et al.  Stabilization of computational procedures for constrained dynamical systems , 1988 .

[6]  D. T. Greenwood.  Principles of dynamics , 1965 .

[7]  R. Rohrer,et al.  Automated Design of Biasing Circuits , 1971 .

[8]  Thomas R. Kane,et al.  Formulation of Equations of Motion for Complex Spacecraft , 1980 .

[9]  J. Baumgarte.  Stabilization of constraints and integrals of motion in dynamical systems , 1972 .

[10]  C. W. Gear,et al.  Simultaneous Numerical Solution of Differential-Algebraic Equations , 1971 .

[11]  C. W. Gear,et al.  Differential-Algebraic Equations , 1984 .

[12]  C. W. Gear,et al.  Automatic integration of Euler-Lagrange equations with constraints , 1985 .

[13]  R. M. Howe,et al.  Space trajectory computation at the University of Michigan , 1966 .

[14]  William W. Hooker,et al.  The Dynamical Attitude Equations for n-Body Satellite , 1965 .

[15]  Sugjoon Yoon.  Real-time simulation of constrained dynamic systems. , 1990 .

[16]  A. Hausner.  Analog and analog/hybrid computer programming , 1971 .

[17]  R. Roberson.  A form of the translational dynamical equations for relative motion in systems of many non-rigid bodies , 1972 .