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。

暂停时,拖动滑块即可即时更新结果。

计算结果(相反性实时验证)
δ_BA:A 点加荷载 P → B 点挠度 (mm)
δ_AB:B 点加荷载 P → A 点挠度 (mm)
相反性差值 |δ_AB − δ_BA| (mm)
δ_BA·P_A(左边)(kN·m)
δ_AB·P_B(右边)(kN·m)
当前荷载 P (kN)
相反定理动画 δ_AB = δ_BA

CASE 1(蓝)=在 A 点施加荷载 P,测量 B 点的挠度 δ_BA;CASE 2(红)=把相同荷载 P 加到 B 点,测量 A 点的挠度 δ_AB。两个黄色标记(测得的交叉挠度)始终相等——这就是 Maxwell-Betti 相反定理 δ_AB = δ_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 \gt 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 相反定理模拟器简介

🙋
我听过"相反定理"这个名字,但实际上什么在"相反"呢?
🎓
简单来说,Betti 相反定理是指:在线性弹性体上分别作用两个荷载系时,它们通过对方作用点处的变位产生的"功"相等。用公式表示为 $\delta_{AB}\,P_A = \delta_{BA}\,P_B$。你可以在上面的模拟器中改变 P_A 或 a、b 的位置。左下和右下的数值卡应该始终相等。
🙋
真的!它们完全一致!这不是巧合吧?
🎓
不是巧合。实际上影响函数 $C(a,b)$ 本身就具有 $C(a,b) = C(b,a)$ 的对称性。这就是 Maxwell 相反定理。看公式,当 a 和 b 交换时,2La·b − a² − b² 部分保持不变。学过 FEM 的人会听说"刚度矩阵总是对称的($K = K^T$)",这本质上就是相反定理在数值上的体现。
🙋
那这个定理对所有结构都适用吗?
🎓
不是。"线性弹性"是最基本的前提条件。一旦涉及塑性变形、接触、大变形等,响应就会依赖于荷载历史,相反性就会破坏。比如,把钢材压到屈服后再加不同荷载,就不会满足这个对称关系了。所以在实务中,这个定理主要用在"微小变形、线性材料"的范围内,比如影响线的制作和 FEM 结果验算。
🙋
实务中具体怎么用呢?
🎓
一个很好的例子是实验中的测量节约。比如你要测量汽轮机叶片的 100 个不同点对某个荷载点的响应。按常规要换 100 次传感器位置。但用相反定理,你可以把传感器固定在一个地方,去叶片的 100 个位置依次施加荷载——这样测量装置的成本大大降低。Mohr 挠度计算法和振型分析的验算中也经常用到这个原理。

常见问题

δ_AB 表示"仅作用系 A 时,在系 B 的作用点 b 处产生的变位",δ_BA 表示"仅作用系 B 时,在系 A 的作用点 a 处产生的变位"。记住:添字的第一个字母是"观测点一侧",第二个字母是"作用点一侧"就容易理解了。本工具从简支梁的影响函数解析求得两者,并验证 Betti 乘积 δ_AB·P_A 和 δ_BA·P_B 相等。
在线性 FEM 中,节点位移 u 和外力 F 满足 F = K·u,Betti 相反定理要求 K 是对称矩阵(K_{ij} = K_{ji})。反过来说,如果 K 中出现非对称项,说明单元刚度或边界条件的组装中混入了非保守力(追踪荷载或阻尼),可以作为检查的指标。商用 FEM 求解器能够优化使用对称求解器(Cholesky 等),也正是由于这种对称性。
点 b 的变位影响线(单位移动荷载对 b 点变位的函数)通过 Maxwell 相反定理 C(a,b) = C(b,a),等于在 b 点作用单位荷载时整条梁的挠度形状。这就是穆勒·布雷斯劳原理的出处。在本工具中改变 a 的位置并观察 δ_BA 时,你实际上在读取"点 a 移动时 b 点变位的影响线纵坐标"。
解析上 δ_AB·P_A 和 δ_BA·P_B 严格相等。本工具采用 IEEE 754 双精度浮点计算,两者的相对误差通常在 10⁻¹⁰ 以下,在显示精度(kN·m 小数后 3 位)上完全相同。这说明"相反定理不仅是理论,在数值上也严格成立",也是 FEM 结果对称性检查的基准。

实际应用

影响线、影响面的制作:桥梁、起重机主梁、楼板等受活荷载结构的设计需要了解某个截面的应力或挠度随荷载位置如何变化。用 Maxwell 相反定理,只需在关注点施加单位荷载得到的挠度形状,就直接是影响线,避免了在所有荷载位置计算。这就是穆勒·布雷斯劳原理,长期以来被桥梁设计师活用。

实验测量装置的节约:大型结构的振动、变位试验中,加振点和测量点可以互换(线性系统的相反性)。汽轮机叶片或航空器结构的振型分析中,可以选择固定一个加速度计,依次叩击多个位置,或固定加振器,多点测量——这种灵活性来自这个定理。

FEM 解的验算:从线性分析输出中检查对角单元 K_{ii} 和其转置 K_{ji} 的对称性,可以简便验证单元刚度组装、边界条件、约束方程的一致性。出现非对称项时,怀疑混入了追踪荷载或非保守阻尼。

反演问题与优化:传感器布局最优化、荷载反演(反问题)中,相反性也很强大。通过交换观测点和"虚拟单位荷载位置"计算灵敏度的伴随法(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 = 10⁴ 手工代入计算,会得到 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 只是不同的表述方式。

使用指南

  1. 设定梁的全长 L(1~5 m 范围),输入 A 点和 B 点位置(A < B)
  2. 荷载系 1 在 A 点施加荷载 P_A,荷载系 2 在 B 点施加荷载 P_B(1~50 kN)
  3. 模拟器自动计算 δ_BA·P_A 和 δ_AB·P_B,用 kN·mm 单位验证相反定理的成立

具体计算示例

长度 L=4 m 的简支梁(E=200 GPa、I=500×10^6 mm^4)在 A=1 m、B=3 m、P_A=20 kN、P_B=15 kN 时,FEM 分析得 δ_BA≈1.28 mm、δ_AB≈1.71 mm。左边:1.28×20=25.6 kN·mm,右边:1.71×15=25.65 kN·mm,在数值误差范围内验证相反定理成立。

实务注意事项