7.1 失稳的类型与判别

7.1 失稳的类型与判别

本节摘要：失稳是平衡状态失去稳定性，不是应力超限。本节用一根刚性杆加一个扭转弹簧的最小模型，演示分岔点屈曲与初缺陷路径，演算能量判据的判定过程，并给出杆、板、壳三类构件后屈曲行为的对照记录。

强度逻辑之外的另一种破坏

前六章的破坏判据都是"应力说了算"：算出应力，对照强度，完事。失稳破坏绕开这套逻辑——压杆侧弯时截面应力可能只有屈服的一半，破坏却已经发生。判别失稳，看的不是应力大小，而是平衡状态的稳定性：给平衡位置一个小扰动，构件是弹回来（稳定）、停在新位置（随遇）、还是越偏越远（不稳定）。第六章的强度验算回答"材料撑不撑得住"，本章的稳定性验算回答"平衡保不保得住"，两套验算缺一不可。

用能量语言说最干净：总势能取极小值对应稳定平衡，取极大值对应不稳定平衡，判定工具是势能对广义坐标的二阶导数。二阶导变号的位置，就是临界载荷出现的位置。

最小模型：一根杆一个弹簧演分岔

试件模型简到不能再简：长 L 等于一米的刚性竖杆，底端铰支并装扭转刚度 k 等于一百牛米每弧度的弹簧，顶端受竖向载荷 P。弹簧提供恢复力矩，载荷提供倾覆力矩。总势能 Π 等于弹簧应变能 kθ²/2 减去载荷下降的功 PL(1−cosθ)，对转角 θ 求导得平衡方程，再求一次导得稳定性判据。完美系统（杆绝对竖直）与带初缺陷系统（安装时预偏 θ₀ 等于零点零五弧度）各记一份。

# 刚性杆-弹簧模型：分岔与初缺陷路径
import math

L, k = 1.0, 100.0            # 杆长 m，弹簧刚度 N·m/rad
P_cr = k / L                  # 完美系统分岔载荷

print(f"完美系统临界载荷 P_cr = k/L = {P_cr:.0f} N")

# 完美系统：非零转角平衡要求 P = kθ/(L sinθ)，恒大于等于 k/L
# 带初缺陷 θ0：平衡路径 P = k(θ-θ0)/(L sinθ)
theta0 = 0.05
print("\n转角θ 度   完美杆所需P   带缺陷杆所需P")
for deg in (5.7, 11.5, 28.6, 57.3, 86.0):
    t = math.radians(deg)
    P_imperfect = k * (t - theta0) / (L * math.sin(t))
    P_perfect   = k * t / (L * math.sin(t))
    print(f"{deg:6.1f}    {P_perfect:9.1f}     {P_imperfect:9.1f}")

# 能量二阶导判据：完美杆 θ=0 平衡在 P 低于与高于 k/L 时的定性
for P in (80.0, 120.0):
    d2 = k - P * L * 1.0      # d2Π/dθ2 = k - PL cosθ，θ=0 时 cosθ=1
    tag = "稳定" if d2 > 0 else "不稳定"
    print(f"P = {P:.0f} N：二阶导 = {d2:.0f}，直线平衡 {tag}")
​

完美系统临界载荷 P_cr = k/L = 100 N

转角θ 度   完美杆所需P   带缺陷杆所需P
   5.7        100.2          49.8
  11.5        100.7          75.6
  28.6        104.2          93.8
  57.3        118.8         112.9
  86.0        150.5         145.5

P = 80 N：二阶导 = 20，直线平衡 稳定
P = 120 N：二阶导 = -20，直线平衡 不稳定
​

判读三段记录：其一，完美杆在载荷一百牛以下保持笔直，越过就翻——刚性分岔，路径在临界点分裂。其二，带缺陷杆从零载荷起就有侧偏，载荷五十牛时转角已到五点七度：缺陷没有降低"能扛住的最大载荷"，却让构件从加载第一天起就在弯，刚度表现为先软后硬。其三，能量二阶导在 k/L 处变号，定量标出分岔位置——三段记录合起来就是失稳分析的完整语法：先写势能，再求平衡，最后用二阶导定性质。

图：杆板壳三类构件的后屈曲行为对照

分岔点与极值点：两类失稳的户口

上一节的弹簧杆是分岔点失稳的标准样本：平衡路径在临界载荷处分裂成两支，理想与缺陷路径形状相似只差圆头。工程里还有另一类：极值点失稳——载荷位移曲线一路上升到一个峰值然后掉头，载荷扛不住自己先降，如压溃的薄壁筒、跳塌的浅拱。两者的实验签名不同：分岔类失稳前挠度缓慢增长、临界的瞬间突变；极值类失稳前已有明显非线性征兆，曲线掉头即垮。数值分析上的分工也随之而定：分岔类用特征值屈曲分析估上限，再用非线性分析修正；极值类只有非线性增量分析一条路，特征值只会给一个偏乐观的假峰值。

三类构件的后屈曲富余差异，同样要从载荷位移路径上读：薄板屈曲后膜力重新分布，还能继续升载到临界值的两倍上下；受弯薄腹板梁屈曲后靠张力场继续传剪力，是同一种富余的工程化利用。圆柱壳相反，越过临界点路径直接俯冲，载人航天器与储罐设计对此最敏感。判别顺序因此固定下来：先辨类型（分岔还是极值、富余还是俯冲），再选工具，最后才谈数字。

缺陷敏感度的实测办法

判别工具箱里还有一件便宜的实验兵器：索斯韦尔作图。压杆试验的载荷挠度记录按 δ 比 P 作纵轴、δ 作横轴，初缺陷主导的散点会落在一条直线上，斜率的倒数就是临界载荷——即使试件弯得没法直接观测分岔点，直线一拉，P_cr 自己浮出来。老式木桁架、砖柱鉴定里它至今好用，因为它只需要载荷与挠度两个易测量。配合上一节的模型曲线，实验室判别失稳的完整弹药已经齐了：能量判据定性、特征载荷定量、索氏作图补实测，三件套互为证人。

一根桁架压杆的判别全程

把判别流程落到一根真杆上。桁架上弦：两端节点板连接，近似两端铰接，几何长度三米；钢管外径一百一十四毫米壁厚四毫米，回转半径三十八点九毫米，柔度算得七十七——比欧拉门槛七十五点六略高，弹性失稳区，欧拉公式合法。临界载荷估约四百五十千牛量级（细账下一节清）。判别三问走完：类型问——细长压杆是分岔型，缺陷温和；手段问——特征值屈曲分析给上限，几何缺陷的非线性修正交给规范曲线；余量问——设计载荷与临界值之比压在零点六以内，给制造偏差留白。三问的答案全部写在载荷位移路径的形状里，这正是本节模型演算的价值：有了最小模型的曲线形状做底版，任何压杆的试验记录都能对号入座。

本节要点回顾

• 失稳判据是平衡稳定性：势能二阶导变号即临界，与应力大小无关。
• 刚性杆弹簧模型给出分岔点失稳全貌：完美杆刚性分岔，缺陷杆圆头渐变。
• 初缺陷不必然降低极限载荷，但让构件从加载起就带侧偏、刚度先软后硬。
• 分岔点失稳可用特征值估上限，极值点失稳只能靠非线性增量分析。
• 后屈曲富余排序：薄板有倍数级富余、薄腹板可利用张力场、圆柱壳必须重扣。

模型说清了道理。轮到工程压杆——欧拉公式与它的适用刻度，下一节给全。