Energy adjusting algorithm for the nonlinear dynamic equation with multi-degrees of freedom

An energy adjusting numerical algorithm for the nonlinear dynamic equation with multi-degree of freedom is proposed. First, by the Taylor's expansion and Duhamel integration, an integral iteration formula with an undetermined parameter is obtained. Second, the original dynamical equation is transformed into an energy adjusting equation to determine the undetermined parameter. Finally, substituting the parameter into the integral iteration formula, an accurate numerical value is obtained. Some examples show that the method can eliminate the algorithm damping and enjoys better stability than the Runge-Kutta method under a large integral step.