4. 遥感载荷热物理模型构建
热物理模型是热控设计与稳定性分析的“数字孪生底座”。遥感载荷在轨工作时,其温度场由太阳辐照、地球红外、内部功耗、辐射散热与结构导热共同决定。若模型构建失真,后续的热控参数优化、热变形预测和稳定性评估都会失去意义。本章围绕几何简化、材料热物性、接触热阻、多层隔热组件等效建模与热网络模型搭建五个关键环节,说明如何从物理对象逐步得到可用于计算的热网络模型。
一、几何简化:从真实结构到热分析模型
遥感载荷结构复杂,包含镜筒、焦面组件、探测器、遮光罩、支撑桁架与电子学箱等。直接对完整CAD模型划分热网格,计算量巨大且大量细节对温度场影响甚微。几何简化的原则是“保热流路径、保辐射面形、保热容主体”。
具体操作上,可将薄壁镜筒简化为等厚度壳体,保留其外表面辐射面积与轴向导热截面;将螺钉、接插件、线缆等小尺度特征合并为等效热质量或忽略;对遮光罩等以辐射换热为主的结构,保留其开口角系数与内表面发射率,而不必保留全部加强筋。简化后应核对三项守恒:总质量偏差一般控制在5%以内,主要散热面积偏差控制在10%以内,关键传热路径的热阻量级不发生数量级变化。
二、材料热物性:温度相关的输入参数
热网络模型中每个节点都需要导热系数λ、比热容c、密度ρ与发射率ε。遥感载荷常用材料包括铝合金、钛合金、碳纤维复合材料、铍与光学玻璃。这些材料的λ和c通常随温度变化,尤其在低温段变化显著。例如某铝合金在20℃时λ约为160 W/(m·K),在-50℃时可降至约120 W/(m·K);碳纤维复合材料沿纤维方向与垂直方向导热系数可相差一个数量级。
建模时应优先采用材料供应商提供的温度相关数据,并以表格或分段线性函数输入求解器。若仅有室温单点值,应在-40℃~+60℃范围内做敏感性分析,确认温度场对λ变化的敏感度。表面发射率ε同样关键:抛光铝ε约0.05,阳极氧化铝约0.8,镀金表面约0.03。多层隔热组件外表面通常取低发射率,内表面取高发射率。
三、接触热阻:被忽视的温降来源
两个固体表面即使压紧,实际接触仅发生在微观凸峰处,其余间隙由空气或真空填充,形成接触热阻Rc。在真空环境中,间隙无气体导热,接触热阻往往成为局部温降的主因。接触热阻可表示为:
R_c = 1 / (h_c · A)
h_c = 1 / (1.25 · σ · (m/σ)^0.95 · (P/H)^0.85 · k_s)
其中hc为接触导热系数,A为名义接触面积,σ为表面粗糙度,m为斜率,P为接触压力,H为较软材料硬度,ks为等效导热系数。工程上更常用经验值:铝-铝干接触hc约1000~5000 W/(m²·K),加导热硅脂后可提升至5000~15000 W/(m²·K),加铟箔可达10000~20000 W/(m²·K)。
建模时,接触热阻应作为独立热阻串联在对应节点之间,而不是简单合并到材料导热中。对于探测器与冷平台、光学镜与镜座等关键界面,建议通过实验或供应商数据标定hc,并在模型中保留其温度依赖性。
四、多层隔热组件等效建模
多层隔热组件(MLI)由多层低发射率反射屏与间隔层交替叠合而成,是遥感载荷最常用的被动热控措施。其等效热阻极大,但并非无限。MLI的等效传热包括辐射、残余气体导热与固体导热三部分。工程上常用等效发射率εeff或等效热导率keff描述:
Q = ε_eff · σ · A · (T_h^4 - T_c^4)
ε_eff ≈ 0.01 ~ 0.05(真空、多层、低发射率)
在热网络模型中,MLI通常被等效为一个高阻值热阻RMLI,连接被包覆结构与外部环境节点。若MLI包覆面积大,还应考虑其边缘和开口处的“热短路”效应,这些位置往往成为漏热主通道。建模时可将MLI分为主体区域与边缘区域,分别赋予不同等效热阻。
五、热网络模型搭建:节点、热阻与热容
热网络模型将连续温度场离散为若干等温节点,节点之间用热阻连接,每个节点配置热容。其数学形式为:
C_i · dT_i/dt = Q_i + Σ_j (T_j - T_i) / R_ij
其中Ci为节点热容,Qi为节点内热源或外热流,Rij为节点间热阻。稳态分析时左侧为零,求解线性方程组即可得到各节点温度;瞬态分析时需时间积分。
节点划分应遵循“热流集中处加密、温度梯度小处合并”的原则。例如探测器焦面、镜体、镜座、遮光罩、电子学箱、散热面应分别设节点;而结构支架若导热良好且温差小,可合并为单一节点。热阻类型包括导热热阻、辐射热阻与接触热阻,三者量级可能相差很大,建模时应分别计算并标注来源。
六、建模流程与结果判读
完整的热物理模型构建可按以下顺序推进:先依据任务剖面确定外热流与内热源,再完成几何简化与节点划分,随后输入温度相关材料热物性,逐对界面计算接触热阻,对MLI区域赋予等效热阻,最后组装热网络并求解。稳态结果应检查各节点温度是否落在合理区间,例如探测器焦面通常控制在-20℃~+40℃,光学镜体温度梯度应小于2℃~5℃。瞬态结果应关注升温/降温速率与时间常数,判断是否满足成像任务对热稳定性的要求。
若某节点温度异常偏高或偏低,应优先排查三类问题:接触热阻是否被低估或高估、MLI边缘漏热是否被忽略、辐射角系数是否因几何简化而失真。通过逐项修正并复算,使模型逐步逼近真实热行为,才能为后续热控设计与稳定性分析提供可靠依据。