2. 轨道控制基础与燃料消耗建模
轨道控制的核心问题可以归结为一句话：用尽可能少的推进剂，把航天器从当前轨道状态送到目标轨道状态。要回答“少到什么程度”，必须先建立三类基础认识——机动有哪些基本类型、速度增量如何换算成推进剂质量、以及推进系统的效率参数如何影响这个换算。本节围绕这三条主线展开，并给出可直接用于估算的建模方法。
一、轨道机动的基本类型
按作用效果划分，轨道机动通常分为以下几类：
- 轨道面内机动：改变轨道半长轴、偏心率或近远地点高度，典型如霍曼转移、共椭圆机动、近地点抬升。速度增量方向与速度矢量共线，效率最高。
- 轨道面改变机动：改变轨道倾角或升交点赤经，速度增量方向垂直于轨道面。在速度较高处改变倾角代价极大，因此工程上常采用“组合机动”，把面内加速与倾角改变合并到一次点火中。
- 轨道保持与位置保持：对抗摄动（大气阻力、J2项、太阳光压）引起的轨道漂移，如地球静止轨道的南北/东西保持、太阳同步轨道的降交点地方时保持。
- 相位调整与交会机动：调整沿迹方向相位角，用于编队飞行、交会对接的远距离导引段。
- 离轨与再入机动：降低近地点进入大气，或直接制动再入。
不同类型的机动，其速度增量方向与速度矢量的夹角不同，直接决定了“同样的Δv能换来多少轨道改变”。面内机动效率最高，面改变机动效率最低，这是后续优化的基本出发点。
二、齐奥尔科夫斯基方程与速度增量—推进剂质量关系
火箭推进的基石是齐奥尔科夫斯基方程：
Δv = Isp · g0 · ln(m0 / mf)
其中：
Δv —— 速度增量（m/s）
Isp —— 比冲（s）
g0 —— 标准重力加速度，取 9.80665 m/s²
m0 —— 点火前总质量（kg）
mf —— 点火后总质量（kg）把推进剂质量 mp = m0 − mf 代入，可改写为：
mp = m0 · [1 − exp(−Δv / (Isp · g0))]
推进剂质量分数：
mp / m0 = 1 − exp(−Δv / (Isp · g0))这个式子揭示了两个关键事实。第一，推进剂质量分数只取决于 Δv 与 Isp·g0 的比值，与航天器绝对质量无关；第二，当 Δv 接近或超过 Isp·g0 时，指数项迅速衰减，推进剂需求急剧上升。例如 Isp = 300 s 时，Isp·g0 ≈ 2942 m/s；若单次机动 Δv = 1500 m/s，则 mp/m0 ≈ 1 − exp(−0.51) ≈ 40%，也就是说近一半的初始质量要作为推进剂。
工程上常用“质量比” R = m0/mf 来评估：R = exp(Δv/(Isp·g0))。R 越大，说明结构效率压力越大。对于多级或多次机动任务，应逐次计算，把上一次的 mf 作为下一次的 m0，而不是把总 Δv 一次性代入。
三、比冲与推力器效率对燃料消耗的影响
比冲 Isp 是推进系统最核心的效率指标，定义为推力与推进剂质量流率之比：
Isp = F / (ṁ · g0)
F —— 推力（N）
ṁ —— 推进剂质量流率（kg/s）比冲越高，产生同样 Δv 所需的推进剂越少。以 1000 kg 航天器、Δv = 200 m/s 为例：
| 推进类型 | 比冲 Isp (s) | 推进剂质量 (kg) | 质量分数 |
|---|---|---|---|
| 冷气推进 | 60 | 约 288 | 28.8% |
| 单组元肼 | 220 | 约 89 | 8.9% |
| 双组元 | 310 | 约 64 | 6.4% |
| 电弧加热 | 500 | 约 40 | 4.0% |
| 离子推进 | 3000 | 约 6.8 | 0.68% |
表中数据为按上述公式估算的示例值，用于说明量级差异，不代表具体型号实测结果。可以看到，从冷气到离子推进，推进剂需求相差超过 40 倍。但高比冲往往伴随低推力，离子推进推力常在毫牛量级，完成同样 Δv 需要长时间连续工作，这对轨道保持类任务可以接受，对快速交会则不合适。
除比冲外，还需关注推力器效率链：电源功率、推力器效率、羽流损失、推力矢量偏斜等。实际可用 Δv 应按下式修正：
Δv_可用 = Isp · g0 · ln(m0/mf) · η
η 为综合效率，包含推力方向损失、重力损失、姿控耦合损失等，
工程估算常取 0.95～0.99，具体依任务剖面确定。四、建模流程与结果判读
把上述内容串成可执行的估算流程：
- 确定任务剖面，列出各次机动的 Δv 及方向类型。
- 选定推进系统，查取 Isp 与推力范围。
- 按顺序逐次代入 mp = m0[1 − exp(−Δv/(Isp·g0))]，更新 m0。
- 累加各次推进剂质量，加入 5%～15% 余量。
- 核对总质量是否满足发射质量约束，若不满足则回到第 1 步优化机动方案。
结果判读的关键是看“推进剂质量分数”和“质量比 R”。若单次机动质量分数超过 30%，说明该次机动代价过高，应考虑分解为多次小机动、利用引力辅助或改用更高比冲的推进方式。若总质量比 R 超过 3，则结构质量压力很大，需要重新审视任务设计。
上图把三类信息放在同一视图中：左侧是机动类型对 Δv 方向与效率的影响，中间是齐奥尔科夫斯基方程及其变形，右侧是比冲、推力、效率与余量参数，底部用条形图直观对比不同推进方式的推进剂需求。实际建模时，先由任务剖面确定 Δv，再由推进系统确定 Isp 与 η，代入方程逐次计算，最后用余量系数修正并核对质量约束。这样得到的推进剂估算值，既可用于方案比较，也可作为后续优化章节的输入基线。