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

.各种模态分析方法总结与比较一、模态分析模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参数的过程。模态分析的理论经典定义：将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵，其每列为模态振型。模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的，则称为计算模记分析；如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数，称为试验模态分析。通常，模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性，就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此，模态分析是结构动态设计及设备的故障诊断的重要方法。模态分析最终目标是在识别出系统的模态参数，为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。二、各模态分析方法的总结.（一）单自由度法一般来说，一个系统的动态响应是它的若干阶模态振型的叠加。但是如果假定在给定的频带内只有一个模态是重要的，那么该模态的参数可以单独确定。以这个假定为根据的模态参数识别方法叫做单自由度(SDOF)法n1。在给定的频带范围内，结构的动态特性的时域表达表示近似为：h tr e rtT2-1QRR而频域表示则近似为：QrtLRrrUR2-2h j2jr单自由度系统是一种很快速的方法，几乎不需要什么计算时间和计算机内存。这种单自由度的假定只有当系统的各阶模态能够很好解耦时才是正确的。然而实际情况通常并不是这样的，所以就需要用包含若干模态的模型对测得的数据进行近似，同时识别这些参数的模态，就是所谓的多自由度 (MDOF)法。单自由度算法运算速度很快，几乎不需要什么计算和计算机内存，因此在当前小型二通道或四通道傅立叶分析仪中，都把这种方法做成内置选项。然而随着计算机的发展，内存不断扩大，计算速度越来越快，在大多数实际应用中，单自由度方法已经让位给更加复杂的多自由度方法。1、峰值检测峰值检测是一种单自由度方法，它是频域中的模态模型为根据对系统极点进行局部估计( 固有频率和阻尼 ) 。峰值检测方法基于这样的事实：在固有频率附近，频响函数通过自己的极值，此时其实部为零( 同相部分最.小 ) ，而虚部和幅值最大( 相移达 90，幅度达峰值 ) 图 1。出现极值的那个固有频率就是阻尼固有频率r 的良好估计。 相应的阻尼比r ，的估计可用半功率点法得到。设1 和2 分处在阻尼固有频率的两侧 (1 < r < 2 ),则：H J 1H j 2H j r2-32212-4r2r2、模态检测模态检测是根据频域中的模态模型对复模态( 或实模态 ) 向量进行局QrtLR2 中略去剩余部估计的一种单自由度方法。在h jrrURjr项则单个频响函数在r 处的值近似为：H tj j rQr1rjrQr 1rjrA1 jr2-5jj rrrrr由此式可见，频响函数在r 处的值乘以模态阻尼因r ，就是留数 ( 的估计值如图 1。利用这种模态检测方法之前，先要估计出r图1 对频响应函数的幅值进行峰值和模态检测3、圆拟合圆拟合是一种单自由度方法，用频域中的模态模型对系统极点和复模.态 ( 或实模态 ) 向量进行局部估计。此方法依据事实是：单自由度系统的速度频响函数 ( 速度对力 ) 在奈奎斯特图 ( 即实部对虚部 ) 上呈现为一个圆。如果把其他模态的影响近似为一个复常数，那么在共振频率r 附近，频响函数的基本公式为：H tj jUjV2-6jR j1r因此，首先要选择共振频率附近的一组频率响应点，通过这些点拟合成一个圆。 阻尼固有频率r 可以看成是复平面上数据点之间角度变化率最大 ( 角间隔最大 ) 的那个点的频率，也可以看成是相位角与圆心的相位角最为接近的那个数据点的频率。 对于分得开的模态而言， 二者的差别是很小。阻尼比 r 估计如下：r212-7rtan1tan222式中1 ，2 ：分居在r 两侧的两个频率点：1， 2 ：分别为频率点在1 和2 得半径与r 得半径之间的夹角。圆的直径和阻尼固有频率点的角位置含有复留数U+jV的信息：U 2V 2U, tan2-8rV式中：圆的直径：园心与固有频率点的连线跟虚轴之间的夹角圆拟合法速度也很快，但为避免结果出错，特别是在模态节点附近，需要操作者参与。.（二）单自由度与多自由度系统粘性阻尼单自由度SDOF系统如图 2的力平衡方程式表示惯性力、阻尼力、弹性力与外力之间的平衡图 2 单自由度系统M x tCx tKx tf t2-9其中 M：质量 C: 阻尼 K： xxx ：加速度，速度，位移f ：外力 t 时间变量，把结构中所呈现出来的全部阻尼都近似为一般的粘性阻尼。把上面的时间域方程变换到拉氏域复变量P，并假设初始位移和初始速度为零，则得到拉氏域方程：Mp 2CpKF p ，或 Z p Xp F pZ：动刚度经过变换可得传递函数的定义，H pZ 1 p 即 X pH p F pH p1/ M2-102C / MpK / Mp上式右端的分母叫做系统特征方程，它的根即是系统的极点是：1, 2C / 2MC / 2M2K / M2-11如果没有阻尼 C=0，则所论系统是保守系统。我们定义系统的无阻尼固有频率为：1K / M2-4临界阻尼 Cc 的定义为使（ 2.3）式中根式项等于零的阻尼值：.Cc 2M K / M2-5而临界阻尼分数或阻尼比1 为： 1=CCc，阻尼有时也有用品质因数即 Q 因数表示：Q1/ 2 12-6系统按阻尼值的大小可以分成过阻尼系统（ 1 >1）、临界阻尼系统（1=1）和欠阻尼系统（1<1）。过阻尼系统的响应只含有衰减成分、没有振荡趋势。欠阻尼系统的响应时一种衰减振动，而临界阻尼系统则是过阻尼系统与欠阻尼系统之间的一种分界。实际系统的阻尼比很少有大于 10%的，除非这些系统含有很强的阻尼机制，因此我们只研究欠阻尼的情形。在欠阻尼的情况下式2-11 两个共轭复根：11j1 ， 1*1j12-7其中1 为阻尼因子1 为阻尼固有频率。有关系统极点的另外一些关系式有：j 122-81111112-9221111 12-10222-111112-2 式写成 如下形式：H p1/ M2-12p*p11在展开成部分分式形式，则有：.A1A1*1/ M2-13H p1p* ，这里 A1j2 1p1这里的 A1 和 A1* 是留数。多自由度系统多自由度系统可以用简单的力平衡代数方程演化成形式相似的一个矩阵的方程。下面是以而自由度系统为例。如图：图 3多自由度系统该系统的运动方程如下：M 1 x1C1C2 x1 t C2 x2 tK1K 2 x1 t K 2 x2 tf1tM 2 x2C2C3 x2 t C2 x1 tK 2K 3 x2 t K 2 x1 t2-14f2 t写成矩阵形式是M 1 0 x1C1 C2C2x1K1 K 2K 2x1f12-150 M 2 x2C2C2 C3x2K 2K 2 K 3x2f2或者 M x C xK xf2-16其中 M、C、 K、 f(t) 和 x(t) 分别为质量矩阵、阻尼矩阵、刚度矩阵、方向量和响应向量。把这个时间域的矩阵方程变换到拉氏域（变量为p）且假定初始位移和初始速度为零，则得：p2 Mp CKX pF p2-17或者是Z pX pF p式中： Z( p) 动刚度矩阵2-18.可以得到传递函数矩阵为：1adj Z p2-19H pZ pZ p式中adj Z p ： Z p 的伴随矩阵，等于Tij Zij；Zij ： Z p去掉第行第列后的行列式ij1如果 ij等于偶数；1如果 ij等于奇数传递函数矩阵含有幅值函数。2-19 式中的分母，即是Z p的韩烈士，叫做系统的特征方程。与单自由度情况一样，系统特征方程的根，即系统极点，决定系统的共振频率。根据特征值问题， 可以求出系统特征方恒的根。为了把系统方程2-17 转化为一般的特征值问题公式，加入下面的恒等式：p M p M X02-20将此式与 2-17 式结合在一起得：p ABYF 2-21其中A0M，BM0，MC0KYp X，F 0。XF如果力函数等于零，那么式2-19 就成了关于实值矩阵的一般特征值问题，其特征值马祖下列方程的p 值：p AB02-22它的根就是特征方程Z p0 的根。对于N 各自由度系统，此方程有2N 个呈复共轭对出现的特征根：.11j 1002-23NNj N*1j1100*Nj NN同单自由度系统一样，多自由度系统的极点的实部r 是阻尼因子，虚部 r 是阻尼固有频率。（三）实模态和复模态按照模态参数（主要指模态频率及模态向量）是实数还是复数，模态可以分为实模态和复模态。对于无阻尼或比例阻尼振动系统，其各点的振动相位差为零或180 度，其模态系数是实数，此时为实模态；对于非比例阻尼振动系统，各点除了振幅不同外相位差也不一定为零或 180 度，这样模态系数就是复数，即形成复模态。1复模态与实模态理论在拟合频段， 实模态理论中传递函数在k 点激励 Z 点响应的留数表达式为.r Rkle j rH kl2r 12vr2rrarctanrk, l1,2, , n(1)vr其中， r Rkl 为留数； r和 vr 构成的复数为系统的复特征值r ： rrjv r拟合频段复模态理论中传递函数在k 点激励 f点响应的留数表达式为H kl (nr Rkl e j rnr Rklej rj)2err 1 2 jrjvrr 1 2vr2rrarctanrvrk， l1， ， n（2）.由（ 1）、（ 2）式中可以看出，传递函数共振峰处复模态的相位与实模态相位的差别在于多出的复留数相位r ，由传递函数的逆变换可以得到脉冲响应函数，由此可以得到物理坐标系中结构的自由响应表达式。对于无阻尼结构， t 时刻第 r 阶模态 k 点的振动为xkrkr Yr sinr tr（3）粘性比例阻尼： t 时刻第 r 阶模态 k 点的振动为xkrkrYr er t sindr tr（4）一般粘性阻尼： t 时刻第 r 阶模态 k 点的振动为xkr2Tr kr emr t cosdr trkr(5)式中， kr 表示振型幅值；表示模态频率；表示相位角。可以看出，无阻尼和比例阻尼系统的初相位与初始条件有关，与物理坐标无关，具有模态 (振型 )保持性；而一般粘性阻尼系统的初相位还与物理坐标k有关，每个物理坐标振动时并不同时达到平衡位置和最大位置，不具备模态保持性，是行波形式但各物理坐标的相位差保持不变，各点的振动周期、衰减率仍保持相同J从物理坐标点的自由响应公式还可看出，即使各测点留数为复数，但如果留数的相位差，即振型的幅角相同，那么还是可以得到振动周期内形状不变且节点固定的振型这样模态虽是复模态，但表现出实模态的性质因此实模态理论的实振型与复模态理论中复模态的差别在于各测点峰值相位差的大小2 实模态提取方法复模态理论中模态参数(特征值和特征向量) 均为复数，在进行结构模型修正时大量采用复数矩阵和复数迭代运算，计算工作量大，效率低；.实模态理论中模态参数为实数，物理概念明确，后续结构模型修正计算公式简单，计算工作量小又节约空间，故实模态得到广泛的应用，实际测试得到的传递函数留数一般都为复数，要由复模态经过实模态提取技术才能得到实模态。复模态提取实模态的方法主要有：根据复模态的实部、虚部或相位确定实模态的传统方法；I b r a h i m的 扩大模型法；C h e n的传递函数提取法等。目前的模态分析软件中普遍使用的为传统方法。由复模态实部或虚部获得实模态向量的方法为：直接取复留数的实部或虚部作为实模态理论中的留数，进行规格化得到实模态振型由复模态相位获得实模态向量的方法为：取复留数的幅值作为实模态理论中的留数，根据 sinr的数值接近1 或-1 ，将留数相位归为90 或 - 90 ，然后尽享振型规格化，得到实模态振型，此振型中各测点相位差即为 0 或 180 。用复模态理论获得的复模态向量，由复振型的周期变化中 t=0 即振动达到最大幅度时的振幅之比表示。三、模态分析的应用与发展模态分析技术的应用可归结为以下几个方面：1) 评价现有结构系统的动态特性；2) 在新产品设计中进行结构动态特性的预估和优化设计；3) 诊断及预报结构系统的故障；4) 控制结构的辐射噪声；5) 识别结构系统的载荷。对于实际的工程， 用有限元软件分析需要的频率段，可查找振动原因，或校核。模态分析可以看出在那些频率段需要防止或避免共振时很有用。.首先，频率和振型是结构的固有特性，任何结构都可以进行模态分析；其次，结构的功能是不同的， 不同结构对应的模态分析的用途是有差别的。对建筑结构，模态分析可以知道结构的避频设计、用于抗震设计计算以及考虑动力荷载的放大作用等。另外，还可以挖掘振型有关的信息。机器、建筑物、 航天航空飞行器、 船舶、汽车等的实际振动千姿百态、瞬息变化。模态分析提供了研究各种实际结构振动的一条有效途径。首先，将结构物在静止状态下进行人为激振，通过测量激振力与胯动响应并进行双通道快速傅里叶变换（FFT ）分析，得到任意两点之间的机械导纳函数（传递函数） 。用模态分析理论通过对试验导纳函数的曲线拟合，识别出结构物的模态参数，从而建立起结构物的模态模型。根据模态叠加原理，在已知各种载荷时间历程的情况下，就可以预言结构物的实际振动的响应历程或响应谱。模态分析软件以美国的MEScopeVES的功能最为全面。 MEScopeVES软 件 的 功 能 包 括 信 号 处 理 (signalProcessing) 、 运 行 挠 曲 振 型(OperatingDeflection Shapes)、模态分析 (Modal Analysi s)、结构改正 (SDM)和声学分析 (Acoustics Analysi S)等，解决和分析机器与结构的振动噪声问题。主要可用于：1、可以显示被测物体的实际工作形态(0DS)、模态、声学分布形态和工程数据的形态等；2、模块化结构便于用户根据自己的需要选择合适的产品；3、强大的图形显示、结构编辑、数据处理及动画显示功能；.4、 件的开放性好，能 与全球的多家厂商的硬件兼容；5、主要 用的 域：航空航天、建筑 梁、汽 制造、 冶金、 工装 等。模 分析与参数辨 作 构 力学中的一种逆 分析方法并在工程 践中 用是从60 年代中、后期开始，至今已有近四十年的 史了。 一技 首先在航空、宇航及汽 工 中开始 展。由于 子技 、信号 理技 与 的 展，到 80 年代末 技 已成 工程中解决 构 性能分析、振 与噪声控制、故障 断等 的重要工具。目前 一技 已 成熟。 二十余年的研究 展，到目前 止模 分析技 已在我国各个工程 域中广泛 用，成 一种解决工程 的重要手段。在工程 用方面模 分析已渗透到我国各个工程 域，并取得了不少成就。例如，某型火箭全装置的 物模 保 了火箭的准确 射与 航，防止了 射的失 ；模 分析与参数 技 曾被成功地用于解决某型航空 机的 重振 故障，取得重大 及社会效益；某型 雷全装置 物水下模 雷的振 与噪声控制确保 航性能提供了技 依据； 第一高塔的上海 方明珠 塔的振 模 ， 高塔的抗 抗地震安全性 提供了技 依据；目前世界上跨度第一的斜拉索 浦大 的振 大 抗 振 的安全性分析与故障 断提供了技 依据；建立在模 分析技 上的 基断裂 技 已在高 建筑施工中广泛 用，提高了 基的 量，确保高 建筑的安全； 等等， 些成就不 枚 。 之，二十余年的 展是迅速的，成就是 著的，回 一 展 程和取得的成就，可更激励我 朝着新的目 前 。.模 分析技 展到今天已 成熟，特 是 性模 理 方面的研究已日臻完善，但在工程 用方面 有不少工作可做。首先是如何提高模 分析的精度， 大 用范 。增加模 分析的信息量是提高分析精度的关 ， 靠增加 感器的 点数目很 ，目前提出的一种激光 描方法是大大增加 点数的有效 法， 点数目的增加随之而来的是增大数据采集与分析系 的容量及提高分析 理速度，在 方法、数据采集与分析方面 有不少研究工作可做。 复 构空 模 的 量分析、 响函数的耦合、高 模 、 抗噪声干 等等方面的研究尚需 一步开展。模 分析当前的一个重要 展 是由 性向非 性 方向 展。非 性模 的概念早在 1960 年就由 Rosenberg提出， 有不少学者 非 性模 理 行了研究，但由于非 性 本身的复 性及当 工程 践中的非 性 并引引起重 ，非 性模 分析的 展受到限制。近年来在工程中的非 性 日益突出，因此非 性模 分析亦日益受到人 的重 。最近已逐步形成了所 非 性模 力学。关于非 性模 的正交性、解耦性、 定性、模 的分叉、渗透等 是当前研究的重点。在非 性建模理 与参数辨 方面的研究工作亦是当今研究的 点。非 性系 物理参数的 、 荷 方面的研究亦已开始。展望未来，模 分析与 技 仍将以新的速度，新的内容向前 展。姓名：徐海 学号：.