JP5822675B2 - 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム - Google Patents
多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム Download PDFInfo
- Publication number
- JP5822675B2 JP5822675B2 JP2011256242A JP2011256242A JP5822675B2 JP 5822675 B2 JP5822675 B2 JP 5822675B2 JP 2011256242 A JP2011256242 A JP 2011256242A JP 2011256242 A JP2011256242 A JP 2011256242A JP 5822675 B2 JP5822675 B2 JP 5822675B2
- Authority
- JP
- Japan
- Prior art keywords
- current
- rocket
- stage
- end time
- 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は、算出した次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、次段燃焼終了時刻と、最終段燃焼終了時刻と、次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、最終段推力加速度プロファイルと、に基づいて最終段燃焼終了時刻の軌道の予測値を積分計算によって算出する。
誘導演算部230は、最終段燃焼終了時刻の軌道の予測値と所定の軌道の目標値との差に基づいて多段式ロケットの新たな推力方向の角度(ピッチ姿勢角)を算出する。
誘導演算部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に進む。
現在の誘導がコースト誘導である場合、S130に進む。
現在の誘導が3段誘導である場合、S140に進む。
S120の後、S150に進む。
S130の後、S150に進む。
S140の後、S150に進む。
S150の後、S151に進む。
軌道予測値と軌道目標値との値差が誘導閾値未満である場合、ロケットは目標の軌道に向けて正しく誘導されているため、今回の周期の3段式ロケット誘導処理を終了する。
軌道予測値と軌道目標値との値差が誘導閾値以上である場合、S160に進む。
独立変数とは、例えばピッチ姿勢角、後述するステアリング係数、最終段燃焼開始時刻である。
誘導演算部230のうちS120およびS140は、新たな独立変数のうちステアリング係数を式(2)に設定し、時間変数tに現在時刻から現段燃焼終了時刻までの時間および3段燃焼開始時刻から3段燃焼終了時刻までの時間を代入して式(2)を計算する。
式(3)は、ロケットの新たなピッチ姿勢角θのタンジェント値tanθを計算するための時間変数tの1次式である。式(3)を「リニアタンジェント則」という。
但し、誘導演算部230は、2段誘導時および3段誘導時に2次タンジェント則とリニアタンジェント則とのどちらを用いてもよい。
実施の形態1における2段誘導演算処理(S120)について、図3に基づいて説明する。
予測演算部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におけるコースト誘導演算処理(S130)について、図4に基づいて説明する。
予測演算部220は、現在の位置ベクトルと、現在の速度ベクトルと、現在時刻と、3段燃焼開始時刻と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、図3のS121と同様に積分計算によって3段燃焼開始時刻の位置ベクトルと速度ベクトルとを算出する。
S131の後、S132に進む。
予測演算部220は、3段燃焼開始時刻の位置ベクトルと、3段燃焼開始時刻の速度ベクトルと、現在時刻と、3段燃焼終了時刻と、3段燃焼開始時刻から3段燃焼終了時刻までの各時刻の推力加速度と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、図3のS121と同様に積分計算によって3段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。
軌道要素予測値の算出方法については別途説明する。
S132により、コースト誘導演算処理(S130)は終了する。
実施の形態1における3段誘導演算処理(S140)について、図5に基づいて説明する。
予測演算部220は、現在の位置ベクトルと、現在の速度ベクトルと、現在時刻と、3段燃焼終了時刻と、現在時刻から3段燃焼終了時刻までの各時刻の推力加速度と、を入力にして予測演算関数を演算する。予測演算関数の演算において、予測演算部220は、図3のS121と同様に積分計算によって3段燃焼終了時刻の位置ベクトルと速度ベクトルとを算出する。
軌道要素予測値の算出方法については別途説明する。
S141により、3段誘導演算処理(S140)は終了する。
実施の形態1における軌道要素について、図6、図7に基づいて説明する。
軌道要素として、近地点高度、遠地点高度、軌道長半径、軌道傾斜角、近地点引数、真近点離角、などが挙げられる。これらの軌道要素について以下に説明する。
ロケット誘導装置200は例えば1つまたは複数の軌道要素を用いる。
近地点は目標軌道内で地球重心から最も近い点であり、地球重心から近地点までの距離を近地点半径という。「近地点高度rp」は近地点半径から地球の赤道半径を引いた距離である。
遠地点は目標軌道内で地球重心から最も遠い点であり、地球重心から遠地点までの距離を遠地点半径という。「遠地点高度ra」は遠地点半径から地球の赤道半径を引いた距離である。
「軌道傾斜角i」は目標軌道と地球の赤道面とが成す角度である。
「近地点引数ω」は昇交点から近地点までの角度(ロケットの進行方向側の角度)である。
図8に示すように、例えば、地球中心(または地球重心)を座標系の「原点o」とし、原点oからグリニッジ子午線と赤道との交点を通る方向を「X軸方向」とし、北極方向(北方向)を「Z軸方向」とし、Z軸とX軸とに直交し右手系を成す方向を「Y軸方向」とする。
図9に示すように、ロケットのピッチ角θは接線方向と角運動量方向とを含んだ平面に対して推力方向が成す角度である。ロケットのピッチ姿勢角θが表す方向は、ロケットの推力方向とみなすことができる。
接線方向は「(位置ベクトル×速度ベクトル)×位置ベクトル」(「×」は外積を表す)で定義される方向である。
軌道半径方向は位置ベクトルの向きに相当する方向である。
角運動量方向は接線方向と軌道半径方向とに対して垂直な方向であり、「位置ベクトル×速度ベクトル」(「×」は外積を表す)で定義される。
図11は、実施の形態1におけるコースト誘導をシミュレーションした際のパラメータを示す表である。
図12は、実施の形態1における3段誘導をシミュレーションした際のパラメータを示す表である。
図13は、実施の形態1におけるシミュレーションパターンを示す表である。
図14は、実施の形態1における誘導アルゴリズムのシミュレーション結果(円軌道)を示すグラフである。
図15は、実施の形態1における誘導アルゴリズムのシミュレーション結果(楕円軌道)を示すグラフである。
図10から図15に示す各記号の意味は以下の通りである。
「A」は実施の形態1における誘導アルゴリズム(中間目標値を設定せずに全段を統合したロジック)を意味する。
「B」は従来の誘導アルゴリズム(各段毎に中間目標値を設定するロジック)を意味する。
「Q」は2次タンジェント則を意味する。
「L」はリニアタンジェント則を意味する。
「CX」はステアリング係数を意味する。
「θ」は3段燃焼中のロケットのピッチ姿勢角を示す。
「T」は図10においては2段燃焼終了時刻から3段燃焼開始時刻までの時間を示し、図11においては現在時刻から3段燃焼開始時刻までの時間を示す。
「rp」は近地点高度を示す。
「ra」は遠地点高度を示す。
「f」は真近点離角を示す。
図11において、近地点高度rpと遠地点高度raとは3段燃焼終了時刻の値であり、真近点離角fは3段燃焼開始時刻の値である。
図12において、従属変数(rp,ra)は3段燃焼終了時刻の値である。
図14は、目標の円軌道(500×500キロメートル)に対する近地点高度、遠地点高度および軌道長半径の誤差(単位:キロメートル)を示している。
図15は、目標の円軌道(500×250キロメートル)に対する近地点高度、遠地点高度および軌道長半径の誤差(単位:キロメートル)を示している。
実施の形態1における誘導アルゴリズム「A」とリニアタンジェント則「リニア」との組み合わせで誤差が大きいケースがあるが、これは、この組み合せでは、誘導アルゴリズム「A」と2次タンジェント則「2次」との組み合わせに比べてロバスト性が落ちるため、精度が悪くなったケースがあったためである。
図16において、ロケット誘導装置200は、CPU911(Central Processing Unit)を備えている。CPU911は、バス912を介してROM913、RAM914、通信ボード915と接続され、これらのハードウェアデバイスを制御する。通信ボード915は、有線または無線でネットワークに接続している。
フローチャートなどに基づいて説明する実施の形態の処理はCPU911やその他のハードウェアを用いて実行される。
実施の形態1により、固体燃料を用いる多段式ロケットを中間の目標地点を経由せずに高い精度で目標軌道に誘導することができる。また、固体燃料を用いることにより発射整備作業の時間を短縮し、中間の目標地点を設定しないことによりミッション解析期間を短縮することができる。
つまり、実施の形態1により、緊急の打ち上げ要請にもタイムリーに応えられるという即応性を有した上で、多段式ロケットを目標軌道に高い精度で誘導することができる。
実施の形態1では2段誘導以降でロケットのピッチ姿勢角(推力方向)を調整するアルゴリズムを例に挙げて説明したが、1段誘導時や1段誘導の終了時から2段誘導の開始時までのコースト誘導時においてもロケットのピッチ姿勢角を調整しても構わない。調整方法は2段誘導時や2段誘導の終了時から3段誘導の開始時までのコースト誘導時と同様である。
Claims (10)
- 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置において、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルと、に基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルと、に基づいて前記最終段燃焼終了時刻の軌道の予測値を積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記最終段燃焼終了時刻の軌道の予測値と所定の軌道の目標値とを比較し、前記最終段燃焼終了時刻の軌道の予測値と前記所定の軌道の目標値との差に基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力する誘導演算部と
を備えることを特徴とする多段式ロケット誘導装置。 - 前記予測演算部は、
前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが次段の固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段のロケットエンジンが前記次段の固体燃料を燃焼し終わる予定時刻として予め記憶する次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記最終段燃焼終了時刻の軌道の予測値を積分計算によって算出する
ことを特徴とする請求項1記載の多段式ロケット誘導装置。 - 前記予測演算部は、
前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する
ことを特徴とする請求項1または請求項2記載の多段式ロケット誘導装置。 - 前記予測演算部は、前記最終段燃焼終了時刻の軌道の予測値を算出し、
前記誘導演算部は、予測値と目標値との差に基づいて、前記多段式ロケットの新たな推力方向の角度を算出する
ことを特徴とする請求項1から請求項3いずれかに記載の多段式ロケット誘導装置。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置を機能させる多段式ロケット誘導プログラムであって、
前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて前記最終段燃焼終了時刻の軌道の予測値を積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記最終段燃焼終了時刻の軌道の予測値と所定の軌道の目標値とを比較し、前記最終段燃焼終了時刻の軌道の予測値と前記所定の軌道の目標値との差に基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力する誘導演算部として前記多段式ロケット誘導装置を機能させる
ことを特徴とする多段式ロケット誘導プログラム。 - 前記予測演算部は、
前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、次段のロケットエンジンが次段の固体燃料を燃焼し始める予定時刻として予め記憶する次段燃焼開始時刻と、に基づいて前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記次段燃焼開始時刻の予測位置ベクトルと予測速度ベクトルと、前記次段のロケットエンジンが前記次段の固体燃料を燃焼し終わる予定時刻として予め記憶する次段燃焼終了時刻と、前記最終段燃焼終了時刻と、前記次段燃焼終了時刻までの推力加速度の予測値を示す次段推力加速度プロファイルと、前記最終段推力加速度プロファイルと、に基づいて前記最終段燃焼終了時刻の軌道の予測値を積分計算によって算出する
ことを特徴とする請求項5記載の多段式ロケット誘導プログラム。 - 前記予測演算部は、
前記ロケット制御装置が前記多段式ロケットの推力方向の角度を前記新たな推力方向の角度に変更した後、前記航法装置から新たな現在の位置ベクトルと新たな現在の速度ベクトルと新たな現在の推力加速度ベクトルとを入力し、前記新たな現在の位置ベクトルと前記新たな現在の速度ベクトルと新たな現在時刻と前記現段燃焼終了時刻と前記現段推力加速度プロファイルとを用いて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを新たに算出する
ことを特徴とする請求項5または請求項6記載の多段式ロケット誘導プログラム。 - 前記予測演算部は、前記最終段燃焼終了時刻の軌道の予測値を算出し、
前記誘導演算部は、予測値との目標値との差に基づいて、前記多段式ロケットの新たな推力方向の角度を算出する
ことを特徴とする請求項5から請求項7いずれかに記載の多段式ロケット誘導プログラム。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットを誘導する多段式ロケット誘導装置によって実行する多段式ロケット誘導方法において、
予測演算部が、前記多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて前記最終段燃焼終了時刻の軌道の予測値を積分計算によって算出し、
誘導演算部が、前記予測演算部によって算出された前記最終段燃焼終了時刻の軌道の予測値と所定の軌道の目標値とを比較し、前記最終段燃焼終了時刻の軌道の予測値と前記所定の軌道の目標値との差に基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記多段式ロケットの推力方向を制御するロケット制御装置に出力する
ことを特徴とする多段式ロケット誘導方法。 - 固体燃料を燃焼して推力を得るロケットエンジンを複数段搭載した多段式ロケットの位置ベクトルと速度ベクトルと推力加速度ベクトルとを計測する航法装置と、
前記多段式ロケットの推力方向を制御するロケット制御装置と、
前記多段式ロケットを誘導する多段式ロケット誘導装置と、を備える多段式ロケット誘導システムにおいて、
前記多段式ロケット誘導装置は、
前記航法装置から現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルとを入力し、入力した現在の位置ベクトルと現在の速度ベクトルと現在の推力加速度ベクトルと、現在時刻と、推力を得ている現段のロケットエンジンが現段の固体燃料を燃焼し終わる予定時刻として予め記憶する現段燃焼終了時刻と、前記現段燃焼終了時刻までの推力加速度の予測値を示す現段推力加速度プロファイルとに基づいて前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルとを積分計算によって算出し、算出した前記現段燃焼終了時刻の予測位置ベクトルと予測速度ベクトルと、前記現段燃焼終了時刻と、最終段のロケットエンジンが最終段の固体燃料を燃焼し終わる予定時刻として予め記憶する最終段燃焼終了時刻と、前記最終段燃焼終了時刻までの推力加速度の予測値を示す最終段推力加速度プロファイルとに基づいて前記最終段燃焼終了時刻の軌道の予測値を積分計算によって算出する予測演算部と、
前記予測演算部によって算出された前記最終段燃焼終了時刻の軌道の予測値と所定の軌道の目標値とを比較し、前記最終段燃焼終了時刻の軌道の予測値と前記所定の軌道の目標値との差に基づいて前記多段式ロケットの新たな推力方向の角度を算出し、算出した前記多段式ロケットの新たな推力方向の角度を前記ロケット制御装置に出力する誘導演算部とを備えることを特徴とする多段式ロケット誘導システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011256242A JP5822675B2 (ja) | 2011-11-24 | 2011-11-24 | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011256242A JP5822675B2 (ja) | 2011-11-24 | 2011-11-24 | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2013107584A JP2013107584A (ja) | 2013-06-06 |
JP2013107584A5 JP2013107584A5 (ja) | 2015-01-08 |
JP5822675B2 true JP5822675B2 (ja) | 2015-11-24 |
Family
ID=48704808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011256242A Active JP5822675B2 (ja) | 2011-11-24 | 2011-11-24 | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5822675B2 (ja) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112179217B (zh) * | 2020-10-27 | 2022-05-31 | 中国运载火箭技术研究院 | 固体运载火箭的制导方法以及装置、存储介质、电子装置 |
CN112231632B (zh) * | 2020-12-08 | 2021-03-19 | 北京星际荣耀空间科技股份有限公司 | 一种运载火箭风估计方法、装置、设备及存储介质 |
CN113834386B (zh) * | 2021-10-29 | 2023-02-28 | 湖北航天技术研究院总体设计所 | 一种固体运载火箭大气层外制导控制方法 |
CN114199085A (zh) * | 2021-11-12 | 2022-03-18 | 大连船舶重工集团有限公司 | 一种补偿落点偏差的海上火箭回收方法 |
CN114739240B (zh) * | 2022-03-30 | 2023-04-07 | 蓝箭航天空间科技股份有限公司 | 运载火箭级间冷分离时序优化方法及装置 |
CN114964264B (zh) * | 2022-07-11 | 2022-10-21 | 东方空间技术(山东)有限公司 | 一种空间目标定位方法、装置、计算机设备及存储介质 |
CN115600050B (zh) * | 2022-12-13 | 2023-04-07 | 东方空间(西安)宇航技术有限公司 | 一种火箭芯级间推冲气缸的充气量确定方法、装置及设备 |
-
2011
- 2011-11-24 JP JP2011256242A patent/JP5822675B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2013107584A (ja) | 2013-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5822676B2 (ja) | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム | |
JP5822675B2 (ja) | 多段式ロケット誘導装置、多段式ロケット誘導プログラム、多段式ロケット誘導方法および多段式ロケット誘導システム | |
US9702674B2 (en) | Method and apparatus for GPS-denied navigation of spin-stabilized projectiles | |
Desai et al. | Entry, descent, and landing performance of the mars phoenix lander | |
Xie et al. | Highly constrained entry trajectory generation | |
Li et al. | Improved artificial potential field based lateral entry guidance for waypoints passage and no-fly zones avoidance | |
CN111368369B (zh) | 运载火箭的飞行性能参数设计方法、装置和电子设备 | |
JP5683224B2 (ja) | ロケット誘導計算装置、ロケット誘導システム、ロケット、ロケット誘導計算プログラムおよびロケット誘導計算装置のロケット誘導計算方法 | |
Dukeman | Atmospheric ascent guidance for rocket-powered launch vehicles | |
Barton et al. | Autolanding trajectory design for the X-34 | |
Slegers et al. | Terminal guidance of autonomous parafoils in high wind-to-airspeed ratios | |
Wolf et al. | Performance trades for Mars pinpoint landing | |
CN105005651A (zh) | 航天器脉冲交会轨迹的梯度分割区间优化设计方法 | |
CN114280934B (zh) | 一种可重复使用运载火箭全程轨迹规划方法 | |
CN109542111B (zh) | 基于分段式的无人机超低空飞行控制方法 | |
CN114491811A (zh) | 一种运载火箭分离体的弹道设计方法 | |
CN104252548A (zh) | 一种燃料最优的火星探测器入轨目标点设计方法 | |
CN108562293B (zh) | 基于干扰观测器的行星着陆有限时间控制方法 | |
Gerth et al. | Guidance for autonomous precision landing on atmosphereless bodies | |
Fuhry | Adaptive atmospheric reentry guidance for the Kistler K-1 orbital vehicle | |
CN108454884A (zh) | 一种动力上升安全制导方法及系统 | |
Hanson et al. | Learning about ares I from Monte Carlo simulation | |
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 | |
Prabhakar et al. | A predictive explicit guidance scheme for ballistic missiles | |
Striepe et al. | An atmospheric guidance algorithm testbed for the Mars Surveyor Program 2001 Orbiter and Lander |
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: 5822675 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 |