热膨胀与材料 — αΔT 到热应力、双金属、CAE热结构耦合分析

分类: 物理基础 | 2026-03-25 | 网站地图
NovaSolver 贡献者
CAE visualization for thermal expansion basics - technical simulation diagram
热膨胀基础

1. 热膨胀与CAE — 无法回避的现象

温度升高时物体膨胀，温度降低时物体收缩——这是再简单不过的现象，但在工程设计中却会变成无法忽视的严重问题。热应力（Thermal Stress）是在温度变化与变形拘束相互作用下产生的。设计失误会导致裂纹、变形、疲劳破坏，严重时甚至导致结构崩溃。

🧑‍🎓 学生 我看过一条新闻说夏天铁轨弯曲导致脱轨事故，那是热膨胀吗？
🎓 博士 完全正确。钢铁轨的线膨胀系数约为 $12\times10^{-6}\,\text{/K}$。一条100米长的铁轨在气温升高30℃时会伸长 $\Delta L = 12\times10^{-6}\times100\times30 = 0.036\,\text{m} = 36\,\text{mm}$。当轨道两端被固定而无法自由膨胀时，就会产生纵向座屈（横向弯曲），最坏情况下导致脱轨事故。
🧑‍🎓 学生 那铁路是如何防止这种情况的呢？
🎓 博士 现代铁路采用"长焊接轨"来最小化接缝，但同时在铺设时对轨道施加人工拉力（预应力）。这样在夏季膨胀时，拉力就会被抵消，轨道不会座屈——这称为"中性温度管理"。这个中性温度设置的任何偏差都可能导致热座屈事故，所以需要用FEM进行详细的热应力评估。

2. 线膨胀系数的定义与主要材料数据

温度变化 $\Delta T$ 引起的线性尺寸变化率称为线膨胀系数（Coefficient of Thermal Expansion: CTE） $\alpha$ [1/K 或 1/°C]。

$$\Delta L = \alpha L_0 \Delta T$$

三维扩展（各向同性材料的体膨胀）：

$$\Delta V = \beta V_0 \Delta T, \quad \beta \approx 3\alpha$$

热应变向量（各向同性材料）：

$$\{\varepsilon_{th}\} = \alpha\Delta T\{1, 1, 1, 0, 0, 0\}^T$$
材料α [×10⁻⁶/K]备注
结构用碳钢11.7铁轨、桥梁、建筑结构
不锈钢（SUS304）17.3食品、化工、医疗设备
铝合金（A6061）23.1钢的约2倍——异种接合时热应力大
铜16.5电气配线、热交换器
CFRP（纤维方向）0～−1可实现负膨胀，精密设备用
CFRP（纤维垂直方向）30～40与纤维方向差异大，易层间剥离
玻璃（硼硅酸盐）3.3实验玻璃器皿、光学部件
硅2.6半导体芯片——与PCB的CTE不匹配是问题
因瓦合金（Fe-Ni 36%）~1.2精密设备、LNG储罐
混凝土10～12与钢筋的匹配性至关重要

3. 自由热膨胀与拘束热应力

当热膨胀自由进行（无拘束）时，会产生热应变 $\varepsilon_{th} = \alpha\Delta T$，但应力为零。

当处于完全拘束状态（伸长被完全阻止）时，热膨胀转化为弹性压缩应变。

$$\sigma_{th} = -E\alpha\Delta T$$

负号表示：温度上升（$\Delta T > 0$）时产生压缩应力（$\sigma < 0$）。

固定在桥墩间的钢梁夏季热应力（例题）

钢梁（$E = 200\,\text{GPa}$、$\alpha = 12\times10^{-6}\,\text{/K}$）在冬季0℃时安装，夏季升温至40℃时的完全拘束热应力：

$$\sigma_{th} = -200\times10^9 \times 12\times10^{-6} \times 40 = -96\,\text{MPa}\,\text{（压缩）}$$

SS400钢的屈服应力为245 MPa，安全系数约为2.5——在这种情况下是安全的。但还需考虑自重和活荷载等外部荷载，需要评估总应力。

部分拘束的情况

