2. 轨道力学基础与约束分析接口

轨道任务约束分析并不是从零开始“造轨道”,而是把轨道力学已经算出的状态,翻译成任务规划、载荷使用、测控、能源、热控等分系统能够直接判断“能不能做、什么时候做、做多久”的输入条件。因此本章的关键不是重复推导二体问题,而是建立一条清晰的接口链:二体模型给出运动规律 → 轨道六根数给出可读的轨道描述 → 星历计算给出任意时刻的位置速度 → 约束分析把位置速度映射为可见性、光照、姿态、通信等判据。理解这条链,才能知道约束分析里每一个“可用窗口”究竟从哪里来。

一、二体问题:约束分析的时间基准与几何基准

二体问题假设航天器只受中心天体引力作用,忽略摄动、大气阻力、太阳光压和第三体引力。它的价值在于给出解析可算、物理意义明确的运动:轨道是圆锥曲线,航天器在轨道上的位置由开普勒方程联系。对约束分析而言,二体模型至少提供三样东西:

  • 时间基准:由平近点角 M 与偏近点角 E 的关系 M = E − e·sinE,把均匀流逝的时间映射为轨道上的真实位置;
  • 几何基准:由轨道半长轴 a、偏心率 e 确定轨道形状,由倾角 i、升交点赤经 Ω、近地点幅角 ω 确定轨道在惯性空间中的朝向;
  • 周期与角速度:由 a 得到轨道周期 T = 2π√(a³/μ),进而得到平均角速度 n = √(μ/a³),这是所有窗口重复性判断的基础。

需要强调的是,约束分析中常用的“轨道周期”“回归周期”“地面轨迹重复周期”都建立在二体近似之上。当任务精度要求较高时,必须明确二体结果只是初值,后续要用摄动模型修正,但接口形式不变:仍然是“给定时刻 → 给出位置速度”。

二、轨道六根数:从数学解到任务可读参数

轨道六根数 (a, e, i, Ω, ω, M₀) 或 (a, e, i, Ω, ω, ν₀) 是二体解的完整参数化。对约束分析来说,六根数不是六个孤立数字,而是六类约束条件的直接来源:

根数物理含义对约束分析的直接输入
a轨道半长轴周期、速度量级、回归特性、覆盖重复性
e偏心率高度变化范围、近远地点差异、测控与载荷几何变化
i轨道倾角可覆盖纬度范围、太阳光照角变化范围
Ω升交点赤经轨道面与太阳方向的相对关系、光照条件
ω近地点幅角近地点在轨道面内的位置、过顶时刻分布
M₀ 或 ν₀历元时刻的相位具体过顶时刻、窗口起点、星下点轨迹相位

例如,太阳同步轨道通过 i 与 Ω 的特定组合,使轨道面与太阳方向保持近似固定夹角,从而让载荷在相近的地方时过境;回归轨道通过 a 与 i 的配合,使星下点轨迹在若干圈后重复。这些都不是约束分析“额外加上的条件”,而是六根数本身已经决定的几何关系。

三、星历计算:把根数变成逐时刻的位置速度

约束分析真正消费的不是根数,而是星历。星历计算的标准流程是:由历元时刻 t₀ 和 M₀ 计算 t 时刻的平近点角 M = M₀ + n(t − t₀);用开普勒方程迭代求解 E;再由 E 得到真近点角 ν 和地心距 r = a(1 − e·cosE);最后在近焦点坐标系中写出位置,再经三次旋转(ω、i、Ω)转到地心惯性系。速度可由 vis-viva 方程和飞行路径角得到,也可由位置对时间求导得到。

在工程实现中,通常不会每次重新推导,而是调用星历模型或两行根数(TLE)配合 SGP4 等传播器。但无论用解析法还是数值法,约束分析拿到的接口都是一致的:时间序列上的 r(t) 和 v(t)。下面这张图给出从二体模型到约束判据的完整接口链,重点标出每一步的输出如何被下一步消费。

