弱形式入门：从偏微分方程到有限元仿真实战

简介：面向COMSOL Multiphysics用户的弱形式入门文档，采用浅显讲解方式，系统比较偏微分方程、能量最小化形式与弱形式三种物理问题描述方法，帮助读者摆脱繁琐数学细节，直接利用弱形式界面处理非线性、多物理场等难以用内置模块或经典PDE模板求解的场景。文档共1个doc文件，压缩包总大小约695KB，篇幅集中、重点突出。内容结合线弹性结构分析中的Navier方程与胡克定律，演示如何将PDE方程转化为积分弱形式，并解释试函数、泛函、弹性张量等关键概念，使读者理解有限元方法离散化求解的本质。文档还梳理了COMSOL手册中的常用术语与矩阵标记方式，便于后续阅读官方算例和模型库示例。该资料已有1177人学习下载，适合希望深入理解COMSOL求解机制、需要自主编写弱形式方程以扩展建模能力的工程师、科研人员及有限元教学者。

1. 先放下恐惧：弱形式不过是换一种写法

很多刚接触仿真软件的人，一听到“弱形式”三个字就开始紧张，觉得这一定是数学系才会的东西。我第一次也一样，翻了几页文档就被积分号劝退。实际上弱形式不是一种新的物理定律，它只是一个“换种写法”的手段：同一个物理问题，可以写成偏微分方程，可以写成能量最小化，也可以写成弱形式。你不需要先成为数学家才能下手。这篇入门就是带你把弱形式这个概念讲懂，再在软件里输入第一个能跑的表达式。它适合所有被内置模块限制住、想自己定义方程的新手。

2. 为什么绕不开弱形式：三种描述同一个物理问题的方式

2.1 偏微分方程、能量最小化、弱形式：三张面孔

先记住一句话：物理问题有三种写法，本质是同一个东西。第一种写法是偏微分方程，简称 PDE，你在课本上最常见到它。第二种是能量最小化，比如一根弹簧被拉伸，它储存的弹性能会努力让自己最小化。第三种就是弱形式，它像是一个更宽容的积分版本，专门用来处理那些不太听话的变量。

为什么非得知道这些？因为仿真软件的内部逻辑就是围绕这三种写法转的。很多模块已经帮你内置了常见的 PDE，比如传热、结构力学，点几个按钮就能算。可一旦问题变成“内置模块没有的方程”，比如某种特殊的扩散、流体和传质耦合，你就得回到最底层，自己把方程输进去。而软件底层的求解器，恰恰是把所有方程都先转成弱形式再做离散。不理解弱形式，你就只能看着高级设置里的一堆选项发愣。

2.2 从能量最小化切入：弹性能泛函是“函数里的函数”

先别急着接触弱形式，我们从更好理解的能量最小化开始。想象一根被拉长的杆子，内部储存了弹性能。物理规律说，系统总是倾向于让总势能最小。这个总势能是一个数，但它不是简单地把几个数字相加，而是把一个“函数”作为输入，输出一个“能量值”。这种输入是函数、输出是数值的东西，在数学上叫泛函。

泛函听着玄，其实可以类比成“函数的函数”。普通函数是输入一个数，得到另一个数；泛函是输入一条曲线或一个函数，得到一个数。比如你测量一段绳子弯曲后的长度，输入的是绳子形状这个函数，输出的是长度这个数。弹性能泛函就是输入位移分布 u，输出总能量 W。有限元里经常出现这种“函数套函数”的东西，理解它是理解弱形式的一半。

2.3 弱形式为什么更“宽容”：连续性要求更低

PDE 形式要求变量在定义域内要足够光滑，导数要存在，因为方程里有散度、梯度这些偏微分算子。可实际问题里，材料可能由不同介质组成，场变量在界面处可能突变。这时候 PDE 形式就显得太苛刻了。

弱形式的做法是：把 PDE 乘上一个试函数，再在整个区域上积分，然后通过分部积分把求导负担分一半给试函数。这样一来，对原变量的可微要求就降低了一阶。很多原本在 PDE 形式下“没法写”的问题，在弱形式下就能正常求解。这也是为什么弱形式特别适合非线性多物理场问题——流体里的 N-S 方程、化工里的传质问题，往往没有能量可以最小化，但弱形式照样能用。

