结构动力方程的辛RK法和辛两步法.PDFVIP

结构动力方程的辛RK 法和辛两步法 侯小平 邢誉峰 （北京航空航天大学 固体力学研究所，北京 100191） 摘要 线性无阻尼系统的单步相位误差由频率和积分步长的乘积决定。只要频率和积分步长相同，任意的无阻尼 线性系统的主振动都有相同的单步相位误差，即多自由度线性系统某一阶频率的相位误差的求解方法与求具有相同 频率的单自由度系统的相位误差的方法相同，这能明显提高相位误差的分析效率。然而，如何利用各阶频率的单步 相位误差，对相位误差进行修正以提高数值积分的精度，有待于更深入的研究。 本文以线性无阻尼动力学方程为例， 研究了辛 Runge-Kutta （RK ）法和辛两步法的相位误差，前者相位滞后，后者相位超前。理论动态响应应该位于两方 法结果之间，由此来改变算法的时间坐标，会使数值模拟精度得到提高。 关键词 结构动力方程，Runge-Kutta 法，辛两步法，相位误差分析 引 言 实验、理论和模拟方法是科学研究的主要方式，计算可以克服理论分析及实验手段的局限，并 且随着计算机的发展其地位不断上升。算法是计算数学领域的重要问题之一, 特别是对于非线性问 题,算法的研究更为重要。 无耗散的物理过程与现象，可以用Hamilton系统来描述。辛性是Hamilton系统的基本特征之一。 然而在用传统数值方法模拟Hamilton系统时，如差分方法（如中心差分法、Houbolt法）、线性加速 度法（如Wilsonθ，Newmark法）存在明显的弊端。这些方法只有一阶或二阶精度，并且这些方法均 有算法阻尼。冯康[1]提出的关于Hamilton系统的辛算法，解决了算法的耗散问题。钟万勰[2,3]将Hamilton 体系引入弹性力学中，开创了Hamilton系统在弹性力学中应用的先例。 辛算法虽然可以保结构，但其相位误差问题依然存在，并最终导致动态响应长期数值模拟的失 败。钟万勰[3]对线性定常结构动力系统提供了精细时程积分法，可以得到计算机精度意义下的解析 解。Wu[4]等用微分求积法求解了具有线性和非线性特性的边值和初值问题，提高了相位误差精度。 [5] [6] Franco 等用新的SDIRKN法对二阶常微分系统进行积分，重点研究了相位误差问题。邢誉峰 等给 出了Euler中点辛差分格式相位误差的精确估计公式，提出了简单的修正方法。 如何进一步提高辛算法的数值模拟精度，本文从相位角度进行了探讨。研究了四阶辛RK法[1]和 四阶辛两步法[7，8]的相位误差，并给出了提高模拟精度的方法。 1 结构动力方程的数值积分 考虑如下动力学方程 Mx + Cx + Kx = r t (1) ( ) 式中 M 、C、K 分别为 n×n 的质量矩阵、反对称陀螺阵和刚度矩阵，x 为位移向量，r(t)为激励向量。 1.1 RK 法及其辛对偶性 单步 s-级 RK 法都具有如下格式： s z k +1 z k =+h∑b f (Y ), i i i 1