第 12 章　能量法初步

12.1应变能与功能原理

拉伸一根弹簧，你做了功；松手，弹簧弹回来——它把你做的功「存」了起来。弹性体也是这样：外力对构件做功 W，这部分功以应变能（弹性变形能）U的形式储存在构件内部。在小变形、缓慢加载（无能量耗散）的条件下，两者相等：

这就是功能原理，能量法的起点。它的妙处在于：我们不直接去解「力-位移」的复杂关系，而是先算能量，再从能量里「读出」位移。这就像「知道花了多少力气（能量），反推物体动了多远（位移）」。

对线弹性体，力与位移成正比（胡克定律），外力做功等于力-位移图线下三角形的面积：W = P·δ/2。这正是应变能里「1/2」系数的由来。

「1/2」的来历：克拉贝隆定理

外力做功为什么要乘 1/2？因为力是「从零慢慢长上去」的。缓慢加载时，外力 P 与其作用点的位移 δ 同步从零增长到终值，且在线弹性范围内两者成正比 P=kδ。加载过程中力做的功等于力-位移曲线下的面积——一条过原点的直线下的三角形：

这个结论叫克拉贝隆定理（Clapeyron's theorem）：线弹性体在缓慢加载下，外力功等于「终态力 × 终态位移」的一半。若荷载一步到位地「砸」上去（突然加载、冲击），外力功是 P·δ 全量，其中只有一半变成应变能，另一半变成结构的动能与振荡——这就是冲击问题比静力问题危险得多的能量根源，也是第 8 章动荷系数的能量法出处。能量法并不只求位移，它还是分析冲击与振动载荷的必经之路。

应变能存在哪里：微观图景

应变能并不是一个抽象的记账符号，它物理地储存在被拉伸的原子键里。杆件受拉时，晶格内原子间距离被拉开，原子间的结合力（本质上是一种保守力）对抗外力做功，把功转化为势能储存起来；卸载时这些势能重新释放，宏观上表现为杆件回弹。这跟拉长一根弹簧、松手弹回是完全相同的机制——整个构件就是天文数字个「微弹簧」的并联体。这个图景还顺带解释了两件事：其一，弹性模量 E 反映的是原子键的刚度，所以它几乎不随合金成分、热处理状态改变（各种钢材的 E 都在 200GPa 上下）；其二，应变能恒为正——无论拉压，键的势能都升高，「内力平方」的结构正是这一物理事实的数学体现。

由此还能引出应变能密度（比能）的算法：单位体积内储存的应变能 u=σ²/(2E)（拉压情形）。把比能对体积积分，就得到构件的应变能——这正是下一节三种基本变形应变能公式统一的推导源头。

功能原理的成立条件

辨析　功能原理 U=W 并不是无条件成立：它要求加载缓慢（构件各处加速度可忽略、无宏观动能）且无耗散（无摩擦热、无塑性功）。工程上把满足这两条的加载称为「准静态加载」。一旦加载很快（冲击、爆炸）或有内耗（塑性、蠕变），能量账本就要分册记账：动能、耗散热都得单列。本章一切结论都建立在准静态、线弹性两个前提上——这与卡氏定理、单位荷载法的适用范围完全一致。

工程落点　应变能思维在设计中的两个经典去处：一是缓冲吸能——汽车的溃缩吸能结构、直升机可折叠起落架，本质上是「设计一个用塑性变形或结构压溃来大量消耗应变能的区域」，把冲击能量挡在乘员与关键设备之外；二是刚度分配——串联结构的总柔度等于各段柔度之和，应变能表达式里「除以刚度」的结构直接告诉你：软的那一段变形最大、储能也最多，想提高整体刚度就要优先加硬「最软的环节」。

刚度的能量语言