实际工程结构几乎都是部分拘束。使用拘束度 $R$（0=自由，1=完全拘束）来表示：

$$\sigma_{th} = -R \cdot E\alpha\Delta T$$

FEM会根据边界条件和周围结构的刚度自动计算出合适的热应力。

4. 双金属的弯曲变形

当两种膨胀系数不同的金属板被粘合在一起时，就形成了双金属（Bimetal）。温度变化时，两层间膨胀差异导致结构自然弯曲。

🧑‍🎓 学生 双金属在温度计和自动温度断路器中使用。CAE也需要分析吗？
🎓 博士 在电子封装领域至关重要——半导体芯片安装在印刷电路板（PCB）上。硅芯片（$\alpha \approx 2.6\times10^{-6}$/K）和FR4 PCB（$\alpha \approx 18\times10^{-6}$/K）的膨胀系数差异巨大。工作时（高温）和室温之间反复的温度循环会在焊点处产生循环剪切应力，导致疲劳破坏（热疲劳）。
🧑‍🎓 学生 手机和汽车电子元件的故障就是这个原因吗？
🎓 博士 正是如此。汽车发动机舱内的电子元件要承受 −40℃（严寒地区冬季）到 +125℃（发动机运行时）的反复温度循环。按照IEC 60068标准进行热循环测试，使用FEM模拟焊点的剪切应变，然后用Coffin-Manson公式预测疲劳寿命，这已成为可靠性设计的标准方法。

Timoshenko双金属曲率公式

两层双金属（厚度 $t_1, t_2$、杨氏模量 $E_1, E_2$）在温度变化 $\Delta T$ 时的曲率 $\kappa$（近似公式）：

$$\kappa \approx \frac{6(\alpha_2 - \alpha_1)\Delta T}{(t_1 + t_2)\left[3(1+m)^2 + (1+mn)\left(m^2 + \frac{1}{mn}\right)\right]}$$

其中 $m = t_1/t_2$（厚度比）、$n = E_1/E_2$（刚性比）。对于简化情形（相同厚度和杨氏模量）：

$$\kappa \approx \frac{3(\alpha_2 - \alpha_1)\Delta T}{2t}$$

曲率半径 $R = 1/\kappa$——温度变化越大、板越薄，弯曲越明显。

5. 异种材料接合部的热应力

当CTE（线膨胀系数）不同的材料被接合在一起时，温度变化会在接合面产生剪切应力和分离应力集中。

接合面剪切应力分布

CTE差异 $\Delta\alpha = \alpha_1 - \alpha_2$。温度变化 $\Delta T$ 时，"自由膨胀差"为 $\Delta\alpha \cdot \Delta T \cdot L$（$L$为长度）。这种差异在接合层产生剪切变形。

剪切应力集中在接合面端部（端部应力集中）。可通过接合面圆角设计或底部填充树脂来缓解应力。

焊点热疲劳寿命（Coffin-Manson公式）

