彈炮結(jié)合武器系統(tǒng)末端反導(dǎo)毀傷評(píng)估方法研究*

謝曉方，王誠(chéng)成，張龍杰，孫濤，曹建

(海軍航空工程學(xué)院 兵器科學(xué)與技術(shù)系，山東 煙臺(tái)　264001)

彈炮結(jié)合武器系統(tǒng)末端反導(dǎo)毀傷評(píng)估方法研究*

謝曉方，王誠(chéng)成，張龍杰，孫濤，曹建

(海軍航空工程學(xué)院 兵器科學(xué)與技術(shù)系，山東 煙臺(tái)264001)

以反艦導(dǎo)彈在彈炮結(jié)合武器系統(tǒng)防御末端艦炮攔截情況下的毀傷評(píng)估為研究方向，提出了一種在小口徑艦炮作用下，基于命中能力及命中后毀傷效能的反艦導(dǎo)彈毀傷評(píng)估方法。從統(tǒng)計(jì)學(xué)角度分析了彈丸從發(fā)射到命中后毀傷過(guò)程，建立了相應(yīng)模型并進(jìn)行了仿真分析。本研究對(duì)反艦導(dǎo)彈在近程反導(dǎo)武器系統(tǒng)作用下的易損性評(píng)估提供了一種新的研究思路。

反艦導(dǎo)彈；毀傷評(píng)估；彈炮結(jié)合；小口徑艦炮

0　引言

針對(duì)典型目標(biāo)的毀傷評(píng)估一直以來(lái)都是各國(guó)軍事專家研究的熱點(diǎn)。國(guó)外就目標(biāo)毀傷評(píng)估研究工作起步較早，美軍ARL及BRL實(shí)驗(yàn)室現(xiàn)已研究并開發(fā)了多種典型目標(biāo)毀傷評(píng)估模型及程序[1-3]，如VISIA，MUVES，COVART，SHOTGEN等，國(guó)內(nèi)該領(lǐng)域的專家、學(xué)者也相繼就不同目標(biāo)的毀傷評(píng)估開展了相關(guān)研究[4-5]。

艦載彈炮結(jié)合武器系統(tǒng)結(jié)合艦空導(dǎo)彈與小口徑艦炮的優(yōu)勢(shì)，是防御反艦導(dǎo)彈的有效手段，一般采用艦空導(dǎo)彈實(shí)施中遠(yuǎn)程攔截，小口徑艦炮實(shí)施末端攔截[6-7]。可以說(shuō)艦炮是水面艦艇近程防御的最后一道屏障，其對(duì)反艦導(dǎo)彈的毀傷評(píng)估對(duì)于攻防雙方而言都具有十分重要的意義，也是當(dāng)前國(guó)內(nèi)外防空反導(dǎo)技術(shù)領(lǐng)域的重要研究方向。本文針對(duì)彈炮結(jié)合武器系統(tǒng)對(duì)反艦導(dǎo)彈的末端防御問(wèn)題，重點(diǎn)研究艦炮作用下對(duì)反艦導(dǎo)彈的毀傷評(píng)估。鑒于實(shí)彈打靶實(shí)驗(yàn)存在的研究成本高、周期長(zhǎng)以及無(wú)法全面反映各種彈目交匯條件等問(wèn)題，利用計(jì)算機(jī)進(jìn)行反艦導(dǎo)彈易損性仿真評(píng)估的作用日益顯著，合理有效的仿真評(píng)估方法是實(shí)現(xiàn)反艦導(dǎo)彈毀傷評(píng)估的關(guān)鍵。

1　反艦導(dǎo)彈毀傷評(píng)估模型