把 U=Pδ/2 与线弹性关系 P=kδ 联立，还能挖出一个漂亮的恒等式：U=kδ²/2=P²/(2k)。反过来，k=2U/δ²=P²/(2U)——刚度可以从能量里读出来。以例题 1 的直杆验证：k=P/δ=40000/0.955×10⁻³≈41.9MN/m；而用 k=2U/δ²=2×19.1/(0.955×10⁻³)²≈41.9MN/m，两算一致。这层等价关系带来一个系统级的直觉：串联结构（力路依次通过各段）的柔度等于各段柔度之和——因为总应变能是各段之和，而各段的「P²/(2k)」按 1/k 累加；并联结构则反过来刚度相加。结构里的「力流」概念，用应变能的语言说得最干净：力流偏爱刚硬的路径（储能少），薄弱的路径以大变形、高储能的方式「分担」——这也是评估传力路径合理性时，工程师爱看应变能分布图的原因。

与第 9 章的握手　比能公式 u=σ²/(2E) 描述的是单向应力状态。推广到复杂应力状态时，应变能密度等于各应力分量贡献的叠加；更妙的分解是把它拆成「体积改变比能」与「形状改变比能（畸变能）」两部分——前者对应单元在静水应力下的均匀胀缩，后者对应单元体的剪切畸变。第 9 章第四强度理论的物理内核「畸变能密度达到单向拉伸屈服时的畸变能即屈服」，正是用这把能量尺子量出来的。能量语言与应力状态语言在此握手：一个从「储存了多少」看问题，一个从「哪个面上多大」看问题，两边说的是同一件事。

12.2三种基本变形的应变能

把前面几章的公式代入「应变能密度」积分，可得到三种基本变形下构件的应变能：

这三个公式形式惊人地统一，可记为「内力² × 长度 / (2 × 刚度)」。内力 N、T、M 分别对应刚度 EA、GI_p、EI。记忆规律：内力的平方、除以两倍对应刚度、再乘长度。

若内力沿杆长变化（如弯曲时 M 是 x 的函数），则需对杆长积分：

三个公式是同一个推导的三次出口

表格里三个公式的统一形式不是巧合，它们都从比能积分这一个源头流出。以拉压为例：杆内各点应力均匀，比能 u=σ²/(2E)=N²/(2E²A²)，乘以体积 V=Al，得 U=N²l/(2EA)。扭转完全同理，把 σ 换成最大切应力对应的比能 u=τ_max²/(4G)……更严格地沿半径积分后得到 U=T²l/(2GI_p)。弯曲时弯矩沿杆长变化，取微段 dx：微段的应变能 dU=M²(x)dx/(2EI)，积分便是 U=∫M²/(2EI)dx。三个出口，一个源头：比能 × 体积，再对不均匀的内力场积分。

「内力的平方」还有一层深意：应变能与内力的符号无关、恒为正。拉力与压力储存同样多的能量，正弯矩与负弯矩也各记各的账、互不抵消。这使应变能成为一个纯粹的「标量账本」——不关心方向、只关心积累，正式这份「不挑方向」的性质，让能量法绕开了矢量力学的方向纠缠，成为处理复杂结构的第一选择。

应变能不可叠加：能量法的第一大坑

一个极易踩中的误区：既然三根杆的应变能可以相加，那两个荷载的应变能也能相加吗？不能。看轴向受拉的展开式——设两荷载单独引起的轴力为 N₁、N₂，共同作用时总轴力是 N₁+N₂，于是：

交叉项 N₁N₂l/(EA) 是「两荷载互相在对方的位移上做功」——只要两荷载引起的内力不同时为零，它就存在。所以「先算各自再相加」必然漏掉这一项。应变能是状态的函数，不是荷载的线性账本；只有「不同变形形式」之间（如一根轴同时受拉与受扭，拉与扭的内力互不交叉）才能直接叠加。这条辨析在用卡氏定理处理多荷载结构时性命攸关：必须先把总应变能用全部荷载写出，再求偏导，绝不许「分荷载算能量再相加」。

从应变能看弯矩图　U=∫M²/(2EI)dx 说明一件事：梁内应变能正比于弯矩图的「平方面积」。平方运算会放大峰值、削弱小值区——弯矩图上 100 与 50 的两段，储能贡献是 4:1 而不是 2:1。这条几何直觉直接指向两条设计经验：其一，削峰最划算——把弯矩峰值从 Pl/4 降到 Pl/6（例如连续梁把支座处负弯矩「借」给跨中正弯矩），总应变能与挠度都显著下降；其二，变截面梁（弯矩大的区段加高）之所以省料，本质是让 M²/(EI) 沿全梁趋于均匀——能量账本上每一毫米都吃得一样多，材料才算物尽其用。

