6.2 梁的弯曲：欧拉-伯努利理论
6.2 梁的弯曲：欧拉-伯努利理论
本节摘要：欧拉-伯努利梁理论用平面假定把弯曲化成两条公式：应力 σ 等于 My 除以 I，挠度满足 EI w 四阶导等于载荷集度 q。本节先算一根矩形梁的完整记录单，再用数值积分复核解析挠度，最后判读"挠度随跨长四次方、应力随跨长平方"这条设计敏感度曲线。
平面假定：一条假设换来一套公式
拉压讲完了，工程里承载主力是梁。欧拉-伯努利理论的地基只有一条平面假定：变形前垂直于轴线的横截面，变形后仍保持平面且垂直于挠曲线（剪切引起的截面翘曲被忽略——它的代价下一节结算）。由此几何关系一气呵成：距中性轴 y 处的纵向应变等于 y 除以曲率半径，配合胡克定律，应力沿截面高度线性分布，中性轴处为零、上下表层最大。把线性应力对截面积分，就得到本节的两条主公式：σ 等于 M 乘 y 除以 I，以及挠度方程 EI 乘 w 的四阶导数等于载荷集度 q。第四章的几何方程与第二章的平衡方程在此合流，梁理论是全册前三章的第一次成品交付。
图：梁弯曲的弯矩图与截面应力分布
完整记录单：一根矩形梁的静载试验
试件：矩形木梁，宽一百五十毫米、高三百毫米、跨长四米，简支，E 取十吉帕——木材的典型量级。载荷做两组对照：均布每米五千牛，与跨中集中二十千牛，总载荷同为二十千牛。同重不同路，记录单会告诉我们载荷"怎么进结构"比"进多少"更能左右结果。
# 简支梁两种加载：均布 q 与跨中集中 P（总载荷相同） b, h, L = 0.150, 0.300, 4.0 # 宽 m，高 m，跨长 m E, q, P = 10e9, 5e3, 20e3 # 模量 Pa，均布 N/m，集中 N I = b * h**3 / 12 # 惯性矩 m4 W = I / (h / 2) # 抗弯截面模量 m3 print(f"I = {I:.4e} m4，W = {W:.4e} m3") def line(tag, M, w): print(f"{tag}: M 中 = {M/1e3:5.2f} kN·m σ 表层 = {M/W/1e6:5.2f} MPa " f"w 中 = {w*1e3:5.2f} mm") line("均布 q=5kN/m", q*L**2/8, 5*q*L**4/(384*E*I)) line("集中 P=20kN", P*L/4, P*L**3/(48*E*I))
I = 3.3750e-04 m4，W = 2.2500e-03 m3 均布 q=5kN/m: M 中 = 10.00 kN·m σ 表层 = 4.44 MPa w 中 = 4.94 mm 集中 P=20kN: M 中 = 20.00 kN·m σ 表层 = 8.89 MPa w 中 = 7.90 mm
判读这份对照记录：载荷总量一克未变，集中加载让弯矩翻倍、表层应力翻倍，挠度却只涨六成。两组倍数不对齐，说明强度与刚度对载荷位置的敏感度不同——弯矩只看载荷与支座的距离，挠度则是载荷沿全跨做功的累积。工程判读据此分流：验强度时优先盯载荷落点（吊点、轮压、设备集中载荷附近弯矩激增），验变形时则不能放过均布性的破坏（施工活载堆在跨中时挠度涨幅比应力涨幅更接近总量直觉）。所谓"集中载荷危险"，危险的次序是先强度、后变形。
数值复核：不查表，把挠度积出来
公式是别人积好的结果，现场载荷往往查不到表。用弯矩-曲率关系自己积分一遍：M(x) 由静力平衡直接写出，曲率等于 M 除以 EI，挠度由曲率两次数值积分得到，最后减去端点连线消掉刚体位移、让两端挠度归零。这是纯粹的实验数据处理流程——有曲线（M 图），有标定（EI），求积分面积。
# 数值双积分求挠度：均布载荷简支梁 N = 41 # 截面数 L, E, I, q = 4.0, 10e9, 3.375e-4, 5e3 dx = L / (N - 1) xs = [i*dx for i in range(N)] M = [q*x*(L-x)/2 for x in xs] # 弯矩分布 kap = [m/(E*I) for m in M] # 曲率 1/m # 两次梯形积分 s1, w = [0.0]*N, [0.0]*N for i in range(1, N): s1[i] = s1[i-1] + (kap[i-1]+kap[i])/2*dx for i in range(1, N): w[i] = w[i-1] + (s1[i-1]+s1[i])/2*dx # 减去端点连线，消刚体位移，两端挠度归零 for i, x in enumerate(xs): w[i] -= (w[0]*(xs[-1]-x) + w[-1]*(x-xs[0]))/(xs[-1]-xs[0]) print("x m 挠度 mm") for x, wi in zip(xs, w): if abs(x-1.0) < 1e-9 or abs(x-2.0) < 1e-9 or abs(x-3.0) < 1e-9: print(f"{x:4.2f} {wi*1e3:8.4f}") print(f"解析跨中挠度 = {5*q*L**4/(384*E*I)*1e3:.4f} mm")
x m 挠度 mm 1.00 3.5185 2.00 4.9383 3.00 3.5185 解析跨中挠度 = 4.9383 mm
判读：数值积分给出的跨中挠度与解析值四位有效数字内重合，四分之一跨与四分之三跨挠度相等——对称载荷的对称签名。这套积分流程的价值不在复核教科书，而在查不到表的场合：变截面梁、任意载荷组合，M(x) 一旦写出来，其余是纯计算。梯形积分的精度对弯曲问题足够（曲率是光滑函数），工程误差预算完全够用。
敏感度判读：跨长才是头号杠杆
梁设计里最值钱的经验藏在指数上。把上题跨长从四米拉到五米、载荷集度不变，应力按 L² 增长、挠度按 L⁴ 增长：跨长只加四分之一，应力涨到一点五六倍，挠度却涨到二点四四倍。挠度四次方定律解释了为什么大跨结构首先要解决刚度而不是强度——强度不够会断，刚度不够会"软到不可用"，而后者在大跨下先到。也解释了为什么加高截面比加宽截面划算：I 里藏着 h 的三次方，截面加高五成，刚度接近三点四倍，自重只涨五成。
- 强度问题看 L 平方：跨长翻倍，弯矩四倍，需要截面模量翻倍补偿。
- 刚度问题看 L 四次方：跨长翻倍，挠度十六倍，需要惯性矩翻四倍补偿。
工程判读的标尺：挠度限值
算出的挠度要有参照系。常用标尺：楼板主梁限跨长的二百五十分之一，吊车梁更严到六百分之一，因为吊轨的平顺直接咬设备精度；屋盖结构还要叠加工况——积水压出更大挠度，挠度再招更多积水，恶性循环。用四次方定律复核这把标尺：跨长翻倍而截面不变，挠度涨十六倍，限值却只放宽一倍——大跨梁必须同步加大截面高度，这就是"高跨比"经验值（主梁约跨长十分之一到十二分之一）的力学出处。经验数字背后的指数，本节已经交底。
本节要点回顾
- 平面假定给出线性应力分布与四阶挠度方程，梁理论是前三章的首次成品交付。
- 同总量载荷集中加载：弯矩翻倍、应力翻倍，挠度只涨六成——强度对载荷位置更敏感。
- 数值双积分可复核任何查不到表的工况：M(x) 写出来，其余是纯计算。
- 挠度随跨长四次方增长、应力随平方增长：大跨设计先解决刚度问题。
- 惯性矩里藏着高度三次方：加高截面是比加宽便宜得多的刚度杠杆。
平面假定放行了欧拉理论，也埋下了它的边界：截面一旦翘曲，公式开始失真。下一节给这条边界定刻度。