网格质量指标与结果收敛之间是什么关系？

网格质量指标全部合格，结果也不一定收敛。本文区分网格质量、求解收敛与网格收敛，梳理常用质量指标与参考阈值、坏单元对结果的影响，并给出一套检查与收敛验证流程。

网格质量与结果收敛是两个相互关联、却不能互相替代的问题：关注区内的网格质量合格，是结果可信的前提之一，但不是充分条件。质量差的单元会拉低精度、拖慢甚至阻止迭代收敛；而质量全部达标的网格，只要分辨率不够，结果仍会随加密继续变化。

本文以 Simcenter Hypermesh 的结构有限元前处理，以及 Simcenter STAR-CCM+ 的 CFD 为背景，梳理常用的网格质量指标及参考阈值、质量差的单元如何影响结果、网格质量与网格收敛的分工，并给出一套可执行的检查与收敛验证流程。边界条件、材料与物理模型带来的误差不在讨论范围内。

本文要点

  • 网格质量衡量单元的形状，网格收敛衡量网格的分辨率，前者合格不代表后者满足。
  • 质量指标没有跨软件的统一标准，阈值必须按所用软件和求解器的定义来理解。
  • 先修复或重新生成网格，消除负体积等无效单元，再处理关注区内的低质量单元，最后做加密对比。
  • 加密后关键结果的变化小于事先约定的阈值，且至少用三套网格确认变化趋势，才能认为网格收敛。

一、先分清三个概念：网格质量、求解收敛与网格收敛

三者常被混着说，回答的却是三个不同的问题：

  • 网格质量：单元的形状是否规整。通过偏斜角、长宽比、雅可比等指标判断，决定求解器能否稳定、精确地处理这个单元。本文统称“单元”，在 CFD 中对应 cell。
  • 求解收敛（迭代收敛）：同一套网格上，求解是否算完。CFD 与非线性分析通过残差下降和监控量是否稳定判断；线性静力分析一次求解完成，不涉及迭代。
  • 网格收敛（网格无关性）：网格加密后，关注的结果是否基本不再变化。目的是确认离散误差已降到可接受的水平。通过关键结果在逐次加密的网格之间的差异是否小于约定阈值判断。

网格质量影响前两者的难易，却不能代替对网格收敛的检验。

还有一个前提：比较各套网格的结果之前，每套网格都必须先完成求解，否则差异里混着未算完的误差。

二、常用网格质量指标与参考阈值

下表列出常用的网格质量指标及参考阈值。同名指标在不同软件、不同求解器里的算法并不相同，阈值只能在对应软件里使用。

指标 衡量什么 参考阈值 软件
偏斜角（Cell Skewness Angle） 面法向与相邻两单元形心连线的夹角，0° 为完全正交；其他软件常称非正交度（non-orthogonality） 大于 85° 视为坏单元；90° 及以上通常引发求解收敛问题 STAR-CCM+
单元质量（Cell Quality） 综合形心分布与面朝向，扁平、面法向与形心连线明显不正交的单元质量低；仅适用于三维网格 1.0 为理想；低于 1.0e-5 视为坏单元 STAR-CCM+
体积变化（Volume Change） 单元体积与其最大相邻单元体积之比 0.01 及以下视为坏单元 STAR-CCM+
面有效性（Face Validity） 各面法向是否都指向单元外侧 低于 1.0 视为坏单元；低于 0.5 表示负体积单元 STAR-CCM+
长宽比 最长边与最短边（或最小归一化高度）之比；三维单元取最差的面 结构单元很少应超过 5:1（不适用于 CFD 棱柱层等边界层网格） Hypermesh
雅可比比值（Jacobian） 最小与最大雅可比行列式之比，反映单元偏离理想形状的程度，范围 0 到 1，1.0 为理想 在高斯点评价时，0.7 及以上一般可接受 Hypermesh
翘曲（Warpage） 单元（或实体单元的面）偏离平面的程度 5° 以内一般可接受 Hypermesh

Hypermesh 会按目标求解器（如 Nastran、Optistruct、Abaqus）分别提供对应的检查方法，所以同一个单元在不同求解器下的评价结果可能不同。负体积、面法向反向等问题属于网格有效性问题，不是“质量好坏”的问题，必须先于质量指标清除。

三、质量差的单元如何影响结果

质量差的网格不一定让求解器报错，却会降低结果的精度和效率，所以“算得出来”不等于“算得对”。影响的路径因分析类型而异：

  • CFD：对流通量要求相邻单元的形心连线尽量穿过公共面的形心，扩散通量要求这条连线尽量垂直于公共面，后者的偏离程度就是偏斜角。偏斜角接近 90° 时，扩散项分母中的点积趋近于零，Simcenter STAR-CCM+ 虽然避免了除零，扩散计算的精度和稳健性仍会下降。体积相差悬殊的相邻单元也会带来误差和不稳定。
  • 结构有限元：单元畸形时，雅可比行列式在单元内部的变化加大，极端情况下出现负值，单元发生翻转。单元内的位移插值和应变计算精度随之下降，应力梯度大的区域受影响最明显。
  • 补救手段的边界：STAR-CCM+ 的单元质量修复模型（Cell Quality Remediation）会标记低质量单元及其邻居，并修正这些单元的梯度计算，以提高求解的稳健性。它让带坏单元的网格能够算下去，但单元本身并没有变好。

