5.2 能量原理与近似解
5.2 能量原理与近似解
本节摘要：能量原理换一条路求解：虚功原理把平衡写成功等式，最小势能原理把求解化成泛函极值，瑞利-里茨法再用少数试探函数把极值问题缩成线性代数。本节完成从原理到近似的全链条演算，验证近似挠度必然偏小的上界性质，并交代它与第九章有限元的传承关系。
换一条路：不解方程，改求极值
5.1 节的应力函数走的是精确路线：把平衡与相容打包进一个函数，解双调和方程。这条路对规整几何有效，一旦几何复杂就断路。本节换思路：彻底放弃"逐点满足方程"，改为在全部许可位移场里挑一个使总势能最小的。真实平衡位移场使总势能取最小值——这是最小势能原理的内容，也是本章从"微分方程"跳到"泛函极值"的那一步。极值问题天然适合近似：试探函数选得粗，误差大但可控；选得细，逐项逼近真解。第九章有限元正是把这里的思想离散化，每个单元的形函数就是一组局部试探函数。
虚功与最小势能：一分钟的推导
一切从虚功原理开始：平衡力系在任意虚位移上做的总功为零。对弹性体，外力虚功等于应变能的变分，δW 外 = δU。把它整理成变分形式， δ(U − W 外) = 0，括号里的量正是总势能 Π = U − W 外。于是平衡条件等价于总势能驻值；再进一步可以证明，在线弹性稳定问题里这个驻值是最小值——任何偏离平衡的许可位移都让势能上升。这就是"解方程"与"求极值"两件事等价的原因。判读时要抓住一个细节：许可位移场只需满足位移边界条件（几何许可），不必预先满足平衡——平衡恰恰是极值条件替你自动完成的。
图：能量原理谱系与近似解路线
瑞利-里茨法：试探函数与单参数演算
方法本体只有三步：把位移场写成试探函数的线性组合 c₁φ₁ + c₂φ₂ + ……，其中 φᵢ 事先选好、只满足几何边界条件；代入总势能 Π(c₁, c₂, ……)；对每个 cᵢ 求偏导置零，得到线性方程组。试探函数一旦选定，无穷维的泛函极值就缩成了有限个未知量的代数问题——缩得多狠，取决于你敢用几个参数。
拿标准试件验一遍。简支梁跨中受集中载荷 P，抗弯刚度 EI，跨长 L。选单参数正弦试探挠度 w = c·sin(πx/L)——它在两端挠度为零、弯矩为零，几何条件全中。应变能 U = EI/2 乘 w'' 平方的积分，积分限零到 L；外力功 W = P 乘跨中挠度 = Pc。代入 Π = U − W，对 c 求导置零，解出 c = 2PL³ 除以 π⁴EI。四阶导为零是正弦函数的天赋，这一步用不着数值积分。
# 瑞利-里茨单参数解：简支梁跨中集中载荷 import math P, L, EI = 10e3, 6.0, 4.5e6 # 载荷 N、跨长 m、抗弯刚度 N·m2 c_ritz = 2 * P * L**3 / (math.pi**4 * EI) # 里茨系数 w_exact = P * L**3 / (48 * EI) # 材料力学精确解 print(f"里茨挠度 = {c_ritz*1e3:.3f} mm") print(f"精确挠度 = {w_exact*1e3:.3f} mm") print(f"偏差 = {(c_ritz/w_exact - 1)*100:+.2f} %") # 应变能对照：里茨应变能应低于精确值（刚度偏大的签名） U_ritz = EI * math.pi**4 * c_ritz**2 / (4 * L**3) U_exact = P**2 * L**3 / (96 * EI) print(f"应变能 里茨 = {U_ritz:8.2f} J，精确 = {U_exact:8.2f} J，" f"比值 {U_ritz/U_exact:.4f}")
里茨挠度 = 9.855 mm 精确挠度 = 10.000 mm 偏差 = -1.44 % 应变能 里茨 = 49.28 J，精确 = 50.00 J，比值 0.9856
判读这份记录：单参数近似就把跨中挠度算到百分之一点五以内的误差里，代价只是一次求导和一次积分。更值得注意的是偏差的方向——挠度偏小，应变能同样偏低约百分之一点五。这不是巧合而是定理：位移型里茨解相当于给结构额外加了约束，结构变"硬"，真实解的挠度与应变能是它的上界。工程判读里这是免费的保险：用里茨挠度做变形验算，结果只会偏安全一侧。
逐项加料：收敛快得像付清全款
单参数解给了方向感，多参数解给精度。正弦试探函数族的妙处在于简支梁上彼此"不正相交也相容"：跨中集中载荷只激发奇数半波，各级系数可以逐个独立解出，挠度级数是 2PL³ 除以 π⁴EI 乘 (1 + 1/81 + 1/625 + 1/2401 + ……)。四阶衰减的级数收敛极快——第二项只补百分之一点二，第三项只补千分之一点六。
# 逐项加试探函数：跨中挠度级数收敛过程 import math P, L, EI = 10e3, 6.0, 4.5e6 w_exact = P * L**3 / (48 * EI) coef = 2 * P * L**3 / (math.pi**4 * EI) total = 0.0 print("项数 n 跨中挠度 mm 相对误差") for n in range(1, 5): total += coef / (2*n - 1)**4 # 奇数半波项 1 3 5 7 err = (total/w_exact - 1) * 100 print(f" {n} {total*1e3:9.3f} {err:+8.4f} %")
项数 n 跨中挠度 mm 相对误差 1 9.855 -1.4445 % 2 9.977 -0.2278 % 3 9.993 -0.0722 % 4 9.997 -0.0306 %
判读：四项之后误差只剩万分之三点一，而每一项的边际贡献缩到前一项的三分之一以下——1/n⁴ 衰减是这类光滑试探函数的典型签名。实践中根本用不到这么多项：变形估算取两项，概念校核取一项，足够。真正记住的应该是模式：近似解的误差单调地从一侧逼近真解，看到误差符号不变，就可以放心停机。
从里茨到有限元：一脉相承
瑞利-里茨法有个先天短板：全局试探函数要求"一段光滑函数逼近整个结构"，几何一复杂就写不出来。有限元的解法朴素而致命——把结构切成小片，每片内用最简单的多项式试探函数（杆用一次、梁用三次），片与片之间只在节点上握手。全局函数写不了，分段多项式永远写得出；参数从"几个全局系数"变成"所有节点的位移"。第五章开篇那句"理解本章是给数值方法打地基"在这里兑现：单元刚度矩阵就是单元总势能对节点位移求导的产物，整体刚度矩阵的组装就是总势能求和。第九章动手装配时，请回头翻这条谱系——你做的每一个矩阵运算，都是本节泛函极值的代数化身。
本节要点回顾
- 虚功原理给出平衡的功等价形式，整理成 δΠ 等于零就得到最小势能原理。
- 里茨法三步：选试探函数、代入总势能、对系数求导，泛函极值缩成线性代数。
- 位移型近似解刚度偏大：挠度、应变能都低于真解，是偏安全方向的上界。
- 奇数半波级数按 1/n⁴ 收敛，四项误差即到万分之三。
- 有限元是里茨思想的分段版本：单元刚度矩阵来自单元势能求导。
解析台的干货到此搬空。构件才是实验室的常客——下一章进入构件试验段。