拉格朗日点

贡献者： addis

本文存在未完成的内容。

预备知识　雅可比常量

　　¹对于大质量天体 A 和 B 围绕质心作匀速圆周运动的限制性三体问题，是否存在若干个特殊位置，使得处于该位置并满足一定的初始条件的第三质点将在运动中保持与天体 A 和 B 的相对位置不变，即第三质点同样围绕系统质心作匀速圆周运动，且角速度与天体 A 和 B 相同？这样的位置称为 “平动点” 或 “拉格朗日点（Lagrangian point）”。

　　提取式 1 中的势能部分并除以质量，记为等效势 $\widetilde{U}$

\begin{equation} \widetilde{U} = -\frac{Gm_1}{r_1} -\frac{Gm_2}{r_2} -\frac{1}{2}\omega^2 r^2~, \end{equation}

其中 $r_i = \sqrt{(x - x_i)^2 - (y - y_i)^2}$。

　　等效势能在 $xoy$ 平面内部分区域的分布如图 1 所示，图中 A、B 两点为等效势能的理论奇点，在其附近形成 “势能陷阱”。

图 1：等效势能的分布（红色代表高势能，蓝色代表低势能）

　　要计算平衡点，即 $\widetilde{U}$ 的极值点，只需令

\begin{equation} \boldsymbol\nabla \widetilde{U} = \boldsymbol{\mathbf{0}} ~. \end{equation}

得到

\begin{equation} \frac{-Gm_1(x-x_1)}{r_1^3} - \frac{Gm_2(x-x_2)}{r_2^3} + \omega^2 x = 0~, \end{equation}

\begin{equation} -\frac{Gm_1 y}{r_1^3} - \frac{Gm_2 y}{r_2^3} + \omega^2 y = 0~. \end{equation}

即

\begin{equation} -\frac{Gm_1( \boldsymbol{\mathbf{r}} - \boldsymbol{\mathbf{r}} _1)}{r_1^3} - \frac{Gm_2( \boldsymbol{\mathbf{r}} - \boldsymbol{\mathbf{r}} _2)}{r_2^3} + \omega^2 \boldsymbol{\mathbf{r}} = \boldsymbol{\mathbf{0}} ~. \end{equation}

显然这相当于寻找两个万有引力和离心力合力为零的点。式 4 有显然的解 $y = 0$，带入式 3 得

未完成：见注释

1. 使用拉格朗日方程推导

未完成：该小节有待整理

　　保持与天体 A 和 B 的相对位置不变，即要求第三质点在旋转坐标系中满足 $\dot{x}=\dot{y}=\dot{z}=\ddot{x}=\ddot{y}=\ddot{z}=0$，代入第三质点的动力学微分方程组（式 2 ，式 3 ，式 4 ）可得方程组

\begin{equation} -\omega^2 x =- \frac{\partial U}{\partial x} =-\frac{\mu GM}{r_1^3}[x+(1-\mu)R]-\frac{(1-\mu)GM}{r_2^3}(x-\mu R)~, \end{equation}

\begin{equation} -\omega^2 y =- \frac{\partial U}{\partial y} =- \left[\frac{\mu GM}{r_1^3}+\frac{(1-\mu)GM}{r_2^3} \right] y~, \end{equation}

\begin{equation} 0 =- \frac{\partial U}{\partial z} =- \left[\frac{\mu GM}{r_1^3}+\frac{(1-\mu)GM}{r_2^3} \right] z~. \end{equation}

观察式 8 可得 $z=0$，即待求点在天体 A、B 的运动平面内。

　　若 $y=0$，则式 7 恒成立。并且注意到 A、B 作匀速圆周运动时，根据外有引力等于向心力的关系，可得

\begin{equation}\omega^2 =\frac{GM}{R^3}~. \end{equation}

此时式 6 便可化简为

\begin{equation}\frac{x}{R^3} =\frac{\mu}{|x+(1-\mu)R|^3}[x+(1-\mu)R]+\frac{1-\mu}{|x-\mu R|^3}(x-\mu R)~. \end{equation}

此方程没有解析解，只能借助数值方法求得其满足一定精度要求的数值解。可以验证在天体 A、B 之间及两边延长线方向上各有一个解，分别记为 $L_1,L_2,L_3$。

　　若 $y\neq 0$，则有

\begin{equation} \left\{\begin{aligned} \frac{x}{R^3} &=\frac{\mu}{r_1^3}[x+(1-\mu)R]+\frac{1-\mu}{r_2^3}(x-\mu R)\\ \frac{1}{R^3} &=\frac{\mu}{r_1^3}+\frac{1-\mu}{r_2^3} \end{aligned}\right. ~. \end{equation}

解得 $r_1=r_2=R$，即平动点与 A、B 形成正三角形，这样的平动点有两个，分别记为 $L_4,L_5$。

　　拉格朗日点的力学本质是，第三质点在该位置所受的引力合力恰好等于维持匀速旋转所需的的向心力。五个拉格朗日点与 A、B（在旋转坐标系内）的相对位置大致如图 2 所示。

图 2：拉格朗日点位置

1. ^ 参考 Wikipedia 相关页面。