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

1、各种模态分析方法总结与比较一、模态分析模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参 数的过程。模态分析的理论经典定义：将线性定常系统振动微分方程组中的物理 坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描 述的独立方程，以便求出系统的模态参数。坐标变换的变换矩阵为模态矩 阵，其每列为模态振型。模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程 振动领域中的应用。模态是机械结构的固有振动特性，每一个模态具有特 定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析 取得，这样一个计算或试验分析过程称为模态分析。这个分析过程如果是 由有限元计算的

2、方法取得的，则称为计算模记分析；如果通过试验将采集 的系统输入与输出信号经过参数识别获得模态参数，称为试验模态分析。 通常，模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整 体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率 围各阶主要模态的特性，就可能预言结构在此频段在外部或部各种振源作 用下实际振动响应。因此，模态分析是结构动态设计及设备的故障诊断的 重要方法。模态分析最终目标是在识别出系统的模态参数，为结构系统的振动特 性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。 二、各模态分析方法的总结（一）单自由度法般来说，一个系统的动态响应是它的若干阶模态

3、振型的叠加。但是(SDOF)法2-12-2如果假定在给定的频带只有一个模态是重要的，那么该模态的参数可以单 独确定。以这个假定为根据的模态参数识别方法叫做单自由度 n1 o在给定的频带围，结构的动态特性的时域表达表示近似为:,ertQR T而频域表示则近似为:Qrh j j单自由度系统是一种很快速的方法，几乎不需要什么计算时间和计算 机存。这种单自由度的假定只有当系统的各阶模态能够很好解耦时才是正 确的。然而实际情况通常并不是这样的，所以就需要用包含若干模态的模 型对测得的数据进行近似，同时识别这些参数的模态，就是所谓的多自由 度（MDOF）法。单自由度算法运算速度很快，几乎不需要什么计算和计

4、算机存，因此 在当前小型二通道或四通道傅立叶分析仪中，都把这种方法做成置选项。 然而随着计算机的发展，存不断扩大，计算速度越来越快，在大多数实际 应用中，单自由度方法已经让位给更加复杂的多自由度方法。1、峰值检测峰值检测是一种单自由度方法，它是频域中的模态模型为根据对系统极点进行局部估计（固有频率和阻尼）。峰值检测方法基于这样的事实：在固有频率附近，频响函数通过自己的极值，此时其实部为零（同相部分最专业资料小），而虚部和幅值最大（相移达90° ,幅度达峰值）图1。出现极值的那个固有频率就是阻尼固有频率r的良好估计。相应的阻尼比,，的估计可用半功率点法得到。设 1和2分处在阻尼固有频率

5、的两侧（1< r< 2）,则:2-32、模态检测212""r2-4模态检测是根据频域中的模态模型对复模态（或实模态）向量进行局部估计的一种单自由度方法。在 h jQrrj rt-URLR中略去剩余项则单个频响函数在 处的值近似为:QrHtj j r-j rjrj rQr 1r jrA jr2-5由此式可见，频响函数在 处的值乘以模态阻尼因r,就是留数（的估计值如图1。利用这种模态检测方法之前，先要估计出图1对频响应函数的幅值进行峰值和模态检测3、圆拟合圆拟合是一种单自由度方法，用频域中的模态模型对系统极点和复模态（或实模态）向量进行局部估计。此方法依据事实是：单

6、自由度系统的速度频响函数（速度对力）在奈奎斯特图（即实部对虚部）上呈现为一个圆。如果把其他模态的影响近似为一个复常数,那么在共振频率r附近，频响函数2-5专业资料的基本公式为：Htj jU jVjrR j12-6因此，首先要选择共振频率附近的一组频率响应点，通过这些点拟合 成一个圆。阻尼固有频率可以看成是复平面上数据点之间角度变化率最 大（角间隔最大）的那个点的频率，也可以看成是相位角与圆心的相位角最 为接近的那个数据点的频率。 对于分得开的模态而言，二者的差别是很小。阻尼比,估计如下:2-72-821;1;rr tan 12tan 2 ?2：分居在r两侧的两个频率点:2：分别为频率点在1和2

