1. 什么是热分析？
热分析（Thermal Analysis）是利用数值方法求解物体内部及其与环境之间的热传递问题，预测温度场、热流密度及由温度变化引起的热应力。在现代工业产品设计中，几乎每一类产品都存在热问题：电子元器件因过热而失效、发动机零件因温度梯度产生热疲劳、建筑围护结构的隔热性能……热分析正是解决这些问题的核心工具。
1.1 稳态热分析 vs. 瞬态热分析
系统达到热平衡后温度不再随时间变化。求解 $\nabla \cdot (k\nabla T) + Q = 0$。适用场景：电子设备正常运行工况、建筑热损失计算、散热器稳定工作状态。计算代价低，是热分析的入门起点。
温度随时间变化的动态过程。求解完整的热扩散方程含时间项。适用场景：焊接过程热循环、发动机启动瞬间温升、PCB热冲击测试、相变凝固过程。计算量远大于稳态。
1.2 线性热分析 vs. 非线性热分析
当材料热导率 $k$、比热容 $c_p$、密度 $\rho$ 不随温度变化，且边界条件为线性时，可使用线性热分析，方程组直接求解，速度极快。然而在大多数真实工程问题中：
- 材料性能随温度显著变化（如钢铁导热率从室温到1000°C变化显著）
- 辐射热流与温度的四次方成正比（高度非线性）
- 相变过程（熔化/凝固）引入潜热，导致非线性
- 接触热阻随接触压力变化
这些情况都需要使用非线性热分析，通过Newton-Raphson迭代求解。
2. 三种传热方式及其控制方程
2.1 热传导（Heat Conduction）— 傅里叶定律
热传导是固体内部（有时也包括静止流体）因温度梯度驱动的热量传递方式。微观上，它是由分子振动或自由电子的碰撞传递能量的过程。
热流密度矢量 $\mathbf{q}$（单位：W/m²）与温度梯度的关系：
$$\mathbf{q} = -k \nabla T$$其中：$k$ 为热导率（W/m·K），$\nabla T$ 为温度梯度（K/m）。负号表示热量从高温流向低温。对于各向异性材料，$k$ 为张量。
由此推导出稳态热传导控制方程：
$$\nabla \cdot (k \nabla T) + Q = 0$$瞬态（含时间项）热传导方程，即热扩散方程：
$$\rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + Q$$其中：$\rho$ 为密度（kg/m³），$c_p$ 为定压比热容（J/kg·K），$Q$ 为体积热源强度（W/m³），$t$ 为时间（s）。这是整个热分析的核心方程，所有CAE热求解器本质上都在求解这个偏微分方程的数值解。
2.2 对流传热（Heat Convection）— 牛顿冷却定律
对流传热发生在固体表面与流动流体之间。在CAE中，纯热分析（不含流体计算）通常将对流作为边界条件处理，而非内部物理过程。
单位面积的对流热流量 $q_{\text{conv}}$（W/m²）：
$$q_{\text{conv}} = h(T_s - T_\infty)$$其中：$h$ 为对流换热系数（W/m²·K），$T_s$ 为固体表面温度（K），$T_\infty$ 为远处流体温度（K）。$h$ 的准确确定是工程难点，通常借助努塞尔数（Nusselt number）关联式：$h = Nu \cdot k_f / L$。
对流换热系数 $h$ 的典型范围（仅供参考）：
| 传热方式 | h 典型值 (W/m²·K) | 典型工程案例 |
|---|---|---|
| 自然对流（空气） | 2 ~ 25 | 无风扇的电子设备散热 |
| 强制对流（空气） | 25 ~ 250 | 散热风扇、风洞 |
| 强制对流（水/液冷） | 500 ~ 10,000 | 汽车发动机冷却水套 |
| 沸腾传热 | 2,500 ~ 100,000 | 核反应堆冷却、蒸发器 |
2.3 辐射传热（Thermal Radiation）— 斯蒂芬-玻尔兹曼定律
热辐射是物体因温度高于绝对零度而发射电磁波（主要是红外线）传递能量的方式。它不需要介质，在真空中也可传热，因此在航天热控、工业炉炉内传热中具有不可替代的重要性。
理想黑体单位面积辐射功率 $E_b$（W/m²）：
$$E_b = \sigma T^4$$对于实际灰体，考虑发射率 $\varepsilon$（0~1）后：
$$q_{\text{rad}} = \varepsilon \sigma (T_s^4 - T_{\text{surr}}^4)$$其中：$\sigma = 5.67 \times 10^{-8}$ W/m²·K⁴ 为斯蒂芬-玻尔兹曼常数，$T_s$ 和 $T_{\text{surr}}$ 均需使用绝对温度（K）。注意温度四次方使辐射在高温下（>500°C）占主导地位。
在有多个辐射表面的密闭空间（如工业炉膛）中，还需要引入视角因子（View Factor） $F_{ij}$，描述表面 $i$ 发出的辐射中有多少比例到达表面 $j$：
$$q_{i \to j} = \varepsilon_i \sigma F_{ij} A_i (T_i^4 - T_j^4)$$教授，我做的是PCB板热分析，客户问我用稳态就够了吗？还是必须算瞬态？
这取决于你要验证什么工况。如果是验证芯片在额定功率下连续工作时散热够不够——稳态就行，算最热的那个状态。但如果客户担心的是"突然上电瞬间温升太快、焊点热冲击"或者"系统开关机循环导致疲劳"，就必须算瞬态，才能看到温度随时间怎么变化。
那对流换热系数 $h$ 怎么确定？我看到有人直接填个10 W/m²K，有人用200，差好多啊。
对，$h$ 的取值是热分析里最容易犯错的地方之一。简单说：没有风扇的密闭机箱内自然对流，$h$ 大概5~15 W/m²K；有强制风冷的话可以到50~200甚至更高；液冷能到几千。准确的做法是：要么用CFD单独算出流场再提取 $h$（这就是共轭传热CHT的思路），要么用经验关联式，比如垂直板自然对流可以用Churchill-Chu关联式计算努塞尔数，再换算出 $h$。随便填一个数字，结果会差十倍都不止。
那辐射呢？PCB板也需要考虑辐射传热吗？
对于大多数电子设备（工作温度80°C以下），辐射可以忽略——低温下辐射量很小。但如果你做的是炉膛、烟道、航天器热控，或者工业高温炉（温度超过500°C），辐射就是主导传热机制，必须计算视角因子。记住一个规律：辐射与温度的四次方成正比，高温下辐射的占比急剧上升。
热传导方程里的 $\nabla \cdot (k \nabla T) + Q = 0$，有限元是怎么求解它的？
和结构分析的原理类似。把连续的温度场 $T(x,y,z)$ 用单元节点的温度值近似。对每个单元用形函数插值，然后用加权余量法（Galerkin法）对控制方程进行弱形式推导，最终得到一个线性代数方程组：$[K_T]\{T\} = \{F_Q\}$。这里 $[K_T]$ 叫热刚度矩阵，$\{F_Q\}$ 是热载荷向量，里面包含体热源、表面热流、对流和辐射边界条件的贡献。对于非线性问题（比如辐射），材料属性随温度变化，$[K_T]$ 随迭代更新。
瞬态热分析里时间步长怎么定？有什么规律吗？
有个简单的估算方法：对于最小单元，热扩散时间尺度约为 $\Delta t_{\text{min}} \approx \rho c_p L_{\text{min}}^2 / k$，时间步长不要超过这个量的数倍，否则容易出现数值振荡（用隐式格式可以放宽这个限制）。实际操作中，Abaqus和ANSYS都支持自动时间步长，让软件自己根据收敛情况调整。如果发现结果出现锯齿状温度波动，第一个要检查的就是时间步长是不是太大。
3. 热分析控制方程总结
以下是CAE热分析中涉及的核心偏微分方程，所有热求解器的理论基础均来源于此：
边界条件类型
| 类型 | 数学表达式 | 工程含义 |
|---|---|---|
| 第一类（Dirichlet） | $T = T_0$ | 指定温度（如冷却水管壁温） |
| 第二类（Neumann） | $-k\partial T/\partial n = q_s$ | 指定热流密度（如芯片功耗） |
| 第三类（Robin/对流） | $-k\partial T/\partial n = h(T_s - T_\infty)$ | 表面对流换热 |
| 辐射 | $-k\partial T/\partial n = \varepsilon\sigma(T_s^4 - T_{\text{surr}}^4)$ | 表面热辐射 |
4. 主要工程应用场景
芯片结温预测、PCB热管理、散热片设计优化。关键指标：热阻 $R_\theta = \Delta T / P$（°C/W）。结温超限是电子产品失效的头号原因。
内燃机缸盖热分析、涡轮叶片冷却通道优化、排气管热疲劳分析。涡轮入口温度高达1700°C，热分析直接影响材料选择和寿命预测。
外墙保温性能、热桥分析、冬夏两季温度分布。与建筑能耗模拟（EnergyPlus等）结合，支撑绿色建筑认证（LEED、被动房标准）。
核燃料棒温度分布、堆芯热工水力分析、乏燃料贮存热管理。安全法规要求极高，V&V验证至关重要。
卫星在轨热环境（太阳辐射、地球红外、深空辐射），热控涂层选择，星上仪器温度维持。真空中只有辐射传热，无对流。
锂电池组在充放电过程中的温升、温度均匀性分析。电芯温度超过45°C会加速老化，低于0°C容量下降。热分析支撑冷却系统（液冷板）设计。
5. 耦合分析：热分析不是孤立存在的
5.1 热-结构顺序耦合（Thermal-Structural Sequential Coupling）
最常见的耦合方案：先用热求解器计算温度场，再将温度场作为载荷施加到结构模型上，计算热应力和热变形。适用条件：变形对温度场影响可以忽略。
热膨胀引起的应变（热应变）：
$$\varepsilon_{\text{th}} = \alpha (T - T_{\text{ref}})$$其中 $\alpha$ 为线膨胀系数（1/K），$T_{\text{ref}}$ 为无应力参考温度。
5.2 共轭传热（Conjugate Heat Transfer, CHT）
同时求解固体导热和流体对流，无需预先指定对流换热系数 $h$，而是通过流固界面的能量守恒自动计算。这是最准确的热分析方法，但计算代价较高。典型应用：CPU散热器优化、涡轮叶片气膜冷却分析。主要使用CFD软件（OpenFOAM、ANSYS Fluent）实现。
5.3 焦耳加热（Joule Heating / Resistive Heating）
电流流过导体产生热量，热量与电流密度和电阻率的乘积成正比：
$$Q = \mathbf{J} \cdot \mathbf{E} = \sigma_e |\mathbf{E}|^2$$其中 $\mathbf{J}$ 为电流密度（A/m²），$\mathbf{E}$ 为电场强度（V/m），$\sigma_e$ 为电导率（S/m）。需要同时求解电场方程和热方程——这是电-热耦合问题。应用：保险丝熔断分析、电阻焊接、高压导线载流量验证。
顺序耦合和完全耦合到底差在哪里？什么时候必须用完全耦合？
顺序耦合是单向的：热→结构，变形不反过来影响热。大多数工程问题——比如你算发动机缸盖热应力——用顺序耦合就够了，因为几毫米的热变形不会明显改变热边界条件。但是，如果你要模拟橡胶轮胎在路面上滚动产生热，橡胶变形改变了接触区域，接触区域的变化又改变摩擦热输入——这就必须完全耦合，热和力同时迭代求解。还有热塑性加工（热锻）、摩擦焊接等，也需要完全耦合。代价是：求解器设置更复杂，收敛更难，计算时间是顺序耦合的几倍甚至十几倍。
我在做电动车电池包热管理，冷却液是水，要用CHT方法吗？
不一定非得用CHT。对于液冷系统，如果冷却液的流动比较规则（层流或已知流态），很多工程师会这样做：用CFD单独算出冷却流道内的换热系数分布，然后把这个 $h$ 分布导入到热分析模型作为边界条件。这样可以把计算量拆分成独立的两步。如果精度要求极高，或者流道形状复杂（比如有湍流器或流道弯折），才用CHT一体化计算。ANSYS Fluent的Conjugate Heat Transfer模块就是专门用来做这件事的。OpenFOAM的 chtMultiRegionFoam 求解器也是这个功能，而且免费。
做热分析时最常见的错误是什么？
我见过的最多的错误有三类：第一，$h$ 值随意取，差了一个数量级结果差很多；第二，忘记给模型所有外表面施加边界条件，导致绝热表面温度虚高（软件默认绝热！）；第三，瞬态分析里时间步太大，导致数值振荡，结果出现负温度或者温度"跳动"。还有一个容易被忽略的：接触热阻。如果两个零件接触，接触面之间是有热阻的（不是完全导热），不考虑接触热阻会低估热点温度，在功率电子封装分析中尤其重要。
6. 主流热分析软件对比
| 软件 | 类型 | 热分析能力 | 优势 | 劣势 |
|---|---|---|---|---|
| Abaqus | 商业FEA | 稳态/瞬态导热、热-结构耦合、辐射、相变 | 非线性能力强，热-结构耦合无缝，材料库丰富 | 价格昂贵，流体分析需配合第三方 |
| ANSYS Mechanical | 商业FEA | 全面热分析 + Steady State Thermal/Transient Thermal模块 | 与ANSYS Fluent集成度高，Workbench工作流完善 | 授权费用高，大模型内存占用大 |
| OpenFOAM | 开源CFD | buoyantSimpleFoam（对流）、chtMultiRegionFoam（共轭传热） | 完全免费，并行计算能力强，CHT一流 | 纯固体传热需要专门配置，学习曲线陡 |
| ANSYS Fluent | 商业CFD | 对流传热、CHT、燃烧辐射传热 | 湍流模型丰富，辐射模型完善（P1、DO、S2S），工业认可度高 | 价格高，纯固体导热不是强项 |
| CalculiX | 开源FEA | 稳态/瞬态热分析，热-结构耦合 | 完全免费，接近Abaqus语法，适合学习 | 并行计算有限，辐射功能较弱 |
7. 热分析子栏目
选择你感兴趣的专题深入学习：
推荐学习路线
| 阶段 | 主题 | 关键概念 | 推荐工具 |
|---|---|---|---|
| 入门 | 稳态导热 | 傅里叶定律、热阻网络、温度边界条件 | CalculiX / Abaqus |
| 进阶 | 瞬态热分析与对流 | 热容、时间步控制、$h$ 值估算、Biot数 | Abaqus / ANSYS Mechanical |
| 高级 | 辐射传热 + 非线性 | 视角因子、温度相关材料属性、相变潜热 | Abaqus / ANSYS Fluent |
| 专家 | 多物理场耦合 | CHT、热-结构完全耦合、焦耳加热、数字孪生 | OpenFOAM / ANSYS Fluent + Mechanical |