2.4 为什么说内置模块也在用弱形式

有个事实可能会让你重新看待这个软件：所有应用模式和 PDE 模式在求解时，都会先把方程式系统转成弱形式，再进行有限元离散。也就是说，弱形式不是“高级用户才碰的偏门功能”，而是这款软件最底层的通用语言。

这意味着，你平时用的固体传热、结构力学模块，界面是友好的，但内核仍然在走弱形式这条路。当你有一天打开模型库，看到一个用弱形式做的特殊算例，你不会再觉得那是天书。反过来，你也会慢慢理解为什么有些求解器设置、网格要求、边界条件写法是那样设计的——因为底层都是同一套积分和离散的逻辑。

3. 第一步动手：在软件里输入一个真正的弱形式

3.1 选对流-扩散方程作为第一个例子

入门第一个例子，我建议你选对流-扩散方程。原因很简单：它没有对应的能量最小化形式，只有弱形式能干净地处理，最能让你体会到弱形式“不可替代”的价值。

方程长这样：

-∇·(c∇u) + P·∇u + a*u = f

这里u是待求的浓度场，c是扩散系数，P是对流速度矢量，a是反应或吸收系数，f是源项。如果你觉得符号陌生，就把它想成一个同时受“扩散”和“平流”影响的浓度分布问题。扩散让浓度从高往低走，对流让浓度随着流场搬走。

3.2 弱形式推导：乘试函数，分部积分

我们现在动手把它转成弱形式。流程一共三步，非常机械。

第一步，给方程乘上一个试函数v。试函数你可以暂时理解为“一个用来测试方程是否成立的任意小扰动”。

第二步，在计算区域上做积分。

第三步，对含有∇·的扩散项做分部积分，把二阶导数拆开。具体操作是：把带散度的项变成边界项和体积项，体积项里原来的梯度现在分了一半给试函数。

最终得到弱形式表达式如下：

c*(ux*test(ux) + uy*test(uy)) + (Px*ux + Py*uy)*test(u) + a*u*test(u) - f*test(u)

这里用二维举例，所以只写了 x 和 y 方向。如果你做三维，还要再加上uz和test(uz)。

3.3 在弱形式界面里的输入：一个代码块

打开软件，新建一个模型，在“模型向导”里选择“二维”，然后在“添加物理场”的搜索框里输入“弱形式”，找到“弱形式 PDE (Weak Form PDE)”，添加它。

在“变量”设置里，把因变量名字保持为u。然后在“弱形式 PDE”节点下的“弱表达式”输入框里，粘贴下面这段表达式：

c*(ux*test(ux) + uy*test(uy)) + (Px*ux + Py*uy)*test(u) + a*u*test(u) - f*test(u)

注意：如果界面里没有ux、uy这些变量，需要先在“变量”节点里定义它们。通常软件会自动识别u对坐标的导数，test(u)表示试函数，test(ux)表示试函数的 x 方向导数。

这段表达式是什么意思呢？我拆开给你看：

• c*(ux*test(ux) + uy*test(uy))：扩散项的弱贡献，也是你用分部积分之后留下来的那个体积积分。
• (Px*ux + Py*uy)*test(u)：对流项。你会发现这里u的梯度被保留了，但试函数test(u)没有梯度，这正是“不可重排”的关键——它让刚度矩阵不对称。
• a*u*test(u)：反应或吸收项。
• - f*test(u)：源项，相当于给系统一个外力。

如果你把 P、a、f 都设成 0，只留扩散项，本质就是求解稳态扩散方程。你可以用这个最简单的设置做第一个算例：左边边界设浓度 1，右边边界设浓度 0，中间区域设扩散系数c = 1，算完后浓度会从 1 线性降到 0。这个结果和手算解析解完全一致，能立刻验证你的设置没毛病。

3.4 区域积分和边界积分分开设置

细心的你可能会问：分部积分得到的边界项跑哪儿去了？答案是：边界项单独放在“边界弱表达式”里。

