3. 航天器热平衡基本方程与能量守恒
航天器在轨运行时，其温度分布本质上由能量守恒关系决定。要判断某个部件是否会过热、某个表面是否会过冷，或者某台单机的温度波动是否在允许范围内，都必须回到热平衡这一基本物理约束上来。本节从热平衡方程的建立出发，讨论稳态与瞬态两种能量守恒形式，并引入节点网络法作为工程求解的基础框架。
一、热平衡方程的物理含义
对航天器上任意一个控制体（可以是整星、一个舱段、一块单机安装板，也可以是一台具体设备），能量守恒可表述为：单位时间内进入控制体的热量，减去离开控制体的热量，等于控制体内部储能的变化率。写成通用形式：
Q_in - Q_out = dU/dt其中 Q_in 包括太阳辐照吸收、地球红外与反照吸收、内部功耗、相邻节点导热或辐射输入等；Q_out 包括向深空辐射散热、向低温面辐射、导热流出等；dU/dt 为控制体热容与温度变化率的乘积。当 dU/dt = 0 时，系统进入稳态热平衡，此时净输入热量为零，温度不再随时间变化。
需要特别强调符号约定：本文统一采用“流入为正、流出为负”的约定，辐射散热项以负号出现在方程左侧。若切换为“散热为正”的写法，所有项的符号需整体调整，但物理结论不变。工程计算中符号混乱是常见错误来源，建立方程前必须先固定约定。
二、稳态与瞬态能量守恒
稳态热平衡用于评估航天器在某一固定工况（如长期对日定向、阴影区稳定运行）下的最终温度水平。此时方程退化为代数方程组：
Σ Q_in,i = Σ Q_out,j例如某单机安装面稳态时，内部功耗 20 W，通过安装板导热导出 12 W，其余 8 W 通过表面辐射排向舱内低温环境。稳态方程即 20 = 12 + 8，对应安装面温度稳定在某一数值。
瞬态热平衡则用于描述进出阴影、变轨、载荷开关机等温度随时间变化的过程。此时必须保留热容项：
C dT/dt = Q_in(t) - Q_out(T)其中 C 为节点热容（J/K），Q_out 通常与 T⁴ 相关，因此瞬态方程是非线性微分方程。工程上常关注两个指标：温度变化速率 dT/dt 和达到新稳态所需的时间常数 τ = C / (dQ_out/dT)。热容越大、散热越弱，温度响应越慢。例如某电池组热容约 8000 J/K，进出阴影时若净热流约 40 W，则初始温变率约 0.005 K/s，即约 0.3 K/min，这一量级可用于初步判断温度波动是否可接受。
三、节点网络法基础
节点网络法（Nodal Network Method）是航天器热分析的核心建模方法。其基本思想是：将航天器离散为有限个等温节点，节点之间通过导热、辐射和流体换热建立热连接，每个节点列出一个热平衡方程，最终形成方程组求解。
节点划分遵循“热物性相近、温度梯度可忽略、热流路径清晰”的原则。例如可将卫星本体划分为：本体结构节点、电池节点、载荷节点、散热面节点、多层隔热材料外表面节点等。节点之间连接关系包括：
- 导热连接：Q = kA(T_i - T_j)/L，适用于安装面、结构梁、热管等；
- 辐射连接：Q = σεFA_i(T_i⁴ - T_j⁴)，适用于舱内表面之间、散热面与空间；
- 外部热源：太阳、地球红外、地球反照，按轨道位置和表面光学参数计算。
对第 i 个节点，稳态方程可写为：
Σ_j k_ij (T_j - T_i) + Σ_j R_ij (T_j⁴ - T_i⁴) + Q_ext,i + Q_int,i = 0其中 k_ij 为导热系数，R_ij 为辐射交换系数，Q_ext,i 为外部热源，Q_int,i 为内部功耗。该方程组通常是非线性的，需迭代求解。
四、热平衡温度求解思路
求解热平衡温度一般按以下思路推进：
- 明确对象与边界：确定分析对象、轨道工况、姿态、表面涂层参数（α_s、ε）和内部功耗。
- 建立节点网络：划分节点，确定节点间导热与辐射连接，形成热阻网络。
- 列写方程：对每个节点写出稳态或瞬态热平衡方程，注意符号统一。
- 迭代求解：由于辐射项非线性，常用牛顿迭代或高斯-赛德尔迭代，先给初值，逐步逼近。
- 结果判读：检查各节点温度是否在允许范围内，热流是否合理，必要时调整涂层、散热面面积或热管布局。
以一个简化示例说明：某载荷节点内部功耗 30 W，仅通过一块 0.2 m² 的辐射面排热，表面发射率 ε = 0.85，环境等效温度约 250 K。稳态时：
30 = σ ε A (T⁴ - T_env⁴)
T⁴ = 30 / (5.67e-8 × 0.85 × 0.2) + 250⁴
T ≈ 322 K若该载荷允许温度为 273～313 K，则 322 K 已超限，需要增大辐射面积或降低功耗。这一判断过程正是热平衡方程的实际价值所在。
从上述过程可以看出，热平衡方程不仅是理论公式，更是热设计迭代的核心工具。建立方程时，节点划分的粗细、连接关系的取舍、外部热源的计算精度，都会直接影响温度求解结果。工程实践中，通常先以少量节点快速估算，再逐步细化模型，最终用专业热分析软件完成全星级稳态与瞬态计算。理解方程本身的结构和每一项的物理来源，是正确使用软件和判断结果合理性的前提。