欲就擔(dān)負(fù)艦載彈炮結(jié)合武器系統(tǒng)末端防御的小口徑艦炮對(duì)反艦導(dǎo)彈的毀傷進(jìn)行合理有效的評(píng)估，需要同時(shí)考慮2個(gè)主要方面因素，分別為艦炮對(duì)反艦導(dǎo)彈的命中能力評(píng)估以及反艦導(dǎo)彈受彈后的毀傷效果評(píng)估，二者缺一不可。反艦導(dǎo)彈易損性評(píng)估總體架構(gòu)如圖1所示。

圖1　毀傷評(píng)估整體架構(gòu)Fig.1　Framework of damage estimation

1.1對(duì)反艦導(dǎo)彈的命中能力評(píng)估

1.1.1命中能力評(píng)估仿真流程

通過(guò)建立小口徑艦炮對(duì)空射擊的火控解算模型以及誤差傳遞模型，對(duì)反艦導(dǎo)彈的受彈情況進(jìn)行統(tǒng)計(jì)分析，并記錄導(dǎo)彈受彈時(shí)刻的有關(guān)信息，在此基礎(chǔ)上統(tǒng)計(jì)分析小口徑艦炮對(duì)反艦導(dǎo)彈的命中能力。

系統(tǒng)整體仿真框架如圖2所示。

圖2　命中能力評(píng)估仿真框架Fig.2　Simulation framing of shooting ability

在小口徑艦炮對(duì)反艦導(dǎo)彈命中能力研究過(guò)程中，首先建立反艦導(dǎo)彈和艦炮彈丸運(yùn)動(dòng)模型，在此基礎(chǔ)上對(duì)艦炮反導(dǎo)射擊各個(gè)環(huán)節(jié)中所涉及的誤差進(jìn)行分析計(jì)算，最后根據(jù)命中判讀模型對(duì)艦炮的反導(dǎo)能力進(jìn)行統(tǒng)計(jì)分析。

1.1.2艦炮彈丸運(yùn)動(dòng)模型

艦炮彈丸發(fā)射后空中運(yùn)動(dòng)過(guò)程中只受到重力和空氣阻力，阻力對(duì)彈丸質(zhì)心速度大小和方向的影響是通過(guò)阻力加速度ax來(lái)體現(xiàn)的，即有

(1)

S=πd2/4.

(2)

結(jié)合式(1)和(2)，有

cH(h)G(cx,vc)vc，

(3)

如圖3所示，在彈丸運(yùn)動(dòng)坐標(biāo)系Sxczc下，彈丸質(zhì)心運(yùn)動(dòng)的矢量方程為

(4)

圖3　彈丸運(yùn)動(dòng)坐標(biāo)系Fig.3　Projectile motion coordinate system

將彈丸質(zhì)心運(yùn)動(dòng)的矢量方程向各坐標(biāo)軸投影，由式(3)，(4)，得到

(5)

式中：Dc為彈丸運(yùn)動(dòng)距離。初始條件

(6)

式中：vc0為彈丸初始飛散速度；θc0為初始射角。

由式(5)和式(6)建立的是無(wú)散布誤差情況下艦炮彈丸的質(zhì)心運(yùn)動(dòng)方程，實(shí)際情況下，由于彈丸發(fā)射過(guò)程中藥溫、初速、大氣擾動(dòng)等隨機(jī)因素的影響，彈丸相對(duì)于射擊線會(huì)存在一個(gè)隨機(jī)散布，具體表現(xiàn)為射擊方位角方向和高低角方向的隨機(jī)散布。在仿真分析中，根據(jù)彈丸隨機(jī)散布誤差的大小，通過(guò)在理論射擊線的基礎(chǔ)上附加一個(gè)方位角誤差和高低角誤差來(lái)實(shí)現(xiàn)對(duì)彈丸隨機(jī)散布的模擬。

