Runge-Kutta方法求解结构动力学方程

将几种具有不同稳定性的Runge-Kutta方法应用到结构动力学方程的数值求解中。针对增量形式的动力学方程,使用改进的Newton-Raphson迭代,研究了减少计算量的两种方法：（1）使用单对角隐式Runge-Kutta方法,（2）应用转化矩阵。采用逼近算子的谱半径分析了稳定性与数值阻尼特性,解释了L-稳定方法抑制高频振荡的原因。数值算例表明在精确解上较小的物理阻尼能有效的抑制高频振荡,但对各种直接积分方法的影响很小,高精度的L-稳定Runge-Kutta方法能在有效抑制高频振荡的同时高精度的求解低频振动。

Abstract：

Several Runge-Kutta methods with the different stability were applied to solve the equations of motion in structural dynamics. For incremental dynamical equations,using the modified Newton-Raphson iteration,two methods to reduce the amount of work were proposed. The first one is the singly diagonally implicit Runge-Kutta methods,and the second one is to apply the transform matrix. Using the spectral radii of approximation operators,the stability analysis and the numerical damping property were studied,and the reason why the L-stability methods could wipe out the high oscillations was explained. Numerical example was solved by several direct integration methods,the result show that the small physical damping can wipe out high oscillations effectively on exact solution,but it has little effect on numerical solution,and the high order L-stability Runge-Kutta methods can wipe out the high oscillation effectively,at the same time,solve the vibration of low frequencies with high accuracy.

Runge-Kutta Methods for Time Integration in Structural Dynamics

Several Runge-Kutta methods with the different stability were applied to solve the equations of motion in structural dynamics. For incremental dynamical equations,using the modified Newton-Raphson iteration,two methods to reduce the amount of work were proposed. The first one is the singly diagonally implicit Runge-Kutta methods,and the second one is to apply the transform matrix. Using the spectral radii of approximation operators,the stability analysis and the numerical damping property were studied,and the reason why the L-stability methods could wipe out the high oscillations was explained. Numerical example was solved by several direct integration methods,the result show that the small physical damping can wipe out high oscillations effectively on exact solution,but it has little effect on numerical solution,and the high order L-stability Runge-Kutta methods can wipe out the high oscillation effectively,at the same time,solve the vibration of low frequencies with high accuracy.