7、得半径与r得半径之间的夹角。圆的直径和阻尼固有频率点的角位置含有复留数U+jV的信息:上上tan U rV式中：圆的直径:园心与固有频率点的连线跟虚轴之间的夹角.圆拟合法速度也很快，但为避免结果出错，特别是在模态节点附近, 需要操作者参与。（二）单自由度与多自由度系统粘性阻尼单自由度SDOF系统如图2的力平衡方程式表示惯性力、 阻尼力、弹性力与外力之间的平衡图2单自由度系统2-91,2 C/2M V C/2M 2 K/M2-11Mx t Cxt Kx t ft其中M:质量C:阻尼K: xxx:加速度，速度，位移f:外力t时间变 量，把结构中所呈现出来的全部阻尼都近似为一般的粘性阻尼。把上面的时

8、间域方程变换到拉氏域复变量P,并假设初始位移和初始速度为零，则得到拉氏域方程：Mp2 Cp K F p ,或Z pXZ:动刚度经过变换可得传递函数的定义，H p Z 1 p即X p2-10一1/MH p -p2 C/M p K/M上式右端的分母叫做系统特征方程，它的根即是系统的极点是:如果没有阻尼C=0 ,则所论系统是保守系统。我们定义系统的无阻尼 固有频率为:2-41 K/M临界阻尼Cc的定义为使（2.3）式中根式项等于零的阻尼值:Cc 2M . K/M c而临界阻尼分数或阻尼比（1为：匕=CCc,阻尼有时也有用品质因数即2-6Q因数表示:Q 1/ 2 i系统按阻尼值的大小可以分成过阻尼系统

9、（1>1 ）、临界阻尼系统（1 = 1 ）和欠阻尼系统（1<1 ）。过阻尼系统的响应只含有衰减成分、没有振荡趋势。欠阻尼系统的响应时一种衰减振动，而临界阻尼系统则是过阻尼系统与欠阻尼系统之间的一种分界。实际系统的阻尼比很少有大于 10%的，除非这些系统含有很强的阻尼 机制，因此我们只研究欠阻尼的情形。在欠阻尼的情况下式2-11两个共辗复根:.*.11 j 1 ,11 j 12-7其中1为阻尼因子1为阻尼固有频率。有关系统极点的另外一些关系式有:2-2式写成如下形式:在展开成部分分式形式，则有:1 j J 1211/M*H p A-,这里 A1p 1 p 11/Mj2 12-82-9

10、2-102-112-122-13这里的A和A*是留数。多自由度系统多自由度系统可以用简单的力平衡代数方程演化成形式相似的一个矩阵的方程。下面是以而自由度系统为例。如图:多自由度系统图3专业资料该系统的运动方程如下:M1x1C1C2x1tM 2乂2C2C3 x? tC?x2 tC2X tK1K2K2 x1 tK3也tK2x1 tf1 tf2 t2-14写成矩阵形式是M10x10M2 x2C2C2C2C2C3x1x2KiK2K2K2K2 K3x2或者M其中M、C、K、f(t)和x(t)分别为质量矩阵、fi2-152-16阻尼矩阵、刚度矩阵、方向量和响应向量。把这个时间域的矩阵方程变换到拉氏域(变量

11、为P)且假定初始位移和初始速度为零，则得:2_p M p C K X p F p2-172-18或者是 Z p X p F p 式中：Z(p)动刚度矩阵可以得到传递函数矩阵为:2-19i adj Z pH P Z PZ PT式中adj Z p : Z p的伴随矩阵，等于 j Zj ;Zij : Z p去掉第行第列后的行列式1 如果i j等于偶数1 如果ij等于奇数传递函数矩阵含有幅值函数。2-19式中的分母，即是 Z p的韩烈士，叫做系统的特征方程。与单自由度情况一样，系统特征方程的根，即系统极点，决定系统的共振频率。根据特征值问题，可以求出系统特征方恒的根。为了把系统方程2-17转 化为一般

