原创力文档- 1、本文档共22页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
特征值解法
《结构动力学》大作业
结构大型特征值问题的求解
0810020035 吴亮秦
1振动系统的特征值问题
1.1实特征值问题
n自由度无阻尼线性振动系统的运动微分方程可表示为：
（1.1）
其中，是位移向量，和分别是系统的质量矩阵和刚度矩阵，都是n阶正定矩阵，是激励向量。
此系统的自由振动微分方程为
（1.2）
设其主振型为： （1.3）
其中，为振幅向量，ω为圆频率，为初相位。将(1.3)代入自由振动微分方程(1.2)，
得：
（1.4）
其中，(1.4)具有非零解的条件是
(1.5)
式(1.4)称为系统的特征方程，由此可以确定方程的n个正实根，称为系统的特征值，称为系统的固有频率，（i=1,2,…..n）为对应于特征值的特征向量或称为系统的振型或模态。
因为矩阵正定，则有Cholesky分解：
(1.6)
其中，是下三角矩阵。引入向量满足：，则：
(1.7)
代入(1.4)，得：
(1.8)
其中，，式(1.8)称为标准实特征值问题。
1.2复特征值问题
多自由度阻尼自由振动系统的运动方程为如下二阶常系数微分方程组：
(1.9)
其中，，分别是n阶的质量、阻尼和刚度矩阵，是n维可微向量函数。用分离变量法，设，其中是与时间无关的常向量，为待定参数。将代入上述齐次方程，得确定参数,的特征方程：
(1.10)
(1.10)具有非零解的条件是
(1.11)
式(1.11)的求解就是复特征值问题。
2 实特征值求解方法
2.1特征方程法
求解实特征值最直接的方法就是特征方程法，即把式(1.5)展开得到特征值多项式：
(2.1)
求解(1.10)即得特征值，特征方程法仅适用于的低阶情况特征值求解，并不是求解特征值的一般方法，实际求解大型结构的实特征值问题的方法很多，归纳起来大致分为两类，即基于矩阵相似变换原理的相似变换法和基于瑞利—李兹（Rayleigh-Ritz）法、迭代方法的瑞利—李兹类方法。
2.2相似变换法
相似变换法基于矩阵的相似变换原理，即相似矩阵与原矩阵有相同的特征值，主要有以下几种方法。
2.2.1雅克比法（Jacobi）
雅克比法是求解实对称矩阵全部特征值的简单有效的方法。它的基本思想是，通过正交
相似变换使矩阵对角化，从而求出矩阵的全部特征值，进行步骤如下：
将系统的自由振动微分方程：
化为标准实特征值形式：；
(2)找矩阵[P]中绝对值最大的非对角线元素，构造正交矩阵，对[P]矩阵作正交变换使元素化成0；
（3）找矩阵中的绝对值最大非对角线元素，构造正交矩阵，对矩阵作正交变换使元素化成0；
（4）依次进行下去，每次找矩阵中的最大非对角线元素，构造正交矩阵，对矩阵作正交变换使元素化成0。
重复使用该变换，每次变换可使矩阵更接近于对角矩阵，若干次变换厚，原矩阵[P]化成对角矩阵，对角线元素即是原矩阵的特征值。
每步变换的关键在于构造正交矩阵，实际采用吉文斯（Givens）旋转矩阵，通过多次坐标系的旋转来实现原矩阵的对角化，取如下形式：
其中，u、v为中最大非对角线元素所在的行列序号，为旋转角，由正交变换使矩阵的元素==0来确定，计算得到：
2.2.2 豪斯厚德三对角化法（Householder）
当矩阵的阶数较高，非对角元素较稠密，且数值较大时，使用雅克比法的收敛
文档评论（0）