2. 轨道力学基础与任务约束接口

轨道选择不是凭经验“挑一条好看的曲线”,而是把任务需求翻译成一组可计算、可验证的力学参数。二体问题、开普勒轨道根数、位置速度转换、轨道周期与能量,正是这套翻译过程的四个接口:前两者回答“轨道长什么样”,后两者回答“卫星此刻在哪里、飞得多快、需要多少能量”。本节围绕这四个接口展开,并给出可直接用于后续轨道选择计算的参数约定与示例。

一、二体问题:把复杂引力场简化为可解模型

在惯性参考系中,忽略其他天体摄动、大气阻力与光压,只考虑中心天体与航天器之间的万有引力,就得到二体问题。设中心天体质量为 M,航天器质量为 m,位置矢量为 r,则运动方程为:

r̈ = -μ · r / |r|³
μ = G·M   (地球 μ ≈ 3.986004418×10¹⁴ m³/s²)

该方程存在解析解,轨道是圆锥曲线:e=0 为圆,0<e<1 为椭圆,e=1 为抛物线,e>1 为双曲线。对绝大多数地球卫星任务,我们工作在椭圆或圆轨道区间。二体模型的价值在于:它给出了轨道选择的“基准解”,后续所有摄动分析、保持策略和燃料预算,都是在这个基准上叠加修正量。因此,本节所有公式与示例均默认采用地心惯性系(ECI),角度以春分点为起点逆时针度量。

二、开普勒轨道根数:六个参数锁定一条轨道

二体问题的解需要六个独立常数,工程上最常用的是开普勒六根数:

  • 半长轴 a:决定轨道尺度与周期,单位 m 或 km;
  • 偏心率 e:决定轨道形状,0≤e<1;
  • 轨道倾角 i:轨道面与赤道面夹角,0°≤i≤180°;
  • 升交点赤经 Ω:从春分点沿赤道量到升交点的角度;
  • 近地点幅角 ω:从升交点沿轨道面量到近地点的角度;
  • 真近点角 ν:从近地点量到航天器当前位置的角度,是随时间变化的量。

前五个根数确定轨道在空间中的几何位置,第六个根数确定航天器在轨道上的瞬时位置。任务约束通常直接映射到这些根数上:覆盖纬度范围约束倾角,回归特性约束半长轴与倾角组合,太阳同步约束升交点赤经的进动率,冻结轨道约束偏心率与近地点幅角。因此,轨道选择的第一步,就是把任务指标写成对 (a, e, i, Ω, ω) 的约束区间。

开普勒轨道根数几何关系与位置速度转换接口 地球 赤道面 轨道面 升交点 近地点 远地点 航天器 r v 2a i Ω ω ν 位置速度转换接口: (a, e, i, Ω, ω, ν) → r, v (由根数求位置速度) r, v → (a, e, i, Ω, ω, ν) (由位置速度反求根数)

三、位置速度转换:从根数到状态矢量

轨道选择过程中,经常需要在“根数”和“位置速度”之间来回转换。例如,用根数描述设计轨道,但用位置速度做数值积分或碰撞预警。转换的核心是先在轨道面内建立近焦点坐标系(PQW),再旋转到地心惯性系。

设近焦点坐标系中位置为 rPQW = [r·cosν, r·sinν, 0]T,其中 r = a(1-e²)/(1+e·cosν)。速度在 PQW 中为:

v_PQW = sqrt(μ/a) · [ -sinν, e + cosν, 0 ]^T / sqrt(1-e²)

然后通过三次旋转得到 ECI 下的状态:先绕 z 轴转 ω,再绕 x 轴转 i,最后绕 z 轴转 Ω。写成矩阵形式为 rECI = Rz(Ω)·Rx(i)·Rz(ω)·rPQW。反向转换则先由 r×v 得到角动量矢量 h,由 h 确定 i 和 Ω,再由偏心率矢量 e 确定 e 和 ω,最后由 ν 确定位置。

以一条近地圆轨道为例:a = 6878 km,e = 0.001,i = 98°,Ω = 120°,ω = 90°,ν = 0°。此时 r ≈ 6871 km,v ≈ 7.61 km/s。若把 ν 改为 180°,位置和速度方向都会改变,但 a、e、i、Ω、ω 保持不变。这说明:前五个根数描述轨道本身,第六个根数只描述“此刻在轨道上的哪里”。

四、轨道周期与能量:任务约束的量化接口

轨道周期由半长轴唯一决定:

T = 2π · sqrt(a³/μ)

对 a = 6878 km,T ≈ 5677 s ≈ 94.6 min。若任务要求每天绕地球 15 圈,则 T ≈ 5760 s,对应 a ≈ 6940 km,高度约 562 km。这就是周期约束直接反推半长轴的过程。

能量方面,单位质量机械能为:

ε = v²/2 - μ/r = -μ/(2a)

圆轨道 ε = -μ/(2a),椭圆轨道 ε 为负且与 a 一一对应。能量决定了轨道“绑定”还是“逃逸”:ε < 0 为椭圆或圆,ε = 0 为抛物线,ε > 0 为双曲线。在轨道选择中,能量差 Δε 直接对应速度增量 Δv 的量级,是燃料预算的起点。

参数约定:本节所有角度均以春分点为起点、逆时针为正;倾角 i 从赤道面量起,0° 为赤道顺行,90° 为极轨道,180° 为赤道逆行。后续章节若切换参考系,会明确说明。

五、把任务约束翻译成根数区间

下面给出一份可直接用于轨道选择初筛的参数映射表。表中“任务约束”是工程语言,“根数约束”是力学语言,两者之间的转换就是本节四个接口的综合应用。

任务约束主要影响根数典型取值示例判断依据
全球覆盖ii ≈ 90° 或太阳同步星下点轨迹覆盖纬度范围
回归轨道a, ia ≈ 6940 km, i ≈ 98°相邻轨迹重复周期
太阳同步a, e, i, Ωi ≈ 98°, a ≈ 7000 kmΩ 进动率 ≈ 0.9856°/天
冻结轨道e, ωe ≈ 0.001, ω ≈ 90°近地点高度长期稳定
快速重访a, ea 减小,e 增大近地点角速度增大

以太阳同步轨道为例,任务要求降交点地方时保持 10:30,则升交点赤经的进动率需约为 +0.9856°/天。由 J2 摄动公式可知,该进动率与 a、e、i 相关。取 e ≈ 0.001、a ≈ 7000 km,可反解出 i ≈ 98.2°。这就是从“地方时约束”到“倾角约束”的完整翻译链:任务指标 → 进动率 → 根数组合 → 位置速度验证。

常见误区:把“轨道高度”直接当作半长轴使用。高度 h = a - Re,其中 Re ≈ 6378 km。若任务要求 500 km 圆轨道,则 a = 6878 km,而不是 500 km。周期、能量、速度全部由 a 决定,不是由高度直接决定。

六、本节小结与后续接口

二体问题提供了可解析的基准轨道;开普勒六根数把轨道几何与瞬时位置分离;位置速度转换让根数与数值计算互通;周期与能量把任务指标量化为半长轴和速度增量。四者共同构成后续轨道选择的力学底座。下一阶段讨论摄动、轨道保持与发射窗口时,所有修正量都将叠加在本节给出的基准解之上。