Stability and accuracy analysis of Baumgarte's constraint violation stabilization method

When Baumgarte's Constraint Violation Stabilization Method (CVSM) is used in the simulation of Lagrange equations of motion with holonomic constraints, it is shown that, with suitable assumptions on the integration step size h and the eigen-values (lambda's) of the linearized system, the constraint variable are effectively integrated by the same algorithm as that used for the state variables. A numerical stability analysis of the constraint violations can be performed using this so-called pseudo-integration equation. A study is also made of truncation errors and their modeling in the continuous time domain. This model can be used to determine the effectiveness of various constraint controls and integration methods in reducing the errors in the solution due to truncation errors. Examples are presented to illustrate the use of a higher-order truncation error model which leads to an accurate quantitative steady-state analysis of th constraint violations.