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

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

一、模态分析

模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参数的过程。

模态分析的理论经典定义：将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵，其每列为模态振型。

模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的，则称为计算模记分析；如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数，称为试验模态分析。通常，模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性，就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此，模态分析是结构动态设计及设备的故障诊断的重要方法。

模态分析最终目标是在识别出系统的模态参数，为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。

二、各模态分析方法的总结 （一）单自由度法

一般来说，一个系统的动态响应是它的若干阶模态振型的叠加。但是如果假定在给定的频带内只有一个模态是重要的，那么该模态的参数可以单独确定。以这个假定为根据的模态参数识别方法叫做单自由度(SDOF)法n1。在给定的频带范围内，结构的动态特性的时域表达表示近似为：

()[]}{}{T R R t r Q e t h r

ψψλ= 2-1

而频域表示则近似为：

()[]}}{

{()[]2

ω

λωψψωLR UR j Q j h r t

r r r -+-= 2-2 单自由度系统是一种很快速的方法，几乎不需要什么计算时间和计算机内存。

这种单自由度的假定只有当系统的各阶模态能够很好解耦时才是正确的。然而实际情况通常并不是这样的，所以就需要用包含若干模态的模型对测得的数据进行近似，同时识别这些参数的模态，就是所谓的多自由度(MDOF)法。

单自由度算法运算速度很快，几乎不需要什么计算和计算机内存，因此在当前小型二通道或四通道傅立叶分析仪中，都把这种方法做成内置选项。然而随着计算机的发展，内存不断扩大，计算速度越来越快，在大多数实际应用中，单自由度方法已经让位给更加复杂的多自由度方法。 1、峰值检测

峰值检测是一种单自由度方法，它是频域中的模态模型为根据对系统极点进行局部估计(固有频率和阻尼)。峰值检测方法基于这样的事实：在固有频率附近，频响函数通过自己的极值，此时其实部为零(同相部分最小)，而虚部和幅值最大(相移达90°，幅度达峰值)图1。出现极值的那个固有频率就是阻尼固有频率r ω的良好估计。相应的阻尼比r ζ，的估计可用半功率点法得到。设1ω和2ω分处在阻尼固有频率的两侧(1ω

()()()2

21r j H j H J H ωωω=

= 2-3

r

r ωωωζ21

2-=

2-4 2、模态检测

模态检测是根据频域中的模态模型对复模态(或实模态)向量进行局

部估计的一种单自由度方法。在()[]}}{

{()[]2ω

λωψψωLR

UR j Q j h r t

r r r -+-=中略去剩余

项则单个频响函数在r ω处的值近似为：

()()()

r

jr r

jr

r r r r r jr

r r r tj A Q j j Q j H σσψψωσωψψω-≈

-≈

+-≈

111 2-5

由此式可见，频响函数在r ω处的值乘以模态阻尼因r σ，就是留数(的估计值如图1。利用这种模态检测方法之前，先要估计出r ω

图1 对频响应函数的幅值进行峰值和模态检测

3、圆拟合

圆拟合是一种单自由度方法，用频域中的模态模型对系统极点和复模态(或实模态)向量进行局部估计。此方法依据事实是：单自由度系统的速度频响函数(速度对力)在奈奎斯特图(即实部对虚部)上呈现为一个圆。如果把其他模态的影响近似为一个复常数，那么在共振频率r ω附近，频响函数的基本公式为：

()()

1j R j jV

U j H r tj ++-+-+=

ωωσω 2-6

因此，首先要选择共振频率附近的一组频率响应点，通过这些点拟合成一个圆。阻尼固有频率r ω可以看成是复平面上数据点之间角度变化率最大(角间隔最大)的那个点的频率，也可以看成是相位角与圆心的相位角最为接近的那个数据点的频率。对于分得开的模态而言，二者的差别是很小。 阻尼比r ζ估计如下：

⎪⎭

⎫ ⎝

⎛+⎪⎭

⎫ ⎝

⎛⎪⎭⎫ ⎝

⎛-=

2tan 2tan 211

2θθωωωζr r 2-7

式中1ω，2ω：分居在r ω两侧的两个频率点：

1θ，2θ：分别为频率点在1ω和2ω得半径与r ω得半径之间的夹角。

圆的直径和阻尼固有频率点的角位置含有复留数U+jV 的信息：

()V

U

V U r =-+=ασφtan ,22 2-8

式中φ：圆的直径

α：园心与固有频率点的连线跟虚轴之间的夹角．

圆拟合法速度也很快，但为避免结果出错，特别是在模态节点附近，需要操作者参与。

（二）单自由度与多自由度系统

粘性阻尼单自由度SDOF 系统如图2的力平衡方程式表示惯性力、阻尼力、弹性力与外力之间的平衡

图2 单自由度系统

()()()()t f t Kx t x C t x

M =++&&& 2-9 其中M ：质量C: 阻尼K ：x

x x &&&&&：加速度，速度，位移 f ：外力 t 时间变量，把结构中所呈现出来的全部阻尼都近似为一般的粘性阻尼。

把上面的时间域方程变换到拉氏域复变量P ，并假设初始位移和初始速度为零，则得到拉氏域方程：()()p F K Cp Mp =++2，或()()()p F p X p Z = Z ：动刚度经过变换可得传递函数的定义，()()p Z p H 1-= 即()()()p F p H p X =

()()()

M K p M C p M

p H ///12

++=

2-10 上式右端的分母叫做系统特征方程，它的根即是系统的极点是：