3. 热平衡方程与能量流网络构建
热平衡分析的本质，是把一个连续的温度场问题离散成一组可求解的能量守恒方程。要做到这一点，第一步不是写方程，而是划定控制体——明确"我们研究的是哪一块、边界在哪里、哪些能量可以穿过边界"。控制体划分决定了后续所有热阻、节点和残差的物理含义，一旦划错，后面的计算再精确也没有意义。
控制体划分：从物理对象到计算域
控制体是热分析的基本单元。划分时遵循三条原则：材料与结构连续（同一控制体内导热路径单一）、边界可描述（每个界面能明确是导热、对流还是辐射）、热源归属清晰（每个发热元件只属于一个控制体，避免重复计入）。
以某星载电子设备机箱为例：内部有功率放大器（发热约 25 W）、电源模块（发热约 12 W）、结构框架和外壳。合理的划分是：功放模块为一个控制体，电源模块为一个控制体，机箱壳体为一个控制体，安装底板为一个控制体。每个控制体有唯一的节点温度 Ti，控制体之间的连接就是热阻路径。
需要特别说明符号约定：本文所有热流 Q 以流入控制体为正，流出为负；温度以摄氏度为统一单位；热阻 R 定义为温差与热流之比，即 R = ΔT/Q，单位 K/W。全文及图示均采用这一约定。
三类传热路径的建模方式
控制体之间的能量传递只有三种机制，对应三种热阻表达：
- 导热路径：Rcond = L/(kA)，其中 L 为路径长度，k 为材料导热系数，A 为截面积。例如铝合金框架 k ≈ 180 W/(m·K)，一段长 0.05 m、截面 0.002 m² 的导热路径，热阻约 0.139 K/W。
- 对流路径：Rconv = 1/(hA)，h 为对流换热系数。自然对流下 h 约 5～10 W/(m²·K)，强制风冷可达 30～100 W/(m²·K)。
- 辐射路径：Rrad = 1/(hrA)，其中辐射换热系数 hr = εσ(T1²+T2²)(T1+T2)，ε 为发射率，σ = 5.67×10⁻⁸ W/(m²·K⁴)。注意辐射热阻随温度变化，这是它与前两者的关键区别，迭代求解时必须更新。
热阻网络与节点方程
把每个控制体画成一个节点，把每条传热路径画成一个热阻，就得到热阻网络。对任意节点 i，稳态能量守恒方程为：
Σ_j (T_j - T_i) / R_ij + Q_i = 0
其中：
T_j 相邻节点温度 (K 或 ℃)
T_i 当前节点温度
R_ij 节点 i 与 j 之间的等效热阻 (K/W)
Q_i 节点 i 的净发热功率 (W)，流入为正对 n 个节点可写出 n 个方程，构成线性方程组（辐射热阻线性化后）。以三节点系统为例，节点 1 为功放、节点 2 为壳体、节点 3 为环境（边界节点，温度已知 25 ℃）：
节点1: (T2 - T1)/R12 + Q1 = 0
节点2: (T1 - T2)/R12 + (T3 - T2)/R23 = 0
节点3: T3 = 25 ℃ (边界条件)
设 Q1 = 25 W, R12 = 0.5 K/W, R23 = 1.2 K/W
由节点1: T1 = T2 + Q1·R12 = T2 + 12.5
由节点2: (T1 - T2)/0.5 = (T2 - 25)/1.2
代入: 12.5/0.5 = (T2 - 25)/1.2 → 25 = (T2 - 25)/1.2
解得: T2 = 55 ℃, T1 = 67.5 ℃这组结果说明：功放节点温度 67.5 ℃，壳体 55 ℃，环境 25 ℃。若功放结温上限为 85 ℃，则当前设计有约 17.5 ℃ 余量。
能量守恒校核与残差分析
求解完成后必须做全局能量守恒校核：所有热源总功率应等于所有边界路径散出的总功率。上例中总发热 25 W，壳体到环境的热流为 (55-25)/1.2 = 25 W，完全平衡，残差为 0。
实际工程中残差不会严格为零，需要设定判据。常用做法是计算相对残差：
残差 = |ΣQ_in - ΣQ_out| / ΣQ_in × 100%
判据建议：
< 1% 可接受，网格与参数合理
1%～5% 需检查边界条件或热阻取值
> 5% 不可接受，重新划分控制体或修正模型残差偏大通常来自三个原因：控制体边界重叠导致热源重复计入、辐射热阻未随温度迭代更新、或某个传热路径被遗漏（如忽略了安装面的接触热阻）。排查时逐节点列出流入流出清单，比整体看总账更容易定位问题。
上图把控制体、热阻路径、节点方程和守恒校核放在同一张网络里，可以直观看到：能量从功放节点经导热+接触热阻流向壳体，再从壳体经对流+辐射热阻散入环境，每一步的热流大小相等，这正是稳态热平衡的物理含义。
完成上述步骤后，一套可求解、可校核的热平衡模型就建立起来了。它的价值不仅在于算出几个温度值，更在于当残差超标时，能沿着热阻网络逐段定位是哪个环节的参数或边界条件出了问题——这正是后续温升控制和热稳定性验证的基础。