此外，在建立彈丸質(zhì)心運(yùn)動(dòng)方程過(guò)程中，假設(shè)彈丸的運(yùn)動(dòng)完全位于二維平面，即彈丸運(yùn)動(dòng)坐標(biāo)系Sxczc內(nèi)，但是彈丸在實(shí)際運(yùn)動(dòng)過(guò)程中是高速旋轉(zhuǎn)的，由于旋轉(zhuǎn)過(guò)程中會(huì)受到空氣阻力的影響，彈丸會(huì)一定程度上偏離Sxczc平面，由于在實(shí)際射擊過(guò)程中系統(tǒng)能夠較好的對(duì)偏移量進(jìn)行修正，為了簡(jiǎn)化分析過(guò)程，不再考慮彈丸運(yùn)動(dòng)的偏移。

1.1.3艦炮對(duì)空射擊誤差分析

彈炮結(jié)合武器系統(tǒng)中小口徑艦炮對(duì)空攔截作戰(zhàn)的整個(gè)過(guò)程主要可由目標(biāo)坐標(biāo)觀測(cè)、目標(biāo)參數(shù)估計(jì)、目標(biāo)提前點(diǎn)解算以及隨動(dòng)系統(tǒng)瞄準(zhǔn)射擊等4個(gè)環(huán)節(jié)組成，在不同的處理環(huán)節(jié)都會(huì)有誤差產(chǎn)生，以γ表示方位角方向的誤差，φ表示高低角方向的誤差，σR表示距離測(cè)量誤差，則各環(huán)節(jié)誤差可分別表示為目標(biāo)坐標(biāo)觀測(cè)誤差e1(γ1,φ1,σR)、目標(biāo)參數(shù)濾波誤差e2、隨動(dòng)系統(tǒng)跟蹤誤差e3(γ3,φ3)(主要包括靜態(tài)瞄準(zhǔn)誤差e3s(γ3s,φ3s)和動(dòng)態(tài)跟蹤誤差e3d(γ3d,φ3d))、目標(biāo)提前點(diǎn)解算誤差e4(σx4,σy4,σz4)以及彈丸散布誤差e5(γ5,φ5)等。

在彈丸射擊過(guò)程中，不同環(huán)節(jié)誤差甚至同一誤差不同時(shí)段所服從分布的形式都各不相同，需要大量樣本進(jìn)行統(tǒng)計(jì)分析，本文采用Monte Carlo仿真方法，以每一次彈丸射擊為樣本，對(duì)各個(gè)環(huán)節(jié)誤差進(jìn)行多航次仿真。

以艦炮對(duì)導(dǎo)彈目標(biāo)的觀測(cè)與估計(jì)為例，在處理目標(biāo)坐標(biāo)觀測(cè)誤差時(shí)，小口徑艦炮對(duì)空攔截射擊過(guò)程中，需要同時(shí)測(cè)量目標(biāo)的斜距、方位角以及高低角。對(duì)于角度測(cè)量誤差γ1和φ1，一般情況下服從正態(tài)分布，對(duì)于距離測(cè)量誤差σR，由于計(jì)數(shù)器的計(jì)數(shù)精度和A/D轉(zhuǎn)換過(guò)程的量化誤差是主要的影響因素，可以采用均勻分布假設(shè)對(duì)其進(jìn)行模擬。

假設(shè)在載體甲板坐標(biāo)系下，以球坐標(biāo)表示的目標(biāo)真實(shí)參數(shù)為(Rw,αw,βw)，則目標(biāo)距離測(cè)量值Rεw的仿真計(jì)算公式為

Rεw=Rw+uσRσR.

(7)

目標(biāo)方位角和高低角觀測(cè)值的仿真計(jì)算公式為

(8)

式中：uσR為-1和1之間服從均勻分布的隨機(jī)數(shù)；uγ1和uφ1為服從N(0,1)分布的隨機(jī)數(shù)。

