熱-構造連成設計（Thermo-Mechanical Design）
理论与物理
什么是热-结构耦合设计
热-结构耦合设计，不就是计算热应力吗？
不，热应力评估只是其中一个方面。更重要的是热变形导致的尺寸精度劣化。例如半导体光刻设备的载物台需要0.1nm的精度，但温度波动哪怕只有0.01°C，也会产生数nm的变形。
0.01°C就产生数nm！？那已经是温度管理的范畴了…
是的。因此，从因瓦合金（Fe-Ni低膨胀合金，CTE ≈ 1.2×10⁻⁶/K）等特殊材料的选型、冷却流道的最优布置、到约束条件的设计，这些都包含在“热-结构耦合设计”中。不仅仅是计算 $\sigma = E \alpha \Delta T$ 就结束了。
从材料选型到冷却设计全都包含啊。具体在哪些产业中使用呢？
主要有这5个领域：
- 半导体制造设备：光刻载物台、晶圆卡盘（nm级精度）
- 航空发动机：涡轮叶片、燃烧室衬套（1600°C环境）
- 电子设备：BGA/CSP封装、功率模块（温度循环可靠性）
- 汽车：排气歧管、刹车盘
- 工厂设备：蒸汽轮机、压力容器（启停循环）
热变形与尺寸精度
热变形问题那么严重吗？感觉铁很硬，不容易变形啊。
看数字就知道有多可怕了。例如，一根长1m的钢材（CTE ≈ 12×10⁻⁶/K）温度升高10°C时：
0.12mm——对于普通机械设计可能没问题，但如果机床的加工精度是μm级，那就是致命的。而且温度分布不均匀，产生的不是简单的膨胀，而是“翘曲”或“扭曲”。这才是热变形真正可怕的地方。
如果是均匀膨胀还能补偿，但翘曲和扭曲就很棘手了…
没错。所以在设计阶段，需要从两方面考虑：“如何减小温度分布的不均匀性”和“如何承受残留的不均匀性”。前者是冷却设计，后者是材料选型和约束设计。
CTE不匹配的力学
CTE不匹配，是指不同材料接合时的问题吧？
是的。例如在电子基板中，硅芯片（CTE ≈ 2.6×10⁻⁶/K）通过焊料连接到FR-4基板（CTE ≈ 14×10⁻⁶/K）上。温度循环时，CTE差异直接转化为界面剪切应力。
这里 $\gamma$ 是焊料的剪切应变，$L$ 是芯片的半宽，$h_{solder}$ 是焊料高度。CTE差异越大、芯片越大、焊料越薄，剪切应变就越大。这就是BGA焊点开裂的根本原因。
车载ECU的温度循环测试有“-40°C〜+125°C，1000次循环”这样的规格，正是在评估这个吧。
很敏锐。涡轮叶片也存在同样的问题。在镍基高温合金上喷涂TBC（热障涂层，陶瓷系，CTE ≈ 10×10⁻⁶/K），但镍基高温合金的CTE约为13×10⁻⁶/K。这种差异在反复启停过程中会导致TBC剥落。
控制方程
数学公式是怎样的？热和结构的方程是耦合的吧？
首先是热传导方程（瞬态）：
这里 $\rho$ 是密度，$c_p$ 是比热容，$k$ 是热导率，$Q$ 是内部热源。然后是结构的平衡方程：
耦合的关键在于本构关系。从总应变中减去热应变得到弹性应变：
也就是说，先求解温度场，然后用这个温度求解结构。温度确定了，热应变就确定了，如果有约束就会产生热应力。
基本是这样。但是，如果材料的弹性模量 $E(T)$ 或屈服应力 $\sigma_y(T)$ 是温度相关的，那么结构侧的行为也会依赖于温度，此时即使是“单向耦合”也能准确反映温度的影响。更进一步，如果结构变形改变了接触状态从而改变了热阻（例如：螺栓连接处的接触热传导），则需要双向耦合。
瞬态热应力 vs 稳态热应力
只看稳态热应力不行吗？稳态不是最严苛的条件吗？
这是一个很大的误区。瞬态应力比稳态应力大的情况非常多。以蒸汽轮机转子为例。启动时高温蒸汽冲击外表面，外表面迅速升温但内部仍然较冷。此时外表面产生压缩应力，内表面产生拉伸应力，形成很大的温度梯度。
达到稳态后，温度梯度缓和，应力下降。也就是说，启动和停止时的瞬态过程才是最严苛的条件。实际上，蒸汽轮机制造商必定会对每种启动模式（冷启动/温启动/热启动）进行瞬态分析，以确定允许的升温速率（°C/min）。
如果只看稳态就判断“没问题”，那可能在启动时就损坏了。真可怕…
所以，在热-结构耦合设计中，对整个生命周期（启动→稳态→停止→再启动）的瞬态分析是必须的。这是它与单纯的“热应力计算”本质上的不同点。
各项的物理意义
- 热应变 $\varepsilon_{th} = \alpha \Delta T$：温度变化引起的自由膨胀。若无约束则不产生应力。铁轨接头留有缝隙就是为了释放热膨胀。
- CTE不匹配应力：异种材料接合处，膨胀量差异受到约束而产生。集中在界面，是剥离或裂纹的起点。
- 温度相关材料特性 $E(T)$, $\sigma_y(T)$：高温下弹性模量和屈服应力会下降。600°C的钢强度只有室温的约60%。
- 瞬态温度梯度 $\partial T / \partial r$：急加热/急冷却时，厚壁部件的内外会产生大的温差。毕渥数越大越显著。
假设条件与适用范围
- 线性弹性假设：高温下发生塑性变形时，需要弹塑性/蠕变本构关系
- 小变形假设：薄壁结构产生大热变形时，需要几何非线性
- 各向同性材料假设：复合材料或CFRP需考虑各向异性CTE（面内与厚度方向可差一个数量级）
- 热与结构的时间尺度分离：通常热响应比结构响应慢，因此准静态结构分析适用
量纲分析与单位制
| 物理量 | SI单位 | 代表值备注 |
|---|---|---|
| 线膨胀系数 $\alpha$ | 1/K（= μm/m/K） | 钢: 12, Al: 23, Cu: 17, Si: 2.6, 因瓦合金: 1.2 |
| 热导率 $k$ | W/(m·K) | Cu: 400, Al: 237, 钢: 50, 陶瓷: 2-30 |
| 比热容 $c_p$ | J/(kg·K) | 钢: 500, Al: 900, Cu: 385 |
| 弹性模量 $E$ | GPa | 钢: 200 (室温), 170 (500°C), 140 (700°C) |
半导体光刻设备 — 0.01°C的较量
EUV光刻设备的晶圆载物台，要求nm级的定位精度。主要结构件采用因瓦合金或SiC陶瓷，正是因为它们的CTE（线膨胀系数）极小。因瓦合金的CTE约为1.2×10⁻⁶/K，是钢的1/10。此外，载物台周围通过温控空气或He气氛控制在0.01°C以下，并由激光干涉仪进行数百Hz、pm分辨率的位移反馈控制。即便如此，热变形依然存在，因此需要用FEM预测瞬态温度分布，并预先计算控制补偿值。这是CAE“直接决定产品动作精度”领域的极致体现。
数值解法与实现
耦合策略：弱耦合 vs 强耦合
热和结构是同时求解吗？还是分开求解？
主要有三种方法：
- 单向耦合（One-way）：热分析→温度传递给结构→结构分析。适用于结构变形不影响温度场的情况。实际工作中80%的情况用这个就够了。
- 弱耦合（Sequentially coupled）：热→结构→（根据结构变形更新热条件）→热→… 循环。适用于接触热传导依赖于变形的情况。
- 强耦合（Fully coupled）：在一个矩阵中同时求解热和结构。适用于结构变形速度对发热有贡献的情况（塑性加工、摩擦发热）。
单向耦合就能解决80%啊。那强耦合在什么时候用呢？
典型例子是锻造或搅拌摩擦焊（FSW）。塑性变形产生的热量 $Q_{plastic} = \eta \cdot \sigma : \dot{\varepsilon}^p$（$\eta$ 是塑性功的热转换系数，通常为0.9）使温度升高，温度改变材料特性进而改变变形行为——形成了完整的反馈。Abaqus的 `*Coupled Temperature-Displacement` 就是处理这种强耦合的代表性分析步。
| 耦合类型 | 精度 | 计算成本 | 典型用途 |
|---|---|---|---|
| 单向 | ○（多数情况足够） | 低 | 电子设备的温度循环、稳态运行机械 |
| 弱耦合 | ◎ | 中 | 螺栓连接处的接触热传导、垫片 |
| 强耦合 | ◎ | 高 | 锻造、FSW、刹车摩擦、切削加工 |
FEM公式化
在FEM中如何引入热应变？
结构的有限元方程是：
这里热载荷向量 $\{\mathbf{F}_{th}\}$ 由每个单元的热应变求得：
$\mathbf{B}$ 是应变-位移矩阵，$\mathbf{D}$ 是弹性矩阵。也就是说，通过“温差 → 热应变 → 等效节点力”这样的转换，热效应以与机械载荷相同的形式进入结构方程的右侧。这样就能用同一个求解器求解，实现起来很方便。
原来如此，把热载荷转换成等效的力加到右边。这样就能直接使用现有的结构求解器了。
瞬态分析的时间积分
瞬态热分析的时间步长怎么确定？
对热传导方程进行FEM离散化，使用后向欧拉法时：
$[\mathbf{C}]$ 是热容矩阵，$[\mathbf{K}_T]$ 是热传导矩阵。后向欧拉法（隐式法）是无条件稳定的，所以可以取较大的时间步长，但要捕捉瞬态温度峰值，需要足够小的时间步长。参考标准是：
- 急加热初期：$\Delta t \leq L_e^2 / (6 \kappa)$（$\kappa = k/(\rho c_p)$ 是热扩散率，$L_e$ 是最小单元尺寸）
- 接近稳态后：增大 $\Delta t$ 以提高计算效率（自动时间步长控制）
初期细密，接近稳态后粗放。很合理。
热-结构间的数据传递
如果热分析和结构分析的网格不同怎么办？
这是实际工作中非常重要的点。热分析网格粗一些也能收敛，但结构分析在应力集中部位需要细密的网格。网格不一致时，需要进行温度数据的传递（映射）：
- 形函数插值：查找结构侧节点位于热侧哪个单元内，用形函数插值。最常用。
- 最近邻法：直接使用最近节点的值。简便但精度低。
- RBF插值：用径向基函数重构平滑的温度场。对网格不一致的鲁棒性强。
なった
詳しく
報告