院試hub

中心力とケプラーの法則有効ポテンシャル・軌道の方程式・惑星の周期

力が常に一点を向き大きさが距離 rr だけで決まる中心力では、角運動量 L\boldsymbol{L} と力学的エネルギーが保存し、運動は一つの平面内に限られる。動径方向の運動は、遠心力の項を加えた有効ポテンシャル Ueff(r)=L2/(2mr2)+U(r)U_{\text{eff}}(r) = L^2/(2mr^2) + U(r) の中の1次元運動として扱える。逆2乗の引力 U=−k/rU = -k/r では軌道が円錐曲線 r=p/(1+ecos⁡θ)r = p/(1 + e\cos\theta) になり、楕円軌道の周期は T2=4π2a3/(GM)T^2 = 4\pi^2a^3/(GM) を満たす。

力学 の第3単元(全4単元)。先に読む単元:運動の法則と保存則

力学を出題する院試の解答・解説PDF

大学・研究科ごとに、過去問の解答と解説を年度別にPDFにしている。収録年度・見本・価格を各商品ページで確認できる。

科目から解答PDFを探す

中心力と保存量

定義

中心力

原点に置いた力の中心からの位置を r\boldsymbol{r} として、F=f(r) r/r\boldsymbol{F} = f(r)\,\boldsymbol{r}/r の形の力を中心力という。f(r)<0f(r) < 0 が引力である。このとき F=−∇U\boldsymbol{F} = -\nabla U、U(r)=−∫f(r) drU(r) = -\int f(r)\,dr と書けるので、中心力は保存力である。

定理・法則

角運動量保存と面積速度一定(ケプラーの第2法則)

条件

質点に働く力が、固定した中心 O に関する中心力だけである。

主張

L=mr×r˙\boldsymbol{L} = m\boldsymbol{r} \times \dot{\boldsymbol{r}} は一定で、運動は L\boldsymbol{L} に垂直な平面内に限られる。その平面の極座標 (r,θ)(r, \theta) で

L=mr2θ˙=一定,dAdt=12r2θ˙=L2mL = mr^2\dot{\theta} = \text{一定}, \qquad \frac{dA}{dt} = \frac{1}{2}r^2\dot{\theta} = \frac{L}{2m}

要点

dL/dt=r×F=0d\boldsymbol{L}/dt = \boldsymbol{r} \times \boldsymbol{F} = \boldsymbol{0} から出る。面積速度一定は逆2乗則に限らず、すべての中心力で成り立つ。

極座標での速度の2乗は r˙2+r2θ˙2\dot{r}^2 + r^2\dot{\theta}^2 なので、θ˙=L/(mr2)\dot{\theta} = L/(mr^2) を使ってエネルギーは

E=12mr˙2+L22mr2+U(r)=12mr˙2+Ueff(r)E = \frac{1}{2}m\dot{r}^2 + \frac{L^2}{2mr^2} + U(r) = \frac{1}{2}m\dot{r}^2 + U_{\text{eff}}(r)

となる。E≥Ueff(r)E \ge U_{\text{eff}}(r) を満たす rr の範囲だけを動ける。UeffU_{\text{eff}} の極小点が円軌道で、EE が極小値より大きく零より小さければ rr は2つの転回点の間を往復する(束縛軌道)。E≥0E \ge 0 で U(∞)=0U(\infty) = 0 なら無限遠へ飛び去る。

軌道の方程式と円錐曲線

定理・法則

軌道の方程式(ビネの公式)

条件

中心力 f(r)f(r) のもとで L≠0L \neq 0。u=1/ru = 1/r を θ\theta の関数とみる。

主張

d2udθ2+u=−mL2u2 f ⁣(1u)\frac{d^2u}{d\theta^2} + u = -\frac{m}{L^2u^2}\,f\!\left(\frac{1}{u}\right)

要点

r˙=−Lmdudθ\dot{r} = -\dfrac{L}{m}\dfrac{du}{d\theta}、r¨=−L2u2m2d2udθ2\ddot{r} = -\dfrac{L^2u^2}{m^2}\dfrac{d^2u}{d\theta^2} を動径方向の運動方程式 m(r¨−rθ˙2)=f(r)m(\ddot{r} - r\dot{\theta}^2) = f(r) に代入して tt を消す。

定理・法則

逆2乗の引力の軌道(ケプラーの第1法則)

条件

f(r)=−k/r2f(r) = -k/r^2(k>0k > 0)、太陽と惑星なら k=GMmk = GMm で M≫mM \gg m とする。