式(7)和式(8)構(gòu)成了載體甲板坐標(biāo)系下目標(biāo)參數(shù)的球坐標(biāo)觀測(cè)結(jié)果(Rεw,αεw,βεw)。目標(biāo)參數(shù)估計(jì)誤差主要產(chǎn)生于對(duì)目標(biāo)坐標(biāo)觀測(cè)值進(jìn)行濾波估計(jì)的過(guò)程中，其誤差值與濾波算法有關(guān)，需要根據(jù)相應(yīng)算法(如最小二乘濾波、Wiener濾波、Kalman濾波等)進(jìn)行動(dòng)態(tài)解算，這里限于篇幅不再展開。主要誤差模型計(jì)算過(guò)程詳見文獻(xiàn)[8]與文獻(xiàn)[9]。

1.1.3反艦導(dǎo)彈命中判讀模型分析

對(duì)于導(dǎo)彈彈體，采用“圓柱+圓錐”的幾何模型進(jìn)行等效，彈體部分采用圓柱體進(jìn)行等效，導(dǎo)彈頭部則采用圓錐體進(jìn)行等效，如圖4所示。

圖4　導(dǎo)彈受彈等效模型Fig.4　Equivalent model of missile body

圖4中：dm為彈體直徑；lm為導(dǎo)彈整體長(zhǎng)度；lh為導(dǎo)彈圓錐頭部的長(zhǎng)度。

則在AB段彈體的橫截圓半徑為

rAB=0.5dm.

(9)

在BD段，導(dǎo)彈的橫截圓半徑隨彈體縱軸坐標(biāo)xd的變化規(guī)律為

(10)

式中：xd∈[0.5lm-lh,0.5lm]。

在艦炮反導(dǎo)作戰(zhàn)能力統(tǒng)計(jì)分析中，需要根據(jù)導(dǎo)彈具體的受彈位置確定受彈部件，按照某型超音速反艦導(dǎo)彈的彈體布局，從整體上將彈體劃分為3個(gè)主要艙段：分別為導(dǎo)引與控制系統(tǒng)艙段、戰(zhàn)斗部系統(tǒng)艙段以及推進(jìn)系統(tǒng)艙段。導(dǎo)引和控制系統(tǒng)艙段主要由導(dǎo)引系統(tǒng)和控制系統(tǒng)的關(guān)鍵部件構(gòu)成；戰(zhàn)斗部系統(tǒng)艙段主要由位于彈體中前部區(qū)域的引信和戰(zhàn)斗部構(gòu)成；推進(jìn)系統(tǒng)主要由位于彈體中部區(qū)域的油箱、輸油管和中后部區(qū)域的發(fā)動(dòng)機(jī)構(gòu)成。3個(gè)艙段的比例近似為3∶1∶4。

假設(shè)某時(shí)刻導(dǎo)彈的理論點(diǎn)坐標(biāo)為(Rj1,αj1,βj1)，實(shí)際射擊點(diǎn)坐標(biāo)為(Rj2,αj2,βj2)，則實(shí)際射擊點(diǎn)坐標(biāo)在艦炮瞄準(zhǔn)線坐標(biāo)系Sxsyszs下的坐標(biāo)為

(11)

1.2反艦導(dǎo)彈受彈后的毀傷效果評(píng)估

在小口徑艦炮命中反艦導(dǎo)彈的前提下，進(jìn)一步分析反艦導(dǎo)彈受彈后的毀傷情況，同時(shí)研究艦炮彈丸與反艦導(dǎo)彈遭遇時(shí)刻的交匯角以及反艦導(dǎo)彈關(guān)鍵部位材料及厚度對(duì)反艦導(dǎo)彈毀傷程度的影響。

以小口徑艦炮彈丸命中反艦導(dǎo)彈時(shí)的入射角、遭遇速度、命中位置等信息作為初始條件，借助顯式動(dòng)力分析軟件ANSYS/LS-DYNA對(duì)小口徑艦炮彈丸侵徹反艦導(dǎo)彈進(jìn)行有限元分析，受彈后毀傷效果評(píng)估總體架構(gòu)如圖5所示。

(1) 建立彈丸及反艦導(dǎo)彈各個(gè)系統(tǒng)所在艙段的有限元模型。模型尺寸及材料盡量與原型一致，以增加仿真數(shù)據(jù)的真實(shí)性。

