物理/惑星の運動(2)

提供: Internet Web School

(版間での差分)
( ケプラーの第2法則の導出)
( 微分方程式を解いて、惑星の軌道形を求める  )
 
(間の34版分が非表示)
7 行: 7 行:
運動の第1、第2法則のもとで、<br/>
運動の第1、第2法則のもとで、<br/>
万有引力の法則と惑星に関するケプラーの3法則が同等であることを導く。
万有引力の法則と惑星に関するケプラーの3法則が同等であることを導く。
-
== 惑星の運動方程式   ==
+
== ニュートンの運動の第一、第二法則による惑星の運動方程式の導出   ==
-
他の惑星の影響を無視して、太陽と一つの惑星の二体からなる系を考える。<br/>
+
他の惑星の影響を無視して、太陽と一つの惑星の二体からなる系をある慣性系から観測する。(注1と注2を参照のこと)<br/>
太陽の質量を$M$、惑星の質量を$m$とする。<br/>
太陽の質量を$M$、惑星の質量を$m$とする。<br/>
この二体を、ともに質点であるとみなし、位置はそれぞれの重心とする。 <br/>
この二体を、ともに質点であるとみなし、位置はそれぞれの重心とする。 <br/>
-
この系には、系外から力が作用しないため、任意の慣性系から観測すると(注1参照)、<br/>
+
この系には、系外から力が作用しないため、、<br/>
[[物理/質点の運動と質点系#質点系の重心とその運動|系の重心は等速直線運動をする。]] <br/>
[[物理/質点の運動と質点系#質点系の重心とその運動|系の重心は等速直線運動をする。]] <br/>
-
ところが、太陽の質量は惑星の質量より圧倒的大きいため二体の重心は太陽の重心とほぼ一致し、 <br/>
+
従って、この重心を原点とする回転しない座標系も慣性系となる。<br/>
-
太陽重心は慣性系からみると等速直線運動しているとみなしてよい。<br/>
+
今後はこの慣性系で観測する。<br/><br/>
-
*[[wikipedia_ja:大きさ順の太陽系天体の一覧#一覧 |ウィキペディア(大きさ順の太陽系天体の一覧)]]
+
-
<br/>
+
(注1)慣性系の存在はニュートン力学(古典力学)の大前提であり、運動法則の第一法則としてまとめられている。<br/>
(注1)慣性系の存在はニュートン力学(古典力学)の大前提であり、運動法則の第一法則としてまとめられている。<br/>
 +
太陽系の重心に原点を取り、恒星系に対して回転しない座標系は慣性系に近いことが現在分かっている。
 +
*[[wikipedia_ja:慣性系 |ウィキペディア(慣性系)]]
慣性系からみた質点系の運動は、[[物理/質点の運動と質点系#質点系の運動|テキストの2.3]] で、その導入部分を説明している。
慣性系からみた質点系の運動は、[[物理/質点の運動と質点系#質点系の運動|テキストの2.3]] で、その導入部分を説明している。
<br/>
<br/>
(注2)他の惑星の影響を無視し、太陽と当該惑星の2体を質点とみなしても、<br/>
(注2)他の惑星の影響を無視し、太陽と当該惑星の2体を質点とみなしても、<br/>
他の惑星からの影響の少ない地球などの惑星の運動は、かなり正確に予測できることが知られている。<br/>
他の惑星からの影響の少ない地球などの惑星の運動は、かなり正確に予測できることが知られている。<br/>
-
正確にするには二体ではなく他の惑星も考慮した多体問題として論じなければならないが、解析は非常に難しくなる。<br/><br/>
+
正確にするには二体ではなく他の惑星も考慮した多体問題として論じなければならないが、解析は非常に難しくなる。
 +
=== 2体の重心は太陽の重心と一致するという近似を行う  ===
 +
 太陽の質量は惑星の質量より圧倒的大きいため二体の重心は太陽の重心とほぼ一致する。 <br/>
 +
*[[wikipedia_ja:大きさ順の太陽系天体の一覧#一覧 |ウィキペディア(大きさ順の太陽系天体の一覧)]]
 +
<br/>
 +
そこで議論を簡単にするため、このテキストでは2体の重心は太陽の重心と一致し、太陽は不動であると仮定する。<br/>
 +
興味のある読者は、この仮定をしないで、惑星と太陽の運動を解析してみてほしい。<br/>
 +
議論は殆ど同じように進められます。<br/><br/>
 +
 
そこで太陽$S$を原点$O$とする慣性座標系を考える。<br/>
そこで太陽$S$を原点$O$とする慣性座標系を考える。<br/>
任意の時刻$\ t\ $の惑星$P_t$の位置ベクトル$\vec{OP_t}$を$\vec{r}(t)$とかく。<br/>
任意の時刻$\ t\ $の惑星$P_t$の位置ベクトル$\vec{OP_t}$を$\vec{r}(t)$とかく。<br/>
28 行: 36 行:
惑星の運動は、ニュートンの第2法則から次の微分方程式で形式的には記述されが、<br/>力が未知なので、この方程式は解を持ち得ない。<br/>
惑星の運動は、ニュートンの第2法則から次の微分方程式で形式的には記述されが、<br/>力が未知なので、この方程式は解を持ち得ない。<br/>
$\ M\frac{d^{2} \vec{r}}{dt^2}(t)=\vec{f} $<br/>
$\ M\frac{d^{2} \vec{r}}{dt^2}(t)=\vec{f} $<br/>
-
== 運動の第1,第2法則と万有引力の法則からケプラーの3法則を導く==
 
-
=== 惑星の運動方程式  ===
 
-
万有引力を仮定すると、位置$\vec{r}$にある惑星に働く力$\vec f$は、<br/>
 
-
$\vec f = -\frac{GMm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}$<br/>となるので、惑星の運動方程式は次式で与えられる。<br/>
 
-
$\ M\frac{d^{2} \vec{r}}{dt^2}(t)=-\frac{GMm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} $<br/>
 
-
両辺をMで除すると<br/>
 
-
$\ \frac{d^{2} \vec{r}}{dt^2}(t)=-\frac{Gm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} \qquad \qquad (1)$<br/>
 
-
ここで、式(1)の右辺中の$G$は万有引力定数である。<br/>
 
-
式(1)を満たす(ベクトル値)関数$\vec{r}$を求めることを、微分方程式(1)を解くという。<br/>
 
-
なお、ある時刻$t_0$(以後、初期時刻と呼ぶ)の位置ベクトルと速度ベクトルが指定されたベクトルになるような<br/>
 
-
微分方程式の解を求めなければならないこともある。<br/>
 
-
指定された位置ベクトルを$\vec{r}_0$、速度ベクトルを$\vec{v}_0$と書くと、<br/>
 
-
速度ベクトル$\vec{v}(t)$は<br/>
 
-
$\vec{v}(t)=\frac{d \vec{r}}{dt}(t)$なので、<br/>
 
-
これ等の条件(初期条件と呼ぶ)は次のように表せる。<br/>
 
-
$\vec{r}(t_0)=\vec{r}_0、\quad \frac{d \vec r}{dt}(t_0)=\vec{v}_0  \qquad \qquad (2)$<br/>
 
-
初期条件をみたす微分方程式の解を求める問題を微分方程式の初期値問題と呼ぶ。<br/>
 
 +
== 運動の第1、第2法則と万有引力の法則からケプラーの3法則を導く ==
=== 惑星は平面上を運動する  ===
=== 惑星は平面上を運動する  ===
不動の太陽の位置$O$,時刻tの惑星の位置$P(t_0)$と$P(t_0)$を始点とした初速ベクトル$\vec{v}_0$を含む<br/>
不動の太陽の位置$O$,時刻tの惑星の位置$P(t_0)$と$P(t_0)$を始点とした初速ベクトル$\vec{v}_0$を含む<br/>
109 行: 101 行:
'''角運動量の保存法則'''<br/>
'''角運動量の保存法則'''<br/>
求心力、ないし遠心力(力の向きが原点向きか、その逆向き)だけを受けて運動する質点の角運動量は保存される。<br/><br/>
求心力、ないし遠心力(力の向きが原点向きか、その逆向き)だけを受けて運動する質点の角運動量は保存される。<br/><br/>
-
==== ケプラーの第2法則の導出  ====
+
=== ケプラーの第2法則の導出  ===
上述したことから惑星は平面上を動くので、<br/>
上述したことから惑星は平面上を動くので、<br/>
位置ベクトルと運動量ベクトルが、xy平面に載る様な正規直交右手系xyzを選べる。<br/>
位置ベクトルと運動量ベクトルが、xy平面に載る様な正規直交右手系xyzを選べる。<br/>
116 行: 108 行:
$=\bigl(0,0,r_{1}(t)p_{2}(t)-r_{2}(t)p_{1}(t)\bigr)$<br/>
$=\bigl(0,0,r_{1}(t)p_{2}(t)-r_{2}(t)p_{1}(t)\bigr)$<br/>
-
他方、ベクトル$\vec{r}(t)$ と ベクトル$\vec{p}(t)$のなす角を$\theta(t)$とすると<br/>
+
他方、ベクトル$\vec{r}(t)$ と ベクトル$\vec{p}(t)$のなす劣角(180度以下のほうの角)を$\theta(t)$とすると<br/>
この2つのベクトルと直交するz方向の$\vec{r}(t)\times \vec{p}(t) $の成分は<br/>
この2つのベクトルと直交するz方向の$\vec{r}(t)\times \vec{p}(t) $の成分は<br/>
-
ベクトル積の性質から、<br/>
+
[[物理/剛体と回転力#ベクトル積の定義|ベクトル積の定義]]から、<br/>
$\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t)$<br/>
$\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t)$<br/>
なので、<br/>
なので、<br/>
-
$\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t)  \qquad \qquad \qquad (1)$<br/>
+
$\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t)  \qquad \qquad \qquad (6)$<br/>
が、時不変であることが分かる。<br/>
が、時不変であることが分かる。<br/>
惑星質量は一定なので
惑星質量は一定なので
-
式(1)から<br/>
+
式(6)から<br/>
-
$\|\vec{r}(t)\| \times \| \vec{v}(t) \| \sin \theta(t)/2  \qquad \qquad \qquad (2)$<br/>
+
$\|\vec{r}(t)\| \times \| \vec{v}(t) \| \sin \theta(t)/2  \qquad \qquad \qquad (7)$<br/>
が、時不変であることが分かる。<br/>
が、時不変であることが分かる。<br/>
-
式(2)は動径$OP(t)$が時刻tの瞬間に単位時間あたりに掃く面積を表すので、惑星の面積速度と呼ばれる。
+
式(7)は動径$OP(t)$が時刻tの瞬間に単位時間あたりに掃く面積を表すので、惑星の面積速度と呼ばれる。
*[[File:GENPHY00010211-02.pdf|right|frame|図 時刻tから、微小時間に動径の掃く面積]]
*[[File:GENPHY00010211-02.pdf|right|frame|図 時刻tから、微小時間に動径の掃く面積]]
<br/>
<br/>
135 行: 127 行:
=== ケプラーの第1法則の導出  ===
=== ケプラーの第1法則の導出  ===
-
 
+
[[物理/惑星の運動(2)#惑星の運動方程式|惑星の運動方程式]](1)を解いて惑星の軌道を求めよう。<br/>
 +
求心力の下での運動を扱うには、極座標が都合がよい。<br/>
==== 惑星運動の微分方程式の極座標表示  ====
==== 惑星運動の微分方程式の極座標表示  ====
-
==== 微分方程式を解いて、惑星軌道の極座標表示を求める  ====
+
そこで惑星の位置ベクトルを、x軸の正の部分のつくる半直線を角度の基準とする極座標で<br/>
-
==== 極座標表示の惑星軌道を、直交座標表示に変換  ====
+
$\vec{r}(t)=r(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$<br/>
 +
と表す。<br/>
 +
この式を時間tで微分すると<br/>
 +
$\dot{\vec{r}}(t)={\dot r}(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
 +
+r(t){\dot \theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$<br/>
 +
$\qquad $但し、$\dot{\vec{r}}(t)\triangleq \frac{d\vec{r}}{dt}(t)$<br/>
 +
さらにtで微分すると<br/>
 +
$\ddot{\vec{r}}(t)$<br/>
 +
$=\ddot {r}(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr)
 +
+\dot{r}(t)\dot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$<br/>
 +
$+\dot{r}(t)\dot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$<br/>
 +
$+r(t)\ddot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)
 +
-r(t)\dot{\theta}^{2}(t)\bigl(\cos \theta (t),\sin \theta (t) \bigr) $<br/>
 +
$=\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$<br/>
 +
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$<br/>
 +
故に<br/>
 +
$\ddot{\vec{r}}(t)$<br/>
 +
$=\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$<br/>
 +
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr) \qquad \qquad \qquad (8)$<br/>
 +
が得られた。<br/>
 +
惑星に働く万有引力を、この極座標で表示すると<br/>
 +
$-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}
 +
=-\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr) \qquad (9)$<br/>
 +
式(8),(9)を惑星の運動方程式<br/>
 +
$\ddot{\vec r}(t)=-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}$<br/>
 +
に代入すると<br/>
 +
$\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$<br/>
 +
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$<br/>
 +
$= -\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr)$<br/>
 +
右辺を左辺に移行すると<br/>
 +
$\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$<br/>
 +
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)\equiv 0$<br/>
 +
$\bigl(\cos \theta (t),\sin \theta (t)\bigr)$と$\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$は直交ベクトルなので、この式から<br/><br/>
 +
$\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\equiv 0 \qquad \qquad (10)  $<br/>
 +
$2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \equiv 0 \qquad \qquad (11) $<br/>
 +
が得られる。<br/><br/>
 +
命題;惑星運動の微分方程式の極座標表示<br/>
 +
太陽の位置Oを原点、x軸の正の部分の半直線を角の基準線とする極座標で
 +
惑星の運動方程式を表現すると<br/>
 +
$\ddot {r}-r\dot{\theta}^{2}=-\frac{GM}{r^{2}}\qquad \qquad (12)  $<br/>
 +
$r^2\dot{\theta} = h \quad (任意の定数) \qquad \qquad (13) $<br/>
 +
となる。(注2 を参照のこと)<br/>
 +
証明<br/>
 +
式(12)は式(10)から明らか。<br/>
 +
式(13)の左辺をtで微分すると
 +
$\frac{d r^2\dot{\theta}}{dt}=r^2\ddot{\theta} + 2r\dot{r}\dot{\theta}
 +
=r\bigl(r\ddot{\theta} + 2\dot{r}\dot{\theta}\bigr)$<br/>
 +
$\qquad $式(11)から<br/>
 +
$= 0$<br/>
 +
左辺の時間微分が零なので、左辺は任意の定数hに等しい。<br/>
 +
証明終わり。 $\qquad \qquad \qquad \qquad \Box$<br/><br/>
 +
(注)式(13)は面積速度が一定であることを表現している。<br/>
 +
先に、このことは証明してあるが、その別証明である。<br/>
 +
(注2)関数の変数tまで書くと、式(12),(13)は、次のように書ける。<br/>
 +
$\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)=-\frac{GM}{r^{2}(t)}\qquad \qquad (12')  $<br/>
 +
$r^2(t)\dot{\theta}(t) = h \quad (任意の定数) \qquad \qquad (13') $<br/>
-
=== ケプラーの第3法則の導出  ===
+
==== 微分方程式を解いて、惑星の軌道形を求める  ====
 +
連立微分方程式(12)、(13)を解いて、$r=r(t),\theta=\theta(t)$を求めれば、<br/>
 +
任意の時刻tにおける惑星の位置がわかる。<br/>
 +
軌道の形は、この2式からtを消去して、rとθの関数関係<br/>
 +
$f(r,\theta)=0$, または$r=f(\theta)$<br/>
 +
を求めれな良い。<br/>
 +
しかし、残念なことにこれは大変困難である。<br/><br/>
-
== 運動の第2法則とケプラーの3法則から万有引力の法則を導く ==
+
軌道の形だけを知りたいときには、<br/>
 +
これを回避して簡単な(しかし技巧的な)直接的な解法がある。<br/>
 +
それは、連立微分方程式の解が<br/>
 +
$r=f\circ \theta ,\qquad$変数まで書けば、$r(t)=f\circ \theta (t)\triangleq f(\theta (t))$  <br/>
 +
という合成関数になることを仮定して、<br/>
 +
$r=f(\theta)$を直接求める問題に変換し、解を求める方法である。<br/><br/>
 +
===== 軌道の形を直接求めるための微分方程式の導出  =====
 +
$\frac{d r(t) }{dt}=\frac{d f\circ \theta (t) }{dt}$<br/>
 +
$\qquad $ 合成関数の微分公式を用いて、<br/>
 +
$=\frac{d f}{d\theta}(\theta (t))\dot{\theta}(t)$<br/>
 +
 
 +
$\qquad $式(13)を$\dot{\theta}=\frac{h}{r^2}$に変形し、上式に代入すると<br/>
 +
$=\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)}$<br/>
 +
故に<br/>
 +
$\frac{d r(t) }{dt}=\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)} \qquad \qquad (b)$<br/>
 +
補題 $g(\theta)\triangleq \frac{1}{f(\theta)}$とおくと、<br/>
 +
$\frac{df(\theta)}{d\theta}= -\frac{1}{g^2(\theta)}\frac{dg(\theta)}{d\theta}$
 +
<br/>
 +
証明は、
 +
[[wikibooks_ja:高等学校数学III/微分法#商の導関数|ウィキブックス高等学校数学III/微分法#商の導関数]]をご覧ください。<br/><br/>
 +
この補題を用いて式(b)を変形すると<br/>
 +
$\frac{d r(t) }{dt}=-\frac{1}{g^2(\theta (t))}\frac{dg}{d\theta}(\theta (t))\frac{h}{r^2(t)} $<br/>
 +
$\qquad r(t)=f(\theta (t))$ なので、$g^2(\theta (t))r^2(t)=1$となり、<br/>
 +
$=-h\frac{dg}{d\theta}(\theta (t)) $<br/>
 +
故に<br/>
 +
$\dot{r}(t)=-h\frac{dg}{d\theta}(\theta (t)) \qquad \qquad (c)$<br/>
 +
これをさらに時間tで微分すると、合成関数の微分の公式から<br/>
 +
$\ddot{r}(t)=-h\frac{d^{2}g}{d\theta ^2}(\theta(t))\dot{\theta}(t)$<br/>
 +
$\qquad $ 式(13)から得られる$\dot{\theta}=\frac{h}{r^2}$を代入して、<br/>
 +
$=-h\frac{d^{2}g}{d\theta ^2}(\theta(t)) \frac{h}{r^2(t)} $<br/>
 +
故に<br/>
 +
$\ddot{r}(t)=-\frac{h^{2}}{r^2(t)}\frac{d^{2}g}{d\theta ^2}(\theta(t))\qquad \qquad (d)$<br/>
 +
が得られる。<br/><br/>
 +
式(d)と$\dot{\theta}=\frac{h}{r^2}$を、式(12)に代入して整頓すると次式が得られる。<br/>
 +
$\frac{d^{2}g}{d\theta ^2}(\theta(t))+\frac{1}{r(t)}=\frac{GM}{h^2}  $<br/>
 +
$\frac{1}{r(t)}=\frac{1}{f(\theta (t))}=g(\theta (t))$なので<br/>
 +
上式から次式が得られる。<br/>
 +
$\frac{d^{2}g}{d\theta ^2}(\theta(t))+g(\theta (t))=\frac{GM}{h^2} $<br/>
 +
従って変数を $\theta$ とする微分方程式<br/>
 +
$\frac{d^{2}g}{d\theta ^2}(\theta)+g(\theta )=\frac{GM}{h^2} \qquad \qquad (14)$<br/>
 +
の解$g=g(\theta)$が分かれば<br/>
 +
$r=f(\theta)=\frac{1}{g(\theta)}$  <br/>
 +
が、惑星の軌道の形を表す式になる。
 +
===== 微分方程式の解法  =====
 +
 
 +
==== 極座標表示の惑星軌道を、直交座標表示に変換  ====

2018年6月25日 (月) 06:11 時点における最新版

目次

 ☆☆惑星の運動(2)

2.3節の惑星運動の項では、
惑星の軌道を円と仮定すると、
運動の第1、第2法則のもとでは、
ケプラーの法則と万有引力の法則が同等であることを示した。
この節では軌道について仮定せず、
運動の第1、第2法則のもとで、
万有引力の法則と惑星に関するケプラーの3法則が同等であることを導く。

 ニュートンの運動の第一、第二法則による惑星の運動方程式の導出

他の惑星の影響を無視して、太陽と一つの惑星の二体からなる系をある慣性系から観測する。(注1と注2を参照のこと)
太陽の質量を$M$、惑星の質量を$m$とする。
この二体を、ともに質点であるとみなし、位置はそれぞれの重心とする。
この系には、系外から力が作用しないため、、
系の重心は等速直線運動をする。
従って、この重心を原点とする回転しない座標系も慣性系となる。
今後はこの慣性系で観測する。

(注1)慣性系の存在はニュートン力学(古典力学)の大前提であり、運動法則の第一法則としてまとめられている。
太陽系の重心に原点を取り、恒星系に対して回転しない座標系は慣性系に近いことが現在分かっている。

慣性系からみた質点系の運動は、テキストの2.3 で、その導入部分を説明している。
(注2)他の惑星の影響を無視し、太陽と当該惑星の2体を質点とみなしても、
他の惑星からの影響の少ない地球などの惑星の運動は、かなり正確に予測できることが知られている。
正確にするには二体ではなく他の惑星も考慮した多体問題として論じなければならないが、解析は非常に難しくなる。

 2体の重心は太陽の重心と一致するという近似を行う

 太陽の質量は惑星の質量より圧倒的大きいため二体の重心は太陽の重心とほぼ一致する。


そこで議論を簡単にするため、このテキストでは2体の重心は太陽の重心と一致し、太陽は不動であると仮定する。
興味のある読者は、この仮定をしないで、惑星と太陽の運動を解析してみてほしい。
議論は殆ど同じように進められます。

そこで太陽$S$を原点$O$とする慣性座標系を考える。
任意の時刻$\ t\ $の惑星$P_t$の位置ベクトル$\vec{OP_t}$を$\vec{r}(t)$とかく。
惑星が太陽からうける力(万有引力の法則を仮定しないと未知)を仮に$\vec f$と書くと、
惑星の運動は、ニュートンの第2法則から次の微分方程式で形式的には記述されが、
力が未知なので、この方程式は解を持ち得ない。
$\ M\frac{d^{2} \vec{r}}{dt^2}(t)=\vec{f} $

 運動の第1、第2法則と万有引力の法則からケプラーの3法則を導く 

 惑星は平面上を運動する

不動の太陽の位置$O$,時刻tの惑星の位置$P(t_0)$と$P(t_0)$を始点とした初速ベクトル$\vec{v}_0$を含む
平面$H$を考える。

この平面と直交する一つのベクトル$\vec h$を考え、
式(1)の両辺の左側からベクトル積として作用させる。
$\ \vec h \times \Bigl(M\frac{d \vec{v}}{dt}(t)\Bigr)=\ \vec h \times \Bigl(-\frac{GMm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}\Bigr) \qquad \qquad (3)$
ベクトル積の性質から
  左辺は
$\ \vec h \times \Bigl(M\frac{d \vec{v}}{dt}(t)\Bigr)=\ M\vec h \times \Bigl(\frac{d \vec{v}}{dt}(t)\Bigr)=\ M\vec h \times \Bigl(\frac{d^2 \vec{r}}{dt^2}(t)\Bigr)$
$=\ M\Bigl(\frac{d^2 \vec h \times\vec{r}}{dt^2}(t)\Bigr)$(注1をみよ)
右辺は
$\ \vec h \times \Bigl(-\frac{GMm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}\Bigr) = -\frac{GMm}{\|\vec{r}(t)\|^{3}}\vec h \times \vec{r}(t)$
なので、 $\ M\Bigl(\frac{d^2 \vec h \times\vec{r}}{dt^2}(t)\Bigr) = -\frac{GMm}{\|\vec{r}(t)\|^{3}}\vec h \times \vec{r}(t) $
が得られる。
故に
$\ \Bigl(\frac{d^2 \vec h \times\vec{r}}{dt^2}(t)\Bigr) = -\frac{Gm}{\|\vec{r}(t)\|^{3}}\vec h \times \vec{r}(t)\qquad \qquad (4)$

$\vec h$をベクトル積として初期条件に作用させると
$\vec h \times \vec{r}(t_0) = \vec h \times \vec{r_0}=0, \frac{d \vec h \times \vec{r}}{dt}(t_0) =\vec h \times \frac{d\vec{r}}{dt}(t_0) =\vec h \times \vec{v_0} = 0$
故に、
$\vec h \times \vec{r}(t_0) = 0, \frac{d \vec h \times \vec{r}}{dt}(t_0) = 0\qquad \qquad (5)$
tの関数$\vec{x}\triangleq \vec h \times \vec{r}$を考えると、
式(4),(5)は
  $\ \Bigl(\frac{d^2 \vec{x}}{dt^2}(t)\Bigr) = -\frac{Gm}{\|\vec{r}(t)\|^{3}}\vec{x}(t)\qquad \qquad (4')$
$\vec{x}(t_0) = 0,\quad \frac{d \vec{x}}{dt}(t_0) = 0\qquad \qquad \qquad (5')$
常微分方程式(4')を初期条件(5')のもとで解くと、
$\vec{x}\triangleq \vec h \times \vec{r}\equiv 0$
が得られる。(注2参照)
これよりベクトル$\vec{r}(t)$は,常にベクトル$\vec{h}$と直交し、 平面H上にあることが証明された。

(注1)

(注2)証明をみれば明らかなように、
惑星に働く力が、太陽と惑星を結ぶ直線と同じ方向(求心力あるいが遠心力)でありさえすれば、
惑星は平面上を運動する。

 惑星の角運動量は保存される

式(1)の両辺に、左側から$\vec{r}(t)$をベクトル積として作用させると、
$\vec{r}(t)\times \ \frac{d^{2} \vec{r}}{dt^2}(t)=-\vec{r}(t)\times \frac{Gm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} \qquad \qquad (a)$
ここで、ベクトル積の性質から
$\quad \frac{d}{dt}\Bigl(\vec{r}\times \frac{d \vec r}{dt}\Bigr)=\frac{d \vec r}{dt}\times \frac{d \vec{r}}{dt}+\vec{r}\times \ \frac{d^{2} \vec{r}}{dt^2} $
$\quad =\vec{r}\times \ \frac{d^{2} \vec{r}}{dt^2}$
$\quad$また右辺は
$\quad -\vec{r}(t)\times \frac{Gm}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} = 0$(注参照)
なので、式(a)は、
$\frac{d}{dt}\Bigl(\vec{r}\times \frac{d \vec r}{dt}\Bigr)=0$
となる。
(注)万有引力でなく、惑星に作用する力の向きが原点向きであるか、それと逆向きという条件だけで良い。

$\frac{d \vec r}{dt}=\vec v$(速度ベクトル)なので、上式は
$\frac{d}{dt}\Bigl(\vec{r}\times \vec v \Bigr)=0\qquad \qquad (b)$
あるいは、惑星の運動量$\vec p=m\vec v$を用いて
$\frac{d}{dt}\Bigl(\vec{r}\times \vec p \Bigr)=0\qquad \qquad (c)$
となる。この2式から、
$\vec{r}\times \vec v $ と$\vec{r}\times \vec p $
は、ともに時不変(時間が経過しても一定値)であることが分かる。

定義(角運動量)
$\vec{r}(t)\times \vec{p}(t) $を、(原点からみた)惑星の時刻tにおける角運動量という。

すると上述したことから、次の法則が証明された。

角運動量の保存法則
求心力、ないし遠心力(力の向きが原点向きか、その逆向き)だけを受けて運動する質点の角運動量は保存される。

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

上述したことから惑星は平面上を動くので、
位置ベクトルと運動量ベクトルが、xy平面に載る様な正規直交右手系xyzを選べる。
この座標系のもとで座標表示すると
$\vec{r(t)}\times \vec p(t)=\bigl(r_{1}(t),r_{2}(t),0)\bigr)\times \bigl(p_{1}(t),p_{2}(t),0\bigr)$
$=\bigl(0,0,r_{1}(t)p_{2}(t)-r_{2}(t)p_{1}(t)\bigr)$

他方、ベクトル$\vec{r}(t)$ と ベクトル$\vec{p}(t)$のなす劣角(180度以下のほうの角)を$\theta(t)$とすると
この2つのベクトルと直交するz方向の$\vec{r}(t)\times \vec{p}(t) $の成分は
ベクトル積の定義から、
$\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t)$
なので、
$\|\vec{r}(t)\| \times \| \vec{p}(t) \| \sin \theta(t) \qquad \qquad \qquad (6)$
が、時不変であることが分かる。
惑星質量は一定なので 式(6)から
$\|\vec{r}(t)\| \times \| \vec{v}(t) \| \sin \theta(t)/2 \qquad \qquad \qquad (7)$
が、時不変であることが分かる。
式(7)は動径$OP(t)$が時刻tの瞬間に単位時間あたりに掃く面積を表すので、惑星の面積速度と呼ばれる。


系 ケプラーの第2法則
求心力、ないし遠心力(力の向きが原点向きか、その逆向き)だけを受けて運動する質点の面積速度は一定である。従って、
「惑星の動径(惑星と太陽とを結ぶ線分)が単位時間に描く面積(面積速度)は、一定である」

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

惑星の運動方程式(1)を解いて惑星の軌道を求めよう。
求心力の下での運動を扱うには、極座標が都合がよい。

 惑星運動の微分方程式の極座標表示  

そこで惑星の位置ベクトルを、x軸の正の部分のつくる半直線を角度の基準とする極座標で
$\vec{r}(t)=r(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$
と表す。
この式を時間tで微分すると
$\dot{\vec{r}}(t)={\dot r}(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr) +r(t){\dot \theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$
$\qquad $但し、$\dot{\vec{r}}(t)\triangleq \frac{d\vec{r}}{dt}(t)$
さらにtで微分すると
$\ddot{\vec{r}}(t)$
$=\ddot {r}(t)\bigl(\cos \theta (t),\sin \theta (t)\bigr) +\dot{r}(t)\dot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$
$+\dot{r}(t)\dot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$
$+r(t)\ddot{\theta} (t)\bigl(-\sin \theta (t),\cos \theta (t)\bigr) -r(t)\dot{\theta}^{2}(t)\bigl(\cos \theta (t),\sin \theta (t) \bigr) $
$=\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$
故に
$\ddot{\vec{r}}(t)$
$=\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr) \qquad \qquad \qquad (8)$
が得られた。
惑星に働く万有引力を、この極座標で表示すると
$-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|} =-\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr) \qquad (9)$
式(8),(9)を惑星の運動方程式
$\ddot{\vec r}(t)=-\frac{GM}{\|\vec{r}(t)\|^{2}}\frac{\vec{r}(t)}{\|\vec{r}(t)\|}$
に代入すると
$\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$
$= -\frac{GM}{r^{2}(t)}(\cos \theta (t),\sin \theta (t)\bigr)$
右辺を左辺に移行すると
$\bigl(\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\bigr)\bigl(\cos \theta (t),\sin \theta (t)\bigr)$
$+\bigl(2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \bigr)\bigl(-\sin \theta (t),\cos \theta (t)\bigr)\equiv 0$
$\bigl(\cos \theta (t),\sin \theta (t)\bigr)$と$\bigl(-\sin \theta (t),\cos \theta (t)\bigr)$は直交ベクトルなので、この式から

$\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)+\frac{GM}{r^{2}(t)}\equiv 0 \qquad \qquad (10) $
$2\dot{r}(t)\dot{\theta} (t) +r(t)\ddot{\theta} (t) \equiv 0 \qquad \qquad (11) $
が得られる。

命題;惑星運動の微分方程式の極座標表示
太陽の位置Oを原点、x軸の正の部分の半直線を角の基準線とする極座標で 惑星の運動方程式を表現すると
$\ddot {r}-r\dot{\theta}^{2}=-\frac{GM}{r^{2}}\qquad \qquad (12) $
$r^2\dot{\theta} = h \quad (任意の定数) \qquad \qquad (13) $
となる。(注2 を参照のこと)
証明
式(12)は式(10)から明らか。
式(13)の左辺をtで微分すると $\frac{d r^2\dot{\theta}}{dt}=r^2\ddot{\theta} + 2r\dot{r}\dot{\theta} =r\bigl(r\ddot{\theta} + 2\dot{r}\dot{\theta}\bigr)$
$\qquad $式(11)から
$= 0$
左辺の時間微分が零なので、左辺は任意の定数hに等しい。
証明終わり。 $\qquad \qquad \qquad \qquad \Box$

(注)式(13)は面積速度が一定であることを表現している。
先に、このことは証明してあるが、その別証明である。
(注2)関数の変数tまで書くと、式(12),(13)は、次のように書ける。
$\ddot {r}(t)-r(t)\dot{\theta}^{2}(t)=-\frac{GM}{r^{2}(t)}\qquad \qquad (12') $
$r^2(t)\dot{\theta}(t) = h \quad (任意の定数) \qquad \qquad (13') $

 微分方程式を解いて、惑星の軌道形を求める  

連立微分方程式(12)、(13)を解いて、$r=r(t),\theta=\theta(t)$を求めれば、
任意の時刻tにおける惑星の位置がわかる。
軌道の形は、この2式からtを消去して、rとθの関数関係
$f(r,\theta)=0$, または$r=f(\theta)$
を求めれな良い。
しかし、残念なことにこれは大変困難である。

軌道の形だけを知りたいときには、
これを回避して簡単な(しかし技巧的な)直接的な解法がある。
それは、連立微分方程式の解が
$r=f\circ \theta ,\qquad$変数まで書けば、$r(t)=f\circ \theta (t)\triangleq f(\theta (t))$
という合成関数になることを仮定して、
$r=f(\theta)$を直接求める問題に変換し、解を求める方法である。

 軌道の形を直接求めるための微分方程式の導出

$\frac{d r(t) }{dt}=\frac{d f\circ \theta (t) }{dt}$
$\qquad $ 合成関数の微分公式を用いて、
$=\frac{d f}{d\theta}(\theta (t))\dot{\theta}(t)$

$\qquad $式(13)を$\dot{\theta}=\frac{h}{r^2}$に変形し、上式に代入すると
$=\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)}$
故に
$\frac{d r(t) }{dt}=\frac{df}{d\theta}(\theta (t))\frac{h}{r^2(t)} \qquad \qquad (b)$
補題 $g(\theta)\triangleq \frac{1}{f(\theta)}$とおくと、
$\frac{df(\theta)}{d\theta}= -\frac{1}{g^2(\theta)}\frac{dg(\theta)}{d\theta}$
証明は、 ウィキブックス高等学校数学III/微分法#商の導関数をご覧ください。

この補題を用いて式(b)を変形すると
$\frac{d r(t) }{dt}=-\frac{1}{g^2(\theta (t))}\frac{dg}{d\theta}(\theta (t))\frac{h}{r^2(t)} $
$\qquad r(t)=f(\theta (t))$ なので、$g^2(\theta (t))r^2(t)=1$となり、
$=-h\frac{dg}{d\theta}(\theta (t)) $
故に
$\dot{r}(t)=-h\frac{dg}{d\theta}(\theta (t)) \qquad \qquad (c)$
これをさらに時間tで微分すると、合成関数の微分の公式から
$\ddot{r}(t)=-h\frac{d^{2}g}{d\theta ^2}(\theta(t))\dot{\theta}(t)$
$\qquad $ 式(13)から得られる$\dot{\theta}=\frac{h}{r^2}$を代入して、
$=-h\frac{d^{2}g}{d\theta ^2}(\theta(t)) \frac{h}{r^2(t)} $
故に
$\ddot{r}(t)=-\frac{h^{2}}{r^2(t)}\frac{d^{2}g}{d\theta ^2}(\theta(t))\qquad \qquad (d)$
が得られる。

式(d)と$\dot{\theta}=\frac{h}{r^2}$を、式(12)に代入して整頓すると次式が得られる。
$\frac{d^{2}g}{d\theta ^2}(\theta(t))+\frac{1}{r(t)}=\frac{GM}{h^2} $
$\frac{1}{r(t)}=\frac{1}{f(\theta (t))}=g(\theta (t))$なので
上式から次式が得られる。
$\frac{d^{2}g}{d\theta ^2}(\theta(t))+g(\theta (t))=\frac{GM}{h^2} $
従って変数を $\theta$ とする微分方程式
$\frac{d^{2}g}{d\theta ^2}(\theta)+g(\theta )=\frac{GM}{h^2} \qquad \qquad (14)$
の解$g=g(\theta)$が分かれば
$r=f(\theta)=\frac{1}{g(\theta)}$
が、惑星の軌道の形を表す式になる。

 微分方程式の解法

 極座標表示の惑星軌道を、直交座標表示に変換  

個人用ツール