跳到主要内容

東京大学 総合文化研究科 広域科学専攻 広域システム科学系 2019年8月実施 総合科目 第3問

Author

Miyake

Description

質量 m1m_1 の質点 11 と質量 m2m_2 の質点 22 が、お互いに万有引力で引き合いながら真空中を運動している。 質点 22 から質点 11 に向かうベクトルを r\vec{r} とし、 二つの質点の重心を原点 O\text{O} とする座標系を用いる。 二つの質点は同一平面内を運動するが、r\vec{r} をその面内の二次元極座標で表し、 動径および方位角成分を (r,ϕ)(r, \phi) とする。万有引力定数を GG とし、

μ=m1m2m1+m2M=m1+m2\begin{aligned} &\mu = \frac{m_1 m_2}{m_1 + m_2} \\ &M = m_1 + m_2 \end{aligned}

として、以下の設問に答えよ。なお必要に応じて、時間微分を r˙\dot{r}ϕ˙\dot{\phi} などと上付きドットを用いて表しても良い。

(1) 質点 11 および質点 22 の位置ベクトルを、m1m_1, m2m_2 および r\vec{r} 用いて表せ。

(2) 二つの質点の運動エネルギーの総和を μ\mu, rr および ϕ\phi を用いて表せ。

(3) この系の重力ポテンシャルエネルギーを μ\mu, MM, rr および GG を用いて表せ。

(4) lr2ϕl \equiv r^2 \phi が保存量となることを示せ。

(5) 動径成分 rr についての運動方程式を、μ\mu, MM, rr, ll および GG を用いて書き下せ。

(6) 系の全エネルギー EEμ\mu, MM, rr, ll および GG で表せ。

(7) EEが最小値を取る時の EE の表式を μ\mu, MM, ll および GG を用いて求め、その時の二つの質点の軌道の様子を述べよ。

(8) r=(r1,ϕ1)\vec{r}=(r_1, \phi_1) を質点 11 の位置ベクトルとする。 二つの質点のが重力的に束縛されお互いの周りを周回する条件と、その時の r1r_1ϕ1\phi_1 の関係を求めよ。

Kai

(1)

題意より、

m1r1+m2r2=0,    r=r1r2\begin{aligned} m_1 \vec{r}_1 + m_2 \vec{r}_2 = \vec{0} , \ \ \ \ \vec{r} = \vec{r}_1 - \vec{r}_2 \end{aligned}

であるから、

r1=m2m1+m2r,    r2=m1m1+m2r\begin{aligned} \vec{r}_1 = \frac{m_2}{m_1+m_2} \vec{r} , \ \ \ \ \vec{r}_2 = - \frac{m_1}{m_1+m_2} \vec{r} \end{aligned}

を得る。

(2)

r\vec{r} のデカルト座標による成分を (x,y)(x,y) とすると、

x=rcosϕy=rsinϕx˙=r˙cosϕrϕ˙sinϕy˙=r˙sinϕ+rϕ˙cosϕr˙2=x˙2+y˙2=r˙2+r2ϕ˙2\begin{aligned} x &= r \cos \phi \\ y &= r \sin \phi \\ \dot{x} &= \dot{r} \cos \phi - r \dot{\phi} \sin \phi \\ \dot{y} &= \dot{r} \sin \phi + r \dot{\phi} \cos \phi \\ \left| \dot{\vec{r}} \right|^2 &= \dot{x}^2 + \dot{y}^2 = \dot{r}^2 + r^2 \dot{\phi}^2 \end{aligned}

である。

よって、求める運動エネルギー KK は、

