第 9–10 章 · 收敛性、误差与单元性能

算出一个彩色应力云图不难，难的是回答"这个结果可信吗、还差多少"。这一节给有限元结果的质量控制体系：收敛理论、误差估计与自适应。

收敛的理论保证

完备 + 协调的位移元，随网格加密按能量范数单调收敛到真解。误差下降速率由解的光滑性决定：解光滑时，$h$ 型加密下能量范数误差 $O(h^p)$（$p$ 为形函数阶次）；解在裂纹尖端有奇异性（$r^{1/2}$ 型）时，全域均匀加密的收敛速率被奇点拖慢——这正是自适应加密的目标场景。

先验误差估计与后验误差估计

• 先验：在计算前从理论上估计误差阶（$\|e\| \le C h^p |u|_{p+1}$），告诉你"该买几阶单元、网格大致多密"。
• 后验：算完后用数值解自身估计误差分布。最经典的是 Zienkiewicz-Zhu（ZZ）恢复法：把积分点应力磨平成连续应力 $\sigma^*$，误差估计 $\|\sigma^* - \sigma_h\|$——因为磨平应力比分片应力更接近真解（超收敛性）。

自适应分析闭环

循环直到目标量（最大应力、某点位移）误差达标。现代实现是 $h$-$p$ 自适应：梯度大的区域加密（$h$），光滑但要求高的区域升阶（$p$）。p 型收敛（指数速率）在光滑问题上远快于 h 型。

单元性能的检验文化

工业软件的单元库是几十年基准测试的沉淀：悬臂梁（锁死）、厚壁圆筒（对称性）、分片检验（收敛资格）、Cook 膜（翘曲敏感）。看到新单元类型，先问它在这套基准上的表现。

老师的提醒

• 网格收敛性研究是最低成本的质检：加密一倍重算，关键量变化 < 2% 才写报告。跳过这步的分析在同行眼里不值钱。
• 奇点处的"最大应力"没有网格无关性——用名义应力法或断裂力学参数评价，不要报一个随网格变化的假峰值。

自检

• 先验与后验误差估计分别在什么时候给出什么信息？
• ZZ 恢复法利用了什么性质？
• 为什么裂纹尖端的应力值不可直接比较？

更深一层：误差估计把有限元从"算数"变成"科学"

没有误差估计的有限元分析是一份"我算了一个数"的报告；有误差估计的分析是一份"这个数的可信区间是……"的报告——两者的科学地位天差地别。误差估计分两个时机：先验（计算前，从理论给误差上界 $Ch^p\|u\|_{p+1}$，指导选阶次与网格规模）与后验（计算后，从数值解自身估出误差分布，指导局部加密）。工程价值主要在后者：它把"网格够不够"从经验问题变成可计算问题——误差云图上红色区域就是下一轮加密的靶子。

ZZ 恢复法的核心洞察值得单独品味：分片应力（导数）在单元间跳变，但它的磨平版本（超收敛应力）比原解更接近真解——于是"磨平应力 − 原始应力"就是误差的免费估计器。不需要真解、不需要加密计算，一次后处理就得到误差场。这个"从数值解自身提取误差信息"的思想，是自适应方法、目标导向误差估计（对特定工程量估误差）的共同祖先。

推导细读：自适应闭环的三个开关

自适应有限元 = 求解 → 估计 → 标记 → 加密 → 重解的循环，工程实现的关键在三个开关。估计开关：能量范数误差（全局指标，稳健但不指向具体工程量）vs 目标导向（DWR 方法，对"你关心的那个数"——如某点应力、裂缝尖端的 $K_I$——估计误差）；标记开关：Dörfler 标记（对误差贡献最大的前 $\theta$ 比例单元加密，$\theta = 0.5$ 是经典默认）比"超过均值的单元"更保证收敛率；加密开关：h 加密（拆单元，梯度大处）vs p 升阶（解光滑但要求高精度处），h-p 混合在理论上最优、实现上复杂。三个开关的配置决定了自适应是"自动达到目标精度"还是"盲目加密烧算力"。

在真实世界里：网格收敛性研究的行业标准做法

结构报告里的网格研究有通行范式：至少三套网格（粗/中/细，相邻加密比 $\ge 1.5$）、报告关键量的变化率、用 GCI（网格收敛指数，Richardson 外推的工程化）给出"数值不确定度"。ASME V&V 10 类标准把这套流程写成了规范——报告没有 GCI 的仿真结论，在审核眼中等同没有结论。一个实用数字感：二阶方法、网格加密比 2、关键量相邻变化 4 倍关系吻合时，外推值可信度高；变化率杂乱无章则说明还没进渐近收敛区，需要更细网格或先修建模问题。

常见误区与辨析

• "误差估计器说哪里错就加密哪里"：估计器自身有可靠性指标 $\theta$（估计误差/真误差），$\theta$ 远离 1 时估计器失效——先验证估计器再信任它。
• "应力云图平滑 = 结果好"：过度磨平恰恰掩盖误差；看未磨平的单元间跳变量才是诚实的诊断。
• "奇点区域加密到天荒地老"：裂纹尖点应力理论发散，加密只会让它逼近无穷——正确做法是换评价对象（应力强度因子、J 积分），这些量才随网格收敛。

专题深入：ZZ 恢复法的"超收敛"原理

Zienkiewicz-Zhu 恢复法的核心是超收敛性：应力在积分点比其他位置更接近真解（超收敛点）。把积分点应力磨平成连续应力 $\sigma^*$（投影到形函数空间），误差估计 $\|\sigma^* - \sigma_h\|$——"磨平应力比原始应力更接近真解"这个不对称性是 ZZ 的全部秘密。有效性指标 $\theta = \|\hat e\|/\|e\|$（估计误差/真误差应接近 1）是估计器的可靠性证书——$\theta$ 远离 1 时估计器失效，先验证估计器再信任它。

数学细节展开：目标导向误差估计（DWR）

目标导向误差估计（Dual Weighted Residual）把误差估计从"全局能量范数"升级到"你关心的那个数"：对某个工程量 $J(u)$（如某点应力、裂缝尖端的 $K_I$），用对偶问题（adjoint）的解加权残差，估计 $J$ 的误差。"我不关心全场误差，只关心这个数准不准"——DWR 是自适应加密的"精确制导"版本：加密只在你关心的量的敏感区域进行，算力效率的数量级提升。

工程案例细节：GCI——网格收敛指数的行业标准

ASME V&V 10 标准的 GCI（网格收敛指数）流程值得会：三套网格（粗/中/细，加密比 ≥1.5）、Richardson 外推得"网格无关解"、GCI 给数值不确定度区间。报告没有 GCI 的仿真结论在工程审核中等同没有结论——网格收敛性研究不是"精益求精"，是行业法定动作。这个标准在航空航天与核电的仿真审核里是硬性要求。