物理/惑星の運動(2)
提供: Internet Web School
(ページの作成: === ☆☆惑星の運動(2) === この節(2.3)ですでに説明した惑星運動の項では、<br/> 惑星の軌道…) |
(→ ☆☆惑星の運動(2)) |
||
1 行: | 1 行: | ||
- | + | = ☆☆惑星の運動(2) = | |
- | + | 2.3節の[[物理/質点の運動と質点系#惑星運動|惑星運動]]の項では、<br/> | |
- | + | 惑星の軌道を円と仮定すると、<br/> | |
+ | 運動の第2法則のもとでは、<br/> | ||
ケプラーの法則と万有引力の法則が同等であることを示した。<br/> | ケプラーの法則と万有引力の法則が同等であることを示した。<br/> | ||
この節では軌道について仮定せず、<br/> | この節では軌道について仮定せず、<br/> | ||
運動の第2法則のもとで、<br/> | 運動の第2法則のもとで、<br/> | ||
- | 万有引力の法則と惑星に関するケプラーの3法則が同等であることを導く。 | + | 万有引力の法則と惑星に関するケプラーの3法則が同等であることを導く。 |
- | + | == 惑星の運動方程式 == | |
他の惑星の影響を無視して、太陽と一つの惑星の2体からなる系を考える。<br/> | 他の惑星の影響を無視して、太陽と一つの惑星の2体からなる系を考える。<br/> | ||
太陽の質量をM、惑星の質量をmとする。<br/> | 太陽の質量をM、惑星の質量をmとする。<br/> | ||
30 行: | 31 行: | ||
ここで、式(1)の右辺中のGは万有引力定数である。<br/> | ここで、式(1)の右辺中のGは万有引力定数である。<br/> | ||
- | + | == 惑星は平面上を運動する == | |
O,P(t0)とP(t0)を始点とした初速→v0を含む<br/> | O,P(t0)とP(t0)を始点とした初速→v0を含む<br/> | ||
平面Hを考える。<br/> | 平面Hを考える。<br/> | ||
66 行: | 67 行: | ||
(注2)<br/> | (注2)<br/> | ||
- | + | == 運動の第2法則と万有引力の法則からケプラーの3法則を導く == | |
- | === | + | === ケプラーの第2法則の導出 === |
第2法則(面積速度一定の法則)、<br/> | 第2法則(面積速度一定の法則)、<br/> | ||
「惑星と太陽とを結ぶ線分が単位時間に描く面積(面積速度)は、一定である」<br/> | 「惑星と太陽とを結ぶ線分が単位時間に描く面積(面積速度)は、一定である」<br/> | ||
73 行: | 74 行: | ||
====== 惑星運動の微分方程式の極座標表示 ====== | ====== 惑星運動の微分方程式の極座標表示 ====== | ||
+ | |||
+ | === ケプラーの第1法則の導出 === | ||
+ | |||
+ | === ケプラーの第3法則の導出 === | ||
+ | == 運動の第2法則とケプラーの3法則から万有引力の法則を導く == |
2018年6月3日 (日) 04:11時点における版
目次[非表示] |
☆☆惑星の運動(2)
2.3節の惑星運動の項では、
惑星の軌道を円と仮定すると、
運動の第2法則のもとでは、
ケプラーの法則と万有引力の法則が同等であることを示した。
この節では軌道について仮定せず、
運動の第2法則のもとで、
万有引力の法則と惑星に関するケプラーの3法則が同等であることを導く。
惑星の運動方程式
他の惑星の影響を無視して、太陽と一つの惑星の2体からなる系を考える。
太陽の質量をM、惑星の質量をmとする。
この2体を、ともに質点であるとみなし、位置はそれぞれの重心とする。
この系には、系外から力が作用しないため、任意の慣性系から観測すると、
系の重心は等速直線運動をする。
ところが、太陽の質量は惑星の質量より圧倒的大きいため2体の重心は太陽の重心とほぼ一致し、
太陽重心は慣性系からみると等速直線運動しているとみなしてよい。
(注)他の惑星の影響を無視し、太陽と惑星を質点とみなしても誤差は小さいことが知られている。
そこで太陽Sを原点Oとする静止座標系を考える。
任意の時刻 t の惑星Ptの位置ベクトル→OPtを→r(t)、
速度ベクトルを→v(t)とかく。
→v(t)=d→rdt(t))である。
指定した初期時刻t0における惑星の
位置ベクトル→r(t0)を、→r(t0)=→r0,
速度ベクトル→v(t0)を、→v(t0)=→v0,
と仮定する。
すると惑星の運動は、ニュートンの第2法則と万有引力の法則より次の微分方程式で記述される。
Md→vdt(t)=−GMm‖→r(t)‖2→r(t)‖→r(t)‖(1)
→r(t0)=→r0、→v(t0)=→v0(2)
ここで、式(1)の右辺中のGは万有引力定数である。
惑星は平面上を運動する
O,P(t0)とP(t0)を始点とした初速→v0を含む
平面Hを考える。
この平面と直交する一つのベクトル→hを考え、
式(1)の両辺の左側から外積として作用させる。
→h×(Md→vdt(t))= →h×(−GMm‖→r(t)‖2→r(t)‖→r(t)‖)(3)
外積の性質から
左辺は
→h×(Md→vdt(t))= M→h×(d→vdt(t))= M→h×(d2→rdt2(t))
= M(d2→h×→rdt2(t))(注1をみよ)
右辺は
→h×(−GMm‖→r(t)‖2→r(t)‖→r(t)‖)=−GMm‖→r(t)‖3→h×→r(t)
なので、
M(d2→h×→rdt2(t))=−GMm‖→r(t)‖3→h×→r(t)
が得られる。
故に
(d2→h×→rdt2(t))=−Gm‖→r(t)‖3→h×→r(t)(4)
→hを外積として初期条件に作用させると
→h×→r(t0)=→h×→r0=0,d→h×→rdt(t0)=→h×d→rdt(t0)=→h×→v0=0
故に、
→h×→r(t0)=0,d→h×→rdt(t0)=0 (5)
tの関数→x≜を考えると、
式(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)
運動の第2法則と万有引力の法則からケプラーの3法則を導く
ケプラーの第2法則の導出
第2法則(面積速度一定の法則)、
「惑星と太陽とを結ぶ線分が単位時間に描く面積(面積速度)は、一定である」
を導出しよう。