📚 课程目录第3章摄动概念与分类

3. 摄动概念与分类

二体问题给出的开普勒轨道是一条理想化曲线：在中心引力场中，航天器沿固定椭圆（或圆、抛物线、双曲线）运动，六个轨道根数全部为常数。但真实航天器并不处在这种理想环境中。地球并非质量均匀的球体，其赤道隆起产生的非球形引力、大气阻力、太阳与月球的第三体引力、太阳光压等都会叠加到中心引力之上，使实际运动偏离开普勒轨道。这种"偏离"就是摄动，而描述偏离的数学工具，就是摄动加速度与摄动运动方程。

一、摄动加速度的定义

设航天器质量为 m，位置矢量为 r。中心天体（如地球）的引力加速度为 a₀ = −μr/r³，其中 μ = GM 为地球引力常数（约 3.986×10¹⁴ m³/s²）。除中心引力之外，作用在航天器上的所有其他加速度之和，统称为摄动加速度：

a_pert = a_total − a_0
       = a_NS + a_3rd + a_drag + a_SRP + a_tide + ...

式中各项分别对应地球非球形引力（主要是 J₂ 项）、第三体引力、大气阻力、太阳光压和固体潮等。摄动加速度的典型量级差异极大：在近地轨道（约 500 km 高度），J₂ 摄动加速度约为 10⁻² m/s² 量级，而大气阻力加速度约为 10⁻⁷～10⁻⁶ m/s²；在 36000 km 的地球同步轨道，J₂ 已衰减到 10⁻⁵ m/s² 量级，太阳光压则升至 10⁻⁷ m/s² 量级。量级差异决定了不同轨道高度上主导摄动源的不同。

二、摄动运动方程

将摄动加速度代入牛顿第二定律，得到摄动运动方程：

r̈ = −μ r / r³ + a_pert(r, ṙ, t)

这是一个二阶非线性非齐次微分方程。直接对位置矢量积分虽然可行，但工程上更关心轨道根数的变化——因为轨道根数直接对应轨道的形状、大小和空间取向。为此，需要把摄动加速度投影到轨道坐标系中，建立轨道根数形式的摄动方程（又称拉格朗日行星方程）。以半长轴 a、偏心率 e、倾角 i、升交点赤经 Ω、近地点幅角 ω、真近点角 f 为例，其一般形式为：

da/dt = (2a²/h)[e sin f · R + (p/r) · S]
de/dt = (1/h)[p sin f · R + ((p+r)cos f + r e) · S]
di/dt = (r cos(ω+f)/h) · W
dΩ/dt = (r sin(ω+f)/(h sin i)) · W
dω/dt = (1/(he))[−p cos f · R + (p+r) sin f · S] − (r sin(ω+f) cos i/(h sin i)) · W

其中 R、S、W 分别是摄动加速度在径向、横向（轨道面内垂直于径向）和法向（轨道面法线方向）的分量，h 为角动量大小，p = a(1−e²) 为半通径。这组方程的意义在于：它把"位置偏离"翻译成了"轨道根数如何随时间变化"，使摄动分析从几何问题转化为参数演化问题。

三、保守力摄动与非保守力摄动

按力学性质，摄动可分为保守力摄动和非保守力摄动。保守力摄动存在势函数，做功与路径无关，航天器的机械能（动能+势能）在摄动过程中守恒。地球非球形引力、第三体引力、固体潮都属于保守力摄动。非保守力摄动不存在势函数，做功与路径有关，会持续改变航天器的总机械能。大气阻力和太阳光压是典型的非保守力摄动。

这一区分在工程上极为关键。对于保守力摄动，轨道半长轴 a 没有长期变化，只有周期性振荡，轨道能量整体保持稳定；对于非保守力摄动，大气阻力会持续消耗轨道能量，使 a 单调减小，轨道逐渐降低，最终再入大气层。太阳光压虽然量级小，但对高面质比航天器（如太阳帆、大型柔性太阳翼）会产生可观的累积效应，且其方向随航天器姿态和日地几何关系变化，处理时需结合姿态模型。

