地球在自轉。現在挖一個隧道,使得一個小球從一端進入,從另一端出來,當中僅受引力作用,不和隧道壁有任何接觸。這樣的隧道是直線嗎?如果不是,是怎樣的曲線?出口和入口是地球上相對的兩個點嗎?

假設地球密度分布均勻,隧道入口與地球表面垂直。
相同問題:
如果地球上挖了一個洞,出口在地球的另一端,一個人跳進洞中,請問他會怎樣運動?
如果有個打通地球的通道人跳下去會怎樣?


陳浩同志在考我大學物理……

以地球為參考系(取直角坐標系,取自轉矢量方向為z軸),則要受到Coriolis力,聯合重力和Coriolis力,我們得到的加速度為:

分量形式為:

前兩式為耦合的微分方程,用差分法進行數值解。

  • 初始條件,地表,緯度45度,經度0度(經度不影響,下面MATLAB代碼中你可以自己改緯度,已注釋)。則

  • 差分解的話,要將3個二階的寫成6個一階的(x、z方向初速度為0,y方向初速度因為地球自轉,是ωRcos(θ)):

  • 北緯45度結果(8萬秒):

    俯視圖:

  • 北緯45度結果(100萬秒):

  • 俯視圖:

————————————————————————————————————————————
萬有引力常量:

地球質量:

地球半徑:

地球自轉角速度:

————————————————————————————————————————————MATLAB代碼,要玩的自己複製去玩吧:
function ball_in_earth( )
% lambda是緯度(正值為北緯,負值為南緯)
lambda = 45;
theta = (90 - lambda)*2*pi/360;
% 緯度lambda轉為球坐標的theta
G = 6.674e-11;
M = 5.9722e24;
R = 6.3781e6;
omega = 7.2921e-5;

x0 = sin(theta)*R;
y0 = 0;
z0 = cos(theta)*R;

vx0 = 0;
vy0 = omega*R*cos(theta);
vz0 = 0;


[t, y] = ode45(@f, [ 0, 8e4], [ x0, y0, z0, vx0, vy0, vz0]);

plot3( y(:,1), y(:,2), y(:,3));
xlabel("x");
ylabel("y");
zlabel("z");

end

function dy = f( t, y)
G = 6.674e-11;
M = 5.9722e24;
R = 6.3781e6;
omega = 7.2921e-5;

dy = [ y(4); y(5); y(6); -G*M/R^3*y(1)+2*omega*y(5); -G*M/R^3*y(2)-2*omega*y(4); -G*M/R^3*y(3)];

end
——————————————————————————————————————————
代碼中倒數第13行的8e4代表的是模擬時間,100萬秒就是1e6了。
8萬秒和100萬秒雖然側視圖不一樣,但是俯視圖一樣。側視圖不一樣的原因可能是MATLAB進行數值計算的精度有限,求解釋、求高精度解法。
————————————————————————————————————————————
如果小球沒有初速度,即你從外太空把它放下去,則代碼中vy0=0,結果如下(北緯45度):

————————————————————————————————————————————
(數值解完了,解析解還沒寫完。)來個解析解,記D為對時間求導的算符,則

兩式相加,得到:

其中,a和b定義為:

由x的特徵方程

得到結果:
r =
0
0.000000381408422 + 0.001758880872482i
0.000000381408422 - 0.001758880872482i
-0.000000762816843
於是,

可見是有衰減項的。(具體解有點崩潰,想解析地解的同志們可以自己玩一下)


當小球在赤道,相對於地球的速度為0的情況下,我用HTML5模擬了一下,於是得到:

圖中的紅色線條就是隧道(or小球的運動軌跡),所有的隧道都是在地球的赤道平面上,這個圖所在的平面就是赤道平面。我只寫了個簡單的二維物理模擬程序,所以暫時只能做赤道上的情況。當小球的初始位置在其它位置時,得到的結果一定是三維空間中的一條複雜曲線,不是簡單的用一張二維的圖片能夠展現的。

另外,求數學大牛用公式分析此運動啊。我嘗試著寫了一下方程,發現很複雜,我微積分太弱了,拜一下牛頓大神。

當然,如果修改不同的模擬參數(比如小球的質量,地球的質量,地球的半徑等等等等),還可以得到很多不一樣的結果,比如:

似乎挺漂亮的:-)PS:請無視某些多餘的直線,謝謝...
PPS:去這裡http://bigbing.org/holybang/html_js_src/earth.html 取代碼或者在線玩


哈哈,居然碰到跟我有一樣想法的人。我提這個問題的時候是初二,那年正好學習引力定律。我就像物理老師提了這個問題,當時的物理老師並沒有回答我。而後來,我又提問了大學物理老師。
如果是太陽為參考系,這個物體的運動軌跡應該是一個以地球公轉軌道為收斂極限的螺旋形狀。


