ランベルト問題シミュレーター 戻る
宇宙工学

ランベルト問題シミュレーター

2点(出発点 r₁ と到達点 r₂)を結ぶ遷移軌道を求める「ランベルト問題」を可視化するツールです。最小エネルギー楕円の遷移時間、Hohmann遷移との比較、各噴射点のΔvと合計Δvがリアルタイムで分かり、LEO→GEOなど軌道遷移の感覚を掴めます。

パラメータ設定
出発点距離 r₁
km
中心天体中心から出発点までの距離(LEO初期軌道)
到達点距離 r₂
km
到達目標までの距離(既定はGEO静止軌道半径)
遷移角 Δθ
°
出発点と到達点の中心角。180°でHohmann遷移
中心天体重力定数 μ
km³/s²
μ=GM。地球:398600、月:4903、太陽:1.327e11
遷移方向
同じ2点間でも2解存在する
計算結果
弦長 c (km)
最小エネルギー長半径 a_min (km)
遷移時間 (h)
Hohmann遷移時間 (h)
1段目 Δv (km/s)
合計 Δv (km/s)
軌道遷移の可視化 — 中心天体・初期軌道・目標軌道・遷移楕円

中央が中心天体。青の円が初期円軌道 r₁、緑の円が目標円軌道 r₂、橙の楕円が遷移軌道。短経路は実線、長経路は点線で描画され、白の探査機が遷移軌道を周回します。

遷移時間 vs 遷移角 Δθ (30〜330°)
合計Δv vs 軌道半径比 r₂/r₁
理論・主要公式

$$c = \sqrt{r_1^{2} + r_2^{2} - 2\,r_1 r_2 \cos\Delta\theta}, \qquad s = \frac{r_1 + r_2 + c}{2}, \qquad a_{\min} = \frac{s}{2}$$

弦長 c、半周長 s、最小エネルギー楕円の長半径 a_min。a_min はランベルト方程式の特殊解で、与えられた2点を結ぶ全楕円の中でエネルギー(−μ/2a)が最小になる軌道。

$$t_{\min} = \frac{1}{\sqrt{\mu}}\,\frac{s^{3/2}}{\sqrt{2}}\,\left[\,\pi - \arcsin\sqrt{1-\frac{c}{s}} + \sqrt{\frac{c}{s}\!\left(1-\frac{c}{s}\right)}\,\right]$$

最小エネルギー遷移軌道の飛行時間(Lagrangeの定理)。μ=GMは中心天体の重力定数。c, s, μ を入れれば一意に決まる。

$$\Delta v_{\text{total}} = \left|v_p - \sqrt{\mu/r_1}\right| + \left|\sqrt{\mu/r_2} - v_a\right|, \qquad v_{p,a} = \sqrt{\mu\!\left(\tfrac{2}{r_{p,a}} - \tfrac{1}{a_H}\right)}$$

Hohmann遷移(180°)の合計Δv。a_H = (r₁+r₂)/2 で、v_p は近点速度、v_a は遠点速度。両端の円軌道速度との差の和が必要噴射量となる。

ランベルト問題と軌道遷移

🙋
「ランベルト問題」って初めて聞きました。要は何を解いてるんですか?
🎓
ざっくり言うと、「中心天体まわりで、出発する点 P1 と到達したい点 P2、そしてそこに着くまでの時間を決めたら、その間どんな軌道を飛べばいいか」を解く問題なんだ。1761年にスイスのランベルトさんがまとめた古典的問題で、今でも月・火星探査の打上げ計算、宇宙ステーションのドッキング、デブリ回避まで、宇宙ミッションの根本にずっと座ってる。要素は3つ:r₁、r₂、そして両者をはさむ遷移角 Δθ。これで遷移軌道(楕円)の長半径 a と離心率 e が一意に決まる。
🙋
2点を結ぶ軌道なら無数にありそうですけど、そんなに簡単に1つに決まるんですか?
🎓
いいツッコミ。実はランベルトが示したのは「飛行時間 t を1つ決めると軌道が一意に決まる」という性質なんだ。だからツールでは時間を直接入れる代わりに、計算しやすい「最小エネルギー解」を表示してる。これは a が最小になる楕円で、s = (r₁+r₂+c)/2 として a_min = s/2、c は弦の長さ。ここから飛行時間が Gauss の公式で出る。スライダーで r₁=6678(LEO)、r₂=42164(GEO)、Δθ=180° にしてみて。a_min=24421 km、遷移時間 5.27時間、合計Δv 3.89 km/s——これは教科書に載っているLEO→GEOホーマン遷移の正解値そのものだよ。
🙋
短経路と長経路って何ですか?ボタンで切り替わるみたいですけど…。
🎓
同じ2点でも、行き方は実は2通りある。Δθ<180°側を通る「短経路」と、反対側を回る「長経路」だ。短経路は早く着くけどΔvが大きく、長経路は時間はかかるけど位相合わせがしやすい。例えば国際宇宙ステーションへの「FastTrack(6時間ドッキング)」は短経路、ソユーズの古い「Standard(2日ドッキング)」は長経路寄り、というイメージ。Δθ=180°のときだけ両者が一致して、それがHohmann遷移という特別な解なんだ。
🙋
そのHohmann遷移はなんで「最効率」と呼ばれるんですか?
🎓
同一平面の円軌道どうしをつなぐ場合、両端で接線方向に2回噴射するHohmannが「合計Δv最小」になることが理論的に証明されているからだよ。下の「合計Δv vs r₂/r₁」グラフを見ると、半径比を上げるほどΔvが増えていくのが分かる。ただし所要時間は r₂ の3/2乗で効くから、火星みたいに遠い目的地だと「Hohmannでも258日かかる→もっと速いがΔvを払う非ホーマン遷移」とトレードオフ判断するんだ。実務ではこの判断のために、出発日×到着日のマトリクスで何千通りもランベルト問題を解いた「ポークチョップ図」を作って最適化する。
🙋
μ(重力定数)を変えると、月や火星まわりの遷移も計算できるんですか?
🎓
そう、μ を切り替えれば「どの天体まわりか」を変えられる。地球 398600、月 4903、火星 42828、太陽 1.327e11 km³/s²。例えば月まわりで r₁=1838 km(低月軌道)、r₂=8000 km、Δθ=180° と入れると、Hohmann遷移時間が約3.4時間、合計Δv 約0.6 km/s と出てくる。アポロ計画の月遷移軌道(地球発・月着)はもっと大規模だけど、まずは中心天体ひとつの「2体問題」として近似するこの計算が出発点になる。実機ミッションでは更に月の重力影響圏や太陽摂動を入れて3体・多体問題に拡張していくんだ。

