姜穎資,王偉力,傅　磊,黃雪峰

(海軍航空工程學(xué)院,山東煙臺　264001)

離散桿和沖擊波對超音速導(dǎo)彈復(fù)合毀傷研究*

姜穎資,王偉力,傅磊,黃雪峰

(海軍航空工程學(xué)院,山東煙臺264001)

摘要:為研究離散桿和沖擊波復(fù)合作用下對超音速導(dǎo)彈的毀傷效應(yīng),對離散桿戰(zhàn)斗部沖擊起爆和結(jié)構(gòu)毀傷超音速導(dǎo)彈進(jìn)行了仿真研究,獲得了離散桿戰(zhàn)斗部沖擊起爆超音速導(dǎo)彈戰(zhàn)斗部的臨界距離以及彈目距離、交匯角對毀傷超音速導(dǎo)彈發(fā)動機(jī)艙的影響規(guī)律。結(jié)果表明:離散桿戰(zhàn)斗部只有在桿條與沖擊波的復(fù)合作用下,才可能引爆超音速導(dǎo)彈戰(zhàn)斗部;離散桿戰(zhàn)斗部桿條和沖擊波疊加作用下對超音速導(dǎo)彈發(fā)動機(jī)艙的結(jié)構(gòu)毀傷效應(yīng)明顯優(yōu)于桿條獨(dú)立作用下的情況。

關(guān)鍵詞:離散桿戰(zhàn)斗部;沖擊波;超音速導(dǎo)彈;沖擊起爆;結(jié)構(gòu)毀傷;復(fù)合毀傷效應(yīng)

0引言

離散桿戰(zhàn)斗部是一種新興的防空反導(dǎo)戰(zhàn)斗部,按桿條飛散姿態(tài)是否可控,分為姿態(tài)可控和無控離散桿戰(zhàn)斗部。離散桿戰(zhàn)斗部殺傷單元相比破片戰(zhàn)斗部,質(zhì)量大、切割目標(biāo)能力強(qiáng);相比連續(xù)桿戰(zhàn)斗部,桿條初速高、殺傷半徑大,近些年已成為戰(zhàn)斗部研制的熱點[1-4]。離散桿戰(zhàn)斗部主要依靠桿條所攜帶的動能毀傷來襲目標(biāo),當(dāng)離散桿戰(zhàn)斗部距離目標(biāo)一定距離內(nèi)爆炸時,爆炸產(chǎn)生的沖擊波和爆轟產(chǎn)物的作用將不可忽略,此時,戰(zhàn)斗部在桿條和爆炸產(chǎn)生的沖擊波、爆轟產(chǎn)物的復(fù)合作用下毀傷目標(biāo),對目標(biāo)毀傷模式同桿條單獨(dú)作用下有很大不同。

目前,速度更快、機(jī)動性和突防能力更強(qiáng)的超音速導(dǎo)彈正廣泛取代亞音速導(dǎo)彈;能否有效攔截、毀傷來襲超音速導(dǎo)彈,已成為制約甚至決定戰(zhàn)爭進(jìn)程和結(jié)局的重要因素。

鑒于以上考慮,文中采用LS-DYNA軟件,研究一種新穎的雙束可控離散桿戰(zhàn)斗部桿條和沖擊波復(fù)合作用下對超音速導(dǎo)彈的毀傷效應(yīng)。

1計算模型及物理參數(shù)

1.1　計算模型

雙束可控離散桿戰(zhàn)斗部,桿條外形為圓柱形,半徑為2.5 mm,長度為10 mm,每束桿條數(shù)量為72根,離散桿戰(zhàn)斗部桿條排列如圖1所示。

雙束可控離散桿戰(zhàn)斗部的兩束桿條,在向外飛散的同時分別向相反的方向自旋,當(dāng)飛散到兩束桿條形成的“殺傷環(huán)”半徑約為1.5 m時,兩束桿條均旋轉(zhuǎn)了90°,形成了兩束準(zhǔn)連續(xù)的圓環(huán)。兩束桿條形成的“殺傷環(huán)”半徑均約為0.75 m和1.5 m時,兩束桿條飛散情況分別如圖2和圖3所示。

