第二章 热力学基础：热传导、热对流、热辐射的基本方程、热边界条件、材料热属性
各位同学，欢迎来到热力耦合仿真的第二站。说实话，热力学这块内容，很多搞结构的人一开始会觉得有点「隔」。我当年刚接触时也一样，觉得不就是个温度场嘛，能有多复杂？直到第一次做电子元器件的热仿真，芯片温度怎么都压不下去，我才意识到——热，其实是个很「狡猾」的东西。
这一章，咱们就把热力学里最核心的几个概念捋清楚。你不需要成为热力学专家，但得知道仿真软件背后在算什么，边界条件该怎么设，材料参数从哪来。
2.1 热传导：最基础的传热方式
热传导，说白了就是热量在固体内部「手拉手」传递。分子振动得厉害，就把能量传给旁边的分子。这个过程的数学描述，就是大名鼎鼎的傅里叶定律。
傅里叶定律（一维形式）：
q = -k * (dT/dx)其中：
- q — 热流密度 (W/m²)，单位面积上流过的热量
- k — 导热系数 (W/(m·K))，材料传导热量的能力
- dT/dx — 温度梯度 (K/m)，温度变化的剧烈程度
注意那个负号。它表示热量永远从高温区流向低温区。嗯，这个方向感很重要，我在做热流路径分析时，经常靠这个判断散热通道是否通畅。
三维情况下，傅里叶定律就变成了向量形式。不过别怕，仿真软件里已经帮你处理好了。你只需要知道：温度梯度越大，热流越猛。
热传导的控制方程
把傅里叶定律和能量守恒结合起来，就得到了热传导方程。我习惯叫它「热扩散方程」：
ρ * Cp * (∂T/∂t) = k * (∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z²) + Q这里：
- ρ — 密度 (kg/m³)
- Cp — 比热容 (J/(kg·K))，材料储存热量的能力
- ∂T/∂t — 温度随时间的变化率
- Q — 内部热源 (W/m³)，比如焦耳热、化学反应热
这个方程左边是「热量积累」，右边是「热量流入」加上「内部产热」。稳态分析时，左边等于零，温度场不再随时间变化。
我的经验：做瞬态热分析时，比热容这个参数特别关键。有一次我仿真一个电机启动过程，温度上升曲线总比实测慢。查了半天，发现供应商给的比热容数据是室温下的，而实际工作温度下比热容会变化。后来我用了分段线性数据，结果就对上了。
2.2 热对流：流体带走的热量
热对流发生在固体表面和流体之间。你想想看，风扇吹散热器，就是典型的对流换热。
描述对流的公式是牛顿冷却定律：
q = h * (T_s - T_f)其中：
- h — 对流换热系数 (W/(m²·K))
- T_s — 固体表面温度 (K)
- T_f — 流体温度 (K)
这个公式看起来简单，但难点全在那个h上。h不是材料属性，它跟流体速度、流体性质、表面形状、流动状态都有关系。
避坑指南：我曾经在仿真一个风冷系统时，直接用了手册上的「典型对流系数」——25 W/(m²·K)。结果仿真结果比实测低了十几度。后来发现，实际风速比手册假设的低，而且散热器翅片间距影响了局部对流。所以，对流系数最好通过经验公式计算或实测标定，别偷懒。
对流分为两种：
- 自然对流 — 靠温差驱动，h一般在5~25 W/(m²·K)
- 强制对流 — 靠风扇或泵驱动，h可以到50~500 W/(m²·K)
2.3 热辐射：不需要介质的传热
热辐射是唯一不需要介质的传热方式。太阳的热量穿过真空到达地球，靠的就是辐射。
描述辐射的核心公式是斯特藩-玻尔兹曼定律：
q = ε * σ * (T₁⁴ - T₂⁴)其中：
- ε — 发射率 (0~1)，黑体的发射率为1
- σ — 斯特藩-玻尔兹曼常数，5.67×10⁻⁸ W/(m²·K⁴)
- T₁、T₂ — 两个表面的绝对温度 (K)
注意，温度是四次方关系。这意味着高温下辐射会急剧增强。我做高温炉仿真时，辐射往往是主要的传热方式，对流反而可以忽略。
关键点：辐射计算是非线性的，而且涉及表面之间的「视角因子」。在仿真中，辐射设置通常比传导和对流复杂得多。我的建议是：如果温差不大（比如几十度以内），辐射可以忽略；如果温差上百度，辐射必须考虑。
2.4 热边界条件：给问题定个调
有了控制方程，还得有边界条件才能求解。热分析中常见的边界条件有四种：
| 类型 | 描述 | 数学形式 | 典型应用 |
|---|---|---|---|
| 第一类 | 给定温度 | T = T₀ | 恒温壁、冰水混合物 |
| 第二类 | 给定热流 | q = q₀ | 加热膜、激光加热 |
| 第三类 | 对流换热 | q = h(T_s - T_f) | 风冷、水冷 |
| 第四类 | 辐射换热 | q = εσ(T⁴ - T_env⁴) | 高温炉、太空环境 |
实际工程中，往往是多种边界条件同时存在。比如一个电子设备，外壳既有对流散热，也有辐射散热，内部还有热传导。
我的习惯：设置边界条件时，先问自己三个问题：
- 这个表面是绝热的还是散热的？
- 散热的主要方式是什么？对流还是辐射？
- 边界条件的数据来源可靠吗？是实测还是估算？
这三个问题想清楚，边界条件就不会设错。
2.5 材料热属性：仿真准确度的基石
仿真结果准不准，很大程度上取决于材料参数准不准。我见过太多人随便从网上扒个导热系数就开始算，结果可想而知。
导热系数 (k)
导热系数衡量材料传导热量的能力。金属的k很高（铜约400 W/(m·K)），气体的k很低（空气约0.026 W/(m·K)）。
注意：导热系数通常是温度的函数。比如不锈钢，从室温到500°C，k可能变化20%以上。做高温分析时，一定要用温度相关的导热系数。
比热容 (Cp)
比热容衡量材料储存热量的能力。水的比热容很大（4180 J/(kg·K)），所以水冷效果很好。金属的比热容相对较小（铜约385 J/(kg·K)）。
比热容也随温度变化，尤其在相变温度附近变化剧烈。做焊接或热处理仿真时，这个参数必须准确。
热膨胀系数 (α)
热膨胀系数描述材料受热后尺寸的变化。这是热力耦合分析中连接温度场和应力场的桥梁。
ε_thermal = α * (T - T_ref)其中：
- ε_thermal — 热应变
- α — 线膨胀系数 (1/K)
- T_ref — 参考温度，通常取无应力状态下的温度
避坑指南：我曾经做一个电子封装的热应力分析，结果应力值大得离谱。查了半天，发现是参考温度设错了。封装工艺中，材料在固化温度下是零应力状态，而我用了室温作为参考温度。这个错误直接导致应力结果翻倍。所以，参考温度一定要跟实际工艺对应。
2.6 知识体系总览
下面这张图，是我自己总结的热力学基础框架。你可以把它当作这一章的「地图」：
这张图把这一章的核心内容串起来了。从三大传热方式出发，到控制方程，再到边界条件和材料属性。你每次做热仿真前，都可以对照这张图检查一下：我的模型里，传热方式考虑全了吗？边界条件设对了吗？材料参数准确吗？
2.7 本章小结
这一章的内容，说白了就是热仿真的「语法」。掌握了这些，你才能读懂仿真软件在干什么，才能判断结果合不合理。
最后说一句：热分析没有捷径，但一定有方法。我做了十几年仿真，最大的体会就是——基础越扎实，遇到问题越不慌。下一章，咱们就开始把这些理论用到实际案例中。
公众号:蓝海资料掘金营，微信deep3321