Stabilization of Baumgarte's method using the Runge-Kutta approach

The dynamic equations of motion of the constrained multibody mechanical system are mixed differential-algebraic equations (DAE). The DAE systems cannot be solved using numerical integration methods that are commonly used for solving ordinary differential equations. 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 Runge-Kutta method is found.