圖1　桿條排列圖

圖2　約0.75 m時桿條飛散情況

圖3　約1.5 m時桿條飛散情況

離散桿戰(zhàn)斗部對超音速導(dǎo)彈主要可造成對戰(zhàn)斗部裝藥的沖擊起爆以及對導(dǎo)彈蒙皮和殼體撕裂、凹陷等結(jié)構(gòu)毀傷。文中以超音速導(dǎo)彈戰(zhàn)斗部和發(fā)動機(jī)艙為毀傷目標(biāo),研究離散桿戰(zhàn)斗部桿條和沖擊波、爆轟產(chǎn)物復(fù)合作用下對超音速導(dǎo)彈的毀傷效應(yīng)。真實超音速導(dǎo)彈結(jié)構(gòu)十分復(fù)雜,研究毀傷效應(yīng)一般對其進(jìn)行簡化,建立等效靶模型。

建立的典型超音速導(dǎo)彈戰(zhàn)斗部和發(fā)動機(jī)艙等效靶分別如圖4和圖5所示。其中,戰(zhàn)斗部等效靶為圓柱形結(jié)構(gòu),由戰(zhàn)斗部殼體和主裝藥兩部分組成,戰(zhàn)斗部殼體外半徑為25 cm,厚度為10 mm;發(fā)動機(jī)艙等效靶為雙圓環(huán)結(jié)構(gòu),外圓環(huán)模擬導(dǎo)彈蒙皮,外半徑為25 cm,厚度為4 mm;內(nèi)圓環(huán)模擬發(fā)動機(jī)殼體,外半徑為23.9 cm,厚度為2.5 mm。

圖4　戰(zhàn)斗部和發(fā)動機(jī)艙等效靶

采用ANSYS/LS-DYNA軟件的流固耦合算法,數(shù)值模擬離散桿戰(zhàn)斗部兩束桿條受爆轟驅(qū)動在空氣中飛散以及對超音速導(dǎo)彈戰(zhàn)斗部和發(fā)動機(jī)艙的毀傷效應(yīng)。計算模型中離散桿戰(zhàn)斗部主裝藥、空氣使用歐拉單元;離散桿戰(zhàn)斗部殼體、桿條、前后端蓋和超音速導(dǎo)彈戰(zhàn)斗部殼體、導(dǎo)彈蒙皮、發(fā)動機(jī)殼體使用拉格朗日單元。在空氣域邊界處施加透射邊界,計算單位采用cm-g-μs。

1.2　物理參數(shù)

離散桿戰(zhàn)斗部主裝藥為PBX9404炸藥,采用High-explosive-burn模型和JWL狀態(tài)方程描述[5]:

(1)

式中:P為爆轟壓力;V是相對體積;E是炸藥單位體積內(nèi)能;ω、A、B、R1、R2為材料常數(shù)。各參數(shù)如表1所示[6]。

表1　PBX9404炸藥材料參數(shù)

離散桿桿條材料為10號鋼,彈體蒙皮材料為2A12鋁合金,戰(zhàn)斗部殼體和發(fā)動機(jī)殼體材料為4340鋼,材料模型均采用Johnson-Cook模型。各材料主要參數(shù)如表2所示[7-8]。

戰(zhàn)斗部主裝藥為B炸藥,其狀態(tài)方程采用Lee-Tarver點火增長方程。反應(yīng)率方程形式為:

(2)

式中:F為反應(yīng)分?jǐn)?shù);t為時間;ρ0為初始密度;ρ為當(dāng)前密度;p為壓力;I、G1、G2、b、x、a、c、d、y、e、g和z都是常數(shù)。B炸藥Lee-Tarver點火增長參數(shù)如表3所示[9-10]。

表3　B炸藥點火增長參數(shù)

2計算結(jié)果及分析

2.1　沖擊起爆計算結(jié)果及分析

