【Demllie航天】二体问题（是平面问题）

前言

我打算做个二维的坎巴拉，其实就是电脑版的航天模拟器，如果你也有这样的想法，那你应该看看这篇文章。我们简单点，只讨论二体问题，多了也不会。
时间系统、参考系、二体问题、三体问题的一些基本概念。是第三章航天器轨道与轨道动力学的两个小节的内容。

---

二体问题

受其他物体摄动影响的两个物体相对运动的一般方程，很容易就可以简化成仅有两个物体的运动方程，即二体问题。首先做两个假设：

1. 物体为球对称的，就可以把物体看作质量集中在其中心
2. 除了沿两物体中心连线作用的引力外，没有其他外力和内力作用

确定一个惯性参考系（无加速度、无转动），其中有两个质量分别是$M$和$m$的物体，它们在惯性系中的位置矢量分别是$r_1、r_2$，两物体相距$r=r_1-r_2$

根据牛顿万有引力定律，对$m、M$有
$$m\ddot{\overrightarrow{r_2}}=-\frac{GMm}{r^3}\vec{r}$$$$M\ddot{\overrightarrow{r_1}}=\frac{GMm}{r^3}\vec{r}$$
两式相减得到二体问题相对运动方程
$$\ddot{\vec{r}}=\ddot{{\vec{r}}_2}-\ddot{{\vec{r}}_1}=-\frac{G(M+m)}{r^3}\vec{r}$$
当两个星球中的一个质量大到$r_1\to 0$，那么其会成为中心天体，有中心天体运动方程
$$\ddot{\vec{r}}=-\frac{GM}{r^3}\vec{r}$$
这说明二体问题与中心天体问题是一致的，为了方便和一般性，定义$\mu =GM$，二体运动方程为$$\ddot{\vec{r}}+\frac{\mu }{r^3}\vec{r}=0（1）$$

对于不同的中心天体，有不同的$\mu$值。对于地球，$\mu =3.986005e5km/s$；对于月球，$\mu =4.902794e3km/s$；对于太阳，$\mu =1.32712438e11km/s$

二体轨道运动常数$h$

1机械能守恒

用$\dot{\vec{r}}$与（1）做点乘，
$$\dot{\vec{r}}\cdot \ddot{\vec{r}}+\frac{\mu }{r^3}\dot{\vec{r}}\cdot \vec{r}=\vec{v}\cdot \dot{\vec{v}}+\frac{\mu }{r^3}\dot{\vec{r}}\cdot \vec{r}=0$$
上式因为结果为零，两项都是标量。于是右边的$\vec{r}$就可以约掉下面一个$r$
$$\vec{v}\cdot \dot{\vec{v}}+\frac{\mu }{r^2}\dot{\vec{r}}=\frac{d}{dt}(\frac{v^2}{2})+\frac{d}{dt}(-\frac{\mu }{r})=0$$
积分后有
$$E=\frac{v^2}{2}+(c-\frac{\mu }{r})（2）$$
E是比机械能（把质量约掉的结果）。它等于从空间一点到另一点克服引力所做的功，常数c依赖于零势能点的选取。如果将中心引力体作为零势能点（例如地球表面），$c=\frac{\mu }{r}$，$r$是星球半径；如果将无穷远处作为零势能点，$c=0$,此时势能始终是负的。
理想情况下，比机械能在卫星沿轨道运行时一直是不变的。

2角动量守恒

用$\vec{r}$与（1）做叉乘，
$$\vec{r}\times \ddot{\vec{r}}+\frac{\mu }{r^3}\vec{r}\times \vec{r}=\vec{r}\times \ddot{\vec{r}}=0$$
所以矢量 $\vec{r}\times \vec{v}$一定是个常数，记为$h$。显然其垂直于$r,v$形成的平面。由此证明航天器的运动被限制在一个固定平面内，称为轨道平面。
观察$r,v$及他们与中心引力体的关系可知，无论航天器处于空间上面位置，当地垂线与矢量$r$的方向一致的，所以当地水平面必定垂直$r$。
定义速度矢量方向与当地垂线间夹角$\gamma$是天顶角。
速度矢量与当地水平面夹角称为航迹仰角，即航迹角$\Phi$。

$$h=|\vec{h}|=|\vec{r}\times \vec{v}|=|r||v|sin\gamma =rvcos\Phi （3）$$

二体轨道几何方程

用$\vec{h}$与（1）做叉乘，
$$\vec{h}\times \ddot{\vec{r}}=-\frac{\mu }{r^3}\vec{h}\times \vec{r}$$
因为$h$是个常数并且有叉乘公式$(a\times b)\times c=b(ac)-a(bc)$
$$\frac{d}{dt}(\vec{h}\times \dot{\vec{r}})=\vec{h}\times \ddot{\vec{r}}+\dot{\vec{h}}\times \dot{\vec{r}}=\vec{h}\times \ddot{\vec{r}}$$$$\vec{h}\times \vec{r}=(\vec{r}\times \vec{v})\times \vec{r}=\vec{v}(\vec{r}\vec{r})-\vec{r}(\vec{v}\vec{r})=\vec{v}{r}^2-r{\vec{e}}_r(\vec{v}\vec{r})$$
其中${\vec{e}}_r是r$方向的单位向量，拿出一个$r$就是为了后面方便约掉
$$-\frac{\mu }{r^3}\vec{h}\times \vec{r}=\frac{\mu }{r^3}\vec{r}\times \vec{h}=-\frac{\mu }{r}\vec{v}+\frac{\mu }{r^2}\vec{v}\vec{r}=-\frac{\mu }{r}\frac{d\vec{r}}{dt}+\frac{\mu }{r^2}\frac{d\vec{r}}{dt}\vec{r}=-\mu \frac{d}{dt}(\frac{\vec{r}}{r})$$$$\frac{d}{dt}(\dot{\vec{r}}\times \vec{h})=\mu \frac{d}{dt}(\frac{\vec{r}}{r})$$
两边积分
$$\dot{\vec{r}}\times \vec{h}=\mu \frac{\vec{r}}{r}+\vec{B}$$
其中B是积分常矢量，方向是从中心引力体指向近拱点。下面画图解释
$f是r$与圆锥曲线上离焦点最近的一点与焦点连线间的夹角

