热与温度基础 — 比热・热传导・对流与CAE热分析
目录
1. 热分析在电子和航空中的重要性
CPU散热器设计需要考虑哪些参数？我知道要用铜或铝，但具体怎么算？
稳态主要看热阻：$R_{th} = \Delta T / Q$，你要把热量Q从芯片导到空气，温差$\Delta T$越小越好。热阻由材料导热系数、接触面积和散热翅片效率共同决定。但瞬态很不一样——比如游戏突发负载时，铜的大热容量$mc$能"吸住"一时的热量，让结温不那么快上升。所以稳态关注导热系数，瞬态关注热容量。
热分析是CAE中仅次于结构分析的第二大应用领域。电子设备的可靠性、航空发动机的寿命、汽车电池的安全性，都直接由热场决定。"热死机"（芯片过热保护）和"热失控"（电池热失控）都是热设计失败的工程代价。
2. 温度与热量的区别
温度（Temperature）是物质热运动强度的度量，是状态量（intensive property）。热量（Heat）是热传递过程中传递的能量，是过程量（extensive property）。
| 概念 | 符号 | 单位 | 本质 |
|---|---|---|---|
| 温度 | $T$ | K 或 °C | 分子平均动能的度量 |
| 热量 | $Q$ | J | 传递的热能（过程量） |
| 热功率（热流） | $\dot{Q}$ | W = J/s | 单位时间传递的热量 |
| 热通量 | $q''$ | W/m² | 单位面积的热功率 |
| 热导率 | $k$ 或 $\lambda$ | W/(m·K) | 材料导热能力 |
温度尺度换算：
3. 比热容与热容量
物体温度升高 $\Delta T$ 所需热量：
其中 $c_p$ 是定压比热容（J/(kg·K)），$m$ 是质量（kg）。
比热容在工程热分析中的作用：
- 瞬态热分析：决定温度响应的时间尺度 $\tau \sim mc_p/hA$
- 热冲击评估：$c_p$ 大的材料对突发热载荷"缓冲"能力强
- 相变材料（PCM）：利用相变潜热 $L$ 存储大量热量（$Q = mL$）
| 材料 | 比热容 $c_p$ (J/(kg·K)) | 密度 ρ (kg/m³) | 体积热容 ρcp (MJ/(m³·K)) |
|---|---|---|---|
| 铜 | 385 | 8960 | 3.45 |
| 铝 | 900 | 2700 | 2.43 |
| 钢 | 490 | 7850 | 3.85 |
| 水 | 4182 | 1000 | 4.18 |
| 空气（常压） | 1005 | 1.2 | 0.0012 |
| 混凝土 | 880 | 2400 | 2.11 |
4. 热传导：Fourier定律
Fourier热传导定律（一维）：
负号表示热量从高温流向低温（与温度梯度方向相反）。
通过平板的热功率（稳态）：
热阻 $R_{th}$（K/W）是电路热模型（Thermal Circuit）的核心概念：多层材料的热阻串联，多通道散热的热阻并联。
三维稳态热传导方程（Fourier 偏微分方程）：
其中 $\dot{q}$ 是体积热源（W/m³），如焦耳热、核反应热。
5. 热对流：Newton冷却定律
Newton冷却定律（对流热边界条件）：
其中 $h$ 是对流换热系数（W/(m²·K)），$T_s$ 是壁面温度，$T_\infty$ 是流体远场温度。
$h$ 值取决于流动状态，差异非常大：
| 对流类型 | h (W/(m²·K)) | 说明 |
|---|---|---|
| 空气自然对流 | 2 - 25 | 无风条件下的散热 |
| 空气强制对流 | 25 - 250 | 风扇冷却 |
| 水强制对流 | 500 - 10,000 | 液冷、水冷散热器 |
| 沸腾传热 | 2,500 - 100,000 | 热管、相变冷却 |
| 液态金属强制对流 | 10,000 - 100,000 | 核反应堆冷却 |
那h值这么不确定，CAE热分析的时候怎么设置？
h值是热分析最大的不确定性来源。通常有三个途径：一是用经验关联式（如Dittus-Boelter、Churchill-Bernstein）根据流速和几何尺寸估算；二是做CFD热流体仿真，精确计算h的空间分布，然后映射到结构热分析模型上（热-流耦合）；三是试验测量（把热电偶数据反推h）。工程中CFD+FEM耦合是主流做法，尤其是发动机散热、电池包热管理这类复杂系统。
6. 热辐射：Stefan-Boltzmann定律
黑体辐射热功率：
其中 $\sigma = 5.67 \times 10^{-8}$ W/(m²·K⁴) 是Stefan-Boltzmann常数。
实际物体的辐射（考虑发射率 $\varepsilon$）：
两表面间的净辐射热交换：
其中 $F_{12}$ 是视角系数（View Factor），取决于两表面的几何关系（0~1）。
辐射传热在以下工程场景中占主导：
- 真空环境（航天器热控——无对流）
- 高温工业炉（$T > 800°C$）
- 屋顶/外墙太阳辐射（建筑节能）
- 发动机辐射热（含视角系数的复杂几何计算）
7. FEM热分析方程
瞬态热传导方程（强形式）：
FEM弱形式离散后，得到热分析矩阵方程：
其中：
- $[C_T] = \int_V \rho c_p [N]^T [N] \, dV$：热容矩阵（类比质量矩阵）
- $[K_T] = \int_V [B]^T k [B] \, dV + \int_S h [N]^T [N] \, dS$：热刚度矩阵（传导 + 对流）
- $\{Q\} = \int_V [N]^T \dot{q} \, dV + \int_S [N]^T q'' \, dS + \int_S h T_\infty [N]^T \, dS$：热载荷向量
稳态热分析（$\dot{T} = 0$）：$[K_T]\{T\} = \{Q\}$
瞬态热分析：隐式Euler法或Crank-Nicolson法时间积分
8. 工程实例：CPU散热器设计
以典型桌面处理器（TDP = 125 W，结温上限 $T_{j,\max} = 95°C$，环境温度 $T_a = 25°C$）为例：
总热阻预算：
热阻分配（热路串联）：
| 热阻段 | 热阻 (K/W) | 材料/手段 |
|---|---|---|
| 芯片结到壳 $R_{jc}$ | 0.1 ~ 0.2 | 硅芯片 + 封装（固定） |
| 壳到散热器 $R_{cs}$ | 0.05 ~ 0.15 | 导热硅脂（TIM） |
| 散热器到空气 $R_{sa}$ | 0.2 ~ 0.35 | 翅片面积 + 风扇流量 |
散热器热阻估算：
其中 $A_{\text{fin}}$ 是翅片总面积，$\eta_{\text{fin}}$ 是翅片效率（0.7~0.9）。
FEM热分析可以：
- 计算散热翅片的温度分布，找出热点
- 优化翅片几何（间距、高度、材料）
- 评估风扇失效时的应急热容量（$mc_p \Delta T$ 能撑多久）
- 进行热-结构耦合分析，预测焊点热疲劳寿命
9. 工程热物性参数对比
| 材料 | k (W/(m·K)) | cp (J/(kg·K)) | ρ (kg/m³) | 热扩散率 α = k/(ρcp) (mm²/s) |
|---|---|---|---|---|
| 纯铜 | 398 | 385 | 8960 | 115.4 |
| 纯铝 | 237 | 900 | 2700 | 97.5 |
| 结构钢 | 50 | 490 | 7850 | 13.0 |
| 不锈钢 316L | 16 | 500 | 7990 | 4.0 |
| 混凝土 | 1.7 | 880 | 2400 | 0.80 |
| 碳纤维（轴向） | 80 | 710 | 1600 | 70.4 |
| 氧化铝（Al₂O₃） | 30 | 880 | 3900 | 8.7 |
| 金刚石 | 2200 | 520 | 3500 | 1209 |
热扩散率 $\alpha$ 决定了温度变化在材料内部传播的速度——$\alpha$ 大的材料温度均匀化快。这就是为什么铜散热器比铝更快地"扩散"热点。
总结
热与温度的基础概念——Fourier热传导、Newton冷却、Stefan-Boltzmann辐射——构成了FEM热分析的物理输入。从简单的热阻模型到复杂的瞬态三维热场，工程师需要根据问题的物理特征选择合适的传热机制建模，并对h值、边界条件等不确定参数进行敏感性分析。