地球扁率的简单推导
By 苏剑林 | 2011-02-08 | 61645位读者 |假如没有自转,单凭物质之间的引力作用,天体应该都会呈现一个很完美(不是绝对完美)的球形。不过绝大多数的天体都存在着自转,因此他们的赤道半径都比极半径长。BoJone粗糙地考虑了一下在引力和惯性离心力的共同作用下,天体所呈现的形状,并与太阳系的一些天体进行对比,发现还是能够吻合到一定程度。现在此和大家分享,供读者参考。
以地球为例,假如没有自转,地球应该是一个球体。这里只在平面考虑,地球方程为$x^2+y^2=R^2$,R是地球半径,即极半径。由于发生了自转,赤道会相应地被“拉长”;越靠近赤道,“拉”得就越长。设地球均匀自转,角速度为$\omega$。没有转动时,地球表面上一点(x,y)(设质量为1)所受到的引力的水平分量为$\frac{GM}{R^2}\cdot \frac{\sqrt{R^2-y^2}}{R}$(我们用含y而不是x的式子写出,因为x轴在旋转前后会变化,而y不变)。这里我们只需要考虑水平分量,因为向心力是水平方向的,大小为$\omega^2 \sqrt{R^2-y^2}$。向心力由引力的水平分量提供,即等价于减少了引力,即旋转后引力的水平分量为$\frac{GM}{R^2}\cdot \frac{\sqrt{R^2-y^2}}{R}-\omega^2 \sqrt{R^2-y^2}$。
接着我们做一个比较主观性的假设:“拉长”的比例与引力的水平分量成反比(稍后我们将会发现这个假设相当正确)。于是我们得到原来的圆形的横坐标被拉长为原来的
$$\frac{\frac{GM}{R^2}\cdot \frac{\sqrt{R^2-y^2}}{R}}{\frac{GM}{R^2}\cdot \frac{\sqrt{R^2-y^2}}{R}-\omega^2 \sqrt{R^2-y^2}}$$
倍。这个结果恰好是一个常数,等于$\frac{1}{1-\frac{\omega^2 R^3}{GM}}$。即旋转后的图形变成了一个椭圆,其方程为$(1-\frac{\omega^2 R^3}{GM})^2 x^2+y^2=R^2$。作近似处理,$\frac{1}{1-\frac{\omega^2 R^3}{GM}}\approx 1+\frac{\omega^2 R^3}{GM}(|\frac{\omega^2 R^3}{GM}| \ll 1)$,即水平方向的长度被拉长了$\frac{\omega^2 R^3}{GM}$。
我们用一些例子来检验一下。就地球来说,查维基可以得到赤道半径6,378km,两极半径6,356km,相差大约22km。用上式计算
$$R\cdot \frac{\omega^2 R^3}{GM}=R\cdot \frac{(\omega R)^2}{v_1^2}=21.76km$$
其中$v_1=7.9km/s$是第一宇宙速度,$R=6356km,\omega=\frac{2\pi}{86400s}$。可见,这是一个相当好的近似。为了证明这不是“纯属巧合”,我们再来用木星的数据进行检验。
木星是太阳系最大的一个行星,赤道半径71,492km,两极半径66,854km,相差4638km(差不多一个地球了,呵呵),自转一周只需要9.9小时。用上式计算
$$R\cdot \frac{\omega^2 R^3}{GM}=R\cdot \frac{(\omega R)^2}{v_1^2}\approx 5100km$$
虽然不是特别精确,不过BoJone认为吻合得还是相当不错了,毕竟这是一个如此粗糙的模型!最后我们得到了一个扁率因子$\frac{\omega^2 R^3}{GM}$!欢迎大家作进一步的修正!
转载到请包括本文地址:https://kexue.fm/archives/1253
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Feb. 08, 2011). 《地球扁率的简单推导 》[Blog post]. Retrieved from https://kexue.fm/archives/1253
@online{kexuefm-1253,
title={地球扁率的简单推导},
author={苏剑林},
year={2011},
month={Feb},
url={\url{https://kexue.fm/archives/1253}},
}
February 9th, 2011
拉长的同时竖直方向会不会变扁呢?
若地球为流体,则极半径会相对于平均半径缩短。貌似这里楼主假设极半径不是原来的平均半径,而是缩短之后的。
是的,你的说法没有错。或者说假设只有赤道拉伸而没有两极收缩。因为求的目的在于一个比例,而比例又很小,可以说各种不确定因素都相互抵消了一部分,累积不明显。
February 9th, 2011
貌似木星自转周期只是9小时多一些哦。
改过来了。
我记得太阳系最扁的是土星,不过也可能是木星。
是土星...
……自己查了一下,发现原来土星更扁。
November 11th, 2019
站长写得真是太精彩了,必须赞一个!不过我要小小地捉一下虫,在这里地球自转的角速度ω应该代入恒星日而不是太阳日,也就是说86400s应该改为86164s(23小时56分4秒)喔。这样算出来的结果就更加接近22km了~~~
哈哈,好的~