主張

r=p1+ecos⁡(θ−θ0),p=L2mk,e=1+2EL2mk2r = \frac{p}{1 + e\cos(\theta - \theta_0)}, \qquad p = \frac{L^2}{mk}, \qquad e = \sqrt{1 + \frac{2EL^2}{mk^2}}

0≤e<10 \le e < 1(E<0E < 0)で楕円、e=1e = 1(E=0E = 0)で放物線、e>1e > 1(E>0E > 0)で双曲線。楕円の長半径は a=p/(1−e2)=k/(2∣E∣)a = p/(1 - e^2) = k/(2|E|) である。

要点

右辺は定数 mk/L2mk/L^2 になり、u=(mk/L2)(1+ecos⁡θ)u = (mk/L^2)(1 + e\cos\theta) が解である。2体問題では mm を換算質量 μ=mM/(m+M)\mu = mM/(m + M) に、相対座標を r\boldsymbol{r} に置き換えれば同じ式になる。

定理・法則

ケプラーの第3法則

条件

逆2乗の引力による楕円軌道で M≫mM \gg m。

主張

T2=4π2GMa3T^2 = \frac{4\pi^2}{GM}a^3

要点

楕円の面積 πab\pi ab(b=a1−e2b = a\sqrt{1 - e^2})を面積速度 L/(2m)L/(2m) で割ると周期になり、L2=GMm2a(1−e2)L^2 = GMm^2a(1 - e^2) を代入すると ee が消える。M≫mM \gg m でないときは GMGM を G(M+m)G(M + m) に置き換える。

軌道計算の公式

円軌道の速さ
v=GMrv = \sqrt{\frac{GM}{r}}

mv2/r=GMm/r2mv^2/r = GMm/r^2 から。UeffU_{\text{eff}} の極小点 r0=L2/(mk)r_0 = L^2/(mk) でもある。

脱出速度
vesc=2GMRv_{\text{esc}} = \sqrt{\frac{2GM}{R}}

12mv2−GMm/R=0\frac{1}{2}mv^2 - GMm/R = 0 から。空気抵抗と天体の自転は無視する。

活力の式(vis-viva)
v2=GM(2r−1a)v^2 = GM\left(\frac{2}{r} - \frac{1}{a}\right)

楕円軌道上の任意の点で成り立つ。E=−GMm/(2a)E = -GMm/(2a) とエネルギー保存から出る。

近点・遠点の距離
rp=a(1−e),ra=a(1+e)r_{\text{p}} = a(1 - e), \qquad r_{\text{a}} = a(1 + e)
近点・遠点の速さ
vp=GMa⋅1+e1−e,rpvp=ravav_{\text{p}} = \sqrt{\frac{GM}{a}\cdot\frac{1 + e}{1 - e}}, \qquad r_{\text{p}}v_{\text{p}} = r_{\text{a}}v_{\text{a}}

近点と遠点では速度が動径に垂直なので、L=mrvL = mrv がそのまま使える。

中心力の計算例

計算例

地球からの脱出速度

地球の GM=3.986×1014 m3/s2GM = 3.986 \times 10^{14}\ \mathrm{m^3/s^2}、平均半径 R=6.371×106 mR = 6.371 \times 10^6\ \mathrm{m} として、地表からの脱出速度を求めよ。

解答

vesc=2×3.986×10146.371×106≈1.12×104 m/sv_{\text{esc}} = \sqrt{\frac{2 \times 3.986 \times 10^{14}}{6.371 \times 10^6}} \approx 1.12 \times 10^4\ \mathrm{m/s}

およそ 11.2 km/s11.2\ \mathrm{km/s} である。

検算

地表すれすれの円軌道の速さ GM/R≈7.9 km/s\sqrt{GM/R} \approx 7.9\ \mathrm{km/s} の 2\sqrt{2} 倍になっている。(m3/s2)/m=m/s\sqrt{(\mathrm{m^3/s^2})/\mathrm{m}} = \mathrm{m/s}。

計算例

静止軌道の半径

地球の自転周期(恒星日)T=86164 sT = 86164\ \mathrm{s} と同じ周期の円軌道の半径を求めよ。

解答

第3法則から

a=(GMT24π2)1/3=(3.986×1014×8616424π2)1/3≈4.216×107 ma = \left(\frac{GMT^2}{4\pi^2}\right)^{1/3} = \left(\frac{3.986 \times 10^{14} \times 86164^2}{4\pi^2}\right)^{1/3} \approx 4.216 \times 10^7\ \mathrm{m}