判断技巧：若摄动力可写成某个标量函数 U 的梯度（a = ∇U），则为保守力；若摄动力方向与速度方向相关且做负功（如阻力），或与光照条件相关（如光压），则为非保守力。保守力摄动不改变轨道平均能量，非保守力摄动会改变。

四、三类轨道根数变化的物理含义

摄动对轨道根数的影响，按时间尺度分为三类：

长期变化（Secular）：随时间线性累积，与轨道周期无关。最典型的是 J₂ 摄动引起的升交点赤经 Ω 退行和近地点幅角 ω 旋进。对于太阳同步轨道，正是利用 Ω 的长期退行速率与地球公转速率匹配（约 0.9856°/天），使轨道面与太阳保持固定夹角。长期变化是轨道寿命分析和星座构型维持的核心关注对象。
长周期变化（Long-period）：周期远大于轨道周期，振幅有限，不随时间无限累积。例如 J₂ 摄动下偏心率 e 和倾角 i 的缓慢振荡，其周期与 ω 或 Ω 的进动周期相关。长周期变化在数天到数月尺度上影响轨道形状，需要在长期预报中加以考虑。
短周期变化（Short-period）：周期等于或接近轨道周期，在轨道周期内均值为零，不产生累积效应。例如半长轴 a 在一个轨道周期内的周期性抖动。短周期变化主要影响瞬时轨道精度，对长期演化趋势影响有限，但在高精度定轨和碰撞预警中不可忽略。

五、教学示例：J₂ 摄动下 Ω 的长期退行

以一颗近地太阳同步轨道卫星为例，设轨道高度 h = 700 km，地球半径 Rₑ = 6378 km，则半长轴 a = 7078 km。取 J₂ = 1.08263×10⁻³，轨道倾角 i = 98.2°（逆行轨道）。升交点赤经的长期变化率近似为：

dΩ/dt ≈ −(3/2) · J₂ · (Rₑ/p)² · n · cos i

其中：
  p = a(1−e²) ≈ a = 7078 km（取 e ≈ 0）
  n = √(μ/a³) = √(3.986×10¹⁴ / 7078³) ≈ 1.048×10⁻³ rad/s
  cos i = cos(98.2°) ≈ −0.1426

代入：
  dΩ/dt ≈ −1.5 × 1.08263×10⁻³ × (6378/7078)² × 1.048×10⁻³ × (−0.1426)
         ≈ +1.99×10⁻⁷ rad/s
         ≈ +0.985°/天

计算结果约为 +0.985°/天，与地球绕太阳公转的平均角速率（360°/365.25 天 ≈ 0.9856°/天）几乎一致，这正是太阳同步轨道的设计条件。该示例说明：通过选择轨道高度和倾角，可以主动利用 J₂ 摄动的长期效应实现特定的轨道面进动，而不是被动地"承受"摄动。这也体现了摄动分析从"误差来源"到"设计工具"的转变。

注意：上述计算为近似估算，忽略了 e 的高阶项、大气阻力和其他摄动源的影响。实际工程中需用数值积分或高精度解析模型进行验证。示例参数为教学假设，不代表任何具体任务的实际参数。

六、摄动分类的工程意义

理解摄动分类，最终要落到工程决策上。保守力摄动（J₂、第三体、潮汐）主要影响轨道面的空间取向和形状的周期性变化，通常可以通过轨道设计加以利用或补偿；非保守力摄动（大气阻力、光压）则持续改变轨道能量，决定航天器的轨道寿命和燃料消耗。在轨道维持策略中，对保守力摄动引起的长期漂移，可通过定期轨道机动修正；对非保守力摄动引起的能量衰减，则需要根据阻力系数和大气密度模型预估衰减速率，提前规划抬轨机动。三类时间尺度变化的区分，则决定了轨道预报的时间步长和模型精度要求：短周期项影响瞬时精度，长周期项影响中期预报，长期项决定长期演化趋势。