馮元偉 吳 強(qiáng)

(92941部隊(duì)94分隊(duì) 葫蘆島 125000)

?

導(dǎo)彈擺式機(jī)動與螺旋機(jī)動突防反導(dǎo)艦炮分析*

馮元偉 吳 強(qiáng)

(92941部隊(duì)94分隊(duì) 葫蘆島 125000)

為探討擺式機(jī)動與螺旋機(jī)動對反導(dǎo)艦炮的突防能力,建立反艦導(dǎo)彈與反導(dǎo)艦炮對抗仿真模型,計(jì)算了多組機(jī)動周期、機(jī)動幅值、艦炮射速組合下,反艦導(dǎo)彈對反導(dǎo)艦炮的突防概率。結(jié)果表明,擺式機(jī)動與螺旋機(jī)動對反艦導(dǎo)彈突防能力均有大幅提升,且螺旋機(jī)動比擺式機(jī)動突防能力更強(qiáng)。而艦炮射速是決定反導(dǎo)艦炮反導(dǎo)能力的最關(guān)鍵因素。

艦炮; 機(jī)動; 仿真; 突防

Class Number TJ391

1 引言

眾所周知,反艦導(dǎo)彈作為海面戰(zhàn)場的主要進(jìn)攻性武器,已成為現(xiàn)代海戰(zhàn)主力。其作戰(zhàn)效能將直接左右海面戰(zhàn)場局面。但面對層層設(shè)防、不斷完善的綜合防御體系,其攻擊能力仍大打折扣。為提高突防能力,反艦導(dǎo)彈近幾年來的發(fā)展可謂突飛猛進(jìn)、異彩紛呈,采用的技術(shù)也五花八門,令人刮目相看。但萬變不離其宗,其突防策略始終離不開兩個(gè)基本點(diǎn):高速突防,機(jī)動突防。目前世界上大部分反艦導(dǎo)彈已具備高速巡航能力,且許多反艦導(dǎo)彈已裝備有機(jī)動突防技術(shù)。如俄羅斯的“寶石”系列反艦導(dǎo)彈就規(guī)劃有蛇形機(jī)動彈道,用以規(guī)避防空導(dǎo)彈的攔截。美國的魚叉、意大利的奧托馬特等反艦導(dǎo)彈則采用躍升俯沖機(jī)動來對付密集陣類近程反導(dǎo)艦炮武器系統(tǒng)[1]。雖然到目前為止,在反艦導(dǎo)彈與反導(dǎo)艦炮的攻防對抗中,除躍升俯沖機(jī)動模型外,尚未見到有其它機(jī)動方式付諸實(shí)施。但是不難想見,如果將蛇形機(jī)動、擺式機(jī)動和螺旋機(jī)動[2～4]等機(jī)動方式用于反艦導(dǎo)彈對反導(dǎo)艦炮的突防,肯定會有意想不到的效果。那么,當(dāng)反艦導(dǎo)彈穿越多重屏障,面對最后一道防線—近程反導(dǎo)艦炮武器系統(tǒng)時(shí),這些機(jī)動行為會為其增加幾成勝算？為一探究竟,本文參考某亞音速反艦導(dǎo)彈與某反導(dǎo)艦炮武器系統(tǒng)之戰(zhàn)術(shù)技術(shù)參數(shù),建立二者攻防對抗仿真模型,系統(tǒng)研究了各種機(jī)動模式下,反艦導(dǎo)彈對反導(dǎo)艦炮的突防能力,得到了一些有價(jià)值的數(shù)據(jù)和結(jié)論[5～6]。本文著重分析比較擺式機(jī)動與螺旋機(jī)動對近程反導(dǎo)艦炮武器系統(tǒng)的突防性能。

2 仿真模型

仿真模型[7～8]包括各種理論值、測量值生成模型,反導(dǎo)艦炮火控解算模型、反導(dǎo)艦炮射擊諸元理論值計(jì)算模型及脫靶量評估模型等,如圖1所示。

圖1 仿真計(jì)算流程

2.1 反艦導(dǎo)彈航跡模擬

2.1.1 坐標(biāo)系

坐標(biāo)系為左手直角坐標(biāo)系,原點(diǎn)O在艦炮回轉(zhuǎn)中心,X、Y軸在水平面內(nèi),X軸指向?qū)椫鬟\(yùn)動方向的反方向,Y軸指向X軸右側(cè),Z軸垂直向上。