被忽略的剪切应变能　细心的读者会问：梁弯曲时剪力 F_S 也做功，为什么表中不储？答案是量级。以跨中受力的简支矩形梁为例，剪切应变能与弯曲应变能之比约为 1.2×(E/G)×(h/l)²：钢梁取 E/G≈2.6，则 h/l=1/10 的细长梁只有约 3%，忽略完全正当；h/l=1/5 时升至约 12%；h/l=1/3 的深梁可达约 35%，剪切变形就必须与弯曲变形一起计入（短跨木梁、深梁腹板属此类）。第 8 章挠度公式默认「细长」、本章应变能公式同样默认「细长」，两条线索在同一个近似上会师——学到这里，你对「每个公式的适用范围」应当已形成条件反射。

12.3互等定理

能量法引出了两个优美而实用的定理：

• 功的互等定理：第一组力在第二组力引起的位移上做的功，等于第二组力在第一组力引起的位移上做的功。
• 位移互等定理：在点 1 作用单位力引起点 2 的位移，等于在点 2 作用单位力引起点 1 的位移（即 δ₁₂ = δ₂₁）。

位移互等定理听起来很「反直觉」——为什么 A 点受力、B 点会位移，反过来 B 点受力、A 点的位移居然相等？但它的本质很简单：因为应变能只取决于最终的受力状态，而与加载的先后顺序无关。两种不同顺序的加载，最终应变能相同，比较两者的功，就得到互等关系。

功的互等定理：两种加载顺序的对账

互等定理的证明只需要一个事实：应变能是最终状态的函数，与加载顺序无关。在弹性体上先加力 F₁（于 1 点）、再加力 F₂（于 2 点）：

• 顺序一（先 F₁ 后 F₂）：外力功 = F₁δ₁₁/2 + F₂δ₂₂/2 + F₁δ₁₂。前两项是各自加载的「三角形功」，第三项是 F₁ 在 F₂ 引起的位移 δ₁₂ 上「全程满力」做的功；
• 顺序二（先 F₂ 后 F₁）：外力功 = F₁δ₁₁/2 + F₂δ₂₂/2 + F₂δ₂₁。

两个顺序的终态完全相同，储存的应变能必须相等，于是交叉项必须相等：

这就是功的互等定理：第一组力在第二组力引起的位移上做的功 = 第二组力在第一组力引起的位移上做的功。取特例 F₁=F₂=1（都是单位力），立即得到位移互等定理 δ₁₂=δ₂₁。注意这里 δ 的第一个下标指位移发生的位置、第二个下标指力的作用位置；「力」应作广义理解——集中力对应线位移，力偶对应角位移（相应地也有「力偶×转角」的互等）。

工程落点　互等定理是结构试验的「换位神器」：桥梁静载试验中，加载车停靠某些位置不便、而测点可以灵活布设——你可以在方便加载的位置加力、去读关心位置的位移，再把结果换算到「荷载作用在关心位置」的工况。用互等定理解释：两种工况的位移响应一一互等，试验数据因此一张表顶两张表用。机身与机翼的全机静力试验中，加载点与测量点的布置同样大量运用这条定理来减少加载作动筒的数量。

反力互等定理　互等家族还有第三位成员：支座 1 发生单位位移引起的支座 2 的反力，等于支座 2 发生单位位移引起的支座 1 的反力（R₂₁=R₁₂）。它同样源自「能量与路径无关」。工程价值在支座沉降分析：连续梁某一中间支座下沉 δ，各支座因此产生的附加反力，可以用「单位沉降」工况的反力按比例放大得到，而单位工况又可以通过互等关系从「支座处加单位力」的位移工况反查——一套表格覆盖全部沉降组合，是桥梁与基础工程分析的常用套路。

