基于matlab的发动机悬置系统的固有频率和主振型计算20542.docVIP

基于Mat lab的发动机悬置系统的固冇频率和主振型计算（二） 3运用MATLAB对动力总成悬置系统固有特性的计算 3. 1理论计算 动力总成系统固有特性的计算，即计算系统的固有频率和振型。动力总成悬直系统无阻 尼的自由振动微分方程： 式中：M——对称正定惯性矩阵； K——对称正定刚度矩阵。 求多自山度振动系统的固有频率，从数学上讲就是求特征值的问题： 设式（13）的解为:X=Xsin （讥+a） 代入式（13）化简后得：KX二3 2MX 左乘吐 1 得：M- 1KX= w2X （14） 令 M- 1K=A,则:AX二32X （15） 32即为A阵的特征值，X为其特征向量。 由于M对称正定，K也是对称阵，因而式（13）是广义特征值问题。可用广义特征值 的方法求得特征值及特征向屋，所求特征值即为系统的固冇频率。 3. 2 MATLAB计算过程 Mat lab是Matrix Laboratory （矩阵实验宗）的缩写，它是山美国Mathwork公司于 1967年推出的软件包，已发展为一种功能强大的计算机语言，特别适合于科学与工程计算。 （1） 将动力总成系统质量参数代入式（6）可得惯性矩阵M。 （2） 将各悬置点的位置参数及悬置块的主刚度代入，可得EiBiDi。再根据式（12）求 得总体的刚度矩阵K。 （3） 编制Matlab程序,由上述（1）、（2）得到矩阵M, K,由式（14）、（15）即可求 得Ac 由式(15),通过Matbb命令eig (A)，即可求出矩阵A的特征值宀2。利用 公式32二2 n f,即可得到悬置系统的各个振动固有频率fo 4振动占优方向的判定 在系统定处标系中，根据系统的质量矩阵[M]及振型矩阵，nJ以求出系统在做各阶主振 动时的能量分布，将它写成矩阵形式，定义为能量分布矩阵[EG] j0当系统以笫j阶固有频 率振动时，此矩阵的(k, j)元素为: 式中： [M]kl——质量矩阵的(k, j)元素； (u ( j) }k——第j阶振型列阵的第k个元素； {u ( j) }1——第j阶振型列阵的第1个元素； 3j——为第j阶固冇频率。 根据此矩阵各行元素值总和的大小便可以判別出在系统以笫j阶固有频率振动时的占 优方向。 在确定了系统的固有频率和固有振型以及相应的振动占优方向后，便可以从避开共振频 率这一理论来初步评价系统的隔振性能。 目录 文尾口口口 m工程地质计算机应用2004年1期 总33期 应用MATLAB计算结构自振频率和振型的一种方法 关文阁杨黎萌魏翠玲(河北工程学院土木工程系河北邯郸056038) 【摘要】木文从质量归一化原理出发，采用MATLAB工具箱屮eig函数，推导出一个新的求 解多质点弹性体系白振频率和振型的方法。算例表明该方法比传统的雅可比方法更简捷更易 于应用。 【关键词】自振频率振型雅可比法特征值法Mat lab 当采用地震反应谱方法计算水平地震作用标准值时，需要求得模型结构的多个主振型及其相 应的自振频率(周期)。因此计算多质点弹性体系的自山振动(包括自振周期、振型等)是 进行结构抗震设计的必要步骤。 1基本理论 多质点弹性体系的无肌尼自山振动方程 (1) 式中:为体系刚度矩阵；为体系质量矩阵。 方程(1)左乘整理后得 (2) 令则有 (3a) 或 (3b) 这是一个求特征值和特征向量问题。该方程非零解的充要条件，是它的系数行列式等于零， 即 ⑷ 式(4)称为方阵 的特征方程。称为方阵 的特征值或特征根。将所求得的 个 逐个回 代到式(3b),解出，称为方阵与相对应的特征向量，也就是所要求解的第振型。 因此，求解体系的自振频率与振型问题也就是求解方阵的特征值和特征向量问题。再利丿IJ 公式(4)求出体系的频率。 (5) 1. 1雅可比法 雅可比(Jacobi)法可以求解实对称阵的特征值和特征向量。对式⑷首先需要把 作对称化处 理。 令 ⑹ 式中 代入式(2),得 等式两边左乘整理后得 即 (8) 式屮 为一实对称矩阵。矩阵 与 有相同的特征值，但它们的特征向量之间存在式⑷的关 系，所以，振型可用式(7)求得，即 (9) 用雅可比方法求解实对称矩阵的特征值和特征向量的基本原理是，寻找一个正交矩阵， 使 (10) 式中，为一对角短阵。这时矩阵的个对角元素就是对称矩阵的个特征值；正交矩 阵 中的第 列就是与对角矩阵 中第 个对角元素对应的特征向呆。所以，雅可比法实际 上是运用平面旋转变换的方法消去矩阵中的非对角元素使其对介化。具体步骤参见文献 [1]。 1. 2本文方法 本文采用MATLAB程序Ieig函数求 特征矩阵 和特征值，应用质量归一-化原理对特征矩 阵 变形就得到体系的振型矩阵。再按(5)求得口振频率。 由质量归一化 (11) 令振型矩阵表示为 (12) 把(⑵带入(11)得