離散桿戰(zhàn)斗部桿條和沖擊波都可對導(dǎo)彈目標(biāo)造成毀傷。桿條對導(dǎo)彈目標(biāo)的毀傷主要表現(xiàn)為對導(dǎo)彈殼體結(jié)構(gòu)的“切割”破壞,沖擊波對導(dǎo)彈目標(biāo)的毀傷主要體現(xiàn)于對導(dǎo)彈殼體結(jié)構(gòu)的擠壓、變形破壞。當(dāng)離散桿戰(zhàn)斗部距離導(dǎo)彈目標(biāo)較近爆炸時,對導(dǎo)彈目標(biāo)造成的毀傷,不是桿條和沖擊波分別對導(dǎo)彈目標(biāo)造成的毀傷的簡單相加,而是存在疊加毀傷效應(yīng)。此時,桿條和沖擊波復(fù)合作用下對導(dǎo)彈目標(biāo)的毀傷模式、毀傷效果,明顯有別于兩種毀傷元單獨(dú)作用下對導(dǎo)彈目標(biāo)的毀傷效應(yīng)。

文中按離散桿戰(zhàn)斗部與超音速導(dǎo)彈彈目交匯角0°、30°、60°,計算離散桿戰(zhàn)斗部沖擊起爆超音速導(dǎo)彈戰(zhàn)斗部的臨界距離。彈目交匯時,離散桿戰(zhàn)斗部運(yùn)動速度為800 m/s,來襲超音速導(dǎo)彈速度為700 m/s。通過數(shù)值模擬計算,0°、30°、60°時,離散桿戰(zhàn)斗部沖擊起爆超音速導(dǎo)彈戰(zhàn)斗部的臨界距離分別為30 cm、37 cm、42 cm(超音速導(dǎo)彈戰(zhàn)斗部表面距離離散桿戰(zhàn)斗部爆心垂直距離)。0°、30°、60°時,離散桿戰(zhàn)斗部臨界沖擊起爆超音速導(dǎo)彈戰(zhàn)斗部的工況,離散桿與超音速導(dǎo)彈戰(zhàn)斗部初始相對位置以及超音速導(dǎo)彈被沖擊起爆后壓力云圖,分別如圖5~圖7所示。

圖5　交匯角0°時,炸藥沖擊起爆后壓力云圖

圖6　交匯角30°時,炸藥沖擊起爆后壓力云圖

圖7　交匯角60°時,炸藥沖擊起爆后壓力云圖

從圖5~圖7可以看出,離散桿戰(zhàn)斗部沖擊起爆超音速導(dǎo)彈戰(zhàn)斗部,是由沖擊波、爆轟產(chǎn)物和桿條的復(fù)合作用實現(xiàn)的。離散桿戰(zhàn)斗部爆炸后,桿條在沖擊波和爆轟產(chǎn)物的推動下高速向外飛散,由于沖擊波的初速遠(yuǎn)大于桿條的飛散速度,當(dāng)離散桿戰(zhàn)斗部距離超音速導(dǎo)彈戰(zhàn)斗部較近時,沖擊波先一步到達(dá)超音速導(dǎo)彈戰(zhàn)斗部,造成戰(zhàn)斗部殼體局部凹陷,凹陷處殼體對裝藥產(chǎn)生擠壓破壞作用;而后,桿條和爆轟產(chǎn)物幾乎同時到達(dá)戰(zhàn)斗部殼體表面,此時殼體凹陷處,在桿條和爆轟產(chǎn)物的切割和擠壓下,殼體發(fā)生嚴(yán)重的撕裂、破壞,桿條和爆轟產(chǎn)物在撕裂殼體后繼續(xù)切割、擠壓裝藥,當(dāng)造成戰(zhàn)斗部裝藥內(nèi)最大壓力超過B炸藥臨界起爆壓力5.63 GPa時,戰(zhàn)斗部裝藥內(nèi)部壓力迅速升高,直至形成穩(wěn)定爆轟。

