Maxwell-Betti 互等定理模拟器 返回
结构力学模拟器

Maxwell-Betti 互等定理模拟器

在简支梁的两组荷载下验证 δ_AB·P_A = δ_BA·P_B 严格成立。直观体会影响函数 C(a,b) 的对称性如何对应线性弹性以及 FEM 刚度矩阵的对称性。

参数设置
梁长 L
m
点 a 位置
m
点 b 位置
m
荷载 P_A
kN

荷载 P_B 固定为 50 kN,弯曲刚度 EI 固定为 1×10⁴ kN·m²。a 与 b 的上限按 L−0.1 m 自动限制。

计算结果
a 点 B 系挠度 δ_BA (mm)
b 点 A 系挠度 δ_AB (mm)
δ_BA·P_A(左侧)(kN·m)
δ_AB·P_B(右侧)(kN·m)
系统 A 与系统 B 的挠度曲线

上:系统 A — a 点作用 P_A,蓝色为挠度曲线,黄色标记为 b 点的 δ_AB;下:系统 B — b 点作用 P_B,红色挠度曲线,黄色标记为 a 点的 δ_BA。

理论与主要公式

等截面简支梁(长度 $L$,弯曲刚度 $EI$ 恒定),在 $x = a$ 处作用集中荷载 $P$ 时,$x = b$ 处的影响函数 $C(a,b)$ 为:

当 $b \le a$:

$$C(a,b) = \frac{b\,(L-a)\,(2La - a^2 - b^2)}{6\,EI\,L}$$

当 $b > a$:

$$C(a,b) = \frac{a\,(L-b)\,(2Lb - b^2 - a^2)}{6\,EI\,L}$$

Maxwell 互等定理:$C(a,b) = C(b,a)$。推广到两组任意荷载即为 Betti 互等定理:

$$\delta_{AB}\,P_A = \delta_{BA}\,P_B$$

其中 $\delta_{BA} = C(a,b)\,P_B$ 为仅系统 B 作用时 a 点的挠度,$\delta_{AB} = C(b,a)\,P_A$ 为仅系统 A 作用时 b 点的挠度。

关于 Maxwell-Betti 互等定理模拟器

🙋
"互等定理"听过这个名字,但到底什么东西在"互等"啊?
🎓
简单说,对线性弹性体施加两组独立荷载时,A 在 B 引起的位移上做的功,与 B 在 A 引起的位移上做的功完全相等,写成式子就是 $\delta_{AB}\,P_A = \delta_{BA}\,P_B$。你在上面调动 P_A 和 a, b 的位置看看,左侧和右侧两个卡片的数字会始终一致。
🙋
真的完全相等!这总不会是巧合吧?
🎓
绝非巧合。影响函数 $C(a,b)$ 本身就满足 $C(a,b) = C(b,a)$,这就是 Maxwell 互等。从公式里能看到 2La·b − a² − b² 关于 a, b 对称。学过有限元的同学应该听过"刚度矩阵必然对称 ($K = K^T$)",那其实就是这一定理在数值层面的化身。
🙋
那是不是任何结构都成立呢?
🎓
前提是"线性弹性"。一旦出现塑性、接触或大变形,响应就依赖荷载历史,互等性就崩了。比如钢梁先屈服后再加另一组荷载,两边就不再相等。所以工程上是在"小变形、线弹性"范围内使用,比如绘制影响线、检验 FEM 结果。
🙋
工程中具体怎么用呢?
🎓
最直观的是节省实验成本。要测涡轮叶片任意点 i 加载、点 j 响应的传递关系,本来要把传感器搬到 100 个 i 上去。但根据互等性,把传感器固定在 j 处,反过来在 100 个 i 上敲击就行——实验装置量级减少。模态分析的检验、伴随法灵敏度计算都是同一思路的延伸。

常见问题

