模态分析与参数识别.docx

1、模态分析与参数识别模态分析方法在发动机曲轴上的应用研究XX（XX 大学 xxxxxxxx 学院，lL1 西 太原 030051 ）摘要：综述模态分析在研究结构动力特性中的应用，介绍模态分析的两大方法：数 值模态分析与试验模态分析。并着重介绍目前的研究热点一一工作模态分析。通过发动机曲轴的模态分析这一具体的实例，综述了运行模态分析国内外研究现状，指出了其尖键技术、存在问题以及研究发展方向。尖键词：模态分析数值模态试验模态工作模态Abstract : Sums up methods of model analysis applied on the research of con figurati

2、on dyn amic;al characteristio. It in troduces two methods of model an alysis: nu merical value model an alysis and experime ntati on model an alysis. The n it stresses the hotspot-working model analysis. Some key techniques, unsolved problems and research direct ions of OMA were also discussed.Key w

3、ords: Model analysis Numerical value model analysis Experimentation model analysis Working model analysis1 -1模态分析的基本概念物体按照某一阶固有频率振动时，物体上各个点偏离平衡位置的位移是满足一定 的比例尖系的，可以用一个向量表示，这个就称之为模态。模态这个概念一般是在振 动领域所用，你可以初步的理解为振动状态，我们都知道每个物体都具 有自己的固有 频率，在外力的激励作用下，物体会表现出不同的振动特性。一阶模态是外力的激励频率与物体固有频率相等的时候出现的，此时物体的 振动形态叫做一阶

4、振型或主振型；二阶模态是外力的激励频率是物体固有频率的两倍 时候出现，此时的振动外形叫做二阶振型，以依次类推。一般来讲，外界激励的频率非常复杂，物体在这种复杂的外界激励下的振动反应 是各 阶振型的复合。模态是结构的固有振动特性，每一个模态具有特定的固有频率、阻尼 比和模态振型。这些模态参数可以由计算或试验分析取得， 这样一个计算或试验分析过程称为模态分析。模态分析经典定义：将线性定常系统振动微分方程组中的物理坐标变换为模 态坐 标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程， 以便求出系统的模态参数。模态分析方法主要分三类，分别是试验模态分析EMA工作模态分析OMA和 工 作变形

