热膨胀与材料 — αΔT热应变与热-结构耦合分析
1. 为什么热膨胀在工程中如此重要
铁路钢轨夏天会膨胀多少？这会有什么影响？
算一下：100m钢轨，温差50°C（冬-15°C到夏+35°C），线膨胀系数α=12×10⁻⁶/°C：伸长量 = α×L×ΔT = 12e-6×100×50 = 60mm。如果两端完全固定不让它膨胀，就会产生约120MPa的压应力，接近钢的屈服强度的一半。旧式钢轨每隔几十米留一个2cm伸缩缝，现代无缝钢轨则通过铺设时施加预拉应力来抵消热压应力，并靠强固道钉约束防止横向屈曲（胀轨）。
热膨胀是工程设计中经常被初学者忽视，却往往导致严重问题的物理现象。桥梁伸缩缝、发动机配合间隙、电子封装可靠性、输油管道补偿器……都是针对热膨胀进行设计的结果。
2. 线膨胀系数（CTE）
线膨胀系数（Coefficient of Thermal Expansion，CTE 或 $\alpha$）定义：
对于各向同性材料，$\alpha$ 在三个方向相同。各向异性材料（CFRP、单晶）则需要三个方向的CTE。
体膨胀系数（体积膨胀）：
3. 热应变与热位移
温度变化 $\Delta T$ 引起的线应变（自由膨胀，无约束）：
长度为 $L$ 的构件自由伸长量：
三维热应变（各向同性材料）：
热应变不产生应力（自由膨胀时），但当热膨胀受到约束时，会产生弹性应变，进而产生热应力。
4. 约束热应力
两端完全固定的杆件，升温 $\Delta T$ 后产生的热压应力：
负号表示升温产生压应力（膨胀被阻止）。对于钢（$E = 210$ GPa，$\alpha = 12 \times 10^{-6}$/°C）：
即温度每升高1°C，产生约2.52 MPa的压应力。升温50°C就是-126 MPa，升温100°C达到-252 MPa（已超过Q235钢屈服强度！）。
5. 双材料系统热应力
当两种CTE不同的材料粘接在一起（如电子芯片贴在PCB上），温度变化会在界面处产生剪切应力，并导致整个结构弯曲（双金属效应）：
其中 $\kappa$ 是曲率，$h_1, h_2$ 是两层厚度，$\alpha_1, \alpha_2$ 是CTE，$E_1, E_2$ 是弹性模量。
界面剪切应力的简化模型（CTE失配引起的剪切应力）：
这种界面剪切应力是焊锡点疲劳失效的根本原因：每次温度循环都是一次剪切疲劳循环。
手机里的芯片那么小，CTE不匹配能有多大问题？平时手机还挺耐用的。
CTE失配是消费电子最大的可靠性挑战之一。手机从低温（-20°C存放）到工作温度（60°C以上）温差达80°C，硅芯片CTE约2.5 ppm/°C，FR4 PCB约18 ppm/°C——差了7倍！每次冷热循环，焊锡球就承受一次剪切应力循环。消费电子一般要求10年内热循环不失效（几千次循环），这就是无铅焊锡（SAC305等）、底部填充（Underfill）工艺开发的动力——填充材料降低了整体CTE失配，大幅提高焊点寿命。
6. 材料热膨胀系数对比
| 材料 | CTE α (ppm/°C = 10⁻⁶/°C) | 适用温度范围 | 工程备注 |
|---|---|---|---|
| 因瓦合金（Invar, Fe36Ni） | 1.2 | -100 ~ 200°C | 精密仪器、激光镜座 |
| 硅（Si） | 2.6 | 室温 | 芯片基底 |
| 氧化铝（Al₂O₃陶瓷） | 7.5 | 室温 | 陶瓷基板、绝缘体 |
| 铸铁 | 10.5 | 室温 | 发动机缸体 |
| 结构钢 | 11 ~ 13 | 室温 | 建筑、桥梁 |
| 不锈钢 304 | 17.2 | 室温 | 化工管道 |
| 铝合金 | 22 ~ 24 | 室温 | 航空、汽车轻量化 |
| 铜 | 17 | 室温 | 电气导体 |
| FR4 PCB（面内） | 14 ~ 18 | 室温 | 电子封装 |
| 碳纤维（轴向） | -0.5 ~ 0 | 室温 | 近零膨胀（精密结构） |
碳纤维复合材料在纤维方向的CTE接近零甚至为负，与基体（CTE约50 ppm/°C）的极端差异，使得CFRP的热分析比金属复杂得多。
7. FEM热-结构耦合分析流程
热-结构耦合分析（Thermal-Structural Coupled Analysis）有两种方式：
顺序耦合（Sequential Coupling）——最常用：
- 第一步：FEM热分析，计算温度场 $\{T\}$
- 第二步：将温度场作为"热载荷"映射到结构FEM模型
- 第三步：FEM结构分析，计算热应力 $\{\sigma\} = [D](\{\varepsilon\} - \{\varepsilon_{th}\})$
热应变从弹性应变中扣除，只有弹性应变产生应力：
完全耦合（Fully Coupled）——适用于热应力影响热传导时：
温度和位移同时求解，适用于橡胶大变形（热-力耦合）、高速撞击（绝热剪切热）等情况。
FEM热应力分析注意事项：
- 参考温度（Stress-free temperature）设置至关重要——通常取制造温度（焊接温度）
- 温度相关的材料属性（$E(T)$、$\alpha(T)$）对高温工况影响显著
- 蠕变（高温长时）需要专门的蠕变本构
8. 工程实例：铁路无缝钢轨
无缝钢轨（CWR, Continuous Welded Rail）消除了传统有缝钢轨的接缝冲击，提高了乘坐舒适性和行驶速度，但必须解决热膨胀问题：
热力分析：轨温年变化约 $\Delta T = 50 ~ 70°C$（中国华北地区），100m钢轨自由膨胀量：
通过铺设时张拉至"锁定轨温"（比最低温高约30°C），使钢轨在最高温时受约束压应力，最低温时受拉应力，但两者均在安全范围内。FEM分析需要计算：
- 夏天最高温时的轴向压应力（防胀轨）
- 冬天最低温时的轴向拉应力（防断轨）
- 曲线段的横向稳定性（弯道更容易胀轨）
9. 电子封装焊点热疲劳
BGA（球栅阵列）焊点热疲劳寿命预测——Coffin-Manson关系的热版本：
其中 $\Delta\gamma$ 是每个温度循环的剪切应变范围，$\varepsilon_f'$ 和 $c$ 是焊锡疲劳参数。
工程设计目标：
| 应用场景 | 温度循环条件 | 要求寿命（循环） |
|---|---|---|
| 消费电子（手机） | 0 ~ 70°C，每天2次 | 2,000 ~ 5,000 次 |
| 汽车电子（车舱） | -40 ~ 85°C，AEC-Q100 | 1,000 次（TS-C级） |
| 汽车发动机舱 | -40 ~ 125°C，AEC-Q100 | 500 次（TS-A级） |
| 航空电子 | -55 ~ 125°C，MIL-STD-810 | 200 次 |
总结
热膨胀是机械工程和电子工程中最普遍的热-机械耦合问题。从铁路钢轨的宏观热变形管理到芯片焊点的微观热疲劳，核心公式 $\varepsilon_{th} = \alpha\Delta T$ 和约束热应力 $\sigma_{th} = -E\alpha\Delta T$ 提供了工程快速估算的基础。FEM热-结构顺序耦合分析是精确评估复杂热应力场的标准方法。