非線形後座屈解析
理论与物理
概要 — 为何要追踪“后屈曲”
线性屈曲分析已经能知道屈曲载荷了，为什么还需要追踪屈曲的“之后”呢？
有三个理由。
1. 实际结构屈曲后未必立即崩溃。 例如飞机机翼的外蒙皮在运行载荷下会发生局部屈曲，但整体结构仍具有足够的强度。这被称为“后屈曲强度的利用”。
2. 特征值屈曲是上界值，会高估实际崩溃载荷。 包含初始缺陷和材料非线性的实际崩溃载荷，只有通过非线性分析才能准确评估。
3. 确定崩溃模式。 了解结构“如何破坏”是安全设计的基础。通过追踪载荷-位移的全路径，还可以评估能量吸收能力和延性。
飞机的外蒙皮屈曲了还能继续使用！？
是的。薄铝蒙皮在加强筋之间会在设计载荷以下发生屈曲，但加强筋会重新分配载荷，使整体结构仍能承载。如果能正确评估这种后屈曲承载能力，就能进一步实现结构轻量化。相反，如果忽略后屈曲而将所有部分都设计成不屈曲，就会导致结构过重。
载荷-位移路径与奇异点
能再稍微从数学角度讲解一下“后屈曲”吗？
在非线性结构力学中，载荷参数 $\lambda$ 与位移向量 $\{u\}$ 的关系称为平衡路径（equilibrium path）。这条路径上存在两种奇异点：
分岔点（bifurcation point） — 平衡路径分叉的点。完美结构中的欧拉屈曲就属于此类。
极限点（limit point） — 载荷达到极大值的点。跳跃屈曲（snap-through buckling）属于此类。
两者都是切线刚度矩阵 $[K_T]$ 变得奇异呢。
是的。但物理意义完全不同。分岔点处，结构“切换到另一种变形模式”。极限点处，结构“无法再承受更多载荷”。后屈曲分析需要正确追踪这两者。
Koiter的后屈曲理论
有人奠定了后屈曲的理论基础吧。
是Warner T. Koiter。他在1945年于代尔夫特理工大学的博士论文（因用荷兰语撰写而长期不为人知）中，系统地对后屈曲路径的稳定性进行了分类。
将分岔点附近的后屈曲路径进行幂级数展开：
其中 $\xi$ 是模态振幅参数。
- $a \neq 0$（非对称分岔） — 后屈曲路径向一侧倾斜。对缺陷极其敏感
- $a = 0, b > 0$（稳定对称分岔） — 后屈曲后载荷上升。平板的屈曲属于此类
- $a = 0, b < 0$（不稳定对称分岔） — 后屈曲后载荷下降。圆柱壳的屈曲属于此类
$b < 0$ 是“不稳定”且对缺陷敏感性高…圆柱壳的实验值只有理论值的两成就是因为这个吧！
正是如此。Koiter证明了可以从“后屈曲路径的形状”预测“对缺陷的敏感性”。这是结构力学中最重要的理论贡献之一。
切线刚度矩阵的构成
非线性分析中的 $[K_T]$ 与线性屈曲的 $[K_0] + \lambda[K_\sigma]$ 有什么不同？
在大变形理论中，切线刚度矩阵由三项构成：
- $[K_0]$ — 小位移的弹性刚度（材料刚度）
- $[K_\sigma]$ — 几何刚度（初始应力的影响）
- $[K_L]$ — 大位移刚度（位移引起的形状变化的影响）
线性屈曲分析忽略了 $[K_L]$ 呢。
没错。线性屈曲基于“屈曲前的变形是微小的”这一假设，因此 $[K_L] \approx 0$。但在后屈曲中，位移变大，$[K_L]$ 就不可忽略了。包含 $[K_L]$ 正是“几何非线性分析”的核心。
如果再加入材料非线性（塑性），$[K_0]$ 本身也会依赖于应力状态…。
是的。在弹塑性后屈曲分析中，所有三项都随位移和应力而变化。因此才需要用增量迭代法（Newton-Raphson法的载荷增量版）一点点地追踪路径。
总结
我来整理一下非线性后屈曲的理论。
要点：
- 后屈曲对于“准确评估崩溃载荷”和“确定崩溃模式”是必要的
- 平衡路径上的分岔点和极限点 — $[K_T]$ 变得奇异的两种点
- Koiter的理论 — 后屈曲路径的形状（$a, b$ 的符号）决定了缺陷敏感性
- $[K_T] = [K_0] + [K_\sigma] + [K_L]$ — 大位移项 $[K_L]$ 是非线性分析的核心
- 有些结构屈曲后仍能承载 — 利用后屈曲强度是轻量化的关键
如果说线性屈曲告诉我们“何时屈曲”，那么非线性后屈曲就是告诉我们“屈曲之后会怎样”呢。
理解得非常完美。下次，我们将详细看看如何实际进行这种路径追踪——Riks法（弧长法）的数值算法。
后屈曲行为的发现：von Kármán的薄板理论
首次将屈曲后仍能承载的“后屈曲（post-buckling）强度”理论化的是Theodore von Kármán（1932年）。薄板在屈曲后会将压缩载荷集中到板边缘，从而能够承受初始屈曲载荷的2〜4倍。这个“有效宽度”概念是现代冷弯型钢设计规范（AISI S100）的基础，是能将薄板结构重量减轻30〜50%的重要见解。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重，越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施加力，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的特性就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形”，强度是“不易破坏”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓紧固力…全都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。弹一下吉他的弦试试。声音会一直响吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦变成了热。汽车的减震器也是同样原理——故意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设定适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项・换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入为mm时，载荷・弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
Riks法（弧长法）的原理
Newton-Raphson法无法通过屈曲点对吧。Riks法有什么不同？
根本区别在于控制量。Newton-Raphson法以载荷为控制量。一点点增加载荷，在每个增量步求解平衡。但在载荷-位移曲线的峰值（极限点）处，同一个载荷对应两个平衡解，导致求解器混乱。
Riks法以弧长为控制量。将载荷 $\lambda$ 和位移 $\{u\}$ 同时作为未知数，沿着平衡路径的“距离”前进：
$\Delta s$ 是弧长（arc length），$\psi$ 是缩放系数…。将载荷和位移作为“一条曲线上的距离”统一处理呢。
Modified Riks法的实现
Abaqus中使用的“modified Riks法”与普通的Riks法有什么不同？
这是对Crisfield（1981）和Riks（1979）方法的改进，主要区别在于弧长约束的应用方式。要点是：
1. 预测步 — 沿前一增量的切线方向进行预测
2. 弧长约束 — 用弧长约束来限制从预测步开始的修正
3. 修正迭代 — 用Newton-Raphson法减小残差
Crisfield（1981）和Riks（1979）方法的改进，主要区别在于弧长约束的应用方式。要点是：
1. 预测步 — 沿前一增量的切线方向进行预测
2. 弧长约束 — 用弧长约束来限制从预测步开始的修正
3. 修正迭代 — 用Newton-Raphson法减小残差
相关主题
なった
詳しく
報告