地球を貫通する穴を通る物体の運動
はじめに
この記事を書くことになったきっかけは先日放送していた頭脳王というテレビ番組である。その中で地球を貫通する穴を通る物体の運動についての問題が出た。このような条件の場合、結論から言うと物体は自由落下せずに単振動をするという興味深い結果が出てくるのだが、この番組はやっぱりお約束の条件写す時間Δt&解説量εだったので、*1物理的に興味深い背景のある内容なのに説明もしないままほったらかしにするのはもったいないと思ったので勝手に説明しようと思った次第である。*2
「高度な物理力」を手に入れたい人は、少なくともこの記事の内容は理解できる必要があるだろう。
流れは以下の通り
レベルとしては1,2が大学初年度でやる力学のレベルであり、3,4は高校レベルの物理だ。*3地球を貫通する穴はたまに入試問題のテーマとなるが、そのようなときは1,2の内容は問題文に前提として与えられている。(ちなみに”数学について解説するブログなのに物理を解説するのはどういうことだ”というツッコミは気にしない方向性で行きます。積分計算とか実質数学でしょう(適当))
制約条件
まず以下の条件を仮定する。
- 地球は密度が場所によらず一定値である。
- 地球は半径がの球であり、質量はとおく。
- つまりである。
- 地球の中心を原点とする座標を取る。
- つまり地球はという領域で表される
- 穴は軸である。つまり上を物体は運動する
- 物体は質点とみなすため穴の大きさについては考えない。物体の質量はとりあえずとおく。ここではに比べて十分に小さいものとする。*4
- 万有引力定数をとおく。
また、空間に対して球面座標を取る。
変換公式としては
$$\left\{\begin{array}{ccl}x&=&r\sin{\theta}\cos{\phi}\\y&=&r\sin{\theta}\sin{\phi}\\z&=&r\cos{\theta}\end{array} \right.$$
である。このとき、となる。
1.地球内部の重力場1 - 球殻周辺の重力
いきなり球面から考えるのは難しいので、半径がまでの球殻から考える。
球殻上のが示す微小領域が物体に及ぼす万有引力を考える。
微小領域は各辺がの直方体に近似できるため、体積はである。よってこの領域の質量は密度をかければいいので
$$\Delta M= \rho r_0^2\sin{\theta_0}\Delta r\Delta \theta \Delta \phi$$
となる。微小領域がある位置は
であり、物体がある位置は
である。
$$\Delta F=\frac{G\cdot \Delta M\cdot m\cdot (\vec{r_1}-\vec{r_2})}{|\vec{r_1}-\vec{r_2}|^3}$$
という式で表される。
ここで、
$$|\vec{r_1}-\vec{r_2}|=\sqrt{|\vec{r_1}|^2+|\vec{r_2}|^2-2\vec{r_1}\cdot\vec{r_2}}$$
$$=\sqrt{r_0^2+z^2-2zr_0\cos{\theta_0}}$$
であることに注意して、球殻全体での万有引力の値を求めると、これはで積分すればよいので、
$$\int_{0}^{\pi}d\theta \int_{0}^{2\pi}d\phi \frac{G \rho mr_0^2\Delta r\sin{\theta}}{\sqrt{z^2+r_0^2-2zr_0\cos{\theta}}^3}(r_0\sin{\theta}\cos{\phi},r_0\sin{\theta}\sin{\phi},r_0\cos{\theta}-z)$$
ここで、力の成分に着目すると、
$$\int_{0}^{\pi}d\theta \int_{0}^{2\pi}d\phi \frac{G \rho mr_0^2\Delta r\sin{\theta}}{\sqrt{z^2+r_0^2-2zr_0\cos{\theta}}^3}r_0\sin{\theta}\cos{\phi}$$
,ごとに分けると
$$\int_{0}^{\pi}\frac{G \rho mr_0^3\Delta r\sin^2{\theta}}{\sqrt{z^2+r_0^2-2zr_0\cos{\theta}}^3}d\theta\cdot\int_{0}^{2\pi}\cos{\phi}d\phi$$
となるが、であるため、成分全体の値は0となる。成分についても同様にして、であるため、成分の値は0となる。よって球殻から受ける万有引力の向きは軸と平行である。
座標についての積分は
$$\int_{0}^{\pi}d\theta \int_{0}^{2\pi}d\phi \frac{G \rho mr_0^2\Delta r\sin{\theta}(r_0\cos{\theta}-z)}{\sqrt{z^2+r_0^2-2zr_0\cos{\theta}}^3}$$
積分の中身はには依存しないため、先にから積分してしまおう。そして定数を外に移動させると以下のようになる。
$$2\pi G \rho mr_0^2\Delta r\int_{0}^{\pi}\frac{\sin{\theta}(r_0\cos{\theta}-z)}{\sqrt{z^2+r_0^2-2zr_0\cos{\theta}}^3}d\theta$$
この積分については、と置換すればよい。するととなり、ならばとなるため、
$$ -2\pi G \rho mr_0^2\Delta r\int_{-1}^{1}\frac{r_0s+z}{\sqrt{z^2+r_0^2+2zr_{0}s}^3}ds$$
簡単のためとでもする。このとき積分
$$\int_{-1}^{1}(r_0s+z)(A+Bs)^{-3/2}ds$$
$$\frac{1}{Bz}\left[(A+Bs)^{1/2}-(z^2-r_0^2)(A+Bs)^{-1/2} \right]_{-1}^{1}$$
となり、であることから、
$$\frac{1}{2z^2r_0}\left\{|r_0+z|-|r_0-z|+(z^2-r_0^2)\left(\frac{1}{|r_0-z|}-\frac{1}{|r_0+z|}\right) \right\}$$
となる。
絶対値が多いのでの値で場合分けをする。
(i)
このとき、より、となる。これを代入してもともとあった定数を掛けると
$$\Delta F_z=4\pi G \rho m r_0^2\Delta r\frac{1}{z^2}$$
となる。
(ii)
このときとなるため、絶対値のところに代入をすると
きれいに0となる。よってこのときとなる。
(iii)
このとき、となるため、代入して計算すると
$$\Delta F_z=-4\pi G \rho m r_0^2\Delta r\frac{1}{z^2}$$
となる。((i)とは符号が逆になる。なぜならが一斉に符号が変換されて全体が-1倍となるからである。)
2.地球内部の重力場2 - 球殻の結果を積分
$$\Delta F(z,r_0)=\left\{\begin{array}{cc}4\pi G \rho mr_0^2\frac{1}{z^2}\Delta r&z\leq -r_0\\0&-r_0\leq z\leq r_0\\-4\pi G \rho mr_0^2\frac{1}{z^2}\Delta r&r_0\leq z\end{array} \right.$$
地球内部での万有引力を考えているので、を仮定する。
これを球殻たちを積分して球にすることを考える。
このとき、球全体で考えたときの万有引力は
$$F(z)=\int_{0}^{R}\Delta F(z,r)=\int_{0}^{R}\frac{\Delta F(z,r)}{\Delta r}dr$$
となる。がよりも内側である場合のみを考えれば良いので
$$\int_{0}^{z}\frac{\Delta F(z,r)}{\Delta r}dr$$
となる。これについて、に先ほど求めたものを当てはめると
$$\int_{0}^{z}-4\pi G \rho mr^2\frac{1}{z^2}dr$$
$$F(z)=-\frac{4}{3}\pi G \rho mz$$
となる。
これに対して、とおいて整理すると(は重力加速度であり、約9.8m/s^2となる。)
$$F(z)=-mg\frac{z}{R}$$
となる。以降は座標を横に回転させて、
$$F(x)=-mg\frac{x}{R}$$と書くことにする。
3.物体の運動方程式を立てる
この系において、物体に働く力は万有引力だけを考えればよいだろう。
つまり、
$$ma=m\frac{d^2x(t)}{dt^2}=-mg\frac{x}{R}$$
となる。物体の移動する領域は一次元の直線であり、物体に働く力は位置に比例している。この時点で†高度な物理力†を持っている人々はあたかも†当たり前のように†単振動であるとわかる。だが高度な物理力を有していていようがいまいが、任意の人々にこの記事を見てほしいので微分方程式の解き方は次章に書く。
4.運動方程式を解く
この章の目的はをの関数で陽に表すことである。ここでの条件を整理してみる。
- を満たす。
- を満たす。(のときから地球表面へ自由落下する)
こんなところだろう。以降は簡単のため、とおき、微分はのようにかく。
ここで、の微分方程式はとなる。
ここで力学的エネルギーの保存量を出して連立微分方程式を解く方法もあるが、あえて特性方程式を使った解き方をしてみよう。
つまりこれの一般解は任意の定数を用いることで、
$$x(t)=C_1e^{i\omega t}+C_2e^{-i\omega t}$$
となる。これをオイラーの公式などを使って変形すると、
$$x(t)=(C_1+C_2)\cos{(\omega t)}+i(C_1-C_2)\sin{(\omega t)}$$
となる。
はもちろん実数でなくてはいけないので(が実数とは限らないことに注意!)とでもすると、(は実数)このとき
$$C_1=\frac{A-iB}{2},C_2=\frac{A+iB}{2}$$
とすれば、
$$x(t)=A\cos{(\omega t)}+B\sin{(\omega t)}$$
となる。あとは制約条件からを具体的に求めてみる。
であることから、となる。
よって、結局特殊解として、
$$x(t)=R\cos{\left(\sqrt{\frac{g}{R}}t \right)}$$
という公式が出てくるのである。これこそ番組の例の問題で出てきた制約条件に出てくる公式である。
あとがき
というように、テレビ画面に数秒しか映っていなかった条件にぱっと出てきた公式も、長い道のりをたどってきた結果ようやく求めることができるものである。公式というものは天から降ってきたものではなく、*7人間が頑張って計算した結果導くことができるものだ。たまには「意味不明な問題が解けるわけがない」と思考停止するのをやめて、問題自体の考察をするのも面白いだろう。もしかしたらそこに高度な物理力を有する者の努力が垣間見えるのかもしれない。
関連記事(2020年3月追記)
頭脳王2020についても記事を作成したので良かったら見てください
記事概要:「東京スタジアムから凱旋門までボールを投げる問題について、地球の丸みを考慮したらどうなるか考察してみた」
*1:測定したところΔtは4秒くらいだった。またεについては河野玄斗とかいう天才の基準にとって簡単だということくらいであった。実際(この記事で導出する)単振動の式が条件に書いてあったし、特定のsinの値が明記してあったのでやることは実質算数の計算だけだった。一回簡単な微分をするだけだが算数の計算の重さに比べたら微々たるものである。この記事を見た人で頭脳王に出るという場合、物理は高校レベルの知識しか要求されないので、筆算の高速化訓練をしたり、ルートを開平法を使うことで計算する訓練をするといいだろう。(筆者は頭脳王に出たことがないので適当なことを言っているだけかもしれないが)
*2:頭脳王自体が番組に問題の解説が少ないという気がします。これを見た人々が問題の解説を自主的にネットに垂れ流すような流れになるのも面白いでしょう。
*3:駿台予備校の物理の参考書「新物理入門」にもこの記事と同様なことが書いてあったのでもしかしたら高校生でも1,2が理解できる人はいるかもしれない
*4:地球が物体から受ける万有引力が無視できるくらい小さいものとする
*5:ベクトルの形式なので高校物理的な視点だと異なるように見える。このベクトルの絶対値を見ればきちんと距離の二乗に反比例している万有引力の式になる。これは元の式に加えて、万有引力は内側の向きに伝わるということを考慮して、距離の逆二乗の式に内側向きの単位ベクトルを掛けたものと考えればよいだろう。
*6:ちなみにz座標の位置を地球外部に置くと見慣れた万有引力の式になる。また、球殻自体の質量はであるため、球殻の内側では各方向からの重力は相殺されるが、外側だと球殻の質量分の質点が点にある場合の万有引力と等しくなる。