上昇軌道設計と最適化

はじめに

ロケットの上昇軌道は、地上の発射台から目標軌道に到達するまでの飛行経路であり、その設計はミッション成功の根幹を成す。上昇軌道の最適化は、限られた推進剤で最大のペイロードを軌道に投入するための重要な工学問題だ。

軌道設計には多数の制約が同時に存在する。構造荷重の制限、空力加熱、飛行安全の地上リスク制約、エンジンのスロットル範囲、さらには再使用ロケットではブースター帰還のための推進剤確保も必要だ。本記事ではロケットの上昇軌道設計の基礎から最適化手法までを解説する。


軌道力学の基礎

ツィオルコフスキーの公式

ロケット推進の基本方程式であるツィオルコフスキーの公式:

ΔV = Isp · g₀ · ln(m₀ / mf)

  • ΔV:速度増分
  • Isp:比推力
  • g₀:標準重力加速度
  • m₀:初期質量(推進剤含む)
  • mf:最終質量(推進剤消費後)

LEO投入に必要なΔVは約9.3〜9.5 km/s(軌道速度7.8 km/s + 重力損失 + 空力抵抗損失)だ。

損失の分類

理想的なΔVに加えて、実際の上昇軌道では以下の損失が発生する。

損失 典型的な値 原因
重力損失 1.0〜1.5 km/s 推力が重力に対抗するために消費
空力抵抗損失 0.05〜0.15 km/s 大気抵抗
ステアリング損失 0.05〜0.1 km/s 推力方向が速度方向と一致しない
背圧損失 ノズル設計に含む 海面での排気膨張不足

上昇軌道の最適化は、これらの損失を最小化して実効ΔVを最大化することに帰着する。


重力ターン(Gravity Turn)

基本概念

重力ターン(Gravity Turn)は、最も基本的かつ効率的な上昇軌道パターンだ。

  1. 垂直上昇フェーズ:打上げ直後、ロケットは垂直に上昇
  2. ピッチオーバーマヌーバ:数秒後にわずかなピッチ角を与える(キックアングル)
  3. 重力ターン:以後、推力方向を速度方向に一致させ、重力の成分でゆるやかに水平方向へ曲がる

推力方向と速度方向が一致するため、ステアリング損失がゼロ(理想的な場合)となる。

ピッチオーバーの設計

ピッチオーバーの角度とタイミングは上昇軌道の最適化において重要なパラメータだ。

パラメータ 典型的な値 影響
ピッチオーバー開始時間 T+5〜15秒 早すぎると動圧過大、遅すぎると重力損失増大
キックアングル 1〜5° 大きすぎるとMax-Q増大
ピッチオーバー継続時間 5〜15秒 緩やかな遷移が好ましい

飛行フェーズ

典型的な上昇プロファイル

2段式ロケットの上昇軌道は以下のフェーズで構成される。

Phase 1:第1段燃焼

サブフェーズ 時間(概算) 内容
垂直上昇 T+0〜10s 発射台クリアランス
ピッチオーバー T+10〜25s 重力ターンの開始
重力ターン T+25s〜Max-Q 速度方向に推力を整列
Max-Q通過 T+60〜80s 最大動圧点の通過
超音速飛行 Max-Q〜MECO 加速度増大(質量減少)
MECO T+150〜180s 第1段エンジンカットオフ

Phase 2:段分離・第2段燃焼

サブフェーズ 時間(概算) 内容
段分離 MECO+数秒 コースト後に分離
第2段点火 分離直後 上段エンジン始動
閉ループ誘導 第2段全期間 IGM等で軌道投入
SECO T+500〜600s 第2段エンジンカットオフ
ペイロード分離 SECO後 目標軌道到達

制約条件

Max-Q制約

最大動圧(Max-Q)はロケットの構造設計を支配する重要な制約だ。

q = ½ρv²

動圧は速度の2乗に比例し、大気密度ρは高度とともに指数関数的に減少する。速度増加と密度減少の競合により、高度10〜15kmで最大値に達する。

ロケット Max-Q(概算)
Falcon 9 約30 kPa
Saturn V 約33 kPa
SLS 約37 kPa

Max-Q通過時の構造荷重を制限するため、一部のロケットはこの区間でエンジンをスロットルダウンする。Falcon 9はMax-Q付近でMerlinエンジンをスロットルバックし、通過後にフルスロットルに復帰する。

加速度制限

推進剤の消費に伴い機体質量が減少するため、一定推力では加速度が増大する。有人ロケットでは乗員の耐G制限(一般に3〜4G)があり、エンジンのスロットリングが必要だ。

飛行安全制約

地上への破片落下リスクを制限するため、上昇軌道は飛行安全制約を満たす必要がある。

  • 排他的飛行範囲:飛行経路下の地上エリアの退避
  • 落下確率制約:エンジン停止時の破片落下範囲が許容範囲内
  • AFTSの判断境界:自律型飛行中断システムの動作範囲

最適制御理論の適用

最適上昇軌道問題

上昇軌道の最適化は、最適制御問題(Optimal Control Problem)として定式化される。

目的関数(最大化): 軌道投入時のペイロード質量(≒最終質量の最大化、≒ΔV損失の最小化)

制御変数: 推力方向(ピッチ角、ヨー角)、エンジンスロットル

