“二体+恒力”问题
By 苏剑林 | 2011-05-14 | 16636位读者 |看完了“双不动中心”问题,我们不妨再来看一个貌似简单一点的力学问题,在一个固定质点的引力吸引的基础上,增加一个恒力作用,研究这样的力场中小天体的运动。
咋看上去这个问题比“双不动中心”简单多了,至少运动方程也显得更简单:
$$\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}},
}
最近评论