(2) 在導(dǎo)彈及彈丸有限元模型建立的基礎(chǔ)上，就導(dǎo)彈在動(dòng)態(tài)條件下彈丸對(duì)導(dǎo)彈各個(gè)艙段侵徹進(jìn)行有限元仿真，各個(gè)艙段仿真分多角度(0°～90°范圍)、多部位進(jìn)行。

(3) 觀察侵徹過(guò)程，分析在相應(yīng)入射部位下不同入射角度及速度對(duì)侵徹效果的影響，對(duì)最終得到的仿真結(jié)果進(jìn)行毀傷判讀。

圖5　受彈后毀傷效果評(píng)估總體架構(gòu)Fig.5　Framework of damage effect estimation after being shot

通過(guò)多次仿真，從統(tǒng)計(jì)學(xué)角度得出導(dǎo)彈相應(yīng)部位不同相對(duì)速度及交匯角度下受彈后的毀傷概率。

需要特別指出的是，對(duì)于侵徹毀傷仿真而言，建立細(xì)致的有限元模型以及選擇正確的材料模型參數(shù)至關(guān)重要，二者關(guān)系到仿真結(jié)果的有效性及真實(shí)性。因而在建模過(guò)程中，一方面要使得物理模型盡量貼近原型，另一方面還要明確反艦導(dǎo)彈及艦炮彈丸的材料參數(shù)、選擇正確的材料模型。由于小口徑艦炮彈丸侵徹反艦導(dǎo)彈對(duì)于二者金屬材料而言是一種大應(yīng)變、高應(yīng)變速率、高溫變形的過(guò)程，因而仿真主要采用Johnson-Cook本構(gòu)模型和GRUNEISEN狀態(tài)方程[10]。

2　仿真分析

假設(shè)小口徑艦炮對(duì)小航路捷徑反艦導(dǎo)彈進(jìn)行攔截射擊，對(duì)象彈參考某大中型超音速反艦導(dǎo)彈技術(shù)參數(shù)：彈長(zhǎng)9.4 m，彈徑0.8 m，末端近似水平勻速運(yùn)動(dòng)，速度Ma數(shù)2.5。取小口徑艦炮誤差參數(shù)：e1=(0.5,0.5)mrad，σR=±5 m，e3s=(1.0,1.0)mrad，e5=(1.3,1.3)mrad，其余誤差動(dòng)態(tài)解算。首先采用Monte Carlo方法對(duì)對(duì)象彈進(jìn)行受彈數(shù)量統(tǒng)計(jì)分析，表1為艦炮對(duì)不同航路捷徑ASM攔截時(shí)發(fā)射的彈丸數(shù)。

仿真進(jìn)行10 000次，初始仿真步長(zhǎng)0.001 s，得到不同航路捷徑下對(duì)象彈不同部位及整體的受彈情況，如圖6所示。

圖6　命中彈數(shù)統(tǒng)計(jì)Fig.6　Statistics of hit counts

從圖6中可以看出，最小航路捷徑不同時(shí)，導(dǎo)彈的受彈數(shù)量也不同。在400～800 m區(qū)間導(dǎo)彈的平均受彈統(tǒng)計(jì)數(shù)量大于6枚，在此區(qū)間內(nèi)艦炮最大可持續(xù)射擊時(shí)間較長(zhǎng)，隨動(dòng)系統(tǒng)動(dòng)態(tài)跟蹤誤差相對(duì)較低，導(dǎo)彈的受彈數(shù)量也相對(duì)較多；當(dāng)最小航路捷徑大于800 m時(shí)，艦炮最大可持續(xù)射擊時(shí)間明顯縮短，并且目標(biāo)距離艦炮位置較遠(yuǎn)，命中彈數(shù)逐漸減小。而從導(dǎo)彈不同部位受彈數(shù)量來(lái)看，在不同最小航路捷徑下，推進(jìn)系統(tǒng)艙段受彈數(shù)量最多，導(dǎo)引及控制系統(tǒng)艙段次之，戰(zhàn)斗部艙受彈數(shù)量最少。