适用边界　互等定理成立的前提与整个能量法一致：线弹性、小变形、准静态。它对「力」没有形态限制——力与位移都可以是广义的（力偶↔转角、一对反向力↔相对位移），互等关系照样成立；但两组力必须作用在同一个弹性体系上（尺寸、约束、材料完全一致），拿「加撑前」与「加撑后」两个体系去谈互等是初学最常见的误解——体系一变，系数全变，对账就再也对不上了。

12.4卡氏定理

能量法最有用的成果是卡氏第二定理（简称卡氏定理）：

意思是：构件某点沿某力方向的位移 δ，等于总应变能对该力求偏导数。这是个「一键求位移」的公式：

1. 写出构件的总应变能 U（作为各外力的函数）；
2. 对「想求位移的那个力」求偏导；
3. 代入实际力值，即得该点该方向的位移。

以弯曲为主的梁为例，代入应变能积分式，卡氏定理变成实用的积分形式：

一个实用技巧：如果要求的点处没有真实外力，可以先假想地在那里加一个「虚拟力」，求完偏导后令虚拟力=0。这个技巧把卡氏定理的适用范围大大扩展了。

卡氏定理为什么成立：一步微扰看穿本质

定理的证明只需要「微扰」一步。设构件在力 P₁、…、P_n 作用下储存应变能 U（状态的函数）。保持其他力不变，只给 P_i 一个微增量 dP_i：外力功的增量按两种口径计算——

• 口径一（能量角度）：dW=dU，而 U 是各力的函数，dU=(∂U/∂P_i)·dP_i；
• 口径二（直接做功）：dP_i 是从零长上去的微量，其他力「全程满力」，但它们作用点的位移没变；只有 dP_i 在自己方向位移 δ_i 上做功，dW=δ_i·dP_i（微量自乘项是二阶小量，略去）。

两个口径必须相等：δ_idP_i=(∂U/∂P_i)dP_i，即 δ_i=∂U/∂P_i。证毕。这一步微扰还顺带解释了求导的「方向含义」：对力求偏导时其他力保持不变，所以多力场合绝不能漏写任何一项内力。求转角时道理相同——对作用在该截面的力偶求偏导即得转角，θ=∂U/∂M₀。

顺带认识它的姊妹命题卡氏第一定理：若把 U 写成位移的函数，则 P_i=∂U/∂δ_i。第二定理「对力求导得位移」，第一定理「对位移求导得力」。第二定理要求线弹性（否则 W≠U 的关系被破坏），第一定理连非线性弹性也成立——两条定理在现代非线性结构分析与有限元里各司其职。

卡氏定理打开超静定之门

能量法真正的威力在超静定结构上。回顾第 8 章：超静定结构的多余约束使静力方程解不出全部未知力——还差「变形协调条件」。用卡氏定理，解题流程机械得几乎可以流水线化：

1. 解除多余约束（如支座 B），代之以未知的约束反力 R——结构化为静定基本体系，内力全部写成「已知荷载 + R」的函数；
2. 写出基本体系的应变能 U（荷载与 R 的函数）；
3. 原结构在 B 处的位移是零（或已知值），用卡氏定理翻译成方程：∂U/∂R=0（或=已知位移）；
4. 解方程得 R，再按静定体系算全部内力与位移。

这里的妙处在于：解除约束后写出的内力自动满足全部平衡条件，卡氏定理的偏导又自动给出变形协调条件——两大类方程各就各位，不需要像力法手工拆解那样反复调整未知量的组合。一次超静定补一个方程，n 次超静定联立 n 个偏导方程；方程系数都是「单位力在彼此位移上做的功」，正是 12.3 节互等定理里的那些量——三节内容在这里会师。后续《结构力学》的力法、位移法，乃至大型有限元软件求解器里的柔度/刚度方程，走的都是这条能量路线。

卡氏定理操作清单　用卡氏定理求位移，把四条纪律背成肌肉记忆：① U 必须包含全部荷载（含待求位移之外的一切外力与多余约束力），否则交叉项漏项；② 「对力求导」中的力是广义力——求线位移对集中力求导、求转角对力偶求导；③ 多个同向力 P 作用在同一结构上时，不能对「名义上相同的 P」求导（它们是独立变量，须先区分 P₁、P₂），或改用单位荷载法避坑；④ 先求偏导、后代入数值——附加力法里「令 F=0」永远是最后一步。四条纪律条条都在防同一个错误：把「状态的函数」当成「荷载的线性账本」。

