仿真圈三大玄学问题
仿真圈有三大玄学问题：
- 为什么网格加密之后，结果总往坏的方向走？
- 为什么同一个模型，同事算出来和我不一样？
- 为什么每次我说”结果没问题”，老板都要反问一句”确定？”
其中最玄学的就是第一条：模型一模一样，加密一次，应力从 198 MPa 变 287 MPa，再加密一次又变 312 MPa。你说，这结果到底是哪个？
答案很简单：你之前算的，从来就不是”这个结构的应力”，而是”这一套网格下的应力”。
这个现象的发生有理论背书：
- 位移有限元基于最小势能原理，求得的位移解是真实解的下界（偏刚、偏小）；
- 应力由位移求导得到，粗网格插值太粗糙，在梯度大的地方”抹平”峰值 → 应力集中系数被低估 → 加密后峰值”长出来”。
以下两张图便是同一个方形板上加同样的力的一个模型，但通过两套不同的网格跑出来的结果。怎么样，结果是不是差了很多？
01 | 什么是网格无关性研究？
好了，谜底揭晓一半，我们来看点正经的。
网格无关性研究（Mesh Independence Study），说白了就一句话：
当你的网格密到一定程度，仿真结果不再随着网格加密而发生明显变化，这时候的结果，才勉强能信。
展开来说的话就是：有限元仿真的本质，是把连续体离散成有限个单元来近似求解。网格越细，近似越逼近真实解——但逼近不是无限的。当网格细到某个程度，继续加密对结果的影响已经小到可以忽略（比如应力只变化 0.1%），这时候结果就”稳了”。
打个比方：网格细度就像手机拍照的像素。
- 像素低 → 照片糊成一团，细节全丢（对应结果误差大）
- 像素高 → 细节清晰（对应结果准确）
- 像素从 5000 万加到 1 亿 → 你肉眼根本看不出区别
最后这个状态，就叫”像素无关”，也就是”网格无关”。
但注意：网格无关 ≠ 网格越细越好。 它指的是找到一个”够用”的网格密度，让结果稳定、误差可接受，而不是无脑加密到电脑风扇开始咆哮、求解器卡到天荒地老。
02 | 为什么必须做网格无关性研究？
① 有限元算的是”近似解”，不是”精确解”
真实结构是一个连续体，上面有无数个点；我们只能算有限个节点。网格把连续体”切”成小块，用一个有限维空间去逼近无限维空间——这一步天然引入离散误差：
数值解 = 真实解 + 离散误差
网格越粗，误差越大；网格越细，误差越小，理论上收敛到真实解。如果不做网格无关性研究，你根本不知道自己的误差有多大，甚至不知道自己离真实解是近还是远。
你以为自己在算应力，其实你算的是”这一套网格下的应力”。
② 结果要靠得住，才能用来做决策
仿真结果一旦用于工程决策——判断部件会不会失效、能不能通过验收、要不要改设计——它就变成了”证据“。证据不靠谱，决策就翻车。
做网格无关性研究，本质上是给结果验明正身：我算出来的这个数，不是我画网格画出来的”巧合”，而是接近真实物理的”必然”。
③ 审稿人 / 评审专家 / 甲方会问
如果你要发论文、做认证报告、交设计评审材料，”你做网格无关性验证了吗？“几乎是必考题。
拿不出来？轻则被质疑，重则被拒稿打回。多少仿真人，都是被审稿人一句轻飘飘的话，逼出了人生第一张网格收敛曲线。
03 | 哪些情况需要做网格无关性研究？
直接上表，对号入座：
| 情况 | 是否要做 | 原因 |
|---|---|---|
| 新模型 / 新材料 / 新工况，第一次计算 | ✅ 必须做 | 没有经验数据兜底，结果完全未知 |
| 结果用于工程决策、认证、评审、发论文 | ✅ 必须做 | 结果要”上法庭”当证据 |
| 关注应力峰值、热流密度、裂纹尖端等局部量 | ✅ 必须做 | 局部量对网格极敏感，尤其应力集中/奇异点 |
| 多物理场耦合、强非线性、瞬态大变形 | ✅ 强烈建议 | 各场误差耦合，会传染、会放大 |
| 与实验 / 文献 / 其他软件结果对比验证 | ✅ 建议做 | 消除”网格差异”这个干扰变量 |
| 常规成熟分析，有验证过的网格模板 | ⚠️ 可不重复做 | 模板已传承验证，直接继承结论 |
| 只关心定性趋势（方案谁优、参数怎么变） | ⚠️ 可简化 | 相对比较对绝对误差不敏感 |
| 超大瞬态模型，算一次要一周 | ⚠️ 用替代方案 | 局部加密 + 自适应网格 + 抽查验证 |
重点提醒：应力峰值这类量是”网格钉子户”——尤其在有尖角、约束点、接触边缘的地方，应力可能永远不收敛（理论上的应力奇异点）。
这时候做网格无关性研究，不是为了让应力”绝对收敛”，而是为了确认你关注的位置和数值的物理意义，并明确误差范围。别跟奇异点死磕，该上断裂力学、子模型（Submodeling）就上，成年人要学会换赛道。
番外 | 换个物理场，剧本还一样吗？
前面说的都是结构力学的”剧本”：位移解偏刚、应力是”导数”，粗网格抹平峰值，加密后峰值”长出来”。
但你可能会想：我只做结构吗？流体、传热、电磁的兄弟们，会不会是另一套剧本？
先上结论，再上表：
不管哪个物理场，”全局量稳、局部峰值变”是通用剧本。 积分量（总力、流量、总热流、总电流）对网格天生钝感； 凡是”带导数/带梯度”的局部量（应力、壁面剪切、热流密度、场强峰值），都是网格敏感户，细化网格后峰值大概率往真实值”长”。
| 物理场 | 典型关心量 | 细化网格后峰值一般会 | 坑点提示 |
|---|---|---|---|
| 结构力学 | 应力峰值、位移 | 应力峰值 ↑（偏刚解被修正） | 奇异点永不收敛，上子模型 / 断裂力学 |
| 流场 CFD | 壁面剪切应力、湍动能、分离点 | 局部梯度量 ↑；近壁湍流还看 y+ | 网格无关 ≠ 湍流无关，还得验湍流模型 |
| 传热 | 热流密度、热点温度 | 热流峰值 ↑（温度梯度被”恢复”） | 纯温度场很稳，热流 / 换热系数很敏感 |
| 电磁 | 电场强度、电流密度、损耗密度 | 尖端 / 棱边场强 ↑（场奇异） | 尖端放电、趋肤效应，局部网格需求极高 |
① 流场（CFD）：不止网格，还有”模型”在作妖
流场里，粗网格最大的问题是数值耗散：网格一粗，速度 / 温度的梯度被”糊”掉，相当于给流体加了隐形黏性。
- 壁面剪切应力：对近壁网格（y+）极敏感，粗网格可以低估一大截——做摩擦阻力、压降分析时尤其要小心；
- 湍动能、分离点位置：随网格变化明显，分离点甚至可能”漂移”；
- 反直觉点：流量、总压降这类全局量往往挺稳，粗网格也能算个大概——所以别只看云图漂亮，先盯你最关心的那个数。
而且 CFD 的网格无关性有个著名的”附加条件”：网格无关 ≠ 湍流无关。换一套湍流模型，收敛后的值可能都不一样。网格无关性研究只是”过了第一关”。
② 传热：热流密度是”温度场的导数”
传热问题里，温度本身（尤其是稳态纯导热）对网格相当友好，粗网格也能算个八九不离十。但——
热流密度 q = -k·∇T，是温度的梯度，也就是”导数”。
熟悉的剧本又来了：粗网格把温度梯度抹平 → 热流峰值被低估；加密后，局部热流”长出来”。接触热阻、尖角附近的局部热点尤其明显。如果关心的是”这个接插件会不会过热””芯片热点在哪”，那就老老实实做网格无关性。
③ 电磁：场强峰值是”场奇异专业户”
电磁场是网格敏感度的”极端案例”：
- 尖端、棱边、电极边缘：电场强度理论上趋近无穷大（电场奇异），加密后场强峰值一路”狂飙”——这和结构里的应力奇异点是一个道理；
- 电流密度、焦耳热损耗密度（J²/σ）：都是局部平方量，粗网格能低估到”看起来安全”，一加密就”现原形”；
- 高频电磁波还有额外的”波长约束”：每个波长内必须有足够多的单元（通常 10 个以上），否则波都”传”不动。
所以电磁仿真的网格无关性研究，重点往往不是”全局要不要加密”，而是”关键奇异区域局部加密够不够”。
小结
一句话总结今天的番外：
网格无关性研究，盯的不是网格本身，而是你关心的那个”量”。 全局量，粗网格够用；局部峰值，加密才见真章。
每个物理场的”剧本”略有不同，但主线都是这一条。所以下次被问”网格无关性做了吗”，先想清楚：你监测的是什么量？ 想清楚了，这个回答就值钱了。
04 | 哪些情况可以不做（或者别较真）？
上面表格列了一些，展开说说”可以不做的底气”：
1. 有传承的常规分析
公司里跑了一百遍的成熟工况，网格模板是”祖传”验证过的——这时候没必要每次都从头做网格无关性。但前提是：工况没变、几何没变、边界没变。变了？那还是乖乖做。
2. 只做相对比较
比如优化设计里比较方案 A、B、C 谁更优，用的是同一套网格、同样的离散误差。误差”约掉”了，相对结论依然有效。这种情况做网格无关性是锦上添花，不是雪中送炭。
3. 结果对网格天生不敏感
某些全局量（总力、总体刚度、远场温度）本身对网格就不太挑剔，粗网格也能算个八九不离十。当然，你得先知道它不敏感——这本身往往需要一次研究来确认。
4. 计算资源实在不允许
超大模型 + 瞬态 + 多物理场，算一次要一周。这种就别硬刚全局加密了，用**局部加密 + 自适应网格细化（Adaptive Mesh Refinement）**把资源用在刀刃上，再挑几个关键位置做验证点抽查。
记住一个原则：网格无关性研究服务于”结论可信度”，不服务于”形式主义”。 谁质疑你的结果，谁就需要这份证明；没人质疑、不影响决策，就别为了流程硬凑工作量。
05 | 三步上手：网格无关性研究实操指南
第一步：选定”监测量”
先想清楚：你这个模型到底关心什么？最大应力？最大位移？温度？流量？
选 1~3 个关键量作为”体检指标”，别啥都监测。全都要 = 全都看不清楚，最后只会收获一张眼花缭乱的表格和一句”算了不看了”。
第二步：逐级加密网格
从粗网格开始，每次加密约 1.5~2 倍（比如单元尺寸减半），跑 3~5 个级别。顺手记录：单元数量、单元质量、计算耗时。
加密一时爽，一直加密一直爽——直到你的电脑风扇开始咆哮，就该停了。
第三步：看曲线，定结论
把监测量 vs 网格密度（或单元数）画成曲线：
- 曲线趋于水平 →恭喜，网格无关，取最细且合理的结果；
- 曲线还在明显上升/下降 →继续加密，还没到站；
- 曲线剧烈震荡 → 有”意外情况”（奇异性、接触、材料失稳），先排查问题，别光靠加密硬扛。
经验判据：相邻两档网格，关键量变化 <1%~5%（视精度要求），基本可以认为网格无关。
写在最后：给结果”验明正身”
网格无关性研究，说白了就是一句话：让你的结果从”我画的网格”变成”物理本身的答案”。
它不浪漫，甚至有点枯燥——反复画网格、反复求解、反复对比，像极了调试程序时那句”再加一个 print”。但正是这种”较真”，才让仿真从”画图游戏”升级为”可信工程”。
最后送大家一句话：
网格无关性研究不会让你的结果变得更好看，但它会让你的结果变得更可信。
散会。去画网格吧，记得先存盘。