第三章 热分析理论：热传导、热对流与热辐射
热分析，说白了就是搞清楚热量怎么来、怎么走、最后去了哪。我刚开始做光机仿真那会儿，总觉得热分析就是加个温度边界条件就完事了。结果有一次，一个高功率激光系统，仿真出来的镜面变形和实测差了快一倍。后来才发现，我把自然对流给忽略了。嗯，从那以后，我再也不敢小看热分析的任何一个环节了。
这一章，我们就把热传递的三种基本方式——传导、对流、辐射，以及稳态/瞬态分析、热源建模，一次性讲透。
3.1 热传导：固体内部的“热量接力”
热传导是固体内部热量传递的主要方式。它的物理本质，是分子或原子之间的动能传递。你想想看，高温区域的分子振动更剧烈，通过碰撞把能量传给相邻的低温分子。
描述热传导的核心定律是傅里叶定律：
q = -k · ∇T其中：
- q — 热流密度 (W/m²)，单位面积上的热流量
- k — 导热系数 (W/(m·K))，材料属性
- ∇T — 温度梯度 (K/m)，负号表示热量从高温流向低温
关键参数：导热系数 k
我个人的习惯是，拿到一个新材料，第一件事就是查它的导热系数随温度的变化曲线。很多金属的 k 值在高温下会明显下降，这个不搞清楚，仿真就是空中楼阁。
| 材料 | 导热系数 k (W/(m·K)) | 典型应用 |
|---|---|---|
| 纯铜 | ~400 | 散热器、热沉 |
| 铝合金 6061 | ~167 | 结构件、镜座 |
| 不锈钢 304 | ~16 | 真空腔体、支撑 |
| 熔石英 | ~1.4 | 光学镜片 |
| 空气（静止） | ~0.026 | 隔热层 |
工程经验：在光机结构中，镜片和镜座之间的接触热阻往往是热传导的瓶颈。我曾经遇到过镜片温度比镜座高十几度的情况，后来在接触面加了导热硅脂，温差立刻降到了 3°C 以内。
3.2 热对流：流体参与的“热量搬运”
热对流发生在固体表面和流动的流体之间。它分为两种：
- 自然对流：流体因温度差导致密度变化，产生浮力驱动流动。比如静止空气中的散热。
- 强制对流：通过风扇、泵等外力驱动流体流动。比如风冷散热器。
牛顿冷却定律是对流换热的计算基础：
q = h · (T_s - T_f)其中：
- h — 对流换热系数 (W/(m²·K))
- T_s — 固体表面温度
- T_f — 流体主体温度
避坑指南：对流换热系数 h 不是材料属性，它取决于流体性质、流速、表面几何形状和温度差。我曾经见过有人把 h 设成固定值 10 W/(m²·K) 就完事了，结果仿真和实测完全对不上。正确的做法是：要么查经验关联式，要么做 CFD 仿真来提取 h 值。
| 对流类型 | 典型 h 值范围 (W/(m²·K)) | 工程场景 |
|---|---|---|
| 自然对流（空气） | 2 ~ 25 | 静止环境散热 |
| 强制对流（空气） | 25 ~ 250 | 风冷散热 |
| 强制对流（水） | 100 ~ 15000 | 水冷板、液冷系统 |
| 沸腾换热 | 2500 ~ 100000 | 高热流密度散热 |
3.3 热辐射：不需要介质的“远距离传热”
热辐射是唯一不需要介质的热传递方式。在真空环境中（比如太空望远镜），辐射是唯一的热量交换途径。
斯蒂芬-玻尔兹曼定律描述了辐射换热量：
q = ε · σ · (T₁⁴ - T₂⁴)其中：
- ε — 发射率（黑度），0~1 之间
- σ — 斯蒂芬-玻尔兹曼常数，5.67×10⁻⁸ W/(m²·K⁴)
- T₁、T₂ — 两个表面的绝对温度 (K)
注意：辐射换热与温度的四次方成正比。这意味着在高温下（比如 >500°C），辐射会迅速成为主导传热方式。但在光机系统常见的 20~80°C 范围内，辐射通常只占一小部分，除非是真空环境。
我做过一个空间光学载荷的热分析，在真空条件下，镜筒内部的辐射换热必须精确建模。当时我用了角系数来计算各表面之间的辐射交换，结果发现镜筒内壁的发射率对镜片温度有显著影响。后来我们给内壁做了发黑处理，把发射率从 0.1 提高到 0.9，镜片温度均匀性提升了 40%。
3.4 稳态热分析与瞬态热分析
这两种分析的区别，说白了就是看热量是否随时间变化。
稳态热分析
系统达到热平衡后，温度场不随时间变化。控制方程为：
∇ · (k ∇T) + Q = 0其中 Q 是内热源项。稳态分析适合：
- 长时间运行的设备（如连续激光器）
- 周期性热负载的平均效果
- 热平衡后的温度分布
瞬态热分析
温度场随时间变化。控制方程增加了时间项：
ρ · c_p · ∂T/∂t = ∇ · (k ∇T) + Q其中：
- ρ — 密度
- c_p — 比热容
- ∂T/∂t — 温度变化率
瞬态分析适合：
- 脉冲式热负载（如脉冲激光）
- 启动和关机过程
- 热冲击分析
我的建议：在做瞬态分析时，时间步长的选择很关键。步长太大，会丢失温度变化的细节；步长太小，计算时间会爆炸。我一般先用稳态分析估算热时间常数 τ = ρ·c_p·L²/k，然后把时间步长设为 τ/10 左右，既能保证精度，又不至于太慢。
3.5 热源建模：把“热量从哪里来”说清楚
热源建模是热分析中最容易出错的环节。常见的几种热源类型：
| 热源类型 | 建模方式 | 典型场景 |
|---|---|---|
| 体热源 | 单位体积发热率 (W/m³) | 激光晶体吸收、电阻发热 |
| 面热源 | 单位面积热流 (W/m²) | 镜面吸收、加热膜 |
| 点热源 | 集中热流量 (W) | 小功率器件、焊点 |
| 热流密度分布 | 空间分布函数 | 高斯光束加热 |
避坑指南：我曾经犯过一个低级错误——把激光晶体的体热源直接设成均匀分布。实际上，激光晶体内部的吸收是呈指数衰减的（遵循比尔-朗伯定律）。均匀热源假设导致仿真温度比实测低了 30%。正确的做法是：用吸收系数和光强分布来计算热源的空间分布。
举个例子，高斯光束加热镜面的热流密度分布：
q(r) = q₀ · exp(-2r²/ω²)其中：
- q₀ — 中心最大热流密度
- r — 径向距离
- ω — 光束半径（1/e² 处）
这个公式在激光热效应仿真中非常常用。我建议你把它记下来，或者直接做成参数化模型，方便后续调用。
知识体系总览
下面这张图，把本章的核心知识点串在了一起。你可以把它当作一个快速索引，随时回来查阅。
本章小结
热分析不是简单的加个温度边界条件。你需要理解：
- 热量通过传导、对流、辐射三种方式传递，每种方式都有其主导场景和关键参数
- 稳态分析适合长期运行设备，瞬态分析适合脉冲或动态过程
- 热源建模必须反映真实的物理分布，均匀假设往往是误差的根源