在“弱形式 PDE”节点下，有一个“边界”子节点。这里可以输入边界上的弱贡献。比如你要在某个边界施加一个进入区域的通量flux，就在边界弱表达式里输入：

-flux*test(u)

负号是因为默认的法向指向外侧。软件会把区域积分和边界积分分开用不同的算法处理，因为二者的几何维度不同。你不需要自己操心积分怎么算，只需要把表达式填进正确的对话框。

3.5 算完怎么确认结果对不对

第一次算完，别急着看云图。先做一个“最小验证”：把对流和反应项全部设为 0，只保留扩散和源项，然后和一个你知道解析解的问题比一比。

比如二维拉普拉斯方程，矩形区域三边设 0，一边设 1，温度分布应该有经典傅里叶级数解。你的弱形式结果应该和内置“传热”模块算出来的结果几乎完全一致。如果不一致，先去检查边界条件是不是加错了，再去检查弱表达式里的正负号。这一步是每个用弱形式的人必经的“排雷”环节，能帮你快速建立信心。

4. 新手避坑指南：弱形式输入最常见的四个错误

4.1 坑一：搞不清 test(u) 到底是什么

现象：输入表达式时报错，提示test未定义，或者发现结果和预期完全不一样。

原因：test(u)表示试函数，也就是你在数学推导里写的那个v。它不是测量值，也不是测试版的意思。有些新手把它当成一个“检查按钮”，或者直接删掉，结果表达式变成了一个普通函数，软件根本不知道这是弱形式。

解决：记住test()是弱形式语言里的试函数标记。test(u)对应试函数本身，test(ux)对应试函数的 x 方向梯度。凡是弱表达式里出现的方程残差项，都要乘上对应变量的test()。没有test，系统无法完成“对任意试函数成立”的变分。

4.2 坑二：正负号跟着直觉走，结果差得很远

现象：表达式看着没问题，但解出来场量整体反了，或者出现明显不合理的振荡。

原因：弱形式的符号约定里有两种常见流派。一种把弱形式写成“∫ 弱项 = 0”，另一种写成“∫ 弱项 + ∫ 边界 = 0”。软件界面里的默认约定常常是“所有项移到等号左侧”，也就是说弱表达式里的每一项最终都是被积函数的一部分。你如果只记住了教科书里的推导步骤，直接把正号负号搬进去，很容易差一个全局负号。

解决：每输入一个新问题，先做“零载荷测试”：把源项设为 0，边界设为简单值，看是否得到物理上合理的平凡解（比如全常数）。再对照一个已知解析解验证符号。我现在的习惯是永远准备一个最简验证模型，符号对错一测便知。

4.3 坑三：忘了弱形式是积分表达式，直接填点值

现象：在弱表达式里写了u + 1，结果完全不起作用。

原因：弱形式在物理上是一个积分式。你输入的表达式是“被积函数”，不是“逐点约束”。比如要施加逐点约束u = 1，这不是弱表达式该做的事，而是需要用到“约束”框，或者通过添加罚项、拉格朗日乘子来实现。把点值塞进弱表达式里，只是给积分贡献了一项，不会强制任何点满足该值。

解决：分清两类输入。弱表达式负责“积分形式的方程”，约束框负责“强制点值”。如果一定要在弱表达式里施加约束，需要了解“弱约束”的用法，但初学者不要一上来就用。先用软件自带“约束”功能把你想要的边界固定住。

4.4 坑四：边界条件不知道放哪里

现象：计算结果边界不对，或者默认边界条件偷偷改变了你的设置。

原因：弱形式把边界项和区域项分离了。如果你只在区域弱表达式里输入了内容，边界弱表达式留空，软件默认边界是“自然边界条件”，通常意味着通量为零。这不是你想要的通量边界，反而可能对结果产生隐性约束。

解决：每次建模前，先明确你要的边界类型。第一类边界（固定值）放在“约束”框，第二类边界（通量/载荷）放在“边界弱表达式”。如果你不确定默认边界是什么，先右键查看节点的“边界条件”设置，别停留在默认值。很多莫名其妙的结果，都是因为边界项被遗漏。

5. 再往前一步：弱形式如何变成有限元求解

5.1 离散化：把无限自由度变成有限个节点值

