第四章_结构固有振动特征值问题的数值解.doc

第四章 结构固有振动特征值问题的数值解 §4.1 概述 根据结构振动的数学模型，即振动微分方程所形成的矩阵特征值问题，求解结构的固有振动特性——固有频率与固有振型，是结构振动分析的一个主要任务。结构的固有振动特性是结构振动的内因。固有振动特性也是进行结构振动响应分析和结构动力学设计的基础。 对于简单的结构，如均匀直梁、均匀直杆等，可以用解析的方法解得其固有振动特性。对于一般结构，如果只需获得结构有限阶的固有振动特性，也可以采用试验测试（模态识别）的方法来获得。 但是对于大型复杂结构，不可能用解析分析的方法得到其固有振动特性，而采用试验测试的方法不仅花费高，而且周期长，对于处于设计状态的结构，显然也无法进行试验。 所以对复杂的工程结构，常用的方法是建立结构的数学模型，用数值求解的方法获得结构的固有振动特性。 随着计算技术飞速发展和特征值计算方法的研究进展，通过矩阵特征值问题的求解来获得结构固有振动特性，是已经被振动工程界普遍接受的一个有效和可靠的途径。从数学理论上也可以证明，许多特征值计算方法具有相当好的精度，并且获得了实践和实验的证明。 由于结构固有振动特性求解与矩阵特征值求解问题的密切关系，在结构振动分析中，矩阵特征值问题已经成为结构固有振动特性分析的一个代名词。 所以在本章中，只要不作说明，一般讲的矩阵特征值问题就是指结构的固有振动特性求解问题。所谓系统的特征值就代指结构的固有频率，特征向量代指结构的固有振型（固有模态） 矩阵特征值问题的数值求解方法可以分为三类：矩阵分解法、迭代法和矩阵变换法。由于矩阵（代数）特征值问题本身就是一个完整的系统，本章只能根据结构固有振动分析问题的需要，介绍一些常用的求解方法。详尽的矩阵特征值问题的数值求解方法可以参考威尔金森的名著《代数特征值问题》。 本章的论述是建立在已经用有限元素法建立了结构振动运动数学模型的基础上。 §4.2 结构振动特征值问题的性质 根据结构振动方程，可以得到结构固有振动的代数特征值问题： （4－1） 或 （） （4－2） 振动特征值问题除了第二章所述的性质外，在特征值问题的数值求解中，还要用到如下一些性质： 移轴特性 对特征值问题 （4－3） 若为一已知实数，则有： （4－4） 新的特征值问题可写为： （4－5） （） （4－6） 显然，上面两个特征值问题具有相同的特征向量，而特征值间的关系为： （4－7） 称为移轴量。 在结构振动分析中，移轴特性常用来消除刚度矩阵的奇异性，也可用来加速迭代求解的收敛速度。 特征值对合同变换的不变性 ［合同变换］：若为一个非奇异矩阵，则变换称为对矩阵的合同变换。 对矩阵特征值问题： （4－8） 中的质量阵和刚度阵，用矩阵作如下合同变换： （4－9） （4－10） 则得到一个新的特征值问题 (4－11) 从而有： (4－12) 由于，故有 （4－13） 显然变换后的特征值不变： （4－14） 且可以容易地证明，变换前后两个特征值问题的特征向量之间具有关系： （4－15） 即合同变换不改变矩阵的特征值，特征向量具有转换关系。 特征值对相似变换的不变性 ［相似变换］：对非奇异阵，变换称为相似变换。 对矩阵特征值问题： （4－16） 中的质量阵和刚度阵，用矩阵作如下相似变换： （4－17） （4－18） 则得到一个新的特征值问题： （4－19） 从而有： （4－20） 由于，故有 （4－21） 显然变换后的特征值不变： （4－22） 且可以证明，变换前后的特征向量间具有关系： （4－23） 即相似变换不改变矩阵的特征值，特征向量具有转换关系。 特征值对正交变换的不变性 ［正交变换］：在相似变换中，若对于非奇异