劉 琛, 李 元, 李燕華, 溫玉全

(北京理工大學(xué)爆炸科學(xué)與技術(shù)國家重點實驗室, 北京 100081)

1 引 言

偏心起爆圓柱形定向戰(zhàn)斗部因其結(jié)構(gòu)簡單,并能明顯提高目標(biāo)方向破片速度,一直以來受到廣泛的關(guān)注[1-5],而棱柱形定向戰(zhàn)斗部采用偏心起爆方式后,破片飛散速度及分布密度均可得到顯著提高[6]。美英的“可編程集成彈藥艙”(PIOS)計劃[7]研究了一種自適應(yīng)六棱柱定向戰(zhàn)斗部,其毀傷效率高于常規(guī)戰(zhàn)斗部。Sam Waggener[8]針對六棱柱定向戰(zhàn)斗部破片飛散存在覆蓋盲區(qū)的問題,研究了起爆方式對破片飛散角的影響,通過調(diào)整起爆方式可使破片束的飛散角度發(fā)生7°～8°的偏移,從而實現(xiàn)360°可打擊范圍。李元等[6]對比研究了棱柱型和圓柱型偏心起爆定向戰(zhàn)斗部的毀傷效能,得出棱柱型偏心起爆定向戰(zhàn)斗部能形成高密度的破片束,棱上偏心八點同時起爆時的平均速度和平均動能增益分別達(dá)到24.13%和54.52%,且在彈靶距離40m時仍有較強(qiáng)的毀傷效果。

雖然國外很早就提出了棱柱型定向戰(zhàn)斗部的概念[7-9],但是其具體性能參數(shù)并無報道,而國內(nèi)對棱柱形定向戰(zhàn)斗部研究較少,其破片飛散規(guī)律尚未完全掌握。本研究利用LS-DYNA軟件對六棱柱定向戰(zhàn)斗部進(jìn)行了數(shù)值模擬分析,結(jié)合爆炸邏輯網(wǎng)絡(luò)定向起爆系統(tǒng)的邏輯功能、爆轟波輸出特點,綜合研究了起爆方式對棱柱形定向戰(zhàn)斗部破片速度及分布密度的影響,意圖為棱柱型定向戰(zhàn)斗部及其起爆系統(tǒng)設(shè)計提供理論依據(jù)。

2 六棱形戰(zhàn)斗部模型建立

2.1 模型結(jié)構(gòu)及相關(guān)參數(shù)

六棱柱定向戰(zhàn)斗部模型結(jié)構(gòu)如圖1所示,主要由炸藥、空氣、襯筒和球形破片四部分組成,模型結(jié)構(gòu)未顯示空氣域結(jié)構(gòu),破片直徑6 mm,戰(zhàn)斗部總外徑122 mm,襯筒厚度2 mm,戰(zhàn)斗部總長度240 mm。有限元模型中空氣域直徑300 mm。裝填破片數(shù)量2577個,破片相互交錯均勻排布于襯筒之上。

圖1 六棱柱定向戰(zhàn)斗部模型結(jié)構(gòu)

Fig.1 The model structure of prismatic aimable warhead

戰(zhàn)斗部主裝藥為TNT/RDX 40/60炸藥,炸藥材料采用高能炸藥燃燒材料模型和JWL狀態(tài)方程[10]描述。材料模型和狀態(tài)方程相關(guān)參數(shù)如表1所示[11]。

空氣采用空材料模型和多項式狀態(tài)方程[10],具體參數(shù)如表2所示[11]。

表1 裝藥性能參數(shù)及狀態(tài)方程參數(shù)

Table 1 Parameters of charge property and equation of state

density/g·cm-3pcr/GPaDcr/m·s-1A/GPaB/GPaR1R2ω1.71729.57890524.237.6784.21.10.34

表2 空氣材料參數(shù)

Table 2 Parameters of air material

density/g·cm-3C0C1C2C3C4C50.0012900000.40.4

戰(zhàn)斗部襯筒采用LY12-CZ硬鋁材料,該材料在爆轟作用下發(fā)生彈塑性變形,所以采用各向同性硬化彈塑性模型描述,具體參數(shù)如表3所示[11]。

破片用93W合金材料制成,采用剛性材料模型,以確保破片在驅(qū)動過程中不發(fā)生變形和損壞,其材料參數(shù)如表4所示[11]。

表3 襯筒材料參數(shù)

Table 3 Parameters of linear material

density/g·cm-3elasticmodulus/GPapoissonratioyieldlimit/GPahardeningindex2.7720.330.310.7

表4 破片材料參數(shù)