12、的特征值问题公式，加入下面的恒等式:2-20将此式与2-17式结合在一起得:其中p ABYF'0MM0,BMC0K2-21如果力函数等于零，那么式2-19就成了关于实值矩阵的一般特征值问题，其特征值马祖下列方程的 p值：2-22它的根就是特征方程 Z p 0的根。对于N各自由度系统，此方程有2N个呈复共辗对出现的特征根:11 j 12-230 NN j N* 11 j 100* Nr是阻尼因子，虚同单自由度系统一样，多自由度系统的极点的实部 部r是阻尼固有频率。（三）实模态和复模态按照模态参数（主要指模态频率及模态向量）是实数还是复数，模态 可以分为实模态和复模态。对于无阻尼或比例阻尼

13、振动系统，其各点的振 动相位差为零或180度，其模态系数是实数，此时为实模态；对于非比例 阻尼振动系统，各点除了振幅不同外相位差也不一定为零或180度，这样模态系数就是复数，即形成复模态。1 复模态与实模态理论在拟合频段，实模态理论中传递函数在 k点激励Z点响应的留数表达式 为HkirRkl r 1 2., r2Vrejr arctan k, l 1,2, , n（1）Vr其中，rRkl为留数；和Vr构成的复数为系统的复特征值r： r r jVr拟合频段复模态理论中传递函数在k点激励f点响应的留数表达式为Hkl()nr12j r jVr2121rVr2ejr arctan -k, l 1, ,

14、 n（2）Vr由（1）、（2）式中可以看出，传递函数共振峰处复模态的相位与实模态相 位的差别在于多出的复留数相位r ,由传递函数的逆变换可以得到脉冲响应函数，由此可以得到物理坐标系中结构的自由响应表达式。对于无阻尼结构，t时刻第r阶模态k点的振动为xkr krYr sin rt r（ 3 ）粘性比例阻尼：t时刻第r阶模态k点的振动为xkr krYre rt sin drt r一般粘性阻尼：t时刻第r阶模态k点的振动为xkr 2Tr kre cos drt r kr（5）式中，k kr表不振型幅值；Q表示模态频率；0表示相位角。可以看出，无阻尼和比例阻尼系统的初相位与初始条件有关，与物 理坐标无

15、关，具有模态（振型）保持性； 而一般粘性阻尼系统的初相位还与物理坐标k有关，每个物理坐标振动时并不同时达到平衡位置和最 大位置， 不具备模态保持性，是行波形式.但各物理坐标的相位差保持不变，各点的振动周期、衰减率仍保持相同 J.从物理坐标点的自由响应公式还可看出，即使各测点留数为复数，但如果留数的相位差，即振型的幅角相同，那么还是可以得到振动周期形状不变且节点固定的振型.这样模态虽是复模态，但表现出实模态的性质.因此实模态理论的实振型与复模态理论中复模态的差别在于各测点峰值相位差的大小.2实模态提取方法复模态理论中模态参数（特征值和特征向量）均为复数，在进行结构 模型修正时大量采用复数矩阵和复

16、数迭代运算，计算工作量大，效率低； 实模态理论中模态参数为实数，物理概念明确，后续结构模型修正计算公 式简单，计算工作量小又节约空间，故实模态得到广泛的应用，实际测试 得到的传递函数留数一般都为复数，要由复模态经过实模态提取技术才能 得到实模态。复模态提取实模态的方法主要有：根据复模态的实部、虚部 或相位确定实模态的传统方法；I b r a h i m 的扩大模型法； C h e n的 传递函数提取法等。目前的模态分析软件中普遍使用的为传统方法。由复 模态实部或虚部获得实模态向量的方法为：直接取复留数的实部或虚部作 为实模态理论中的留数，进行规格化得到实模态振型.由复模态相位获得实模态向量的方

