2. 轨道动力学基础与评估关联

轨控任务评估的核心问题可以概括为:在给定控制量、给定时间窗口和给定测量条件下,航天器能否到达目标状态,以及到达的精度和代价是否可接受。回答这个问题,必须先把轨道动力学模型与评估指标之间的映射关系理清楚。本节从二体问题出发,逐层引入摄动、轨道根数、状态递推与误差传播,说明每类模型在评估中的使用条件与局限。

一、二体问题:评估的基准模型

二体问题假设航天器只受中心天体引力作用,运动方程为

r̈ = -μ · r / |r|³

其中 μ 为中心天体引力常数,r 为地心惯性系下的位置矢量。该模型给出开普勒轨道:轨道为圆锥曲线,机械能、角动量和轨道面法向守恒,轨道根数在无控条件下保持不变。在轨控评估中,二体模型主要承担三类任务:一是作为控制量初值设计的基准,例如用高斯摄动方程或兰伯特问题求解脉冲速度增量;二是作为快速评估的解析工具,用开普勒方程由平近点角反解偏近点角,得到任意时刻的位置速度;三是作为误差传播的参考轨道,把摄动和测量误差视为叠加在二体解上的小量。

二体模型的局限同样明确。它无法描述地球非球形引力、大气阻力、太阳光压和第三体引力引起的轨道根数长期漂移,因此不能直接用于评估长时间在轨保持任务的燃料消耗。若评估周期只有几个轨道周期、且任务对位置精度要求不高,二体模型可以作为一阶近似;一旦评估周期超过一天或精度要求进入百米量级,就必须引入摄动模型。

二、摄动模型:从定性到定量的评估依据

摄动模型在二体方程右端叠加扰动加速度。工程评估中最常用的是地球非球形引力中的 J2 项,其加速度量级约为 10⁻³ m/s² 量级,会引起升交点赤经和近地点幅角的长期变化。对于近地轨道,大气阻力在 300~600 km 高度范围内不可忽略,其加速度与面质比和大气密度直接相关,是轨道衰减评估的主要输入。太阳光压对高面质比航天器影响显著,第三体引力则在深空或高轨任务中占主导。

在评估中使用摄动模型,需要明确三个条件。第一,摄动项的选择要与任务轨道和评估周期匹配:近地短周期评估可只保留 J2,长期轨道寿命评估必须加入大气阻力,高轨任务则要考虑日月引力。第二,摄动模型的参数本身带有不确定性,例如大气密度随太阳活动变化可达数倍,评估时应给出参数区间而非单点值。第三,摄动模型通常以数值积分实现,积分步长和积分器阶数会影响结果,评估报告中应记录积分设置,避免不同工具之间结果不可比。

三、轨道根数:评估指标的载体

轨道根数把位置速度映射为六个具有明确物理意义的参数:半长轴 a、偏心率 e、倾角 i、升交点赤经 Ω、近地点幅角 ω 和平近点角 M。轨控任务评估中,不同指标对应不同根数分量。半长轴决定轨道周期和地面轨迹重复特性,是相位保持评估的核心;偏心率影响近地点高度和地面覆盖均匀性;倾角和升交点赤经决定轨道面朝向,是太阳同步轨道和回归轨道评估的关键;近地点幅角在冻结轨道设计中用于抑制偏心率漂移。

使用轨道根数评估时要注意两点局限。一是根数在 e 接近零或 i 接近零时存在奇异,此时应改用春分点根数或位置速度矢量描述。二是根数变化量与控制量之间不是线性关系,大控制量下必须用非线性映射或数值打靶,不能简单用线性敏感矩阵外推。

四、状态递推与误差传播:评估精度的来源

状态递推是把初始状态和控制量沿时间推进,得到评估时刻的状态。数值递推用积分器直接求解含摄动的运动方程,精度高但计算量大;解析递推用平均根数法或小参数展开,速度快但适用范围受限。评估中常采用“数值递推为主、解析递推做快速筛选”的组合方式。

误差传播回答的是“评估结果有多可信”。初始状态误差、控制执行误差和模型误差会沿轨道动力学放大或缩小。工程上用状态转移矩阵 Φ(t, t₀) 描述一阶传播:δX(t) = Φ(t, t₀) · δX(t₀)。对于二体问题,Φ 有解析形式;含摄动时通常用数值方法求变分方程。评估中应重点关注沿迹方向误差,因为它随时间增长最快,直接决定相位保持和交会窗口的可用性。

轨道动力学模型在轨控任务评估中的分层使用与误差传播 二体问题 r̈ = -μr/|r|³ 开普勒轨道 · 根数守恒 控制量初值 · 快速评估 摄动模型 J2 非球形引力 · 大气阻力 太阳光压 · 第三体引力 数值积分 · 参数区间评估 轨道根数 a · e · i · Ω · ω · M 指标载体 · 物理意义明确 奇异点需换用非奇异根数 状态递推 数值/解析 积分器设置 步长与阶数 误差传播链:δX(t) = Φ(t, t₀) · δX(t₀) 初始状态误差 位置/速度偏差 控制执行误差 推力大小/方向/时长 模型误差 摄动参数不确定性 状态转移矩阵 Φ(t, t₀) 评估时刻 状态不确定度 评估判断依据 沿迹方向误差增长最快 → 决定相位保持与交会窗口可用性 模型选择与评估周期、精度要求、轨道类型匹配 → 否则结果不可信 二体局限 无摄动长期漂移 长周期评估失效 仅作一阶近似 根数局限 e≈0 或 i≈0 奇异 大控制量非线性 需数值打靶验证

五、模型选择与评估结论的对应关系

把上述模型放到评估流程中,可以得到一条清晰的判断链。若评估周期小于一个轨道周期、精度要求为公里级,二体模型加解析递推即可给出可信结论。若评估周期为一天到数天、精度要求为百米级,应使用含 J2 的数值递推,并用状态转移矩阵给出误差包络。若评估周期超过一周或涉及轨道寿命,必须加入大气阻力并考虑太阳活动参数区间,此时评估结论应以概率形式给出,例如“在给定参数区间内,半长轴衰减量落在 1.2~2.8 km 的概率不低于 90%”。

需要特别强调的是,误差传播不是评估的附属步骤,而是评估结论的一部分。同一组控制量,在忽略误差传播时可能显示“满足要求”,在加入沿迹方向误差增长后可能显示“窗口边缘不可用”。因此,规范的轨控任务评估报告应同时给出标称结果和误差包络,并说明所用动力学模型的适用条件与已知局限。

使用轨道根数描述状态时,务必先声明参考系和角度定义。本文所有根数均定义在地心惯性系下,升交点赤经从春分点起算,近地点幅角从升交点起算。若切换到地固系或采用不同的角度起算点,必须重新说明,否则评估结果会出现系统性偏差。
在快速评估阶段,可以先用二体模型加解析递推筛选出可行控制方案,再用含 J2 的数值递推对候选方案做精确校核。这种“粗筛—精算”的组合能在保证评估可信度的同时显著降低计算量。