4. 仿真模态分析全流程
仿真模态分析的目标,是求解结构在无外载荷、无阻尼(或忽略阻尼)条件下的固有频率与振型。它与试验模态对标时,真正决定结果可信度的往往不是求解器本身,而是求解之前的建模决策。本节按“几何—网格—材料—边界—求解—后处理”这条主线,说明每个环节中真正影响模态结果的关键动作与判断依据。
一、几何清理:为网格准备“干净”的实体
CAD 模型面向制造与装配,常包含对模态分析无意义甚至有害的细节。几何清理的核心是判断哪些特征保留、哪些简化。典型处理包括:删除直径小于 3~5 mm 的工艺孔、倒角、圆角、螺纹与铭牌凸台;修补导入时产生的碎面、缝隙与重复面;将多个零件合并为单一实体或保留为装配体;对薄壁件抽取中面,对实体件保留三维几何。
判断标准是特征尺寸与关注频段波长的关系。若某细节尺寸远小于所关心最高频率对应的弯曲波长,它对整体刚度与质量的贡献通常可忽略。例如关注 0~2000 Hz 的汽车副车架,2 mm 倒角可安全删除;但若关注局部高频模态,则需保留。清理后应检查几何是否“水密”,否则网格划分会失败或产生自由边。
二、网格划分:单元类型、尺寸与质量
网格决定离散误差。模态分析常用单元包括:实体单元(四面体、六面体)、壳单元、梁单元与质量点单元。薄板类结构优先用壳单元,实体类用六面体或二阶四面体。单元尺寸应保证在关注频率范围内每个弯曲波长至少有 6~10 个单元。经验上,若最高关注频率为 f_max,材料波速为 c,则单元尺寸 L ≤ c/(10·f_max)。
网格质量指标需重点检查:雅可比、翘曲度、长宽比、内角与自由边。模态分析对网格质量比静力学更敏感,因为局部畸形单元会引入虚假局部模态。下表给出常用参考阈值(示例,具体以求解器手册为准):
| 质量指标 | 推荐范围 | 影响 |
|---|---|---|
| 雅可比 | ≥ 0.6 | 过低导致刚度矩阵病态 |
| 长宽比 | ≤ 5(壳) | 过大影响弯曲刚度 |
| 翘曲度 | ≤ 5° | 壳单元翘曲引入虚假刚度 |
| 自由边 | 0(除设计边界) | 非设计自由边意味着连接缺失 |
三、材料属性:弹性模量、密度与泊松比
模态频率对刚度与质量分布高度敏感。频率近似与 √(K/M) 成正比,因此弹性模量 E 与密度 ρ 的误差会直接传递到频率。工程中常用材料示例:钢 E=210 GPa、ρ=7850 kg/m³、ν=0.3;铝合金 E=70 GPa、ρ=2700 kg/m³、ν=0.33;塑料 E=2~3 GPa、ρ=1200 kg/m³、ν=0.35。若材料为复合材料或橡胶,需使用正交各向异性或超弹性模型,并注意温度与频率相关性。
对标试验模态时,材料参数常是首要修正对象。若仿真频率整体偏高,可能是 E 偏大或 ρ 偏小;若整体偏低则相反。建议先核对材料牌号与实测密度,再考虑刚度修正。
四、边界条件:约束、连接与接触
边界条件决定模态的“边界”。自由模态不加任何约束,用于与自由悬挂试验对标;约束模态需按实际安装状态施加。常见约束包括:固定支撑、弹性支撑、对称/反对称边界。连接方式包括:共节点、刚性单元(RBE2)、柔性单元(RBE3)、螺栓连接、胶接与接触。连接刚度不足会显著降低整体频率,是仿真与试验差异的常见来源。
实操中建议:先用共节点或 RBE2 建立刚性连接,快速获得基准模态;再逐步替换为柔性连接或接触,观察频率下降量。若某阶模态频率对连接刚度敏感,说明该连接是薄弱环节,需在试验中重点测量。
五、求解设置:提取阶数、频率范围与求解器
模态求解本质是求解特征值问题 (K − ω²M)φ = 0。常用方法有 Lanczos、子空间迭代与 Block Lanczos。设置要点:提取阶数应覆盖关注频率范围并留 20%~30% 余量;频率范围下限可设 0 或 1 Hz,上限按 f_max 设定;若模型含预应力,需先做静力分析再求解预应力模态。求解后应检查是否有零频刚体模态(自由模态应有 6 阶),以及是否有异常高频局部模态。
六、后处理:频率、振型与 MAC 对标
后处理不只是看频率列表。关键动作包括:查看前若干阶频率与振型,确认振型是否物理合理(如弯曲、扭转、局部呼吸);用 MAC(模态置信准则)矩阵对比仿真与试验振型,MAC 接近 1 表示振型一致;绘制频率误差表,逐阶判断偏差来源。若某阶频率误差超过 5%,需回溯几何、网格、材料或边界条件。
七、实操闭环示例:以某钢制支架为例
假设某钢制支架,材料 E=210 GPa、ρ=7850 kg/m³、ν=0.3,关注 0~1500 Hz。几何清理删除 2 mm 倒角与 M6 螺纹孔;网格采用二阶四面体,单元尺寸 5 mm,检查雅可比 0.72、自由边 0;边界条件为底面 4 个安装孔固定;求解提取前 20 阶,频率范围 1~2000 Hz。求解后得到一阶弯曲 312 Hz、一阶扭转 487 Hz。若试验测得对应频率为 298 Hz 与 470 Hz,误差分别为 4.7% 与 3.6%,MAC 分别为 0.94 与 0.91,可判断该模型在关注频段内可信。若某阶误差达 8%,则应优先检查该阶振型对应的连接区域与材料参数。
综上,仿真模态分析的全流程是一条“建模决策链”:几何清理决定分析对象,网格决定离散精度,材料决定刚度质量比,边界决定边界刚度,求解设置决定提取范围,后处理决定结果判读。每个环节都需留下可追溯的判断依据,才能在与试验模态对标时快速定位差异来源。