3. 能量原理与变分法：虚功原理、最小势能原理、变分法在有限元中的应用

各位同学，今天我们来聊聊有限元法的“灵魂”——能量原理与变分法。说实话，我刚接触这部分内容时，也觉得它太抽象了。什么虚功、什么泛函，听着就头大。但后来我发现，不理解这些原理，你永远只是个“软件操作工”。

你想想看，有限元法本质上在干什么？它是在找一个“最接近真实”的近似解。那怎么判断“最接近”？靠的就是能量原理。说白了，大自然很“懒”，它总是选择能量最小的那个状态。我们有限元就是模拟这个“偷懒”的过程。

3.1 虚功原理：力学世界的“平衡天平”

虚功原理，我个人觉得是理解有限元最直观的入口。它讲的是：如果一个结构处于平衡状态，那么对于任意微小的、满足边界条件的虚位移，外力做的虚功等于内力做的虚功。

用公式表达就是：

δW_ext = δW_int

其中，δW_ext 是外力虚功，δW_int 是内力虚功（也就是虚应变能）。

我在项目中遇到过这样一个案例：一个大型桥梁的吊杆，用传统方法算受力，怎么算都跟实测对不上。后来我改用虚功原理去推导，才发现是边界条件处理错了。嗯，这里要注意，虚功原理对边界条件特别敏感。

3.2 最小势能原理：大自然的选择

最小势能原理是虚功原理的“好兄弟”。它说：在所有可能的位移中，真实位移使系统的总势能取最小值。

总势能 Π 可以写成：

Π = U - W

U 是应变能，W 是外力势能。真实解就是让 Π 对位移的变分为零：

δΠ = 0

为什么会这样？你想想看，一根悬臂梁，你用手压它，它总会停在一个稳定的位置。这个位置就是势能最小的位置。我曾经用这个原理去验证一个复杂结构的有限元结果，发现只要总势能曲线是凹的，结果就一定是收敛的。

3.3 变分法：寻找“最优”的数学工具

变分法，说白了就是“微积分”的升级版。微积分找的是函数的最值，变分法找的是泛函的最值。泛函是什么？简单说，就是“函数的函数”。

在有限元中，我们经常要解这样的问题：

求 u(x)，使得 Π[u] = ∫ F(x, u, u') dx 取最小值

这个问题的解必须满足欧拉-拉格朗日方程：

∂F/∂u - d/dx (∂F/∂u') = 0

我记得刚开始学变分法时，总觉得它跟有限元没什么关系。直到有一次做热-力耦合分析，发现直接解偏微分方程根本解不出来。后来改用变分法，把问题转化成求泛函极值，再用有限元离散，一下子就通了。

3.4 变分法在有限元中的具体应用

这部分是重点，我建议大家认真看。有限元的核心步骤，其实就是把变分法“离散化”。

具体流程是这样的：

1. 建立泛函：根据物理问题，写出总势能 Π
2. 假设位移模式：用形函数 N 和节点位移 u_e 表示单元内位移：u = N u_e
3. 代入泛函：把 u 代入 Π，得到关于 u_e 的表达式
4. 求变分：令 δΠ = 0，得到单元刚度方程：K_e u_e = f_e
5. 组装求解：组装全局刚度矩阵，求解节点位移

举个例子，一维杆单元的推导：

应变能 U = 1/2 ∫ E A (du/dx)² dx
外力势能 W = ∫ f u dx + P u|x=L

代入 u = N1 u1 + N2 u2
得到 K_e = (EA/L) [1 -1; -1 1]

你看，这个经典的单元刚度矩阵，就是从变分法一步步推出来的。我个人习惯在推导新单元时，永远从变分法开始，而不是直接抄书上的公式。这样即使遇到复杂的耦合问题，心里也有底。

3.5 知识体系框架

下面这张图，是我自己总结的本章知识体系。它把虚功原理、最小势能原理和变分法的关系，以及它们在有限元中的位置，都串起来了。

3.6 本章小结

好了，这一章的内容就到这里。我最后再啰嗦几句：

• 虚功原理是“平衡”的另一种表达，它不依赖材料本构，适用范围最广
• 最小势能原理是虚功原理在弹性问题中的特例，它给出了“最优解”的判据
• 变分法是连接物理原理和数值计算的桥梁，没有它，有限元就无从谈起

我个人觉得，这三者就像“三兄弟”：虚功原理是大哥，告诉你“平衡”是什么；最小势能原理是二哥，告诉你“最优”是什么；变分法是三弟，负责把前两位的想法“算出来”。

下次当你用有限元软件算出一个结果时，不妨想想：这个结果是不是让总势能最小了？虚功是不是平衡了？如果答案是肯定的，那你的模型大概率是对的。

---

公众号:蓝海资料掘金营，微信deep3321