|
kuing
发表于 2019-2-26 22:55
至于压力,我以前也没思考过,刚才试着自己推,竟然得出中心压强无穷大
这感觉不太合常理,但我想来想去也不知道问题出在哪,大家帮忙瞧瞧:
先研究简化情形,假设液体密度均匀,取距离球心 `x` 至 `x+\Delta x` 的一层球壳,将它记为 `K`,那么 `K` 的两边的压力之差就等于它所受到的引力(吧?问题最大可能出在这里)
设距离球心 `x` 处的压强为 `p(x)`,则压力之差显然是
\[4\pi x^2p(x)-4\pi(x+\Delta x)^2p(x+\Delta x);\]
由链接中的结论知 `K` 的外部对它是没引力的,所以只需计算 `K` 以内的球体对它的引力。
记半径为 `x` 的球的质量为 `m(x)`,由密度均匀知 `m(x)=(x/R)^3M`,其中 `R`, `M` 为地球半径和质量,则 `K` 的质量为 `m(x+\Delta x)-m(x)`,故此引力为
\[\frac{Gm(x)(m(x+\Delta x)-m(x))}{r^2},\]其中 `r\in (x,x+\Delta x)`。
以上两式相等,再令 `\Delta x\to0`,就有
\[-(4\pi x^2p(x))'=\frac {Gm(x)m'(x)}{x^2},\]把 `m(x)` 的表达式代进去化简后就是
\[2p(x)+xp'(x)=-\frac {3GM^2}{4\pi R^6}x^2,\]解这个微分方程得
\[p(x)=-\frac {3GM^2}{16\pi R^6}x^2+\frac C{x^2},\]由 `p(R)=0` 确定常数 `C` 后,最终得
\[p(x)=\frac {3GM^2(R^4-x^4)}{16\pi R^6x^2},\]这样 `p(0)=\infty`,即球心处的压强无穷大……
到底问题出在哪?是“压力差=引力”不对,还是根本不能假设密度均匀?还是其他……? |
|