中心力を受ける粒子の軌道方程式
📂古典力学中心力を受ける粒子の軌道方程式
中心力を受ける粒子の軌道方程式
中心力を受ける質量がmの粒子の運動方程式を極座標系で表すと以下の通りだ。
mr¨=F(r)r^
F(r)は粒子に作用する中心力だ。極座標系での加速度は以下のようになる。
r¨=a=(r¨−rθ˙2)r^+(2r˙θ˙+rθ¨)θ^
したがって、運動方程式(eq1)を成分ごとに分けて書くと以下のようになる。
m(r¨−rθ˙2)m(2r˙θ˙+rθ¨)=F(r)=0
この時、r1dtd(r2θ˙)=2r˙θ˙+rθ¨なので、上の二番目の式を以下のように書き換えることができる。
dtd(r2θ˙)=0
これはr2θ˙が定数であるという意味だ。今、その定数をlとしよう。
r2θ˙=constant=l
一方、中心力によって運動する粒子の角運動量の大きさは以下の通りである。
L=mr2θ˙(4)
したがって、(eq3), (eq4)により次の式を得る。
∣l∣=mL=∣r×v∣
が成立する。したがって、lは粒子の単位質量あたりの角運動量と解釈できる。lが定数であるということは、中心力によって動く粒子の角運動量が保存されるというすでに知られている結果だ。そこで、時間に関係ない空間での粒子の軌道を求めるために、以下のように新しい変数を導入しよう。
u=r1
すると、(eq3)によりθ˙=lu2が成り立つ。これで、r˙、r¨を求めると以下のようになる。
r˙=dtd(u1)=dud(u1)dtdu=−u21dtdu=−u21dθdudtdθ=−u21θ˙dθdu=−ldθdu
r¨=−ldtd(dθdu)=−ldθd(dθdu)dtdθ=−lθ˙dθ2d2u=−l2u2dθ2d2u
これらを(eq2)に代入すると、以下の式を得る。
⟹−ml2u2dθ2d2u−ml2u3dθ2d2u+u=F(u1)=−ml2u21F(u1)
ここで、力Fが重力のように距離の逆二乗に比例する場合、軌道は楕円になり、これを惑星の運動に適用したものがケプラーの第一法則である。