2.1.2 擺式機(jī)動

圖2 擺式機(jī)動原理圖

作擺式機(jī)動[9]的導(dǎo)彈,沿X軸反方向運(yùn)動(速度為Vm)的同時(shí),在YOZ面內(nèi)像鐘擺一樣繞某點(diǎn)(Om)擺動。圖2為其運(yùn)動規(guī)跡于YOZ面內(nèi)的投影。其中,Om是擺式運(yùn)動的擺心(等同于鐘擺的擺軸),R為擺長(等同于鐘擺的擺臂長度),ξ為擺角(即擺心-導(dǎo)彈連線與鉛垂線夾角的瞬時(shí)值),Am為擺幅(最大擺動幅度),Cm為擺深(導(dǎo)彈在一個(gè)機(jī)動周期內(nèi)最大高度與最小高度之差),hm為導(dǎo)彈高度最小值。

由圖2可知,導(dǎo)彈質(zhì)心坐標(biāo)(x,y,z)為

x=-Vm·t
y=Rsinξ
z=hm+R-Rcosξ

(1)

式中,t為導(dǎo)彈運(yùn)動時(shí)間。擺角ξ與導(dǎo)彈X軸坐標(biāo)x,x0有關(guān),即

(2)

式中,ξmax是擺角最大值,l是節(jié)距,即導(dǎo)彈在一個(gè)機(jī)動周期內(nèi)沿X軸飛行的距離,x0為起點(diǎn)坐標(biāo)。

顯然,導(dǎo)彈擺式機(jī)動軌跡與擺長R、最大擺角ξmax、節(jié)距l(xiāng)等有關(guān),在仿真過程中調(diào)整這些參數(shù),即能考察它們對艦炮命中概率的影響。圖3是l=2800.0m、R=10.0m、Am=3.0m時(shí),擺式機(jī)動軌跡的模擬圖。

圖3 模擬擺式機(jī)動航跡

2.1.3 螺旋機(jī)動

螺旋機(jī)動[9]也稱滾筒機(jī)動。如圖4所示,Xl是螺旋軸心線,導(dǎo)彈沿Xl軸運(yùn)動(速度為Vm)的同時(shí)繞Xl軸作等速率圓周運(yùn)動。圖中l(wèi)是螺距(導(dǎo)彈在一個(gè)機(jī)動周期內(nèi)沿X軸飛行的距離),R是螺旋半徑(導(dǎo)彈到直線Xl的距離),hm為螺旋軌跡最低點(diǎn)的高度,機(jī)動所需過載以n表示。顯然,導(dǎo)彈螺旋機(jī)動軌跡與螺距、螺旋半徑有關(guān),在仿真過程中調(diào)整這些參數(shù),即能考察它們對反艦導(dǎo)彈突防概率的影響。圖5為某一螺旋機(jī)動軌跡模擬圖(l=1400.0m、R=10.0m)。

圖5 模擬螺旋機(jī)動軌跡

2.1.4 導(dǎo)彈坐標(biāo)探測模型

艦炮武器系統(tǒng)對導(dǎo)彈目標(biāo)的跟蹤探測設(shè)備為光電跟蹤儀或跟蹤雷達(dá),一般情況下,其探測誤差為平穩(wěn)正態(tài)隨機(jī)過程[10]。

2.2 艦艇運(yùn)動姿態(tài)模擬

假設(shè),艦炮所在艦艇理論上作等速直線運(yùn)動。搖擺量、升沉量等姿態(tài)參數(shù)是具有一定周期、一定幅值的正弦量[11],且這些參數(shù)的測量誤差服從正態(tài)分布。

2.3 氣象參數(shù)模擬

彈道氣象參數(shù)理論值,按均勻分布抽樣產(chǎn)生。彈道氣象參數(shù)包括彈丸初速、藥溫、氣溫、空氣密度、風(fēng)速風(fēng)向等,其測量誤差一般服從正態(tài)分布。其值在每個(gè)航次中抽樣設(shè)置一次。彈道氣象參數(shù)在求解艦炮射擊諸元預(yù)測值(火控解算)與理論值時(shí)參與計(jì)算。

2.4 反導(dǎo)艦炮火控解算模型

