第 2–3 章 · 从强形式到弱形式

连续体没有"节点"和"弹簧"，怎么把上一节的框架搬过来？答案是弱形式。这一节是有限元的数学心脏：把微分方程改写成积分形式，离散化就变成了"在有限维空间里找最优逼近"。

强形式为什么不好离散

以一维杆的平衡方程（强形式）为例：

它要求解二阶可微——但工程里的集中力、材料突变处导数根本不连续。要求太高的光滑性，是强形式的死穴。

弱化：乘上试函数再积分

把方程乘以任意试函数 $v$（满足 $v=0$ 于位移给定边界），分部积分：

这一步有两个奇迹：导数阶数从二降到一（对光滑性要求减半），且边界力项自然出现（自然边界条件自动包含在积分里，不必单独施加）。分部积分不是技巧，是把物理平衡从"逐点成立"弱化成"加权平均成立"。

伽辽金法：有限维逼近

取试函数与检验函数同属一个有限维空间（形函数张成的空间），令残差在检验函数方向上为零：

这就是伽辽金有限元。对称的微分算子给出对称的 $\mathbf{K}$——上一节的对称性在数学上有了出处。弱形式的另一个别名是"虚功原理"：$v$ 是虚位移，方程就是内力虚功 = 外力虚功。

老师的提醒

• 记住对称性来源：对称算子 + 检验空间 = 试空间 → 对称刚度阵。非对称格式（如对流占优的迎风有限元）正是破坏了这个等式。
• 弱形式的边界条件分两种：位移给定的本质边界（须强制施加于试探空间）与力给定的自然边界（分部积分自动带走）。混用是初学者第一大错。

自检

• 分部积分带来哪两大好处？
• 本质边界与自然边界分别怎么施加？
• 伽辽金法与虚功原理是什么关系？

更深一层：弱形式是"对光滑性的讨价还价"

强形式要求解二阶可微，弱形式只要求一阶（能量范数有定义即可）。这次"讨价还价"换来的是三重自由：允许集中力、材料界面、裂纹这类导数不连续的解进入框架；允许使用低次多项式这类不光滑的试探函数做逼近；让自然边界条件自动嵌入积分式。换句话说，强形式描述的是"理想光滑世界的平衡"，弱形式描述的是"真实工程世界的平衡"——后者才是工程师需要的。

弱形式还有一个被低估的身份：它是数值方法的统一语言。有限差分基于强形式（逐点近似导数），有限元基于伽辽金弱形式，有限体积可理解为"对偶权函数取 1 的弱形式"（每格守恒），谱方法用全局正交函数做试探空间。四种主流离散方法的差异，本质上是"用哪组检验函数与试探空间"的差异——抓住弱形式，四大家族就在同一张地图上了。

推导细读：分部积分的边界项里藏着什么

对 $-\int (EAu')'v\,dx$ 分部积分，边界项为 $-[EAu'v]_0^L$。逐端分析：

• 固定端（$u = 0$ 已知）：检验函数满足 $v(0) = 0$，边界项自动消失——本质边界通过试探空间消化；
• 自由端（力未知）：边界项保留 $EAu'(L)v(L)$，令其等于外加力 $Pv(L)$——自然边界条件在积分里现形；
• 弹簧支撑端：边界项 $= -kv(L)v(L)$，移入左端变成刚度贡献——弹性边界变成了对称刚度项。

三种端部条件，一条分部积分全部安放妥当。更深一层：边界项的处理方式决定了弱形式是否对称、刚度阵是否正定——弹性支撑让刚度阵"更正定"，而自由端加上刚体模态则让它奇异（欠约束）。边界条件从来不是"套公式"，而是弱形式里活的结构。

在真实世界里：加权残量法族的实战分工

第 2 章把配点法、最小二乘法、伽辽金法统一在加权残量框架下，工程上的分工大致是：伽辽金（检验=试探）是结构/传热有限元的默认，对称问题的最优选择；最小二乘（检验=残量对参数的导数）在对流占优与非自伴算子（如双曲方程、一阶系统）上比伽辽金稳定，近年 LSFEM 有复兴之势；配点法在实时场景（模型降阶后的数字孪生）因其"不求积分、只算点值"的速度而被重新启用。选哪种加权方式，取决于算子的对称性与你愿意为积分付出的代价——这条判据比任何教材的排序都实用。

常见误区与辨析

• "弱形式的解就是'差点儿'的解"：在能量范数下，伽辽金解是有限维空间里的最优逼近（投影定理）——不是凑合，是子空间里的冠军。
• "本质边界也可以用罚函数糊弄"：罚函数可行但引入 $O(1/\alpha)$ 误差与条件数恶化；高精度分析（断裂、接触）建议拉格朗日乘子或 Nitsche 方法。
• "检验函数取什么无所谓"：伽辽金之所以给出对称刚度，正因为检验空间 = 试探空间；换成 Petrov-Galerkin（两者不同）就换来稳定性（迎风）——无免费午餐。

专题深入：Petrov-Galerkin 与迎风的"稳定性改造"

伽辽金（检验 = 试探）在对流占优问题上失效——中心差分式的伽辽金离散对流项会振荡（与 CFD 的中心差分不稳定同源）。Petrov-Galerkin（检验空间 ≠ 试探空间）的迎风有限元（SUPG——Streamline Upwind Petrov-Galerkin）在检验函数里加迎风偏向项，恢复稳定性而不牺牲精度（一致性保持）。对称格式 ⟹ 对称刚度阵（SPD，CG 可解）；非对称格式 ⟹ 非对称刚度阵（GMRES/BiCG）——稳定性改造的代价是失去对称性，这是"迎风"在有限元世界的精确代价。迎风的本质是"让信息来源决定检验方向"——与 CFD 的迎风有限差分共享同一个物理直觉。

数学细节展开：误差估计的"Céa 引理"——伽辽金的最优性证书

伽辽金逼近的最优性由 Céa 引理给出：能量范数下，伽辽金解是有限维子空间中的最优逼近（误差 ≤ 常数 × 最优逼近误差）。证明只有两行（弱形式的正交性 + 强制性/连续性），但它的地位是"有限元的数学宪法"：收敛性与精度被归约为"子空间逼近能力"——形函数阶数、网格密度、解的正则性全部通过"逼近能力"这一个量进入误差估计。这就是为什么形函数阶数与解的正则性是收敛率的两个旋钮：解光滑时 $h$ 型加密 $O(h^p)$、解有奇异性（裂纹）时全域均匀加密被奇点拖慢——自适应加密的目标场景。Céa 引理把"算法好不好"翻译成"子空间大不大"。

工程案例细节：弱形式在"非位移型"问题里的变体

弱形式在非位移型问题里的变体清单：传热（温度标量场，$B$ 换成温度梯度形函数——最干净的第一个非线性/场问题练习）；渗流-固结（Biot 耦合：位移 + 孔压双场——岩土工程的核心）；电磁（棱边单元防伪解——$H(\mathrm{curl})$ 空间的特殊插值）；声学（无限元/完美匹配层截断开放域）。每个场问题的"新东西"几乎都是：单元插值空间的选择与稳定化——$H(\mathrm{div})$、$H(\mathrm{curl})$ 的插入空间理论是混合与电磁有限元的数学基础。

常见误区与辨析（补充）

• "弱形式比强形式'弱'所以不精确"：恰恰相反——弱形式允许更广泛的解（包括强形式无定义的间断解），它是"更诚实"而不是"更模糊"。
• "本质边界可以用大刚度模拟"：罚法可行但引入 $O(1/\alpha)$ 误差与条件数恶化——高精度分析建议乘子或 Nitsche。