北朝鮮の弾道ミサイルに備えて

弾道の特性

戻る: Top Page へ  「はじめに」 へ  Nodong Missile の概要 へ 

進む: 標準的な弾道 へ  パラメータ分析 へ  パラメータ分析(Isp) へ  パラメータ分析(Payload) へ  

我が国に飛来する Nodong の弾道 へ  首都圏に飛来する弾道(続き) へ  三陸沖に飛来する弾道

沖縄諸島に飛来する弾道 へ  Taepodong_Missile  へ  Taepodong-2 Missile 

Taepodong-2 Missile 検証 へ  Taepodong-3 の弾道特性を推定する へ  「おわりに」 へ

Nodong Missile

2  Nodong Missile 弾道

2.1 弾道計算のためのパラメータ設定

 弾道計算(computer simulation)のために、基本的な計算パラメータを下表のように設定する。

表1 弾道計算のためのパラメータ
パラメータ 記号 備      考
 弾頭質量 mw  700 -1,200 [kg]  基準値 1,000 [kg]
 ブースター質量 推進剤質量 mp   12,960 [kg]  推進剤質量比 = 0.85
空虚質量 md    2,290 [kg]  構造質量比    = 0.15
 比推力 Isp       296 [s]  UDMH/IRFNA
 燃焼時間 tb       110 [s]  
 抗力係数 機体全部 Cd       0.8 
弾頭部 Cd     0.15 - 0.20
 基準断面積 A     1.37 [m2]  ミサイル径から計算
 発射場位置 緯度 φ  北朝鮮の主要射場
 発射方向 λ  東~南
 経路角
 (上下角)
最終値 θ  ブースト終了時の経路角
角速度 θdot
 空気密度 ρ  近似式による  高度に依存

<弾頭>

 弾頭質量については、例えば、「missile.index」(資料1)によれば、700 - 1,200 [kg] としている。ここでは、弾頭質量をパラメータとして、それに対応する射程の変化をシミュレートする。他のパラメーターを変化させる場合は、弾頭質量の基準値を 1,000 [kg] とする。

<発射質量>(Launch Weight; L.W)

 発射質量については、16,250 [kg] としている資料が多いが、16,500 [kg] としている資料(例えば、「MissileThreat.com」)あるいは、弾頭とブースター(ロケットモーター)の質量のみを示して、あえて、発射質量を示していない資料もある。例えば、「An Analysis of the North Korean Nodong Missile」(資料2)が、その例である。

 「missile.index」では、「発射重量」を、「発射時に満載したミサイルの重量。kgで示す。」(注:「重量」と「質量」の混同がある)としてしている。ここで、「発射質量」=「弾頭質量」+「ブースター質量」として、この定義に従えば、発射質量=16,250 [kg]弾頭質量 = 1,200 [kg] の場合、ブースター質量 = 15,050 [kg] になる。また、発射質量=16,500 [kg] の場合、ブースター質量 = 15,300 [kg] になる。

<ブースター質量>

 ブースター質量 m は、推進剤質量 mp 空虚(dry)質量 md の和 m = mp + md である。そして、λ = mp / m を「推進剤質量比」、ε = md / m を「構造質量比」と呼ぶ。

 ブースター質量を直接示している資料は見当たらない。資料2では、推進剤質量を「Scud B 改」の質量(4,000 [kg])の4倍として、16,000 [kg] と見積もっている。また、空虚質量については、ブースター本体の外壁、フィン及び支持具等がアルミ合金で作られていると仮定して、推進剤タンク、エンジン(4個)+ポンプ等及び誘導制御装置を含めて、2,800 - 3,000 [kg] と見積もっている。したがって、この資料によれば、ブースター質量は、18,800 - 19,000 [kg] になる。

 「北朝鮮のミサ汁脅威と我が国の対応」(金ギョンス韓国国防研究院責任研究委員、国防ジャーナル1998.10月号)(資料3)では、推進剤(燃料:TM-185, 酸化剤:AK-271)質量12,912 [kg] と見積もっている。

 「Federation of American Scientists」(North Korean Missile Nodong)では、空虚質量1,780 - 2,180 [kg] と見積もっている。

 ここでは、発射質量16,250 [kg]、弾頭質量1,000 [kg] に抑えて、これから、ブースター質量を、16,250 - 1,000 = 15,250 [kg] と見積もって、この質量を推進剤質量比0.85)及び構造質量比0.15)をもって、それぞれ、推進剤質量12,963[kg])及び空虚質量2,287[kg])に配分することにする。なお、発射質量 = 16,500 [kg] とすると、これらの値は、それぞれ、13,175 [kg] 及び  2,325 [kg] となる。