假定艦炮武器系統(tǒng)火控解算采用卡爾曼濾波預(yù)測目標(biāo)運(yùn)動參數(shù)、未來點(diǎn)及射擊諸元[12]。

火控解算所用數(shù)據(jù),是各測量模型產(chǎn)生的測量值。

2.5 反導(dǎo)艦炮射擊諸元理論值

反導(dǎo)艦炮射擊諸元的理論值,用拉格朗日插值算法在導(dǎo)彈理論航路上求解[13]。參與計(jì)算的導(dǎo)彈坐標(biāo)、艦艇運(yùn)動參數(shù)(如航速、航向、搖擺、升沉等)、彈道氣象參數(shù)等是各仿真模塊生成的理論值。

2.6 炮彈脫靶量

反導(dǎo)艦炮射擊總誤差是火控解算誤差(以Δγ、Δφ表示)、艦炮隨動誤差、炮彈散布誤差三部分的合成。

考慮到隨動誤差與散布誤差的不相關(guān)性,可先綜合再抽樣。記作γb、φb:

(3)

其中,v1、v2是正態(tài)分布隨機(jī)數(shù)[14]。σγ0、σφ0與σγs、σφs分別是炮彈散布與艦炮隨動誤差在方位、高低兩方向的均方差。

則航路上某射擊點(diǎn)角脫靶量(F,E)為

方位:F=Δγ+γb

高低:E=Δφ+φb

角脫靶量與該射擊點(diǎn)未來點(diǎn)斜距離的乘積,即是該發(fā)炮彈相對于導(dǎo)彈幾何中心的線脫靶量,在X坐標(biāo)系(Q平面)[15]內(nèi)以(X1,X2)表示。

2.7 導(dǎo)彈被命中判斷模型

根據(jù)反艦導(dǎo)彈幾何特征及其與反導(dǎo)艦炮的相對位置、相對姿態(tài)等,求得反艦導(dǎo)彈的命中面積及圓形等效命中域半徑r[15]。

將線脫靶量(X1,X2)換算為炮彈至導(dǎo)彈幾何中心的最小距離rx:

(4)

那么,若rx≤r,反艦導(dǎo)彈被炮彈命中;反之未被命中。

2.8 突防概率

設(shè)n為仿真航次數(shù),n1為導(dǎo)彈被命中的航次數(shù),則突防概率P為

(5)

3 計(jì)算方法

計(jì)算方法步驟如下:

1) 參數(shù)賦值初始化。包括反艦導(dǎo)彈的幾何參數(shù);彈道氣象參數(shù)的有效取值范圍,艦艇運(yùn)動特征量(搖擺周期、幅值等),各參數(shù)測量誤差特征量(均值、均方差等),艦炮隨動誤差特征量,炮彈散布特征量,艦炮射擊區(qū)段,以及仿真航次數(shù)等。

2) 設(shè)定機(jī)動方式(擺式機(jī)動或螺旋機(jī)動)、機(jī)動參數(shù)、艦炮射速等,循環(huán)仿真計(jì)算。

3) 每個(gè)航次開始前,集中進(jìn)行一批隨機(jī)抽樣,以產(chǎn)生某些在一個(gè)航次中只需設(shè)置一次的參數(shù)值,如某些運(yùn)動參數(shù)(導(dǎo)彈坐標(biāo)、艦艇姿態(tài)等)的起始值,某些參數(shù)的系統(tǒng)誤差以及氣象參數(shù)等。然后開始一個(gè)航次的模擬。

4) 每個(gè)航次計(jì)算過程中,在艦炮的每個(gè)射擊點(diǎn),各理論值生成模塊模擬產(chǎn)生相應(yīng)參數(shù)(導(dǎo)彈坐標(biāo)、艦艇運(yùn)動姿態(tài)等)的理論值。這些理論值一方面參與計(jì)算射擊諸元理論值;一方面經(jīng)各測量模塊模擬產(chǎn)生相應(yīng)測量值,進(jìn)一步驅(qū)動火控解算模塊濾波求取射擊諸元,并求該射擊點(diǎn)的火控解算誤差。

與此同時(shí),艦炮射擊誤差生成模塊模擬該發(fā)炮彈射擊誤差,經(jīng)脫靶量計(jì)算模塊計(jì)算出脫靶量,最后由命中判別模型判斷該航次命中情況并予統(tǒng)計(jì)。

