このページを印刷される方はこちらのバージョンをご利用下さい。ブラウザーでは見にくいのですが印刷は鮮明です。

惑星探査機の軌道と飛行速度

　惑星探査機の航行速度の計算例を示します。惑星探査機の飛行は、出発時に所定の速度を得るためにロケットエンジンを用いるが、それ以後はロケットエンジンを使用しない慣性飛行であることに注意。その様に最初にロケット燃料を使い切る飛行をするのは、できるだけ速い時期に最大速度・最高エネルギー値に到達した方が飛行時間を短縮できるからです。それはここでの説明を読めば了解できる。以下の計算では、

とする。また、簡単化して惑星は全て円軌道を回転するとする。そして各惑星の重力ポテンシャルの寄与は全て無視して、太陽の重力ポテンシャル場の中での飛行を考える。惑星の質量は太陽の質量の1000分の1以下だから、ここでの考察のようにおおざっぱな値を見積もるときには無視しても問題ない。そのため計算値は実際の値とは少し異なる。

１．各惑星の公転半径位置で円運動する軌道速度ｖ₁

高校物理�Uで習うように公転軌道速度ｖ₁は

となる。

これは当然の事ですが各惑星の公転軌道速度と同じです。

２．各惑星軌道上に於ける太陽系脱出速度ｖ₂

　惑星公転軌道上から出発して太陽系を脱出するための最低速度ｖ₂は、無限遠点での速度がゼロになるとしたエネルギー保存則を解いて求めればよい。

となる。これは前項で求めたｖ₁の２^1/2倍になる。

この速度で各惑星を出発した探査機の飛行速度ｖを太陽からの距離ｒで表すと

となる。これは当然のことであるが、上記の脱出速度ｖ₂(ｒ)と全く同じ関数形になる。

　このとき、各惑星を上記の太陽系脱出速度で出発するとき、その出発方向はどの方向でも良いことに注意、どの方向に出発しても太陽を焦点とする放物線軌道を描いて無限遠の彼方へ到達できる。
　この速度以下で出発した場合に描く楕円軌道の場合は、出発の方向によって太陽中心に対して持つ角運動量が異なってくる。そのとき、動径方向の最大値に達したときにも動径に垂直な方向の速度成分を最初持っていた角運動量成分に応じて持たねばならないので、太陽からの最遠到達距離が異なってくる。
　これに対して、上記の放物線軌道の場合は、到達動径距離の最大値は太陽から無限遠の彼方だから、出発時に持っていた角運動量成分が異なっていて、無限の彼方での角運動量が有限でも、無限遠での話しだから動径に垂直な速度成分はゼロと見なせる。もちろん動径方向の速度成分もゼロに成って良いのだから、最初持っていた運動エネルギーを全て位置エネルギーの獲得に使えることになる。そのため放物線軌道や双曲線軌道が取れる初速度で出発した場合は、太陽に対して任意の方向に向かって出発しても最終的に無限遠の彼方へ到達できる。次に上記脱出速度（42.1km/s）で様々な角度に打ち出したときに、たどる放物線軌道の様子を図示する。

次に、太陽系脱出速度よりも大きな速度45km/sで打ち出したときに取る双曲線軌道を示す。

　上記の軌道の中で、ロケット燃料が最も節約できるのは、地球の公転運動速度を最大限利用できるθ₀＝90°の場合です。そのため実際には、θ₀＝90°の方向に出発して、惑星軌道に到着する時に目的の惑星がちょうどその位置にいるような地球−惑星の位置関係を選んで出発する事になる。

３．地球-外惑星ホーマン軌道の地球出発速度ｖ₃と外惑星軌道到達時の速度ｖ₄

　ホーマン軌道は、地球軌道から出発するとき地球の公転運動の速度を最大限利用できるために、ロケットエンジンの燃料が最も節約できる軌道です。ホーマン軌道を取るときの、地球出発速度ｖ₃と惑星到着時の速度ｖ₄を角運動量保存則（面積速度一定の法則）とエネルギー保存則から求める。ホーマン軌道については別稿「質点の二次元運動」４．（３）（ｂ）で述べたように以下の式が成り立つ。ただしそこでの比例定数ｋは今の場合ＧｍＭになる。

また、ホーマン軌道で飛行中の速度ｖ₅を太陽からの距離ｒの関数で表すと

となる。このとき当然　ｖ₅（ｒ_地）＝ｖ₃（ｒ_地）、ｖ₅（ｒ_惑）＝ｖ₄（ｒ_惑）　が成り立つ。

４．地球から外惑星軌道に到達する為の最小出発速度ｖ₆

　ホーマン軌道は探査機の惑星間移動に良く利用されるが、外惑星軌道の位置に到達するためだけの軌道としては、エネルギー最小の軌道でもなければ、最短時間の軌道でもない。エネルギー最小で最短時間で外惑星軌道位置まで到達する軌道は、外惑星軌道半径の位置と太陽の位置が両端になるような、離心率が１に近い平たい楕円軌道です。つまり、角運動量の成分＝ゼロで打ち出す。惑星探査機が地球を出発するときの速度ｖ₆は、そのような軌道の地球軌道半径位置での速度になる。これは、前項のエネルギー保存則で　ｖ₃＝ｖ₆、ｖ₄＝0　とした場合である。

途中の飛行速度ｖ₇を太陽からの距離ｒの関数で表すと

となる。　

　ただしこの軌道では、地球の公転方向に対して垂直に出発する事になるので地球の公転運動が探査機の加速に利用できない。そのため、出発速度まで加速するのにより多くのロケット燃料が必要になり、実用的ではない。
　実際、ボイジャー2号は地球の公転軌道に沿った方向に出発している。そして、木星に近づいたときには木星の引力を利用して木星に到達してスイングバイで増速した。