從各彈目交匯角下,離散桿戰(zhàn)斗部沖擊起爆超音速導(dǎo)彈戰(zhàn)斗部的臨界距離來看,60°的工況最遠(yuǎn),這是由于彈目交匯角60°時,超音速導(dǎo)彈具有約693 m/s迎向離散桿戰(zhàn)斗部的運(yùn)動速度。三種彈目交匯角下,對超音速導(dǎo)彈戰(zhàn)斗部最遠(yuǎn)的臨界沖擊起爆距離僅為42 cm,說明離散桿戰(zhàn)斗部沖擊起爆導(dǎo)彈戰(zhàn)斗部的能力較弱,只有在近距離時桿條與沖擊波復(fù)合作用下才有可能沖擊起爆超音速導(dǎo)彈戰(zhàn)斗部。

2.2　結(jié)構(gòu)毀傷計算結(jié)果及分析

離散桿桿條在空氣中飛行其速度不斷衰減,且桿條姿態(tài)、分布密度也時刻發(fā)生變化;沖擊波在空氣中衰減非常迅速,對不同距離處的目標(biāo)的毀傷差異非常明顯;因此,離散桿戰(zhàn)斗部距離超音速導(dǎo)彈目標(biāo)不同距離爆炸時對其毀傷效應(yīng)存在差異。

文中按離散桿戰(zhàn)斗部與超音速導(dǎo)彈發(fā)動機(jī)艙兩者相距0.5 m、1 m、1.5 m、2 m(超音速導(dǎo)彈發(fā)動機(jī)艙表面距離離散桿戰(zhàn)斗部爆心垂直距離),彈目交匯角0°、30°、60°,研究桿條和沖擊波對超音速導(dǎo)彈發(fā)動機(jī)艙的復(fù)合毀傷效應(yīng)。相距0.5 m和1.5 m,彈目交匯角0°、30°、60°時,超音速導(dǎo)彈發(fā)動機(jī)艙彈體蒙皮和發(fā)動機(jī)殼體毀傷效應(yīng)應(yīng)力云圖如圖8~圖13所示(每幅圖左圖為彈體蒙皮毀傷效應(yīng)應(yīng)力云圖,右圖為發(fā)動機(jī)殼體毀傷效應(yīng)應(yīng)力云圖);各工況下超音速導(dǎo)彈發(fā)動機(jī)艙彈體蒙皮和發(fā)動機(jī)殼體破壞情況如表4所示。

對比圖8~圖10與圖11~圖13可以看出,離散桿戰(zhàn)斗部與超音速導(dǎo)彈發(fā)動機(jī)艙相距0.5 m時,對超音速導(dǎo)彈發(fā)動機(jī)艙的毀傷效應(yīng)明顯強(qiáng)于1.5 m時。0.5 m時,各彈目交匯工況,發(fā)動機(jī)艙彈體蒙皮和發(fā)動機(jī)殼體,都發(fā)生了大面積的破裂和變形;而1.5 m時,各彈目交匯工況,發(fā)動機(jī)艙彈體蒙皮被離散桿桿條切割出長條狀切口,切口長度和寬度明顯小于0.5 m時,發(fā)動機(jī)殼體局部凹陷,但均未破裂。

圖8　距離0.5 m,交匯角0°時,發(fā)動機(jī)艙毀傷應(yīng)力云圖

圖9　距離0.5 m,交匯角30°時,發(fā)動機(jī)艙毀傷應(yīng)力云圖

圖10　距離0.5 m,交匯角60°時,發(fā)動機(jī)艙毀傷應(yīng)力云圖

圖11　距離1.5 m,交匯角0°時,發(fā)動機(jī)艙毀傷應(yīng)力云圖

圖12　距離1.5 m,交匯角30°時,發(fā)動機(jī)艙毀傷應(yīng)力云圖

圖13　距離1.5 m,交匯角60°時,發(fā)動機(jī)艙毀傷應(yīng)力云圖

