子空间迭代法的计算加速研究.pdfVIP

中文科技期刊数据库（引文版）工程技术

中文科技期刊数据库（引文版）工程技术

子空间迭代法的计算加速研究

林庚浩

中国建筑科学研究院有限公司北京构力科技有限公司，北京100013

摘要：采用Wilson移频策略对子空间迭代法进行了加速，给出了详细的移频子空间迭代求解特征值的步骤。为

加速高阶特征值的收敛，对Wilson移频策略进行了改进。讨论了提高子空间迭代计算稳定性的措施和给定有效

质量系数下由程序自动确定振型数的方法。从对比算例可看出，采用移频算法后子空间迭代法求解特征值加速明

显。本文研发的特征值求解器理论上可以求解出所有阶的特征值，且计算时间随特征值数的增加近似呈线性增长

关系，对求解大阶数特征值问题有一定的意义。

关键词：子空间迭代法；移频算法；计算稳定性；自动确定振型数

中图分类号：TU201

广义特征值问题在建筑、机械、电磁学、声学、阵阶数比较小的特征值问题。在这些加速方案中，移

量子化学等领域有着广泛的应用。结构系统的频率和频算法被认为是子空间迭代法用来计算大规模固有振

振型是结构动态分析的重要特性。通常将频率和振型动特征值问题的最好选择。本文采用移频算法，使得

称为特征对，因此，求解结构的广义特征值问题就是求解大型稀疏矩阵的广义特征值问题的效率得以提高。

要求解所需的特征对。在结构有限元分析中，由于实1移频加速

际工程需要，经常需要求解大规模自由度的广义特征

[1]1.1移频算法

值问题，计算规模已达到了30万自由度甚至更大。

移频又称谱变换(SpectrumTransfor-mation)或

随着研究问题的深入，求解数个低阶特征值已经

移位，相当于在频率谱上对频率作了一维坐标变换，

不能完全满足要求，如在建筑结构设计中，为了达到

从而加速待求的特征值和特征向量的收敛速度。

规范90%的有效质量系数的要求，需要计算数十到数百

广义特征值问题如下：

个特征值。因此需要寻求高效、稳定的特征值求解算

Kf-lMf=0(1)

法来提高单机上的计算规模和速度。

移频后，特征值问题变为

子空间迭代法、Ritz向量法和Lanczos法被公认

¢

(K-mM)j=lMj(2)

为求解部分低阶特征对的有效方法。通常认为Ritz向

其中，l¢=l-m，特征值移频后，特征向量不

量法和Lanczos法具有比子空间迭代法更高的计算效

变。

率，Ritz向量法和Lanczos法比子空间迭代法平均快

在标准的子空间迭代法中，第i个特征值的收敛

4-10倍。大量算例表明，在允许模态误差的意义之下，

l/ll/l