Abstract: A new method is developed to correct the violations in constraint multibody systems. Compared with any conventional methods, this method has clear physical meaning, less calculation and obvious correction effect. Besides, this method has little influence on the quality of the original dynamic equations. Using this method, the violations of constraint equations can be limited efficiently in the given accuracy. Suggest that the integration is at the step t=t k , and the generalized coordinates k do not...