12.5单位荷载法

与卡氏定理等价、但更灵活的另一种方法是单位荷载法（莫尔定理）：

其中 M 是真实荷载下的弯矩，M̄ 是在「所求位移处、沿所求位移方向」作用一个单位力产生的弯矩。思路非常直观：

1. 在要求位移的点，加一个单位力（方向就是要求位移的方向）；
2. 分别算真实荷载下的弯矩 M 和单位力下的弯矩 M̄；
3. 两者相乘、除以 EI、沿杆长积分，结果就是该点的位移。

单位荷载法：卡氏定理的「标准化改造」

卡氏定理对每个待求位移都要重写一次 U 并求一次偏导，工程师们很快发现：对线弹性结构，「∂M/∂P」这一步可以标准化——在待求位移的位置沿待求方向加单位力，它单独引起的弯矩正是 M̄，且有恒等式 ∂M/∂P=M̄。于是卡氏定理改写成 δ=∫M M̄/(EI)dx，这就是单位荷载法（莫尔定理）。它把「求导」这道工序替换成「另算一个单位力状态」——两套弯矩各自独立、互不干扰，特别适合多工况、多测点的批量计算，也天然适合编程：M 与 M̄ 不过是两个函数表，相乘积分即可。图乘法（结构力学中把两个弯矩图「面积×形心对应值」相乘以代替积分）是它在手算场合的加速器，此处按下不表。

三个常用变体一并记住：求转角时加单位力偶；求两点相对位移时加一对反向单位力（如桁架某杆两端，得杆长变化）；求相对转角时加一对反向单位力偶。单位荷载从不改变结构——它只是一个「度量工具」，加在哪个自由度上，就量出哪个自由度的位移。

通往有限元的血统　单位荷载法把「求一个位移」组织成「两个状态的函数相乘再积分」。把结构离散成许多小单元、把单位力状态换成「形函数」，这套「两个场相乘积分」的模式就长成了有限元方法中的虚功方程与刚度矩阵组装——现代结构分析软件每一次求解，内核里都在批量执行本章 12.5 节的运算。可以说，能量法是手算时代的高峰，也是机算时代的起点：学会本章，你不但掌握了工程师的「一把尺」，还提前摸到了数字仿真的门把手。

三法选用速查　功能原理直接算能量，适合单一荷载、求储能或校核；卡氏定理适合「几个力同时作用、要其中某个力方向的位移」，特别是超静定问题（对约束反力求偏导=变形协调）；单位荷载法最通用——任意点、任意方向（含转角、相对位移），公式形态统一、便于编程。三者共享同一个地基（U=W 与应变能状态函数性质），解题时按「结构简单程度 + 待求量类型 + 是否编程」三问选刀即可。

桁架形式：单位荷载法的求和版

桁架杆件只受轴力，单位荷载法退化为一个简单的求和：

其中 N 是真实荷载下各杆轴力、N̄ 是单位力状态下各杆轴力、l 为杆长。逐杆列表计算、求和即得位移——航空与桥梁桁架的手算校核至今仍用它，因为一张表格就能让两位工程师独立复算对表。

给全书收个尾　回望十二章：内力与应力是「看清构件在承受什么」，强度理论是「判断它何时破坏」，刚度与变形计算是「量出它变形多少」，压杆稳定是「确认它站得住」，而能量法把这一切重新装订成一本「能量的账」。两百多年前欧拉用变分的眼光看平衡，克拉贝隆把功写成三角形，卡斯蒂利亚诺把能量变成可求导的函数，莫尔把求导标准化成两图相乘——今天每一座桥梁、每一架飞机的结构分析软件里，依然运行着这几条定理的代码化身。掌握了这一章，你手里已经不只是材料力学的知识点清单，而是一整套「从力学现象到数学结构、再到工程判断」的思维方式——带着它去学结构力学、有限元与任何一门工程分析课程，都会事半功倍。