17、法为：取复留数的幅值作为实模态理论中的留数，根据sin r的数值接近1或-1 ,将留数相位归为 90 c或-90。，然后尽享振型规格化，得到实模态振型，此振型中各测点相位差 即为0 °或 180。用复模态理论获得的复模态向量，由复振型的周期变化中t=0即振动达到最大幅度时的振幅之比表示。三、模态分析的应用与发展模态分析技术的应用可归结为以下几个方面：1）评价现有结构系统的动态特性；2）在新产品设计中进行结构动态特性的预估和优化设计；3）诊断及预报结构系统的故障；4）控制结构的辐射噪声；5）识别结构系统的载荷。对于实际的工程，用有限元软件分析需要的频率段，可查找振动原因，或校核。模态分

18、析可以看出在那些频率段需要防止或避免共振时很有用。首先，频率和振型是结构的固有特性，任何结构都可以进行模态分析；其次，结构的功能是不同的，不同结构对应的模态分析的用途是有差别的。 对建筑结构，模态分析可以知道结构的避频设计、用于抗震设计计算以及 考虑动力荷载的放大作用等。另外，还可以挖掘振型有关的信息。机器、建筑物、航天航空飞行器、船舶、汽车等的实际振动千姿百态、 瞬息变化。模态分析提供了研究各种实际结构振动的一条有效途径。首先，将结构物在静止状态下进行人为激振，通过测量激振力与胯动 响应并进行双通道快速傅里叶变换(FFT)分析，得到任意两点之间的机械导纳函数(传递函数)。用模态分析理论通过对

19、试验导纳函数的曲线拟 合，识别出结构物的模态参数，从而建立起结构物的模态模型。根据模态 叠加原理，在已知各种载荷时间历程的情况下，就可以预言结构物的实际 振动的响应历程或响应谱。模态分析软件以美国的ME' ScopeVES的功能最为全面。ME'ScopeVES软件的功能包括信号处理(signal Process ing)、运行挠曲振型 (OperatingDeflection Shapes) 、模态分析(Modal Analysi s)、结构改正 (SDM)和声学分析(Acoustics Analysi S) 等，解决和分析机器与结构的振 动噪声问题。主要可用于：1、可以显示被

20、测物体的实际工作形态 (0DS)、模态、声学分布形态和 工程数据的形态等；2、模块化结构便于用户根据自己的需要选择合适的产品；3、强大的图形显示、结构编辑、数据处理及动画显示功能；4、软件的开放性好，能够与全球的多家厂商的硬件兼容；5、主要应用的领域：航空航天、建筑桥梁、汽车制造、钢铁冶金、 军工装备等。模态分析与参数辨识作为结构动力学中的一种逆问题分析方法并在 工程实践中应用是从60年代中、后期开始，至今已有近四十年的历史了。 这一技术首先在航空、宇航及汽车工业中开始发展。由于电子技术、信号 处理技术与设备的发展，到 80年代末这项技术已成为工程中解决结构动 态性能分析、振动与噪声控制、故障

21、诊断等问题的重要工具。目前这一技 术已渐趋成熟。经过二十余年的研究发展，到目前为止模态分析技术已在 我国各个工程领域中广泛应用，成为一种解决工程问题的重要手段。在工程应用方面模态分析已渗透到我国各个工程领域，并取得了不少 成就。例如，某型火箭全装置的实物模态试验保证了火箭的准确发射与导 航，防止了发射的失败；模态分析与参数识别技术曾被成功地用于解决某 型航空发动机的严重振动故障，取得重大经济及社会效益；某型鱼雷全装 置实物水下模态试验为鱼雷的振动与噪声控制确保导航性能提供了技术 依据；远东第一高塔的上方明珠电视塔的振动模态试验，为高塔的抗风抗 地震安全性设计提供了技术依据；目前世界上跨度第一的斜拉索浦大桥的 振动试验对大桥抗风振动的安全性分析与故障诊断提供了技术依据；建立 在模态