５．地球から外惑星軌道に到達する最小出発速度ｖ₆＝ｖ₀で、打ち出し角度を変えて出発した場合

　エネルギー最小で最短時間で外惑星軌道位置まで到達する軌道は角運動量成分＝ゼロで打ち出されたときの軌道で、外惑星軌道半径の位置と太陽の位置が両端になるような、離心率が１に近い平たい楕円軌道になるが、出発時に太陽中心に対して角運動量成分を持って出発すると、動径に垂直な速度成分ｖ_Φが必ず残るので動径方向の到達距離ｒが小さくなってしまう。
　下図は木星に到達するための最低速度37.8km/sで出発し、角度θ₀を様々に変えたときの軌道の様子を示す。出発時の角度θ₀が増えるにつれて、太陽からの最遠到達距離ｒ_maxが減少してくることに注意。図の描き方は別稿「質点の二次元運動」３．(３)(ｄ)で説明している。数式のグラフ描画ソフトがあれば６．ｒとφで述べる極座標表示の関数をそのまま描いても良い。

　地球Ｅの軌道半径＝ｒ₀の位置から角運動量成分＝0（つまりθ₀＝０）、出発速度ｖ₀で打ち出されたとき、太陽から最も離れるときの距離をｒ₁(今の場合木星の軌道半径＝5.2天文単位）とする。θ₀＝0の時の最遠到達距離ｒ₁は、エネルギー保存則から明らかなようにｖ₀、ｒ₀と次式で関係づけられる。

　次に、同じ出発速度ｖ₀だか出発の角度θ₀を変えたときに描く楕円軌道の、太陽Ｓから最も離れたときの距離をｒ_max、太陽Ｓと楕円の焦点Ｆとの距離をｒ_F、太陽に最も近づいたときの距離を離ｒ_minとして、これらの値を求めてみる。探査機の飛行中には、角運動量とエネルギーが保存されるので

が成り立つ。別稿「質点の二次元運動」３．（３）で説明したように最遠到達距離ｒ_max、及び最近距離ｒ_minではｖ_r＝0だから、ｖ_r＝0として上式を連立させて解くと、ｒ_max、ｒ_minはＬ、Ｅを用いて

と表される。また、三角形ＥＳＦに三角関数の余弦定理を用いると太陽Ｓと楕円の焦点Ｆとの距離をｒ_Fは

と表される。太陽Ｓから見た出発地点Ｅと焦点Ｆとのなす角φ_FEは次式から求まる。

　また、上記の二つの保存則をｖ_φ(ｒ)、ｖ_r(ｒ)の連立方程式と見なして解けば

が得られ、初期値(ｒ₀、ｖ₀、θ₀)及びｒとの関係式が得られる。

以下にθ₀＝67度の場合を例示す。

６．極座標表示に於ける動径ｒと方位角φの関係

　５．で用いた楕円軌道図を用いて説明する。

別項「二次曲線の性質」で説明した楕円軌道に関して成り立つ一般的な性質を用いると、標準形表現におけるａ、ｃはｒ_max、ｒ_minと以下の関係で結びつけられる。

また、楕円の離心率ｅと焦点Ｓと準線との距離ｆは、上記ａ、ｃと次のように結びつれられる。

楕円の極座標表示をもちいると、太陽からの距離ｒは方位角φ(上図参照)を用いて以下の様に表される。

φ測定の基準線SOは探査機を打ち出した方位角SEから（π−φ_FS）だけずれているが、φ_FSの値は

から求めることができるから、打ち出し点から測定した方位角で表すこともできる。

　以上の話では、軌道上での時間の関係は説明しませんでした。時間の関係を取り扱うには高校数学では習わない積分公式を利用しなければならないからです。
　この当たりは、例えばランダウ、リフシッツの理論物理学教程「力学」東京図書などに見通しよく解説されている。　

７．ボイジャー2号の飛行速度ｖ

　ボイジャーの軌道速度を太陽からの距離の関数として表すにはエネルギー保存則を用いればよい。このとき、重要なことは飛行速度の大きさは進行する方向や軌道の形によらず太陽からの距離のみの関数になることと、各惑星近傍でのスイングバイによる増速時以外は全て慣性飛行と見なせる事です。
　その為、各惑星の位置でスイングバイにより得た速度を初速度として計算したグラフを繋いでゆけばよい。つまり、各惑星によるスイングバイで得た速度を初速度ｖ₀とし、それ以後次にスイングバイする惑星に到達するまでの速度ｖの変化を太陽からの距離ｒによってあらわすと、エネルギー保存則より

となる。(別項「惑星探査機のスイングバイ航法」参照)。

　ちなみにNASAのHP(http://www2.jpl.nasa.gov/basics/bsf4-1.html)にボイジャー２号の実際の飛行速度が掲載されている。（ただし、この図のSOLAR SYSTEM ESCAPE VELOCITYのグラフは間違っている）

　NASAによると36km/sで地球を出発したことになっているが、これは地球の引力圏を脱出するときの値であるから、実際にはもっと速い速度（下図は39.2km/s）で出発したであろう。また、各惑星の近傍では惑星の重力のために速度が大きく変動している。NASAの図に最もマッチするような数値を与えて上記の方法で飛行速度のグラフを書いてみると次のようになる。

この図は以下の値を与えて描いたものであるが、実際の飛行速度と良く一致している。

このグラフに前述３〜４の結果ｖ₃(ｒ)〜ｖ₆(ｒ)を重ねると下図になる。