東京大学 総合文化研究科 広域科学専攻 広域システム科学系 2019年8月実施 総合科目 第3問
Author
Miyake
Description
質量 m1 の質点 1 と質量 m2 の質点 2 が、お互いに万有引力で引き合いながら真空中を運動している。
質点 2 から質点 1 に向かうベクトルを r とし、
二つの質点の重心を原点 O とする座標系を用いる。
二つの質点は同一平面内を運動するが、r をその面内の二次元極座標で表し、
動径および方位角成分を (r,ϕ) とする。万有引力定数を G とし、
μ=m1+m2m1m2M=m1+m2
として、以下の設問に答えよ。なお必要に応じて、時間微分を r˙ や ϕ˙ などと上付きドットを用いて表しても良い。
(1) 質点 1 および質点 2 の位置ベクトルを、m1, m2 および r 用いて表せ。
(2) 二つの質点の運動エネルギーの総和を μ, r および ϕ を用いて表せ。
(3) この系の重力ポテンシャルエネルギーを μ, M, r および G を用いて表せ。
(4) l≡r2ϕ が保存量となることを示せ。
(5) 動径成分 r についての運動方程式を、μ, M, r, l および G を用いて書き下せ。
(6) 系の全エネルギー E を μ, M, r, l および G で表せ。
(7) Eが最小値を取る時の E の表式を μ, M, l および G を用いて求め、その時の二つの質点の軌道の様子を述べよ。
(8) r=(r1,ϕ1) を質点 1 の位置ベクトルとする。
二つの質点のが重力的に束縛されお互いの周りを周回する条件と、その時の r1 と ϕ1 の関係を求めよ。
Kai
(1)
題意より、
m1r1+m2r2=0, r=r1−r2
であるから、
r1=m1+m2m2r, r2=−m1+m2m1r
を得る。
(2)
r のデカルト座標による成分を (x,y) とすると、
xyx˙y˙r˙2=rcosϕ=rsinϕ=r˙cosϕ−rϕ˙sinϕ=r˙sinϕ+rϕ˙cosϕ=x˙2+y˙2=r˙2+r2ϕ˙2
である。
よって、求める運動エネルギー K は、
K=21m1r1˙2+21m2r2˙2=21(m1+m2)2m1m22+m12m2r˙2=21m1+m2m1m2(r˙2+r2ϕ˙2)=21μ(r˙2+r2ϕ˙2)
となる。
(3)
求める重力ポテンシャルエネルギー U は、
U=−Grm1m2=−rGμM
である。
(4)
この系のラグランジアン L は、
L=K−U=21μ(r˙2+r2ϕ˙2)+rGμM
であり、
∂ϕ∂L∂ϕ˙∂L=0=μr2ϕ˙=μl
であるから、 ϕ に関する運動方程式より、
l が保存量であることがわかる。
(5)
上で求めたラグランジアン L から、
∂r∂L∂r˙∂L=μrϕ˙2−r2GμM=r3μl2−r2GμM=μr˙
であるから、 r に関する運動方程式は、
μr¨∴ r¨=r3μl2−r2GμM=r3l2−r2GM
となる。
(6)
E=K+U=21μ(r˙2+r2ϕ˙2)−rGμM=21μ(r˙2+r2l2)−rGμM
(7)
V(r)=21μr2l2−rGμM
とすると、
V′(r)=−μr3l2+r2GμM=μr3−l2+GMr
となるので、
r0=GMl2
とすると、 V(r) が最小値をとるのは r=r0 のときで、
V(r0)=−21l2G2μM2
となる。
E が最小値 E0 をとるのもこのときで、
E0=−21l2G2μM2
である。
(8)