📚 课程目录第4章太阳翼发电模型

4. 太阳翼发电模型

太阳翼是卫星在轨运行的主电源，其发电能力直接决定了整星可用功率的上限，也决定了任务排程中“能安排多少负载、能持续多久”的边界。要做出可信的能源约束排程，第一步就是把太阳翼的发电能力从物理过程抽象成一个可计算、可复现的模型。本节围绕理想发电功率、衰减因子、温度效应、光照角度损失四个核心环节，给出工程上常用的简化计算方法，并说明每个参数的物理含义与取值依据。

一、理想发电功率：从太阳常数到电池阵输出

理想发电功率指太阳翼在标准测试条件（AM0 光谱、光照强度 1367 W/m²、电池温度 28℃）下，不考虑任何衰减、温度与角度损失时的理论输出。其计算起点是太阳常数与电池阵有效面积：

P_ideal = S × A × η_cell × η_pack × η_mppt
其中：
  S        = 1367 W/m²（AM0 太阳常数）
  A        = 电池阵有效受照面积（m²）
  η_cell   = 单体电池光电转换效率（如 GaAs 三结约 0.28~0.30）
  η_pack   = 布片系数（有效布片面积 / 阵面总面积，典型 0.85~0.92）
  η_mppt   = MPPT 最大功率点跟踪效率（典型 0.95~0.98）

举例说明：某卫星太阳翼单翼阵面面积 6.0 m²，采用 GaAs 三结电池，η_cell 取 0.29，布片系数 0.90，MPPT 效率 0.97。则单翼理想功率约为 1367 × 6.0 × 0.29 × 0.90 × 0.97 ≈ 2077 W。这个数值是后续所有修正的基准，实际在轨输出会在此基础上逐项打折。

二、衰减因子：把“新电池”折算成“在轨电池”

电池阵在寿命初期（BOL）与寿命末期（EOL）的输出差异，主要由三类衰减构成：紫外辐照与粒子辐照引起的电池性能退化、微流星与空间碎片造成的局部损伤、以及长期热循环导致的接触电阻上升。工程上通常用一个综合衰减因子 D 表示：

D = D_rad × D_uv × D_thermal × D_contamination
典型取值（15 年 GEO 任务）：
  D_rad          ≈ 0.85 ~ 0.90（辐照衰减）
  D_uv           ≈ 0.97 ~ 0.98（紫外衰减）
  D_thermal      ≈ 0.98 ~ 0.99（热循环）
  D_contamination≈ 0.98 ~ 0.99（污染）
  综合 D         ≈ 0.80 ~ 0.86

需要强调的是，衰减因子是随任务时间单调下降的，排程时若按 EOL 取值会过于保守，按 BOL 取值又会高估能力。常用做法是取任务中期（MOL）或按任务阶段分段取值，并在能源裕度中单独留出衰减余量。

三、温度效应与光照角度损失

电池效率随温度升高而下降，GaAs 三结电池的温度系数约为 -0.20%/℃ ~ -0.25%/℃。若在轨工作温度比标准 28℃ 高出 40℃，则效率损失约 8%~10%。温度修正可写为：

k_T = 1 + α_T × (T_cell - 28)
α_T ≈ -0.0022 /℃（GaAs 三结典型值）

光照角度损失来自太阳矢量与太阳翼法线之间的夹角 θ。理想情况下 θ=0，输出最大；当 θ 增大时，有效受照面积按 cosθ 减小。此外，太阳翼对日定向机构的跟踪误差、卫星姿态偏置也会引入额外角度偏差。工程上把角度损失统一写成：

k_θ = cosθ × k_track
k_track ≈ 0.98 ~ 0.99（对日跟踪精度损失）

当 θ 超过 60° 时，cosθ 已降至 0.5 以下，此时太阳翼输出急剧下降，排程中应尽量避免让卫星长时间处于大角度偏置姿态。

四、工程简化计算方法与参数关系

把上述环节串联，得到工程上最常用的太阳翼发电简化模型：

P_out = P_ideal × D × k_T × k_θ × k_shadow
其中：
  k_shadow = 1（全日照）或 0（地影区）
  地影区占比由轨道 β 角与轨道高度决定

该模型把复杂的物理过程压缩为五个可查表、可标定的系数，便于在任务排程中快速迭代。实际使用时，建议按轨道周期分段计算：光照段用 k_shadow=1，地影段用 k_shadow=0，再对光照段内的 θ 与 T_cell 取平均值。若需要更高精度，可把 θ 与 T_cell 按时间步长离散，逐点计算后积分得到每圈发电量。

工程经验：在初步排程阶段，可先用 P_out ≈ P_ideal × 0.75 作为光照段平均输出估算值，再根据任务阶段和姿态计划做细调。这个 0.75 综合了衰减、温度、角度与跟踪损失的典型中值，仅用于快速估算，不能替代正式计算。

注意：衰减因子 D 与温度系数 α_T 的符号约定必须全文统一。本文约定 D 为小于 1 的正系数，α_T 为负值，k_T = 1 + α_T×(T-28)。若采用其他约定，需在计算表中明确标注，避免排程结果出现方向性错误。

把模型落到排程上，关键是把“每圈光照段可用发电量”算清楚。以 6.0 m² 单翼、P_ideal=2077 W 为例，取 D=0.83、k_T=0.91、k_θ=0.85，则光照段瞬时输出约 2077×0.83×0.91×0.85 ≈ 1333 W。若轨道周期 95 min、地影占比 35%，则每圈光照时长约 61.75 min，单圈发电量约 1333 W × 61.75/60 h ≈ 1372 Wh。这个数值就是该圈可分配给载荷、热控、数传等负载的总能量上限，再扣除蓄电池充放电效率与母线损耗后，才是真正可用于任务排程的净能量。

最后需要提醒：上述模型是工程简化模型，适用于方案论证与初步排程。若进入详细任务规划阶段，应结合太阳翼实际 I-V 曲线、MPPT 工作点、蓄电池 SOC 约束与姿态计划，做逐时步的联合仿真，才能得到更贴近在轨的发电能力估计。