2. 有限元离散的基本误差机制

有限元分析的结果与真实解之间总存在差距，理解这个差距从哪里来、有多大、如何度量，是进行误差估计与网格收敛判断的前提。本节从强形式与弱形式的关系出发，梳理有限元离散中两类本质不同的误差来源——投影误差与插值误差，并给出能量范数与L2范数的定义及其物理含义，为后续的误差估计器与自适应网格策略奠定概念基础。

一、从强形式到弱形式：离散误差的起点

以线弹性静力问题为例，强形式要求平衡方程在求解域内逐点成立，同时满足位移边界条件与力边界条件。强形式对解的光滑性要求较高，而工程中的几何、载荷往往只具有较低的正则性，逐点满足既不现实也无必要。弱形式通过引入试探函数与检验函数，将平衡方程乘以检验函数后在域上积分，并利用分部积分把二阶导数转移到检验函数上，从而把对解的光滑性要求降低一阶，只要求解属于某个Sobolev空间（如H¹）。

弱形式的问题提法是：求 u ∈ V（满足本质边界条件的试探空间），使得对任意 v ∈ V₀（检验空间），有 a(u, v) = l(v)。其中 a(·,·) 为双线性型（对应应变能），l(·) 为线性泛函（对应外力功）。有限元方法的核心操作，是把无限维空间 V 替换为有限维子空间 V_h，在 V_h 中寻找 u_h 使弱形式对 V_h 中的检验函数成立。于是真实解 u 与有限元解 u_h 之间的差异，本质上就是“用一个有限维子空间去逼近无限维空间中的解”所产生的误差。

二、投影误差与插值误差：两种不同性质的逼近

有限元解 u_h 具有一个关键性质——Galerkin 正交性：对任意 v_h ∈ V_h，有 a(u − u_h, v_h) = 0。这意味着误差 u − u_h 在能量内积意义下与整个有限元子空间正交。因此 u_h 是 u 在 V_h 中关于能量内积的正交投影，由此产生的误差称为投影误差。投影误差是有限元方法固有的、由离散空间选择决定的误差，它只依赖于子空间 V_h 与真实解 u，而与具体如何构造基函数无关。

插值误差则是另一回事。给定网格节点上的自由度，可以构造一个插值函数 I_h u ∈ V_h，使其在节点上与 u 一致。插值误差 u − I_h u 衡量的是“用分片多项式在节点信息下能多好地表示 u”。经典插值估计给出：对 p 次分片多项式，在网格尺寸 h 下，H¹ 半范数意义下的插值误差为 O(h^p)，L2 范数意义下为 O(h^{p+1})。而由 Céa 引理与 Galerkin 正交性可推出，有限元投影误差在能量范数下不超过插值误差乘以一个与网格无关的常数，即有限元解是“准最优”的。因此工程上常用插值误差的收敛阶来预测有限元误差的收敛阶。

需要强调的是：投影误差 ≤ C × 插值误差，二者量级相同但性质不同。投影误差是实际计算得到的误差，插值误差是可先验估计的上界工具。当网格发生畸变、单元形状变差时，插值常数 C 会增大，投影误差可能明显偏离理想的收敛阶，这正是畸变网格导致精度下降的机制。

三、能量范数与L2范数：两种误差度量及其物理意义

误差需要范数来度量，不同范数回答不同的工程问题。对线弹性问题，能量范数定义为 ‖v‖_E = √(a(v,v))，其中 a(v,v) = ∫ ε(v)ᵀ D ε(v) dΩ。由于 a(u−u_h, u−u_h) 恰好等于真实应变能与有限元应变能之差的两倍，能量范数直接度量了应变能误差，也即应力误差在能量意义下的综合大小。它是有限元误差估计中最自然、最常用的范数，因为 Galerkin 正交性正是在能量内积下成立的，误差估计器（如 Zienkiewicz–Zhu 恢复型估计器）通常也以能量范数为目标量。

L2 范数定义为 ‖v‖_0 = √(∫ v² dΩ)，度量的是位移场本身的整体幅值偏差。它不直接反映应力精度，但对关心位移结果的场合（如刚度、柔度、变形量）更有意义。由于 L2 范数对高频误差不敏感，其收敛阶通常比能量范数高一阶（p 次单元时 L2 为 O(h^{p+1})，能量范数为 O(h^p)），这也解释了为什么位移收敛看起来很快、而应力往往收敛较慢。

工程判断上，若关注应力集中处的峰值应力，应重点看能量范数误差或局部应力恢复误差；若关注整体变形，L2 范数更贴切。两者不可互相替代，报告中应明确所用范数。

四、一个可操作的示例：用插值阶预测收敛

设一维杆问题，长度 L = 1 m，均布轴向载荷，采用 p = 1 的线性单元，网格尺寸 h。理论预测能量范数误差按 O(h) 收敛，L2 范数误差按 O(h²) 收敛。若在 h = 0.1、0.05、0.025 三套网格上计算，得到能量范数误差约为 0.8×10⁻²、0.4×10⁻²、0.2×10⁻²，则相邻网格误差比接近 2，与一阶收敛一致；若某次比值明显偏离 2（如变为 1.3），说明网格畸变或边界处理引入了额外误差，需要检查单元质量与约束施加方式。这一“误差比—收敛阶”的核对方法，是判断离散是否健康的实用手段。