應該不是直線
如果是直線,水池放水時就不會產生漩渦了
然後我在想,這個題目,應近似等價於——
小球在空間中繞與地球等重的質點做「軌道」運動


「地球在自轉。現在挖一個隧道,使得一個小球從一端進入,從另一端出來,當中僅受引力作用,不和隧道壁有任何接觸。這樣的隧道是直線嗎?如果不是,是怎樣的曲線?出口和入口是地球上相對的兩個點嗎?
地球密度分布是均勻的
隧道入口與地球表面垂直」
——————————————————————————
1.至少存在這樣的兩個點,即南北極點,連線是直線。
2.如果入口在赤道。由於小球不受切向的力,而只受到徑向的地心引力,小球初始時刻具有地球自轉表面的切向線速度(日行八萬里)。問題可以部分等效為一個小型天體被一個大型天體俘獲到近軌道然後又逃逸。「隧道的軌跡是沿地球自轉方向的對稱的雙螺旋線」(我這一段不太像對的,因為這個圖形很彆扭,是個w形狀,如果是s形的話轉動方向又錯了。可以看@金秉文 的答案,明顯美),對稱點是地心,「出口關於地心對稱」(恩,這裡也存疑,同上一個括弧里的原因)。
3.入口在地表面其他區域:我猜的,是1和2的疊加。

感謝@陳浩


跟@祁育 一樣,一下子就想到了劉慈欣的《地球大炮》,非常贊!寫了很多相關的細節。
這裡摘一段:

這是一條貫穿地球的隧道。
「它是走切線,還是穿過地心?」沈華北問,只是思維以語言的形式冒了一下頭。
  「聰明的頭腦,這麼快就想到了!」鄧洋驚嘆道。
  「很像他兒子。」有人跟著說,聽上去可能是中部斷裂災難留下的孤兒。
  「是穿過地心,由中國的漠河穿過地球到達南極大陸的最東端南極半島。」鄧洋回答沈華北說。
  「剛才那座城市是漠河?!」
  「是的,它因作為地球隧道起點而繁榮起來。」
  「據我所知,從那裡貫穿地球應該到達阿根廷南部。『』」不錯,但隧道有輕微的彎曲。「
  「既然隧道是彎曲的,我會不會撞上井壁呢?」
  「如果隧道筆直地直達阿根廷,你倒是肯定會撞上,那種筆直的地球隧道只有在貫穿兩極之間的地軸上才能實現,這種與地軸成一定角度的隧道必須考慮地球的自轉因素,它的彎曲正好能讓你平滑地通過。」
  「呵,偉大的工程!」沈華北由衷地讚歎道。


如果不碰到隧道壁,離地心越來越近的時候,角速度會發生變化,所以隧道是彎曲的


為什麼不是空間振子加勻速轉動。
高票沒有計算離心加速度


0、前提假設
不考慮地球自轉,地球密度為
ho ,將跳將跳到洞里的人視為質點,質量為m,挖的通道壁光滑,萬有引力常數為G。任意時刻人到達的位置距離球心r,並以球心O為球心做一個半徑為r的球面(圖中紅色的線),分兩種情況,如圖所示。

1、挖的通道通過球心
建立如圖XOY軸的坐標系,人從A點跳入通道,到達C點,此時r=x,紅線以外的球殼對人的萬有引力為零,所以人僅受到紅線以內的球體的萬有引力

式中k為常數,

由於物體在與位移成正比的恢復力作用下,在其平衡位置附近按正弦規律作往複的運動,即簡諧運動,由F=kx的形式可知,此時人滿足簡諧運動的條件。2、挖的通道不通過球心
建立如圖X"O"Y軸的坐標系,OO"距離為d,人從B點跳入通道,到達D點,此時

和1中同理,紅線以外的球殼對人的萬有引力為零,所以人僅受到紅線以內的球體的萬有引力,但注意此時將人受到的萬有引力分解為平行於通道和垂直於通道兩個方向的分力,y方向的分力為垂直於通道的力,與通道壁給人的正壓力平衡,不再考慮,x"方向的分力為

式中k用1中相同,所以此時人也滿足簡諧運動的條件。綜上,在前提假設條件下,人將做簡諧運動(不用解釋什麼叫簡諧運動了吧?)


個人靠著中學物理僅有的記憶認為,在不考慮阻力什麼的情況下,這個物體應該是在這個洞里來回擺動。


推薦閱讀:

是什麼使星球有地心引力?是什麼使原子核吸引電子?
當今的物理學體系或者化學體系存在崩塌的可能嗎?
複數的物理意義是什麼?
太陽的內部有可見光嗎?
鞅過程與馬爾科夫過程是什麼關係?

TAG:數學 | 物理學 | 冷知識 |