二体问题之2：积分常数

注：笔记，肯请批评指正。
积分常数的推导借助于二体运动方程]：
$\ddot{\mathbf{r}}+\mu \frac{\mathbf{r}}{r^3}=0$

1. 角动量积分$\mathbf{h}$

角动量定义为$$\mathbf{h}=\mathbf{r}\times \mathbf{v}$$

对它求导得，$$\dot{\mathbf{h}}=\mathbf{v}\times \mathbf{v}+\mathbf{r}\times \ddot{\mathbf{r}}$$
其中，前一部分相同向量叉乘得0；且由二体运动方程可知，和同向，二者叉乘同样的0。

由此可知，$$\mathbf{h}=\mathbf{r}\times \mathbf{v}=c$$
根据[[航天器轨道方程]] 的定义可以得到角动量积分和轨道集合参量半通径$p$的关系，即$h=\sqrt{p\mu }$。

2. 能量积分

在二体运动方程右乘$\mathbf{v}$ , 得$$\dot{\mathbf{r}}\cdot \ddot{\mathbf{r}}+\mu \frac{\mathbf{r}\cdot \dot{\mathbf{r}}}{r^3}=0$$
根据[[矢量计算]]公式1：$\mathbf{a}\cdot \dot{\mathbf{a}}=a\dot{a}=(\dot{\frac{1}{2}{a}^2)}$得$$(\dot{\frac{1}{2}{v}^2)}+\frac{\mu }{r^3}(r\dot{r})=0$$
$$(\dot{\frac{1}{2}{v}^2)}+\dot{(-\frac{\mu }{r})}=0$$
得到积分常量$\mathbf{E}$
$$\mathbf{E}=\frac{1}{2}{v}^2-\frac{\mu }{r}$$

根据[[航天器轨道方程]] ，远近地点速度和位置矢量垂直的特性，可以得到能量积分与轨道几何参量的关系。

带入$r_p=\frac{P}{1+e}$ , $v_p=\sqrt{p}(1+e)$ ， 且$a=\frac{P}{1-e^2}$。
$$\begin{array}{rl}& E=\frac{v^2}{2}-\frac{\mu }{r}\\ & =\frac{V_p^2}{2}-\frac{\mu }{r_p}\\ & =\frac{\mu }{2p}\left(1+2e+e^2\right)-\frac{\mu }{2P}(2+2e)\\ & =\frac{\mu }{2p}\left(e^2-1\right)\\ & =-\frac{\mu }{2}\frac{1-e^2}{p}\\ & =-\frac{\mu }{2a}\end{array}$$
即
$$E=-\frac{\mu }{2a}$$

3. 拉普拉斯积分$L$

在二体运动方程左侧叉乘角动量$h=\mathbf{r}\times \mathbf{v}$，得$$\mathbf{h}\times \ddot{\mathbf{r}}+\frac{\mu }{r^3}(\mathbf{r}\times \mathbf{v}\times \mathbf{r})=0$$

第一部分$$\vec{h}\times \ddot{\vec{r}}=\frac{d}{dt}(\vec{h}\times \vec{v})=-\frac{d}{dt}(\vec{v}\times \vec{h})$$
第二部分，根据矢量计算公式
$$\begin{array}{rl}& \frac{\mu }{r^3}(\vec{r}\times \dot{\vec{r}}\times \vec{r})\\ =& \frac{\mu }{r^3}[(\vec{r}\vec{r})\cdot \dot{\vec{r}}-(\dot{\vec{r}}\vec{r})\cdot \vec{r}]\\ =& \frac{\mu }{r^3}\left(r^2\dot{\vec{r}}-r\dot{r}\vec{r}\right)\\ =& \mu \left(\frac{1}{r}\dot{\vec{r}}-\frac{1}{r^2}\dot{r}\vec{r}\right)\\ =& \mu \frac{d}{dt}\left(\frac{1}{r}\vec{r}\right)\end{array}$$
综合可得
$$-\frac{d}{dt}(\vec{v}\times \vec{h})+\frac{d}{dt}\left(\frac{\mu }{r}\vec{r}\right)=0$$
即推导出拉普拉斯常量
$$\vec{L}=\vec{v}\times \vec{h}-\frac{\mu }{r}\vec{r}$$

由[[航天器轨道方程]]$r=\frac{p}{1+e\cos \theta }$可知，当$\theta =0,\vec{L},\vec{r}$同向时，$r=\frac{p}{1+e}=r_p$为近地点矢量，即拉普拉斯积分指向近地点。同时，由圆锥曲线方程可知，其大小$L=e\mu$。

4. 耦合关系

1. 方向上，拉普拉斯常量位于轨道面内，垂直于角动量积分。向量r，v，h相互垂直，关系如下图
   

2. 大小上，拉普拉斯常量有约束
   $$L^2=2E{h}^2+{\mu }^2$$
   在近地点有$h=rv$，则
   $$L^2=(vh-\mu {)}^2=v^2{h}^2-2vh\mu +{\mu }^2$$

$$\begin{array}{rl}2E{h}^2& =2\left(\frac{1}{2}{v}^2-\frac{\mu }{r}\right)h^2\\ & =v^2{h}^2-\frac{2\mu h^2}{r}\\ & =v^2{h}^2-\frac{2\mu rvh}{r}\end{array}$$