弱形式写好了，接下来就是计算机真正干活的部分。你面对的方程仍然有无限多个可能解，计算机只能处理有限个未知数。于是有限元方法出场了。

它做的事很简单：把整个计算区域切成很多小单元，比如二维的三角形、三维的四面体。每个单元由节点连接而成。原来需要求一个连续的未知函数 u，现在只需要求每个节点上的数值。整个连续场被表示成这些节点值的组合，这就是“离散化”的含义。

弱形式在这里的优势是：它是一个积分表达式，每个小单元上的积分可以独立计算，再累加到一起。没有弱形式，你很难把 PDE 自然地切成小单元来处理。

5.2 基函数：像搭积木一样搭出场函数

离散化需要一个工具，叫基函数。每个节点对应一个基函数，它的形状像一个尖顶山包：在自身节点处值为 1，在相邻节点处值为 0，中间平滑过渡。用这些基函数做线性组合，就可以近似表达任意连续场。

弱形式里还有一个试函数。在前面推导里我们说过，试函数是任意可取的。在有限元离散中，我们就“任意”地让每个基函数轮流充当试函数。每代入一个试函数，弱形式就变成一个代数方程。所有节点都试一遍，就得到一组方程，写成矩阵形式就是：

K*U = L

U 是所有节点值组成的向量，K 是刚度矩阵，L 是载荷向量。

5.3 刚度矩阵：为什么里面这么多零

你可能听说过刚度矩阵是稀疏的。为什么？因为每个基函数只在它所属的那些单元里非零，在别处都是零。当试函数和第 j 个节点相关、解函数只在自己周围的节点非零时，积分只有在两者有重叠区域时才不为零。因此矩阵里大量元素是零。

稀疏矩阵非常利于存储和求解。软件内置的求解器专门针对这种结构做了优化。所以你不必担心矩阵太大，只要网格不是离谱地密集，通常都能跑得动。这也是为什么弱形式看起来复杂，但实际计算效率并不低。

5.4 几何解释：把解投影到有限维空间

如果非要给有限元一个几何直觉，可以说它是“投影”。所有可能的连续解构成一个无限维的函数空间，而你的基函数张成了一个有限维的子空间。有限元方法就是在这个子空间里寻找一个“最佳近似解”，使得弱形式残差相对于子空间里的所有试函数都为零。这个性质和“把一个向量投影到另一个方向上”很像，区别只是空间维度变高了。

当你理解了这一点，再看求解器里的高级设置就不会心虚了。比如为什么增加网格单元就可以提高精度？因为基函数张成的子空间更接近真实解所在的空间。为什么在某些对流主导的问题里会出现振荡？因为离散后刚度矩阵不对称，近似解可能会在局部出现过冲。这些现象都不是玄学，而是弱形式和有限元这两个概念共同作用的结果。

6. 下一步学什么：三条路径帮你继续深入

如果你已经跑通了最简扩散模型的弱形式，接下来可以按三条路径继续深化。下面的表格可以帮你快速定位：

入门的话建议先做路径一，因为结构力学里的能量最小化与弱形式天然对应，你能更直观地理解试函数。具体做法是：找一道一维杆拉伸的题目，写出它的弹性能泛函，再手动求一阶变分，对比软件里的弱表达式，你会瞬间明白许多过去死记硬背的公式。

路径二适合你已经敢改方程。可以把时间项du/dt加进弱形式，COMSOL 里用d(u,t)*test(u)表示。注意，加了时间项之后，你需要选择“瞬态研究”，而不是“稳态”。很多初学者卡在这一步，其实只要在“研究”里切换类型就行。

路径三更像是一种高阶修炼。等你见过几个非线性不收敛的报错后，就会懂得看求解器日志。那个看似黑匣子的“非线性方法”其实是牛顿迭代，每一次迭代都在求解一个类似 K*U=L 的线性系统。

从那以后我每次新建弱形式模型，都强制自己先做一个最简验证：先关掉所有复杂项，只留一项，算一个能手动验算的结果，然后再逐步打开其余项。这个习惯替我避开了至少一半的符号错误和边界设置问题。希望帮到你。

本文还有配套的精品资源，点击获取