轨道动力学程序基本测试方法
1. 引言
轨道动力学程序的正确性是计算轨道特性和验证姿轨控程序的必要条件，而轨道动力学程序也需要严格可信的验证和测试。
在轨道与位置动力学仿真中，“程序能跑”不等于“动力学是对的”。
由于数值积分、坐标系、符号方向、力模型等问题，错误的动力学程序往往可以给出“看起来合理”的轨迹，但其物理意义是错误的。
因此，在引入复杂摄动力、姿态耦合或控制算法之前，必须建立一套基于最基本物理规律的位置动力学验证方案。
本文给出一套工程可落地、结论可信、易于自动化的轨道动力学基本性测试方法。
2. 验证对象与基本假设
2.1 验证对象
- 质点位置动力学仿真程序
- 状态变量为(r,v)(\mathbf{r},\mathbf{v})(r,v)
- 数值积分方式不限（RK4、RKF、ODE113 等）
2.2 基本假设
- 惯性参考系（如 ECI）
- 地心引力，两体问题
- 无摄动力、无推力、无非惯性力
动力学模型为：
r˙=v \dot{\mathbf{r}}=\mathbf{v}r˙=v
v˙=−μrr3 \dot{\mathbf{v}}=-\mu\frac{\mathbf{r}}{r^3}v˙=−μr3r
其中r=∣r∣r=|\mathbf{r}|r=∣r∣，μ=GM\mu=GMμ=GM。
3. Level 0：数值与物理一致性检查
3.1 单位与量纲验证
这是最容易忽略但最致命的一步。
- r\mathbf{r}r：m
- v\mathbf{v}v：m/s
- a\mathbf{a}a：m/s²
- μ\muμ：m³/s²
加速度模型应为：
a=−μrr3 \mathbf{a}=-\mu\frac{\mathbf{r}}{r^3}a=−μr3r
判断条件
- 数值量级合理（LEO 轨道∣a∣≈8∼9m/s2|\mathbf{a}|\approx8\sim9\text{m/s}^2∣a∣≈8∼9m/s2）
- 不出现异常放大或衰减
4. Level 1：圆轨道基准测试
4.1 圆轨道初始条件构造
选择最简单、解析解明确的情况：
∣r0∣=R |\mathbf{r}_0|=R∣r0∣=R
∣v0∣=μ/R |\mathbf{v}_0|=\sqrt{\mu/R}∣v0∣=μ/R
r0⋅v0=0 \mathbf{r}_0\cdot\mathbf{v}_0=0r0⋅v0=0
并令轨道位于 XY 平面。
4.2 位置与速度分量行为
理论解为：
x(t)=Rcos(ωt),y(t)=Rsin(ωt) x(t)=R\cos(\omega t),\quad y(t)=R\sin(\omega t)x(t)=Rcos(ωt),y(t)=Rsin(ωt)
vx(t)=−Rωsin(ωt),vy(t)=Rωcos(ωt) v_x(t)=-R\omega\sin(\omega t),\quad v_y(t)=R\omega\cos(\omega t)vx(t)=−Rωsin(ωt),vy(t)=Rωcos(ωt)
其中：
ω=μ/R3 \omega=\sqrt{\mu/R^3}ω=μ/R3
测试方法
- 将数值结果与正弦函数对比
- 检查幅值、相位、周期稳定性
判断条件
- 无相位漂移
- 周期与理论一致
- 分量连续、光滑
该测试对符号错误、坐标轴混乱极为敏感，是最有效的“低级错误过滤器”。
4.3 模长守恒
对于圆轨道：
∣r(t)∣=R |\mathbf{r}(t)|=R∣r(t)∣=R
∣v(t)∣=μ/R |\mathbf{v}(t)|=\sqrt{\mu/R}∣v(t)∣=μ/R
判断条件
- 相对误差小于给定阈值
- 无随时间累积的漂移
5. Level 2：守恒律验证（物理正确性的核心）
5.1 角动量守恒
角动量定义为：
h=r×v \mathbf{h}=\mathbf{r}\times\mathbf{v}h=r×v
在中心引力场中：
- 大小守恒
- 方向守恒
判断条件
- ∣h(t)∣|\mathbf{h}(t)|∣h(t)∣为常数
- h\mathbf{h}h方向不随时间旋转
5.2 机械能守恒
机械能为：
E=12v2−μr E=\frac12 v^2-\frac{\mu}{r}E=21v2−rμ
判断条件
- E(t)E(t)E(t)为常数
- 允许数值振荡，但不允许单调漂移
能量是对积分器稳定性最敏感的指标，建议作为自动化回归测试项。
5.3 关于线动量的说明
在地心引力场中：
- 单个质点的线动量不守恒
- 因此外力场下不将其作为验证指标
6. Level 3：椭圆轨道一致性验证
6.1 Vis-Viva 方程验证
椭圆轨道满足：
v2=μ(2r−1a) v^2=\mu\left(\frac{2}{r}-\frac{1}{a}\right)v2=μ(r2−a1)
测试方法
- 数值计算r(t)r(t)r(t)与v(t)v(t)v(t)
- 验证是否满足上述关系
判断条件
- 曲线吻合
- 反算半长轴aaa为常数
6.2 拉普拉斯–龙格–楞次向量
定义：
e=v×hμ−rr \mathbf{e}=\frac{\mathbf{v}\times\mathbf{h}}{\mu}-\frac{\mathbf{r}}{r}e=μv×h−rr
其性质：
- 方向指向近地点
- 模长等于偏心率eee
判断条件
- ∣e∣|\mathbf{e}|∣e∣恒定
- e\mathbf{e}e方向不旋转
这是对“是否严格为1/r21/r^21/r2引力”的强验证。
7. Level 4：轨道根数稳定性
由(r,v)(\mathbf{r},\mathbf{v})(r,v)反算轨道根数：
- 半长轴aaa
- 偏心率eee
- 倾角iii
- 升交点赤经Ω\OmegaΩ
- 近地点幅角ω\omegaω
判断条件
- a,e,i,Ω,ωa,e,i,\Omega,\omegaa,e,i,Ω,ω为常数
- 真近点角ν\nuν单调变化
8. Level 5：时间反演一致性（终极数值验证）
8.1 测试流程
- 正向积分时间TTT
- 速度取反v→−v\mathbf{v}\rightarrow-\mathbf{v}v→−v
- 再积分时间TTT
8.2 判断条件
∣rend−r∗0∣<ε |\mathbf{r}_{end}-\mathbf{r}*0|<\varepsilon∣rend−r∗0∣<ε
∣v∗end+v0∣<ε |\mathbf{v}*{end}+\mathbf{v}_0|<\varepsilon∣v∗end+v0∣<ε
该测试可有效暴露积分器非对称误差与隐藏不一致。
9. 推荐的自动化测试执行顺序
- 单位与加速度量级
- 圆轨道正弦分量
- |r|、|v| 守恒
- 角动量守恒
- 机械能守恒
- vis-viva 方程
- LRL 向量守恒
- 轨道根数稳定
- 时间反演
10. 总结
轨道动力学仿真的正确性不能依赖“轨迹看起来对”。
只有从最简单物理事实出发，逐层引入不变量与几何约束，才能获得可信的结论。
一句话概括：
- 圆轨道验证结构正确，
- 守恒律验证物理正确，
- 几何不变量验证模型严格，
- 时间反演验证数值可靠。
