看完了“双不动中心”问题,我们不妨再来看一个貌似简单一点的力学问题,在一个固定质点的引力吸引的基础上,增加一个恒力作用,研究这样的力场中小天体的运动。

咋看上去这个问题比“双不动中心”简单多了,至少运动方程也显得更简单:
$$\ddot{vec{r}}=-GM\frac{\vec{r}}{|\vec{r}|^3}+\vec{F}$$

其中$\vec{F}$是一个常向量。不过让人比较意外的是,这个问题本质上和“双不动中心”是一样的,它可以看作是双不动中心问题的一个极限情况。而且它们的解法也是惊人地相似,下面我们就来分析这一个过程。

首先很容易写出这个方程的能量守恒积分:
$$1/2 \dot{vec{r}}^2-GM\frac{1}{|\vec{r}|}-\vec{F}\cdot \vec{r}=h$$

同样我们只考虑平面情况。选择适当的参考点,使得$\vec{F}=(F,0)$,于是可以记该系统的动能为$T=1/2 (\dot{x}^2+\dot{y}^2)$,势能为$U=-\frac{GM}{\sqrt{x^2+y^2}}-Fx$。能量是
$$T+U=h$$
拉格朗日函数是$L=T-U$。

为了消去势能函数U中的根号,我们首先会想到的是使用极坐标。的确,这个坐标系能够消去根号,但不能让我们进一步做下去了。那么,还有没有其他方法能够消去根号的呢?

一个忠实的数学粉丝会知道,在2000多年前,伟大的毕达哥拉斯就想出了一种“勾股数”的表达形式:$(a^2-1)^2+(2a)^2=(a^2+1)^2$,后来,人们进一步发现了恒等式$(a^2-b^2)^2+(2ab)^2=(a^2+b^2)^2$。我们发现,左端的形式正对应着$x^2+y^2$,于是我们可以考虑利用这个公式来进行变换。

现在问题来了,在$x^2+y^2$,x,y是对称的,而在$(a^2-b^2)^2+(2ab)^2$,括号内的两个量不是对称的,于是我们要判断究竟是$x=a^2-b^2$还是$x=2ab$适合?正确的答案是
$$\begin{aligned}x=a^2-b^2 \\ y=2ab\end{aligned}$$

为什么要这样取?BoJone并没有很好的判断方法,一开始BoJone就把两个给调换过来了,结果计算到最后化简不了,就把两者换过来了。所以说,这是一种“事后判断法”,

变换之后:$T=2(a^2+b^2)(\dot{a}^2+\dot{b}^2),U=-\frac{GM}{a^2+b^2}-F(a^2-b^2)$

现在看上去这个变换还是挺不错的,它保持了方程原有的对称性。下面的处理方法几乎和“双不动中心”问题一模一样了。请看——

拉格朗日方程$\frac{d}{dt}(\frac{\partial L }{\partial \dot{a}})=\frac{\partial L }{\partial a}$

即$\frac{d}{dt}[4\dot{a}(a^2+b^2)]=4a(\dot{a}^2+\dot{b}^2)-\frac{\partial U }{\partial a}$

各项乘以$\dot{a}(a^2+b^2)$得
$$\frac{d}{dt}[2\dot{a}^2 (a^2+b^2)^2]=4a\dot{a}(\dot{a}^2+\dot{b}^2)(a^2+b^2)-\frac{\partial U }{\partial a}\dot{a}(a^2+b^2)$$

利用能量守恒积分进行代换
$$\frac{d}{dt}[2\dot{a}^2 (a^2+b^2)^2]=2a\dot{a}(h-U)-\frac{\partial U }{\partial a}\dot{a}(a^2+b^2)$$

右端变成
$$\begin{aligned}\dot{a}\frac{\partial }{\partial a}[(h-U)(a^2+b^2)] \\ =\dot{a}\frac{\partial }{\partial a}[(h(a^2+b^2)+GM+F(a^4-b^4)] \\ =\dot{a}\frac{\partial }{\partial a}(h a^2+F a^4) \\ =\frac{d}{dt}(h a^2+F a^4)\end{aligned}$$

至此,我们应该明白为什么要$x=a^2-b^2$而不是$x=2ab$了,如果选取后者的话,就无法将a的项和b的项相互分离了。

于是我们有积分:$2\dot{a}^2 (a^2+b^2)^2=h a^2+F a^4+C$

由能量积分变形:$2(\dot{a}^2+\dot{b}^2) (a^2+b^2)^2=h (a^2+b^2)+GM+F (a^4-b^4)$

对比两个式子,有$2\dot{b}^2 (a^2+b^2)^2=h b^2+GM-F b^4-C$

就有$(\frac{d b}{d a})^2=\frac{hb^2-F b^4+(GM-C)}{ha^2+F a^4+C}$

积分得:$\int \frac{da}{\sqrt{ha^2+F a^4+C}}=\int \frac{1}{\sqrt{hb^2-F b^4+(GM-C)}}$

看到了吧?整个过程完全就是“双不动中心”问题求解过程的“翻版”!!

而事实上,坐标变换$x=a^2-b^2,y=2ab$在数学上叫做“抛物线坐标变换[引用]。这是否会让我们想起圆锥曲线的相关知识?抛物线其实是椭圆的一种极限情况(离心率趋向于1),于是抛物线坐标其实也就是椭圆坐标的一个极限情况,而本题也是“双不动中心”问题的一个极限情况了。至于怎么对应起来,话说BoJone还没有弄清楚细节,所以不便多说了。

转载到请包括本文地址:https://kexue.fm/archives/1358

更详细的转载事宜请参考:《科学空间FAQ》

如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。

如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!

如果您需要引用本文,请参考:

苏剑林. (May. 14, 2011). 《“二体+恒力”问题 》[Blog post]. Retrieved from https://kexue.fm/archives/1358

@online{kexuefm-1358,
        title={“二体+恒力”问题},
        author={苏剑林},
        year={2011},
        month={May},
        url={\url{https://kexue.fm/archives/1358}},
}