制約条件: – 終端条件:目標軌道要素(高度、速度、軌道傾斜角等) – 経路制約:Max-Q、加速度制限、熱流束制限 – 制御制約:ジンバル角制限、スロットル範囲

間接法

ポントリャーギンの最大原理に基づく間接法は、必要条件(随伴方程式とハミルトニアンの最小化)を解く。

間接法の特徴: – 最適性条件が厳密 – 随伴変数の初期値推定が困難(BVP: 境界値問題) – 制約のアクティブ/非アクティブの切り替えが複雑

ロケット上昇軌道の間接法解析は理論的研究では重要だが、実用では直接法が主流だ。

直接法

直接法は制御変数を離散化し、非線形計画問題(NLP)に変換する。

直接コロケーション法 – 状態と制御を離散点で表現 – 力学方程式をコロケーション条件として拘束 – 大規模NLPをIPOPT等の内点法で解く

直接シューティング法 – 制御変数のみを離散化 – 初期値問題を前方積分 – 終端条件の残差をNLPで最小化


実用的な誘導手法

反復誘導モード(IGM)

IGM(Iterative Guidance Mode)は、Saturn V の上段で使用された閉ループ誘導であり、多くの現代ロケットの上段誘導の基礎となっている。

IGMのアルゴリズム: 1. 現在の航法状態から目標軌道条件までの必要ΔVを計算 2. ΔVを達成するための推力方向と燃焼時間を反復計算 3. 最適推力方向を線形時間関数で近似(Linear Tangent Steering) 4. 各制御サイクルで再計算(閉ループ)

IGMの利点: – 推進剤質量予測を使わず、リアルタイムの推力実測値に基づく – 打上げウィンドウの変動、推力のばらつきに自動適応 – 計算負荷が比較的低い

PEG(Powered Explicit Guidance)

PEG(Powered Explicit Guidance)はスペースシャトルで使用された上段誘導であり、IGMの発展形だ。

PEGの特徴: – 推力加速度の関数としてステアリング角を解析的に表現 – 複数のエンジンカットオフ条件に対応 – コースト弧(自由飛行区間)を含む軌道の最適化

凸最適化ベースの誘導

凸最適化を上昇軌道の誘導にも適用する研究が進んでいる。非線形力学を逐次線形化し、リアルタイムで最適軌道を再計算する。着陸誘導での実績を上昇誘導に展開する方向性だ。


再使用ロケットの軌道設計

RTLS(Return to Launch Site)の影響

Falcon 9のようなRTLS(発射地点帰還)型の再使用ロケットでは、ブースター帰還のための推進剤を確保する必要がある。

これは上昇軌道に以下の制約を追加する: – MECO速度の制限:帰還に必要なΔVを残すため、第1段の到達速度が制限される – MECO高度の調整:帰還軌道の設計に応じた分離高度 – 飛行方位角の制約:帰還飛行経路が陸地上を通過する制約

結果として、RTLS構成ではペイロード能力が15〜30%低下する。

ASDS(ドローン船着陸)

GTO等の高エネルギーミッションでは、第1段の速度が大きすぎてRTLSが不可能になる。この場合、ダウンレンジのASDS(Autonomous Spaceport Drone Ship)に着陸する。

ASDS構成では帰還のためのブーストバックバーンが不要であり、RTLS構成よりペイロード能力が高い。


風への対応

風バイアス軌道

高層風は上昇軌道に大きな影響を与える。特にMax-Q近傍での横風は構造荷重を増大させる。

風バイアス(Wind Biasing)は、事前に測定した風プロファイルに基づいて軌道を微調整する手法だ。打上げ前に上層気象観測(ラジオゾンデ)で取得した風データを誘導に反映する。

ロードリリーフとの連携

ロードリリーフ制御と組み合わせ、風荷重に対する構造応答を管理する。上昇軌道の風バイアスは迎角を最小化するマクロ的な対策であり、ロードリリーフは迎角の微小変動に対するリアルタイム制御だ。


設計最適化の実践

トレードスタディ

上昇軌道設計は多くのトレードオフを含む。

トレードオフ 一方 他方
重力損失 vs 空力損失 水平飛行→重力損失大 垂直飛行→空力損失小
ペイロード vs 帰還推進剤 帰還推進剤減→ペイロード増 帰還マージン低下
Max-Q vs 飛行時間 スロットルダウン→Max-Q低減 重力損失増大
飛行安全 vs 性能 安全制約厳しい→軌道制限 ペイロード能力低下

オフライン最適化 + オンライン誘導

実用的なアプローチは、オフラインで最適基準軌道を設計し、オンラインの閉ループ誘導で追従/再計算する2段階法だ。

オフラインの軌道最適化で公称条件のベスト軌道を設計し、オンライン誘導(IGM/PEG等)で実飛行中の偏差に対応する。


まとめ

上昇軌道設計は、ロケット工学の中でも最も古典的かつ重要な問題だ。重力ターンの基本概念から最適制御理論の適用まで、理論と実践が高度に統合された分野である。

再使用ロケットの登場により、上昇軌道設計には帰還制約という新たな次元が加わった。限られた推進剤をペイロードと帰還に最適配分する問題は、今日の軌道設計の最前線の課題だ。

オンライン凸最適化による軌道の実時間再計算が実現すれば、事前設計への依存度を低減し、より柔軟な上昇誘導が可能になるだろう。