# 三体问题：拉格朗日等边三角形解的两种推导 ## 1 背景与设定 在经典三体问题中，拉格朗日（Lagrange）证明了：**任意三个天体在万有引力作用下，可以保持等边三角形的构型进行整体运动**。这个解是天体力学中少数几个能用初等函数严格表达的周期解之一，也是拉格朗日点 $L_4$、$L_5$ 存在的理论基础。 ### 1.1 基本假设 - 三个质点 $P_1, P_2, P_3$，质量分别为 $m_1, m_2, m_3$，均为正实数。 - 任意时刻，三者构成等边三角形，边长为 $l$。 - 整个系统绕其质心做匀速旋转，角速度 $\omega$ 为常数。 - 仅考虑牛顿万有引力，引力常数 $G$。 ### 1.2 符号约定 - $\vec{r}_i$：质点 $i$ 在惯性系中的位置矢量。 - $\vec{r}_{ij} = \vec{r}_j - \vec{r}_i$：从 $i$ 指向 $j$ 的矢量。 - $r_{ij} = |\vec{r}_{ij}| = l$：因等边三角形假设，所有边长相等。 - $\vec{r}_i'$：质点 $i$ 相对于系统质心的位置矢量。 - $\vec{\omega}$：旋转角速度矢量，垂直于运动平面。 --- ## 惯性参考系中的推导（先求合力，再证圆周运动） 这个版本的核心思路是：**不预先假设质点做圆周运动，而是直接计算任意时刻质点受到的引力合力，然后分析这个合力会导致什么样的运动**。 ### 3.1 初始设定 - 三个质点 $P_1, P_2, P_3$ 在**初始时刻**构成一个边长为 $l$ 的等边三角形。 - 设质心 $O$ 静止（或匀速运动，不影响相对运动）。 - 质点 $i$ 在惯性系中的位置矢量记为 $\vec{r}_i$（以质心 $O$ 为原点）。 - 三个质点的初始速度**待定**——我们希望通过计算发现，只有特定速度才能维持等边三角形。 ### 3.2 计算质点 $P_1$ 受到的合力 质点 $P_1$ 受到 $P_2$ 和 $P_3$ 的引力： $ \vec{F}_1 = \vec{F}_{12} + \vec{F}_{13} = \frac{G m_1 m_2}{r_{12}^3} \vec{r}_{12} + \frac{G m_1 m_3}{r_{13}^3} \vec{r}_{13} $ 在等边三角形假设下，$r_{12} = r_{13} = l$，所以： $ \vec{F}_1 = \frac{G m_1}{l^3} \left( m_2 \vec{r}_{12} + m_3 \vec{r}_{13} \right) \tag{3.1} $ **关键技巧**：引入质心关系。由质心定义： $ m_1 \vec{r}_1 + m_2 \vec{r}_2 + m_3 \vec{r}_3 = 0 $ 且 $\vec{r}_{12} = \vec{r}_2 - \vec{r}_1$，$\vec{r}_{13} = \vec{r}_3 - \vec{r}_1$。代入式 (3.1)： $ \begin{aligned} \vec{F}_1 &= \frac{G m_1}{l^3} \left[ m_2 (\vec{r}_2 - \vec{r}_1) + m_3 (\vec{r}_3 - \vec{r}_1) \right] \\ &= \frac{G m_1}{l^3} \left[ m_2 \vec{r}_2 + m_3 \vec{r}_3 - (m_2 + m_3) \vec{r}_1 \right] \end{aligned} $ 利用质心关系 $m_2 \vec{r}_2 + m_3 \vec{r}_3 = -m_1 \vec{r}_1$，代入得： $ \vec{F}_1 = \frac{G m_1}{l^3} \left[ -m_1 \vec{r}_1 - (m_2 + m_3) \vec{r}_1 \right] = -\frac{G m_1}{l^3} (m_1 + m_2 + m_3) \vec{r}_1 $ 记总质量 $M = m_1 + m_2 + m_3$，则有： $ \vec{F}_1 = -\frac{G m_1 M}{l^3} \vec{r}_1 \tag{3.2} $ **结果解读**：式 (3.2) 是一个**极其重要**的结论——在任意质量分布的等边三角形构型中，每个质点受到的引力合力**恰好等于一个指向质心的、大小与到质心距离成正比的有心力**。 ### 3.3 分析合力的物理含义 式 (3.2) 的形式是： $ \vec{F}_1 = -k m_1 \vec{r}_1, \quad \text{其中 } k = \frac{G M}{l^3} $ 这是一个典型的**线性恢复力**（胡克定律形式），但它是由万有引力在特定几何构型下产生的。 将牛顿第二定律 $\vec{F}_1 = m_1 \ddot{\vec{r}}_1$ 代入： $ m_1 \ddot{\vec{r}}_1 = -k m_1 \vec{r}_1 \quad \Rightarrow \quad \ddot{\vec{r}}_1 = -k \vec{r}_1 \tag{3.3} $ **方程 (3.3) 是二维各向同性谐振子的运动方程**！ ### 3.4 求解运动方程 方程 $\ddot{\vec{r}}_1 = -k \vec{r}_1$ 的通解为： $ \vec{r}_1(t) = \vec{A} \cos(\omega t) + \vec{B} \sin(\omega t) $ 其中 $\omega = \sqrt{k} = \sqrt{\frac{G M}{l^3}}$，$\vec{A}$、$\vec{B}$ 是由初始位置和初始速度决定的常矢量。 **两种可能的运动模式**： 1. **椭圆轨道**：一般情况下，质点在 $xy$ 平面上做椭圆运动，椭圆中心在质心。 2. **匀速圆周运动**：当 $\vec{A} \perp \vec{B}$ 且 $|\vec{A}| = |\vec{B}| = R_1$（$R_1$ 为 $P_1$ 到质心的距离）时，轨道为正圆。 ### 3.5 等边三角形构型与圆周运动的匹配 式 (3.3) 虽然表明每个质点都在做简谐振动，但**要维持三个质点始终构成等边三角形，它们的轨道必须满足严格的相位关系**。 可以证明（通过代入另外两个质点的方程并利用质心关系），维持等边三角形的唯一可能是： - 三个质点都做**匀速圆周运动**（而非椭圆）。 - 角速度相同，均为 $\omega = \sqrt{\frac{G M}{l^3}}$。 - 相位相差 $120^\circ$（或 $240^\circ$），即始终构成等边三角形。 若尝试椭圆轨道，则三个质点到质心的距离 $R_i$ 会随时间变化，导致边长 $l$ 变化，但 $l$ 在万有引力公式中是以 $1/l^3$ 出现的，会进一步改变 $k$ 值，破坏谐振子的线性性质。因此，椭圆轨道不能自洽地维持等边三角形（除非特殊情形如三质点同时径向运动——即相似收缩/膨胀解，但那是另一类解，不是拉格朗日等边三角形解）。 ### 3.6 圆周运动解的验证 当三个质点以相同角速度 $\omega$ 做匀速圆周运动时： - 向心加速度大小：$a_c = \omega^2 R_1$ - 由式 (3.2) 产生的加速度大小：$a = \frac{F_1}{m_1} = \frac{G M}{l^3} R_1$ 两者相等给出： $ \omega^2 R_1 = \frac{G M}{l^3} R_1 \quad \Rightarrow \quad \omega^2 = \frac{G M}{l^3} $ 这与之前的结果完全一致。$R_1$ 在等式中消去，说明**无论质点到质心的距离是多少，只要角速度取这个值，圆周运动就能维持**。 ### 3.7 几何自洽性 还需要验证：三个质点以角速度 $\omega$ 做圆周运动时，是否能始终保持等边三角形？ - 由几何学可知，三个点绕同一点（质心）做同频圆周运动，只要半径 $R_1, R_2, R_3$ 和初始相位满足三角形的边角关系，就能保持形状不变。 - 对于等边三角形，三个半径由质量决定（质心位置决定），初始相位互差 $120^\circ$，形状不变性自然满足。 ### 3.8 推导结论 通过"先求合力"的思路，我们得到： 1. 在等边三角形构型下，任意质点受到的引力合力**必然指向质心**，且大小与到质心的距离成正比（线性恢复力）。 2. 因此，每个质点的运动方程是二维谐振子方程。 3. 在谐振子解中，只有**匀速圆周运动**才能自洽地维持等边三角形的几何不变性。 4. 圆周运动的角速度必须满足 $\omega^2 = \frac{G M}{l^3}$，与质量分布无关。 这个推导的优点在于： - **物理图像清晰**：先认识到"合力是指向质心的线性力"，然后自然地引出圆周运动解。 - **不需要预先假设运动形式**，而是从力的表达式反推出可能的运动。 - 顺便揭示了等边三角形解与二维各向同性谐振子的深刻联系。 --- [三体运动 RKF45 高精度模拟](https://www.myphysics-lab.com/physics_classic/santi) | 质点 | 质量 \(m\) | 初始位置 \((x, y)\) | 初始速度 \((v_x, v_y)\) | | :--- | :--- | :--- | :--- | | **P₁** | 3 | (0.8660254035, -0.1666666667) | (0.1443375673, 0.7500000000) | | **P₂** | 2 | (-0.8660254035, 0.8333333333) | (-0.7216878365, -0.7500000000) | | **P₃** | 1 | (-0.8660254035, -1.1666666667) | (1.010425336, -0.7500000000) |