命中后的毀傷效能評(píng)估涉及彈丸及反艦導(dǎo)彈有限元模型建立、反艦導(dǎo)彈受侵徹區(qū)域劃分以及在相應(yīng)的相對(duì)速度下對(duì)不同部位、不同入射角度的侵徹毀傷進(jìn)行有限元仿真，需要進(jìn)行大量的仿真統(tǒng)計(jì)得出相應(yīng)的毀傷數(shù)據(jù)，本文限于篇幅選取兩種典型侵徹仿真過(guò)程圖示分析，不對(duì)具體模型的建立展開詳細(xì)敘述。

仿真條件設(shè)定：艦炮彈丸彈芯材料為鎢合金，導(dǎo)彈蒙皮為鈦合金，厚度為4 mm；導(dǎo)引與控制系統(tǒng)艙內(nèi)部系統(tǒng)等效為壁厚為4 mm的硬鋁立方箱體。戰(zhàn)斗部殼體采用鈦合金，殼體厚度為18 mm，整流罩艙與后部戰(zhàn)斗部艙之間有5 mm厚硬鋁隔板。導(dǎo)彈與彈丸接觸瞬間，導(dǎo)彈速度Ma為2.5(850 m/s)，彈丸速度為1 200 m/s。計(jì)算中采用的主要材料參數(shù)[11-12]如表2所示。

(1) 沿導(dǎo)彈縱軸45°方向侵徹戰(zhàn)斗部艙

由于導(dǎo)彈為軸對(duì)稱結(jié)構(gòu)，這里取1/4導(dǎo)彈戰(zhàn)斗部艙模型進(jìn)行仿真分析。圖7～8分別為時(shí)間T=0，400 μs時(shí)的彈丸侵徹有限元仿真圖。

從圖中可見看出，彈丸穿透戰(zhàn)斗部艙鈦合金蒙皮后進(jìn)入殼體內(nèi)部，沖擊戰(zhàn)斗部側(cè)壁殼體，仿真顯示彈丸擊中戰(zhàn)斗部側(cè)壁后對(duì)戰(zhàn)斗部殼體造成擊穿毀傷。

(2) 彈丸沿導(dǎo)彈縱軸10°方向侵徹頭部導(dǎo)引與控制系統(tǒng)艙

圖9～11分別為彈丸沿導(dǎo)彈縱軸10°方向侵徹時(shí)，時(shí)間T=0，500及1 000 μs時(shí)的彈丸侵徹有限元

表1　艦炮發(fā)射彈丸數(shù)Table 1　Number of projectiles that naval gun fires

表2　主要材料模型參數(shù)Table 2　Main parameters of material model

仿真圖。從圖中可以看出，彈丸穿透導(dǎo)彈蒙皮后進(jìn)入殼體內(nèi)部，對(duì)內(nèi)部器件等效殼體進(jìn)行侵徹毀傷，彈丸在侵徹過(guò)程中彈體發(fā)生一定程度的變形，運(yùn)動(dòng)軌跡也受到導(dǎo)彈內(nèi)部器件等效殼體的影響。仿真顯示彈丸在擊穿導(dǎo)引與控制系統(tǒng)后繼續(xù)穿透隔板沖擊戰(zhàn)斗部，對(duì)導(dǎo)引與控制系統(tǒng)艙及戰(zhàn)斗部造成擊穿毀傷。

圖7　T=0 μs侵徹仿真Fig.7　Penetrating simulation as T=0 μs

圖8　T=400 μs侵徹仿真Fig.8　Penetrating simulation as T=400 μs

圖9　T=0 μs侵徹仿真Fig.9　Penetrating simulation as T=0 μs

