4. 几何模型与网格离散误差
在有限元、有限体积或时域有限差分等数值仿真中，连续物理场被离散为有限个自由度，这一过程本身就会引入误差。几何模型与网格离散误差，是仿真误差链中最靠前、也最容易被低估的一环。它不同于求解器迭代误差或浮点舍入误差：即使求解器收敛到机器精度，只要几何表达或网格划分与真实物理对象存在偏差，结果仍然会偏离真解。本节围绕几何清理、特征简化、网格类型与密度四个来源，说明它们如何贡献离散误差，并给出网格收敛性的初步判断方法。
一、几何清理：把“看起来对”的模型变成“算得对”的模型
CAD 模型面向制造与装配，常包含大量对仿真无意义甚至有害的细节。几何清理的目标是消除导致网格失败或局部畸变的几何缺陷，而不是改变结构的主要传力路径。常见缺陷与处理方式包括：
- 碎面与短边：由布尔运算或格式转换产生，长度远小于目标网格尺寸。应合并共面碎面、删除长度小于最小网格尺寸的短边，否则会在该处强制生成极小单元，拉低整体时间步或造成矩阵病态。
- 缝隙与重叠：装配体导入后零件间存在微米级缝隙或面片重叠。若缝隙远小于网格尺寸，可缝合或施加绑定约束；若缝隙是真实装配间隙且影响传力，则必须保留并单独处理。
- 小圆角、倒角与孔：半径小于网格尺寸的圆角无法被网格正确捕捉，反而产生大量畸形单元。通常按“圆角半径小于目标单元尺寸的 1/2 时可抑制”的经验处理，但若圆角位于应力集中区，则需保留并局部加密。
- 面法向与自由边：清理后应检查是否存在自由边、非流形边和法向翻转，这些是壳单元和接触分析中常见的报错来源。
几何清理的判断依据不是“模型是否好看”，而是“清理前后关键响应量是否稳定”。建议在清理前后各做一次粗网格试算，比较总质量、质心、关键截面惯性矩等几何等效量，偏差一般应控制在 1% 以内。
二、特征简化：简化的是几何，不是物理
特征简化与几何清理的区别在于：清理针对缺陷，简化针对“真实但可忽略”的细节。简化的合法性取决于该特征对目标响应量的敏感度。例如分析整体弯曲刚度时，直径 2 mm 的螺栓孔对截面惯性矩影响很小，可以忽略；但分析孔边应力集中时，同一孔洞就是核心对象，必须保留并加密。
常用的简化策略有：以等效质量、等效刚度替代小特征；用梁单元或质量点替代细长杆件和集中质量；用二维轴对称或平面应变模型替代三维模型；用对称/反对称边界截取一半或四分之一模型。每一次简化都应记录“简化了什么、依据是什么、预期影响量级”，便于后续误差溯源。
三、网格类型与密度：离散误差的两个主控旋钮
网格类型决定单元对物理场的逼近能力。以结构分析为例，一阶四面体单元（C3D4）在弯曲问题中偏刚，需要极密网格才能达到与二阶六面体单元相当的精度；二阶四面体（C3D10）和二阶六面体（C3D20）对弯曲和应力梯度捕捉更好，但自由度更多。壳单元适合薄壁结构，但要求厚度方向应力假设成立；实体单元通用性强，但薄壁结构需要厚度方向多层单元才能捕捉弯曲。
网格密度则决定离散误差的量级。离散误差随单元尺寸减小而下降，但下降速率取决于单元阶次和问题类型。对光滑解，二阶单元误差约为 O(h²)，一阶单元约为 O(h)。在应力集中、接触边缘、裂纹尖端等梯度剧烈区域，误差下降更慢，必须局部加密。
下面这张图把几何清理、特征简化、网格类型与密度对离散误差的贡献，以及网格收敛性判断的流程组织在一起，便于对照理解。
四、网格收敛性初步判断：用三套网格看趋势
网格收敛性判断的核心思想是：如果网格足够细，继续加密后目标响应量应趋于稳定。工程上常用“三套网格法”做初步判断，具体做法如下。
- 选定目标响应量：只选一个或少数几个关键量，如最大等效应力、最大位移、第一阶固有频率。不要用“整体云图看起来差不多”作为判断依据。
- 生成三套网格：粗、中、细三套，单元尺寸比约为 1.5，例如 6 mm、4 mm、2.7 mm。三套网格的几何清理和简化方案必须完全一致，否则比较的是几何误差而非离散误差。
- 计算相对变化：设粗、中、细网格的响应量分别为 S₁、S₂、S₃，计算 |S₂−S₁|/S₂ 和 |S₃−S₂|/S₃。若后者明显小于前者，说明正在收敛；若两者接近，说明尚未进入收敛区。
- 与阈值比较：工程初步分析常用 5% 作为阈值，高精度分析可取 2%。若 |S₃−S₂|/S₃ 小于阈值，可认为当前网格基本收敛；否则继续加密，并优先加密应力集中、接触边缘、材料界面等梯度剧烈区域。
需要特别提醒：网格收敛并不等于结果正确。如果几何清理时误删了关键特征，或简化时忽略了敏感细节，那么即使网格收敛，收敛到的也是“错误模型的真解”。因此，几何与网格误差控制应遵循“先几何、后网格、再收敛”的顺序，并在每一步保留可追溯的判断依据。
五、一个可复用的检查清单
| 阶段 | 检查项 | 判断依据 |
|---|---|---|
| 几何清理 | 碎面、短边、缝隙、自由边、法向 | 清理前后质量/质心/惯性矩偏差 ≤ 1% |
| 特征简化 | 小圆角、小孔、倒角、细长杆、对称性 | 目标响应量敏感度分析，记录简化依据 |
| 网格类型 | 单元阶次、单元形状、壳/实体选择 | 弯曲问题优先二阶单元；薄壁优先壳单元 |
| 网格密度 | 全局尺寸、局部加密、单元质量 | 雅可比 > 0.6，长宽比 < 5，翘曲角 < 15° |
| 收敛判断 | 三套网格、目标响应量、相对变化 | 细-中网格相对变化 < 5%（或 2%） |
把上述检查项落实到具体项目时，建议在仿真报告中保留三套网格的响应量对比表、几何清理前后的等效量对比，以及简化决策记录。这样，当后续结果出现异常时，可以快速定位误差是来自几何、网格还是求解设置，而不是笼统地归因于“仿真不准”。