各种模态分析方法总结及比较

..各种模态分析方法总结与比较一、模态分析模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参数的过程。模态分析的理论经典定义：将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵，其每列为模态振型。模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的，则称为计算模记分析；如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数，称为试验模态分析。通常，模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性，就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此，模态分析是结构动态设计及设备的故障诊断的重要方法。模态分析最终目标是在识别出系统的模态参数，为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。二、各模态分析方法的总结..（一）单自由度法一般来说，一个系统的动态响应是它的若干阶模态振型的叠加。但是如果假定在给定的频带内只有一个模态是重要的，那么该模态的参数可以单独确定。以这个假定为根据的模态参数识别方法叫做单自由度(SDOF)法n1。在给定的频带范围内，结构的动态特性的时域表达表示近似为：TRRtrQethr2-1而频域表示则近似为：2LRURjQjhrtrrr2-2单自由度系统是一种很快速的方法，几乎不需要什么计算时间和计算机内存。这种单自由度的假定只有当系统的各阶模态能够很好解耦时才是正确的。然而实际情况通常并不是这样的，所以就需要用包含若干模态的模型对测得的数据进行近似，同时识别这些参数的模态，就是所谓的多自由度(MDOF)法。单自由度算法运算速度很快，几乎不需要什么计算和计算机内存，因此在当前小型二通道或四通道傅立叶分析仪中，都把这种方法做成内置选项。然而随着计算机的发展，内存不断扩大，计算速度越来越快，在大多数实际应用中，单自由度方法已经让位给更加复杂的多自由度方法。1、峰值检测峰值检测是一种单自由度方法，它是频域中的模态模型为根据对系统极点进行局部估计(固有频率和阻尼)。峰值检测方法基于这样的事实：在固有频率附近，频响函数通过自己的极值，此时其实部为零(同相部分最..小)，而虚部和幅值最大(相移达90°，幅度达峰值)图1。出现极值的那个固有频率就是阻尼固有频率r的良好估计。相应的阻尼比r，的估计可用半功率点法得到。设1和2分处在阻尼固有频率的两侧(1r2),则：221rjHjHJH2-3rr2122-42、模态检测模态检测是根据频域中的模态模型对复模态(或实模态)向量进行局部估计的一种单自由度方法。在2LRURjQjhrtrrr中略去剩余项则单个频响函数在r处的值近似为：rjrrjrrrrrrjrrrrtjAQjjQjH1112-5由此式可见，频响函数在r处的值乘以模态阻尼因r，就是留数(的估计值如图1。利用这种模态检测方法之前，先要估计出r图1对频响应函数的幅值进行峰值和模态检测3、圆拟合圆拟合是一种单自由度方法，用频域中的模态模型对系统极点和复模..态(或实模态)向量进行局部估计。此方法依据事实是：单自由度系统的速度频响函数(速度对力)在奈奎斯特图(即实部对虚部)上呈现为一个圆。如果把其他模态的影响近似为一个复常数，那么在共振频率r附近，频响函数的基本公式为：1jRjjVUjHrtj2-6因此，首先要选择共振频率附近的一组频率响应点，通过这些点拟合成一个圆。阻尼固有频率r可以看成是复平面上数据点之间角度变化率最大(角间隔最大)的那个点的频率，也可以看成是相位角与圆心的相位角最为接近的那个数据点的频率。对于分得开的模态而言，二者的差别是很小。阻尼比r估计如下：2tan2tan2112rr2-7式中1，2：分居在r两侧的两个频率点：1，2：分别为频率点在1和2得半径与r得半径之间的夹角。圆的直径和阻尼固有频率点的角位置含有复留数U+jV的信息：VUVUrtan,222-8式中：圆的直径：园心与固有频率点的连线跟虚轴之间的夹角．圆拟合法速度也很快，但为避免结果出错，特别是在模态节点附近，需要操作者参与。..（二）单自由度与多自由度系统粘性阻尼单自由度SDOF系统如图2的力平衡方程式表示惯性力、阻尼力、弹性力与外力之间的平衡图2单自由度系统tftKxtxCtxM2-9其中M：质量C:阻尼K：xxx：加速度，速度，位移f：外力t时间变量，把结构中所呈现出来的全部阻尼都近似为一般的粘性阻尼。把上面的时间域方程变换到拉氏域复变量P，并假设初始位移和初始速度为零，则得到拉氏域方程：pFKCpMp2，或pFpXpZZ：动刚度经过变换可得传递函数的定义，pZpH1即pFpHpXMKpMCpMpH///122-10上式右端的分母叫做系统特征方程，它的根即是系统的极点是：MKMCMC/2/2/22,12-11如果没有阻尼C=0，则所论系统是保守系统。我们定义系统的无阻尼固有频率为：MK/12-4临界阻尼Cc的定义为使（2.3）式中根式项等于零的阻尼值：..MKMCc/22-5而临界阻尼分数或阻尼比ζ1为：ζ1=CCc，阻尼有时也有用品质因数即Q因数表示：12/1Q2-6系统按阻尼值的大小可以分成过阻尼系统（ζ11）、临界阻尼系统（ζ1=1）和欠阻尼系统（ζ11）。过阻尼系统的响应只含有衰减成分、没有振荡趋势。欠阻尼系统的响应时一种衰减振动，而临界阻尼系统则是过阻尼系统与欠阻尼系统之间的一种分界。实际系统的阻尼比很少有大于10%的，除非这些系统含有很强的阻尼机制，因此我们只研究欠阻尼的情形。在欠阻尼的情况下式2-11两个共轭复根：111j，11*1j2-7其中1为阻尼因子1为阻尼固有频率。有关系统极点的另外一些关系式有：121111j2-82121112-91112-10212112-112-2式写成如下形式：*11/1ppMpH2-12在展开成部分分式形式，则有：..*1*111pApApH，这里112/1jMA2-13这里的1A和*1A是留数。多自由度系统多自由度系统可以用简单的力平衡代数方程演化成形式相似的一个矩阵的方程。下面是以而自由度系统为例。如图：图3多自由度系统该系统的运动方程如下：tftxKtxKKtxCtxCCxMtftxKtxKKtxCtxCCxM212232122322212212122121112-14写成矩阵形式是212132222121322221212100ffxxKKKKKKxxCCCCCCxxMM2-15或者fxKxCxM2-16其中[M]、[C]、[K]、{f(t)}和{x(t)}分别为质量矩阵、阻尼矩阵、刚度矩阵、方向量和响应向量。把这个时间域的矩阵方程变换到拉氏域（变量为p）且假定初始位移和初始速度为零，则得：pFpXKCpMp22-17..或者是pFpXpZ式中：[Z(p)]动刚度矩阵2-18可以得到传递函数矩阵为：pZpZadjpZpH12-19式中pZadj：pZ的伴随矩阵，等于TijijZ；ijZ：pZ去掉第行第列后的行列式等于奇数如果等于偶数如果jijiij11；传递函数矩阵含有幅值函数。2-19式中的分母，即是pZ的韩烈士，叫做系统的特征方程。与单自由度情况一样，系统特征方程的根，即系统极点，决定系统的共振频率。根据特征值问题，可以求出系统特征方恒的根。为了把系统方程2-17转化为一般的特征值问题公式，加入下面的恒等式：0XMpMp2-20将此式与2-17式结合在一起得：'FYBAp2-21其中CMMA0，KMB00，XXpY，FF0'。如果力函数等于零，那么式2-19就成了关于实值矩阵的一般特征值问题，其特征值马祖下列方程的p值：0BAp2-22它的根就是特征方程0pZ的根。对于N各自由度系统，此方程有2N个呈复共轭对出现的特征根：..NNNNNNjjjj0000\\1111**112-23同单自由度系统一样，多自由度系统的极点的实部r是阻尼因子，虚部r是阻尼固有频率。（三）实模态和复模态按照模态参数（主要指模态频率及模态向量）是实数还是复数，模态可以分为实模态和复模态。对于无阻尼或比例阻尼振动系统，其各点的振动相位差为零或180度，其模态系数是实数，此时为实模态；对于非比例阻尼振动系统，各点除了振幅不同外相位差也不一定为零或180度，这样模态系数就是复数，即形成复模态。1复模态与实模态理论在拟合频段，实模态理论中传递函数在k点激励Z点响应的留数表达式为..1222rjrrklrklrevRHrrrvarctannlk,,2,1,(1)其中，klrR为留数；r和rv构成的复数为系统的复特征值r：rrrjv拟合频段复模态理论中传递函数在k点激励f点响应的留数表达式为nrjrrjklrnrrrjklrklrrreveRvjjeRH122122)(..rrrvarctannlk，，，1（2）由（1）、（2）式中可以看出，传递函数共振峰处复模态的相位与实模态相位的差别在于多出的复留数相位r，由传递函数的逆变换可以得到脉冲响应函数，由此可以得到物理坐标系中结构的自由响应表达式。对于无阻尼结构，t时刻第r阶模态k点的振动为rrrkrkrtYxsin（3）粘性比例阻尼：t时刻第r阶模态k点的振动为rdrtrkrkrteYxrsin（4）一般粘性阻尼：t时刻第r阶模态k点的振动为krrdrtkrrkrteTxmr'cos2(5)式中，φkr表示振型幅值；Ω表示模态频率；θ表示相位角。可以看出，无阻尼和比例阻尼系统的初相位与初始条件有关，与物理坐标无关，具有模态(振型)保持性；而一般粘性阻尼系统的初相位还与物理坐标k有关，每个物理坐标振动时并不同时达到平衡位置和最大位置，不具备模态保持性，是行波形式．但各物理坐标的相位差保持不变