Table 4 Parameters of fragment material

density/g·cm-3Young'smodulus/GPaPoissonratio17.511170.22

2.2 起爆方式設(shè)定

爆炸邏輯網(wǎng)絡(luò)定向起爆系統(tǒng)可以實現(xiàn)側(cè)向一線或多線多點同步輸出、序貫輸出等模式,因此在起爆方式設(shè)定時,首先選定五種起爆方式,即面偏心一線起爆(圖2a)、面偏心兩線起爆(圖2b)、面偏心三線起爆(圖2c)、棱偏心一線起爆(圖2d)、棱偏心兩線起爆(圖2e),起爆點位置如圖2所示。先在同步起爆方式下進(jìn)行數(shù)值計算,再考慮延時起爆方式。延時起爆時,起爆點由上至下依次延時相同的時間間隔。研究中還同時設(shè)置了端面中心一點起爆方式(圖2f),以便和偏心起爆方式作對比,分析毀傷性能參數(shù)變化情況。

a. eccentric one line b. eccentric two lines initiation on face initiation on face

c. eccentric three lines d. eccentric one line initiation on face initiation on arris

e. eccentric two lines f. center one point initiation on arris initiation on end face

圖2 六棱柱定向戰(zhàn)斗部起爆點布置

Fig.2 Initiator layouts of prismatic aimable warhead

3 計算結(jié)果及數(shù)據(jù)分析

利用LS-DYNA計算180 μs,統(tǒng)計此時刻多點同步起爆和延時起爆下定向區(qū)域內(nèi)破片沿軸向及徑向的飛散速度及方向角,將其作為評價定向戰(zhàn)斗部威力的性能參數(shù)。

3.1 同步起爆3.1.1 破片飛散規(guī)律

圖3為該六棱戰(zhàn)斗部起爆180 μs時分別采用不同起爆方式所形成的典型破片飛散圖,由圖3可知，破片飛散匯聚成6個破片束。與端面中心一點起爆相比,采用偏心起爆后,破片沿徑向飛散呈現(xiàn)明顯的不對稱性,如圖3a、圖3c、圖3e所示。起爆點對側(cè)破片速度明顯提高,但各方向破片數(shù)量并未發(fā)生改變。采用棱偏心一線起爆或面偏心兩線起爆后,左側(cè)上方破片速度明顯高于右側(cè),如圖3b、圖3d所示。

a. eccentric one line b. eccentric two lines initiation on face initiation on face

c. eccentric three lines d. eccentric one line initiation on face initiation on arris

e. eccentric two lines f. center one point initiation on arris initiation on end face

圖3 180 μs不同起爆方式下破片飛散情況

Fig.3 Fragment dispersion situation of different initiation ways at 180 μs

3.1.2 破片飛散速度

由于每個面上排布的破片個數(shù)較多,為了便于量化比較,取起爆點對側(cè)平面內(nèi)中間一行的破片作為研究對象。棱柱戰(zhàn)斗部對側(cè)平面上共布置了78行破片,中間行為第39、40行,從左到右依次編號為1-11,其中中間6號破片所在位置稱為中間區(qū)域,兩邊1、11號破片所在位置稱為邊緣部位,破片排列位置如圖4所示。不同起爆方式下的破片飛散速度統(tǒng)計如圖5所示。

圖4 破片排布示意圖

Fig.4 Fragment configuration Diagram

圖5 不同起爆方式下的徑向破片飛散速度

Fig.5 Fragment dispersion velocities in radial direction under different initiation ways

由圖5可以看出,起爆方式對破片飛散速度影響較大。采用端面中心一點起爆方式起爆時,破片最大飛散速度為1935.49 m·s-1; 采用面偏心一線同步起爆時,其破片最大飛散速度為2326.06 m·s-1,比端面中心一點起爆時提高20.18%; 采用面偏心兩線同步起爆時,其破片最大飛散速度為2231.91 m·s-1,提高15.31%; 采用面偏心三線同步起爆時,其破片最大飛散速度為2371.33 m·s-1,提高22.52%; 采用棱偏心一線同步起爆時,其破片最大飛散速度為2272.25 m·s-1,提高17.40%; 采用棱偏心兩線同步起爆時,其破片最大飛散速度為2510.83 m·s-1,提高29.73%。

從以上數(shù)據(jù)可看出,采用偏心起爆時的破片速度明顯高于采用中心起爆時的破片速度,其中采用棱偏心兩線同步起爆時的破片速度增益最大。從圖5可以看出,不論采取何種起爆方式,都是中間區(qū)域的破片飛散速度高于邊緣部位的破片飛散速度。這是因為中間區(qū)域有爆轟波匯聚,破片在爆轟波的匯聚作用下飛散速度會高一些,邊緣部分因為有稀疏波的作用,破片速度會降低。此外,由于爆轟波對破片的加載過程不一致,使得破片沿徑向的飛散速度分布存在較大差異。