圖10　T=500 μs侵徹仿真 Fig.10　Penetrating simulation as T=500 μs

圖11　T=1 000μs侵徹仿真Fig.11　Penetrating simulation as T=1 000 μs

通過(guò)以上2種典型侵徹毀傷仿真可以看出，艦炮彈丸命中反艦導(dǎo)彈頭部后，在入射角度及部位適當(dāng)?shù)那闆r下在毀傷導(dǎo)引與控制系統(tǒng)后可能繼續(xù)對(duì)其后面的戰(zhàn)斗部系統(tǒng)造成毀傷，因而在對(duì)某一系統(tǒng)進(jìn)行毀傷評(píng)估時(shí)在獨(dú)立建模仿真的同時(shí)應(yīng)注意彈丸侵徹運(yùn)動(dòng)的連貫性導(dǎo)致的對(duì)后續(xù)系統(tǒng)可能造成的毀傷，使得毀傷評(píng)估仿真數(shù)據(jù)更加真實(shí)、可靠。

3　結(jié)束語(yǔ)

近程反導(dǎo)武器對(duì)反艦導(dǎo)彈毀傷評(píng)估是一項(xiàng)系統(tǒng)而復(fù)雜的基礎(chǔ)性工作，在各國(guó)海軍越來(lái)越重視發(fā)展反艦導(dǎo)彈及近程反導(dǎo)武器系統(tǒng)的情況下，建立一套合理有效的反艦導(dǎo)彈毀傷評(píng)估方法無(wú)疑具有十分重要的實(shí)際意義。本文提出一種基于命中能力及命中后毀傷效能綜合評(píng)估的反艦導(dǎo)彈易損性評(píng)估方法，對(duì)有效評(píng)價(jià)彈炮結(jié)合武器系統(tǒng)對(duì)反艦導(dǎo)彈的毀傷效果以及導(dǎo)彈整體的易損性分析提供了理論參考。

[1]盧永剛,錢立新,楊云斌.目標(biāo)易損性/戰(zhàn)斗部威力評(píng)估方法[J].彈道學(xué)報(bào)，2005,17(1):46-49. LU Yong-gang, QIAN Li-xin, YANG Yun-bin. Vulnerability/Lethality Assessment Method Based on Virtual Model[J]. Journal of Ballistics，2005,17(1):46-49.

[2]ABELL J M, ROACH L K, STARKS M W. Degraded States Vulnerability analysis [R]. BRL-TR-3010,U.S. Army Ballistic Research Laboratory , Aberdeen Proving Ground, MD, June 1989.

[3]DEITZ P H, Ozolins A. Computer Simulations of the Abrams Live-Fire Field testing [R]. BRL-MR-3750, 1988.

[4]裴揚(yáng).飛機(jī)易損性建模方法研究及DMECA軟件開發(fā)[D].西安：西北工業(yè)大學(xué)，2003. PEI Yang. Modeling Method Research of Airplane Vulnerability & Software development of DMECA[D].Xi′an: NWPU,2003.

[5]李向東, 蘇義嶺, 韓永要. 導(dǎo)彈目標(biāo)在破片式戰(zhàn)斗部作用下的易損性評(píng)估[J]. 爆炸與沖擊，2007，27(5):468-472. LI Xiang-dong, SU Yi-ling , HAN Yong-yao. Vulnerability Assessment of the Missile Subjected to the Fragment Warhead[J]. Explosion and Shock Waves, 2007，27(5):468-472.

[6]彭峰生, 廖振強(qiáng). 彈炮結(jié)合防空武器系統(tǒng)毀傷概率分析與仿真[J]. 火力與指揮控制，2007,32(3):62-64. PENG Feng-sheng, LIAO Zhen-qiang. Kill Probability Analysis and Simulation of Integrated Missile and Anti-aircraft Gun Weapon System [J]. Fire Control and Command Control，2007,32(3)：62-64.

