1. 课程导入与误差认知

在有限元分析的实际工程中，一个常见的困惑是：同一个零件、同一套边界条件，不同工程师算出来的应力却相差百分之二三十。问题往往不在软件操作，而在于对误差来源缺乏系统认知。有限元从来不是"输入模型、点击求解、读取结果"的线性动作，它是一条从物理现实到离散方程、再到数值解的近似链条，每一环都会引入误差。本课程的目标，就是让这条链条上的误差变得可识别、可量化、可控制。

误差从哪里来：三类误差的层次关系

有限元计算的总误差可以按来源划分为三个层次，它们并非并列关系，而是依次叠加、逐级传递的。

模型误差位于最上游，指数学模型与真实物理现象之间的偏差。例如把三维实体简化为平面应力问题、忽略接触面的摩擦、用线弹性本构替代塑性行为、把分布载荷等效为集中力。这类误差在网格加密时不会减小，反而可能因为计算精度提高而更"忠实"地呈现一个错误模型的结果。工程上判断模型误差是否可接受，通常依靠解析解对比、试验数据标定或同类问题的经验区间。

离散误差是有限元特有的核心误差，来源于用有限维形函数空间逼近连续解。它包含两部分：一是单元尺寸导致的逼近误差，网格越粗，单元内位移场越难贴合真实变形；二是单元形状与阶次导致的插值误差，畸变单元、低阶四面体会显著放大这一项。离散误差随网格加密单调下降，这正是网格收敛研究的对象。

求解误差发生在离散方程组的数值求解阶段，包括迭代求解器的收敛容差、病态刚度矩阵引起的舍入误差累积、接触迭代未完全收敛等。它通常远小于前两者，但在大变形、强非线性或极细网格问题中可能上升为不可忽视的量。

为什么必须先建立误差认知

很多初学者把"网格加密"当成提高精度的万能手段，结果在应力集中处反复加密，计算量翻了几倍，关键部位的应力却仍在漂移。原因在于：如果模型误差占主导，加密网格只会让你更精确地算错；如果加密后结果仍不收敛，往往说明存在奇异点（尖角、点载荷、刚性约束边），此时应力理论上趋于无穷，任何网格都无法收敛，需要改用局部细化加应力外推或改用结构应力法。

因此，误差估计与网格收敛不是仿真流程末尾的"可选检查"，而是贯穿建模决策的验证环节。它回答三个问题：当前结果离真实解还有多远？加密网格还能带来多少改善？在给定工程精度要求下，多粗的网格就够用？

一个可操作的收敛判断示例

以某支架在给定载荷下的最大等效应力为例（以下数值为教学假设，非实测）。采用二阶四面体单元，依次取全局尺寸 8 mm、4 mm、2 mm、1 mm 四套网格，得到最大应力分别为 182 MPa、205 MPa、216 MPa、221 MPa。相邻两级的相对变化为 12.6%、5.4%、2.3%，呈明显递减趋势，说明解正在收敛。若工程允许误差为 5%，则 2 mm 网格已满足要求；若要求 2%，则需继续加密或改用外推估计。

网格尺寸(mm)   最大等效应力(MPa)   相邻级相对变化
    8              182.0              —
    4              205.0            12.6%
    2              216.0             5.4%
    1              221.0             2.3%

这个简单表格体现了网格收敛研究的基本逻辑：不看单次结果，看结果随离散参数变化的趋势。趋势平稳，才谈得上可信。

课程目标与学习路径

本课程围绕"识别误差—估计误差—控制误差"这条主线展开。前半部分建立理论基础：误差分类、收敛阶、先验与后验误差估计、h/p 收敛与奇异点处理；中段进入工程方法：网格收敛指数（GCI）、Richardson 外推、局部细化策略、应力奇异与应力线性化的判别；后段结合实操：在主流求解器中设置收敛研究、读取误差指标、编写批量算例并整理收敛曲线。学习路径建议按"概念—判据—操作—复盘"推进，每章都配一个可复现的小算例，最终形成一套适用于自己产品线的网格验证规范。