二体模型 → 轨道六根数 → 星历 → 任务约束判据 接口链 二体问题 μ、a、e 决定运动 M = E − e·sinE T = 2π√(a³/μ) n = √(μ/a³) 轨道六根数 a, e, i, Ω, ω, M₀ 形状 + 朝向 + 相位 历元 t₀ 与参考系约定 约束条件的参数来源 星历计算 M → E → ν → r 三次旋转到惯性系 输出 r(t)、v(t) 序列 TLE/SGP4 或数值传播 约束判据 可见性 光照 测控 姿态/能源 约束分析从星历中提取的关键量 几何可见性 r_sat 与地面站 r_gs 仰角 ≥ 5° 判据 光照条件 太阳矢量与地影 β 角、进出影时刻 姿态与指向 v 方向、天底/天顶 目标矢量夹角 能源与热控 光照时长、阴影时长 太阳入射角变化 接口约定:时间系统(UTC/TAI/TT)、参考系(J2000/ECEF)、角度单位与正负号方向必须统一 二体结果作为初值,摄动修正后仍以 r(t)、v(t) 形式进入约束判据 窗口是否可用 = 几何可见 + 光照满足 + 姿态可达 + 能源/热控不越限

四、从力学量到约束输入:四类典型映射

第一类是几何可见性。给定地面站位置 r_gs 和航天器位置 r_sat,可计算站星矢量与当地水平面的夹角,即仰角。当仰角大于等于最小仰角(例如 5°)时,认为可见。可见窗口的起止时刻就是仰角穿越阈值的时刻,通常用星历序列逐点判断再插值细化。

第二类是光照与地影。太阳方向可由太阳星历给出,判断航天器是否处于地球本影或半影,常用圆柱地影模型或圆锥地影模型。输出是进出影时刻和持续时长,直接决定能源分系统的供电与蓄电池放电安排。太阳矢量与轨道面法线的夹角 β 角,则决定太阳翼受照几何和热控外热流。

第三类是姿态与指向。由 r 和 v 可以构造轨道坐标系,得到天底、天顶、飞行方向等参考方向;再结合目标点位置,计算目标矢量与载荷视轴的夹角,判断是否满足指向和侧摆能力。这里 v(t) 的作用常被低估:没有速度方向,就无法确定偏航和俯仰的参考基准。

第四类是能源与热控。光照时长、阴影时长、太阳入射角随时间的变化,都是 r(t) 与太阳方向的函数。约束分析通常把它们整理成时间轴上的“可用/不可用”区段,再与测控窗口、载荷任务窗口求交集。

五、接口实现中的关键约定与检查点

实际工程中最容易出问题的不是公式,而是约定不一致。建议在接口层面固定以下检查点:

  • 时间系统:明确历元是 UTC、TAI 还是 TT,传播步长与插值方式是否一致;
  • 参考系:位置速度在 J2000 惯性系还是 ECEF 地固系,地面站坐标是否已做岁差章动和极移修正;
  • 角度约定:i、Ω、ω 的参考方向和正负号,仰角、方位角的起算基准;
  • 单位:km 与 m、度与弧度、秒与分钟必须统一,避免量级错误;
  • 精度与步长:低轨任务星历步长常取 10~60 s,窗口边界再用二分或插值细化到秒级。
教学示例:设某低轨任务 a = 6878 km、e = 0.001、i = 97.8°、Ω = 120°、ω = 90°、M₀ = 0°,μ = 398600.4418 km³/s²。可算得 T ≈ 5677 s ≈ 94.6 min,n ≈ 0.001107 rad/s。若地面站最小仰角取 5°,用 30 s 步长生成一轨星历,即可得到若干段可见窗口;这些窗口再与地影区段求交,就得到“可见且光照”的载荷可用时段。上述数值为教学假设,实际任务需用精确摄动星历复核。

把这条接口链固定下来后,后续的约束分析就不必反复回到轨道力学细节,而是面向统一的时间序列输入做判断。二体问题提供规律,六根数提供参数,星历提供逐时刻状态,约束分析负责把这些状态翻译成任务可执行的窗口与限制条件。