4、热学基础:热传导、对流、辐射基本定律、热膨胀系数、热弹性效应
光机热一体化设计的核心挑战在于:温度场变化会引起结构变形,结构变形又改变光学系统的面形精度与装调状态,最终导致成像质量或光束指向性能退化。因此,热学基础不是孤立的热力学复习,而是建立“热载荷→温度场→结构响应→光学性能”分析链路的起点。本章聚焦热传导、对流、辐射三大传热机制的基本定律,以及热膨胀系数与热弹性效应的工程含义,为后续光机结构的热-力耦合分析奠定物理基础。
4.1 热传导:傅里叶定律与热阻网络
核心原理
热传导是固体内部或相互接触的固体之间依靠分子振动和自由电子运动传递热量的过程。傅里叶定律给出了热流密度与温度梯度之间的线性关系:
q = -k · ∇T
其中 q 为热流密度(W/m²),k 为导热系数(W/(m·K)),∇T 为温度梯度(K/m)。负号表示热量从高温向低温传递。对于光机结构中的一维稳态导热(如镜座到镜体的热传导路径),可简化为:
Q = k·A·(T₁-T₂)/L
式中 A 为垂直于热流方向的截面积,L 为导热路径长度。这一形式与欧姆定律高度相似,因此工程上常将导热路径等效为热阻:
R_cond = L/(k·A)
热阻网络法将复杂结构分解为串联、并联的多个热阻单元,是光机系统热分析中最常用的手算与初步仿真验证手段。
判断方法或操作流程
在光机结构热传导分析中,建议按以下流程操作:
- 识别主要导热路径:从热源(如探测器、电机、太阳辐照)到光学元件(镜片、镜座)之间的所有固体接触链。
- 计算各段热阻:对每一段均质材料,按 R=L/(k·A) 计算;对接触界面,需查接触热阻经验值(通常为 10⁻⁴~10⁻³ m²·K/W 量级)。
- 建立热阻网络:串联路径热阻相加,并联路径按倒数相加。
- 求解节点温度:结合边界条件(热流输入、恒温边界、对流散热),用节点法或热路方程求解。
- 校核关键温度:重点校核光学镜片中心与边缘温差、镜座与镜体接触面温度差。
案例拆解:反射镜镜体与镜座间的热传导设计
某空间相机主镜直径为 300mm,采用微晶玻璃(k=1.5 W/(m·K)),镜座为钛合金(k=7.5 W/(m·K)),两者通过 3 个柔性支撑连接。在轨时镜体吸收太阳辐照约 5W,需将热量传导至镜座(恒温 20°C)。初始设计采用平面接触,接触面积 1200mm²,接触热阻 5×10⁻⁴ m²·K/W。计算镜体温度:
接触热阻 R_c = 5×10⁻⁴ / (1200×10⁻⁶) = 0.417 K/W
镜体自身导热热阻 R_m ≈ L/(k·A) = 0.02/(1.5×0.07) ≈ 0.19 K/W
总热阻 R_total = 0.417 + 0.19 = 0.607 K/W
镜体温升 ΔT = Q × R_total = 5 × 0.607 ≈ 3.04 K
镜体温度约 23°C,尚可接受。但若改用 3 个柔性片支撑(每个接触面积仅 80mm²),接触热阻将增大至 2.08 K/W,镜体温升达 11.4°C,将引起显著热变形。因此设计时需在柔性支撑的力学解耦能力与热传导效率之间权衡,必要时在接触面填充导热脂(导热系数 2~5 W/(m·K))以降低界面热阻。
4.2 热对流:牛顿冷却定律与自然/强迫对流
核心原理
热对流是流体(空气、冷却液等)与固体表面之间因温度差引起的热量交换。牛顿冷却定律描述了这一过程:
Q = h·A·(T_s - T_∞)
其中 h 为对流换热系数(W/(m²·K)),A 为换热面积,T_s 为固体表面温度,T_∞ 为流体远场温度。对流换热系数 h 是核心参数,其典型取值范围如下表:
| 对流类型 | h 范围 (W/(m²·K)) | 典型应用场景 |
|---|---|---|
| 自然对流(空气) | 2~10 | 地面光学平台、无风环境 |
| 强迫对流(风扇) | 10~100 | 机载光电吊舱、主动风冷 |
| 强迫对流(液冷) | 100~1000 | 高功率激光器、探测器制冷 |
| 相变换热(沸腾/冷凝) | 1000~10000 | 热管、环路热管 |
判断方法或操作流程
在光机系统对流散热设计中,需判断以下关键问题:
- 自然对流还是强迫对流? 若系统无风扇且体积受限,只能依赖自然对流,h 取 3~8 W/(m²·K);若允许加装风扇或泵,可显著提高散热能力。
- 对流面积是否充足? 计算所需散热面积 A = Q/(h·ΔT),若面积不足需增加散热翅片。
- 气流路径是否顺畅? 光学系统常因密封要求限制通风,需评估密封腔内自然对流的有效性。
- 重力方向影响:自然对流时,热表面朝上换热系数约为朝下的 1.5~2 倍,设计时需考虑安装方位。
实战场景:机载光电吊舱的散热设计
某机载光电吊舱内部有红外探测器(热耗 8W)和伺服电机(热耗 15W),吊舱外壳为铝合金密封结构,外部飞行风速 50m/s。初始设计仅依靠外壳自然对流散热,外壳表面积 0.12m²,内部空气温度可达 65°C,超出探测器工作上限 55°C。分析发现:
- 自然对流散热能力:Q = 5 × 0.12 × (65-20) ≈ 27W,不足以带走 23W 总热耗(且内部温升已超标)。
- 改进方案:在吊舱外壳增加纵向翅片(面积增至 0.35m²),利用飞行时的强迫对流(h≈40 W/(m²·K)),散热能力提升至 Q = 40 × 0.35 × 30 ≈ 420W,远大于热耗。
- 实际实施:翅片方向与气流方向一致,翅片高度 15mm、间距 8mm,内部热源通过导热板传导至外壳。
改进后内部温度稳定在 42°C,满足探测器要求。此案例说明,在高速飞行环境下,强迫对流是光机系统最有效的散热途径,但需注意翅片增加的风阻和重量代价。
4.3 热辐射:斯蒂芬-玻尔兹曼定律与辐射换热
核心原理
热辐射是物体因温度而向外发射电磁波的能量传递方式,不需要介质,在真空环境中(如太空)是唯一传热途径。斯蒂芬-玻尔兹曼定律描述黑体的辐射力:
E_b = σ·T⁴
其中 σ = 5.67×10⁻⁸ W/(m²·K⁴) 为斯蒂芬-玻尔兹曼常数。实际物体表面辐射力为 E = ε·σ·T⁴,ε 为表面发射率(0~1)。两个表面之间的净辐射换热量为:
Q_rad = ε_eff·σ·A₁·(T₁⁴ - T₂⁴)
ε_eff 为系统有效发射率,与两表面的发射率、角系数有关。辐射换热与温度的四次方成正比,因此在高温或低温(太空)场景中辐射占主导;在常温地面环境中,辐射与对流往往量级相当。
判断方法或操作流程
光机系统辐射换热分析的关键步骤:
- 判断辐射是否重要:若系统处于真空环境、或温差大于 50°C、或表面发射率大于 0.5,辐射必须考虑。
- 确定表面辐射属性:查阅材料发射率(抛光铝约 0.05,黑漆约 0.9,阳极氧化铝约 0.8)。
- 计算角系数:对简单几何(平行平板、同心圆柱)查表;复杂结构用辐射仿真软件。
- 评估太阳辐照:在轨光学载荷需计算太阳辐射吸收(吸收率 α_s 与发射率 ε 的比值决定平衡温度)。
案例拆解:空间望远镜遮光罩的辐射热设计
某空间望远镜遮光罩为铝制圆筒(外径 400mm,长 600mm),外表面喷涂黑漆(ε=0.9,α_s=0.9),内表面抛光(ε=0.05)。在轨时太阳辐照 1367 W/m²,遮光罩轴线与太阳光夹角 30°。计算遮光罩平衡温度:
太阳吸收功率:P_abs = α_s · S · A_proj = 0.9 × 1367 × (0.4×0.6×sin30°) ≈ 147.6 W
辐射散热面积:A_rad ≈ 2π×0.2×0.6 + 2π×0.2² ≈ 0.754 + 0.251 = 1.005 m²
平衡方程:P_abs = ε·σ·A_rad·T⁴
T⁴ = 147.6 / (0.9 × 5.67×10⁻⁸ × 1.005) ≈ 2.88×10⁹
T ≈ 231 K = -42°C
若将外表面改为低吸收率白漆(α_s=0.2,ε=0.85),则 P_abs = 32.8W,T ≈ 160K = -113°C。可见表面辐射属性对空间光机系统的温度水平有决定性影响。设计时需根据光学元件的工作温度要求,选择合适的热控涂层。
4.4 热膨胀系数与热弹性效应
核心原理
热膨胀系数(CTE,α)描述材料在温度变化时的尺寸变化率:
ΔL/L₀ = α·ΔT
对于各向同性材料,体膨胀系数约为线膨胀系数的 3 倍。光机结构中,不同材料的热膨胀失配是热应力的主要来源。热弹性效应则描述了温度变化引起的应力-应变耦合关系,其本构方程为:
ε_total = ε_mechanical + ε_thermal = (1/E)·σ + α·ΔT
当结构受到约束时,温度变化会产生热应力:σ = -E·α·ΔT(完全约束时)。光机设计中,热弹性效应直接导致镜面曲率半径变化、镜间距改变、支撑结构变形。
| 材料 | CTE (×10⁻⁶/K) | 导热系数 (W/(m·K)) | 弹性模量 (GPa) | 光机应用 |
|---|---|---|---|---|
| 微晶玻璃 Zerodur | 0.02~0.05 | 1.5 | 90 | 高稳定反射镜基底 |
| ULE 钛硅玻璃 | 0.03 | 1.3 | 67 | 空间望远镜镜片 |
| 铝合金 6061 | 23.6 | 167 | 69 | 镜座、结构件 |
| 钛合金 Ti-6Al-4V | 8.8 | 7.5 | 114 | 柔性支撑、低热导结构 |
| 殷钢 Invar 36 | 1.2 | 10 | 141 | 精密装调基准 |
| 碳纤维复合材料 | -1~2(可设计) | 5~30 | 100~200 | 轻量化结构、遮光罩 |
判断方法或操作流程
热弹性效应分析在光机设计中的标准流程:
- 确定温度载荷:从热分析获得结构温度场(均匀或梯度分布)。
- 识别约束条件:明确结构的边界约束(固定点、柔性支撑、自由膨胀)。
- 计算热变形:对自由膨胀结构,直接计算 ΔL = α·L·ΔT;对约束结构,需用有限元求解热应力。
- 评估光学性能影响:将镜面变形(曲率半径变化、面形 PV/RMS)代入光学仿真,判断是否满足指标。
- 优化设计:若变形超差,采用低 CTE 材料、无热化设计(匹配镜筒与镜片 CTE)、或柔性支撑释放热应力。
实战场景:无热化镜头设计中的热弹性匹配
某红外镜头工作温度范围 -20°C~+50°C,镜筒为铝合金(α=23.6×10⁻⁶/K),镜片为锗(α=6.1×10⁻⁶/K),镜片焦距 f=100mm。若不采取补偿措施,温度升高 70°C 时:
镜筒伸长:ΔL = 23.6×10⁻⁶ × 100 × 70 = 0.165mm
镜片焦距变化:Δf ≈ f·α_ge·ΔT = 100 × 6.1×10⁻⁶ × 70 = 0.043mm
像面位移:Δ = ΔL - Δf = 0.165 - 0.043 = 0.122mm
0.122mm 的像面位移远超红外焦平面允许的离焦量(通常 ±0.02mm)。解决方案:
- 机械被动无热化:在镜筒中插入殷钢补偿环(α=1.2×10⁻⁶/K),利用不同材料的膨胀差抵消像面位移。
- 光学被动无热化:选用负热差光学材料(如 ZnSe,α=7.1×10⁻⁶/K,dn/dT 为负)与锗组合,使焦距变化与镜筒伸长匹配。
- 主动补偿:通过温度传感器驱动电机微调镜片位置。
最终采用机械被动方案,在镜筒中设计殷钢补偿筒,计算其长度 L_invar 满足:
α_al·L_al - α_invar·L_invar = α_ge·f
23.6×10⁻⁶ × 100 - 1.2×10⁻⁶ × L_invar = 6.1×10⁻⁶ × 100
L_invar ≈ 1458mm(不现实)
单段殷钢补偿长度过大,实际采用多级放大机构(杠杆结构)将殷钢的微小膨胀放大,实现紧凑设计。此案例说明热弹性匹配需要结合结构创新,而非简单材料替换。
4.5 三种传热机制的工程权衡与综合应用
核心原理
实际光机系统中,三种传热机制往往同时存在。工程判断的关键是识别主导机制并合理简化。以下给出典型场景的传热主导机制判断表:
| 应用场景 | 主导传热机制 | 次要机制 | 工程简化建议 |
|---|---|---|---|
| 太空轨道(真空) | 辐射 | 传导(结构内部) | 忽略对流,辐射+传导耦合分析 |
| 地面实验室(无风) | 自然对流+辐射 | 传导 | 两者量级相当,不可忽略任一 |
| 机载/车载(高速气流) | 强迫对流 | 传导、辐射 | 以对流为主,辐射可简化 |
| 高功率激光器内部 | 传导(液冷) | 辐射 | 液冷通道设计决定热性能 |
| 密封光学模组 | 传导(外壳) | 内部辐射 | 内部对流极弱,以传导+辐射为主 |
综合设计流程
光机热一体化设计中的热学分析应遵循以下闭环流程:
- 热载荷定义:明确内部热耗(探测器、电机、电路)和外部热环境(太阳辐照、气动加热、冷黑空间)。
- 传热路径规划:设计低热阻导热路径(高热导率材料、大面积接触、热管)将热量引导至散热面。
- 温度场计算:采用热网络法或有限元法求解温度分布,重点关注光学元件的温度梯度和温度水平。
- 热变形评估:将温度场作为载荷施加到结构模型,计算镜面变形和位置偏移。
- 光学性能验证:将变形结果导入光学软件,评估波前误差、MTF 或指向精度。
- 迭代优化:若性能不达标,调整材料选择、散热方案或支撑结构,重复上述步骤。
案例复盘:星敏感器热设计全过程
某星敏感器(光学头部+电子学箱)在轨工作,要求光学头部温度稳定在 20±2°C。初始设计将电子学箱(热耗 5W)与光学头部通过铝合金支架连接,支架截面积 200mm²,长度 80mm。分析发现:
支架导热热阻:R = 0.08/(167×200×10⁻⁶) = 2.4 K/W
电子学箱到光学头部的传热:Q = 5W,ΔT = 12K
光学头部温度 = 电子学箱温度 - 12K
若电子学箱温度 35°C,则光学头部仅 23°C,尚可接受;但若电子学箱温度升至 40°C(高温工况),光学头部达 28°C,超出指标。改进措施:
- 在支架与光学头部之间增加热管(热阻 0.1 K/W),将 ΔT 降至 0.5K。
- 在光学头部外表面包覆多层隔热材料(MLI),减少对冷黑空间的辐射散热,保持温度稳定。
- 增加主动加热器(2W),在低温工况下补偿辐射损失。
最终光学头部温度控制在 19.5~21.5°C,满足指标。此案例展示了热传导路径优化、辐射隔热和主动控温三种手段的综合运用。