过滤高阶失稳振型满足显式算法的稳定性

摘要

以中心差分显式算法为例,证明了算法传递矩阵的特征向量与结构自振向量的一致性,以及显式算法步长过大时的失稳,是传递矩阵的特征值和特征向量两个因素共同作用的结果.在证明的基础上得到了过滤高阶失稳振型来满足算法稳定性的方法,也就是在每步按通常的显式算法得到计算结果后,增加了从结果位移向量中过滤振型参与系数很小的高阶失稳振型的计算步骤,使得计算步长即使取在稳定域外也能使算法不会失稳.这种方法,既大大提高了通常采用的显式方法的计算效率,同时对求解的精度也影响不大.算例证明了这种方法的有效性.为求解刚性常微分方程组提供一个思路.