$$N_f = C \left(\frac{\varepsilon_{f}'}{2\Delta\varepsilon_p/2}\right)^{1/c}$$

$\Delta\varepsilon_p$：单个循环的塑性剪切应变振幅（由FEM获得）
$\varepsilon_f'$、$c$：材料疲劳常数（铅焊料 $c \approx -0.6$～-0.7）

SAC305（Sn-Ag-Cu）无铅焊料在约125℃循环条件下的寿命目标为数百到数千个循环。FEM计算出一个循环的塑性应变分布后，用此公式评估每个焊点的寿命。

6. 热膨胀的温度依赖性与相变

线膨胀系数 $\alpha$ 并非常数，而是随温度变化。高精度的FEM分析需要使用温度相关的材料数据表。

铁的相变与体积变化

铁和钢在特定温度下发生晶体结构转变（相变），体积会急剧变化。

  • $\alpha$-Fe（体心立方: BCC）↔ $\gamma$-Fe（面心立方: FCC）转变（~912℃）：约1%的体积收缩
  • $\gamma$-Fe → 马氏体转变（淬火）：体积膨胀（~0.2～0.4%）

淬火时的不均匀冷却（表面和内部温度差异）与相变的不均匀相结合，会产生巨大的残余应力。弯曲和开裂成为淬火件的常见问题，使用包含相变的热结构耦合FEM（如Abaqus、DEFORM等）来预测变形成为必须。

混凝土的温度依赖性

混凝土的碱硅反应（ASR）会因吸水产生膨胀，高温时的膨胀行为也很重要。核电站格纳容器混凝土在150℃以上会出现强度下降和膨胀增大——极限事故分析时需要考虑这种非线性特性。

7. FEM热结构耦合分析

当温度分布不均匀（存在温度梯度）时，热应力计算必须分两步进行。

逐步耦合分析的步骤

  1. 热分析（第1步）：设定边界条件（热源、冷却、绝热）计算温度场 $T(x, y, z)$。FEM热分析方程：$[K_T]\{T\} = \{Q\}$
  2. 结构分析（第2步）：将第1步的温度分布作为输入，热应变 $\varepsilon_{th} = \alpha(T - T_{ref})\delta_{ij}$ 作为"荷载"输入到结构分析中，计算位移和应力。

两个步骤可以使用不同的有限元（热单元DC3D8和结构单元C3D8等），但使用相同网格可避免温度数据插值，提高精度。

参考温度（Stress-free Temperature）的设置

最容易被忽视的设置是"参考温度（$T_{ref}$）"——这是制造后应力为零的温度，包括焊接温度、焊后冷却温度、成形温度等。Abaqus中通过 *INITIAL CONDITIONS, TYPE=TEMPERATURE 设置。$T_{ref}$ 的错误设置会导致热应力计算严重偏差。

温度相关材料性质的设置

高温分析时，杨氏模量、屈服应力、热导率、比热都是温度的函数。Abaqus中杨氏模量的温度表输入示例：

*ELASTIC
200e9, , 20.0
190e9, , 200.0
170e9, , 400.0
140e9, , 600.0

8. 热膨胀的控制设计

精密设备（测量仪器、光学系统、空间设备）中，最小化温度变化导致的尺寸变化是性能的关键。

低膨胀材料

  • 因瓦合金（Fe-36%Ni）：$\alpha \approx 1.2\times10^{-6}$/K（钢的1/10）。用于LNG储罐、测长机、望远镜筒。
  • Zerodur（玻璃陶瓷）：$\alpha \approx 0.05\times10^{-6}$/K。用于天文望远镜主镜。
  • CFRP（纤维方向）：通过纤维配向设计可达到 $\alpha = 0$ 或负值。用于宇宙空间结构部件和精密光学架台。

精密设备的热膨胀补偿

三坐标测量机（CMM）的测量臂在温度变化时长度会改变，影响测量精度。使用因瓦合金和温度传感器的实时补偿相结合，可在±0.001℃的温度管理下实现亚微米级精度。

9. 实践：排气歧管的热应力疲劳

汽车和摩托车的排气歧管（排气多歧管）在发动机启动和停止时经历最剧烈的温度循环，是热疲劳最恶劣的工作环境之一。

热工况

  • 发动机启动时：常温（约20℃）快速升温到排气温度（700～900℃）
  • 发动机停止后：在数十分钟内自然冷却至常温
  • 运行中：平均600～800℃，局部可超过1000℃

材料与破损模式

材料最高使用温度 [℃]主要破损模式
SUS430（铁素体不锈钢）~900热疲劳裂纹（焊缝部）
SUS304（奥氏体不锈钢）~870高温蠕变+热疲劳
铸铁（Si-Mo铸铁）~800裂纹扩展、氧化减薄
Inconel 625（镍合金）~1100成本高但高性能，用于赛车

FEM分析流程

  1. 输入启动～定常～停止的温度历程（CFD或实测数据）
  2. 热结构耦合分析计算各时刻的应力应变分布
  3. 评估塑性应变累积（高温下材料易屈服）
  4. 从塑性应变振幅 $\Delta\varepsilon_p$ 用Coffin-Manson公式估算疲劳循环数
  5. 与目标寿命（例如10万km行驶）比较，优化形状和材料

交叉主题

热传递 热分析 耦合分析 应力和应变 热和温度 结构分析