Numerical Solution for Stiff Dynamic Equations of Flexible Multibody System

A nonlinear numerical integration method,based on forward and backward Euler integration methods,is proposed for solving the stiff dynamic equations of a flexible multibody system,which are transformed from the second order to the first order by adopting state variables.This method is of A0 stability and infinity stability.The numerical solutions violating the constraint equations are corrected by Blajer’s modification approach.Simulation results of a slider-crank mechanism by the proposed method are in good agreement with ones from other literature.