4、热分析基础：热传导、热对流、热辐射三种传热方式，热边界条件，稳态与瞬态热分析概念
光机热分析的核心在于理解热量如何在光学系统内部及系统与外界环境之间传递。无论是星载相机的轨道外热流变化，还是高功率激光光学元件的热沉积，其本质都遵循热力学三大定律。本章作为后续光机热耦合分析的理论基石，将系统梳理三种基本传热方式的物理本质与数学描述，明确热边界条件的工程分类，并厘清稳态与瞬态分析在工程应用中的适用场景与求解逻辑。只有透彻掌握这些基础概念，才能在实际工程中正确选择热分析工具、设定边界条件并合理解读温度场与热变形结果。
4.1 三种基本传热方式的物理本质与工程描述
4.1.1 热传导：固体内部热量传递的主导机制
热传导是光机结构中最重要的传热方式。当物体内部存在温度梯度时，热量会从高温区域向低温区域传递，其微观机制包括晶格振动（声子）和自由电子迁移。在光学系统设计中，镜体材料（如微晶玻璃、碳化硅、ULE）的导热系数直接决定了温度梯度的大小和热变形程度。傅里叶定律给出了热流密度与温度梯度的定量关系：q = -k·∇T，其中 k 为导热系数（W/m·K）。对于各向异性材料（如碳纤维复合材料），导热系数为张量形式，需分别定义三个方向的分量。
在工程计算中，热传导路径的截面积和长度是影响热阻的关键因素。例如，一个直径 100mm、厚度 20mm 的反射镜，若采用微晶玻璃（k≈1.6 W/m·K），其轴向热阻约为 0.1 K/W；若改用碳化硅（k≈180 W/m·K），热阻降至约 0.001 K/W，两者相差两个数量级。这解释了为什么高功率激光反射镜普遍采用碳化硅或金刚石等高导热材料。
4.1.2 热对流：流体与固体表面的热量交换
热对流发生在固体表面与流动流体之间，其基本规律由牛顿冷却定律描述：q = h·(T_w - T_∞)，其中 h 为对流换热系数（W/m²·K），T_w 为壁面温度，T_∞ 为流体主流温度。对流换热系数 h 不是材料属性，而是取决于流动状态（层流/湍流）、流体物性、流速和几何形状的综合性参数。
在光机系统中，对流换热主要出现在以下场景：地面测试时的自然对流（h≈5-25 W/m²·K）、风冷强迫对流（h≈20-200 W/m²·K）、以及液冷系统（h≈500-10000 W/m²·K）。对于星载设备，由于空间接近真空，对流换热可以忽略，热辐射成为唯一的外部换热途径。这一差异导致地面热平衡测试与在轨热状态存在显著区别，必须通过热控涂层和辐射散热面进行补偿。
4.1.3 热辐射：真空环境下的主导传热方式
热辐射是唯一不需要介质即可传播的传热方式，其能量传递遵循 Stefan-Boltzmann 定律：q = ε·σ·(T₁⁴ - T₂⁴)，其中 ε 为表面发射率（0-1），σ 为 Stefan-Boltzmann 常数（5.67×10⁻⁸ W/m²·K⁴）。辐射传热与温度的四次方成正比，因此在高温或大温差场景下占据主导地位。
在光机热分析中，辐射换热的建模最为复杂，涉及表面间的角系数计算、多次反射吸收以及光谱特性。工程上通常采用以下简化策略：
- 灰体假设：假设表面发射率与波长无关，适用于大多数金属和涂层材料。
- 角系数计算：对于简单几何（平行平板、同心圆柱），可查表获得；复杂结构需使用蒙特卡洛法或辐射度算法。
- 辐射-传导耦合：在真空环境中，辐射与传导往往同时存在，需建立辐射网络与传导网络的耦合方程。
| 传热方式 | 控制方程 | 关键参数 | 光机工程典型场景 | 分析复杂度 |
|---|---|---|---|---|
| 热传导 | q = -k·∇T | 导热系数 k (W/m·K) | 镜体内部热扩散、结构件传热 | 低（线性） |
| 热对流 | q = h·(T_w - T_∞) | 对流系数 h (W/m²·K) | 地面风冷测试、液冷系统 | 中（需CFD） |
| 热辐射 | q = ε·σ·(T₁⁴-T₂⁴) | 发射率 ε、角系数 | 星载散热、低温光学系统 | 高（非线性） |
4.2 热边界条件的工程分类与施加方法
热边界条件描述了热分析模型与外部环境之间的相互作用，是求解热传导方程的必要条件。在光机热分析中，边界条件的准确性直接影响温度场计算结果的可靠性。根据数学形式的不同，热边界条件分为三类：
4.2.1 第一类边界条件（Dirichlet条件）
给定边界上的温度分布，即 T_w = f(x,y,z,t)。在光机分析中，这类边界常用于模拟与恒温平台接触的安装面，或已知温度的热沉。例如，光学平台底部与温控板接触，可设定为固定温度边界。需要注意的是，实际工程中完全恒温的边界很少存在，施加第一类边界条件时应确认接触面热阻足够小，否则会人为引入过大的热流。
4.2.2 第二类边界条件（Neumann条件）
给定边界上的热流密度，即 q_w = -k·∂T/∂n = g(x,y,z,t)。这类边界适用于已知热载荷的情况，如激光辐照在镜面产生的热流密度、加热片施加的恒定功率等。在光机热分析中，第二类边界条件常用于模拟：
- 高功率激光束在反射镜表面的热沉积（通常为高斯分布热流）
- 加热片或热电制冷器的热流输入
- 绝热边界（q_w = 0），用于对称面或隔热良好的表面
4.2.3 第三类边界条件（Robin条件）
给定边界上的对流换热条件，即 -k·∂T/∂n = h·(T_w - T_∞)。这是工程中最常见的边界类型，适用于所有与流体接触的表面。在光机系统中，第三类边界条件用于模拟：
- 镜面与周围空气的自然对流
- 风冷或液冷系统的强迫对流
- 恒温罩内壁与光学元件之间的辐射-对流复合换热
4.3 稳态与瞬态热分析的概念与工程选择
4.3.1 稳态热分析：平衡状态的温度场求解
稳态热分析求解的是系统在恒定热载荷和边界条件下达到热平衡后的温度分布，其控制方程为拉普拉斯方程（无内热源）或泊松方程（有内热源）：∇·(k∇T) + Q = 0。稳态分析不考虑时间项，计算量较小，适用于以下工程场景：
- 连续工作状态下的光学系统热平衡分析（如长时间工作的红外热像仪）
- 恒定环境温度下的地面测试工况
- 热控设计初步评估，快速比较不同方案的优劣
稳态分析的关键在于判断系统是否真正达到热平衡。工程上通常以温度变化率小于 0.1°C/h 作为稳态判据。对于大热惯性系统（如大质量镜体），达到稳态可能需要数小时甚至数天。
4.3.2 瞬态热分析：时间演化的温度场求解
瞬态热分析求解温度场随时间的变化，控制方程为：ρc·∂T/∂t = ∇·(k∇T) + Q，其中 ρ 为密度，c 为比热容。瞬态分析需要给定初始温度场和随时间变化的边界条件，计算量显著增大，但能揭示以下关键信息：
- 温度梯度随时间的变化速率，评估热冲击效应
- 达到热平衡所需的时间（热时间常数）
- 周期性热载荷（如卫星进出地影）引起的温度波动
- 瞬态热变形对光学性能的动态影响
在光机热分析中，瞬态分析尤为重要。例如，卫星在轨道上经历周期性太阳照射和地影遮挡，光学系统温度呈现周期性波动，导致镜面面形随之变化。这种动态热变形无法通过稳态分析获得，必须进行瞬态热-结构耦合分析。
4.3.3 稳态与瞬态的选择判据
| 判断维度 | 稳态分析 | 瞬态分析 |
|---|---|---|
| 时间依赖性 | ∂T/∂t = 0 | ∂T/∂t ≠ 0 |
| 热载荷特征 | 恒定或缓慢变化 | 快速变化或周期性 |
| 计算资源 | 低（单次求解） | 高（多时间步迭代） |
| 输出结果 | 最终平衡温度场 | 温度场时间序列 |
| 典型应用 | 热控方案初选、恒定工况评估 | 热冲击分析、轨道热环境模拟 |
| 工程判据 | 温度变化率 < 0.1°C/h | 热时间常数 τ = L²/α 与载荷周期比较 |
4.4 案例拆解：星载相机主镜的瞬态热分析
为综合理解上述概念，以下通过一个典型工程案例展示热分析方法的实际应用。
背景与约束：某低轨遥感卫星搭载一台口径 600mm 的卡塞格林望远镜，主镜材料为微晶玻璃（k=1.6 W/m·K，ρ=2530 kg/m³，c=821 J/kg·K），镜体通过 3 个钛合金柔性支撑与碳纤维镜筒连接。卫星轨道高度 500km，轨道周期 95 分钟，其中日照期 65 分钟，地影期 30 分钟。太阳辐射强度 1350 W/m²，地球反照率 0.3，地球红外辐射 240 W/m²。
初始现象与目标：在轨测试发现，当卫星从日照期进入地影期时，主镜面形 RMS 值从 12nm 增大至 28nm，超出光学设计容差（20nm）。需要分析温度场变化规律，找出面形退化的热致因素。
关键输入数据：
- 主镜正面镀高反射膜（ε=0.05），背面及侧面镀黑漆（ε=0.88）
- 镜筒内壁涂黑（ε=0.92），温度维持在 20±2°C
- 主镜初始温度 22°C，均匀分布
- 热时间常数估算：τ = ρcL²/k = 2530×821×(0.03)²/1.6 ≈ 1170s ≈ 19.5min
分析判断过程：由于轨道周期（95min）与热时间常数（19.5min）处于同一量级，必须采用瞬态热分析。建立有限元模型，施加以下边界条件：
- 主镜正面：辐射边界（吸收太阳辐射+地球反照+地球红外，同时向深空辐射）
- 主镜背面及侧面：辐射边界（与镜筒内壁辐射换热）
- 柔性支撑处：传导边界（与镜筒热传导）
- 镜筒外表面：辐射边界（与空间环境换热）
采用瞬态求解器，时间步长 60s，总时长 3 个轨道周期（285min），确保温度场进入周期性稳态。
实际操作步骤：
1. 建立几何模型：主镜、柔性支撑、镜筒（简化）
2. 定义材料属性：k, ρ, c, ε（分表面）
3. 设置辐射组：定义所有辐射面及角系数
4. 施加热载荷：太阳辐射（随时间变化）、地球反照、地球红外
5. 设置求解控制：时间步长 60s，迭代收敛准则 1e-4
6. 后处理：提取主镜前表面温度分布及面形变化结果对比：瞬态分析结果显示，进入地影后 10 分钟，主镜前表面中心温度下降 3.2°C，边缘温度下降 1.8°C，形成径向温度梯度约 0.14°C/cm。该温度梯度导致镜面产生热变形，面形 RMS 从 12nm 增大至 25nm，与实测值 28nm 基本吻合（差异主要来自支撑结构简化）。
风险与代价：若采用稳态分析（取日照期平均热流），计算得到的温度梯度仅为 0.05°C/cm，面形 RMS 为 15nm，严重低估了实际热变形。这说明瞬态效应在周期性热载荷下不可忽略。
复盘经验：
- 对于轨道周期与热时间常数接近的系统，必须进行瞬态分析
- 辐射边界条件的角系数计算精度对结果影响显著，需仔细建模
- 主镜正面高反射膜（低发射率）有效降低了太阳辐射吸收，但背面黑漆导致散热过快，形成温度梯度
- 优化方向：在镜体侧面增加热控涂层梯度，或增加热惯性补偿环