图

雅可比常量

预备知识 限制性三体问题,拉格朗日方程

   在“限制性三体问题”中,我们推导了第三天体 C 的动能、势能、角动量表达式.

   这里对限制性三体问题进一步特殊化,假设原二体系统围绕质心作圆周运动,即大质量天体 A 和 B 的距离 $R$ 在运动过程中不变,则在旋转坐标系中,A、B 是相对静止的.并且天体 A、B 将作匀速圆周运动(参考“轨道参数-时间变量”中关于圆轨道的讨论),即角速度 $\omega$ 也恒定不变.

动力学方程

   三体系统在空间中自由度为 9,但由于天体 A 和 B 的运动模式已知,故限制性三体系统的自由度为 3,可用质点 C 在旋转坐标系中的广义坐标 $x$,$y$,$z$ 描述.将质点 C 的动能 $T$(式 3 )和势能 $U$(式 4 )代入拉格朗日方程

\begin{equation}\dv{t} \pdvTwo{(T-U)}{\dot q_i} = \pdvTwo{(T-U)}{q_i} \quad (q_i\; \rightarrow \; x,y,z) \end{equation}
可得质点 C 的动力学微分方程组
\begin{align} -\omega^2 x - 2\omega\dot{y} + \ddot{x} &=- \pdvTwo{U}{x}=-\frac{\mu GM}{r_1^3}[x+(1-\mu)R]-\frac{(1-\mu)GM}{r_2^3}(x-\mu R)\\%2 -\omega^2 y +2\omega\dot{y} + \ddot{y} &=- \pdvTwo{U}{y}=-\frac{\mu GM}{r_1^3}y-\frac{(1-\mu)GM}{r_2^3}y\\%3 \ddot{z} &=- \pdvTwo{U}{z} =-\frac{\mu GM}{r_1^3}z-\frac{(1-\mu)GM}{r_2^3}z\end{align}

雅可比常量

   由于质点 C 的角动量在惯性坐标轴上的分量和机械能分别守恒,故这些物理量的线性组合也是守恒量.为统一量纲,可以令角动量分量乘以角速度.在所有这些线性组合中,有一种组合具有实际的物理意义,记这个组合量为 $E$

\begin{equation}\begin{aligned} E &=(T+U)-\omega L_z\\ &=\frac{1}{2}m\qtyRound{\dot{x}^2-\dot{y}^2+\dot{z}^2}-\frac{1}{2}m\omega^2\qtyRound{x^2+y^2}-\frac{\mu GMm}{r_1}-\frac{(1-\mu)GMm}{r_2} \end{aligned} \end{equation}
物理量 $E$ 具有能量的量纲.其中 $\frac{1}{2}m\qtyRound{\dot{x}^2-\dot{y}^2+\dot{z}^2}$ 是质点 C 相对旋转坐标系的动能,$E$ 的其余部分只与位置相关,故可解释为势能.在势能部分,引力势能 $-\mu GMm/r_1-(1-\mu)GMm/r_2$ 与坐标系的选取无关,而 $-\frac{1}{2}m\omega^2\qtyRound{x^2+y^2}$ 是由于坐标系 $xOy$ 的旋转引起的,可将其解释为质点 C 所受的向心力由于坐标轴自身的旋转而产生的势能.守恒量 $E$ 称为雅可比常量,可看作是第三天体相对于旋转坐标系具有的能量.

推导

   将式 2 乘以 $\dot{x}$,式 3 乘以 $\dot{y}$,式 4 乘以 $\dot{z}$,再将三个式子相加,可得

\begin{equation} \begin{aligned} &\dot{x} \ddot{x}+\dot{y}\ddot{y}+\dot{z}\ddot{z}-\omega^2\qtyRound{x\dot{x}+y\dot{y}}\\ &=-\frac{\mu GM}{r_1^3}[x\dot{x}+(1-\mu)R\dot{x}+y\dot{y}+z\dot{z}]-\frac{(1-\mu)GM}{r_2^3}(x\dot{x}-\mu R\dot{x}+y\dot{y}+z\dot{z}) \end{aligned} \end{equation}
其中,注意到
\begin{align} \dot{x}\ddot{x}+\dot{y}\ddot{y}+\dot{z}\ddot{z} &=\frac{1}{2}\dv{t}\qtyRound{\dot{x}^2+\dot{y}^2+\dot{z}^2}\\%7 \omega^2\qtyRound{x\dot{x}+y\dot{y}} &=\frac{1}{2}\omega^2\dv{t}\qtyRound{x^2+y^2}\\%8 -\frac{\mu GM}{r_1^3}[x\dot{x}+(1-\mu)R\dot{x}+y\dot{y}+z\dot{z}] &=\dv{t}\qtyRound{\frac{\mu GM}{r_1}}\\%9 -\frac{(1-\mu)GM}{r_2^3}(x\dot{x}-\mu R\dot{x}+y\dot{y}+z\dot{z}) &=\dv{t}\qtyRound{\frac{(1-\mu) GM}{r_2}}\end{align}
将以上四式代入式 6 可得
\begin{equation}\frac{1}{2}\dv{t}\qtyRound{\dot{x}^2+\dot{y}^2+\dot{z}^2}-\frac{\omega^2}{2}\dv{t}\qtyRound{x^2+y^2}-\dv{t}\frac{\mu GM}{r_1}-\dv{t}\frac{(1-\mu) GM}{r_2}=0 \end{equation}
于是
\begin{equation} \frac{1}{2}\qtyRound{\dot{x}^2+\dot{y}^2+\dot{z}^2}-\frac{1}{2}\omega^2\qtyRound{x^2+y^2}-\frac{\mu GM}{r_1}-\frac{(1-\mu) GM}{r_2}=Const \end{equation}
式 12 就是单位质量的雅可比常量.

致读者: 小时物理百科一直以来坚持所有内容免费且不做广告,这导致我们处于日渐严重的亏损状态。长此以往很可能会最终导致我们不得不选择商业化,例如大量广告,内容付费,会员制,甚至被收购。因此,我们鼓起勇气在此请求广大读者热心捐款,使网站得以健康发展。如果看到这条信息的每位读者能慷慨捐助 10 元,我们几天内就能脱离亏损状态,并保证网站能在接下来的一整年里向所有读者继续免费提供优质内容。感谢您的支持。
—— 小时(项目创始人)

编辑词条 返回目录 返回主页 捐助项目 © 小时物理百科 保留一切权利