■ケプラーと正多面体(その14)

 極座標(r,θ)を用いて

  r=r(θ)=a/(1+εcosθ)   (a>0,0<ε<1)

で表される曲線は原点を1つの焦点とする楕円で,εはその離心率となる.なぜなら,分母を払ってr+rεcosθ=a.また,r^2=x^2+y^2,rcosθ=xであるから

  x^2+y^2=a^2−2εax+ε^2x^2

  (1−ε^2)(x+εa/(1−ε^2))^2+y^2=a^2/(1−ε^2)

これは楕円の方程式である.

 2次曲線は円(e=0)として生まれ,楕円に育ち,放物線(e=1)で相転移して双曲線になる.漸近線のなす角度は最初鋭角だがだんだん大きくなり,180°になった時点(e=∞)で虚領域に入る.そして再び円に生まれ変わる.楕円の面積は有限,放物線の面積は∞,この考え方を延長すると双曲線の面積は虚数ということになる・・・というのがケプラーの原理である.射影平面上では,円錐曲線はただ1種類しかなく,双曲線・放物線・楕円などの区別はなく,どれも同種の曲線となる.

===================================

【1】ケプラーの法則

 16世紀前半,コペルニクスは地動説を提案しました.この地動説に基づき,天文学者ケプラーは,ティコ・ブラーエが資料として残してくれた20年にわたる膨大で精確無比な天文観測記録を試行錯誤的に数値計算し,25年も費やして火星の軌道を執拗に模索しました.ケプラーは後年,この苦心談を「マルス(火星)との悪戦苦闘」と表現していますが,とてつもない忍耐力と信念に支えられた偉大な業績です.

 しかし,その結果は「惑星の軌道は太陽を中心とする円軌道である」とするコペルニクスの地動説に反するものでした.そこで,彼は円軌道という前提に疑問をいだき,これに合う理論を求めてケプラーの法則に到達します.

<第1法則>惑星は太陽を焦点のひとつとする楕円軌道上を動く(1609年).

<第2法則>面積速度は一定である(角運動量保存則)(1609年).

<第3法則>公転周期の2乗は平均距離の3乗に比例する(1619年).

すなわち,惑星の軌道は完全無欠な円ではなく楕円であり,太陽はその一つの焦点の位置にあるとすることによって矛盾が解決されることを導き出したのです.日本ではまだ江戸時代が始まったばかりの時代です.

===================================

[Q1]ケプラーの3法則をもとにして力学を確立したのがニュートンですが,ニュートンの法則からケプラーの第2法則を導け(ニュートンからケプラーへ).

[A2]まず,準備としてxy平面の極座標(r,θ)を用いて

  r=r(θ)=a/(1+εcosθ)   (a>0,0<ε<1)

で表される曲線は原点を1つの焦点とする楕円で,εはその離心率となります.なぜなら,分母を払ってr+rεcosθ=a.また,r^2=x^2+y^2,rcosθ=xであるから

  x^2+y^2=a^2−2εax+ε^2x^2

  (1−ε^2)(x+εa/(1−ε^2))^2+y^2=a^2/(1−ε^2)

これは楕円の方程式です.

 この楕円は

  r(θ)=(r(θ)cosθ,r(θ)sinθ)

とパラメータ表示されますから,点P=r(θ)における接ベクトルは

 dr/dθ=(−asinθ/(1+εcosθ)^2,a(ε+cosθ)/(1+εcosθ)^2)

となることがわかります.以上の準備の下でケプラーの法則について考えてみましょう.

 rやθを時間tの関数として,’=d/dt,”=d^2/dt^2とすれば

  r(t)=(x(t),y(t))=(r(t)cosθ(t),r(t)sinθ(t))

  v(t)=(vx(t),vy(t))=(r’cosθ−rθ’sinθ,r’sinθ+rθ’cosθ)

  v’(t)=((r”−rθ’^2)cosθ−(rθ”+2r’θ’)sinθ,(r”−rθ’^2)sinθ+(rθ”+2r’θ’)cosθ+rθ’cosθ)

 ここで,

  er=(cosθ,sinθ)

  eθ=(−sinθ,cosθ)

すなわち,erはr方向の単位ベクトル,eθはそれと直交する単位ベクトルとすると

  v=r’er+rθ’eθ

  v’=(r”−rθ’^2)er+(rθ”+2r’θ’)eθ

 ニュートンの運動法則:F=mv’と万有引力の法則より

  m{(r”−rθ’^2)er+(rθ”+2r’θ’)eθ}=−GMm/r^2er

ですから,

  r”−rθ’^2=−GMm/r^2

  rθ”+2r’θ’=0

が成り立っていなければなりません.

 角運動量Lは位置ベクトルrと運動量mvの外積で定義され,その大きさは|r×mv|=mr^2θ’で表されるのですが,

  d(mr^2θ’)/dt=mr(rθ”+2r’θ’)=0

より,

  mr^2θ’=c(一定)

これは角運動量保存則にほかなりません.

 一方,面積速度は

  S(θ)=1/2∫{r(θ)}^2dθ,dS/dθ=1/2{r(θ)}^2

より

  dS(θ)/dt=1/2r^2dθ/dt=1/2r^2θ’

角運動量が時間によらないわけですから,面積速度dS/dtも時間によらず一定

  c/(2m)

となるわけです.

−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−

[Q2]ニュートンの法則からケプラーの第3法則を導け(ニュートンからケプラーへ).

[A2]近似的に,太陽の周囲を公転する惑星の軌道が円であると仮定する.公転軌道半径R,公転角速度ω,公転周期Tとすると,

  mRω^2=GMm/R^2

  ω=√GM/R^3

  T^2=(2π/ω)^2=4π^2R^3/GM

===================================