快速子空间迭代法、迭代Ritz向量法与迭代Lanczos法的比较.pdf

第18卷第2期 振 动 工 程 学 报 V01．18No．2 ofVibration Jun．2005 2005年6月 Journal Engineering 快速子空间迭代法、迭代Ritz向量法与 迭代Lanczos法的比较 宫玉才，周洪伟，陈 璞，袁明武 (北京大学力学与工程科学系，北京100871) 摘要：以高效的细胞稀疏直接快速解法为核心步骤，实现了快速的固有振动广义特征值问题解法。并在相同的允 许模态误差的意义下检验了三种常用的大型矩阵特征模态算法——子空间迭代法、迭代Ritz向量法和迭代Lanc— ZOS法的计算效率。迭代Ritz向量法平均最快。子空间迭代法最慢，三种解法效率相差不是太大。与ANSYS的子空 较大规模的例题显示，本文对特征值算法的改进是十分有效的，算法的健壮性，通用性都达到了高水平。 关键词：结构振动；特征值；子空间；迭代法 文献标识码：A 文章编号：1004-4523(2005)02一0227—06 中图分类号：0327；0151．21 1振动特征值问题 向量法和Lanczos方法对收敛的判定是相对含糊 的，在实际工程计算中可能造成漏根或多根。 在工程有限元分析中常常要求解广义代数特征 传统上，子空间迭代用特征值的两次迭代的相 值问题 对误差不等式 K妒一剐订P一0 (1) 的部分低阶特征值与特征向量。对于矩阵阶数超过 监丢I掣≤b】(卜P1)l ＼。^ 、o(3)7 1000的大型问题，子空间迭代法、Ritz向量法和 控制收敛，而Lanczos法用其过程中的不等式 Lanczos法被公认为求解部分低阶极端特征值和特 I酊1一可j I％ (4) l≤I成s。i 征向量的有效方法。尽管国内外的有限元软件都提 控制收敛。在大量的工程计算中，发现在允许误差 供广义代数特征值问题(1)的多种解法，但结果仍然 e。=￡。=10叫的情形下，除最低的10多阶模态之外， 不能令人完全满意，漏根与多根、自由模态误判都时 子空间迭代与Lanczos法所得到的特征向量精度都 有发生。 可能不令人满意。这一现象对Lanczos方法尤为严 传统上，低端特征值问题求解过程极度依赖于 重，原因是采用逆迭代技术时，高阶的、密集的特征 谱变换的线性方程组 值不易分离。 (X—fYM)x=LDLTz=My(2) 关于特征模态的收敛，不同的算法往往采用不 的解法，移轴矩阵x一卢M的LDLl三角分解是计算 同的标准，相对速度的比较不是很客观。对于特征模 量最大的。在以变带宽解法为核心步骤的特征值解法 态的近似(瓦，磊)，在各种算法中可以统一用模态误 中，它常常占到特征值问题计算时间的70％～90％。差(5)代替特征值误差作为收敛判据，来衡量算法 本文采用了文Ef]提出的一个效率非常高的有限元解 的效率。 法——细胞稀疏直接快速解法(简称细胞解法)替换 等≈锗≤％㈣ 变带宽解法，极大地提高了三角分解的效率。 如果要求不太