K=12m1r1˙2+12m2r2˙2=12m1m22+m12m2(m1+m2)2r˙2=12m1m2m1+m2(r˙2+r2ϕ˙2)=12μ(r˙2+r2ϕ˙2)\begin{aligned} K &= \frac{1}{2} m_1 \left| \dot{\vec{r}_1} \right|^2 + \frac{1}{2} m_2 \left| \dot{\vec{r}_2} \right|^2 \\ &= \frac{1}{2} \frac{m_1 m_2^2 + m_1^2 m_2}{(m_1+m_2)^2} \left| \dot{\vec{r}} \right|^2 \\ &= \frac{1}{2} \frac{m_1 m_2}{m_1+m_2} \left( \dot{r}^2 + r^2 \dot{\phi}^2 \right) \\ &= \frac{1}{2} \mu \left( \dot{r}^2 + r^2 \dot{\phi}^2 \right) \end{aligned}

となる。

(3)

求める重力ポテンシャルエネルギー UU は、

U=Gm1m2r=GμMr\begin{aligned} U = - G \frac{m_1 m_2}{r} = - \frac{G \mu M}{r} \end{aligned}

である。

(4)

この系のラグランジアン LL は、

L=KU=12μ(r˙2+r2ϕ˙2)+GμMr\begin{aligned} L = K - U = \frac{1}{2} \mu \left( \dot{r}^2 + r^2 \dot{\phi}^2 \right) + \frac{G \mu M}{r} \end{aligned}

であり、

Lϕ=0Lϕ˙=μr2ϕ˙=μl\begin{aligned} \frac{\partial L}{\partial \phi} &= 0 \\ \frac{\partial L}{\partial \dot{\phi}} &= \mu r^2 \dot{\phi} = \mu l \end{aligned}

であるから、 ϕ\phi に関する運動方程式より、 ll が保存量であることがわかる。

(5)

上で求めたラグランジアン LL から、

Lr=μrϕ˙2GμMr2=μl2r3GμMr2Lr˙=μr˙\begin{aligned} \frac{\partial L}{\partial r} &= \mu r \dot{\phi}^2 - \frac{G \mu M}{r^2} = \frac{\mu l^2}{r^3} - \frac{G \mu M}{r^2} \\ \frac{\partial L}{\partial \dot{r}} &= \mu \dot{r} \end{aligned}

であるから、 rr に関する運動方程式は、

μr¨=μl2r3GμMr2    r¨=l2r3GMr2\begin{aligned} \mu \ddot{r} &= \frac{\mu l^2}{r^3} - \frac{G \mu M}{r^2} \\ \therefore \ \ \ \ \ddot{r} &= \frac{l^2}{r^3} - \frac{G M}{r^2} \end{aligned}

となる。

(6)

E=K+U=12μ(r˙2+r2ϕ˙2)GμMr=12μ(r˙2+l2r2)GμMr\begin{aligned} E = K + U &= \frac{1}{2} \mu \left( \dot{r}^2 + r^2 \dot{\phi}^2 \right) - \frac{G \mu M}{r} \\ &= \frac{1}{2} \mu \left( \dot{r}^2 + \frac{l^2}{r^2} \right) - \frac{G \mu M}{r} \end{aligned}

(7)

V(r)=12μl2r2GμMr\begin{aligned} V(r) = \frac{1}{2} \mu \frac{l^2}{r^2} - \frac{G \mu M}{r} \end{aligned}

とすると、

V(r)=μl2r3+GμMr2=μl2+GMrr3\begin{aligned} V'(r) &= - \mu \frac{l^2}{r^3} + \frac{G \mu M}{r^2} \\ &= \mu \frac{- l^2 + G M r}{r^3} \end{aligned}

となるので、

r0=l2GM\begin{aligned} r_0 = \frac{l^2}{GM} \end{aligned}

とすると、 V(r)V(r) が最小値をとるのは r=r0r=r_0 のときで、

V(r0)=12G2μM2l2\begin{aligned} V(r_0) = - \frac{1}{2} \frac{G^2 \mu M^2}{l^2} \end{aligned}

となる。

EE が最小値 E0E_0 をとるのもこのときで、

E0=12G2μM2l2\begin{aligned} E_0 = - \frac{1}{2} \frac{G^2 \mu M^2}{l^2} \end{aligned}

である。

(8)