5) 算完所有航次,統(tǒng)計(jì)突防概率。

6) 調(diào)整機(jī)動方式、機(jī)動參數(shù),或加倍艦炮射速,計(jì)算相應(yīng)條件下的突防概率。

如此循環(huán)反復(fù),直至計(jì)算完所有預(yù)設(shè)條件下的突防概率。

4 計(jì)算結(jié)果

4.1 擺式機(jī)動

1) 令擺深Cm=0.1m保持不變,節(jié)距l(xiāng)分別取350m、700m、1400m、2800m,計(jì)算突防概率隨擺幅Am的變化情況,計(jì)算結(jié)果以圖線表示(圖6)。

圖6 突防概率與擺式機(jī)動之?dāng)[幅之間的關(guān)系曲線

2) 令擺幅Am=3m保持不變,節(jié)距l(xiāng)分別取350m、700m、1400m、2800m,計(jì)算突防概率隨擺深Cm的變化情況(結(jié)果見圖7)。

圖7 突防概率與擺式機(jī)動之?dāng)[深之間的關(guān)系曲線

3) 令擺深Cm=0.1m,節(jié)距l(xiāng)=700m,將艦炮射速設(shè)為原來的2倍,計(jì)算突防概率隨擺幅Am的變化情況,計(jì)算結(jié)果也在圖6中表示。

4) 計(jì)算了擺式機(jī)動過載(如圖8所示)。

圖8 導(dǎo)彈擺式機(jī)動所需過載

4.2 螺旋機(jī)動

1) 計(jì)算了螺距l(xiāng)=700m、1400m、2800m時(shí),突防概率隨螺旋半徑R的變化情況(如圖9所示)。

2) 令螺距l(xiāng)=1400m,艦炮射速加倍,其它條件不變,計(jì)算突防概率隨螺旋半徑R的變化情況(結(jié)果一并示于圖9)。

3) 計(jì)算了螺距l(xiāng)=700m時(shí)螺旋機(jī)動過載(見圖9)。

圖9 突防概率與螺旋機(jī)動之螺旋半徑之間的關(guān)系曲線

5 結(jié)語

1) 無論擺式機(jī)動還是螺旋機(jī)動,均能大幅提升反艦導(dǎo)彈對反導(dǎo)艦炮的突防概率。機(jī)動周期(節(jié)距或螺距)越小,機(jī)動幅度(擺幅或螺旋半徑)越大,突防概率越高。如圖9所示,反艦導(dǎo)彈作螺距700m、螺旋半徑2m的螺旋機(jī)動時(shí),其對反導(dǎo)艦炮的突防概率可達(dá)90%以上。但無機(jī)動時(shí)的突防概率只有20%左右。

2) 由圖5、圖9可見,機(jī)動幅度較小時(shí),隨著機(jī)動幅度的增大,螺旋機(jī)動突防概率增加迅速,而擺式機(jī)動突防概率增加緩慢。因此,螺旋機(jī)動比擺式機(jī)動突防效果更好,所需過載更小,且機(jī)動幅度也小。這在機(jī)動周期較小時(shí)尤為明顯。如節(jié)距和螺距同為700m時(shí),欲達(dá)到80%的突防概率,螺旋機(jī)動所需過載不到1.2g,機(jī)動幅度(螺旋半徑)也在2m以內(nèi)(參加圖9),而擺式機(jī)動所需過載卻大于2g,機(jī)動幅度在3m左右(參加圖6、圖8)。進(jìn)而也容易推知,以螺旋機(jī)動突防時(shí),由于導(dǎo)彈機(jī)動幅度比艦艇輪廓小得多,所以,即使導(dǎo)彈在接近艦艇時(shí)始終保持機(jī)動,也不易丟失或錯(cuò)過艦艇。

3) 由圖7可見,在擺式機(jī)動中,突防概率對擺弧深度很不敏感。而蛇形機(jī)動就是擺弧深度為0時(shí)的擺式機(jī)動,不難想象,擺式機(jī)動與蛇形機(jī)動的突防效果應(yīng)該不相上下。

4) 可以看出,反艦導(dǎo)彈突防概率在很大程度上受制于艦炮射速。如導(dǎo)彈作螺距1400m、螺旋半徑6m的螺旋機(jī)動時(shí),艦炮射速提高1倍可使反艦導(dǎo)彈的突防概率降低1倍以上(參見圖9)。因此艦炮射速是決定反導(dǎo)艦炮反導(dǎo)能力的最關(guān)鍵因素。

