An automatic constraint violation stabilization method for differential/ algebraic equations of motion in multibody system dynamics

A new automatic constraint violation stabilization method for numerical integration of Euler-Lagrange equations of motion in dynamics of multibody systems is presented. The parameters a, β used in the traditional constraint violation stabilization method are determined according to the integration time time step size and Taylor expansion method automatically. The direct integration method, the traditional constraint violation stabilization method and the new method presented in this paper are compared finally. Communicated by Ma Xingrui Foundation item: the National Natural Science Foundation of China (19902006); the Natural Science Foundation of Sandong Province (Y97F06152), China Biography: Zhao Weijia (1955∼)