從表4計算結(jié)果可以看出,離散桿戰(zhàn)斗部與超音速導(dǎo)彈發(fā)動機(jī)艙距離相同的各工況對比:彈目交匯角0°的工況,發(fā)動機(jī)艙彈體蒙皮破裂長度最長,但破裂寬度最小;彈目交匯角60°的工況,發(fā)動機(jī)艙彈體蒙皮破裂長度最小,但破裂寬度最大,這是由于,彈目交匯角60°時,離散桿戰(zhàn)斗部與超音速導(dǎo)彈發(fā)動機(jī)艙平行于發(fā)動機(jī)艙橫軸方向相對速度最大。對發(fā)動機(jī)殼體的毀傷,離散桿戰(zhàn)斗部與超音速導(dǎo)彈發(fā)動機(jī)艙距離1.5 m、2 m的各工況,發(fā)動機(jī)殼體均未破裂;距離0.5 m、1 m的工況,相同距離的各工況對比,發(fā)動機(jī)殼體破裂長度和寬度均是彈目交匯角30°的工況最大,彈目交匯角0°的工況次之,彈目交匯角60°的工況最小。對比離散桿戰(zhàn)斗部相對超音速導(dǎo)彈發(fā)動機(jī)艙不同距離的各工況,可以發(fā)現(xiàn),離散桿戰(zhàn)斗部對超音速導(dǎo)彈發(fā)動機(jī)艙彈體蒙皮和發(fā)動機(jī)殼體的毀傷隨著兩者相對距離的增大而迅速變小,毀傷效應(yīng)由離散桿桿條和沖擊波的復(fù)合作用,逐漸變?yōu)闂U條的獨(dú)立切割作用。

表4　發(fā)動機(jī)艙蒙皮和發(fā)動機(jī)殼體裂口長度和寬度

綜合上述分析結(jié)果,離散桿戰(zhàn)斗部對超音速導(dǎo)彈發(fā)動機(jī)艙結(jié)構(gòu)破壞能力較強(qiáng),特別是在近距離情況下,離散桿桿條和沖擊波的復(fù)合作用下,可造成超音速導(dǎo)彈發(fā)動機(jī)艙彈體蒙皮和發(fā)動機(jī)殼體大面積撕裂、變形,超音速導(dǎo)彈發(fā)動機(jī)將會完全喪失其功能。離散桿桿條和沖擊波的疊加毀傷能力,與離散桿戰(zhàn)斗部與目標(biāo)的相對距離密切相關(guān)。離散桿戰(zhàn)斗部爆炸后產(chǎn)生的沖擊波峰值超壓為[11]:

(3)

從上式可以看出,沖擊波超壓與爆點距離的平方成反比,隨著爆點距離的增大,沖擊波的毀傷能力將迅速衰減。離散桿戰(zhàn)斗部與超音速導(dǎo)彈發(fā)動機(jī)艙相距一定距離以后,沖擊波的能量將可以忽略,離散桿戰(zhàn)斗部主要依靠桿條所攜帶的動能切割毀傷發(fā)動機(jī)艙,對發(fā)動機(jī)艙的毀傷效應(yīng)明顯弱于沖擊波與桿條的復(fù)合毀傷效應(yīng)。

3結(jié)論

通過數(shù)值模擬計算以及結(jié)果分析,可以得出如下結(jié)論:

1)離散桿戰(zhàn)斗部對超音速導(dǎo)彈的毀傷主要體現(xiàn)為結(jié)構(gòu)毀傷,很難引爆超音速導(dǎo)彈戰(zhàn)斗部。只有在離散桿戰(zhàn)斗部與超音速導(dǎo)彈戰(zhàn)斗部相距很近時,沖擊波與離散桿戰(zhàn)斗部桿條的復(fù)合作用下,才有可能引爆超音速導(dǎo)彈戰(zhàn)斗部。

