Clohessy-Wiltshire 交会轨道仿真器 返回
航天工程·轨道力学

Clohessy-Wiltshire 交会轨道仿真器

用于设计追踪航天器(追踪器)与目标卫星(目标)进行交会对接的轨道工具。通过改变目标高度、初始相对位置、飞行时间,可从 CW 方程的闭解析解实时计算双脉冲机动的 Δv 预算,帮助选择燃料匹配的接近路线。

参数设置
目标轨道高度
km
目标卫星的圆形轨道高度(距地表)
初始相对位置 X(径向)
km
地心-目标方向。正值表示轨道外侧(高轨道)
初始相对位置 Y(沿迹向)
km
目标速度方向。正值表示前方
初始相对位置 Z(法向)
km
轨道平面外方向。作为独立振动求解
飞行时间
min
从第一脉冲到第二脉冲的飞行时间
机动类型
接近方向选择(参考显示)
计算结果
目标轨道周期(min)
平均运动 n(rad/s)
初始距离(km)
第一脉冲 Δv₁(m/s)
第二脉冲 Δv₂(m/s)
总 Δv(m/s)
LVLH 坐标系 — 追踪器相对轨迹

x = 径向(外向)、y = 沿迹向(速度方向)。黄点●为初始位置,绿点●为目标(原点),红线为接近轨迹,蓝箭头为 Δv 方向。

相对距离 vs 时间
Δv vs 飞行时间(优化曲线)
理论和主要公式

$$\ddot{x} - 3n^2 x - 2n\dot{y} = a_x,\quad \ddot{y} + 2n\dot{x} = a_y,\quad \ddot{z} + n^2 z = a_z$$

Clohessy-Wiltshire 方程。n = 目标轨道平均运动(rad/s)、x = 径向、y = 沿迹向、z = 法向。面内(x,y)耦合,面外(z)独立振动。

$$n = \sqrt{\frac{\mu}{a^3}},\quad a = R_\oplus + h,\quad T = \frac{2\pi}{n}$$

μ = 398600.4418 km³/s²(地球重力常数)、R_⊕ = 6378.137 km、h = 轨道高度。由此确定轨道周期 T 和 n。

$$\Delta v_1 = -\Phi_{rv}^{-1}\,\Phi_{rr}\,\mathbf{r}_0,\qquad \Delta v_{\text{total}} = |\Delta v_1| + |\Delta v_2|$$

状态转移矩阵的闭解析解用于双脉冲交会。第一脉冲 Δv₁ 将初始位置 r₀ 转移到目标位置 r_f=0,第二脉冲 Δv₂ 消去到达时的剩余速度。

Clohessy-Wiltshire 方程 — 航天器交会轨道

