Exact Nonlinear Attitude Dynamics of Gravity-Gradient Spacecraft: Elliptic Libration Orbits, Separatrices and Energy-Based Stability Thresholds

技术摘要：重力梯度空间载具的精确非线性姿态动力学

问题陈述
重力梯度稳定是地球指向型空间载具的一种基础被动姿态控制概念。然而，对大角度动力学的评估往往依赖于小角度近似或低阶多项式模型（例如扭矩的三次或五次截断）。本文指出，这些局部模型无法保留全局相空间结构，导致在空间载具发生显著姿态偏移时，对稳定性边界、摆动周期及逃逸阈值的预测出现偏差。具体而言，多项式近似可能会引入人工平衡点或消除必要的屏障，从而导致将有界运动误判为旋转运动，或反之亦然。

研究方法
本研究聚焦于在圆轨道运行的刚体空间载具，并将分析限制在绕轨道平面法轴的平面转动。研究方法如下：

1. 精确公式化： 从重力梯度扭矩 $\tau_n = -\frac{3}{2}n^2(I_t - I_r)\sin(2\theta)$ 出发，推导出平面运动方程。通过变量代换 $\phi = 2\theta$，将该非线性微分方程转化为可积的单摆形式。
2. 解析解： 利用雅可比椭圆函数求解该系统。解由无量纲能量参数 $q$ 参数化，该参数将动力学划分为三个区域：
   • 摆动区 ($q < 1$)： 在天底方向附近的有界振荡。
   • 分界线 ($q = 1$)： 连接不稳定平衡点的异宿轨道。
   • 旋转区 ($q > 1$)： 相对于局部垂线的连续旋转。
3. 对比分析： 将精确的周期性扭矩与三次（Duffing）及五次泰勒级数截断进行对比，以量化拓扑结构的差异。
4. 数值验证： 针对一个典型的低地球轨道案例（高度 500 km）进行仿真。通过高精度八阶龙格-库塔法（DOP853）对解析椭圆解进行验证，验证时长覆盖五个大振幅周期。

主要贡献与结果

• 精确解析解： 本文利用雅可比椭圆函数 ($\text{sn}, \text{cn}, \text{dn}$) 提供了任意摆动初始条件的闭式表达式。这得到了一个与振幅相关的周期 $T(k) = 4K(k)/\omega_g$，其中 $K(k)$ 是第一类完全椭圆积分。
• 周期发散： 分析证实，当振幅趋近于分界线（$90^\circ$）时，摆动周期呈对数发散。这一特征在精确模型中得到了捕捉，但在近似模型中却无法体现。
• 基于能量的稳定性阈值： 推导出了一个紧凑且精确的有界运动判据：$\dot{\theta}^2 + \omega_g^2 \sin^2 \theta < \omega_g^2$。这定义了 $(\theta, \dot{\theta})$ 平面上的“安全包络线”。此外，本文还推导了在固定姿态下，将状态从旋转区域带回摆动区域所需的最小角冲量 ($J_{\min}$)。
• 多项式近似的失效：
  • 三次模型： 三次 Duffing 近似预测了一个虚假的分界线，位于 $\approx 70.17^\circ$ ($\sqrt{3/2}$ rad) 处，而真实的重力梯度屏障位于 $90^\circ$。这低估了稳定摆动区域，误差接近 $20^\circ$。
  • 五次模型： 五次截断无法产生任何非零实数扭矩平衡点，实际上消除了有限角度的分界线。
• 数值精度： 精确椭圆解与高精度数值积分的匹配度极高，在五个大振幅周期内的最大绝对误差约为 $2.9 \times 10^{-10}$ rad，证实了解析推导的有效性。

意义与主张
本文主张，虽然局部多项式近似在小扰动分析中仍然有效，但在全局稳定性评估、捕获区域定义及大角度恢复裕度方面是不可靠的。所提供的精确解具有以下作用：

1. 透明基准： 为验证非线性传播算法和控制设计提供严谨的标准。
2. 运行阈值： 提供明确的公式，用于判断扰动是否会导致脱离摆动势阱，以及防止脱离所需的最小控制量。
3. 建模警示： 展示了增加泰勒级数阶数并不保证能保留全局相空间拓扑结构。

作者强调，对于可能发生大幅度姿态偏移的空间载具，在评估分界线和逃逸阈值时，必须保留精确的周期性重力梯度扭矩。所推导的保守解提供了一个精确的“骨架”，可在未来工作中进一步扩展，以包含轨道偏心率、阻尼或柔性附件等扰动。