JP5822676B2 - 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム - Google Patents
多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム Download PDFInfo
- Publication number
- JP5822676B2 JP5822676B2 JP2011256243A JP2011256243A JP5822676B2 JP 5822676 B2 JP5822676 B2 JP 5822676B2 JP 2011256243 A JP2011256243 A JP 2011256243A JP 2011256243 A JP2011256243 A JP 2011256243A JP 5822676 B2 JP5822676 B2 JP 5822676B2
- Authority
- JP
- Japan
- Prior art keywords
- current
- rocket
- predicted
- stage
- vector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Description
空中発射システムは、飛行体(飛行機など)に搭載したロケットを空中で発射するシステムであり、ロケットを地上から発射する地上発射システムに比べて発射地点の自由度が高い。
一方、空中発射システムでは、安全な地点からロケットを発射することができるため中間の目標地点を設ける必要が無く、ロケットを発射したときから最終的な目標地点に誘導することができる。
一方、固体ロケットの場合、固体燃料の装填に液体燃料ほどの時間はかからないため発射整備作業に要する時間を短縮することができる。但し、固体燃料の燃焼を途中で止めることはできない。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する。
前記多段式ロケット誘導装置は、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。
前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが次段の固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段のロケットエンジンが前記次段の固体燃料を燃焼し終わる予定時刻として予め記憶する次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出する。
前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する。
前記多段式ロケット誘導装置は、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる。
前記多段式ロケット誘導プログラムは、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させる。
前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが次段の固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段のロケットエンジンが前記次段の固体燃料を燃焼し終わる予定時刻として予め記憶する次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出する。
前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる。
前記多段式ロケット誘導プログラムは、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させる。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する。
前記多段式ロケット誘導方法は、
予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出し、
誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する。
前記多段式ロケット誘導方法は、
予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出し、
誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置と、
前記多段式ロケットの推力方向を制御するロケット制御装置と、
前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える。
前記多段式ロケット誘導装置は、
前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。
固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置と、
前記多段式ロケットの推力方向を制御するロケット制御装置と、
前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える。
前記多段式ロケット誘導装置は、
前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部とを備える。
多段式ロケットをできるだけ遠地点高度の大きな楕円軌道、もしくは、できるだけ近地点高度の大きな楕円軌道に誘導するシステムの形態について説明する。
実施の形態1における多段式ロケット誘導システム100の構成について、図1に基づいて説明する。
例えば、多段式ロケット誘導システム100は多段式ロケットに搭載される。但し、多段式ロケット誘導システム100の構成の少なくとも一部を多段式ロケットの外部(例えば、地上の施設)に設けても構わない。
例えば、航法装置110は慣性計測装置(IMU)とも呼ばれ、ジャイロと加速度計とを備える。航法装置110は、ジャイロによって計測される3軸方向の角速度と加速度計によって計測される3軸方向の加速度(加速度ベクトル)とを積分して3軸方向の速度(速度ベクトル)の変化量と3軸方向の位置(位置ベクトル)の変化量とを算出する。航法装置110は、速度ベクトルの変化量を前回の速度ベクトルに加算して現在の速度ベクトルを算出し、位置ベクトルの変化量を前回の位置ベクトルに加算して現在の位置ベクトルを算出する。
例えば、ロケット制御装置120は、多段式ロケットの機体やロケットエンジンの噴射向きを制御して多段式ロケットの推力方向を制御する。
推力加速度プロファイルは、ロケットエンジンの各時刻(経過時間)の推力加速度の予測値を示すデータである。例えば、推力加速度プロファイルは、各段のロケットエンジンが固体燃料を燃焼し始める予定時刻(推力加速度がゼロから上昇する時刻)や各段のロケットエンジンが固体燃料を燃焼し終わる予定時刻(推力加速度がゼロになる時刻)を示す。
予測演算部220は、現在の位置ベクトルと現在の速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる現段燃焼終了時刻と、現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出する。
予測演算部220は、算出した現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルと、に基づいて最終段燃焼終了時刻に多段式ロケットが飛行する予測の楕円軌道の軌道長半径(軌道要素予測値の一例)と近地点高度(または遠地点高度)(軌道要素予測値の一例)とを積分計算によって算出する。
予測演算部220は、現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、現段燃焼終了時刻と、次段のロケットエンジンが固体燃料を燃焼し始める次段燃焼開始時刻と、に基づいて次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出する。
予測演算部220は、算出した次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、次段燃焼終了時刻と、最終段燃焼終了時刻と、次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、最終段推力加速度プロファイルと、に基づいて予測の楕円軌道の軌道長半径と近地点高度(または遠地点高度)とを積分計算によって算出する。
予測演算部220は、新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在時刻と現段燃焼終了時刻とを用いて現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出し、最終段燃焼終了時刻に多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度(または遠地点高度)とを新たに算出する。
誘導演算部230は、算出した多段式ロケットの新たな推力方向の角度をロケット制御装置120に出力することによって、目標の近地点高度に位置する近地点を通って予測の楕円軌道より遠地点高度が大きい楕円軌道へ多段式ロケットを誘導する。
誘導演算部230は、算出した多段式ロケットの新たな推力方向の角度をロケット制御装置120に出力することによって、目標の遠地点高度に位置する遠地点を通って予測の楕円軌道より近地点高度が大きい楕円軌道へ多段式ロケットを誘導する。
推力加速度プロファイル、位置ベクトル、速度ベクトル、推力加速度ベクトル、軌道要素予測値、軌道要素目標値およびピッチ姿勢角は、誘導装置記憶部290に記憶するデータの一例である。
但し、多段式ロケット誘導システム100が誘導する多段式ロケットは2段式ロケットであっても4段以上の多段式ロケットであっても構わない。
実施の形態1における多段式ロケット誘導システム100による3段式ロケット誘導処理(多段式ロケット誘導方法、多段式ロケット誘導プログラムの一例)について、図2に基づいて説明する。
以下、第n段(n:1以上の整数)のロケットエンジンが固体燃料を燃焼しているときに行うロケットの誘導を「n段誘導」という。
また、前段のロケットエンジンが固体燃料を燃焼し終わって切り離されてから次段のロケットエンジンが固体燃料を燃焼し始めるまでのロケットの誘導、特に、2段誘導の終了時から3段誘導の開始時までのロケットの誘導を「コースト誘導」という。
図示省略する1段誘導(1段誘導の終了時から2段誘導の開始時までのコースト誘導を含む)において、多段式ロケット誘導システム100は、所定の誘導処理(例えば、従来の誘導処理)によってロケットのヨー姿勢角(方位角)を調整してロケットを目標の軌道面に誘導する。目標の軌道面は、目標の楕円軌道(遠地点高度または近地点高度が大きい楕円軌道)を含んだ平面である。
多段式ロケット誘導システム100は、図2に示す2段誘導以降の処理によってロケットのピッチ姿勢角を調整してロケットを目標の楕円軌道に誘導する。
つまり、多段式ロケット誘導システム100は、ロケットのピッチ姿勢角(推力方向)を周期的に正しい角度に調整する。
現在の誘導がいずれの誘導であるかを示すデータはロケットエンジンの切り離し時や固体燃料の燃焼開始時に誘導装置記憶部290に記憶されるものとする。
現在の誘導が2段誘導である場合、S120に進む。
現在の誘導がコースト誘導である場合、ロケットのピッチ姿勢角を変更せず、今回の周期の3段式ロケット誘導処理を終了する。
現在の誘導が3段誘導である場合、S140に進む。
S120の後、S150に進む。
S140の後、S150に進む。
オプションA「投入誤差の最小化」は、目標の遠地点と目標の近地点とを通る楕円軌道(目標軌道)へロケットを誘導するための誘導オプションである。
オプションB「遠地点高度の最大化」は、目標の近地点を通って遠地点高度ができるだけ大きい楕円軌道へロケットを誘導するための誘導オプションである。
オプションC「近地点高度の最大化」は、目標の遠地点を通って近地点高度ができるだけ大きい楕円軌道へロケットを誘導するための誘導オプションである。
誘導オプションの種類がオプションB(遠地点高度の最小化)である場合、S152に進む。
誘導オプションの種類がオプションC(近地点高度の最小化)である場合、S153に進む。
オプションAのΔyは、遠地点高度の目標値(軌道要素目標値の一例)と予測値との値差を用いた成分と、近地点高度の目標値(軌道要素目標値の一例)と予測値との値差を用いた成分と、を含んだベクトルである。オプションAのΔyの詳細については別途説明する。
S151の後、S154に進む。
オプションBのΔyは、軌道長半径の予測値の逆数を用いた成分と、近地点高度の目標値と予測値との値差を用いた成分と、を含んだベクトルである。オプションBのΔyの詳細については別途説明する。
S152の後、S154に進む。
オプションCのΔyは、軌道長半径の予測値の逆数を用いた成分と、遠地点高度の目標値と予測値との値差を用いた成分と、を含んだベクトルである。オプションCのΔyの詳細については別途説明する。
S153の後、S154に進む。
誘導演算部230は、誘導オプションの種類(S150)と現在の誘導(S110)との組み合わせに基づいて、図3に示すΔy(更新ベクトル)の計算式を計算してΔyを算出する(S151からS153)。
「Re」は赤道半径である。
「a」は軌道長半径である。
「rp」は近地点高度の予測値であり、「rpf」は近地点高度の目標値である。
「ra」は遠地点高度の予測値であり、「raf」は遠地点高度の目標値である。
「ω」は近地点引数の予測値であり、「ωf」は近地点引数の目標値である。
誘導オプションがオプションA「投入誤差の最小化」である場合(S151)、誘導演算部230は図3の項番1に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。
誘導オプションがオプションB「遠地点高度の最大化」である場合(S152)、誘導演算部230は図3の項番3に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。
誘導オプションがオプションC「近地点高度の最大化」である場合(S153)、誘導演算部230は図3の項番5に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は異なる。
誘導オプションがオプションA「投入誤差の最小化」である場合(S151)、誘導演算部230は図3の項番2に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。
誘導オプションがオプションB「遠地点高度の最大化」である場合(S152)、誘導演算部230は図3の項番4に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。
誘導オプションがオプションC「近地点高度の最大化」である場合(S153)、誘導演算部230は図3の項番6に示すΔyを算出する。2段誘導と3段誘導とでΔyの計算式は同じである。
誘導演算部230は、二乗和Δy2と誘導装置記憶部290に予め記憶する誘導閾値とを比較する。
二乗和Δy2が誘導閾値未満である場合、ロケットは目標の軌道に向けて正しく誘導されているため、今回の周期の3段式ロケット誘導処理を終了する。
二乗和Δy2が誘導閾値以上である場合、S160に進む。
独立変数とは、例えばピッチ姿勢角、後述するステアリング係数、最終段燃焼開始時刻である。誘導演算部230は、算出したピッチ姿勢角をロケット制御装置120に出力する。
関係式(1)において感度行列Jは、軌道要素予測値を独立変数で偏微分して得られるヤコビ行列である。
係数λは、繰り返し計算においてΔy2が前回値より小さい場合に小さく、Δy2が前回値より大きい場合に大きく設定する値である。
重み付け行列Wは、誘導オプションの種類に応じて予め定義する行列である。
誘導演算部230は、新たな独立変数のうちステアリング係数を式(2)に設定し、時間変数tに現在時刻から現段燃焼終了時刻までの時間および3段燃焼開始時刻から3段燃焼終了時刻までの時間を代入して式(2)を計算する。
式(3)は、ロケットの新たなピッチ姿勢角θのタンジェント値tanθを計算するための時間変数tの1次式である。式(3)を「リニアタンジェント則」という。
但し、誘導演算部230は、2段誘導で2次タンジェント則を用いて、3段誘導でリニアタンジェント則を用いてもよい。また、誘導演算部230は、各段の誘導時に2次タンジェント則とリニアタンジェント則とのいずれを選択しても構わない。
実施の形態1における2段誘導演算処理(S120)について、図4に基づいて説明する。
予測演算部220は、現在の位置ベクトルと、現在の速度ベクトルと、現在時刻と、2段燃焼終了時刻と、現在時刻から2段燃焼終了時刻までの各時刻の推力加速度と、を入力にして所定の予測演算関数を演算し、2段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。
S121の後、S122に進む。
予測演算部220は、2段燃焼終了時刻の位置ベクトルと、2段燃焼終了時刻の速度ベクトルと、現在時刻と、3段燃焼開始時刻と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、S121と同様に積分計算によって3段燃焼開始時刻の位置ベクトルと速度ベクトルとを算出する。
S122の後、S123に進む。
予測演算部220は、3段燃焼開始時刻の位置ベクトルと、3段燃焼開始時刻の速度ベクトルと、現在時刻と、3段燃焼終了時刻と、3段燃焼開始時刻から3段燃焼終了時刻までの各時刻の推力加速度と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、S121と同様に積分計算によって3段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。
軌道要素および軌道要素予測値の算出方法については別途説明する。
S123により、2段誘導演算処理(S120)は終了する。
実施の形態1における3段誘導演算処理(S140)について、図5に基づいて説明する。
予測演算部220は、現在の位置ベクトルと、現在の速度ベクトルと、現在時刻と、3段燃焼終了時刻と、現在時刻から3段燃焼終了時刻までの各時刻の推力加速度と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、図3のS121と同様に積分計算によって3段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。
軌道要素予測値の算出方法については別途説明する。
S141により、3段誘導演算処理(S140)は終了する。
実施の形態1における軌道要素について、図6、図7に基づいて説明する。
軌道要素として、近地点高度、遠地点高度、軌道長半径、軌道傾斜角、近地点引数、真近点離角、などが挙げられる。これらの軌道要素について以下に説明する。
近地点は目標軌道内で地球重心から最も近い点であり、地球重心から近地点までの距離を近地点半径という。「近地点高度rp」は近地点半径から地球の赤道半径を引いた距離である。
遠地点は目標軌道内で地球重心から最も遠い点であり、地球重心から遠地点までの距離を遠地点半径という。「遠地点高度ra」は遠地点半径から地球の赤道半径を引いた距離である。
「軌道傾斜角i」は目標軌道と地球の赤道面とが成す角度である。
「近地点引数ω」は昇交点から近地点までの角度(ロケットの進行方向側の角度)である。
図8に示すように、例えば、地球中心(または地球重心)を座標系の「原点o」とし、原点oからグリニッジ子午線と赤道との交点を通る方向を「X軸方向」とし、北極方向(北方向)を「Z軸方向」とし、Z軸とX軸とに直交し右手系を成す方向を「Y軸方向」とする。
図9に示すように、ロケットのピッチ角θは接線方向と角運動量方向とを含んだ平面に対して推力方向が成す角度である。ロケットのピッチ姿勢角θが表す方向は、ロケットの推力方向とみなすことができる。
接線方向は「(位置ベクトル×速度ベクトル)×位置ベクトル」(「×」は外積を表す)で定義される方向である。
軌道半径方向は位置ベクトルの向きに相当する方向である。
角運動量方向は接線方向と軌道半径方向とに対して垂直な方向であり、「位置ベクトル×速度ベクトル」(「×」は外積を表す)で定義される。
図11は、実施の形態1における3段誘導をシミュレーションした際のパラメータを示す表である。
図12は、実施の形態1におけるシミュレーションパターンを示す表である。
図13は、実施の形態1における誘導アルゴリズム(遠地点高度の最大化)のシミュレーション結果を示すグラフである。
図14は、実施の形態1における誘導アルゴリズム(近地点高度の最大化)のシミュレーション結果を示すグラフである。
図10、図11に示す各記号の意味は以下の通りである。
「A」は実施の形態1における誘導アルゴリズム(中間目標値を設定せずに全段を統合したロジック)を意味する。
「Q」は2次タンジェント則を意味する。
「L」はリニアタンジェント則を意味する。
「CX」はステアリング係数を意味する。
「θ」は3段燃焼中時刻のロケットのピッチ姿勢角を示す。
「T」は2段燃焼終了時刻から3段燃焼開始時刻までの時間を示す。
「rp」は近地点高度を示す。
「ra」は遠地点高度を示す。
「a」は軌道長半径を示す。
「ω」は近地点引数を示す。
500×250キロメートルの楕円軌道(遠地点高度:500km、近地点高度:250km)を当初の予定軌道として誘導オプションB「遠地点高度の最大化」を実行した結果、標準偏差「+1σ」と標準偏差「ノミナル」と標準偏差「−1σ」のRSS値はいずれも遠地点高度を200キロメートル以上大きくすることができることを示した(図13参照)。
同様に、誘導オプションC「近地点高度の最大化」を実行した結果、標準偏差「+1σ」と標準偏差「ノミナル」と標準偏差「−1σ」のRSS値はいずれも近地点高度を75キロメートル以上大きくすることができることを示した(図14参照)。
実施の形態1において、例えば、図15に示す誘導オプションに従ってロケットを目標軌道へ誘導する誘導アルゴリズムについて説明した。
図15において、「○」は軌道要素(近地点高度、遠地点高度、近地点引数)の目標値に対する予測値の誤差を最小にすることを意味し、「max」は軌道要素(近地点高度、遠地点高度)をできるだけ大きくすることを意味する。
図16において、ロケット誘導装置200は、CPU911(Central Processing Unit)を備えている。CPU911は、バス912を介してROM913、RAM914、通信ボード915と接続され、これらのハードウェアデバイスを制御する。通信ボード915は、有線または無線でネットワークに接続している。
フローチャートなどに基づいて説明する実施の形態の処理はCPU911やその他のハードウェアを用いて実行される。
実施の形態1により、固体燃料を用いる多段式ロケットを中間の目標地点を経由せずに高い精度で目標軌道に誘導することができる。また、固体燃料を用いることにより発射整備作業の時間を短縮し、中間の目標地点を設定しないことによりミッション解析期間を短縮することができる。
つまり、実施の形態1により、緊急の打ち上げ要請にもタイムリーに応えられるという即応性を有した上で、多段式ロケットを目標軌道に高い精度で誘導することができる。
例えば、ペイロード(人工衛星、探査機など)を搭載したロケットをペイロードの目標軌道またはそれに近い軌道へ投入することができる。これにより、ペイロードの燃料を節約することができ、余った燃料をペイロードの軌道修正などに利用することもできる。
実施の形態1では2段誘導以降でロケットのピッチ姿勢角(推力方向)を調整するアルゴリズムを例に挙げて説明したが、1段誘導においてもロケットのピッチ姿勢角を調整しても構わない。調整方法は2段誘導と同様である。
Claims (12)
- 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置において、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部と
を備えることを特徴とする多段式ロケット誘導装置。 - 前記予測演算部は、
前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが次段の固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段のロケットエンジンが前記次段の固体燃料を燃焼し終わる予定時刻として予め記憶する次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出する
ことを特徴とする請求項1記載の多段式ロケット誘導装置。 - 前記予測演算部は、
前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する
ことを特徴とする請求項1または請求項2記載の多段式ロケット誘導装置。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置において、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部と
を備えることを特徴とする多段式ロケット誘導装置。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる多段式ロケット誘導プログラムであって、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させる
ことを特徴とする多段式ロケット誘導プログラム。 - 前記予測演算部は、
前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが次段の固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段のロケットエンジンが前記次段の固体燃料を燃焼し終わる予定時刻として予め記憶する次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記予測の楕円軌道の前記軌道長半径と前記近地点高度とを積分計算によって算出する
ことを特徴とする請求項5記載の多段式ロケット誘導プログラム。 - 前記予測演算部は、
前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する
ことを特徴とする請求項5または請求項6記載の多段式ロケット誘導プログラム。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる多段式ロケット誘導プログラムであって、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部として前記多段式ロケット誘導装置を機能させる
ことを特徴とする多段式ロケット誘導プログラム。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する多段式ロケット誘導方法において、
予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出し、
誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する
ことを特徴とする多段式ロケット誘導方法。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する多段式ロケット誘導方法おいて、
予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出し、
誘導演算部が、前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する
ことを特徴とする多段式ロケット誘導方法。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルとを計測する航法装置と、
前記多段式ロケットの推力方向を制御するロケット制御装置と、
前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える多段式ロケット誘導システムにおいて、
前記多段式ロケット誘導装置は、
前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と近地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の近地点高度と目標の近地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の近地点高度に位置する近地点を通って前記予測の楕円軌道より遠地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部と
を備えることを特徴とする多段式ロケット誘導システム。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルとを計測する航法装置と、
前記多段式ロケットの推力方向を制御するロケット制御装置と、
前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える多段式ロケット誘導システムにおいて、
前記多段式ロケット誘導装置は、
前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度を示す最終段推力加速度プロファイルとに基づいて、前記最終段燃焼終了時刻に前記多段式ロケットが飛行する予測の楕円軌道の軌道長半径と遠地点高度とを積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記予測の楕円軌道の軌道長半径と、前記予測演算部によって算出された前記予測の楕円軌道の遠地点高度と目標の遠地点高度との差とに基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力することによって前記目標の遠地点高度に位置する遠地点を通って前記予測の楕円軌道より近地点高度が大きい楕円軌道へ前記多段式ロケットを誘導する誘導演算部と
を備えることを特徴とする多段式ロケット誘導システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011256243A JP5822676B2 (ja) | 2011-11-24 | 2011-11-24 | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011256243A JP5822676B2 (ja) | 2011-11-24 | 2011-11-24 | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2013107585A JP2013107585A (ja) | 2013-06-06 |
JP2013107585A5 JP2013107585A5 (ja) | 2015-01-08 |
JP5822676B2 true JP5822676B2 (ja) | 2015-11-24 |
Family
ID=48704809
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011256243A Active JP5822676B2 (ja) | 2011-11-24 | 2011-11-24 | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5822676B2 (ja) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6256196B2 (ja) * | 2014-05-28 | 2018-01-10 | 株式会社Ihi | ロケットの制御装置 |
CN110031003B (zh) * | 2019-02-21 | 2023-03-24 | 中国西安卫星测控中心 | 一种火箭上面级最优可达轨道快速规划计算方法 |
CN112179217B (zh) * | 2020-10-27 | 2022-05-31 | 中国运载火箭技术研究院 | 固体运载火箭的制导方法以及装置、存储介质、电子装置 |
CN113761670B (zh) * | 2020-10-27 | 2024-02-13 | 深圳三零三防务科技有限公司 | 一种运载火箭的平衡飞行理论及在线入轨能力评估方法 |
CN112462794B (zh) * | 2020-11-09 | 2024-03-26 | 航天科工火箭技术有限公司 | 一种演示验证火箭悬停制导方法及系统 |
CN112270108B (zh) * | 2020-11-13 | 2023-07-18 | 北京航天自动控制研究所 | 一种运载火箭椭圆轨道在线规划方法 |
CN113189870A (zh) * | 2021-04-12 | 2021-07-30 | 大连理工大学 | 一种火箭推力下降故障下椭圆救援轨道的轨迹重规划方法 |
CN113467497B (zh) * | 2021-07-08 | 2023-09-19 | 北京星途探索科技有限公司 | 一种满足载荷投放点多约束条件下的能量管理制导方法 |
CN113602532B (zh) * | 2021-08-25 | 2023-05-16 | 重庆零壹空间科技集团有限公司 | 一种固体运载火箭入轨修正方法 |
CN114216376B (zh) * | 2021-12-09 | 2023-11-14 | 北京航天自动控制研究所 | 运载火箭的多载荷分级优化方法 |
CN114413691A (zh) * | 2021-12-24 | 2022-04-29 | 北京航天自动控制研究所 | 一种运载火箭推力下降故障的跨滑行段解析制导重构方法 |
CN114858011B (zh) * | 2022-07-05 | 2022-10-14 | 星河动力(北京)空间科技有限公司 | 运载火箭 |
-
2011
- 2011-11-24 JP JP2011256243A patent/JP5822676B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2013107585A (ja) | 2013-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5822676B2 (ja) | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム | |
JP5822675B2 (ja) | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム | |
CN110471450B (zh) | 在高度速度剖面内直接规划再入轨迹的方法 | |
CN109398762B (zh) | 一种固体火箭基于椭圆转移轨道的入轨式弹道设计方法 | |
Xie et al. | Highly constrained entry trajectory generation | |
CN110908396B (zh) | 可重复使用运载器全阶段再入返回制导方法 | |
US9702674B2 (en) | Method and apparatus for GPS-denied navigation of spin-stabilized projectiles | |
US9711052B2 (en) | Flight path search device and flight path search program | |
CN109062241B (zh) | 基于线性伪谱模型预测控制的自主全射向再入制导方法 | |
Li et al. | Improved artificial potential field based lateral entry guidance for waypoints passage and no-fly zones avoidance | |
CN104007665A (zh) | 一种固液动力飞行器飞行仿真测试系统 | |
CN111368369A (zh) | 运载火箭的飞行性能参数设计方法、装置和电子设备 | |
Wolf et al. | Performance trades for Mars pinpoint landing | |
JP5683224B2 (ja) | ロケット誘導計算装置、ロケット誘導システム、ロケット、ロケット誘導計算プログラムおよびロケット誘導計算装置のロケット誘導計算方法 | |
CN114280934B (zh) | 一种可重复使用运载火箭全程轨迹规划方法 | |
CN106371312A (zh) | 基于模糊控制器的升力式再入预测‑校正制导方法 | |
Von der Porten et al. | Powered explicit guidance modifications and enhancements for space launch system Block-1 and Block-1B vehicles | |
US11505324B2 (en) | Methods and systems for controlling supersonic flight entry/exit of a vehicle | |
CN108298110A (zh) | 一种两级入轨空天飞行器上升段轨迹及设计方法 | |
CN114491811A (zh) | 一种运载火箭分离体的弹道设计方法 | |
Gerth et al. | Guidance for autonomous precision landing on atmosphereless bodies | |
CN104020678A (zh) | 一种基于月表地形的动力下降初始点参数优化方法 | |
Fuhry | Adaptive atmospheric reentry guidance for the Kistler K-1 orbital vehicle | |
De Vita et al. | Assessment of hypersonic flights operation scenarios: Analysis of launch and reentry trajectories, and derived top level vehicle system and support infrastructure concepts and requirements | |
Belokonov et al. | Passive stabilization systems for CubeSat nanosatellites: general principles and features |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20141118 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20141118 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150805 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150908 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20151006 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5822676 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |