ケプラーの第3法則
質量 M の天体のまわりを円軌道(半径 r)で回る物体の公転周期を T とすると、万有引力定数 G を使って
r3T2=GM4π2右辺は中心の天体で決まる定数なので、周期の2乗は軌道半径の3乗に比例する。
この記事の要点
- 向心力の正体は万有引力。円運動の運動方程式 mrv2=Gr2Mm を立てる
- 速さを v=T2πr と書き直すと、v が消えて T と r の関係だけが残る
- 比例定数は GM4π2。周期と半径から、中心の天体の質量が求まる
ケプラーの第3法則は、「T2∝r3」と暗記されがちです。この記事では、万有引力の法則と円運動の運動方程式から出発して導き、比例定数の意味まで確かめます。使う定数は、万有引力定数 G=6.67×10−11N⋅m2/kg2 です。円運動の運動方程式は 力学的エネルギー保存則の導出 の発展例題でも使っています。
ケプラーの第3法則は何を言っているか
ケプラーは惑星の観測データから、惑星ごとに「公転周期の2乗」と「軌道の大きさの3乗」の比が同じになることを見つけ、1619年の『世界の調和』で発表しました。その約70年後の1687年、ニュートンが『プリンキピア』で、万有引力の法則から同じ関係が出ることを示しました。
太陽系の惑星で、横軸を r3(単位 AU³)、縦軸を T2(単位 年²)にとると、点は原点を通る傾き1の直線に乗ります。
惑星の点は、原点を通る直線の上に並ぶ(AU は地球の軌道半径、年は地球の公転周期)ケプラーの第3法則の導出
軌道を円とします。中心の天体(質量 M)は動かないものとし、まわりを回る物体の質量を m とします。
万有引力 F が中心を向いて、向心力のはたらきをする1. 運動方程式を立てる
物体に働く力は中心向きの万有引力 F=Gr2Mm だけです。一方、半径 r・速さ v の円運動の加速度は中心向きに rv2 です。運動方程式 ma=F より
mrv2=Gr2Mm2. 質量 m を消して v2 を出す
両辺の m を約分し、両辺に r を掛けます。
v2=rGMこれは、物体の質量によらず、半径 r で決まる速さです。
3. 速さを周期で書き直す
知りたいのは周期 T なので、v を消します。1周は長さ 2πr を時間 T で進むので v=T2πr です。これを2乗して代入します。
(T2πr)2T24π2r2=rGM=rGM4. T2 と r3 の比にまとめる
分母を払うため、まず両辺に rT2 を掛けます。
4π2r3=GMT2次に両辺を GMr3 で割ります。
r3T2=GM4π2右辺は G と中心の天体の質量 M だけで決まる定数です。したがって T2∝r3 で、比例定数は GM4π2 です
。物体の質量 m は式に残らないので、どの惑星でも同じ法則になります。
楕円軌道のときは
実際の惑星の軌道は楕円で、ケプラーの第3法則も「半長軸 a の3乗と周期の2乗」で書かれます。円は半長軸が半径に等しい楕円なので、上の結果で r を a に置き換えた a3T2=GM4π2 が成り立ちます。楕円の場合の証明は、面積速度の保存などを使うため、この記事では結果だけ紹介します。
楕円軌道では、半径のかわりに半長軸 a を使う使い方の例題
例題
地球は太陽のまわりを、半径 1.50×1011m の円軌道で、周期 3.16×107s(1年)で回っているとみなす。太陽の質量を求めよ。G=6.67×10−11N⋅m2/kg2 とする。
導いた式を M について解き(逆向きに使い)ます。
M=GT24π2r3=6.67×10−11×(3.16×107)24π2×(1.50×1011)3≒6.66×1041.33×1035≒2.0×1030 kg検算:地球の速さは v=T2πr≒3.0×104m/s。手順2の式 M=Gv2r に入れても 2.0×1030kg で一致します。単位は、N⋅m2/kg2⋅s2m3=N⋅s2kg2⋅m=kg(N=kg⋅m/s2 を使う)で質量になります。
ケプラーの第3法則の発展的な使い方
地上での重力加速度 g と地球の半径 R を使うと、G や地球の質量を知らなくても人工衛星の問題が解けます。地表の物体では、重力が万有引力に等しいので mg=GR2Mm、つまり GM=gR2 です。
発展例題
地球の自転と同じ周期 8.64×104s(24時間)で赤道上空を円運動し、地上から見て静止して見える静止衛星の、地表からの高さを求めよ。地球の半径を R=6.4×106m、地表の重力加速度を g=9.8m/s2 とする。
1. GM を g と R で置き換える
導いた式 r3T2=GM4π2 の GM に gR2 を入れ、r3 について解きます。
r3=4π2gR2T22. 数値を入れる
r3r=4π29.8×(6.4×106)2×(8.64×104)2≒7.59×1022 m3≒4.23×107 m3. 地表からの高さにする
r は地球の中心からの距離なので、半径 R を引きます。
h=r−R=4.23×107−0.64×107≒3.6×107 m検算(別の方法):地表すれすれの円軌道の周期は T0=2πgR≒5.08×103s です。T2∝r3 より Rr=(T0T)2/3≒6.62 なので、r≒4.23×107m で一致します。
この問題では、ケプラーの第3法則の GM を gR2 に置き換え、周期から軌道半径を逆算しました
。
注意
r は中心からの距離で、地表からの高さではありません。「高さ」を問われたら、最後に地球の半径を引くのを忘れないでください。
物理の参考書・問題集は テキスト購入(物理) で、レベル別に紹介しています。
この記事のまとめ
- 万有引力が向心力になる:mrv2=Gr2Mm から v2=rGM。
- v=T2πr を入れて整理すると r3T2=GM4π2。質量 m は消え、比例定数は中心の天体で決まる。
- 楕円軌道では r を半長軸 a に置き換える。周期と半径から中心天体の質量が、GM=gR2 と組み合わせると衛星の軌道が求まる。