欧拉屈曲分析
理论与物理
什么是屈曲
屈曲和拉伸破坏不一样，是突然发生的吗？我看到踩铝罐时它一下子就瘪了，很好奇。
那个铝罐正是屈曲的绝佳例子。从正上方慢慢按压罐子，在达到某个力之前它都安然无恙。但一旦超过临界点，侧面就会突然向外凸出，一下子瘪掉。材料并没有破裂，但形状失去了稳定性——这就是屈曲。
材料没事但结构却崩溃了，这不可怕吗？仅靠应力计算无法预测吗？
这正是与拉伸、弯曲决定性的不同点。拉伸可以判断“应力是否超过屈服点？”，但屈曲是几何稳定性问题，所以即使应力低于屈服点的一半也可能发生。1757年，莱昂哈德·欧拉在世界上首次用数学公式对此进行了表述。所以至今仍称之为“欧拉屈曲”。
1757年！？ 270多年前的理论现在还在CAE中使用吗？
欧拉是当时数学的巅峰。他使用变分法解析地求出了“柱子在哪一个载荷下会变得不稳定”。这个公式成为了当今CAE中屈曲特征值分析的起点。我们先来推导一下看看。
控制方程的推导
欧拉公式是怎么推导出来的？
首先，考虑承受轴向力 $P$ 的柱子稍微向横向弯曲的状态。设挠度为 $y(x)$，则由于该挠度会产生附加弯矩 $M = -Py$。将其代入欧拉-伯努利梁的弯曲方程，得到：
$EI$ 是弯曲刚度对吧。式子看起来很简单……
这是一个二阶线性常微分方程，令 $k^2 = P/(EI)$，则 $y'' + k^2 y = 0$。通解为：
正弦和余弦！和简谐振动的公式形式一样呢。
正是如此。现在代入边界条件。如果是两端铰支（$y(0) = 0$, $y(L) = 0$），则 $B = 0$，得到 $A\sin(kL) = 0$。我们想要求有挠度的解（$A \neq 0$），所以：
$n = 1$ 对应最小的载荷对吧？在实际结构中，会在那里发生屈曲吗？
没错。将 $k^2 = P/(EI)$ 代回，得到最小的临界载荷：
这就是著名的欧拉屈曲载荷。$n = 2, 3, \ldots$ 是高阶模态，但在实际结构中，最低阶模态（$n=1$，半波正弦曲线）是主导的。
边界条件与有效屈曲长度
两端铰支的情况明白了。但实际的柱子，有的是用螺栓固定在基础上，有的是上端连接到楼板，对吧？
这时就出现了有效屈曲长度系数 $K$。它将实际边界条件下的屈曲长度换算为“等效的两端铰支柱的长度”。一般公式如下：
| 边界条件 | K值 | $P_{cr}$ 比率 | 实际工程中的例子 |
|---|---|---|---|
| 两端铰支 | 1.0 | 1.0（基准） | 桁架构件的理想化 |
| 一端固定·一端自由（悬臂） | 2.0 | 0.25 | 旗杆、独立柱 |
| 两端固定 | 0.5 | 4.0 | 混凝土墙之间的柱子 |
| 一端固定·一端铰支 | 0.7 | 2.04 | 钢框架的柱脚固定·梁连接 |
$K = 2.0$ 的悬臂，和相同长度的两端铰支柱相比，$P_{cr}$ 只有四分之一！？差别这么大吗？
因为 $P_{cr}$ 与 $(KL)^2$ 成反比。悬臂柱的设计需要特别注意。相反，如果是两端固定，$K = 0.5$，所以同样柱子的 $P_{cr}$ 是4倍。例如，建筑钢结构的柱子，通过横向连接梁或支撑约束端部，可以降低 $K$ 值——也就是增强抗屈曲能力。
但现实中的螺栓连接，既不是“完全固定”也不是“完全铰支”吧？
问得好。实际的连接节点是“半刚性连接”，$K$ 值与理论值不同。所以设计规范（如AIJ或Eurocode 3）指导要从安全角度出发，取较大的 $K$ 值。在CAE中，也有使用旋转弹簧单元来表现半刚性连接的技术。这里如果考虑不足会很危险。
长细比与适用范围
又粗又短的柱子不会屈曲吗？我很难想象压缩一个粗短的柱子会向侧面弯曲。
这是个核心问题。用临界应力重写一下就容易理解了：
这里 $r = \sqrt{I/A}$ 是截面的回转半径，$KL/r$ 是长细比（slenderness ratio）。长细比越大，$\sigma_{cr}$ 就越低——也就是说，柱子越细长越容易屈曲。
那反过来，非常粗的柱子（长细比小的柱子），$\sigma_{cr}$ 会超过屈服应力吧？这在物理上可能吗？
这正是欧拉公式的适用范围限制。对于 $\sigma_{cr} > \sigma_Y$（屈服应力）这样的短柱，材料会先于屈曲发生屈服。欧拉屈曲有效的条件是：
对于钢材（$E = 200$ GPa, $\sigma_Y = 250$ MPa），大约是 $KL/r > 89$。对于长细比小于此值的“中长柱”，需要使用Johnson抛物线公式或切线模量理论。
在CAE中进行屈曲分析时，不注意这个区分会很糟糕吗？
非常糟糕。线性屈曲分析（特征值分析）的前提是欧拉型弹性屈曲。如果只对长细比小的构件进行线性屈曲分析，会得到比实际更高的屈曲载荷，导致偏危险的评估。这种情况下需要进行非线性屈曲分析（几何非线性+材料非线性）。
初始缺陷的影响
教科书上的柱子是完美笔直的，但实际工厂制造的柱子都有点弯曲吧？这会影响屈曲吗？
影响非常大。实际的柱子必然存在制造误差、焊接残余应力、运输变形导致的初始缺陷（initial imperfection）。完美的柱子在理论上直到 $P_{cr}$ 都能保持笔直，但有初始缺陷时，挠度会随着载荷增加而逐渐增大。
假设初始挠度 $y_0 = a_0 \sin(\pi x / L)$，则在载荷 $P$ 作用下的挠度为：
当 $P$ 接近 $P_{cr}$ 时分母趋近于零……挠度发散！
是的。这就是缺陷敏感性。在实际结构中，在达到 $P_{cr}$ 之前就会因过大的挠度或塑性而破坏。所以设计中会乘以安全系数，在CAE中，引入初始缺陷的非线性分析也很重要。
具体是怎么引入初始缺陷的呢？
典型的方法有两种：
- 屈曲模态法 —— 将线性屈曲分析的一阶模态形状按微小振幅（构件长度的 $L/1000$ 左右）缩放后赋予初始形状
- 设计规范法 —— 根据制造公差直接反映缺陷模式（例如：AIJ的 $L/1000$、Eurocode 的 $L/500$）到形状上
在Abaqus中，使用 *IMPERFECTION 关键字就可以直接将一阶模态作为初始缺陷导入，非常简单。
理论总结
我想整理一下到目前为止理论中需要掌握的重点。
有四点：
- 屈曲是稳定性问题而非材料破坏 —— 发生在弹性范围内的几何失稳现象
- $P_{cr} = \pi^2 EI / (KL)^2$ —— 边界条件（$K$值）支配屈曲载荷
- 适用范围由长细比决定 —— 对短柱使用欧拉公式会导致偏危险的过高估计
- 初始缺陷对实际行为影响很大 —— 不能盲目相信完美柱子的理论值 $P_{cr}$
屈曲原来是“形状是否稳定”的问题，而不是“材料能否承受”的问题啊。如果只盯着应力值，似乎很容易忽略这一点。
这种直觉很重要。屈曲分析不仅仅是输出数值，其本质是物理上理解为什么该结构会变得不稳定。
欧拉的天才与屈曲理论的诞生
1744年，莱昂哈德·欧拉发表的屈曲载荷公式 $P_{cr}=\pi^2 EI/(KL)^2$ 震惊了当时的学者。因为“临界载荷完全不依赖于材料强度，仅由截面形状（$I$）和长度（$L$）决定”这一点，无人能够预料。据说欧拉本人对这个公式的实用性也缺乏信心，但它作为细长柱设计的基础，至今已被使用了270多年。
数值解法与实现
基于FEM的屈曲分析公式化
明白了欧拉的理论公式，那么在FEM中是如何求解的呢？和普通的 $[K]\{u\} = \{F\}$ 不一样吗？
本质上不同。普通的静力分析只是求解联立方程组，但屈曲分析是特征值问题：
这里 $[K]$ 是通常的弹性刚度矩阵，$[K_\sigma]$ 是几何刚度矩阵（应力刚度矩阵），$\lambda$ 是屈曲载荷系数，$\{\phi\}$ 是屈曲模态形状。
特征值问题，在振动分析中也出现过 $([K] - \omega^2 [M])\{\phi\} = \{0\}$ 呢。形式很像……
在数学上是完全相同的结构。振动的质量矩阵 $[M]$ 在屈曲中被几何刚度矩阵 $[K_\sigma]$ 替换了。掌握这个对应关系，在阅读求解器手册时理解速度也会快很多。
$[K_\sigma]$ 到底是什么？从名字看就和普通的刚度矩阵不一样。
通常的 $[K]$ 由材料的弹性模量和形状决定——与载荷无关。而 $[K_\sigma]$ 依赖于当前的应力状态。存在压应力时，$[K_\sigma]$ 的贡献为负，会降低整体刚度。
啊，明白了！随着压缩载荷增加，结构变得越来越“柔软”，在某个点整体刚度变为零——那就是屈曲点对吧！
理解得非常完美。当 $[K] + \lambda [K_\sigma]$ 变为奇异（行列式为零）时的 $\lambda$ 就是屈曲载荷系数。参考载荷 $\{F_{ref}\}$ 乘以 $\lambda$ 就得到临界载荷。
几何刚度矩阵的构建
なった
詳しく
報告