有限元理论中散度定理的物理意义
在有限元理论相关的著作中，当涉及变分运算时，常常出现散度定理
或者矩阵形式
式中符号说明：
• ：体积 • ：表面积 • ：应力分量 • ：位移分量 • ：边界法向量的方向余弦 • ：应变矩阵（微分算子矩阵） • 、 ：应力和位移的向量表示 • 下标逗号表示偏导数（如 ）
散度定理的物理意义可以从“通量与源”的关系来理解。它的核心思想是：
一个封闭曲面上的通量，等于该曲面所包围的体积内所有源的代数和。
具体到有限元分析中，它起到了连接“边界积分”与“体积积分”的桥梁作用，这是推导弱形式和控制方程离散化的关键工具。
1. 从数学定义看
对于一个向量场 ，其散度定义为：
它就是一个标量，表示每一点上向量场的“源”强度。正值表示向外发散（有源），负值表示向内汇聚（有汇）。
例 1 对于向量场 有
这个向量场如图1所示。在 中的任何一点，尾部在该点的箭头比头部在该点的箭头更长。因此，从每一点流出的流体多于流入的；也就是说， 在在每一点都在"发散"。（因此，我们看到了"散度"一词的由来。）
▲图1
例 2 对于向量场 指向与图1中向量场 相反 在的方向（见图2），并且应该很清楚 的散度 是如何在图中体现出来的。
▲图2
例 3 常数向量场 ， 处处成立，则称该向量场是不可压缩的。直观上，我们可以看到 的每一点都有一条尾部在该点的箭头表示 ，还有另一条头部在该点的箭头也表示 ，这样在任何一点处的向内和向外流量都是相同的（见图3）。
▲图3
2. 在弹性力学中的核心应用：平衡方程
弹性力学里，我们最常处理的是二阶张量——应力张量 。对张量求散度，就是对它的每一行（或每一列）分别求散度，得到一个向量。
最关键的关系就是平衡方程（无惯性项时的静力平衡）：
这里：
• 是应力张量的散度，代表单位体积内由应力引起的净力。 • 是单位体积上的体力（如重力）。
这个方程说的就是：某一点上，应力散度产生的力与外力相平衡。
直观理解：
想象物体内部的一个微小立方体。六个面上都有应力作用。如果某个方向上的应力随位置变化（比如左面压力小、右面压力大），那么这个微元体就会在该方向上受到一个净力。这个净力的大小，正是应力在该方向上的偏导数之和，也就是散度。
▲图4
所以，应力散度 ≠ 0 的地方，意味着存在不平衡的内力，必须由外力来补偿，否则物体会加速运动。
可以把应力散度想象成“力的净流出率”：
• 如果一个点的应力散度为正，说明从这个微元体“流出”的力大于“流入”的力，相当于该点有一个向外的“力源”，需要施加一个向内的体力才能平衡。 • 如果为负，则相反，需要施加向外的体力。
所以，当你看到 时，可以把它理解为：“看看这一点周围，应力是怎么变化的，这些变化加起来会产生多大的净力”。
3. 在有限元中的作用：从强形式到弱形式
有限元方法通常从控制方程的微分形式（强形式）出发，例如：
但直接求解这个偏微分方程很困难。利用高斯公式，可以将其转化为积分形式（弱形式）：
1. 将方程乘以一个测试函数（权函数）并在整个域上积分。 2. 利用高斯公式将包含散度的项从体积积分转化为边界积分。 3. 得到弱形式：
其中 是边界上的面力。
这一步骤的物理意义：
• 原本需要满足每一点都平衡（强形式），现在只要求在加权平均意义上平衡（弱形式）。 • 高斯公式将“内部散度”转化为“边界通量”，使得自然边界条件（如给定面力）自动进入方程，无需额外强制施加。
在有限元中，高斯公式的物理本质是：
将局部守恒律（微分方程）转化为全局守恒律（积分方程），从而允许用分片多项式近似来求解连续介质问题。