3.1.3 破片徑向飛散角

不同起爆方式下的破片徑向飛散角如表5所示,表中參數(shù)定義見圖6,其中β1為中心起爆時的徑向飛散角,β2為面偏心兩線同步起爆或棱偏心一線同步起爆時的徑向飛散角,α為偏轉(zhuǎn)角。

表5 不同起爆方式對破片徑向飛散角的影響

Table 5 Influence of different initiation ways on fragment dispersion angle in radial direction

initiationwayβ/(°)α/(°)centerinitiation27.330eccentriconelineinitiationonface24.920eccentrictwolinesinitiationonface26.630eccentricthreelinesinitiationonface24.360eccentriconelineinitiationonarris22.820eccentrictwolinesinitiationonarris25.580

Note:βis the radial dispersion angle,αis the deflection angle.

圖6 徑向飛散角

Fig.6 Dispersion angles in radial direction

由表5可得出,采用端面中心單點起爆時,破片束的徑向飛散角在27°左右,偏轉(zhuǎn)角0°; 采用棱偏心一線起爆時,破片束的徑向飛散角在23°左右,偏轉(zhuǎn)角是0°; 采用其它偏心起爆方式時,破片束的徑向飛散角在25°左右,偏轉(zhuǎn)角也是0°。這表明起爆方式對破片束的徑向飛散角有一定的影響,但對偏轉(zhuǎn)角沒有影響。這與采用偏心起爆可以明顯提高定向區(qū)域內(nèi)的破片飛散速度,但對目標(biāo)方向的破片密度并無影響的研究結(jié)果一致[6]。

3.2 棱偏心兩線序貫延時起爆

由于采用棱偏心兩線同步起爆時破片速度增益最大,故在棱偏心兩線序貫延時起爆方式下,研究序貫延時時間對破片飛散的影響。設(shè)相鄰兩起爆點之間的爆轟波傳播時間為t=L/D,L為相鄰兩起爆點間距,D為裝藥爆速,則設(shè)定起爆延時時長分別為t/4、2t/4、3t/4,對不同延時時長下的計算效果進(jìn)行比較。

3.2.1 爆轟波壓力及破片飛散規(guī)律

圖7為棱偏心兩線序貫延時t/4起爆時爆轟波傳播情況。

圖8為棱偏心兩線序貫延時t/4起爆180 μs時破片飛散圖,可以看出,經(jīng)過延時處理后,起爆點對側(cè)的破片受到較大影響,在飛散過程中會產(chǎn)生一定角度的傾斜,這是因為隨起爆點起爆時間的不同,其對側(cè)破片向外飛散時間及距離也不同,最終使破片飛散角發(fā)生改變。

圖7 序貫延時爆轟波傳播

Fig.7 Detonation wave propagation under sequential delay detonation

圖8 180 μs序貫延時破片飛散情況

Fig.8 Fragment dispersion situation under sequential delay at 180 μs

3.2.2 延時起爆對破片飛散速度的影響

為便于分析,取起爆點對側(cè)面中間一列的破片速度作為比較對象。對側(cè)面共排11列破片,取中間第6列,依次從上到下提取破片速度。中間列在不同延時時長下的破片飛散速度分布如圖9所示。

圖9 序貫延時起爆對軸向破片飛散速度的影響

Fig.9 Effect of sequential delay initiation on fragment dispersion velocities in axial direction

由圖9可看出,對棱偏心兩線起爆進(jìn)行序貫延時處理后,其破片飛散速度會發(fā)生一定的改變。采用棱偏心兩線同步起爆時,其破片最大飛散速度為2510.83 m·s-1; 對其進(jìn)行t/4序貫延時處理后,其破片最大飛散速度為2444.82 m·s-1; 對其進(jìn)行2t/4序貫延時處理后,其破片最大飛散速度為2420.89 m·s-1; 對其進(jìn)行3t/4序貫延時處理后,其破片最大飛散速度為2260.71 m·s-1。以上結(jié)果表明,經(jīng)過序貫延時處理后,破片飛散速度與棱偏心兩線同步起爆時相比明顯減小,且破片飛散速度隨延時時長增加而減小,這是因為延時處理使得爆轟波之間的碰撞條件發(fā)生了改變。

3.2.3 延時起爆對破片軸向方向角的影響

