3. 航天器动力学与轨道要素摄动

轨道控制的燃料优化，本质上是在“轨道会自己变”和“我们必须让它按需变”之间做权衡。要理解这个权衡，先要建立两个基础：一是理想情况下轨道如何描述，二是真实环境中轨道如何被摄动改变。前者是二体问题与轨道六根数，后者是J2项等主要摄动，而两者之间的差值，正是需要燃料去补偿的部分。

二体问题与轨道六根数

二体问题假设航天器只受中心天体引力作用，忽略其他天体、大气阻力、太阳光压和非球形引力。在这个假设下，轨道是严格的圆锥曲线，且六个轨道根数保持不变。这六个根数分别是：半长轴 a、偏心率 e、轨道倾角 i、升交点赤经 Ω、近地点幅角 ω、真近点角 ν（或平近点角 M）。

其中，a 和 e 决定轨道形状与大小，i、Ω、ω 决定轨道在空间中的姿态，ν 决定航天器在轨道上的瞬时位置。前五个根数描述轨道本身，第六个描述航天器沿轨道的相位。理解这一点很关键：轨道控制通常改变的是前五个根数中的某一个或某几个，而相位调整则更多与任务时序相关。

在二体模型下，轨道能量 ε = -μ/(2a)，角动量 h = √(μa(1-e²))。这两个量决定了轨道的基本性质。半长轴越大，能量越高；偏心率越大，轨道越扁。这些关系是后续分析摄动和燃料代价的起点。

J2项摄动：最主要的地球非球形影响

真实地球不是理想球体，赤道略鼓、两极稍扁。这种扁率主要由二阶带谐项 J2 描述，其值约为 1.08263×10⁻³。J2项虽然量级不大，但长期累积效应显著，是低地球轨道和太阳同步轨道设计中最关键的摄动源。

J2项对轨道根数的影响具有明显的选择性：它不改变半长轴 a 和偏心率 e 的长期平均值，但会引起升交点赤经 Ω、近地点幅角 ω 和平近点角 M 的长期漂移。其中最重要的两个漂移率是：

dΩ/dt = -3/2 · J2 · (Re/p)² · n · cos i
dω/dt =  3/4 · J2 · (Re/p)² · n · (5cos²i - 1)

这里 Re 是地球赤道半径，p = a(1-e²) 是半通径，n = √(μ/a³) 是平均角速度，i 是轨道倾角。这两个公式揭示了几个重要规律：

• dΩ/dt 与 cos i 成正比。当 i < 90°（顺行轨道）时，升交点赤经向西漂移；当 i > 90°（逆行轨道）时，升交点赤经向东漂移。太阳同步轨道正是利用这一特性，选择特定倾角使 dΩ/dt 恰好等于地球绕太阳公转的平均角速度（约 0.9856°/天）。
• dω/dt 的符号由 5cos²i - 1 决定。当 i ≈ 63.43° 时，dω/dt = 0，这就是著名的“临界倾角”，常用于冻结轨道设计，使近地点幅角长期稳定。
• 漂移率与 a 的 -7/2 次方相关（因为 n ∝ a^(-3/2)，(Re/p)² ∝ a^(-2)）。轨道越低，J2漂移越快。例如，500 km高度的太阳同步轨道，dΩ/dt 约为 1°/天量级；而地球静止轨道（约35786 km）的J2漂移则小得多。

其他主要摄动及其影响

除J2项外，还有几类摄动需要在轨道控制中考虑：

• 大气阻力：主要影响低地球轨道（通常低于800 km）。阻力使半长轴 a 持续衰减，轨道高度下降，进而影响任务寿命。阻力加速度约为 a_drag = -½ ρ v² C_d A/m，其中 ρ 是大气密度，C_d 是阻力系数，A/m 是面质比。太阳活动高年，大气密度上升，衰减可加快数倍。
• 第三体引力：日月引力对高轨（如地球静止轨道、地月转移轨道）影响显著，主要引起轨道倾角和偏心率的长周期变化。
• 太阳光压：对大面积、低面质比航天器影响较大，尤其在同步轨道附近，可引起偏心率矢量缓慢旋转。

这些摄动的共同特点是：它们不消耗燃料，却持续改变轨道。如果任务要求轨道保持稳定，就必须用推进系统定期修正，这就是摄动补偿。

摄动补偿的燃料代价

摄动补偿的燃料代价，可以用速度增量 Δv 来度量。根据齐奥尔科夫斯基公式，Δv = Isp · g₀ · ln(m₀/m_f)，其中 Isp 是比冲，g₀ 是标准重力加速度，m₀ 和 m_f 分别是初始和最终质量。因此，Δv 直接对应燃料消耗。

以太阳同步轨道为例，若J2引起的升交点赤经漂移率为 0.9856°/天，而任务要求轨道面保持固定，则每年需要补偿约 360° 的赤经漂移。但实际太阳同步轨道正是利用J2漂移来匹配太阳公转，所以不需要补偿 Ω。真正需要补偿的往往是大气阻力引起的半长轴衰减。假设500 km轨道，大气阻力使半长轴每天衰减约 100 m，一年约 36.5 km。为恢复高度，需要施加沿速度方向的 Δv，量级约为 Δv ≈ (Δa/a) · v/2。取 a ≈ 6878 km，v ≈ 7.6 km/s，则每年约需 20 m/s 的 Δv。若卫星质量1000 kg，Isp = 300 s，则每年燃料消耗约 6.8 kg。这个量级说明：低轨轨道保持的燃料代价主要来自大气阻力补偿，而非J2项本身。

对于地球静止轨道，J2项引起的经度漂移需要定期修正。地球静止轨道在J2作用下，升交点赤经和近地点幅角漂移，导致星下点经度发生周期性摆动。典型的经度保持精度为 ±0.1°，每年需要约 50 m/s 的 Δv 用于东西向和南北向保持。其中南北向保持（倾角修正）消耗最大，约占总量的 80% 以上。这也是为什么电推进等高效推进方式在地球静止轨道卫星上越来越重要。

从图中可以清晰看到，二体问题提供理想轨道基准，J2项和其他摄动使轨道根数发生漂移，而燃料补偿只针对那些与任务要求冲突的漂移量。理解每个摄动对哪个根数、以什么速率、朝什么方向产生影响，是制定轨道保持策略和估算燃料预算的前提。在实际任务中，通常先用轨道递推模型（如STK、GMAT或自编的J2解析模型）预测未来一段时间的根数变化，再根据任务约束确定需要补偿的量和时机，最后折算成Δv和燃料质量。