ケプラーの第3法則の導出:万有引力と円運動から導く方法

物理力学公式の成立 更新
目次

ケプラーの第3法則

質量 MM の天体のまわりを円軌道(半径 rr)で回る物体の公転周期を TT とすると、万有引力定数 GG を使って

T2r3=4π2GM \frac{T^2}{r^3}=\frac{4\pi^2}{GM}

右辺は中心の天体で決まる定数なので、周期の2乗は軌道半径の3乗に比例する。

この記事の要点

  • 向心力の正体は万有引力。円運動の運動方程式 mv2r=GMmr2m\dfrac{v^2}{r}=G\dfrac{Mm}{r^2} を立てる
  • 速さを v=2πrTv=\dfrac{2\pi r}{T} と書き直すと、vv が消えて TT と rr の関係だけが残る
  • 比例定数は 4π2GM\dfrac{4\pi^2}{GM}。周期と半径から、中心の天体の質量が求まる

ケプラーの第3法則は、「T2∝r3T^2\propto r^3」と暗記されがちです。この記事では、万有引力の法則と円運動の運動方程式から出発して導き、比例定数の意味まで確かめます。使う定数は、万有引力定数 G=6.67×10−11 N⋅m2/kg2G=6.67\times10^{-11}\,\mathrm{N\cdot m^2/kg^2} です。円運動の運動方程式は 力学的エネルギー保存則の導出 の発展例題でも使っています。

ケプラーの第3法則は何を言っているか

ケプラーは惑星の観測データから、惑星ごとに「公転周期の2乗」と「軌道の大きさの3乗」の比が同じになることを見つけ、1619年の『世界の調和』で発表しました。その約70年後の1687年、ニュートンが『プリンキピア』で、万有引力の法則から同じ関係が出ることを示しました。

太陽系の惑星で、横軸を r3r^3(単位 AU³)、縦軸を T2T^2(単位 年²)にとると、点は原点を通る傾き1の直線に乗ります。

惑星の点は、原点を通る直線の上に並ぶ(AU は地球の軌道半径、年は地球の公転周期)

ケプラーの第3法則の導出

軌道を円とします。中心の天体(質量 MM)は動かないものとし、まわりを回る物体の質量を mm とします。

万有引力 F が中心を向いて、向心力のはたらきをする

1. 運動方程式を立てる

物体に働く力は中心向きの万有引力 F=GMmr2F=G\dfrac{Mm}{r^2} だけです。一方、半径 rr・速さ vv の円運動の加速度は中心向きに v2r\dfrac{v^2}{r} です。運動方程式 ma=Fma=F より

mv2r=GMmr2 m\frac{v^2}{r}=G\frac{Mm}{r^2}

2. 質量 mm を消して v2v^2 を出す

両辺の mm を約分し、両辺に rr を掛けます。

v2=GMr v^2=\frac{GM}{r}

これは、物体の質量によらず、半径 rr で決まる速さです。

3. 速さを周期で書き直す

知りたいのは周期 TT なので、vv を消します。1周は長さ 2πr2\pi r を時間 TT で進むので v=2πrTv=\dfrac{2\pi r}{T} です。これを2乗して代入します。

(2πrT)2=GMr4π2r2T2=GMr \begin{aligned} \left(\frac{2\pi r}{T}\right)^2&=\frac{GM}{r}\\[4pt] \frac{4\pi^2r^2}{T^2}&=\frac{GM}{r} \end{aligned}

4. T2T^2 と r3r^3 の比にまとめる

分母を払うため、まず両辺に r T2r\,T^2 を掛けます。

4π2r3=GM T2 4\pi^2r^3=GM\,T^2

次に両辺を GM r3GM\,r^3 で割ります。

T2r3=4π2GM \frac{T^2}{r^3}=\frac{4\pi^2}{GM}

右辺は GG と中心の天体の質量 MM だけで決まる定数です。したがって T2∝r3T^2\propto r^3 で、比例定数は 4π2GM\dfrac{4\pi^2}{GM} です 。物体の質量 mm は式に残らないので、どの惑星でも同じ法則になります。

楕円軌道のときは

実際の惑星の軌道は楕円で、ケプラーの第3法則も「半長軸 aa の3乗と周期の2乗」で書かれます。円は半長軸が半径に等しい楕円なので、上の結果で rr を aa に置き換えた T2a3=4π2GM\dfrac{T^2}{a^3}=\dfrac{4\pi^2}{GM} が成り立ちます。楕円の場合の証明は、面積速度の保存などを使うため、この記事では結果だけ紹介します。

楕円軌道では、半径のかわりに半長軸 a を使う

使い方の例題

例題

地球は太陽のまわりを、半径 1.50×1011 m1.50\times10^{11}\,\mathrm{m} の円軌道で、周期 3.16×107 s3.16\times10^{7}\,\mathrm{s}(1年)で回っているとみなす。太陽の質量を求めよ。G=6.67×10−11 N⋅m2/kg2G=6.67\times10^{-11}\,\mathrm{N\cdot m^2/kg^2} とする。

導いた式を MM について解き(逆向きに使い)ます。

M=4π2r3GT2=4π2×(1.50×1011)36.67×10−11×(3.16×107)2≒1.33×10356.66×104≒2.0×1030 kg \begin{aligned} M&=\frac{4\pi^2r^3}{GT^2}\\[4pt] &=\frac{4\pi^2\times\left(1.50\times10^{11}\right)^3}{6.67\times10^{-11}\times\left(3.16\times10^{7}\right)^2}\\[4pt] &\fallingdotseq\frac{1.33\times10^{35}}{6.66\times10^{4}}\fallingdotseq 2.0\times10^{30}\ \mathrm{kg} \end{aligned}

検算:地球の速さは v=2πrT≒3.0×104 m/sv=\dfrac{2\pi r}{T}\fallingdotseq3.0\times10^{4}\,\mathrm{m/s}。手順2の式 M=v2rGM=\dfrac{v^2r}{G} に入れても 2.0×1030 kg2.0\times10^{30}\,\mathrm{kg} で一致します。単位は、m3N⋅m2/kg2⋅s2=kg2⋅mN⋅s2=kg\dfrac{\mathrm{m^3}}{\mathrm{N\cdot m^2/kg^2\cdot s^2}}=\dfrac{\mathrm{kg^2\cdot m}}{\mathrm{N\cdot s^2}}=\mathrm{kg}(N=kg⋅m/s2\mathrm{N}=\mathrm{kg\cdot m/s^2} を使う)で質量になります。

ケプラーの第3法則の発展的な使い方

地上での重力加速度 gg と地球の半径 RR を使うと、GG や地球の質量を知らなくても人工衛星の問題が解けます。地表の物体では、重力が万有引力に等しいので mg=GMmR2mg=G\dfrac{Mm}{R^2}、つまり GM=gR2GM=gR^2 です。

発展例題

地球の自転と同じ周期 8.64×104 s8.64\times10^{4}\,\mathrm{s}(24時間)で赤道上空を円運動し、地上から見て静止して見える静止衛星の、地表からの高さを求めよ。地球の半径を R=6.4×106 mR=6.4\times10^{6}\,\mathrm{m}、地表の重力加速度を g=9.8 m/s2g=9.8\,\mathrm{m/s^2} とする。

1. GMGM を gg と RR で置き換える

導いた式 T2r3=4π2GM\dfrac{T^2}{r^3}=\dfrac{4\pi^2}{GM} の GMGM に gR2gR^2 を入れ、r3r^3 について解きます。

r3=gR2T24π2 r^3=\frac{gR^2T^2}{4\pi^2}

2. 数値を入れる

r3=9.8×(6.4×106)2×(8.64×104)24π2≒7.59×1022 m3r≒4.23×107 m \begin{aligned} r^3&=\frac{9.8\times\left(6.4\times10^{6}\right)^2\times\left(8.64\times10^{4}\right)^2}{4\pi^2}\\[4pt] &\fallingdotseq7.59\times10^{22}\ \mathrm{m^3}\\[4pt] r&\fallingdotseq4.23\times10^{7}\ \mathrm{m} \end{aligned}

3. 地表からの高さにする

rr は地球の中心からの距離なので、半径 RR を引きます。

h=r−R=4.23×107−0.64×107≒3.6×107 m h=r-R=4.23\times10^{7}-0.64\times10^{7}\fallingdotseq3.6\times10^{7}\ \mathrm{m}

検算(別の方法):地表すれすれの円軌道の周期は T0=2πRg≒5.08×103 sT_0=2\pi\sqrt{\dfrac{R}{g}}\fallingdotseq5.08\times10^{3}\,\mathrm{s} です。T2∝r3T^2\propto r^3 より rR=(TT0)2/3≒6.62\dfrac{r}{R}=\left(\dfrac{T}{T_0}\right)^{2/3}\fallingdotseq6.62 なので、r≒4.23×107 mr\fallingdotseq4.23\times10^{7}\,\mathrm{m} で一致します。

この問題では、ケプラーの第3法則の GMGM を gR2gR^2 に置き換え、周期から軌道半径を逆算しました 。

注意

rr は中心からの距離で、地表からの高さではありません。「高さ」を問われたら、最後に地球の半径を引くのを忘れないでください。

物理の参考書・問題集は テキスト購入(物理) で、レベル別に紹介しています。

この記事のまとめ

  • 万有引力が向心力になる:mv2r=GMmr2m\dfrac{v^2}{r}=G\dfrac{Mm}{r^2} から v2=GMrv^2=\dfrac{GM}{r}。
  • v=2πrTv=\dfrac{2\pi r}{T} を入れて整理すると T2r3=4π2GM\dfrac{T^2}{r^3}=\dfrac{4\pi^2}{GM}。質量 mm は消え、比例定数は中心の天体で決まる。
  • 楕円軌道では rr を半長軸 aa に置き換える。周期と半径から中心天体の質量が、GM=gR2GM=gR^2 と組み合わせると衛星の軌道が求まる。