传热学基础 — 导热、对流、辐射与CAE热分析完整解析
老师，我做芯片散热仿真，软件里有好多热边界条件选项——"热流"、"对流系数"、"辐射"，这些分别对应什么物理？
这就是传热学的三种基本方式：导热（固体内部的热传递）、对流（流体流动带走热量，就是你说的"对流系数"）、辐射（电磁波形式的热传递，不需要介质）。芯片散热通常导热和对流都有，高温设备（如炉子）还得加辐射。
1. 传热三种方式总览
| 传热方式 | 基本定律 | 驱动力 | 需要介质？ | 关键参数 |
|---|---|---|---|---|
| 导热（Conduction） | 傅里叶定律 | 温度梯度 | 是（固体/静止流体） | 导热系数 k [W/m·K] |
| 对流（Convection） | 牛顿冷却定律 | 温差 + 流动 | 是（流体） | 对流系数 h [W/m²·K] |
| 辐射（Radiation） | 斯特藩-玻尔兹曼定律 | 温度（绝对温度⁴） | 否（真空中可传播） | 发射率 ε，视因子 F |
2. 热传导——傅里叶定律与热传导方程
导热系数k具体代表什么？铜和不锈钢差多少？
导热系数k表示材料传导热量的能力——k越大，相同温差下传递的热流越大。铜的k约385 W/m·K，不锈钢约16 W/m·K，差了约24倍！这就是为什么散热器用铜，而保温杯用不锈钢。空气的k约0.026 W/m·K，比固体低很多——所以双层玻璃窗靠中间的空气层隔热。
2.1 傅里叶导热定律
热流密度（单位面积热流量，W/m²）：
负号表示热流从高温流向低温（热力学第二定律）。对于一维问题：
通过面积 $A$ 的总热流率 $\dot{Q} = q_x \cdot A$（单位：W）。
2.2 热传导方程（热扩散方程）
含内热源 $\dot{q}$（W/m³）的瞬态热传导方程：
各项含义：
- $\rho c_p \frac{\partial T}{\partial t}$：材料热容量（温度随时间变化的惯性）
- $\nabla \cdot (k \nabla T)$：热扩散（来自周围的净热流）
- $\dot{q}$：内热源（电阻发热、化学反应热、辐射吸收等）
热扩散系数 $\alpha = k / (\rho c_p)$（单位：m²/s）衡量温度变化传播的速度。$\alpha$ 越大，材料"热响应"越快。
2.3 热阻类比
热传导类比于电路的欧姆定律：$\dot{Q} = \Delta T / R_{th}$，热阻的定义：
多层结构的总热阻 = 各层热阻串联（平板）或并联。这一概念在电子散热（芯片→焊料→基板→散热器→空气的热阻链）中广泛应用。
| 材料 | 导热系数 k (W/m·K) | 密度 (kg/m³) | 比热容 cp (J/kg·K) | 热扩散系数 α (m²/s) |
|---|---|---|---|---|
| 纯铜 | 385 | 8960 | 385 | 1.12×10⁻⁴ |
| 铝合金（6061） | 167 | 2700 | 896 | 6.9×10⁻⁵ |
| 结构钢 | 50 | 7800 | 500 | 1.28×10⁻⁵ |
| 不锈钢（316） | 16 | 8000 | 500 | 4.0×10⁻⁶ |
| 碳化硅（SiC） | 120 | 3210 | 750 | 5.0×10⁻⁵ |
| 硅（半导体） | 148 | 2330 | 700 | 9.1×10⁻⁵ |
| 热界面材料（TIM） | 1～10 | — | — | — |
| 空气（20°C） | 0.026 | 1.205 | 1007 | 2.1×10⁻⁵ |
3. 对流换热——Nu、Re、Pr关系
软件里对流边界条件要输入一个h值，这个对流系数怎么确定？能直接查表吗？
h值有两种获取方式：① 从经验公式（Nu关系式）估算；② 做CFD仿真直接算出来（这样最准确）。h的量级参考：自然对流空气约5～25 W/m²·K，强制对流空气约25～250，水冷约500～10000，沸腾传热可达10000以上。
3.1 牛顿冷却定律
$h$ 是对流换热系数（W/m²·K），$T_s$ 是壁面温度，$T_\infty$ 是流体远场温度。
3.2 努塞尔数（Nu）——无量纲对流强度
$L$ 是特征长度，$k_f$ 是流体导热系数。$Nu = 1$ 相当于纯导热；$Nu \gg 1$ 说明对流大幅增强传热。
3.3 强制对流——Dittus-Boelter公式（管内湍流）
适用条件：$Re > 10000$，$0.6 < Pr < 160$，$L/D > 10$（充分发展）。
其中普朗特数 $Pr = \mu c_p / k_f$（动量扩散/热量扩散比）：空气 $Pr \approx 0.71$，水 $Pr \approx 6.9$（20°C），润滑油 $Pr \approx 100$～1000。
3.4 自然对流——Churchill-Chu公式（垂直平板）
其中瑞利数 $Ra_L = Gr_L \cdot Pr = \frac{g\beta\Delta T L^3}{\nu^2}\cdot Pr$，$\beta$ 是热膨胀系数（理想气体：$\beta = 1/T$）。适用于整个层流和湍流范围（$10^{-1} < Ra_L < 10^{12}$）。
3.5 典型对流系数参考值
| 传热场景 | h (W/m²·K) | 备注 |
|---|---|---|
| 自然对流（空气） | 5～25 | 电子设备被动散热 |
| 强制对流（空气，低速） | 25～100 | 风扇冷却，v=1～5 m/s |
| 强制对流（空气，高速） | 100～300 | 风扇冷却，v>10 m/s |
| 强制对流（水） | 500～10,000 | 水冷散热器 |
| 沸腾传热（水，大气压） | 2,000～100,000 | 池沸腾/流动沸腾 |
| 冷凝传热（水蒸气） | 5,000～50,000 | 冷凝器 |
4. 辐射换热——视因子与封闭系问题
辐射换热为什么跟温度的四次方有关？什么时候辐射不能忽略？
辐射遵循斯特藩-玻尔兹曼定律——黑体辐射功率与绝对温度的四次方成正比。所以温度越高，辐射占的比例越大。一般来说：300K以下（室温）辐射可以忽略；500K以上辐射不可忽略；1000K以上辐射常常是主要传热方式。汽车排气管、炉子内衬、卫星热控——都必须考虑辐射。
4.1 黑体辐射——斯特藩-玻尔兹曼定律
黑体（理想辐射体，发射率=1）单位面积辐射功率：
实际物体（灰体假设：发射率与波长无关）的辐射功率：
发射率 $\varepsilon \in (0,1)$：高度抛光金属 $\varepsilon \approx 0.05$，氧化铁 $\varepsilon \approx 0.8$，人体皮肤 $\varepsilon \approx 0.95$，黑漆 $\varepsilon \approx 0.97$。
4.2 两灰体表面间的净辐射热流
对于两个面积分别为 $A_1$、$A_2$ 的平行大平板（$F_{12}=1$）：
线性化（当 $T_1 \approx T_2$ 时）有时用到：$T_1^4 - T_2^4 \approx 4T_m^3(T_1-T_2)$，其中 $T_m = (T_1+T_2)/2$。
4.3 视因子（角系数）
视因子 $F_{ij}$ 表示从面 $i$ 出发的辐射中，到达面 $j$ 的比例（$0 \leq F_{ij} \leq 1$）：
互换关系：$A_i F_{ij} = A_j F_{ji}$（互反律）
封闭关系：$\sum_{j=1}^N F_{ij} = 1$（封闭系，能量守恒）
4.4 辐射视因子计算方法
- 解析公式：简单几何（平行平板、同轴圆柱）有闭合公式，查辐射换热手册
- Hottel弦法：二维问题的图解法
- 蒙特卡洛法：复杂几何的随机射线追踪，FEM软件内置
- 有限元法：与导热方程耦合求解，商业软件（ANSYS、Abaqus）支持
5. FEM热分析
热分析的有限元方程是怎么建立的？跟结构分析的Ku=f有什么相似之处？
热分析的FEM方程和结构分析很像，只是未知量变成了温度而不是位移。稳态热分析对应线性方程组 $\mathbf{K}_T \mathbf{T} = \mathbf{f}_T$，其中 $\mathbf{K}_T$ 是热传导矩阵，$\mathbf{T}$ 是节点温度向量，$\mathbf{f}_T$ 是热载荷向量。
5.1 稳态热分析方程
热传导矩阵：
热载荷向量（体热源 + 对流 + 热流）：
5.2 瞬态热分析——热容矩阵
瞬态问题需要加上热容矩阵 $\mathbf{C}_T$：
$\mathbf{C}_T = \int_\Omega \rho c_p \mathbf{N}^T \mathbf{N}\, d\Omega$。时间积分用向后差分（隐式，无条件稳定）：$\dot{\mathbf{T}} \approx (\mathbf{T}^{n+1} - \mathbf{T}^n)/\Delta t$。
5.3 热边界条件类型
| 边界条件类型 | 数学形式 | 软件中的输入 | 典型应用 |
|---|---|---|---|
| 温度约束（第一类） | $T|_\Gamma = T_0$ | 固定温度 | 与恒温体接触的边界 |
| 热流（第二类） | $-k\partial T/\partial n = q_0$ | 热流密度 (W/m²) | 电加热膜、已知热流 |
| 对流（第三类，Robin） | $-k\partial T/\partial n = h(T-T_\infty)$ | h值 + 环境温度 | 风冷、水冷散热 |
| 辐射 | $-k\partial T/\partial n = \varepsilon\sigma(T^4-T_{rad}^4)$ | 发射率 + 辐射环境温度 | 高温炉子、太空热控 |
| 绝热（零热流） | $\partial T/\partial n = 0$ | 默认（不施加其他条件） | 对称面、保温面 |
5.4 热-结构耦合分析
温度场会产生热应变，导致热应力和变形（热-结构耦合）：
其中 $\alpha_{CTE}$ 是热膨胀系数（CTE，单位：1/K）。热应力的有限元分析流程：
- 进行热分析，得到温度场 $T(x,y,z)$
- 以温度场为载荷，计算热应变 $\varepsilon_{th} = \alpha(T-T_{ref})$
- 热应变作为"初始应变"加入结构方程，求解热应力和变形
PCB板在焊接后冷却过程中，因铜（$\alpha \approx 17$ ppm/K）和FR4基板（$\alpha_{xy} \approx 16$，$\alpha_z \approx 60$ ppm/K）的CTE差异，焊点会承受很大的热应力——这就是焊点疲劳失效的主要原因。
6. 实际热分析技巧与常见错误
做芯片热分析时，模型要多精细？全部建实体模型工作量太大了。
芯片热分析的关键是热阻网络。可以先用简化模型（R_θJC + R_θCS + R_θSA）快速估算结温，再对最热路径建精细模型。一般来说：焊料层和TIM的热阻最关键（厚度薄、热阻大），优先精细建模；PCB的过孔热导等效处理（等效各向异性导热系数）就够了。
热分析常见错误
- h值选取不合理：自然对流空气5 W/m²·K和强制对流200 W/m²·K差40倍，设错了结果差很多
- 忽略接触热阻：两个固体面接触，接触界面的热阻（取决于表面粗糙度和接触压力）不可忽略，一般需要实验测量或查厂商数据
- 低温时加了辐射：低于200°C时辐射换热量很小，加了反而增加了收敛难度（辐射方程高度非线性）
- CTE方向搞错：各向异性材料（PCB板、复合材料）的CTE在x/y/z方向不同，坐标系设置一定要对
- 不检查能量守恒：求解后对比输入热功率和边界导出热量，误差应 <1%
- 稳态分析但热源是瞬态的：芯片功耗随工作状态变化，稳态分析只给出最差情况，不能反映峰值脉冲热效应
6.1 热分析精度目标
| 分析场景 | 主要热阻路径 | 可接受精度 | 验证方法 |
|---|---|---|---|
| 电子芯片结温 | 硅→焊料→基板→散热器→空气 | ±5°C | 热电偶测壳温 + 结温推算 |
| 发动机冷却系统 | 燃烧壁面→冷却液 | ±10°C | 热像仪 + 冷却液出口温度 |
| 建筑外墙热桥 | 金属框架导热 | ±0.1 W/m²·K | 热箱法 |
| 航天器热控 | 辐射+导热+内热源 | ±3°C | 热真空试验 |