🙋
「交会」就是航天船与空间站对接的时候,对吧?那直接加速追上去不就行了吗?为什么要这么复杂的计算?
🎓
这是轨道力学最神奇的地方。在地面上的直觉完全不适用。如果你对前面的目标向前加速,能量增加,你的轨道膨胀,反而会变慢,结果是离目标越来越远。反过来减速会让轨道下降,速度反而加快。所以直觉的「追上去」做法在太空中会失败。这就是为什么需要 Clohessy-Wiltshire 方程——它用数学语言告诉你,在目标卫星的参考框架里怎样运动才能可靠地交会。
🙋
那为什么不直接用轨道计算,非要"线性化"?「线性化」是什么意思?
🎓
因为完整的两体问题求解非常复杂,需要数值计算。但如果追踪器离目标足够近——比如说距离只有几百公里,而轨道半径有六千多公里——那么相对距离可以看作「小扰动」。这时我们用泰勒展开,只保留一阶项,方程变成线性的,就能用矩阵求出闭形式解析解。这就是「线性化」。结果就是 CW 方程:x'' - 3n²x - 2ny' = a_x 这样简洁的形式。在 LVLH 坐标系里——也就是在目标卫星的坐标系里——能给出完全解析的轨迹。
🙋
LVLH 坐标系って何ですか?Local-Vertical/Local-Horizontal ですか?
🎓
そう、LVLH 坐標系(局所竪・局所水平)は目標衛星に張り付いた回転座標系です。x 軸は地心-目標方向(radial、外向きが正)、y 軸は速度方向(along-track、前方が正)、z 軸は軌道面外(cross-track)。この座標系で見ると、地球からの二体力場の影響が「力の勾配」として現れ、それが CW 方程式の 3n² x や 2ny' という項になるわけです。
🙋
「双脉冲交会」是什么?为什么有第一脉冲和第二脉冲?
🎓
最基本的交会策略就是「双脉冲」。第一脉冲时,点燃发动机,从初始位置「打向」目标位置的轨道。然后无推力漂移,让航天器沿着这条轨道飞向目标。到达目标时,相对速度还有剩余,所以第二脉冲再点燃发动机把相对速度消为零。这样就完成对接准备了。CW 方程是线性的,所以状态转移矩阵的逆矩阵 Φ_rv⁻¹ を使うと、第一脉冲 Δv₁ = -Φ_rv⁻¹·Φ_rr·r₀ という公式が出ます。この式は Clohessy と Wiltshire (1960年) が紙と鉛筆で導出したもので、Gemini 計画のランデブー設計で使われた歴史的成果です。
🙋
时间越长,Δv 越小吗?下面的图有很多"谷"……怎么选?
🎓
好眼力!确实 τ = n·t 接近 π、2π、3π... 时,状态转移矩阵的行列式接近零,第一脉冲的大小会发散。这不是计算错误,而是一个物理现象——当传递时间恰好等于轨道周期的半倍数时,起点和终点的几何关系变得「奇异」了。所以不是「时间长 Δv 就小」,而是谷和峰交替出现。实务中一般选择「轨道周期的一半到一个周期」范围内的最低点,比如空间站的交会通常用 1~2 个轨道周期(90~180 分钟),分几个阶段慢慢接近。SpaceX 龙飞船也是这样做的。
🙋
V-bar、R-bar、Natural Drift 这三个模式怎么选啊?
🎓
是接近方向和策略的不同选择。V-bar Hop 是沿着目标速度方向(前后)接近,这是标准做法,姿态稳定好操控。R-bar Hop 是从下方(径向)上升接近,有个优点是万一发动机失灵,重力会自然把你拉开,更安全。国际空间站规定新来的航天器必须用 R-bar 接近。Natural Drift 用的是 CW 自由解中的椭圆轨道(足球形轨道),基本不用推力,让轨道自己转圈改变位相。燃料最省但要花时间。HTV 这样的大型货运飞船通常用 Natural Drift + V-bar 的组合:前期无推力漂移节省燃料,最后阶段用 V-bar 精准靠近。

常见问题

Clohessy-Wiltshire (CW) 方程(1960年)是用于描述追踪航天器在目标卫星周围相对运动的线性化微分方程组。在目标卫星固定的 LVLH 坐标系中(x=径向, y=沿迹向, z=法向),由三个耦合方程 x'' - 3n²x - 2ny' = ax、y'' + 2nx' = ay、z'' + n²z = az 表示。对于短距离(数百公里以内)和近圆形轨道,CW 方程拥有闭解析解,可用于快速设计交会 Δv 预算。联盟号-空间站、SpaceX 龙飞船、JAXA HTV 等所有对接航天器的最终接近阶段都采用 CW 方程作为标准。
设状态转移矩阵 Φ(t) = [Φ_rr Φ_rv; Φ_vr Φ_vv],初始相对位置 r0 在时刻 t_f 要到达目标相对位置 r_f=0,第一脉冲为 Δv1 = -Φ_rv^(-1)·Φ_rr·r0。到达目标时的剩余速度通过第二脉冲 Δv2 消去,总 Δv = |Δv1| + |Δv2| 构成燃料预算。本工具根据 τ = n·t_f 计算矩阵元素,面内方向(x,y)作为二元联立方程求解,面外方向(z)独立求解,最后返回 Δv 值。
V-bar Hop 是沿着速度方向(y 轴,沿迹向)的接近方式,R-bar Hop 是从径向方向(x 轴,地心-卫星方向)上升的接近方式。V-bar 是标准最终接近方式,便于姿态稳定和相对速度控制。而对于国际空间站这样的大型目标,为了安全考虑常采用 R-bar 接近(以避免排气羽流冲击船体)。Natural Drift 是使用轨道自由运动消除位相差的无推力模式,燃料最经济但耗时较长。
CW 方程适用于以下条件:(1) 目标轨道近似圆形(离心率 e<0.01);(2) 相对距离远小于轨道半径(通常数百公里以内);(3) 飞行时间较短(数周期以内);(4) J2 摄动和大气阻力可忽略。对于离心率较大的轨道需采用 Tschauner-Hempel 方程;对于长距离和长时间需考虑 Hill-Clohessy-Wiltshire 高阶项或数值积分;低轨道受空气阻力影响时需加入 LEO 修正项。本工具专为 200~2000 km 高度范围的低地球轨道标准交会设计。

实际应用

国际空间站补给任务:SpaceX 龙飞船、诺格公司的天鹅座、JAXA 的 HTV/HTV-X、俄罗斯进步号、ESA ATV 等所有补给飞船都采用基于 CW 方程的最终接近阶段。地面计算出轨道高度和位相差后,用 Δv 预算来规划接近路线,在多个保持点(距离几公里处)逐步降低相对速度和距离,最终以 200m / 30m / 10m 三个阶段靠近。本工具计算的 Δv 预算与实际任务规划所用公式完全相同。

