COMSOL有限元模拟：从几何建模到网格收敛的完整闭环

做有限元模拟的人，第一次面对网格收敛性验证的时候，往往会觉得这是多此一举——我的网格已经够密了，计算结果看起来也合理，为什么还要再加密一倍去验证？

这个问题背后有一个残酷的现实：看起来合理的计算结果，可能是错的，而且错得很有说服力。

有限元的核心误差来源

有限元计算的误差来自三个地方：离散化误差（单元尺寸太大）、截断误差（形函数阶次太低）、数值积分误差（高斯积分点不够）。这三个误差里，唯一可控而且可以量化的是离散化误差——通过网格加密把它压下去。

COMSOL默认用二阶拉格朗日单元（Quadratic Lagrange），形函数本身是二阶精确的，对光滑解来说收敛速度是O(h³)——单元尺寸减半，误差降到1/8。但对有奇点的解（比如应力集中处的尖角），收敛速度会掉到O(h)，甚至完全不收敛。

这就是为什么网格收敛性验证不能省——你不知道你的解里有没有奇点，也不知道当前网格下的误差有多大。

网格收敛性验证的标准流程

第一步：做一个粗网格的计算，记录关心的结果量（比如最大位移、最大Mises应力、最高温度）。

第二步：把全域网格尺寸减半（或者把关心区域的网格局部加密），重新计算，再记录结果量。

第三步：重复第二步，直到连续两次计算的相对误差小于预设阈值（通常3-5%）。

COMSOL里有个很实用的功能叫”自适应网格细化”（Adaptive Mesh Refinement），可以自动做这件事——它根据误差估计（后验误差估计，基于残差或梯度指标）自动在误差大的区域加密网格，比你手动加密更高效。

但要注意：自适应细化的误差估计不一定可靠，特别是对非线性问题。更稳妥的做法是做系统的网格收敛研究（Mesh Convergence Study），手动加密至少三次，画一条结果量随单元数的收敛曲线。

结构力学问题的网格策略

结构力学里，应力集中区域的网格密度直接决定峰值应力的精度。一个经典的例子是带圆孔的平板受拉，孔边的应力集中因子Kt=3，但有限元算出来的峰值应力如果网格不够密，可能只有真实值的60%。

COMSOL里的处理方法：在应力集中区域（圆孔、缺口、接触边角）加”尺寸”节点，把最大单元尺寸设为全局尺寸的1/4-1/8。同时加”边界层”网格，沿应力梯度大的方向（通常是法向）加密。

另一个要点是单元类型选择。结构力学默认用四面体单元，但在薄壁结构里，六面体单元（砖块单元）的精度高得多。COMSOL的”扫掠”网格（Sweep）可以把薄壁区域划成全六面体网格，代价是几何清理的工作量大。

热传导问题的网格策略

热传导方程本身是椭圆型PDE，解的光滑性好，网格收敛通常比结构力学快。但如果有对流边界条件（牛顿冷却），边界层附近的温度梯度大，需要边界层网格。

多物理场耦合的热电问题（比如焦耳热），电流密度和温度场互相耦合，网格需要在电流密度大的区域（接触点、尖角）和温度梯度大的区域同时加密。这两套要求不一定重合，需要分别检查。

自适应网格细化的触发条件

COMSOL的”自适应”功能有两个模式：

• 几何自适应：根据解的梯度指标加密，适合光滑解。
• 残差自适应：根据方程残差加密，适合有激波或者不连续的问题。

实际使用中，几何自适应更稳健。残差自适应容易在奇点附近无限加密（因为残差永远降不下来），需要手动设置最大细化次数和最大单元数限制。

常见坑点

坑一：只在初始网格上跑一遍就出结果。这是最危险的用法。你不知道当前结果的误差有多大，也不知道网格是否足够。

坑二：自适应细化次数不够。默认的自适应细化通常只做2-3次，对有些问题（特别是有应力奇点的问题）不够。手动做网格收敛研究更可靠。

坑三：把不同网格上的结果直接比较而不做插值。粗网格和细网格的节点位置不一样，直接比数值没有意义。正确的做法是把细网格的解插值到粗网格的节点上再比较，或者用积分量（比如总应变能、总热流）来做比较。

计算规模的控制策略

网格收敛性和计算成本是一对矛盾。全域加密的话，单元数按h^(-d)增长（d是维度），三维问题三次加密单元数就涨了27倍。

实际策略是局部加密——只在关心区域和误差大的区域加密，其他地方用粗网格。COMSOL里可以用”尺寸”节点定义多个尺寸域，分别赋给不同的几何域或边界。

另一个策略是用二阶单元。二阶单元的精度比一阶单元高很多，代价是每个单元的计算成本高2-3倍，但总单元数可以少5-10倍，综合成本是更低的。

COMSOL有限元模拟的精度，最终是由你愿意在网格上投入多少计算资源决定的。网格收敛性验证给你的不是一个”精确”的答案，而是一个误差估计——让你知道当前结果的可靠程度在什么范围内。这个误差估计，比一个精确到小数点后三位的数字更有价值。