よくある質問

ランベルト問題とは、中心天体まわりの2点(出発点P1と到達点P2)と、その間の飛行時間が与えられたときに、両点を結ぶ唯一の遷移軌道(速度ベクトルと軌道要素)を求める問題です。1761年にJ.H.ランベルトが定式化し、現代の宇宙ミッション設計の基礎になっています。月・惑星探査の打上げウインドウ計算、ランデブー、軌道変更など、ほぼ全ての軌道計画問題がこの形式に帰着します。本ツールでは特に最小エネルギー解(a_min=s/2)と、180°遷移の特殊ケースであるHohmann遷移を扱います。
ランベルト問題は遷移角 Δθ に対し原則として2解持ち、それぞれ「短経路(short way、Δθ<180°側)」と「長経路(long way、Δθ>180°側)」と呼ばれます。短経路は飛行距離が短く到達が速いですが、要求Δvが大きくなる傾向があります。長経路は迂回するため時間がかかる一方、出発時の方向自由度が大きく、ランデブー時の位相合わせに使われます。Δθ=180°の場合は両者が一致し、この特殊解がHohmann遷移です。本ツールの「遷移方向」セレクタで切り替えられます。
Hohmann遷移は、同一平面内で2つの円軌道(半径 r₁ と r₂)を最小Δvでつなぐ遷移楕円で、遷移軌道の近点を r₁ に、遠点を r₂ に一致させた180°遷移です。両端2回の接線方向加速のみで完結し、ほぼ任意の Δv 配分の中で合計が最小になることが証明されています。代表例のLEO→GEO(r₁=6678km, r₂=42164km)では合計Δv≈3.89km/s、遷移時間約5.27時間と本ツールでも確認できます。ただし所要時間は長いため、火星探査などでは時間を節約する非ホーマン遷移が選ばれることもあります。
ミッション設計では、出発日 t₁ と到着日 t₂ を多数組み合わせ、各ペアに対してランベルト問題を解いてΔvを求める「ポークチョッププロット(豚カツ図)」を作成します。地球出発時のC3(脱出エネルギー)と火星到着時のV∞(双曲線過剰速度)を等高線で表示し、ロケット能力・探査機質量・科学観測ウインドウの制約をすべて満たす最適な打上げ日を選定します。NASAのMariner・Voyager・Mars Reconnaissance Orbiterなど全ての惑星間探査でこの手法が標準的に使われています。

実世界での応用

静止衛星投入(GTO→GEO):商業通信衛星はロケットで一旦GTO(Geostationary Transfer Orbit、近点LEO・遠点GEO)に投入され、衛星自身の遠点キック(アポジエンジン噴射)で円形のGEO軌道に乗り換えます。これがまさに本ツールの既定条件(r₁=6678→r₂=42164、Δθ=180°)で示される LEO→GEO Hohmann遷移で、典型的Δvは約3.89 km/s、遷移時間は約5.3時間です。打上げプロバイダ(H-IIA、Falcon 9、Ariane 5)の能力曲線はこのΔvを基準に作られています。