载人龙飞船和联盟号对接:NASA 的载人龙飞船和俄罗斯联盟号都有自动交会系统(龙飞船用 DragonEye LIDAR,联盟号用 KURS 微波系统)可执行 CW 解算轨迹。但最终 200 米以内航天员可接管控制,采用「手动模式」。地面预测的 CW 轨迹与实时传感器数据差异超过许可范围时,自动系统进入 Hold 状态。设计时会在最终 Δv 预算基础上增加 20~30% 的裕度来应对导航误差和应急中止。

卫星维修和空间碎片清除:对于故障或无合作的目标卫星(燃料耗尽、控制失效)的维修和碎片清除任务,地面可通过雷达测量轨道参数,然后用 CW 方程计划接近路线。Northrop Grumman 的 MEV(任务延期飞行器)已成功为商业通信卫星 Intelsat 进行了在轨加注和轨道抬升,其轨道计算采用的是 CW 方程的扩展形式(Yamanaka-Ankersen 解)以适应离心轨道。

大学和小型卫星编队飞行:NASA MMS(四星编队)、CanX-4/5、PRISMA 等多卫星编队保持相对距离和队形的运行,都基于 CW 方程的自由解(齐次解)。这种「足球形轨道」(football orbit)可以周期性地描绘出椭圆,使用几乎零推力就能在编队内维持相对位置。Natural Drift 模式正是这一原理的实现。

常见误区和注意事项

最大的陷阱是忽视「目标轨道必须是圆形」这个前提。CW 方程在推导时假设离心率 e≈0,对于国际空间站和太阳同步卫星这类 e<0.001 的完美圆轨道很精确,但对于椭圆轨道的目标(如轨道转移中的上面级或周期彗星)就不适用。当 e>0.01 时误差明显,e>0.1 时 Δv 预测可能偏离实际 10~50%。椭圆轨道交会要用 Tschauner-Hempel 方程或 Yamanaka-Ankersen 解析解,或者用数值轨道积分。本工具仅适用低地球轨道的圆形轨道。

第二个误解是「飞行时间越长 Δv 越低」。下面的「Δv vs 飞行时间」曲线显示,当 τ = n·t 接近 π 的整数倍时,状态转移矩阵行列式趋于零,第一脉冲所需 Δv 会爆炸。这不是工具的 bug,而是物理上的真实现象——轨道初值和终值的配置在某些时间点变成奇异的。航天器设计早期就要规避这些「禁忌时间」,这就是为什么历史上 Gemini 和 Apollo 项目有详细的「ΔV vs 时间」查表。实际选择通常是局部最小值,并留出安全裕度。

最后,切勿将 Δv 等同于燃料消耗。本工具计算的是理想瞬时脉冲(impulsive maneuver)的数值,实际飞行中需要加上:(1) 有限燃烧时间的重力损失(gravity loss);(2) 姿态调整用的 RCS 推进剂独立消耗;(3) 导航误差和应急中止用的裕度(20~30%);(4) 最终接近阶段(相对速度<0.3 m/s)的精细控制。还要通过 Tsiolkovsky 方程 Δm/m = 1 - exp(-Δv/(Isp·g)) 换算实际质量——同样 10 m/s 的 Δv,在联盟号的 hypergolic 推进剂(比冲 280 s)下比 SpaceX Dragon 高效推进剂(比冲 340+ s)消耗更多燃料。

使用指南

  1. 在 300~2000 km 范围内设置目标轨道高度。国际空间站约 408 km,静止轨道约 35786 km。
  2. 输入初始相对位置(X、Y、Z 分量),定义追踪器距目标多远的起始状态。
  3. 指定所需交会时间(分钟),工具即实时计算双脉冲机动的第一、第二 Δv。
  4. 验证总 Δv 预算与航天器搭载推进剂的匹配。通常追踪航天器携带 50~200 m/s 的推进剂。

具体计算示例

对于空间站轨道(高度 408 km)的交会,初始相对位置 X=50 km、Y=10 km、Z=5 km,飞行时间 120 分钟:目标轨道周期为 92.68 分钟,平均运动 n=0.001107 rad/s。第一脉冲 Δv₁=8.3 m/s,第二脉冲 Δv₂=7.9 m/s,总 Δv=16.2 m/s。而对于静止轨道(35786 km)的同样相对距离,n=7.29×10⁻⁵ rad/s,第一脉冲仅需 Δv₁≈0.58 m/s。

实务操作要点