1. 模态分析到底在解决什么问题
1.1 模态不是“算一个频率”那么简单
做结构仿真的人,迟早都会碰到模态分析。有限元模态分析这个名头听起来挺学术,但它本质上是回答一个非常朴素的问题:这个结构在什么频率下容易振动,以及它振动起来是什么形状。
先说一个最常见的误解——很多人以为模态分析就是“算一下结构的固有频率”,算完拿几个频率数字就交差。实际项目里,频率只是结果的一部分,甚至不是最重要的部分。真正要关心的是每一阶频率对应的振型(Mode Shape),以及这个振型在真实激励下会不会被激发出来。一个结构可能有几十阶模态,但真正影响产品性能的,往往只有头几阶,或者与激励频率重合的那几阶。
我接触模态分析是从一件很普通的支架开裂问题开始的。那是一个安装在设备内部的薄板支架,说实话结构刚度不算差,但设备一运行起来,支架根部焊缝就开裂。静态强度校核算了很多遍,应力水平也不高,怎么都想不通。后来做了模态分析,发现设备工作转速对应的激振频率正好落在支架某一阶模态的固有频率附近,这一阶振型又是支架末端大幅摆动,根部焊缝附近变形集中。激励频率和固有频率重合,直接引发共振,焊缝反复受载最终开裂。这个案例让我彻底明白,模态分析不是用来“验证强度够不够”的,而是用来判断“这个结构在真实工况下会不会被激励玩坏”的。
所以,对于做结构设计、仿真分析和故障诊断的工程师来说,模态分析是绕不开的基础功。它适合所有人——不管你是做消费电子、工程机械、汽车零部件还是家电产品,只要你的结构会被振动、冲击或交变载荷影响,模态分析就应该出现在你的仿真流程里。
1.2 谁需要模态分析:从设计到故障诊断
模态分析的应用场景比大多数人想的宽得多。
设计阶段,模态分析用于避频。设计工程师关心的是,产品工作时有没有旋转部件,转频是多少,齿轮啮合频率是多少,这些已知的激励频率必须避开结构固有频率。通常的做法是要求结构的最低阶固有频率高于主要激励频率的1.2到1.5倍,或者让固有频率落在激励频率范围之外。
故障诊断阶段,模态分析用于锁定问题源头。结构出现异常振动或异响,光靠感觉很难判断是哪块板子在振,通过模态分析可以得到结构的振型,观察哪个区域位移最大、应变能最集中,基本就能锁定问题区域。
还有一类场景是跟试验对标。样机做出来了,锤击试验或激振器试验能得到试验模态,仿真模态分析可以和试验结果相互印证,修正有限元模型中的材料参数和边界条件。很多企业把仿真模态与试验模态的误差控制在5%以内作为模型可信的标准。
从行业来看,汽车行业的白车身模态、发动机悬置模态,航空航天领域的机翼和舱段模态,电子行业的PCB和散热器模态,家电行业的箱体和压缩机模态,医疗设备里的关键结构件模态,这些都是模态分析的高频应用场景。可以说,凡是会动、会被振动的结构,都有模态分析的用武之地。
2. 理论底子:方程、特征值和三个核心概念
2.1 无阻尼自由振动方程与特征值问题
做模态分析之前,不需要把整个弹性力学从头啃一遍,但有一个方程必须搞清楚,就是无阻尼自由振动方程:
[ [M]{\ddot{u}} + [K]{u} = 0 ]
[M]是质量矩阵,[K]是刚度矩阵,{u}是位移向量。这个方程的物理含义很直白:结构在没有外力和阻尼的情况下,惯性力和弹性力互相平衡,结构做自由振动。
解这个方程,假设位移是简谐运动,也就是{u} = {φ}sin(ωt),代进去就能得到一个经典的特征值问题:
[ ([K] - \omega^2[M]){\phi} = 0 ]
要让这个方程有非零解,系数矩阵的行列式必须等于零,解出来的ω就是固有圆频率,对应的{φ}就是振型向量。这里注意,ω是圆频率,单位是rad/s,平时我们说的固有频率f = ω/(2π),单位是Hz。
为什么要强调这个方程?因为模态分析的所有核心概念——固有频率、振型、模态质量、模态刚度——全都是从这个方程推导出来的。理解了这个方程,后面看任何一款软件的模态分析设置都不会发怵,因为万变不离其宗,求解器做的事情本质上就是组装[M]和[K],然后求解特征值问题。
有限元软件做的具体事情是:把连续体离散成有限个单元,每个单元有自己的质量矩阵和刚度矩阵,组装成整体矩阵后,用数值方法求解特征值。这里有个细节值得注意,质量矩阵有两种常见形式:一致质量矩阵(Consistent Mass Matrix)和集中质量矩阵(Lumped Mass Matrix)。一致质量矩阵由形函数积分得到,更精确但计算量大;集中质量矩阵把质量分布到节点上,是对角矩阵,计算效率高但在某些问题上精度略差。大多数软件默认用一致质量矩阵,对于壳单元和梁单元,这个选择通常没有太大争议。
2.2 固有频率、振型、模态参与质量
模态分析的结果里,你最少要看三样东西。
第一个是固有频率。这个大家都懂,但要注意它跟刚度和质量的关系。刚度越高,频率越高;质量越大,频率越低。这个关系可以直观地用单自由度系统的公式理解:
[ f = \frac{1}{2\pi}\sqrt{\frac{k}{m}} ]
所以,如果模态分析得到的频率比目标值低,你要么加刚度,要么减质量。如果在结构上加筋肋、增加壁厚这些手段改造成本太高,有时候靠调整材料分布也能达到目的。
第二个是振型。振型描述的是结构在某个固有频率下以什么形状振动。第一阶振型往往是整体弯曲,第二阶可能是扭转,第三阶可能是更复杂的弯扭耦合。振型至关重要,因为如果激励力的方向跟振型的运动方向正交,那这个激励就不会激发这一阶模态。举个例子,一个平板的第一阶振型是上下弯曲,如果你在平板平面内施加面内拉力激励,这一阶模态根本不会被激发。
第三个是模态参与质量(Modal Participation Mass / Effective Modal Mass)。这个概念很多人容易忽略,但它才是判断“哪一阶模态重要”的关键指标。模态参与质量表示每一阶模态在某个方向上有多少有效质量参与振动,它跟外部激励方向的分布有关。一个模态虽然频率落在激励范围内,但如果它的参与质量非常小,说明这个模态在对应方向上对整体响应的贡献几乎可以忽略,可以大胆舍弃。
在工程上有个常用经验:在某个方向上,所有有效模态参与质量之和应该达到总质量的90%以上,否则说明你提取的模态阶数不够,模态截断误差太大,需要多提取几阶。
2.3 为什么前几阶很重要:截断与模态叠加
一个连续体结构理论上拥有无穷多阶模态,但有限元模型离散之后,模态数量等于自由度数,依然非常庞大。实际分析中不可能把每一阶模态都算出来,也没有必要。
结构在低频段的表现主要由前几阶模态控制,因为低频模态的振型通常是大范围的整体变形,对结构响应贡献大。高频模态则表现出明显的局部特征,能量分散,对整体响应贡献小,而且高频激励在工程中往往被系统阻尼消耗掉了。
所以模态分析实践中的标准做法是:只提取一定频率范围内的模态,或者提取前N阶模态。对于大多数工程结构,前10到20阶模态已经足够覆盖主要的动态响应。如果你要做模态叠加法(Modal Superposition)的后续分析,比如谐响应分析或瞬态分析,提取的模态阶数是否足够,直接决定后续计算的精度。
这引出一个关键问题:提取多少阶才算够?除了看频率范围,更重要的是看有效模态质量。我在做汽车零部件分析时,客户经常要求提取前20阶模态,但20阶不一定够。稳妥的做法是:先提取前30阶,查看各方向有效模态质量之和占总质量的百分比,如果不足90%,继续增加提取阶数,直到满足要求为止。
3. 实操流程:以某支架类零件为例
3.1 前处理:网格、材料、单位制(最容易翻车的点)
用一个实际例子来说,假设现在要做的是一个设备内部安装支架的模态分析。支架是钣金件,厚度2mm,材料是Q235钢,通过四个螺栓固定在设备机架上。
第一件事不是画网格,而是定单位制。这是我见过最坑的问题之一,很多新手在建模时模型尺寸用的是毫米,但材料参数里的密度用的是kg/m³,弹性模量用的是Pa,最后算出来的频率单位"看起来"是Hz,实际上差了三个数量级都不止。
有限元软件本身不识别单位,只认数字。所以单位制必须自洽,推荐使用一套统一的单位体系:
| 物理量 | 推荐单位制 | 备选单位制 |
|---|---|---|
| 长度 | mm | m |
| 质量 | t(吨) | kg |
| 时间 | s | s |
| 力 | N | N |
| 应力 | MPa | Pa |
| 密度 | t/mm³(即7.85e-9) | kg/m³(即7850) |
| 弹性模量 | MPa(即2.06e5) | Pa(即2.06e11) |
| 频率 | Hz | Hz |
用mm-t-s-MPa这套单位制,钢的密度要输入7.85e-9,弹性模量输入2.06e5。用m-kg-s-Pa这套单位制,密度输入7850,弹性模量输入2.06e11。两套算出来的频率结果完全一致,但混用就废了。
几何清理方面,钣金件的圆角、小孔、翻边、加强筋这些特征对模态结果有影响,但影响程度不同。大孔和加强筋会明显改变刚度分布,建议保留。小圆角和小倒角对低阶模态影响很小,可以简化掉。焊点和螺栓孔要保留,因为后面定义边界条件要用。如果模型里有非常细小的特征导致网格质量很差,优先简化的顺序是:不影响刚度的小圆角、不影响连接关系的小台阶、不影响质量分布的工艺孔。
网格划分的密度对模态结果有直接影响。网格太粗,模型偏刚硬,计算出来的频率偏高;网格太细,计算时间暴增,但对精度的提升趋于饱和。一个实用的判断方法:用两套网格(一套粗、一套细)分别计算第一阶固有频率,如果结果差异在2%以内,说明网格密度已经合适。对于壳单元,单元尺寸建议为厚度的3到5倍以上;对于实体单元,至少保证厚度方向有2到3层单元,否则弯曲刚度偏大,频率会明显偏高。
边界条件建模是另一个重点。支架通过螺栓固定在机架上,在实际建模时要不要把螺栓画出来?通常不需要。简化的做法是在螺栓孔周围建立刚性连接(RBE2)或梁单元模拟螺栓,把螺栓孔中心的节点约束六个自由度。有些情况下,螺栓连接界面的接触刚度会影响模态结果,但这种影响主要在高频段,低频段影响较小。对于初步分析,刚性连接足够;对于需要跟试验精确对标的分析,可以考虑用弹簧单元模拟螺栓连接刚度,弹簧刚度根据螺栓规格和预紧力估算。
3.2 求解设置:频率范围、阶数、提取算法
在求解器里设置模态分析,核心参数有三个:提取方法、频率范围、模态阶数。
提取方法方面,主流有限元软件(Abaqus、ANSYS、Nastran)最常用的是Lanczos法,也叫Block Lanczos法。这个方法对大型稀疏矩阵的特征值求解效率高,而且能很好地处理刚体模态和重频模态。如果你的模型里有刚体模态(比如自由-自由边界条件),Lanczos法依然可以稳定求解。对于中小规模模型,Subspace法也是一种选择,但效率不如Lanczos。我的习惯是,除了极小的教学模型,一律用Lanczos法。
频率范围如何设置?关键看激励来源。设备转速3000rpm,对应的转频是50Hz,那么结构的模态频率至少要覆盖到50Hz的2到3倍,也就是150Hz以上。但更稳妥的做法是,不设置上限或者设置一个偏高的上限,比如500Hz或1000Hz,然后看有效模态质量来决定保留哪些模态。频率范围设得太窄会漏掉重要模态,设得太宽会增加计算时间但通常不会带来精度问题。
模态阶数和频率范围是配合使用的。Nastran里通常设置为提取前20阶,Anonymous里设置为提取频率0到500Hz之间的模态。我的建议是:如果开始不清楚应该提取多少阶,先设置一个较大的阶数,比如前30阶或前50阶,然后查看有效模态质量汇总,如果各方向有效质量之和达到90%以上,说明够了;如果不够,增加阶数重新计算。
求解器的输出设置里,除了位移振型,建议同时输出每阶模态的参与系数和有效模态质量,这对接下来的后处理和后续动力响应计算非常重要。具体输出关键字因软件而异,但都可以在模态分析步的输出请求中直接勾选。
3.3 后处理:看振型、验频率、校核模态质量
模态分析的后处理不像静力分析那样直接看应力云图,而是看三件事。
第一件事是逐一查看各阶振型动画。不要只看频率数字,要看结构在每阶频率下如何变形。检查有没有明显的应力集中区域、变形过大的部位、振型是否连续。如果某个区域振型发生突变,比如出现不自然的尖角变形、穿越、锯齿状振型,通常说明网格有问题或者有未连接的部件。振型动画还有一个作用,就是快速判断它是不是“局部模态”——如果一个支架整体刚度很好,但某个小耳朵或某个薄片在单独剧烈振动,那就是局部模态,这种模态对整体结构响应影响小,但可能恰好被某个局部激励激发,需要留意。
第二件事是验证频率的合理性。例如,一个2mm厚的钣金支架,第一阶固有频率通常在几十到几百Hz之间。如果你算出来第一阶频率是几千Hz,要么是单位制错了,要么是边界条件过约束;如果算出来只有几Hz,要么是约束不足,要么是结构刚度被严重低估。跟手头上类似产品的历史数据对比,是排查频率量级是否合理的最有效方法。
第三件事是检查有效模态质量。以支架为例,正常情况下第一阶模态在X或Y方向(弱刚度方向)的有效模态质量占比最高,可能在30%到60%之间,前几阶模态累加起来,各方向的有效模态质量占比应该达到90%以上。如果在某个方向,所有模态的累计有效质量都不到50%,说明这个方向的约束可能过于刚性,或者提取的模态阶数不足。
这些后处理操作看起来简单,但很多人就是只看一个频率云图就交差了。我第一次做模态分析的时候也是这样,后来被一个有经验的同事提醒,才知道模态分析的关键输出根本不在第一阶频率本身,而在振型、有效质量和频率分布规律。
4. 边界条件与建模细节:决定结果可信度的隐藏因素
4.1 约束与刚体模态
模态分析中最容易出现的问题,就是边界条件设置错误。
一个自由体结构(没有任何约束)拥有六个刚体模态,对应三个平动和三个转动方向,频率理论上为零。如果你的模型设置了自由-自由边界条件,软件会算出六个接近零频率的刚体模态,这是正常现象,不代表模型有问题。
但如果你本意是模拟固定边界,结果却出现了负频率或接近零的频率,那就说明约束不够,结构有机构运动自由度。比如两个零件之间没有建立连接,或者螺栓孔中心节点的自由度没有约束到位。排查的方法是查看零频模态的振型,它能直接告诉你结构在哪个方向发生了刚性位移,从而找到遗漏的约束位置。
边界条件过强和过弱都会造成误差。实际结构中的“固定”几乎不可能是完全刚性的。设备机架本身有柔性,螺栓连接有接触刚度,所以仿真中设置全固定约束,算出来的频率通常会比试验值偏高。如果发现仿真频率普遍比试验值高10%以上,第一个要怀疑的就是边界条件过于刚性。这时候可以用弹簧单元或弹性约束替代刚性固定,弹簧刚度根据实际支撑结构估算。
4.2 螺栓连接、焊接与接触要不要考虑
很多人在模态分析时纠结:螺栓要不要建出来?接触要不要定义?焊缝要不要建模?我的建议是按分析目的区分。
如果是设计初期的方案对比,模态分析只用来比较不同方案的固有频率高低,那么一切连接都可以简化。刚性连接、共节点、MPC都可以,关键是所有方案用同样的简化方式,保证对比的一致性。
如果是跟试验对标、验证模型,那么连接方式的建模就要谨慎了。螺栓连接的接触刚度会降低结构的整体刚度,从而降低固有频率。焊接连接同理,焊点或焊缝的实际刚度小于母材,通过共节点处理焊接连接,会高估局部刚度。
对于需要精确模拟螺栓连接的场景,常见的做法是用弹簧单元(如CBUSH单元)模拟螺栓,弹簧刚度包含拉伸刚度和剪切刚度,可以通过经验公式或微观有限元模型标定。一个M8级别的螺栓连接,其轴向刚度一般在1e7到1e8 N/mm量级。这个值看起来大,但在模态分析中,它对整体结构的低频模态影响有限,对高频局部模态影响显著。
焊接连接的处理跟螺栓类似。如果焊缝是连续焊,可以用共节点或绑定接触(Tie)处理,误差可控。如果焊缝是断续焊或点焊,焊点间距较大,最好用焊点单元(如CWELD或ACM)建模,每个焊点用一根弹簧单元连接两个部件的对应节点,这样能反映焊点间距对刚度的影响。
4.3 阻尼、预紧力、温度等附加影响
标准的模态分析默认是无阻尼的自由振动。但实际结构一定有阻尼,阻尼会影响模态频率和振型,只是影响幅度通常很小。低频结构阻尼比一般在1%到5%之间,对固有频率的影响在0.5%以内,所以绝大多数模态分析不考虑阻尼是合理的。
但有一个例外:如果结构中含有橡胶、聚氨酯等高阻尼材料,或者结构主要通过摩擦界面吸收能量,那么阻尼对频率的影响就不能忽略。橡胶材料的弹性模量本身就随频率变化,而且橡胶配方不同,模量差异很大,直接用钢的材料参数思维套用橡胶,误差会非常大。对于含橡胶的结构,建议用试验数据或者参考相近材料数据库来定义橡胶的弹性模量,并开启复模态分析(Complex Eigenvalue Analysis)以考虑阻尼的影响。
预紧力对模态的影响取决于结构类型。螺栓预紧力对低频模态影响通常较小,但对高频局部模态有影响。更典型的例子是皮带和张紧轮结构,预紧力直接影响结构的等效刚度。对于这类结构,模态分析前需要先做静力分析,把预应力效应考虑进去,使用预应力模态分析(Pre-stressed Modal Analysis)。
温度的影响也值得注意。温度升高导致材料弹性模量下降,结构固有频率下降。金属材料从室温升高到100度,弹性模量可能下降5%左右,固有频率下降约2.5%。如果产品工作环境温差大,模态分析时应该取工作温度下的材料参数,或者分别计算高温和低温工况,看频率漂移是否越过了激励频率。
5. 模态分析常见问题和排查思路
5.1 频率异常低或高:先查单位制和材料
模态分析结果出来,第一阶频率低得离谱(比如0.01Hz)或高得离谱(比如几百万Hz),95%以上是单位制问题。
典型场景:几何模型单位是mm,材料密度输入7850(对应kg/m³),弹性模量输入2.06e11(对应Pa)。这种情况下质量矩阵和刚度矩阵的量级不匹配,导致频率异常。排查方法是把密度和弹性模量做成一个表格,核对每个数值在当前单位制下是否正确。
单位制没问题但频率仍然偏高,常见原因有三个:一是网格太粗导致模型刚硬,细化网格后频率下降;二是边界条件过约束,把应该柔性的连接做成了刚性固定;三是遗漏了附加质量,比如支架上的电机、线束、传感器没有建模,实际质量没有反映到质量矩阵里。
频率偏低则通常是对结构刚度估计不足,比如忽略了某些加强结构,或者网格太细但局部特征(加强筋、翻边)没有体现。还有一种容易被忽略的情况:如果模型中存在未连接的零件,它们以“自由零件”的形式存在,会贡献额外的低阶模态,但这些模态振型只在该零件上出现,对整体结构没有意义,需要检查装配关系。
5.2 振型不连续或出现局部突变
振型云图不光滑、出现锯齿状变化,通常是网格质量或连接定义有问题。
第一种情况是网格过渡不合理,局部单元严重畸变或尺寸突变过大,导致局部刚度分布失真。解决办法是优化网格质量,尽量保证单元尺寸渐近过渡。壳单元要注意单元的翘曲角、长宽比;实体单元要注意Jacobian和偏斜度。
第二种情况是部件之间只有节点重合但节点没有合并,或者连接定义(MPC、绑定的)存在遗漏,导致振型传递断裂。这种情况在装配体中经常发生,排查方法是查看模型的连接关系,确认每个接触面都定义了合适的连接类型。
第三种情况是在同一个位置出现了模态“穿越”,即两阶频率非常接近(重频模态),振型互相叠加,看起来不光滑。这在对称结构中非常常见,比如对称支架的前两阶模态频率相同、振型正交。这不是错误,是结构的对称性导致的物理现象,后处理时分别查看两阶模态即可。
5.3 大量零频或负频模态:约束问题
零频模态的数量和形式能反映约束状态。完全自由的结构有6个刚体模态。如果你的模型只约束了3个平动自由度,没有约束转动,模态分析依然能算,但结果里会出现刚体转动模态,这属于模型设置问题。
出现负频率通常是质量矩阵或刚度矩阵奇异,由过度约束或约束冲突导致。需要检查是否有重复的约束、相互矛盾的边界条件,比如同一节点既被固定又被施加了非零位移。
如果模型中有弹簧单元且刚度为零,也会导致刚体模态或异常低频模态,检查弹簧单元的参数是否合理,是否有悬空的、未连接的弹簧单元。
5.4 与试验对不上:边界条件和传感器附加质量
仿真模态和试验模态对不上,这是做对标项目时最头疼的问题。
先检查试验本身:锤击试验中用到的加速度传感器,即使是几克重的小型传感器,贴在薄壁结构上也会显著改变局部质量和局部模态。一个5g的传感器贴在2mm厚的钣金件上,可能导致局部模态频率偏移5%到10%。你需要在仿真中把传感器质量加到对应位置,重新计算后对比。
再检查边界条件:试验中结构用橡胶绳悬挂,模拟自由-自由边界,这种情况下仿真应该用自由边界条件,而不是固定约束。试验中结构通过工装固定在刚性平台上,工装本身不是绝对刚性,平台的刚度会影响测量结果。仿真中可以考虑把工装也建出来,或者用弹性边界替代理想固定。
还有一个常见问题是材料参数不对。供应商提供的弹性模量数据往往是室温静载下的,实际试验中结构可能处于不同温度或不同应变率下,导致模量有差异。这时候可以用试验频率反推材料弹性模量,对模型进行标定,再用标定后的模型预测其他工况。
5.5 排查顺序速查表
| 现象 | 首要排查项 | 次要排查项 |
|---|---|---|
| 频率普遍偏高 | 单位制错误、网格过粗 | 边界过约束、漏建附加质量 |
| 频率普遍偏低 | 材料弹性模量输入错误 | 连接刚度不足、网格质量差 |
| 出现零频/负频 | 约束不够 | 弹簧单元零刚度 |
| 振型不连续 | 网格畸变 | 部件连接遗漏 |
| 出现重频模态 | 结构对称性(正常现象) | 网格不对称导致虚假重频 |
| 仿真与试验偏差大 | 传感器附加质量 | 试验边界条件与仿真不一致 |
这个表格看起来简单,但在实际项目里,每一次排查都涉及大量的重复计算和结果比对。我的习惯是先把单位制和材料参数检查完,再看边界条件,最后才动网格——因为改网格的代价最大,而且经常不是网格的问题。
6. 从模态到更真实的动力学:延伸方向
6.1 谐响应分析、随机振动与响应谱
模态分析本身是一个基础分析,它的价值很大程度上体现在为后续动力学分析提供模态空间。
谐响应分析(Harmonic Response Analysis)计算结构在正弦激励下的稳态响应。如果激励频率正好靠近某阶模态频率,响应幅值会显著放大。在模态叠加法中,谐响应分析需要先完成模态分析,然后基于模态坐标叠加求解,计算效率很高。工程上常用来评估设备在旋转机械激励下的振动响应,确定危险转速区间。
随机振动分析(Random Vibration,也叫PSD分析)处理的是宽带随机激励,比如道路颠簸、风载、湍流激励。这类分析同样基于模态叠加,在模态空间内计算功率谱密度响应,输出结构在频域内的RMS应力、加速度响应。随机振动分析对模态提取阶数非常敏感,模态截断不够会导致响应被严重低估。
响应谱分析(Response Spectrum Analysis)主要用在抗震和冲击场景。它利用模态分析得到的频率和振型,从响应谱中读出各阶模态对应的最大加速度,再通过模态组合方法(如SRSS、CQC)合成总响应。这个方法效率极高,是工程上评估结构在地震或冲击载荷下响应的标准做法。
可以看出,模态分析是所有频域动力学分析的基石。很多人说“模态分析没什么用,只是算个频率”,其实是因为他们只做了模态分析这一步,没有把它接上后续的响应计算。
6.2 模态测试与有限元对标
模态分析做完之后,如果条件允许,用试验模态做一次对标,是提升模型可信度的最佳路径。
试验模态分析(Experimental Modal Analysis,EMA)的基本流程是:用激振器或力锤激励结构,用加速度传感器测量各测点的响应,通过频响函数(FRF)识别出结构的固有频率、阻尼比和振型。锤击法适合轻小结构,激振器法适合大型结构。
对标时,先用频响函数峰值识别固有频率,跟仿真频率对比。如果偏差小(5%以内),说明模型基本可信;如果偏差大,需要逐项排查前述的因素。振型对标最常用的是MAC(Modal Assurance Criterion)指标,MAC值在0.8以上说明振型相关性好,0.6以下说明振型不匹配。
试验对标还有一个好处,就是能获取真实的阻尼比。仿真模态分析不包含阻尼数据,而后续的谐响应分析和随机振动分析又需要阻尼比,这个参数只能从试验中获取。在仿真中常用的做法是取试验测得的阻尼比的平均值,应用到后续所有模态上。
6.3 拓扑优化与模态约束
模态分析还有一个重要应用场景,就是跟拓扑优化结合。
很多结构优化项目的目标是在减重的同时保证低阶固有频率不低于某个下限。在拓扑优化中,约束条件里可以加入模态频率约束,优化器会在每次迭代中重新计算模态频率,在材料分布上寻找满足频率约束的最优解。
这个场景下,模态分析的精度直接决定优化结果的有效性。优化的中间迭代步中,网格可能产生大量低密度单元,这些单元的存在会导致虚假的局部模态——低密度区域刚度极低,产生很多无关紧要的低频模态。处理方法是设置最小密度阈值,并把频率下限设在合理范围内,避免优化器把注意力放在虚假模态上。
我做过一个电机支架的拓扑优化项目,目标是把第一阶固有频率从250Hz提升到350Hz以上,同时减重15%。优化迭代过程中,第一阶频率从250Hz逐步爬升到380Hz,但减重效果不理想。后来把频率约束改为“前五阶模态的频率均大于300Hz”,同时增加了对称约束,优化收敛的结果才更合理。这说明模态约束的数量和形式对优化结果影响很大,不能只盯着第一阶频率。
这也算是我做了多年模态分析后最深的体会:模态分析不是一个孤立的工作,它的边界条件、提取阶数、质量分布,全都跟后续的响应计算、试验验证、结构优化环环相扣。项目前期多花一点时间把模态分析做扎实,后续的动力学分析会顺畅很多。现在的习惯是,只要拿到一个新结构,我第一件事就是先跑一遍模态分析,哪怕最终不需要做动力学评估——因为它能最快地告诉我这个结构的基本力学性格从哪里来,哪里弱,哪里强,哪里是下一步优化的切入点。