惑星間探査(Mars Reconnaissance Orbiter等):火星探査では地球公転軌道(r₁≈1.496e8 km)と火星公転軌道(r₂≈2.279e8 km)を結ぶ太陽中心のホーマン遷移を基本とし、ランベルト解で出発日・到着日の組み合わせを多数評価して打上げウインドウを決定します。Hohmann解で約258日(μ=1.327e11 km³/s²の太陽中心問題)、Δv約5.6 km/sが理論最小値。実機では火星のスイングバイ角度や着陸サイトの可視時間を加味した非ホーマン解が選ばれることが多いです。

ISSランデブー・ドッキング:ソユーズ宇宙船やCrew Dragonが国際宇宙ステーション(高度約400 km)へ追いつく際、6時間ドッキング・2日ドッキングなど複数の遷移シナリオがあり、いずれもランベルト解で位相と接近速度を設計します。短時間ドッキングではISSの位置を高精度に予測し、追跡誤差をリアルタイムで補正する複数回ランベルト解が用いられます。

軌道上サービス・デブリ除去:故障衛星の救援、燃料補給、デブリ捕獲ミッションでは、サービス衛星が複数の目標を巡回する必要があり、各目標へのアクセスは個別のランベルト問題として解かれます。多目標巡回ではTSP(巡回セールスマン問題)的な最適化と組み合わせ、総Δvと所要時間を同時最小化する経路を遺伝的アルゴリズムやAIで探索する研究が活発です。

よくある誤解と注意点

まず最も多い誤解は、「Δθ=180°のときがいつも最効率」というもの。これは円軌道→円軌道(共面・同心)の場合に限った話で、軌道が傾いていたり離心率を持つ場合、最適遷移角は180°から外れます。さらに、目標との位相合わせ(出発時に到着点が正しい位置に来るタイミング)を考えると、Hohmannは数年に1回しか機会がない(火星なら約26ヶ月周期の「火星会合周期」)ことが多い。実機ではこの「打上げウインドウ」の制約が決定的で、Δv最適性を多少犠牲にしても日程を優先することがほとんどです。本ツールの計算結果は「理想的なΔv下限」として参照してください。

次に、「ランベルト解は2体問題」であることを忘れる落とし穴。本ツールも含めて、ランベルト方程式は中心天体ひとつの重力場だけを考えます。実際には地球→月遷移なら月の重力影響圏(SOI)内では月中心の2体問題に切り替えるパッチドコニック近似、惑星間遷移では出発惑星・太陽・到着惑星の3区間でそれぞれランベルト解を継ぎ足すのが一般的です。さらに高精度ミッションでは太陽輻射圧、J2扁平度項、第三体摂動(月・木星)まで考慮した数値積分(一般摂動法)が必要で、ランベルトはあくまで「初期推定値」を与えます。

最後に、「Δv=0なら燃料は要らない」と思い込むこと。本ツールが出すΔvは「速度を変えるのに必要な瞬時的増分」で、燃料消費はツィオルコフスキー式 Δv = Iₛₚ·g·ln(m₀/m₁) で質量比に変換しなければなりません。例えばΔv=3.89 km/s をIₛₚ=320秒(ヒドラジン2液系)のスラスタで実現すると、初期質量の約71% が燃料になります。商業衛星の打上げ重量の半分以上が「軌道変更のための推進剤」なのはこのため。Δvを1割減らせるミッション設計は、それだけで搭載機器を倍にできる級の価値があります。

使い方ガイド

  1. 初期軌道半径(km)と最終軌道半径(km)を入力します。例えばLEO(高度200km、地心距離6578km)からGEO(高度35786km、地心距離42164km)への遷移を想定してください。
  2. 遷移角度(度)を0~360の範囲で指定します。180度はHohmann遷移、180度未満は短弧遷移で到達時間が短くなりますがΔvが増加します。
  3. 天体の重力定数μ(km³/s²)を設定します。地球は3.986×10⁵、月は4.905×10³です。シミュレーターが弦長、最小エネルギー長半径、遷移時間、必要なΔvを自動計算します。

具体的な計算例

LEO(6578km)からGEO(42164km)へ遷移角度180度で移行する場合:弦長c≈47300km、最小エネルギー楕円の長半径a_min≈24371km、遷移時間≈5.3時間、Hohmann遷移時間≈5.3時間(同一)、1段目Δv≈2.45km/s、合計Δv≈3.88km/s。遷移角度120度の短弧遷移なら遷移時間≈2.8時間に短縮されますが、合計Δv≈4.20km/sに増加します。

実務での注意点

  1. Hohmann遷移は最小エネルギー需要ですが、遷移時間が長いため、衛星寿命やミッション期間に制約がある場合は短弧遷移(180度未満)を選択する必要があります。
  2. 重力定数μの値を誤入力するとすべての計算結果が無効になります。地球3.986×10⁵、火星4.283×10⁴を確認してください。
  3. 遷移角度0度や360度に近い値ではランベルト問題が特異点に近づき、数値計算の精度が低下します。実務では5~355度の範囲内で設計してください。