用$\vec{r}$点乘上式可得到标量方程
$$\vec{r}(\dot{\vec{r}}\times \vec{h})=\vec{r}\frac{\mu \vec{r}}{r}+\vec{r}\vec{B}$$$$左边=\vec{r}(\vec{v}\times (\vec{r}\times \vec{v}))=(\vec{r}\times \vec{v})(\vec{r}\times \vec{v})=h^2$$$$h^2=\mu r+rBcosf$$
上面第二步利用了向量混合积的概念，混合积结果是三个向量形成的立体的体积，所以满足结合律。解出$r$得轨道的几何方程为
$$r=\frac{h^2/\mu }{1+(B/\mu )cosf}=\frac{p}{1+ecosf}$$
$p$是半正焦弦，常数$e$是偏心率。e=0时轨迹是圆，0<e<1是椭圆，e=1是抛物线，e>1是双曲线。
综上，航天器轨道运动可以总结如下，

• 圆锥曲线族为二体问题中的航天器唯一可能的运动轨道。
• 中心引力体中心必定是圆锥曲线轨道的一个焦点。
• 当航天器沿同一条圆锥曲线轨道运动时，比机械能保持不变。动能与势能可相互转换。当航天器高度增加时，速度必定变慢，反之变快。
• 航天器绕中心引力体运动，当$\vec{r}、\vec{v}$沿轨道变化时，比角动量$\vec{h}$保持不变——比的意思是质量为1kg，角动量（动量矩）公式$\vec{L}=\vec{r}m\vec{v}$
• 轨道运动总是处于一个固定于惯性空间的平面内

二体轨道几何性质

$$偏心率e=\frac{c}{a}$$$$半正焦弦p=\frac{a^2-c^2}{a}=a(1-e^2)$$

1轨道的近拱点与远拱点

轨道长轴的两个端点称为拱点，离主焦点近的称为近拱点（或者称近X点，如近日点），远的称为远拱点（远X点）。对于圆轨道，这些点不是确定的，对于开曲线（抛物线与双曲线），远拱点没有物理意义。
主焦点至近拱点距离
$$r_p=\frac{p}{1+e}=a(1-e)$$
主焦点至远拱点的距离
$$r_a=\frac{p}{1-e}=a(1+e)$$
在任何圆锥曲线轨道的近拱点或远拱点处，总有位置矢量与速度矢量垂直
$$h=r_p{v}_p=r_a{v}_a$$

2轨道形状

对所有的圆锥曲线轨道都成立这个简单的关系式表明，轨道长半轴a仅与航天器比机械能E有关。进一步说，E仅与轨道上任一点的r\v有关，即

• 圆和椭圆轨道a>0,E<0
• 抛物线轨道a=$\infty$(抛物线的a,c值都是无穷),E=0
• 双曲线轨道a<0,E>0

因为默认以无穷远点为零势能点，航天器向心力与引力平衡，椭圆平均半径等于半长轴，推出E与a的关系如下，
E = 势 能 + 动 能 = { m v 2 2 − G M m r } ∣ m = 1 , v = G M / r = − μ 2 r = − μ 2 a E=势能+动能=\{\dfrac{mv^2}{2} - \dfrac{GMm}{r}\}|_{m=1,v=\sqrt{GM/r}}=-\dfrac{\mu}{2r}=-\dfrac{\mu}{2a} E=势能+动能={2mv2​−rGMm​}∣m=1,v=GM/r ​​=−2rμ​=−2aμ​

所以仅由航天器比机械能的符号就可以确定航天器处于哪种类型的圆锥曲线轨道内。由于$p=h^2/\mu ，E=-\mu /2a$，h单独决定p,E单独决定a，而E\h决定e,即确定了圆锥曲线轨道的具体形状。
$$e=\sqrt{1+\frac{2E{h}^2}{{\mu }^2}}$$

圆型限制性三体问题

只考虑两个相互作用的质点所组成的质点系在惯性系中的运动问题，称为二体问题。三体问题则是研究3个引力体之间的运动问题。限制性三体问题是低三个质量较之于其他两个可以忽略，即$m对M_1,M_2没有引力作用，只考虑M_1和M_2的相互吸引力和对m的引力$
$M_1{M}_2$按照圆锥曲线轨道做无摄动运动，按运动曲线可分为双曲线型限制性三体问题，抛物线型限制性三体问题，圆型限制性三体问题。

假设$Sun$和$planet$围绕其公共质心做圆周运动，则由$Sun、planet、satellite$组成的系统属于原型限制性三体问题。为了描述卫星的轨道运动，建立如下图所示的质心旋转坐标系C-XYZ,M1是$Sun$，M2是$planet$，m是$satellite$。

• 圆心C是质心，$Sun和planet$围绕质心旋转，两者都处在X轴上，$satellite$
  位置随意。旋转时C-XYZ坐标系以Z轴垂直纸面向内做顺时针旋转。

平动点（拉格朗日点）

---

参考《航天器动力学与控制》