3. 航天器热平衡基本方程

航天器热分析的一切计算，最终都归结为对能量守恒的求解。无论结构多复杂、轨道多极端，热控工程师真正反复书写的，只有一条核心关系：某个对象在某一时刻储存能量的变化率，等于它接收的所有热流减去它排出的所有热流。把这条关系写成数学形式，就是热平衡方程；把它按部件离散成一组方程，就是节点热平衡方程组；再根据是否保留时间项，就得到稳态与瞬态两种描述。

一、能量守恒原理与控制体

热平衡方程的物理基础是热力学第一定律。取航天器上一个有限区域作为控制体（例如一块蜂窝板、一台行波管、一组蓄电池），在时间间隔 dt 内，控制体内能的变化等于净流入的热量与内部产热之和。写成微分形式：

dU/dt = Q_in - Q_out + Q_gen

其中 U 为控制体内能，Q_in 为外部输入热流（太阳辐射、地球反照、行星红外、其他部件辐射、加热器功率），Q_out 为向外排出的热流（自身辐射、向其他部件辐射、向冷空间辐射），Q_gen 为内部产热（电子设备功耗、电池充放电损耗、推进剂化学反应热等）。工程上通常把 Q_in 与 Q_gen 合并为“得热项”，把 Q_out 称为“排热项”。

需要特别强调符号约定：本课程统一采用“流入控制体为正、流出控制体为负”的约定。太阳辐射、反照、行星红外、加热器功率、设备功耗均为正；向冷空间的辐射、向低温部件的净辐射、制冷机制冷量均为负。若中途切换参照系或改变符号约定，必须在计算文件中显式声明，否则节点方程会出现整体符号错误。

二、从连续体到节点热平衡方程

真实航天器是连续体，温度分布是空间和时间的连续函数。工程上无法直接求解连续偏微分方程，于是采用集总参数法：把每个等温化良好的部件或区域抽象为一个节点，节点内部温度均匀，节点之间通过热导（导热、辐射、对流）连接。这样，连续问题就离散为一组常微分方程。

对第 i 个节点，其热平衡方程可写为：

C_i · dT_i/dt = Q_i^env + Q_i^gen + Σ_j G_ij·(T_j - T_i) - ε_i·σ·A_i·T_i^4

式中 C_i 为节点 i 的热容（J/K），T_i 为节点温度（K），Q_i^env 为外部环境直接投射到该节点的热流，Q_i^gen 为节点内部产热，G_ij 为节点 i 与节点 j 之间的热导（W/K），ε_i、A_i 分别为节点 i 的发射率和有效辐射面积，σ 为斯特藩-玻尔兹曼常数（5.67×10⁻⁸ W/(m²·K⁴)）。

这个方程包含三类项：热容项（瞬态时保留，稳态时为零）、线性导热/辐射耦合项、非线性四次方辐射项。正是四次方项的存在，使热平衡方程成为非线性方程组，通常需要迭代求解。

三、稳态热平衡的数学描述

当航天器处于轨道热环境周期性稳定、或工作模式长时间不变时，各节点温度不再随时间变化，dT_i/dt=0，瞬态项消失，方程变为：

Q_i^env + Q_i^gen + Σ_j G_ij·(T_j - T_i) - ε_i·σ·A_i·T_i^4 = 0

稳态方程描述的是“得热与排热恰好相等”的平衡状态。它没有时间信息，只回答“最终稳定在多少度”。稳态分析常用于：确定散热面面积、校核高温工况和低温工况下的温度包络、验证加热器功率是否足够。求解时通常给定环境热流和内部功耗，未知量为各节点温度，采用牛顿-拉夫逊迭代或高斯-赛德尔迭代逐次逼近。

四、瞬态热平衡的数学描述

当航天器进出阴影、设备开关机、变轨或姿态机动时，热环境与内热源随时间剧烈变化，必须保留热容项：

C_i · dT_i/dt = Q_i^env(t) + Q_i^gen(t) + Σ_j G_ij·(T_j - T_i) - ε_i·σ·A_i·T_i^4

瞬态方程是一组一阶非线性常微分方程，需要给定初始温度场 T_i(0)，然后沿时间步进求解。常用方法包括显式欧拉、四阶龙格-库塔以及隐式向后差分。显式方法受稳定性限制，时间步长需满足 Δt < C_i / ΣG_ij 量级；隐式方法虽每步需迭代，但允许更大步长，适合长周期轨道分析。

瞬态分析的核心价值在于捕捉温度变化速率和极值。例如某设备在阴影期最低温度、出阴影后升温速率、加热器启动后的响应滞后，都只能由瞬态方程给出。稳态结果只是瞬态过程的渐近终点。

五、稳态与瞬态的对比与选用

实际工程中两者配合使用：先用稳态分析快速确定设计点和散热能力，再用瞬态分析验证极端工况下的温度是否越限、加热器功率是否足够、温度变化率是否超出部件允许范围。稳态是瞬态的极限，瞬态是稳态在时间轴上的展开。

理解这张图，就抓住了热平衡分析的主线：左侧是得热，右侧是排热，中间是节点自身的热容。稳态时热容项为零，方程退化为代数方程；瞬态时热容项决定了温度变化的快慢。节点之间的热导 G_ij 把单个节点方程耦合为整个航天器的方程组，求解这组方程，就得到了各节点的温度分布与时间历程。