2)離散桿戰(zhàn)斗部對超音速導(dǎo)彈發(fā)動機(jī)艙的結(jié)構(gòu)毀傷能力與兩者的相對距離密切相關(guān),在沖擊波與離散桿桿條疊加效應(yīng)明顯的距離內(nèi),毀傷效應(yīng)明顯優(yōu)于離散桿桿條獨(dú)立作用下的情況。在離散桿戰(zhàn)斗部與超音速導(dǎo)彈發(fā)動機(jī)艙距離相同時,是彈目交匯角30°的工況對超音速導(dǎo)彈發(fā)動機(jī)毀傷效果最佳。

參考文獻(xiàn):

[1]李龍俊, 董素榮, 陳秀文, 等. 不同起爆方式下離散桿戰(zhàn)斗部爆炸驅(qū)動桿條的數(shù)值研究 [J]. 彈箭與制導(dǎo)學(xué)報, 2009, 29(4): 103-105.

[2]孫傳杰, 路中華, 盧永剛, 等. 可控旋轉(zhuǎn)離散桿空間運(yùn)動分析 [J]. 爆炸與沖擊, 2008, 28(4): 378-382.

[3]宋浦, 余建斌, 梁安定, 等. 桿式戰(zhàn)斗部對柱殼裝藥的毀傷實驗研究 [J]. 含能材料, 2004, 12(4): 249-251.

[4]呂紅超, 張新偉. 可控離散桿式戰(zhàn)斗部設(shè)計技術(shù) [J]. 航空兵器, 2011(3): 36-38.

[5]陳軍, 曾代朋, 孫承緯. JB-9014炸藥超壓爆轟產(chǎn)物的狀態(tài)方程 [J]. 爆炸與沖擊, 2010, 30(6): 583-586.

[6]章冠人, 陳大年. 凝聚炸藥起爆動力學(xué) [M]. 北京: 國防工業(yè)出版社, 1991.

[7]申志強(qiáng). 穿甲子彈偏心入射陶瓷復(fù)合靶板的侵徹機(jī)理研究 [D]. 國防科學(xué)技術(shù)大學(xué), 2006.

[8]杜茂華, 王偉力, 盧明章. EFP對反艦導(dǎo)彈發(fā)動機(jī)艙的侵徹毀傷機(jī)理研究 [J]. 工程爆破, 2011, 17(4): 9-13.

[9]賈憲振, 陳松, 楊建, 等. 雙破片同時撞擊對B炸藥沖擊起爆的數(shù)值模擬研究 [J]. 高壓物理學(xué)報, 2011, 25(5): 469-474.

[10]陳衛(wèi)東, 張忠, 劉家良. 破片對屏蔽炸藥沖擊起爆的數(shù)值模擬和分析 [J]. 兵工學(xué)報, 2009, 30(9): 1187-1191.

[11]陳文, 張慶明. 爆炸沖擊波下導(dǎo)彈結(jié)構(gòu)動力分析模型初步研究 [J]. 爆炸與沖擊, 2009, 29(2): 190-204.

收稿日期:2014-12-27

作者簡介:姜穎資(1982-),男,山東威海人,博士研究生,研究方向:武器系統(tǒng)作戰(zhàn)效能與毀傷能力。

中圖分類號:TJ410.3

文獻(xiàn)標(biāo)志碼:A

Research on Composite Damage Effect on Supersonic Missile by
Discrete Rod Warhead and Blast Wave

JIANG Yingzi,WANG Weili,FU Lei,HUANG Xuefeng

(Naval Aeronautical and Astronautical University, Shandong Yantai 264001, China)

Abstract:To study the composite damage effect on supersonic missile under discrete rod and blast wave, numerical simulation studies about impact initiation and structural damage supersonic missile warhead by discrete rod warhead was conducted. Both the critical distance that discrete rod warhead impact initiating supersonic missile warhead and the rules that the relative distance and the intersection angle of the missile and the target damaged the engine cabin of supersonic anti-ship missile were obtained. The results show that only under co-action of bar and blast wave, the supersonic missile warhead can be detonated. The structural damage effect of the discrete rod warhead on the engine cabin of the supersonic missile by co-action of bar and blast wave is much better than action of bar only.

Keywords:discrete rod warhead; blast wave; supersonic missile; impact initiation; engine cabin; composite damage effect