地心から約 42,200 km42{,}200\ \mathrm{km}、地表からは約 35,800 km35{,}800\ \mathrm{km} の高さである。

検算

このときの速さ GM/a≈3.07 km/s\sqrt{GM/a} \approx 3.07\ \mathrm{km/s} と 2πa/T≈3.07 km/s2\pi a/T \approx 3.07\ \mathrm{km/s} が一致する。周期に太陽日 86400 s86400\ \mathrm{s} を使うと半径がわずかにずれるので、恒星日を使う。

計算例

逆2乗則で軌道の方程式を解く

f(r)=−k/r2f(r) = -k/r^2 のとき、軌道の方程式を解いて r(θ)r(\theta) を求めよ。

解答

f(1/u)=−ku2f(1/u) = -ku^2 なので

d2udθ2+u=mkL2\frac{d^2u}{d\theta^2} + u = \frac{mk}{L^2}

定数の特解 mk/L2mk/L^2 と斉次解 Ccos⁡(θ−θ0)C\cos(\theta - \theta_0) を足して u=mkL2{1+ecos⁡(θ−θ0)}u = \dfrac{mk}{L^2}\{1 + e\cos(\theta - \theta_0)\}(e=CL2/(mk)e = CL^2/(mk))。逆数をとって r=p/{1+ecos⁡(θ−θ0)}r = p/\{1 + e\cos(\theta - \theta_0)\}、p=L2/(mk)p = L^2/(mk)。

検算

e=0e = 0 では r=L2/(mk)r = L^2/(mk) の円で、UeffU_{\text{eff}} を微分して零とおいた半径と一致する。そのとき E=Ueff(r0)=−mk2/(2L2)E = U_{\text{eff}}(r_0) = -mk^2/(2L^2) で、ee の式に代入すると e=0e = 0 になる。

計算例

近点と遠点の速さの比

離心率 ee の楕円軌道で、近点の速さ vpv_{\text{p}} と遠点の速さ vav_{\text{a}} の比を求めよ。

解答

近点と遠点では速度が動径に垂直なので、角運動量保存は mrpvp=mravamr_{\text{p}}v_{\text{p}} = mr_{\text{a}}v_{\text{a}} となる。

vpva=rarp=1+e1−e\frac{v_{\text{p}}}{v_{\text{a}}} = \frac{r_{\text{a}}}{r_{\text{p}}} = \frac{1 + e}{1 - e}

検算

活力の式から vp2=GMa1+e1−ev_{\text{p}}^2 = \frac{GM}{a}\frac{1 + e}{1 - e}、va2=GMa1−e1+ev_{\text{a}}^2 = \frac{GM}{a}\frac{1 - e}{1 + e} で、比をとると同じ結果になる。e=0e = 0 では比が 11 で円軌道に一致する。

答案でよくある誤り

  1. 1

    軌道上の任意の点で L=mrvL = mrv とする。L=mrvsin⁡αL = mrv\sin\alpha(α\alpha は r\boldsymbol{r} と v\boldsymbol{v} のなす角)で、mrvmrv になるのは速度が動径に垂直な近点・遠点と円軌道だけである。

  2. 2

    有効ポテンシャルに遠心力の項 L2/(2mr2)L^2/(2mr^2) を入れ忘れる。θ˙\dot{\theta} を LL で消すと動径方向の運動にこの項が現れる。

  3. 3

    第3法則の aa に近点距離や平均距離を入れる。aa は楕円の長半径 (rp+ra)/2(r_{\text{p}} + r_{\text{a}})/2 である。

  4. 4

    2体問題で中心の天体を固定したまま扱う。質量が同程度なら換算質量 μ\mu と相対座標を使い、周期の式では G(M+m)G(M + m) にする。

  5. 5

    面積速度一定を逆2乗則の性質と考える。角運動量保存から出るので、すべての中心力で成り立つ。楕円軌道と第3法則の形は逆2乗則に固有である。

力学の前後の単元

このページの本文を数式のLaTeXのまま読むには Markdown版。最終更新 2026-11-16。誤りに気づいたら 運営者 へ連絡してほしい。

力学を出題する院試の解答・解説PDF

大学・研究科ごとに、過去問の解答と解説を年度別にPDFにしている。収録年度・見本・価格を各商品ページで確認できる。

科目から解答PDFを探す