三种常用固有振动特征值解法的比较.pdf

2005 全国结构动力学学术研讨会 海南省海口市，2005.12.19-20 中国振动工程学会结构动力学专业委员会 三种常用固有振动特征值解法的比较 宫玉才1 周洪伟 陈 璞 袁明武 （北京大学力学与工程科学系 北京，100871） Email：yuanmw@pku.edu.cn 摘要： 本文以高效的细胞稀疏直接快速解法为核心步骤，实现了快速的固有振动广义特征值问题解法, 并 在相同的允许模态误差的意义下检验了三种结构动力学中常用的大型矩阵特征模态算法——子空间迭代法、 迭代 Ritz 向量法和迭代 Lanczos 法的计算效率。迭代 Ritz 向量法平均而言最快，子空间迭代法最慢，三种 解法效率相差不是太大。与 ANSYS 的子空间迭代和 Lanczos 法相比，本文的子空间迭代比 ANSYS 的效率 高很多，Lanczos 法和 ANSYS 的差不多 。大量较大规模的例题显示，本文对特征值算法的改进是十分有效 的，算法的健壮性，通用性都达到了高水平。 关键词：特征值，结构振动，迭代法，高效能计算 引言 在工程有限元分析中常常要求解广义代数特征值问题 KM? ?=λ?0 (1) 的部分低阶特征值与特征向量。对于矩阵阶数超过1000的大型问题，子空间迭代法、Ritz向量法和 Lanczos法被公认为求解部分低阶极端特征值和特征向量的有效方法。尽管国内外的有限元软件都提 供广义代数特征值问题(1)的多种解法，但结果仍然不能令人完全满意，漏根与多根、自由模态误判 都时有发生。 传统上，低端特征值问题求解过程极度依赖于谱变换的线性方程组 ()K ?=μMx LDLxT = My (2) 的解法，移轴矩阵 K ? μM 的LDLT三角分解是计算量最大的主要步骤。在以变带宽解法为核心步骤 的特征值解法中，它常常占到特征值问题计算时间的70%到90%。本文采用了文[1]提出的一个效率 非常高的有限元解法-细??稀疏直接快速解法（简称细胞解法）替换变带宽解法，极大地提高了三角 分解的效率。 如果要求不太多的特征模态，例如10个，通常认为Ritz向量法和Lanczos法具有比子空间迭代法 更高的计算效率，Ritz向量法和Lanczos法比子空间迭代法平均快4～10倍[2]。但是，标准的Ritz向量 法和Lanczos方法对收敛的判定是相对含糊的，在实用的工程计算中可能造成漏根或多根。 1高等学校博士学科点专项科研基金资助项目 （编号：20030001112） 172 传统上，子空间迭代用特征值的两次迭代之相对误差不等式 ()ll (+ 1) ||λλkk? (1)l+ ≤ ε λ (3) ||λk 控制收敛，而Lanczos法用其过程中的不等式 ?1 ||||λki?≤νβ qqis εν (4) ?5 控制收敛。在大量的工程计算中，发现在允许误差εελν==10 的情形下，除最低的十几阶模态之 外，子空间迭代与Lanczos法所得到的特征向量精度都可能不令人满意。这一现象对Lanczos方法尤