MENU

【古典力学】第12講 ケプラー問題1

  • このエントリーをはてなブックマークに追加

第12講の導入

今回は万有引力の下での惑星の運動、いわゆるケプラー問題について考えてみましょう。歴史的には、ニュートンはケプラーの法則などを基にしてニュートンの万有引力とニュートンの運動法則を見出しました。ここでは、その逆を辿り、ニュートンの万有引力の法則とニュートンの運動法則からケプラーの法則が導くことを目指してお話を進めて行きます。

中心力問題の運動は平面内に限られるということ

質量mの惑星(例えば地球など)が質量Mの太陽の周りを回っている状況を考えてみましょう。ニュートンの万有引力の法則によれば、惑星は太陽に向かう方向に、以下の式で与えられるような距離の2乗に反比例した力Fを受けています。
F=mHr2


但し、HH=GMという定数です。今、太陽の質量は惑星の質量よりも十分大きい(Mm)ということを考えると、太陽は近似的に動かないとみなすことが出来ます。そこで、太陽を座標系の原点を取れば、この力は中心力であるとみなせるため、原点のまわりの角運動量が保存され、惑星の運動は質点の速度ベクトルと原点を含むような平面内で留まり続けることになります。このことは、第07講で扱いましたので、忘れてしまっている人はそこの項目を復習してください。以下ではそのような平面内に限定して惑星運動を論じることにします。

軌道の方程式を解くこと

さて、惑星に生じる加速度を平面極座標成分に分けて
a=arer+aθeθ


とあらわすと、引力が中心力であることからθ成分の加速度は0となる必要があります。すなわち、
aθ=r¨θ+2˙r˙θ=1rddt(r2˙θ)=0

となります。ここで角運動量L=mr2˙θは変化しないので、Lを定数であると考えて、
˙θ=Lmr2

と書くことが出来ます。これはケプラーの第2法則、すなわち面積速度一定の法則に対応しています。これを動径方向の運動方程式に代入すると、
Hr2=ar=¨rr˙θ2=d2rdt2L2m2r3

となります。これは距離rの時間tに関する常微分方程式ですから、これを解くとその時間変化を追うことが出来ます。しかし、ここでは距離の時間依存性を直接見るのではなく、軌道の方程式を出したいので、距離rを時間tではなく角度θについての常微分方程式に書き直すことを考えることになります。

まず、時間に関する1階微分は合成関数の微分の関係式を利用して
˙r=drdt=dθdtdrdθ=Lmr2drdθ=Lmddθ(1r)


と書くことが出来ます。同様に、時間に関する2階微分は
¨r=ddt˙r=dθdtddθ˙r=Lmr2ddθ˙r=L2m2r2ddθddθ(1r)

となります。これを先ほどの常微分方程式に代入すると、
L2m2{1r2d2dθ2(1r)+1r3}=H1r2

が得られます。これを解きやすくするために、両辺に(mr/L)2を掛けてu=1/rとおくと、
d2udθ2+u=Hm2L2

となります。これで、太陽からの距離rの時間に対する依存性をあらわす常微分方程式は、距離の逆数uの角度に対する依存性をあらわす常微分方程式へと書き換えられました。これを解くと、
u=Hm2L2+Acos(θ+δ)=Hm2L2+Acosθ

となります。但し、途中で、初期位相δは座標系を見ているだけの文字なので、δ=0としても一般性を失わないとして0と採りました。これをrの戻すと、
r=1Hm2/L2+Acosθ=L2/(Hm2)1+AL2/(Hm2)cosθ

と書くことが出来ます。

第12講のまとめ

最後に一言コメントをして第12講を終わりにしましょう。ここで求めた軌道の式は|AL2/(Hm2)|の値によって異なる曲線を表すことになります。これがケプラーの第1法則と関連することになります。特に軌道を楕円軌道と限定することでケプラーの第3法則も導くことが出来ます。これらの解説を次回の記事で行いたいと思います。

  • このエントリーをはてなブックマークに追加

SNSでもご購読できます。