四、网格质量达标，结果就一定收敛吗

不一定。质量指标只看单元自身的形状，不看单元够不够小。评价 CFD 网格时，STAR-CCM+ 把网格密度（高梯度区是否足够密）、网格分布和近壁层，与偏斜角、体积比等单元级指标并列为不同的考察项，说明分辨率和单元质量本来就是两件事。两个维度交叉，有四种情况：

网格质量 \ 分辨率 分辨率足够 分辨率不足
质量达标 可以进入收敛对比，确认收敛后结果才可信 最隐蔽：求解稳定、云图平滑，结果却随加密持续变化
质量不达标 能算出结果，但精度、效率和稳定性受影响，关注区内的坏单元必须处理 既不稳定，也不准确

反过来，坏单元也不必全部消灭。STAR-CCM+ 文档说明，取决于所选的物理模型，单元质量较低时仍可能得到有效解。

工程经验坏单元的位置比数量更重要。远离关注区的少量坏单元通常影响有限，关注区、壁面附近和载荷传递路径上的坏单元应优先处理。这一判断来自工程经验，不是软件的官方判据。

五、可执行流程：从网格检查到收敛验证

下面的流程适用于结构和 CFD 分析：先确定看什么、差多少算收敛，再检查网格，最后用加密对比验证。下面的选择器可以按当前结果给出下一步去向，完整步骤见其后的文字版。

交互工具

收敛判断选择器

每套网格先完成求解再比较；回答两个问题，得到下一步去向。

1　加密后，关键结果的变化小于事先约定的阈值吗？
查看完整流程（文字版）
  1. 明确关注量与阈值确定要对比的关键结果，并约定可接受的变化阈值。
  2. 生成基准网格修复负体积等无效单元。
  3. 检查质量处理关注区内的坏单元。
  4. 求解CFD 与非线性分析需达到迭代收敛；每套网格都先完成求解，再比较结果。
  5. 加密并重新求解网格收敛至少用三套网格（两次加密）。
  6. 判断变化是否小于阈值小于阈值，采用该网格，并记录阈值、网格规模与结果。
  7. 否则判断是否奇异点是奇异点，修正几何或载荷施加方式后回到质量检查；不是奇异点，继续加密并重新求解。
工程经验阈值应按分析目的事先约定，概念阶段可以放宽，设计验证应更严。单元尺寸每次缩小约 1.5 到 2 倍，差异才容易与数值噪声区分开；加密后还要重新检查质量，过渡区可能产生新的坏单元。网格收敛至少用三套网格（两次加密）才能判断变化趋势。这些做法来自工程经验。

六、不同分析类型的关注重点

同一套流程，在不同分析中盯的指标和结果各不相同：

分析类型 重点检查的网格质量 收敛验证时重点对比
结构静力（应力） 关注区内的雅可比、翘曲、长宽比；三维单元按最差的面评价 峰值应力、位移、反力；峰值应力最难收敛，需要局部加密
稳态 CFD 偏斜角、体积变化 压降、受力、壁面换热量等积分量，而不是某一点的速度；同时确认近壁层分辨率（y+）满足所选湍流模型的要求
高保真 CFD（LES/DES、声学） 在上述基础上增加 STAR-CCM+ 的面翘曲质量（Warpage Quality）：低于 0.15（对应翘曲角约 50°）视为坏单元。它与第二章 Hypermesh 的翘曲角不是同一指标，阈值不可互相套用 时间平均量与脉动统计量
工程经验位移和反力通常比应力更早稳定，所以只对比位移就判定网格收敛，容易高估网格的可靠程度。这是工程经验。

七、三个常见误区

  • 所有指标都合格，网格就收敛了。质量指标不度量分辨率。合格的网格仍然需要加密对比（至少三套网格），才能确认结果不再随网格变化。
  • 仅靠单元质量修复模型或自动移除无效单元把计算跑通，问题就解决了。这些手段提高的是求解的稳健性，关注区的网格并没有变好。应回到前处理工具中修复，再重新检查。
  • 应力随加密一直上升，就继续加密。点约束、点载荷、凹角（重入角）等引起的奇异点，其处的应力不会收敛，加密只会让数值越来越大。应先判断是真实的应力集中还是奇异点，奇异点要靠修正几何（如加圆角）或载荷施加方式解决，并且在远离奇异点的位置做收敛对比。

八、结语

网格质量和网格收敛各管一件事：先用质量指标保证每个单元可用、关注区没有明显的坏单元，再用加密对比确认分辨率足够，最后把阈值、网格规模和对比结果写进分析记录。

单元类型的选择同样影响收敛行为，可参考《壳单元、四面体与六面体网格应该如何选择？》；网格之外的误差来源，可参考《仿真与测试结果不一致时，应该如何定位差异？》。

如果您正在使用 Simcenter Hypermesh 检查和改善网格，或希望评估现有模型的网格与收敛性，可以了解数领科技的 工程仿真与物理测试 服务。

相关文章

继续阅读