2. 轨道动力学基础回顾
轨道保持与漂移控制的全部推导,最终都要落回一个基本问题:航天器在某一时刻处于哪里、速度指向何方、未来会漂到何处。要回答这个问题,必须先建立一套自洽的动力学描述体系。本节把二体问题、轨道根数、摄动运动方程与参考坐标系四件事串成一条链:二体问题给出理想运动的解析解,轨道根数把解的参数化结果变成可读的几何量,摄动运动方程描述真实力场下这些几何量如何随时间变化,参考坐标系则规定所有量的度量基准。四者缺一不可,且参照约定必须前后一致。
二体问题:一切轨道计算的基准解
二体问题只考虑两个质点之间的万有引力,忽略其他天体、非球形引力、大气阻力、太阳光压等一切扰动。设中心天体质量为 M,航天器质量为 m,引力常数 G,则引力加速度为 μ/r²,其中引力参数 μ = G(M+m)。对地球轨道,工程上通常取 μ ≈ 3.986004418×10¹⁴ m³/s²,这一数值直接决定周期、速度与半长轴之间的换算关系。
二体问题的核心结论有三条。第一,轨道是圆锥曲线,偏心率 e 决定其类型:e=0 为圆,0<e<1 为椭圆,e=1 为抛物线,e>1 为双曲线;航天器在轨运行对应椭圆或圆轨道。第二,面积速度为常数,等价于角动量守恒,这是开普勒第二定律的动力学本质。第三,轨道周期 T 与半长轴 a 满足 T = 2π√(a³/μ),即开普勒第三定律。以半长轴 a = 7000 km 的近地轨道为例,T ≈ 2π√((7.0×10⁶)³/3.986×10¹⁴) ≈ 5828 s ≈ 97.1 min,这是后续判断漂移速率是否合理的量级基准。
二体解还给出活力公式 v² = μ(2/r − 1/a),它把某点的速度大小与当地地心距、半长轴直接联系起来。轨道保持中判断一次机动是否"够用",往往先用活力公式估算速度增量,再进入摄动修正。
轨道根数:把解析解翻译成几何量
二体解的参数化结果就是六根数:半长轴 a、偏心率 e、轨道倾角 i、升交点赤经 Ω、近地点幅角 ω、真近点角 f(或用平近点角 M、偏近点角 E 替代)。前五个决定轨道在空间中的形状与朝向,第六个决定航天器在轨道上的瞬时位置。对轨道保持而言,a 决定周期与漂移速率,e 决定地心距的周期性起伏,i 与 Ω 决定星下点轨迹的南北与东西走向,ω 决定近地点的方位。地球静止轨道保持中,最关心的通常是 a(东西向漂移)、e(南北向摆动)与 i(倾角漂移)三个量的长期演化。
需要特别强调:轨道根数在 e→0 或 i→0 时会出现奇异,Ω 与 ω 的物理意义退化。工程实现中常改用非奇异根数,例如用 ecos(ω)、esin(ω) 或直角坐标分量代替,避免数值计算中的除零与跳变。这一点在静止轨道保持的数值积分中尤为关键。
摄动运动方程:真实力场下的根数演化
真实航天器受到的力远不止中心引力。地球非球形引力(主要是 J₂ 项)、大气阻力、太阳与月球三体引力、太阳光压等都会使轨道根数缓慢变化。把这些扰动加速度代入二体运动方程,再用常数变易法把根数视为随时间变化的量,就得到摄动运动方程,典型形式为拉格朗日行星方程:每个根数的时间导数由扰动势函数对相应根数的偏导给出。
以 J₂ 摄动为例,它引起升交点赤经的长期漂移率约为 −(3/2)J₂(Re/p)²ncos i,其中 p=a(1−e²),n 为平均角速度。对 i≈0 的静止轨道,cos i≈1,Ω 漂移率约为 −0.013°/天量级,这意味着若不控制,升交点会持续西退。太阳光压与日月引力则主要引起 e 与 i 的周期性及长期变化,是静止轨道南北保持的主要扰动源。理解这些量级,才能判断哪些项必须建模、哪些项可以忽略。
参考坐标系:所有量的度量基准
轨道计算中最容易出错的地方不是公式,而是坐标系约定。常用坐标系包括:地心惯性系(ECI,如 J2000),用于描述轨道根数与动力学积分;地心固连系(ECEF),用于描述星下点、地面站可见性与静止轨道定点位置;轨道坐标系(径向—横向—法向),用于描述相对运动与机动方向;以及卫星本体坐标系,用于姿态与推力方向。切换坐标系时必须明确旋转顺序与时间尺度(UTC、TAI、TT 等),否则会出现方向反号或量级错误。
在轨道保持任务中,一个典型约定是:ECI 下积分得到瞬时根数,再转换到 ECEF 计算星下点经度与纬度,最后用轨道坐标系给出东西/南北机动的方向定义。所有正负号、角度起算点、旋转方向必须在任务文档中固定,并在图示与代码中保持一致。
四者如何协同工作
把上述四件事串起来看:先用二体问题给出基准轨道与周期量级,再用轨道根数把基准轨道参数化,接着用摄动运动方程描述根数在真实力场下的长期与周期变化,最后用参考坐标系把根数、星下点、机动方向统一到同一度量基准上。轨道保持的每一次判断——是否需要机动、机动多少、机动后是否达标——都依赖这条链的完整性。若跳过摄动直接按二体推算,静止轨道的定点经度会在数周内明显漂移;若坐标系约定不统一,东西/南北机动的方向可能完全反号。