5、分析ODS(1 )试验模态分析(Experimental Modal Analysis，EMA，也称为传统模态 分析或经典模态分析，是指通过输入装置对结构进行激励，在激励的同时测量结构的 响应的一种测试分析方法。输入装置主要有力锤和激振器，因此，实验模态分析又 分为力锤激励EMA技术和激振器激励EMA技术。(2)工作模态分析(Operational Modal Analysis 、OMA，也称为只有输 出 的模态分析，而在土木桥梁行业，工作模态分析又称为环境激励模态分析。这类分 析最明显的特征是对测量结构的输出响应，不需要或者无法测量输入。当受传感器 数量和采集仪通道数限制时，需要分批次进行测

6、量。(3)工作变形分析(Operational Deflection Shape, ODS，也称为运行 响应 模态。这类分析方法也只测量响应，不需要测量输入。但是它跟OMA勺区别在于， OMA寻到的是结构的模态振型，而ODS得到的是结构在某一工作状态下的变形形 式。此时分析出来的ODS振型已不是我们常说的模态振型了，它实际是结构模态振 型按某线性方式叠加的结果。模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程振动领域中 的应用。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比 和模态振型。这些模态参数可以由计算或试验分析取得， 这样一个计算或试验分析过程称为模态分

7、析。这个分析过程如果是由有限元计算的方法取得的，则称 为计算模态分析；如果通过试验将采集的系统输入与输出信号经过参数识别 获得模态 参数，称为试验模态分析。振动模态是弹性结构固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内，各阶主要 模态的特性，就可能预知结构在此频段内，在外部或内部各种振源作用下实际振动响 应。因此，模态分析是结构动态设计及设备故障诊断的重要方法。近十余年以来，模态分析的理论基础，已经由传统的线性位移实模态、复模 态 理论发展到广义模态理论，并被进一步引入到非线性结构振动分析领域，同时模态 分析理论汲取了振动理论、信号分析、数据处理、数理统计以

8、及自动控制的相尖理 论，结合自身的发展规律，形成了一套独特的理论体系，创造了更加广泛的应用前 景。这一技术已经在航空、航天、造船、 机械、建筑、交通运输和兵器等工程领域得到广泛应用。1.2数值模态分析与试验模态分析现状及局限性模态分析过程如果是由有限元计算的方法取得的，则称为数值模态分析 ；如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数，称为试验 模态分析。两种方法各有利弊，目前的发展趋势是把有限元方法和试验模态分析技术 有机地结合起来，取长补短，相得益彰。利用试验模态分析结果检验、补充和修 正原始有限元动力模型；利用修正后的有限元模型计算结构的动力特性和响应，进行 结构的优化

9、设计。数值模态分析主要采用有限元法，它是将弹性结构离散化为有限数量的具体质 量、弹性特性单元后，在计算机上作数学运算的理论计算方法。它的优点是可以在结 构设计之初，根据有限元分析结果，便预知产品的动态性能，可以在产 品试制出来之 前预估振动、噪声的强度和其他动态问题，并可改变结构形状以消除或抑制这些问 题。只要能够正确显示出包含边界条件在内的机械振动模型， 就可以通过计算机改变机械尺寸的形状细节。有限元的不足是计算繁杂，耗资费时5。这种方法，除要求计算者有熟练的技巧与经验外，有些参数如阻尼、结合 面特 征等）目前尚无法定值，并且利用有限元法计算得到的结果，只能是一个近 似值 L6】o试验模态分

10、析是模态分析中最常用的，它与有限元分析技术一起成为解决现代复杂结构动力学问题的两大支柱。利用试验模态分析研究系统动态性能是一种更经济、更有实效的方法。首先，根据已有的知识和经验，在老产品基础上试制 出一 台新的模型；其次，用试验模态分析技术，对样机作全面的测试与分析，获 得产品的 动力特性，由此识别出系统的模态参数，建立数学模型，进而了解产品在实际使用中 的振动、噪声、疲劳等现实问题；再次，在计算机上改变产品的结 构参数，了解动 态性能可能获得的改善程度，或者反过来，设计者事先指定好动力特性，由计算机 来回答所需要的结构参数质量、刚度、阻尼）的改变量 （7）。传统的试验模态分析方法是建立在系统

11、输入输出数据均已知的基础上，利用 激励和响应的完整信息进行参数识别。将结构物在静止状态下进行人为激振，通过 测量激振力与响应并进行双通道快速傅里叶变换分析，得到任意两点之间的机械导纳 函数即传递函数。用模态分析理论通过对试验导纳函数的曲线拟合，识别出结构物的 模态参数，从而建立起结构物的模态模型。根据模态叠加原理，在已知各种载荷 时间历程的情况下，就可以预知结构物的实际振 动的响应历程或响应谱。在试验模态分析中大致可以分为四个步骤：1 ）动态数据的采集及频响函 数或脉冲响应函数分析；2 ）建立结构数学模型;3 ）参数识别；4）振形动画。与有 限元方法相比，结构动力修改的问题在试验模态分析基础上

12、要容易。传统的模态分析方法已经在桥梁、汽车和航空航天工程等几乎所有和结构动态分 析有尖的领域中得到广泛应用，数值模态分析与试验模态分析的方法在 理论上已经趋于完善，然而这些方法在具体应用时还是存在局限性，因为对于某些实 际工程结构，要获得输入激励的完整信息是难以实现的，或者根本就没有 获得任何输入信息，具体表现为：1）海洋平台、建筑物以及桥梁等在风、浪以及大地脉动作用下引起的振动；导 弹以及航天器在飞行运输过程中所产生的振动等，这些结构在实际工作时所承受的载 荷往往是不可测量或很难测得，因而无法获得结构系统的激励输入信息；2 ）某些结构待识别的自由度很多，并且所受载荷的空间分布复杂，往往没有足

13、够的传感器，无法得到完整的输入信息；3）所需要的载荷测试量力）与能够测试的加速度不是同一类信号，所需 要的量不能直接测试，不能满足识别方法的要求；4)很多实际工作中，例如武器结构的振动试验，已经得到大量的振动响应数 据，但却没有输入数据。然而目前根据实测振动响应数据往往只能得到诸如共振频 率、最大峰值、总均方根值等特征量，而不能进一步用于对产品结构的动力特性分 析，这就难以对产品的整体变形、响应特性作完整了解，大量的试验结构不能得到充 分利用。针对传统的试验模态分析方法的局限性，发展仅基于响应数据的工作模态分析技 术显得尤其重要。采用工作模态分析技术可以避免对输入信息的采集， 这样也就解决了传

14、统分析方法中很多状况下输入不可测的问题(8)o1.3工作模态分析现状及发展趋势工作模态分析常称为环境激励下的模态分析、只有输出或激励未知条件下的模 态分析，正是近年来模态分析领域发展活跃，新理论、新技术的应用层出不穷的一 个研究方向，被视为对传统试验模态分析方法的创新和扩展。 工作模态分析的优点是：仅需测试振动响应数据，由于这些数据直接来源于结构实际所经受的振动 工作环境，因而识别结果更符合实际情况和边界条件；无需对输入激励进行测试，节 省了测试费用；利用实时响应数据进行模态参数识别，其结果能够直接 应用于结构的 在线健康监测和损伤诊断何0因此工作模态试验技术使试验模态分析，由传统的主要 针对

15、静止结构被扩展到处于现场运行状态的结构， 不仅可以实现对那些无法测得载荷的工程结构进行所谓在线模态分析，而且利用实际工作状 态下的响应数据识别的模态参数，能更加准确地反映结构的实际动态特性， 已经在桥梁、建筑、机械领域取得实质性的进展。工作模态分析的理论和思想的提出早在20世纪70年代初期就已开始。工作 模 态的主要手段都是基于响应信号的时域参数辨识技术。 随机减量技术最早被用来处理环境激励下的结构响应数据，这一技术 主要是将结构的随机响应转化为结构的自由响应Mo以此为基础基于时域的辨识方法Ibrahim时域法被提出，极 大推 动了工作模态分析技术的发展。随着控制理论和计算机技术的发展，多输入

16、、多输 出、参数辨识技术也被相继推出，广泛运用的时域模态辨 识方法有多参考 点复指数方法、特征系统实现算法等。目前工作模态辨识的其他主要方法还有功 率谱密度函数的峰值提取方法、建立自回归滑动平均模型的时间序列分析法、结 合时 域参数识别的随机减量技术等I。1965年Clarkson和Mercer提出使用互相尖函数估计承受白噪声激励下结 构的 频响特性，从而提出了当激励未知时使用相尖函数替代脉冲响应函数的思想框架。20 世纪90年代以来，美国Sandia国家实验 室结合时域模态辨识方法，提出了 NExT 技术，利用结构在环境激励下的响应的相尖函数进行工作模态识别形成上述技术思路 后，美国Sand

17、ia国家实验室已经将此分析成果成功运用于航天涡轮机、地面载重、 高速公路大桥和濒海建筑的工况信号测量和结构分析中13O在国内，南京航空航天大学振动工程研究所也一直从事着模态分析的研究工 作， 从传统的模态分析到工作模态分析，也包括只利用响应数据进行系统模态参 数识别方 法的研究，并且发表了多篇矣于环境激励下工作模态参数识别的文章。中国振动协 会，上海交通大学振动、冲击、噪声国家重点实验室，哈尔滨工业大 学等也致力于研 究工作模态参数识别方法2。现有和各种工作模态参数识别方法虽然都有一些很好的应用，但在理论上还 需要完善。而且各种工作模态分析方法还有着各自的局限性，如时域法通常要求 激励是平稳白

18、噪声，结构系统具有线性时不变特性，其中Thrihim法不易剔除噪声和 虚假模态；而时间序列法的模型阶次较难确定；基于响应相尖函数的最小二乘复指数 法和特征系统实现法要求数据样本长、平均次数多；随机子空间法模型阶次的确定较 为繁琐，在测点较多时 Hankel矩阵阶次很高，所需要的数据采样量较大。频域法的弊病是要求频率分辨率高、样本长，结构是小阻尼的。时频分析法 利用的响应信息太少，是一种局部识别法。工作模态分析不仅在方法求解上还存在局限性，而且在方法考证中，针对的 是较为简单的结构，但即使对于简单结构，现有的方法也不能说解决了所有的问题， 当响应测试数据不完整或者测试数据信噪比较低，现有的方法将

19、会遇到困难。同时 工作模态识别方法同传统的模态识别方法相比， 无论是理论模型、分析手段，还是计算方法都更为复杂，这就可能带来求解上的困难，因此在模型自由度较多 时如何保证数值分析的稳定性，是值得进一步考虑的问题。再者如何将它与有限元分析相结合，以获得更加准确的反映结构在实际运行时的动态特性模型， 也是目前没有解决的问题佝。其中的难点在于，现有的工作模态识别方法所得到的振 型只 是一个相对量，不以质量和刚度阵归一化，同有限元计算结果进行比较有困难。因此 工作模态分析的方法还有待于进一步完善。2 曲轴有限元模型的建立2.1曲轴模态计算要求曲轴是一异形转轴类零件，具有轴线不连续、长径比大、结构、复杂

20、等特点，从 总体上看，曲轴不是对称或反对称体，对曲轴进行有限元模态分析时必须取整 体为研 究对象阴。根据曲轴的结构形状特点，考虑曲轴的实际使用条件，以有限元计算的数据准 备工作量、求解时间及精度等为基本尺度，曲轴模态计算时要求 口：(1)曲轴结构形状复杂，建模时应以不影响其结构动特性为简化原则；(2)有限元模喷网格要有足够的密度，以保证计算结果的精度，真实反映曲轴的模态特性；对曲轴进行模态分析，求解其前9阶模态。2.2曲轴模型的简化曲轴结构形状复杂，小圆角(倒角和细油孔(油道较多，如果在建模时考虑这些因素则会使有限元模型网格非常密，大大增加数据准备工作最和求解时间，并造成单元形状不理想，求解时

21、累积误差增大，运算精度降低 阴冏。在对曲轴进行简化时，参考一些机械结构的模态计算经验，认为小圆角和细油孔对曲轴整体结构动特性影响很小，在建模时忽略半径小于5mnm勺圆角及直径 小于12mm勺油道，协调好计算精度与计算工作量之间的矛盾。2.3曲轴的建模过程该曲轴由6个连杆轴颈与7个主轴颈组成6个单拐，形成曲轴的主体；第1与第 7主轴颈分别向外延伸，形成前端与后端。由于曲轴的很多部分是相似的，所以整体 曲轴模型可用子结构模塑组装建立。本文的曲轴模是在UG中建立的，曲轴的主体模型是根据曲轴的结构特点 先建立 一个单拐模型，再通过旋转、平移、合并而成；最后，将主体模型与前、后端模型合并 成曲轴的整体模

22、型， 如图1所示。图1曲轴的实体模型2.4单元类型的选择根据曲轴的结构形状特点，结合有限元分析系统的单元库中各种单元的优缺点，选取了求解精度比较厲的10结点单元对曲轴进行模态分析。2.5单元划分及节点设置由于只对曲轴作模态分析，综合考虑数据准备最、计算工作量与计算精度，进行网格划分得到：结合点数为89099,单元数为49403。2.6曲轴有限元模型将建立的曲轴模型导ANSYS中，进行有限元模型的建立，得到曲轴的有限元模型，如图2所示。图2曲轴的有限元模型3曲轴模态分析及求解主要是模仿曲轴在磨床上刚刚被夹持的情形，因此要设定曲轴有限元模甩的并提取出前9阶模态，如表1所示，前9阶模态振型，如图35

23、所示表1曲轴前9阶模态参数阶数123456789频率(HZ)169.38181.84397.53472.38513.84760.94796.29832.72917.88HrHnPT5G图3曲轴1阶振型f卜m J! r如%Z1M1OG 袖那屮他心图4曲轴3阶振型如船沁丿小1碱活I神腳6产倾冃2174图5曲轴5阶振型由上面的计算结果和曲轴模态振型可知，在曲轴的前9阶阶次中，其最低的 频率为169.38H乙随着阶次上升，其频率也相应增加。从图3时以看出曲轴的第1阶模态振型(Z69.38HZ)是曲轴在垂直平面YOZ 和XOZ内的一弯振动；从2阶振锢图(图略可以看出，曲轴的第2阶模态振 型 (f=181

24、.84HZ)S曲轴在平面YOZ和XOZ内的一弯振动；从振型图4可以看出，曲轴 的第3阶模态振型(397.53HZ)是曲轴在XOZ平面内的一弯振动，第三缸和第四缸变 形最大；从4阶振型图可以看出，曲轴的第4阶模态振型(f=472.38HZ)是曲轴在平 面YOZ和XOZ内的二弯振动，第三缸和第四缸变形最大；从振型图5可以看出，曲 轴的第5阶模态振型(f=513.84HZ)是曲轴在平面YOZ内的二弯，以 及在平面XOZ 内的二弯振动；从6阶振型图可以看出，曲轴的第6阶模态振型(J760.94HZ)是曲 轴在平面YOZ内的二弯振动，以及Z向的振动。由于曲轴在第7阶、8阶、9阶的 振型变化相对前几阶要小

25、，所以在这里就不再详细叙述。总体来看，在曲轴的9阶 模态振型中，变形最大基本上都出现在第三和第四缸，以及他们之间的连杆轴颈 处。通过上面曲轴的有限元模态分析，从振型图上可以看出，曲轴的低阶频率下 主要 以弯曲振动为主，并且弯曲变形最大的部位出现在连杆轴颈与曲柄臂和平衡块的结 合处。由此可预见，弯曲裂纹是最容易出现在这些地方的血。4 总结数值模态分析与试验模态分析是目前研究结构动力学特性的两大方法，已经成 为解决现代复杂结构动态特性设计的相辅相成的重要手段，在以后的研究过程中， 它们在应用领域会得到进一步的发展。工作模态分析作为目前模态分析领域中一个研究热点，尽管存在某些不足，但 由于它所固有的、在工程应用上的巨大前景和优势，以技术创新和发展为基础，相 信工作模态技术将会有更加广泛的发展和应用。15李德葆振动模态分析及应用北京：宇航出版社，1989:10211816郑敏，申凡，鲍明利用互相尖函数进行环境激励下的模态分析航空学报，200017王利，向阳，卢艳辉，刘文法曲轴的建模及其自由模态分析洛阳高等专科学校学报71 ,2004.14(1):8 1018张勋基UG的结构特征建模方法研究D兰州理工大学学报，200619鱼春燕曲轴分析及其优化设计D.江苏大学学报，200920刘波曲轴参数化建模和有限元分析D:武汉理工大学出版社，2006