3. 有限元建模流程与NVH分析类型
在车辆NVH仿真中，有限元模型的质量直接决定后续所有分析结果的可信度。一个完整的NVH有限元模型，需要依次完成几何清理、网格划分、材料属性赋予、连接模拟和边界条件设定五个核心环节。这五个环节并非孤立操作，而是相互制约：网格尺寸影响连接点的选取方式，材料属性决定模态频率的量级，边界条件则直接改变刚度矩阵的奇异性。理解它们之间的耦合关系，比记住操作步骤更重要。
网格划分：NVH分析的精度基础
网格划分的核心矛盾是计算成本与精度的平衡。对于车身结构，NVH关注的中高频段通常在20～200Hz（整车模态）到200～1000Hz（局部板件模态）之间。按照“每个弯曲波长至少6个单元”的经验准则，若某板件在500Hz处出现局部模态，其弯曲波长约为λ = c/f，其中薄板弯曲波速与频率相关，实际工程中常直接采用目标频率上限对应最小波长除以6来估算单元尺寸。例如，某地板在200Hz以下无局部模态，采用20～30mm的壳单元即可；若需捕捉到500Hz，则单元尺寸应缩小到8～12mm。
单元类型选择同样关键。车身覆盖件和梁结构优先使用四边形壳单元，避免三角形单元因刚度过大导致频率偏高。对于厚度方向尺寸远小于其他两向的实体（如悬置支架），可使用六面体单元；若几何复杂无法映射，则用四面体单元但需注意其弯曲刚度偏硬，通常需要更密的网格补偿。网格质量检查中，翘曲度应小于15°，雅可比大于0.6，长宽比控制在5以内，否则模态频率可能出现5%以上的偏差。
材料属性与连接模拟
材料属性方面，NVH分析至少需要弹性模量E、密度ρ、泊松比ν。对于各向同性钢材料，常用E=210GPa，ρ=7.85e-9 t/mm³，ν=0.3；铝合金则取E=70GPa，ρ=2.7e-9 t/mm³。需要注意的是，单位制必须统一，推荐使用mm-t-s体系，此时力单位为N，应力单位为MPa。阻尼在模态分析中通常不直接输入，而在谐响应和瞬态分析中通过模态阻尼比或Rayleigh阻尼定义，车身结构常用模态阻尼比0.02～0.05。
连接模拟是NVH建模中最容易引入误差的环节。点焊在车身中数量庞大，常用CWELD单元或ACM单元模拟，前者通过刚性梁连接两层壳单元节点，后者则生成一个小型实体块。对于粘接结构，可采用实体单元或弹簧单元模拟胶层，弹簧刚度根据胶层剪切模量和面积计算。螺栓连接在NVH中若预紧力足够，可简化为刚性连接或梁单元；若需考虑局部柔度，则用CBUSH单元定义六向刚度。连接点位置必须与几何模型一致，否则会引入虚假局部模态。
边界条件：约束与载荷
边界条件分为约束和载荷两类。模态分析中，自由-自由状态用于识别整车或总成的固有特性，此时需施加惯性释放或至少约束6个刚体自由度；约束模态则用于模拟实际安装状态，例如车身与悬置连接点处施加SPC约束。谐响应分析中，载荷通常为单位力或单位加速度，施加在悬置安装点或方向盘、座椅导轨等响应关注点。瞬态分析则输入时域载荷，如路面激励通过轮胎传递的加速度时间历程。随机振动分析输入功率谱密度PSD，常见于路面谱或发动机振动谱。
分析类型的选择依据
四种NVH分析类型并非随意选用，而是由激励性质和关注目标共同决定。模态分析用于获取结构的固有频率和振型，是其他分析的基础。当关注“某频率下振动有多大”时，若激励为稳态正弦（如发动机定转速阶次激励），选择谐响应分析，输入单位力扫频，输出频响函数，观察峰值频率是否与模态频率重合。若激励为短时冲击（如过减速带、发动机启停），选择瞬态分析，输入时域加速度或力，输出响应时程，判断峰值是否超标以及衰减是否迅速。若激励为随机过程（如路面不平度、湍流），选择随机振动分析，输入功率谱密度，输出响应均方根值，用于疲劳寿命评估或舒适性评价。
建模与分析的实际衔接
以某白车身NVH分析为例：首先对车身几何进行清理，去除直径小于5mm的孔和小于3mm的圆角；然后以10mm壳单元划分，焊点用CWELD单元模拟，材料为钢；约束前后悬置安装点，进行约束模态分析，提取前20阶模态。若一阶扭转频率为32Hz，而发动机怠速二阶激励频率为26.7Hz（800rpm/60×2），则频率间隔仅5.3Hz，存在共振风险。此时可调整加强件厚度或增加加强筋，将一阶扭转频率提升至38Hz以上。随后在悬置点施加单位力，进行20～200Hz谐响应分析，观察方向盘安装点加速度响应，若在32Hz附近出现峰值，则需进一步优化。
需要强调的是，所有分析结果的可信度都建立在建模质量之上。网格过粗会遗漏局部模态，连接刚度过大或过小会改变整体刚度，边界条件与实际不符会导致频率偏移。因此，在完成建模后，建议先进行自由模态检查，确认前6阶为刚体模态且频率接近零；再进行网格收敛性验证，将单元尺寸缩小一半，若关键频率变化小于2%，则认为网格已收敛。只有通过这两项检查，后续的谐响应、瞬态和随机振动分析才具有工程参考价值。