精确自由落体运动定律的讨论(二)
By 苏剑林 | 2010-01-09 | 62223位读者 |之前在这篇文章中,我们使用过一个牛顿引力场中的自由落体公式:
\begin{equation}t=\sqrt{\frac{r_0}{2GM}}\left(r_0 \cdot \arctan \sqrt{\frac{r_0 -r}{r}}+\sqrt{r(r_0 -r)}\right)\label{eq:1}\end{equation}
我们来尝试一下推导出这个公式来。
同时,站长在逐渐深入研究的过程中,发现微分方程极其重要。以前一些我认为不可能解决的问题,都用微分方程逐渐解决了。在以后的文章里,我们将会继续体验到微分方程的伟大魔力!因此,建议各位有志研究物理学的朋友,一定要掌握微分方程,更加深入的,需要用到偏微分方程!
首先,质量为$m$的物理在距离地心$r$处的引力为$\frac{GMm}{r^2}$,根据牛顿第二定律$F=ma$,自然下落的物体所获得的加速度为$\frac{GM}{r^2}$。假设物体从距离地心$r$开始向地心自由下落,求位移$s$关于$t$的函数$s=s(t)$.
根据根据加速度的定义,我们有:$\frac{d^2 s}{dx^2}=a=\frac{GM}{(r-s)^2}$,于是问题实质就是解常微分方程
\begin{equation}s''=\frac{GM}{(r-s)^2}\end{equation}
接着我们令$s'=v$,则$s''=v(\frac{dv}{ds})$,代入上式:$GM(r-s)^{-2} ds=vdv$,两端积分:
\begin{equation}\begin{gathered}\int vdv = \int GM(r - s)^{ - 2} ds = - \int GM( r - s)^{ - 2}d(r - s) \\
\Downarrow\\
\frac{1}{2} v^2 = GM[(r - s)^{ - 1} + C_1]\end{gathered}\end{equation}
根据实际情况,当$t=0$时,$v=s=0$,推出$C_1=-r^{-1}$,即
\begin{equation}\frac{1}{2}v^2 = GM[(r - s)^{ - 1} - r^{ - 1}] \quad\Rightarrow\quad \frac{ds}{dt} = v = \sqrt {\frac{2GM}{r}} \sqrt {\frac{s}{r - s}}\end{equation}
两端积分:
\begin{equation}\int dt = \sqrt {\frac{r}{2GM}} \int (\sqrt {\frac{r - s}{s}} )ds =2\sqrt {\frac{r}{2GM}} \int (\sqrt {r - s} )d(s^{0.5}) \end{equation}
令$s^{0.5}=P$,则
\begin{equation}2\int (\sqrt {r - s} )d(s^{0.5}) =2\int (\sqrt {r - p^2} )dp\end{equation}
根据积分公式:
\begin{equation}\int \sqrt{a^2-x^2} dx =\frac{a^2}{2} \arcsin \frac{x}{a} +x/2 \sqrt{a^2-x^2}+C\end{equation}
得出
\begin{equation}2\int (\sqrt {r - p^2} )dp=r\cdot \arcsin \frac{p}{\sqrt{r}}+p \sqrt{r-p^2}+C\end{equation}
换回$s$,有:
\begin{equation}2\int (\sqrt {r - s} )d(s^{0.5}) =r\cdot \arcsin \sqrt{\frac{s}{r}}+\sqrt{s(r-s)}+C\end{equation}
则:
\begin{equation}t=\sqrt{\frac{r}{2GM}}\left(r\cdot \arcsin\sqrt{\frac{s}{r}}+\sqrt{s(r-s)}+C\right)\end{equation}
当$s=0$时,$t=0$,则$C=-\pi r / 2$,得出:
\begin{equation}t=\sqrt{\frac{r}{2GM}}\left(r\cdot \arcsin\sqrt{\frac{s}{r}}+\sqrt{s(r-s)}\right)\end{equation}
而根据反三角函数公式:$\arcsin \frac{a}{b} = \arctan\sqrt{\frac{a^2}{b^2-a^2}}$,即得
\begin{equation}t=\sqrt{\frac{r}{2GM}}\left(r\cdot\arctan\sqrt{\frac{s}{r-s}}+\sqrt{s(r-s)}\right)\label{eq:2}\end{equation}
不难看出,$\eqref{eq:1}$和$\eqref{eq:2}$两式是等价的,$\eqref{eq:2}$的$r$实际上是$\eqref{eq:1}$的$r_0$,而$s$则相当于$r_0 - r$。
转载到请包括本文地址:https://kexue.fm/archives/337
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Jan. 09, 2010). 《精确自由落体运动定律的讨论(二) 》[Blog post]. Retrieved from https://kexue.fm/archives/337
@online{kexuefm-337,
title={精确自由落体运动定律的讨论(二)},
author={苏剑林},
year={2010},
month={Jan},
url={\url{https://kexue.fm/archives/337}},
}
February 19th, 2010
微积分已经忘光了。
高一学生就懂2微积分啊。真不简单。
February 19th, 2010
你是生物专业的吧?或许哪天可能用到微积分的哦(在考虑那些什么结构之类的)...
我其实也是自学的而已,我身边没有一个同好。
February 20th, 2010
I
呵呵,我的专业是金属材料,属于大的机械类。大学时学的普通物理、高等数学等都忘光了。大学时学的东西,在现在最有用的是英语和机械制图。呵呵。
II
身边没有同好没有问题,网上有。
@I
我也得学好英语...有一句经典话:当你研究下去的时候,你会发现只有英语的资料可以看
@II
这话好,网络把四海之内的同好都联系起来了
February 20th, 2010
I
是的,很多前沿的资料都是英文的。对于一个搞(学术)研究的人,英文是不可或缺的。
II
这是互联网的神奇之处。
February 21st, 2010
(抒情中)觉得同站长经历挺像的,同是天涯沦落人啊!
或许这是人生的必经阶段,呵呵
May 16th, 2010
爱压,看了后我自己也了证一遍,一样耶。你好强
不能算什么,只是一个积分的推导而已
June 8th, 2015
同高一 同正在学微积分 同没有同好……