熱座屈解析
理论与物理
什么是热屈曲
老师，会因为热而发生屈曲吗？即使没有施加外力？
问得好。热屈曲（thermal buckling）是由温度变化产生的压缩应力引起的屈曲。如果结构可以自由膨胀，就不会产生温度应力，但当膨胀受到约束时，就会产生压缩应力。
铁轨在炎热天气下弯曲也是热屈曲吗？
正是如此。铁轨被枕木约束了轴向膨胀。温度上升 $\Delta T$ 产生的轴向应力为：
其中 $\alpha$ 是线膨胀系数。对于钢材（$\alpha \approx 12 \times 10^{-6}$ /°C, $E = 200$ GPa），当 $\Delta T = 40$ °C时，$\sigma_{th} = 96$ MPa。当这个应力超过屈曲应力时，铁轨就会横向弯曲。
40°C的温差就能产生96 MPa…比想象的要大呢。
是的。热应力与约束度成正比。完全约束时为 $E\alpha\Delta T$，但如果允许部分膨胀，应力就会降低。正确评估实际结构的约束度是热屈曲分析的第一步。
热屈曲的控制方程
热屈曲的公式化与机械屈曲不同吗？
特征值屈曲的框架是相同的。但参考载荷变为温度载荷：
1. 给定温度分布 $\Delta T(x,y,z)$ 进行热应力分析
2. 用得到的热应力构成 $[K_\sigma]$
3. 求解特征值问题 $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$
$\lambda$ 是屈曲温度载荷系数。$\lambda \cdot \Delta T$ 即为临界温升。
也有温度分布不均匀的情况吧。
在实际结构中，温度分布通常是不均匀的。例如，火灾时的梁，底面高温，顶面相对低温。板单面受热时，还会产生温度梯度引起的弯矩。
温度的影响有两个：
- 膜应力（面内的压缩/拉伸） — 屈曲的驱动力
- 弯矩（板厚方向的温度梯度） — 诱发附加变形
板的热屈曲
板的热屈曲有理论解吗？
四边约束的矩形板均匀受热时的临界温升：
板厚 $t$ 与宽度 $b$ 之比的平方包含在内。板越薄，在越低的温升下就会屈曲。
而且这个公式里没有杨氏模量 $E$！在理想的约束条件下，临界温度仅取决于几何参数（$t/b$）和线膨胀系数 $\alpha$，而与材料的刚度无关。
与材料无关…这有点反直觉。
完全约束时，温度应力为 $E\alpha\Delta T$，屈曲应力为 $E \cdot f(t/b)$，两者相除 $E$ 就消掉了。但实际结构并非完全约束，所以 $E$ 的影响会体现出来。
实际结构中的热屈曲问题
请告诉我一些热屈曲成为问题的实际例子。
代表性的例子：
| 结构 | 情况 | 特征 |
|---|---|---|
| 铁路钢轨 | 夏季温度上升 | 连续焊接钢轨的横向屈曲 |
| 管道 | 高温流体输送 | 海底管道的侧向屈曲（lateral buckling） |
| 飞机蒙皮 | 超音速飞行时的气动加热 | 薄板面板的屈曲 |
| 空间结构 | 向阳/背阴温差 | 太阳能电池板、天线 |
| 焊接结构 | 焊接时的局部加热 | 焊接变形（角变形、屈曲变形） |
| 火灾时的建筑物 | 钢梁受热 | 端部受约束梁的轴力屈曲 |
海底管道的屈曲很有意思呢。
海底管道的热屈曲是石油和天然气行业的一个重要问题。输送高温原油时，管道试图膨胀，但受到海底摩擦力的约束。当温度上升到一定程度以上时，管道会向侧向蛇行（lateral buckling）。有一种称为“designed lateral buckling”的设计方法，旨在有意识地控制这种屈曲。
总结
我来整理一下热屈曲的理论。
要点：
- 热屈曲是由“受约束的热膨胀”引起的压缩应力导致的 — 自由膨胀则不会屈曲
- $\sigma_{th} = E\alpha\Delta T$ — 完全约束时的热应力
- 约束度的评估是关键 — 实际结构既非完全约束也非完全自由
- FEM中采用热应力分析→特征值屈曲分析的两步法 — 即使温度分布不均匀也可应对
- 温度梯度同时产生膜应力和弯矩 — 单面加热尤其危险
“即使不施加力也会屈曲”是热屈曲可怕的地方。温度管理直接关系到结构安全性。
没错。热屈曲的“载荷”是看不见的，因此在设计中容易被忽视。在温度变化大的环境中使用的结构，应始终考虑热屈曲的可能性。
热屈曲与双金属片
热引起的屈曲发生在压缩热应力达到屈曲载荷时。1780年代由Johann Christoph Baugarten（德国）实验的“双金属片”（热膨胀系数不同的两层金属板）正是热屈曲的原型。基于此原理工作的恒温器在1900年代工业化，至今仍用于电力和空调系统的自动切断。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
基于FEM的热屈曲分析
请告诉我用FEM分析热屈曲的步骤。
分为2步（或3步）流程。
Step 1: 确定温度场
通过稳态热传导分析或瞬态热传导分析求得温度分布 $T(x,y,z)$。如果是简单的均匀加热，也可以手动给定温度。
Step 2: 热应力分析
将得到的温度场输入结构分析，计算热膨胀引起的应力分布。此阶段的边界条件（约束条件）对结果影响很大。
Step 3: 特征值屈曲分析
根据Step 2的热应力构成几何刚度矩阵，求解特征值问题。$\lambda$ 即为临界温度载荷系数。
与机械载荷复合的情况呢？
当机械载荷和热载荷同时作用时，将机械载荷作为预加载荷施加，评估其对附加温度载荷的屈曲。当温度和机械载荷独立变化时，需要注意特征值屈曲的“比例载荷”假设可能不成立。
Nastran
```
SOL 105
CEND
SUBCASE 1
TEMPERATURE(LOAD) = 100
SPC = 1
METHOD = 10
```
使用TEMPERATURE卡片指定温度分布。使用TEMP/TEMPD卡片定义节点温度。正确设置MAT1中的热膨胀系数（A）。
相关主题
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
なった
詳しく
報告