熱座屈解析
理论与物理
热屈曲是什么
热屈曲是温度升高就会屈曲吗？即使没有载荷也会？
是的。简单来说，受约束的结构体被加热时，想要膨胀却无法膨胀，因此会积累压缩应力。当该压缩应力超过临界值的瞬间，即使外部载荷为零，也会发生屈曲。
诶，载荷为零时屈曲不是有违直觉吗？实际在什么情况下会发生呢？
举三个身边的例子：
- 铁路轨道的“太阳弯折”：夏季轨道受热，但由于被螺栓固定在枕木上，无法沿纵向伸长。压缩力不断积累，超过某个温度时就会横向蜿蜒屈曲。这是导致列车停运的原因之一。
- 飞机的外蒙皮面板：超音速飞行时的气动加热使外蒙皮温度升高，但由于受到框架（骨架）的约束而产生压缩应力，导致面板屈曲。
- 管道的热膨胀：化工厂的高温管道在支撑点之间受到约束而被压缩，有时会发生管道因屈曲而弯曲的事故。
这些都是“约束 + 温度升高 → 压缩应力 → 屈曲”的相同机制。
原来如此！“约束”是关键啊。如果能自由伸长就不会屈曲了。
没错。如果是自由端，只会伸长 $\Delta L = \alpha L \Delta T$ 就结束了。但如果两端固定，伸长量被约束为零，那么这部分就会完全转化为压缩应力 $\sigma = E \alpha \Delta T$。这就是热屈曲的起点。
热应力的产生机制
能再详细讲解一下热应力的公式吗？
考虑一个单轴约束的杆。温度升高 $\Delta T$ 时，自由热应变为：
但如果两端固定，总应变就为零（$\varepsilon_{\text{total}} = 0$），因此会产生抵消热应变的弹性应变：
根据胡克定律，约束结构中产生的压缩应力为：
负号表示压缩。例如钢材（$E = 200\,\text{GPa}$，$\alpha = 12 \times 10^{-6}\,/\text{K}$）在 $\Delta T = 50\,\text{K}$ 时...
$\sigma = 200 \times 10^3 \times 12 \times 10^{-6} \times 50 = 120\,\text{MPa}$。仅仅50度的温升就会产生120 MPa的压缩应力。对于薄板来说，这个应力足以引起屈曲。
临界温差的理論公式
那么，具体“升高多少度会屈曲”是怎么求的呢？
考虑最基本的情况，即四边简支的矩形板均匀受热。这是欧拉屈曲理论的热版本。
临界屈曲应力由板的屈曲公式给出：
其中 $k$ 是屈曲系数（四边简支、单轴压缩时 $k = 4$），$t$ 是板厚，$b$ 是短边宽度。令其与热应力 $\sigma = E \alpha \Delta T$ 相等，得到：
哦，杨氏模量 $E$ 消掉了！
很好的观察点。因为热应力是 $E \alpha \Delta T$，而临界应力也与 $E$ 成正比，所以 $E$ 被约分消掉了。也就是说，临界温差仅由板厚、宽度、热膨胀系数和泊松比决定，与杨氏模量无关。这是热屈曲特有的性质。
我想通过具体数值来感受一下。例如铝板（$\alpha = 23 \times 10^{-6}$，$\nu = 0.33$），$t = 2\,\text{mm}$，$b = 200\,\text{mm}$ 会怎样？
取 $k = 4$（四边简支）计算：
$\Delta T_{cr} = \dfrac{4 \times \pi^2 \times (0.002)^2}{12 \times 23 \times 10^{-6} \times (1 - 0.33^2) \times (0.2)^2} \approx 1.6\,\text{K}$
仅仅1.6度就会屈曲。可以看出薄铝板对热屈曲有多么敏感。这就是飞机外蒙皮设计中重视热屈曲的原因。
各项的物理意义
- $\alpha$（线膨胀系数）：温度每升高1K的伸长率。铝的值约为钢的2倍，因此铝结构对热屈曲更敏感。
- $t/b$（板厚/短边宽度比）：相当于板的细长比。这个比值越小（板越薄越宽），$\Delta T_{cr}$ 越低，越容易屈曲。与 $t^2$ 成正比，所以板厚加倍，$\Delta T_{cr}$ 变为4倍。
- $\nu$（泊松比）：$(1-\nu^2)$ 项是双轴应力状态的影响。泊松比越大，$\Delta T_{cr}$ 会略微减小。
- $k$（屈曲系数）：随边界条件和长宽比变化。四边固定时 $k \approx 10.67$，是简支的约2.7倍。
假设条件与适用范围
- 线弹性材料（在屈服前范围内有效）
- 小挠度理论（挠度相对于板厚足够小）
- 均匀的温度分布（板面内和板厚方向均均匀）
- 无初始缺陷（假设为理想平板）
- 忽略温度相关的材料属性
- 实际结构需要应用折减系数（约0.5〜0.8）
量纲分析与单位制
| 变量 | SI单位 | 代表值・备注 |
|---|---|---|
| $\alpha$（线膨胀系数） | 1/K | 钢: 12×10⁻⁶、铝: 23×10⁻⁶、钛: 8.6×10⁻⁶ |
| $E$（杨氏模量） | Pa | 钢: 200 GPa、铝: 70 GPa（对 $\Delta T_{cr}$ 无直接影响） |
| $t$（板厚） | m | 飞机外蒙皮: 1〜3 mm、压力容器: 10〜50 mm |
| $\sigma_{cr}$（临界屈曲应力） | Pa | 可通过 $E \alpha \Delta T_{cr}$ 反算 |
作为特征值问题的公式化
在FEM中是如何公式化的？
线性屈曲分析被公式化为特征值问题。使用热应力产生的应力刚度矩阵 $[K_\sigma]$：
各符号的含义如下：
- $[K_0]$：通常的弹性刚度矩阵
- $[K_\sigma(T)]$：由温度 $T$ 引起的热应力计算得到的应力刚度矩阵（也称为几何刚度矩阵）
- $\lambda$：特征值 = 临界温度载荷倍数
- $\{\phi\}$：屈曲模态形状
$\lambda = 1$ 表示在给定的温度载荷下刚好屈曲。$\lambda = 0.5$ 表示在温度的一半时屈曲，所以危险。$\lambda > 1$ 则表示安全。
原来如此，$\lambda$ 起到了类似安全系数的角色。
边界条件与临界温度的关系
改变边界条件，临界温度会变化多少？
屈曲系数 $k$ 随边界条件变化很大，因此临界温度也成比例变化。以正方形板（$a/b = 1$）为例进行比较：
| 边界条件 | 屈曲系数 $k$ | $\Delta T_{cr}$ 的倍数（以SS为基准） |
|---|---|---|
| 四边简支 (SSSS) | 4.0 | 1.0 |
| 四边固定 (CCCC) | 10.67 | 2.67 |
| 2边固定+2边自由 (CCFF) | ≈2.2 | 0.55 |
| 全边自由端（有面内约束） | 1.0 | 0.25 |
仅仅变成四边固定，临界温度就变为2.7倍。设计中增加固定端是热屈曲对策的基本。反之，自由端多的结构则较弱。
仅改变固定条件，临界温度就有近4倍的差异，这在设计上影响很大啊。
数值解法与实现
FEM 两步分析流程
用FEM分析热屈曲时，具体是什么步骤？
基本按两步进行：
Step 1：热应力分析（带热载荷的静力分析）
- 输入温度分布 $\Delta T(\mathbf{x})$ 作为载荷
- 将热应变 $\varepsilon_{\text{th}} = \alpha \Delta T$ 转换为等效外力载荷
- 求解 $[K]\{u\} = \{F_{\text{th}}\}$ 得到应力场 $\sigma_{ij}$
Step 2：线性屈曲特征值分析
- 根据 Step 1 的应力场构建应力刚度矩阵 $[K_\sigma]$
- 求解特征值问题 $([K_0] + \lambda [K_\sigma]) \{\phi\} = \{0\}$
- 最小的正特征值 $\lambda_1$ 即为临界温度倍数
Step 1 中输入温度100度，如果得到 $\lambda_1 = 0.6$，那就是在60度时屈曲吗？
没错。因为 $\Delta T_{cr} = \lambda_1 \times \Delta T_{\text{ref}}$，所以 $0.6 \times 100 = 60$ 度就是临界温差。实际工作中会乘以安全系数1.5〜2.0，所以设计允许温差 $\Delta T_{\text{allow}} = 30 \sim 40$ 度。
温度分布不均匀的情况
如果不是均匀加热，而是单面加热或局部加热，则需要先在 Step 1 之前执行稳态/非稳态的热传导分析以求得温度分布，然后将该结果作为温度载荷映射到 Step 1 的结构分析中。也就是说实际上是三步（热传导 → 热应力 → 屈曲特征值）。
单元选择与公式化
热屈曲分析应该用什么单元？
根据结构形状区分使用：
| 结构类型 | 推荐单元 | 注意事项 |
|---|---|---|
| 薄板面板 | 壳单元（S4R, SHELL181等） | 板厚方向5〜7个积分点以捕捉弯曲梯度 |
| 厚板・块体 | 实体单元（C3D20R, SOLID186等） | 板厚方向至少3〜4层网格 |
| 轨道・梁 | 梁单元（B31, BEAM188等） | 无法捕捉薄壁部分的局部屈曲，需与壳单元并用 |
| 管道 | 管/壳单元 | 要捕捉周向屈曲模态，需要足够的周向分割 |
壳单元用得很多啊。板厚方向的积分点数有什么影响？
当板厚方向存在温度梯度时，积分点太少就无法准确捕捉弯曲应力的分布。例如单面加热时，板厚方向会产生线性温度分布，但如果只有3个积分点，可能会低估弯曲热应力。实际工作中推荐至少5点，最好7点。
非线性后屈曲分析
仅靠线性屈曲分析就足够了吗？如果想了解屈曲后的行为该怎么办？
问得好。线性屈曲能告诉我们“何时屈曲”，但“屈曲后如何”则不清楚。要追踪后屈曲行为，需要进行非线性静力分析（几何非线性）。
步骤如下：
- 通过线性屈曲分析获取屈曲模态形状 $\{\phi\}$
- 将该模态形状按比例缩放后作为初始缺陷添加到几何形状中（典型值为板厚的10〜50%）
- 执行温度载荷逐步增加的非线性静力分析（NLGEOM=ON）
- 追踪载荷-位移曲线（Riks法或Arc-length法）
如果不加初始缺陷会怎样？
对于完美的几何形状，在分岔点解不唯一，求解器无法捕捉屈曲。在非线性分析中，基于“现实中不存在完美的结构”这一前提，加入初始缺陷是常规做法。实际结构也存在制造误差引起的微小凹凸，所以加入初始缺陷的模型反而更真实。
求解器设置要点
求解器设置有什么需要注意的要点吗？
| 设置项目 | 推荐值 | 理由 |
|---|---|---|
| 求解的特征值数量 | 10〜20个 | 不仅确认最低阶模态，也确认高阶模态的影响 |
| 特征值求解器 | Lanczos法 | 对于大规模模型稳定且快速 |
| 温度参考温度 | 无应力温度 | 设置错误导致全部应力偏移的事故很多 |
| NLGEOM（非线性） | ON | 追踪后屈曲时必须开启 |