Abstract: The independent coordinates are defined applying singular value decomposition of the constraint Jacobian matrix. Reduction of the dynamic equations of multibody systems are thus established. The process of numerical dynamic analysis is also given.