ケプラーの法則とは

力学

ケプラーの法則とは

ケプラーの法則とは、17世紀初頭に天文学者のケプラーにより発見された、惑星の運動に関する法則です。ケプラーは、天文学者のブラーエの膨大な観測記録から、太陽に対する火星の運動を推定し、3つの法則としてまとめました。

  • 第1法則:楕円の法則
  • 第2法則:面積速度一定の法則
  • 第3法則:周期と長径の法則

そして、これらのケプラーの法則は、ニュートンの運動方程式から導くことができます。

運動方程式を解く

ここでは、惑星の運動法的式を導きます。まず、極座標のr成分でのニュートンの運動方程式は以下になります。

$$F_r=mA_r=m(\ddot{r}-r\dot{\theta}^2)$$

ページが見つかりませんでした |
数理の散策路

恒星と惑星の万有引力を以下で表されます。

$$F_r=-\frac{GmM}{r^2}$$

これをr方向の運動方程式に代入すると以下になります。ここで、$l=mr^2\dot{\theta}$ で角運動量を表します。

$$\frac{d^2r}{dt^2}+\frac{GM}{r^2}-\frac{l^2}{mr^3}=0  -①$$

運動方程式を解くと以下の解が得られます。この式の導出はこの記事の最後で行います。

$$r=\frac{\lambda}{1+\epsilon\cos{\theta}}  -②$$

$$\lambda\equiv\frac{l^2}{GMm^2}  , \epsilon\equiv\frac{Cl^2}{GMm^2}$$

ケプラーの法則を説明する

第1法則:楕円の法則

第1法則とは、惑星は恒星の周りを楕円軌道を描いて公転し、恒星は楕円の焦点に1つに位置することです。

惑星の運動方程式を解いて得られる以下の式は、惑星の軌道が楕円であることを表しています。

$$\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$$

$$a\equiv\frac{\lambda}{1-\epsilon^2} , b\equiv\frac{\lambda}{\sqrt{1-\epsilon^2}}$$

第2法則:面積速度一定の法則

第2法則とは、面積速度$\Delta$が一定となることで、角運動量の保存により導くことが得きます。 ここで$\phi$は速度Vのr方向とのなす角度です。

$$\Delta=\frac{1}{2}rV_\theta=\frac{1}{2}rV\sin{\phi}=\frac{l}{2m}$$

第3法則:周期と長径の法則

第3法則は、公転周期の2乗は長半径の3乗に比例するという法則です。公転周期は楕円の面積Sを面積速度で割ると得られる。

$$T=\frac{S}{\Delta}=\frac{\pi ab}{\Delta}=\frac{2\pi m\lambda^2}{l(1-\epsilon^2)^{3/2}}$$

$$T^2=\frac{4\pi^2}{GM}a^3$$

惑星の軌道を導く

以下で①から②の導出を行います。まず、①で時間微分を θ の微分に置換えます。$\dot{\theta}=l/mr^2$ であることから、

$$\frac{d}{dt}=\frac{d\theta}{dt}\frac{d}{d\theta}=\frac{l}{mr^2}\frac{d}{d\theta}$$

であるため、以下のように置き換えられます。

$$\frac{d^2r}{dt^2}=\frac{l^2}{m^2r^2}\frac{d}{d\theta}\left(\frac{1}{r^2}\frac{dr}{d\theta}\right)$$

次に $u=1/r$ と置くと、

$$\frac{1}{r^2}\frac{dr}{d\theta}=-\frac{du}{d\theta}$$

となるため、以下のように書き換えられます。

$$\frac{d^2r}{dt^2}=-\frac{l^2u^2}{m^2}\frac{d^2u}{d\theta^2}$$

これを①に代入すると、以下になります。

$$\frac{d^2u}{d\theta^2}=-u+\frac{GMm^2}{l^2}$$

Cと $\alpha$ を積分定数とすると、この微分方程式の解は以下になります。

$$u=C\cos{(\theta+\alpha)}+\frac{GMm^2}{l^2}$$

最後に、$r=1/u$、$\alpha=0$ とすると②になります。

タイトルとURLをコピーしました