4. 网格划分对模态结果的影响

模态分析求解的是广义特征值问题 (K − ω²M)φ = 0，固有频率由刚度矩阵 K 与质量矩阵 M 共同决定。网格划分之所以会改变结果，根本原因在于：离散后的 K 和 M 只是连续体的近似，网格越粗，离散误差越大，频率偏差也越大。理解这一点，比记住"网格越密越好"更重要——因为不同单元类型、不同阶次、不同密度对 K 和 M 的逼近能力并不相同，误差方向甚至可能相反。

一、网格密度：低频收敛快，高频收敛慢

以悬臂梁为例，设长度 L=500 mm、截面 20 mm×20 mm、材料钢（E=2.0×10¹¹ Pa，ρ=7850 kg/m³，ν=0.3），一端固支。用实体单元（SOLID186）沿长度方向分别划分 5、10、20、40、80 个单元，其余方向保持足够密度，提取前 5 阶弯曲频率，可观察到典型规律：

（以上为按欧拉-伯努利梁理论估算的示例量级，用于说明收敛趋势，非实测数据。）可以看到：1 阶频率在 10～20 个单元时已基本稳定，而 3 阶频率要到 40 个单元才收敛。原因是高阶模态振型波长短，需要更多单元才能描述一个完整正弦半波。工程经验是：若只关心前 2～3 阶，每半个波长至少布置 3～5 个单元；若关心高阶，需按最高阶振型的波长反推单元尺寸。

二、单元类型：实体、壳、梁的差异

同一根梁，用梁单元（BEAM188）、壳单元（SHELL181）和实体单元（SOLID186）建模，结果并不完全一致：

• 梁单元：基于梁理论，直接给出解析刚度与质量分布，自由度少、计算快，对细长梁的弯曲频率精度最高。但它无法捕捉局部变形、截面翘曲和厚度方向模态。
• 壳单元：适合薄壁结构，能反映面内与面外变形，对悬臂梁的弯曲、扭转模态都能较好逼近，但需要合理设置厚度与积分点。
• 实体单元：最通用，能捕捉三维应力状态和局部效应，但自由度最多，且对弯曲问题存在"剪切自锁"风险，需足够多的厚度方向单元。

对细长悬臂梁，梁单元与壳单元的前几阶弯曲频率通常非常接近，实体单元在网格足够密时也趋于一致；但若实体单元在厚度方向只划 1 层，弯曲刚度会被高估，频率偏高。因此，单元类型的选择应服从结构特征：细长杆件优先梁，薄板薄壳优先壳，复杂三维实体才用实体。

三、高阶与低阶单元

低阶单元（如 SOLID185、SHELL181 的线性选项）只有角节点，位移场线性，弯曲刚度偏刚，收敛慢；高阶单元（如 SOLID186、SHELL281）含中间节点，位移场二次，对弯曲和曲率描述更好，同样精度下可用更少单元。代价是自由度和计算量增加。实践中常见策略是：用高阶单元配合中等网格，往往比低阶单元配极密网格更经济，且数值更稳定。

四、网格收敛性检查

收敛性检查不是"看网格好不好看"，而是看结果是否随网格加密趋于稳定。推荐做法：固定单元类型与阶次，逐级加密网格，记录目标阶次频率，计算相邻两级之间的相对变化：

Δ = |f_i − f_(i−1)| / f_i × 100%

判据示例：连续两级 Δ < 1%（工程常用 0.5%～2%）
且频率变化单调趋稳，即可认为该阶已收敛。

若加密后频率反而明显跳动，通常说明：网格畸变过大、厚度方向单元不足、或约束/连接方式不合理，应先排查这些因素，而不是继续盲目加密。

五、悬臂梁对比的完整判断逻辑

把上述因素串起来，可以得到一条清晰的判断链：先按结构特征选单元类型，再按关心阶次的波长定初始网格，然后用高阶单元或加密网格做收敛性检查，最后确认目标阶次频率的相对变化小于判据。下图给出这一逻辑的关系示意。

回到悬臂梁案例：若用 BEAM188 划分 20 个单元，前 3 阶弯曲频率与理论解偏差通常在 1% 以内；若改用 SOLID186 且厚度方向仅 1 层单元，1 阶频率可能偏高 3%～5%，此时应增加厚度方向单元或改用高阶单元。可见，网格划分不是"越密越好"的单一维度问题，而是单元类型、阶次、密度与关心阶次之间的匹配问题。只有通过收敛性检查确认目标阶次稳定，模态结果才具备工程可信度，也才能作为后续谐响应分析的可靠输入。