<比推力>

 ここでは、推進剤として IRFNA/UDMH (混合比 3.12)を使用するものとして、比推力の値として、「Sutton, George P. : "Rocket Propulsion Elements, An Introduction to the Engineering of Rocket, 6th Edition", John Wiley & Sons, 1992」に載せられている値: 329 [s]真空中)を基準に、ブースト期間が海面高度(Sea level)から比較的低高度の範囲であることを考慮して、真空中の値の 90 % 、すなわち、296 [s] (329 x 0.9)とする。ただし、資料によっては、Sucad の推進剤(燃料TM-185:Gasoline 20%, Kerosene 80%, 酸化剤AK-271:N2O4 27%. HNO3 73%, Isp 264[s](真空中))をそのまま使っている場合もあるので、ここでは、250、 270、 290 [s] の場合について示す。

<燃焼時間>

 多くの資料が、推進剤の「燃焼時間」(burn time)を 110 [s] と見積もっている。なお、このとき、推進剤の平均燃焼率 mpdotは、12,963 [kg] / 110 [s] = 117.8 [kg/s]推力 T = mpdot x Isp x g = 117.8 x 296 x 9.8 = 341,714 [N] = 34,869 [kgf] となる。

<抗力係数>

 空気抵抗が問題になるのは、(1) ブースト間機体全体の抗力係数 Cdマッハ数 M に応じて変化する。(M<1): 0.5 → (M=1): 1.1 →(M>1): 0.5 ここでは、平均的に 0.8 用いる。

(2) 再突入の間、ブースター部が分離された弾頭部(円錐形、最大径1.32 [m]、長さ 2.5 [m])の抗力係数 Cd : 0.15 - 0.20 (資料2 による)。この抗力係数から、弾頭質量 mw 1,000 [kg]基準断面積 A(=1.37 [m2] として、弾道係数(ballistic coefficent)β は次のように計算される。

 β = mw /(Cd x A) = 1000/((0.15 ~ 0.20 )x 1.37) =  4,866 3,650 [kg/m2]

<発射位置・発射方向>

 地球の東向きの自転は、ミサイルの発射時に初速を与える。自転による地球の円周速度は、赤道で最大、465 [m/s] であり、極に向かうに従て、緯度の余弦に比例して減少する。また、円周速度がミサイルに与える速度は、北極を基準時計回り(clockwise)にとった角度の正弦比例する。

 例えば、北朝鮮東倉里発射基地の場合、緯度は 39°39’36”であるから、自転円周速度は 465 x cos(39.660°)= 358 [m/s] になる。この基地から南に向けて発射した場合、北極から測った角度は 180 °であるから、sin( 180°)= 0 となり、自転速度はミサイルに寄与しない

 舞水端里発射基地(北緯 40°51’20”)の場合:自転円周速度 352 [m/s]、東向きに発射する場合、352 x sin( 90°) = 352 [m/s] が初速として与えられる。

<経路角(上下角)>

  ミサイルの弾道を垂直面内に固定してみた場合、飛翔方向水平面となすを「経路角上下角)」と呼ぶ。経路角 θ は、最初、90 °に設定されて打上げられ、その後、ブースト間に、目標方向に向けて 90°から経路角を減少させ、ブースト終了時に定められた打上角度 θLloft angle)で慣性飛翔に移行する。

 進行方向に垂直な面での運動方程式は、迎角を0と仮定すると、飛翔速度を v として

   v・ dθ/dt = - g cos θ → dθ/dt = - g cos θ / v ・・・・・・・(1)

と表現できるから、重力の損加速度 g によって、上昇につれて、機首は下向きに回転していく。この現象を 「gravity turn」(重力ターン)と呼ぶ。

 ここでは、打上角度 θLと経路角 θ の変化率 θdot を予め設定してシミュレートする。

<空気密度>

 下図に示した「高度対大気密度」の近似式をを使用する。

図1 高度対大気密度
air_density

ρ = 1.225・exp( -0.1051・h ) [kg/m3] ; 0 ≦ 高度 h [km] <10 [Km] ・・・(2)

ρ = 2.020・exp( -0.1558・h ) [kg/m3] ; 10 ≦ 高度 h [km] ・・・・・・・・・・・(3)

 

 

 

戻る:

「はじめに」 へ  Nodong Missile の概要 へ 

進む:

 標準的な弾道 へ  パラメータ分析 へ  パラメータ分析(Isp) へ  パラメータ分析(Payload) へ

 我が国に飛来する Nodong の弾道 へ  首都圏に飛来する弾道(続き) へ  三陸沖に飛来する弾道

 Taepodong_Missile  へ  Taepodong-2 Missile  へ  Taepodong-2 Missile 検証 

 Taepodong-3 の弾道特性を推定する へ  「おわりに」 へ

Top Page へ