[7]嚴(yán)代彪, 鄧建輝, 王航宇, 等. 彈炮結(jié)合武器系統(tǒng)效能評(píng)估[J]. 火力與指揮控制，2003,28(4):38-41. YAN Dai-biao, DENG Jian-hui, WANG Hang-yu. Effectiveness Evaluation of Gun System with Missile [J]. Fire Control and Command Control,2003,28(4):38-41.

[8]張龍杰, 謝曉方, 孫濤, 等. 艦載CIWS對(duì)小航路捷徑超音速導(dǎo)彈的攔截射擊分析[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2013, 30(1): 123-130. ZHANG Long-jie, XIE Xiao-fang, SUN Tao, et al. Intercept Analysis of Shipborne CIWS to Supersonic Missiles on a Shortcut Flight Course [J]. Journal of Harbin Engineering University, 2013, 30(1): 123-130.

[9]張龍杰, 謝曉方, 孫濤, 等. 艦艇搖擺對(duì)小口徑艦炮射擊精度的影響分析[J]. 火炮發(fā)射與控制學(xué)報(bào), 2013, 34(1): 11-15. ZHANG Long-jie, XIE Xiao-fang, SUN Tao, et al. Effects Analysis of Ship Swaying on Small Naval Gun Firing Accuracy [J]. Journal of Artillery Launch and Control, 2013, 34(1): 11-15.

[10]JOHNSON G R, COOK W H. A Constitutive Model and Data for Metals Subjected to Large Strains, High Strain Rates and High Temperatures[C]∥Proceedings of the Seventh International Symposium on Ballistics. The Hague, Netherlands, 1983:541-547.

[11]范亞夫,段祝平.Johnson-Cook材料參數(shù)的實(shí)驗(yàn)測(cè)定[J].力學(xué)與實(shí)踐,2003,25(5):40-43. FAN Ya-fu, DUAN Zhu-ping. Cylinder Explosive Test and Material Model of Johnson-cook[J].Mechanics in Engineering, 2003,25(5):40-43.

[12]蘭彬，文鶴鳴.鎢合金長(zhǎng)桿彈侵徹半無(wú)限鋼靶的數(shù)值模擬及分析[J].高壓物理學(xué)報(bào)，2008,25(3):245-249. LAN Bin, WEN He-ming. Numerical Simulation and Analysis of the Penetration of Tungsten-Alloy Long Rod into Semi-Infinite Armor Steel Targets[J].Chinese Journal of High Pressure Physics，2008,25(3):245-249.

Research on Terminal Antimissile Damage Estimate Method of Integrated Missile and Antiaircraft Gun

XIE Xiao-fang, WANG Cheng-cheng, ZHANG Long-jie,SUN Tao, CAO Jian

(NAEI, Department of Ordnance Science and Technology, Shandong Yantai 264001,China)

The damage estimation of anti-ship missile caused by naval gun under the circumstance of end defense of integrated missile and antiaircraft gun weapon system is taken as research direction. A kind of damage estimation method of anti-ship missile in the effect of close-in warship weapon system is put forward based on the abilities of hitting the target and damage efficiency after being hit. The damage process is analyzed from the perspective of statistical and relevant models are built and simulated. The research offers a new research direction to the vulnerability estimation of anti-ship missile in the effect of close-in weapon system (CIWS).

anti-ship missile; damage estimation; integrated missile and antiaircraft gun; small caliber naval gun

2014-09-26；

2014-10-20

中國(guó)博士后科學(xué)基金(2013T60923;20110491883)

謝曉方(1962-)，男，河北承德人。教授，研究方向?yàn)槲淦飨到y(tǒng)建模與仿真，虛擬現(xiàn)實(shí)技術(shù)。

通信地址：264000山東煙臺(tái)2馬路188號(hào)2系1007室E-mail：497990518@qq.com

10.3969/j.issn.1009-086x.2015.04.008

TJ761；N945.11

A

1009-086X(2015)-04-0043-07