δ_AB 表示仅系统 A 作用时在系统 B 的作用点 b 产生的位移;δ_BA 表示仅系统 B 作用时在系统 A 的作用点 a 产生的位移。可以这样记忆:下标的第一个字母代表"观察点",第二个字母代表"作用点"。本工具利用简支梁影响函数解析地求出两者,并验证 δ_AB·P_A = δ_BA·P_B 始终成立。
线性 FEM 中节点位移 u 与外力 F 满足 F = K·u,Betti 互等等价于 K_{ij} = K_{ji},即 K 是对称阵。一旦组装出的刚度矩阵出现非对称项,往往说明有非保守力(跟随荷载)或非对称阻尼被引入,这成为查错的重要线索。商用 FEM 求解器普遍使用对称 Cholesky 类高效算法,也依赖于这种对称性。
点 b 位移的影响线(单位移动荷载在不同位置时点 b 位移的函数)依据 Maxwell 互等 C(a,b) = C(b,a),正好等于在 b 点施加单位荷载时整梁的挠曲线——这就是 Mueller-Breslau 原理的本质。本工具中移动 a 并读取 δ_BA,等同于在读取"点 a 处的 b 点位移影响线"的纵距。
理论上 δ_AB·P_A 与 δ_BA·P_B 严格相等。本工具采用 IEEE 754 双精度浮点运算,两侧相对误差通常在 1e-10 以下,按显示精度(kN·m 取三位小数)两个卡片的数字完全一致。这说明互等定理不仅是理论结论,在任何一致的数值实现中也精确到舍入误差成立——这是检验 FEM 对称性的有力依据。

实际工程应用

影响线与影响面绘制:桥梁、起重机轨道、楼板等承受移动荷载的结构,设计时需要知道某关键截面的应力或位移随荷载位置变化的规律。利用 Maxwell 互等,只需在响应点施加单位荷载,得到的挠曲线即为该响应的影响线——这就是被桥梁工程师沿用一百多年的 Mueller-Breslau 原理。

节省实验测点:大型结构的振动与位移试验中,激励点与测点互换后传递函数不变。涡轮叶片、飞机结构的模态测试既可固定激振器、移动加速度传感器,也可固定传感器、移动激振点——选择由可达性与成本决定,物理上完全等效。

FEM 结果检验:线性分析的最简检查就是验证 K_{ij} = K_{ji}。若组装后的刚度矩阵出现明显的非对称分量,应立即排查跟随荷载、非对称耦合或边界条件 BUG。

反问题与伴随灵敏度:传感器布置优化、荷载识别、拓扑优化都借助互等性。大规模梯度优化中常用的伴随法(adjoint method)本质上就是 Betti 定理的伪装:通过把荷载点与响应点对调来得到梯度,只需多解一次方程。

常见误解与注意事项

最常见的误解是把互等定理当作"任何力学问题都适用的普适规律"。它实际上要求"线性、弹性、小变形、保守力"四项前提同时满足。出现塑性、损伤、接触、摩擦、屈曲、大变形中的任何一项时一般都不成立。本工具固定 EI、保持挠度在小变形范围内,正是为了保证线性性。当工程中出现"互等性破坏"时,首要怀疑就是潜藏的非线性效应。

另一个常见困惑是"两侧完全相等是不是恰好这样的巧合"。在本工具中无论你怎么调动 L、a、b、P_A,δ_BA·P_A 与 δ_AB·P_B 始终完全相等,这是因为影响函数式子本身关于 a 与 b 对称构造——是数学恒等式而非偶然。代入 a = 1.5、b = 3.5、L = 5、EI = 1e4 手算,可得 C(a,b) = 1.5375×10⁻⁴ m/kN,δ_BA·P_A = δ_AB·P_B = 0.2306 kN·m,数值上完全吻合。

最后是把 Maxwell 与 Betti 混为一谈或彻底分开看待。Maxwell 定理(1864)描述了单位荷载下影响函数的对称性 C(a,b) = C(b,a);Betti 定理(1872)将其推广到两组任意荷载系统。两者本质相同,现代工程中常合称为 Maxwell-Betti 互等定理。核心要点是"线性弹性条件下,对调观测点与作用点不改变响应的对应关系",影响函数公式与 FEM 中的 K = K^T 不过是同一事实的两种表达。