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

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

§4.1 概述

根据结构振动的数学模型，即振动微分方程所形成的矩阵特征值问题，求解结构的固有振动特性——固有频率与固有振型，是结构振动分析的一个主要任务。结构的固有振动特性是结构振动的内因。固有振动特性也是进行结构振动响应分析和结构动力学设计的基础。

对于简单的结构，如均匀直梁、均匀直杆等，可以用解析的方法解得其固有振动特性。对于一般结构，如果只需获得结构有限阶的固有振动特性，也可以采用试验测试（模态识别）的方法来获得。

但是对于大型复杂结构，不可能用解析分析的方法得到其固有振动特性，而采用试验测试的方法不仅花费高，而且周期长，对于处于设计状态的结构，显然也无法进行试验。

所以对复杂的工程结构，常用的方法是建立结构的数学模型，用数值求解的方法获得结构的固有振动特性。

随着计算技术飞速发展和特征值计算方法的研究进展，通过矩阵特征值问题的求解来获得结构固有振动特性，是已经被振动工程界普遍接受的一个有效和可靠的途径。从数学理论上也可以证明，许多特征值计算方法具有相当好的精度，并且获得了实践和实验的证明。

由于结构固有振动特性求解与矩阵特征值求解问题的密切关系，在结构振动分析中，矩阵特征值问题已经成为结构固有振动特性分析的一个代名词。

所以在本章中，只要不作说明，一般讲的矩阵特征值问题就是指结构的固有振动特性求解问题。所谓系统的特征值就代指结构的固有频率，特征向量代指

结构的固有振型（固有模态）

矩阵特征值问题的数值求解方法可以分为三类：矩阵分解法、迭代法和矩阵变换法。由于矩阵（代数）特征值问题本身就是一个完整的系统，本章只能根据结构固有振动分析问题的需要，介绍一些常用的求解方法。详尽的矩阵特征值问题的数值求解方法可以参考威尔金森的名著《代数特征值问题》。

本章的论述是建立在已经用有限元素法建立了结构振动运动数学模型的基础上。

§4.2 结构振动特征值问题的性质

根据结构振动方程，可以得到结构固有振动的代数特征值问题：

}]{[}]{[2x M x K ω=

（4－1）

或 }]{[}]{[x M x K λ= （2ωλ=） （4－2）

振动特征值问题除了第二章所述的性质外，在特征值问题的数值求解中，还要用到如下一些性质： 1． 移轴特性

对特征值问题

}]{[}]{[x M x K λ= （4－3）

若μ为一已知实数，则有：

}]{)[(}]){[]([x M x M K μλμ-=- （4－4）

新的特征值问题可写为：

}]{[}]{ˆ[x M x K

ρ= （4－5） ][][]ˆ[M K K

μ-= （μλρ-=） （4－6） 显然，上面两个特征值问题具有相同的特征向量，而特征值间的关系为：

μρλ+=i i （4－7）

μ称为移轴量。

在结构振动分析中，移轴特性常用来消除刚度矩阵的奇异性，也可用来加速迭代求解的收敛速度。 2． 特征值对合同变换的不变性

［合同变换］：若][A 为一个非奇异矩阵，则变换]][[][A M A T 称为对][M 矩阵的合同变换。

对矩阵特征值问题：

}]{[}]{[x M x K λ= （4－8）

中的质量阵和刚度阵，用矩阵][A 作如下合同变换：

]][[][]ˆ[A M A M

T = （4－9） ]][[][]ˆ[A K A K

T = （4－10） 则得到一个新的特征值问题

}]{ˆ[}]{ˆ[φμφM K

= (4－11) 从而有：

][][][]

[]][[][]][[][]ˆ[]ˆ[=-=-=-A M K A A M A A K A M K T

T T μμμ (4－12)

由于0][≠A ，故有

0][][=-M K μ （4－13）

显然变换后的特征值不变：

μλ= （4－14）

且可以容易地证明，变换前后两个特征值问题的特征向量之间具有关系：

}]{[}{φA x = （4－15）

即合同变换不改变矩阵的特征值，特征向量具有转换关系。

3． 特征值对相似变换的不变性

［相似变换］：对非奇异阵][A ，变换]][[][1A K A -称为相似变换。 对矩阵特征值问题：

}]{[}]{[x M x K λ= （4－16）

中的质量阵和刚度阵，用矩阵][A 作如下相似变换：

]][[][]ˆ[1A M A M

-= （4－17） ]][[][]ˆ[1A K A K

-= （4－18） 则得到一个新的特征值问题：

}]{ˆ[}]{ˆ[φμφM K

= （4－19） 从而有：

][][][]

[]][[][]][[][]ˆ[]ˆ[1

11=-=-=----A M K A A M A A K A M K μμμ （4－20）

由于0][≠A ，故有

0][][=-M K μ （4－21）

显然变换后的特征值不变：

μλ= （4－22）

且可以证明，变换前后的特征向量间具有关系：

}]{[}{φA x = （4－23）

即相似变换不改变矩阵的特征值，特征向量具有转换关系。

4． 特征值对正交变换的不变性

［正交变换］：在相似变换中，若对于非奇异矩阵][A ，有][][][I A A T =，即

1][][-=A A T ，则用][A 进行的相似变换称为正交变换。正交变换不仅不改变矩阵

的特征值，而且正交变换后矩阵的对称性不变。

5． 特征值对旋转变换的不变性

［旋转变换］：取变换矩阵为：

⎥⎦

⎤

⎢

⎣⎡-=θθθθ

cos sin sin cos ][R （4－24） 可以证明][R 为正交矩阵，用][R 矩阵进行的相似变换，称为旋转变换。 通过一系列的旋转变换，可以使对称阵变为对角阵，使一般矩阵变为三角阵。 利用上述特征值问题的性质，可以得到不同的特征值求解方法。 结构振动特征值的数值求解中常用的方法主要有： 【多项式迭代法】

多项式迭代法是利用特征值使特征多项式等于零的性质，即

0][][)(=-=M K p λλ （4－25）

【利用特征多项式的Sturm 序列性质的方法】

通过Sturm 序列的性质，分离出每个特征值所在的区间，通过二分法和其它加速查找方法，逐步缩小特征值所在的分隔区间，最后得到所需的特征值。 【矢量迭代法】

矢量迭代法利用的基本关系式为：

}]{[}]{[i i i x M x K λ= （4－26）

具体的迭代方法有：正迭代法、逆迭代法、子空间迭代（同时迭代）法。 【矩阵变换法】

矩阵变换法利用的是“矩阵特征值对合同变换的不变性”这一性质。利用具有某种性质的变换矩阵，经过一系列变换，得到一个形式简单的、容易进行特