4. 热分析基础与光学系统热载荷

光机热结构光学协同仿真的第一步,是把光学系统在真实环境中所受的热作用,准确翻译成热分析模型能够接受的边界条件与载荷。热分析求解的本质是能量守恒:进入系统的热量减去离开系统的热量,等于系统内能的变化。对光学载荷而言,热量通过传导、对流、辐射三种方式进入或离开结构,而太阳辐照、地球反照、内部热源和制冷机冷指则是最典型的四类热载荷来源。理解它们的物理机制与建模方式,是后续热变形与光学性能耦合分析的基础。

三种传热边界条件的物理含义与建模

传导发生在固体内部或相互接触的固体之间,遵循傅里叶定律,热流密度与温度梯度成正比,比例系数为材料导热率。在有限元热分析中,传导由材料导热率、几何尺寸和网格自动处理,用户需要关注的是接触热阻。两个名义上贴合的零件,实际只有微观凸起处接触,真实接触面积远小于名义面积,接触热阻往往不可忽略。工程上常用接触导热系数(TCC)描述,单位是 W/(m²·K),装配良好的金属-金属界面可取 1000~5000 W/(m²·K),而通过导热硅脂或铟箔的界面可显著提高。若忽略接触热阻,仿真会高估散热能力,导致温度预测偏低。

对流是固体表面与流体之间的换热,遵循牛顿冷却定律,热流密度等于对流换热系数乘以表面与流体的温差。对流换热系数 h 的取值高度依赖流动状态:自然对流在空气中通常为 2~10 W/(m²·K),强迫风冷可达 20~100 W/(m²·K),液体冷却则更高。在轨光学载荷处于真空环境,对流几乎不存在,因此空间热分析中通常不设对流边界,除非载荷内部存在主动流体回路。地面试验阶段则必须考虑对流,且地面与在轨的换热路径差异巨大,这是热试验与在轨预测偏差的主要来源之一。

辐射是真空环境中最重要的换热方式,遵循斯特藩-玻尔兹曼定律,辐射热流与绝对温度的四次方成正比。辐射换热需要定义表面发射率 ε 和吸收率 α,两者对波长和表面状态敏感。抛光金属发射率可低至 0.03,而黑漆或阳极氧化表面可达 0.85 以上。在热分析软件中,辐射通常通过角系数和辐射换热矩阵处理,需要定义辐射面、封闭腔体以及空间环境温度。对于在轨光学载荷,辐射是散热的主要途径,表面光学属性直接决定温度水平。

光学载荷热载荷与传热路径示意 光学载荷本体 镜筒 / 镜座 / 探测器 T₀ 初始温度 20℃ 太阳辐照 1367 W/m² · α/ε 关键 辐射加热 地球反照 反照率 0.3 · 轨道角 短波辐射 地球红外 237 W/m² · 长波 辐射加热 内部热源 探测器 5W · 电子学 8W 焦耳热 制冷机冷指 80K · 吸热 2W 冷量输入 向空间辐射 ε·σ·T⁴ · 4K 背景 散热 传导:接触热阻 TCC 1000~5000 W/(m²·K) | 对流:地面 h=2~100 W/(m²·K),在轨真空可忽略 稳态:能量平衡,温度不随时间变化 | 瞬态:ρc∂T/∂t = ∇·(k∇T) + Q,关注时间常数与温度变化率

四类典型热载荷的工程处理

太阳辐照是空间光学载荷最主要的加热源。太阳常数约 1367 W/m²,但实际到达表面的热流取决于太阳入射角与表面吸收率。建模时通常将太阳辐照作为表面热流载荷施加,方向由轨道和姿态确定。关键参数是吸收率 α 与发射率 ε 的比值 α/ε:比值越大,表面在阳光下升温越高。例如白漆 α/ε 约 0.2/0.85,平衡温度较低;而抛光铝 α/ε 约 0.2/0.05,虽然吸收少但发射更少,反而可能升温。因此选择表面处理时,必须同时考虑 α 和 ε,而非只看吸收率。

地球反照是地球反射的太阳短波辐射,热流大小与轨道高度、太阳-地球-载荷夹角以及地球反照率有关。反照率通常取 0.3,但随云层和地表变化,工程上常按 0.3 设计、0.35 校核。地球红外辐射是地球自身发出的长波辐射,约 237 W/m²,对低轨载荷影响显著。这两项在热分析中通常作为轨道热流边界条件,通过轨道参数和姿态角计算得到随时间变化的热流曲线。

内部热源包括探测器焦耳热、电子学功耗、电机发热等。建模时将其作为体热源或面热源施加在对应部件上。例如某红外探测器功耗 5 W,电子学组件 8 W,需分别赋予对应体或面。内部热源通常相对稳定,但在工作模式切换时可能阶跃变化,瞬态分析中需定义功率随时间的变化曲线。

制冷机冷指是低温光学系统的关键边界。冷指将探测器或光学元件冷却到 80 K 甚至更低,建模时通常作为定温边界或负热流边界。若作为定温边界,需注意冷指与载荷之间的接触热阻和热链接设计;若作为热流边界,则需输入制冷量随温度的变化曲线。冷指温度波动会直接传递到光学元件,因此瞬态分析中常需考虑制冷机启停或负载变化引起的温度波动。

稳态与瞬态热分析的设置要点

稳态热分析假设温度不随时间变化,能量平衡方程为 ∇·(k∇T) + Q = 0。它用于确定载荷在长期工作或极端工况下的温度分布,是热变形分析的输入。设置时需注意:边界条件必须完整,若存在辐射,需定义辐射封闭腔和空间背景温度;若存在接触,需定义接触热阻;若存在内部热源,需确认功率单位与分布方式。稳态求解对初始条件不敏感,但对边界条件敏感,因此边界条件的合理性直接决定结果可信度。

瞬态热分析考虑温度随时间的变化,控制方程为 ρc∂T/∂t = ∇·(k∇T) + Q。它用于分析开机升温、轨道周期温度波动、制冷机启停等过程。设置时需定义初始温度场、时间步长和总时长。时间步长应小于系统最小时间常数的 1/5~1/10,否则可能丢失温度变化细节。例如某光学镜筒时间常数约 600 s,时间步长可取 60~120 s。总时长应覆盖至少 3~5 个轨道周期,以观察温度是否进入周期性稳态。结果判断时,关注关键光学元件的温度变化率、最大温差和温度梯度,这些量直接决定热变形和光学性能漂移。

工程提示:在轨真空环境中,对流可忽略,辐射是主要散热途径。地面试验时若存在空气对流,需在仿真中增加对流边界以匹配试验条件,否则地面与在轨预测会出现系统性偏差。建议分别建立在轨和地面两套热模型,用于试验设计与在轨预测。
常见误区:忽略接触热阻会导致温度预测偏低;辐射表面属性只给发射率不给吸收率,会导致太阳辐照加热计算错误;瞬态分析时间步长过大,会平滑掉真实的温度波动,使热变形预测偏小。

综合来看,热分析建模的核心在于:明确传热路径,正确施加三类边界条件,准确输入四类热载荷,并根据分析目标选择稳态或瞬态求解。只有边界条件与载荷输入可信,后续的热变形与光学性能耦合分析才有意义。