模态分析理论.docVIP

模态分析指的是以振动理论为基础、以模态参数为目标的分析方法。首先建立结构的物理参数模型，即以质量、阻尼、刚度为参数的关于位移的振动微分方程；其次是研究其特征值问题，求得特征对（特征值和特征矢量），进而得到模态参数模型，即系统的模态频率、模态矢量、模态阻尼比、模态质量、模态刚度等参数。 特征根问题 以图3所示的三自由度无阻尼系统为例，设，， 图 SEQ 图 \* ARABIC 3 三自由度系统 其齐次运动方程为： mz+kz=0 （ 其中m，k分别为系统的质量矩阵和刚度矩阵，，，则运动方程展开式为： （9） 定义主振型 由于是无阻尼系统，因此系统守恒，系统存在振动主振型。主振型意味着各物理坐标振动的相位角不是同相（相差0o）就是反相位（相差180o），即同时达到平衡位置和最大位置。主振型定义如下： （10） 其中zi为第阶频率下，各自有度的位移矢量，zmi为第个特征矢量，表示第阶固有频率下的振型，为第阶频率下的第个特征值，为初始相位。 对于三自由度系统，在第阶频率下，等式可以写成 （11） 表示第个自由度在第阶模态下的模态矩阵。 特征值 对式（10）二次求导，得 （12） 代入齐次运动方程得 m-ωi2 去除sinωi k-ωi2 以矩阵的形式展开得： （15） zmi有非零解 （16） 即 （17） 方程解如下：，，。三个解对应该系统的前三阶固有频率，每一个特征根对应一个特征矢量zi，表示对应模态下该系统的振型。 特征矢量 由式 k （18） 展开第一行和第二行，忽略下脚标m和i，得 （19） 得 （20） 如果设定了值，则就可以求出三个特征根值下，和相对于的位移。假设， 一阶模态，：，，即z1=11 二阶模态，：，，即z2=10 三阶模态，：，，即z3=1- 模态矩阵 所谓模态矩阵就是指各列由各阶模态特征矢量构成的矩阵，如图4所示。 图 SEQ 图 \* ARABIC 4 模态矩阵 对于前面提到的三自由度系统，模态矩阵如下： z 运动方程的解耦 对于一个复杂的系统，在物理坐标系统中建立的运动方程之间存在耦合关系，因此求解起来比较麻烦，因此需要进行坐标系转化，将耦合的运动方程变为非耦合的运动方程，再将求得的结果转化为物理坐标系下的结果，运动方程解耦过程如下图5： 图 SEQ 图 \* ARABIC 5 运动方程解耦过程 在进行坐标变换之前需对刚度矩阵和质量矩阵进行归一化。任意上面的三自由度系统为例，由式m-ω ωi2mz ωj2mzm 对式（21）左乘zmj ωi2zmj 又因为 ω 因为系统对称所以，mT ωj2zmjT 对式（24）右乘z ωj2zmjT 则式（23）—式（25）得 (ωi2- 当(ω zmjTmz 当ωi2-ωj2=0 zmjTmz 则对质量矩阵和刚度矩阵的归一化结果如下： mn=zmT kn=zmT 特征矢量的归一化 由于特征矢量只是位移之比，而不是绝对振幅，因此可以对其进行归一化处理。令zn zni=z qi=j=1 对于对角质量矩阵 qi=j=1 则三自由度系统： zm=11 （35） 则归一化的质量矩阵为 （36） 同理归一化后的刚度矩阵为 （37） 可以看出归一化后的刚度矩阵对角线上的各项就是各阶模态固有频率的平方。 运动方程解耦 将物理坐标系下的运动方程按照前面介绍的归一化方法转化为主坐标系下的运动方程，其结果如下： 38） 可以看出在主坐标系中的运动方程之间没有耦合关系，分别单独描述各阶模态的运动特性。 初始条件和激励的坐标转换 物理坐标系中的非齐次运动方程为 （39） 做如下变形 （40） 其中，就是前面介绍的质量和刚度矩阵的对角化。 令，主坐标质量矩阵； ，主坐标刚度矩阵； ，主坐标系加速度矢量； ，主坐标系位移矢量； ，主坐标系激励矢量。 同样的关系也适用于初始位移和速度：