仍取對側(cè)面中間第6列作為研究對象,量取不同延時時長下該列上每個破片的軸向方向角,如圖10所示。

圖10 延時起爆對破片軸向方向角的影響

Fig.10 Effect of delay initiation on fragment direction angles in axial direction

由圖10可看出,對棱偏心兩線起爆進(jìn)行序貫延時處理后,其軸向方向角整體增大,且隨延時時長增加,角度改變量逐漸增大。將序貫延時起爆時各個位置破片軸向飛散角與棱偏心兩線同步起爆時同一位置的破片軸向飛散角相減,得出破片軸向飛散角的改變量,將每種起爆方式下飛散角改變量的最大值作為對比統(tǒng)計量,如表6所示。

表6 不同序貫延時下的破片軸向方向角改變量

Table 6 Change quantity in axial direction angle of fragment under different sequential delays

delaytime0t/42t/43t/4changeofaxialdirectionangle/(°)04.475.616.77

從表6可以看出,隨著延時時長的增加,軸向方向飛散角改變量增大,最大可達(dá)6.77°,表明序貫延時起爆能夠改變破片束飛散方向,可用于精確控制破片的飛散方向,有助于提高對目標(biāo)的毀傷概率。

4 結(jié) 論

(1) 采用偏心起爆后,破片沿徑向飛散呈現(xiàn)明顯的不對稱性,破片束基本分布在25°左右的徑向范圍內(nèi);

(2) 相對于中心起爆,采用棱偏心兩線同步起爆時對目標(biāo)方位破片最大飛散速度增益可達(dá)29.73%;

(3) 序貫延時起爆對破片的軸向飛散角有一定的影響,隨著延時時長的增加,軸向方向角最大可達(dá)6.77°,可用于精確控制破片的飛散方向,提高對目標(biāo)的毀傷概率。

參考文獻(xiàn):

[1] Fairlie G E , Marriott C O. Computer modeling of full size fragmenting aimable warheads using AUTODAN-3D[C]∥ 17th International Symposium on Ballistics,Midrand, South Africa, 1998: 389-396.

[2] Resnyansky A. Directional fragmentation warheads:A theoretical and experimental investigation[C]∥18th International Symposium on Ballistics,San Antonio,TX, 1999: 543-550.

[3] 梁爭峰, 袁寶慧, 孫興昀, 等.定向戰(zhàn)斗部破片能量增益的數(shù)值模擬[J]. 火炸藥學(xué)報, 2005(2): 52-55.

LIANG Zheng-feng, YUAN Bao-hui, SUN Xing-yun, et al. Numerical simulation on the kinetic energy iincrease of aimed warhead[J].ChineseJournalofExplosives&Propellants, 2005(2): 52-55.

[4] 李翔宇,盧芳云.三種類型戰(zhàn)斗部破片飛散的數(shù)值模擬[J].火炸藥學(xué)報, 2007(1): 44-48.

LI Xiang-yu, LU Fang-yun. Numerical Simulation on fragments dispersion of three type warheads[J].ChineseJournalofExplosives&Propellants, 2007(1): 44-48.

[5] 葉小軍, 韓玉, 陳慶寶. 偏心起爆戰(zhàn)斗部速度增益的數(shù)值模擬及實驗[J]. 火炸藥學(xué)報, 2009(3): 29-34.

YE Xiao-jun, HAN Yu, CHEN Qing-bao. Numerical simulation and experiment of velocity gains on the non-central detonation warhead[J].ChineseJournalofExplosives&Propellants, 2009(3): 29-34.

[6] LI Yuan, WEN Yu-quan. Simulation on damage effectiveness of hexagonal prism aimable warhead with multi-point synchronous initiations[J].JournalofBeijingInstituteofTechnology, 2014(1): 1-7.

[7] Cunard D A Thomas K A. Programmable Integrated Ordnance Suite(PIOS)[C]∥1st Annual International Missiles & Rockets Exhibition, 1992

[8] Waggener Sam, Relative performance of anti-air missile warheads[C]∥19th International Symposium on Ballistics, Interlaken, Switzerland, 2001: 623-630.

[9] Borje Nyquist. Abraham to c-ram[R]. Technical report, 2007.

[10] Hallquist J O. LS-DYNA Keyword User′s Manual[M]. Livermore Software Technology Corporation, 2007, 970.

[11] 魏繼峰. 戰(zhàn)斗部多點起爆模擬研究[D]. 北京: 北京理工大學(xué), 2005.

WEI Ji-feng. Numerical research of multi-point initiated warhead[D]. Beijing: Beijing Institute of China, 2005.