![张朝阳的物理课(第二卷)](https://wfqqreader-1252317822.image.myqcloud.com/cover/50/47793050/b_47793050.jpg)
二、小质量物体在大质量天体引力作用下的轨道
要想详细地厘清引力弹弓效应,我们需要计算出小质量物体在大质量天体引力作用下的轨道方程。其实在《张朝阳的物理课》第一卷中就已经推导过行星的轨道方程,不过当时的主要焦点是椭圆轨道。在这里,我们将使用更便捷的方法得到一般性的轨道方程。
简单起见,我们忽略小质量物体自身的引力,并将大质量天体固定于坐标原点。设大质量天体的质量为M,小质量物体的质量为m。在大质量物体的引力场中,物体的动力学方程可以写为
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_1.jpg?sign=1739260957-o8KqmGTrK3CdUxlHPx8xZ2pS5OVBs7Er-0-22868d06dae39e9f11962866a473ab6c)
(1)
我们并不打算直接求解这个矢量微分方程,而是先考虑这个体系的守恒定律。首先,我们有角动量守恒,这一点可以从有心力场引起的力矩上看出:
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_2.jpg?sign=1739260957-98CT2d9uGYr57SUHPZgxKLBbzzr1RBpc-0-101f63e76916777a4fcb8e0970d71aac)
这导致角动量是一个不随时间改变的量。根据向量叉乘的性质,
与
都垂直于
,因此
与
两者总在同一个平面上,这说明物体绕天体的运动是一个平面运动。在这个平面上,以大质量天体的中心为原点建立极坐标系(r,θ),于是物体的速度可以表示为
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_9.jpg?sign=1739260957-dErNQN7txMPfIBQjUz2v3CbHLYFQ3kwN-0-83d7c5a714170089b4ab98dd96b47154)
其中,上方标有一点的r和θ分别表示它们对时间的导数,即与
。
接下来要考察的守恒定律是能量守恒定律,为此,我们在式(1)等号左边点乘物体在时间dt内的位移:
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_13.jpg?sign=1739260957-fnKYB47RgaxUFhFBFM0i3UuykzLm2SEb-0-585af368bdaafd221b7d2e74abbce6b4)
然后,在式(1)右边点乘:
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_15.jpg?sign=1739260957-HLB3C1RrZylsUVfPiGYl3mCdRp38vjeE-0-658e8419744bee43837cda9cbaa3f5cc)
上式使用了标量场的微分关系以及作为保守力场的引力场所满足的条件:
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_17.jpg?sign=1739260957-R5ur9V7Ku6q8mbSGHn56fXeAjvnIRq0E-0-fca5a4233e3722afa0fe9611b35eb11c)
由这些结果我们能够写出
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_18.jpg?sign=1739260957-9HeUWZPunGpxwUACZ6NbhKfbIdX7AVsp-0-a343237841cb55ba280de8e6e5be02b6)
于是
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_19.jpg?sign=1739260957-DVz0OSSyV76UTWUnHmBnOjufgCtbQotf-0-8617fd43950f47a2b074e5a619f020c9)
在物体运动过程中保持不变。事实上,E正是物体的能量。
将速度表达式代入前述两个守恒量中,我们得到
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_20.jpg?sign=1739260957-BuTUubJgTUeq5ZqsN4HCqs98FimtSYqQ-0-e2f52d362aa29361b20c8343014d85aa)
定义常数α=L/m,于是,借此消去E中的
可得
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_23.jpg?sign=1739260957-Td30SoGEoc3aKd10Aawxshl4VmqQwZ24-0-279fa4fc6b9853c1203fce43b34244a6)
(2)
我们只关心轨道方程,即作为θ的函数的r(θ),因此我们使用链式法则将dr/dt进行改写:
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_24.jpg?sign=1739260957-KzCi0F8LgWIuiESgfmlfNZbIQSMLlH7P-0-67a752fb385cddfe1dab90661f4073e7)
其中,撇号表示对θ的一阶导数。将上式代入式(2),我们得到
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_25.jpg?sign=1739260957-ezIBoXeCU7ssmAlM0cMgJ65PBr7lMna7-0-8c01d0243bfbb17c29614be2881b1352)
接着,引入y=1/r,可以将上式改写为
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_26.jpg?sign=1739260957-JEnQu19ydwYYmMNNEDf3cBbtHs87onZ5-0-2b298a90366a936f382baec8bdec8fbd)
这是一个非线性微分方程,在通常情况下是难以求解的。在这里,我们可以通过猜测解的形式来进行求解。方程中同时出现了y与y′ 的同次幂项,这启发我们使用三角函数来试探解的形式。我们不妨设解拥有形式
y=Acosθ+B
将此形式代入方程并整理,可得
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_27.jpg?sign=1739260957-IdwWH79JJwbSNmO3Lxyors2cisZPXHvu-0-7aca0221f1fc9484c56712e2718ecb06)
要想这个方程对任意θ成立,必须让cosθ前面的系数为零,由此得到
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_28.jpg?sign=1739260957-ErigPsLsEvguNejzKISC2BHwPz1pVaoz-0-d4d079387aa512d011a77de41bb9008c)
将其代入前式可以解出A为
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_29.jpg?sign=1739260957-ZMNWjbql2btnlH8TRmIHT5wsObP0da8i-0-7688b7c7170b1b66ec2630a4c87ec5a1)
在上式的开方中只保留了正根,这是因为我们总可以选择θ的起点使得A>0。由这些结果可以得到轨道方程为
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_30.jpg?sign=1739260957-fkklKbmt4hf7ZfnN6YgatOXyzbD6DO3L-0-65a641c8af8d28db1a60c4df536a5646)
此结果表明物体的轨道是圆锥曲线。在θ=0 处,r(θ)取得其最小值,因此在最终的结果中,我们将极轴取在了圆锥曲线的对称轴上,而坐标原点是圆锥曲线的一个焦点。
上述圆锥曲线的类型取决于A与B的相对大小,有如下三种可能性。
(1)A<B,此时Acosθ+B总是大于零的,这使得无论θ取任何值,r(θ)都取有限值,这对应椭圆轨道。
(2)A=B,此时在θ=π 处有Acosθ+B=0。定性地说,这意味着当θ靠近π 时。r将变得任意大,这对应抛物线轨道。
(3)A>B,此时在θ<π 处,Acosθ+B就有机会取到零值,取零值时的θ由A与B决定,其代表了渐近线的方向,这对应双曲线轨道。
对于双曲线轨道,我们能够给出渐近线对应的角方向。如果令其补角为β,则有
![](https://epubservercos.yuewen.com/239267/27207915207960106/epubprivate/OEBPS/Images/txt005_31.jpg?sign=1739260957-juDtdjhCTrdVfqz6nVx3rwZ3vQUi3xO1-0-110083da37d739cc92d5e855df6ae302)
注意,2β是双曲线的两条渐近线在双曲线这一侧的夹角。