軌道力学

軌道力学（きどうりきがく）は、弾道学と天体力学の応用で、ロケットや宇宙船の軌道に関する現実的な問題を解決するための学問である。これらの物体の軌道は、ニュートン力学と万有引力から計算することができる。軌道力学は、宇宙探査ミッションの設計や制御の基本原理である。天体力学は、重力の下での、軌道力学よりも広範な領域を扱い、宇宙船も恒星系、惑星、衛星、彗星等を含めた天体も、どちらも対象となる。軌道力学は、軌道マヌーバ、軌道平面の変更、惑星間移動も含めた宇宙船の軌道に対象を絞っており、ミッションの計画者が宇宙機の推進を予測するために用いられる。一般相対性理論は、ニュートンの法則より正確に軌道を計算し、高い精確さが必要な場面や太陽近傍等の重力が非常に強い環境では必須である。
歴史.
20世紀に宇宙飛行が実現するまで、軌道力学と天体力学の間には、ほとんど差異がなかった。時間の関数として位置を決定するケプラーの問題を解くために用いるような基本技術は、どちらの学問領域でも同じだった。さらに、この学問分野の歴史はほぼ共有していた。
ヨハネス・ケプラーは、高い正確性で惑星の軌道のモデル化に成功した最初の人物であり、1605年にケプラーの法則を発表した。アイザック・ニュートンは、天体の運動のより一般的な法則を1687年の著書『自然哲学の数学的諸原理』の中で発表した。
実用的な技術.
経験則.
以下の経験則は、標準的な前提の下で天体力学で近似できる状況にとって有用である。議論されている特定の例は、惑星の周囲を公転する衛星であるが、この経験則は恒星の周囲の小天体のような他の状況にも適用することができる。
軌道力学の法則の結果は、時として直観と相容れないことがある。例えば、同じ円軌道上の2機の宇宙船がドッキングしようとする場合、その位置が非常に接近していない時に、後ろの宇宙船は速度を速めるために単純にエンジンを吹かすことはできない。そうすると軌道の形が変化し、ターゲットと出会うことができない。ドッキングするための1つの方法は、速度を下げるために逆向きにエンジンの噴射を行い、その後、低い円軌道に戻すために再度噴射を行う。低軌道は高軌道よりも速度が速いため、後ろの宇宙船は追いつくことが出来る。3度目の噴射で、先行する宇宙船の軌道と交わり、後ろから接近できるような楕円の軌道に移行する。
標準的な前提が適用できないようなレベルであれば、実際の軌道は計算したものからずれることになる。例えば、大気の抗力は、地球軌道にある物体について複雑化要因になり得る。これらの経験則は、連星系等の同程度の質量の2つかそれ以上の物体に適用する際には、不正確になる。惑星のような大きな物体にとっては、古典力学と一般相対性理論の差異も重要になる。
天体力学の法則.
天体力学の基本法則は、ニュートンの万有引力の法則とニュートンの運動の法則であり、ニュートンの開発した微分積分学がその計算のための重要な数学的ツールになる。
ケプラーの法則は、周回する天体が中心の天体からの重力のみを受けていると見なせる場合には、ニュートンの法則から導くことができる。推力が働く場合、ニュートンの法則は適用できるが、ケプラーの法則は成り立たなくなる。推力が止まると、結果として軌道は変わるが、再びケプラーの法則が適用できるようになる。ケプラーの3法則は、次のとおりである。
脱出速度.
脱出速度の公式は、以下のように簡単に導くことができる。全ての宇宙船の単位質量当たりのエネルギーは、単位質量当たりの位置エネルギーと単位質量当たりの運動エネルギーという2つの成分から成り立っている。単位質量当たりの位置エネルギーは、惑星の質量"M"と関係があり、次の式で与えられる。
formula_1
一方、単位質量当たりの運動エネルギーは、次の式で与えられる。
formula_2
エネルギー保存の法則から、単位質量当たりの合計の軌道エネルギー
formula_3
は、中心の天体から宇宙船までの距離formula_4には依存しない。これより、天体は、この値が負ではない時に無限大formula_4に達し、これは、以下を意味する。
formula_6
地球表面からの脱出速度は、約11km/sであるが、太陽の重力があるため、物体を無限の距離まで送るにはこの速度では不十分である。地球から太陽までの距離で地球の近傍以外の場所から太陽系の外への脱出速度は、約42km/sであるが、地球から打ち上げられる宇宙船は、方向が同じであれば地球の公転速度も利用することができる。
自由軌道の公式.
軌道は円錐曲線であるため、所与の角度の天体までの距離の公式は、極座標系での曲線の公式と一致し、次のようになる。
ここで、μは標準重力パラメータ、Gは重力定数、m1とm2は天体1と天体2の質量、hは天体1に対する天体2質量当たりの角モーメントである。パラメータθは、真近点角、"p"は半通径、"e"は軌道離心率として知られる値であり、全て6つの軌道要素から導くことができる。
円軌道.
中心の天体の重力が独占的な全ての軌道は、楕円軌道となる。この特殊なケースは円軌道であり、これは軌道離心率が0の楕円軌道である。質量"M"の重力中心からの距離"r"における円軌道上の天体の速度は、次の公式で表される。
ここで、formula_11は重力定数であり、次の値と等しい。
この公式を適切に使うためには、単位は一貫していなければならず、例えば"M"はkg、"r"はmで表されていなければならない。答えはm/sの単位で得られる。
"GM"の大きさは、しばしば標準重力パラメータと呼ばれ、太陽系のそれぞれの惑星や衛星で異なった値である。
円軌道の速度が既知になると、2の平方根をかけることで簡単に脱出速度が求められる。
楕円軌道.
軌道離心率が0<e<1であれば、自由軌道の方程式の分母が真近点角θに伴って変化するが、0になることはなく、必ず正の値である。そのため、相対的な位置ベクトルは有界であり、近点rpにおける最小値は、以下の式で与えられる。
最大値rは、θ=180°の時である。この点は遠点と呼ばれ、raで示される動径座標は、次のようになる。
近点Pから遠点Aまでの距離を2aとすると、以下の方程式が成り立つ。
上記の方程式を展開して、次の値が得られる。
aは、楕円軌道の軌道長半径である。rについて解くと、次のようになる。
軌道周期.
標準的な前提の下で、楕円軌道を公転する物体の軌道周期 (formula_18) は、次のように計算される。
ここで、
ここから次の結論が得られる。
速度.
標準的な前提の下で、楕円軌道を公転する物体の軌道速度(formula_23)は次のように計算される。
ここで、
双曲線軌道の速度の方程式は、+ formula_28となるか、"a"が負の場合は上式となる。
エネルギー.
標準的な前提の下で、楕円軌道を公転する物体の質量当たりの軌道エネルギー(formula_29)は負となり、この軌道の軌道エネルギー保存の方程式は、次のような形になる。
ここで、
ここから次の結論が得られる。
ビリアル定理を用いると、次のことが分かる。
放物線軌道.
軌道離心率が1と等しい場合、軌道の方程式は次のようになる。
ここで、
真近点角θが180°に近づくと、分母は0に近づき、"r"は無限大に発散する。従って、"e"=1の時の軌道エネルギーは0であり、次の式で表される。
ここで、
言い換えると、放物線軌道上の任意の点の速度は、次のようになる。
双曲線軌道.
軌道離心率がe>1であれば、双曲線軌道となり、軌道方程式は次のようになる。
この系は2つの対称な曲線で構成されている。周回する天体がそのうちの1つを占め、もう1つはその空の数学的な像になる。明らかに、上記の方程式の分母は、cos"θ" = -1/"e"となると0になり、この時の真近点角の値を次のように示す。
真近点角が"θ"∞に近づくと、半径方向距離が無限大になる。"θ"∞は、「漸近線の真近点角」として知られる。"θ"∞を90°から180°の間とすると、sin2"θ"+cos2"θ"=1という三角関数の性質から、次のように書ける。
エネルギー.
標準的な前提の下で、双曲線軌道の質量当たりの軌道エネルギー(formula_29)は0より大きく、軌道エネルギー保存の法則により、次の式が成り立つ。
ここで、
双曲線過剰速度.
標準的な前提の下で、双曲線軌道を公転する天体は、無限遠点で双曲線過剰速度(formula_50)と呼ばれる軌道速度に到達する。この速度は、以下の式で表される。
ここで、
双曲線過剰速度は、質量当たりの軌道エネルギーに関係がある。
軌道の計算.
ケプラーの方程式.
主に歴史的に使われてきた軌道を計算する方法の1つは、ケプラーの方程式である。
ここで、"M"は平均近点角、"E"は離心近点角、formula_56は軌道離心率である。
ケプラーの公式では、近点から真近点角formula_57に至るまでの時間は、2つのステップによって求められる。
逆に与えられた時間の離心近点角を求めるのはより難しい。ケプラーの方程式はformula_59に対して超越的で、つまりformula_59について代数的に解くことはできない。ただし、反転させて解析関数的に解くことはできる。
全ての実数formula_64に対して適用できるケプラーの方程式の解は、以下のとおりである。
formula_65
この値を求めることで、次の式が出る。
formula_66