[1] 宋立,李勇,馬良.反艦導(dǎo)彈對近程反導(dǎo)艦炮的機(jī)動突防模型分析[J].艦船電子工程,2009,29(9):122-124.

[2] 顧文錦,趙紅超.超聲速反艦導(dǎo)彈末端非平面機(jī)動研究[J].飛行力學(xué),2003,21(3):36-40.

[3] 蓋玉華,顧文錦,趙紅超,等.反艦導(dǎo)彈的主要突防技術(shù)及其效能[J].海軍航空工程學(xué)院學(xué)報(bào),2008,23(2):168-170.

[4] 顧文錦,王鳳蓮.反艦導(dǎo)彈末端螺旋機(jī)動的控制研究[J].戰(zhàn)術(shù)導(dǎo)彈控制技術(shù),2003(2):2-7.

[5] 馮元偉,徐國亮,劉奎永.目標(biāo)擺式機(jī)動對反導(dǎo)艦炮射擊精度的影響[J].指揮控制與仿真,2010,32(2):101-103.

[6] 馮元偉,金振中,郜彪.導(dǎo)彈蛇形機(jī)動與擺式機(jī)動突防反導(dǎo)艦炮仿真[J].現(xiàn)代防御技術(shù),2010,38(4):22-25.

[7] 邱志明.艦炮武器系統(tǒng)分析[M].北京:兵器工業(yè)出版社,1999:113-129.

[8] 程健慶.仿真建模技術(shù)在新型艦炮武器系統(tǒng)試驗(yàn)中的應(yīng)用[J].計(jì)算機(jī)仿真,2002,19(1):54-58.

[9] 張翼飛,鄧方林,顧文錦.反艦導(dǎo)彈末端機(jī)動彈道設(shè)計(jì)及彈道仿真[J].計(jì)算機(jī)仿真,2005,22(2):19-22.

[10] 楊文安,韓成標(biāo).基于白噪聲的跟蹤雷達(dá)量測誤差建模與仿真[J].吉林大學(xué)學(xué)報(bào),2005,23(6):621-627.

[11] 孫東彥,潘紅華,徐明.某型艦炮火控系統(tǒng)建模與仿真[J].火力與指揮控制,2006,31(12):37-40.

[12] 王基組.艦炮火控系統(tǒng)原理[M].北京:海潮出版社,1999:83-193.

[13] 鄭樹敏,馮元偉,霍文志,等.GJB254A-2000艦炮武器系統(tǒng)設(shè)計(jì)定型試驗(yàn)規(guī)程[S].北京:國家軍用標(biāo)準(zhǔn),2000:39-42.

[14] 朱元昌.火控系統(tǒng)的數(shù)字仿真[M].北京:國防工業(yè)出版社,1992:39-48.

[15] 李小元,吳宇瓊,周衛(wèi)平,等.GJB 592.6-89,艦炮武器系統(tǒng)射擊效力評定確定目標(biāo)命中面積的方法[S].北京:國家軍用標(biāo)準(zhǔn),1988:1-5.

Analysis of Penetration of Cycloid Maneuver and Spiral Maneuver Missile Against Anti-missle Naval Gun

FENG Yuanwei WU Qiang

(Unit 94, No. 92941 Troops of PLA, Huludao 125000)

In order to discuss the penetration ability of cycloid maneuver and spiral maneuver against anti-missle naval gun, the combat engagement models of missile and gun were established,and the penetration probability in different conditions of two maneuver forms were computed respectively. The results show the penetration ability of anti-ship missile can be significantly increased by cycloid maneuver and spiral maneuver. Compared with the cycloid maneuver, spiral maneuver penetration ability is better. While the effective way to improve the anti-missile capability of the short range anti-missle naval gun weapon system is to increase its firing rate.

naval gun, maneuver, simulation, penetration

2014年10月3日,

2014年11月27日

馮元偉,男,高級工程師,研究方向:艦炮武器系統(tǒng)試驗(yàn)。吳強(qiáng),男,工程師,研究方向:艦炮武器系統(tǒng)試驗(yàn)。

TJ391

10.3969/j.issn1672-9730.2015.04.031