王哲， 蔣建偉， 王樹有， 任芮池

(北京理工大學(xué) 爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室， 北京 100081)

雙層藥型罩形成的串聯(lián)爆炸成型彈丸速度計(jì)算模型

王哲， 蔣建偉， 王樹有， 任芮池

(北京理工大學(xué) 爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室， 北京 100081)

針對(duì)典型雙層藥型罩爆炸成型彈丸(EFP)裝藥結(jié)構(gòu)，基于微元爆轟驅(qū)動(dòng)及碰撞理論，考慮藥型罩軸向驅(qū)動(dòng)及徑向壓垮特性，建立了雙層藥型罩串聯(lián)EFP速度分析模型。計(jì)算結(jié)果表明：罩微元徑向壓垮引起的內(nèi)外罩EFP軸向速度差是串聯(lián)EFP分離的核心因素；外罩EFP速度隨外罩曲率及外罩所占厚度比增加而增加；內(nèi)罩EFP速度隨外罩曲率增加而減小，隨外罩所占厚度比增加而增加。結(jié)合數(shù)值模擬與脈沖X光成型試驗(yàn)對(duì)所推導(dǎo)理論進(jìn)行了驗(yàn)證，三者所得到的串聯(lián)EFP速度吻合較好。

兵器科學(xué)與技術(shù)； 爆轟； 串聯(lián)爆炸成型彈丸； 雙層藥型罩； 速度

0 引言

隨著各種目標(biāo)的防護(hù)能力不斷加強(qiáng)，新型裝甲的出現(xiàn)使傳統(tǒng)反裝甲聚能戰(zhàn)斗部受到嚴(yán)峻挑戰(zhàn)。研究人員嘗試應(yīng)用新的技術(shù)途徑和措施應(yīng)對(duì)這種挑戰(zhàn)，從而涌現(xiàn)出了若干新型反裝甲彈藥。串聯(lián)爆炸成型彈丸(EFP)就是其中之一。同軸放置的兩層藥型罩，通過合理設(shè)計(jì)其形狀和材料組合，在炸藥爆轟驅(qū)動(dòng)下能形成兩個(gè)同軸度很高的分離EFP或是單一大長徑比EFP，可應(yīng)付復(fù)合裝甲、爆炸反應(yīng)式裝甲等新型防護(hù)結(jié)構(gòu)的挑戰(zhàn)，提高對(duì)裝甲目標(biāo)的侵徹能力。國內(nèi)外已有不少學(xué)者開展了相關(guān)理論、數(shù)值模擬及試驗(yàn)研究[1-8]，鄭宇[9]基于瞬時(shí)爆轟理論，雖建立了串聯(lián)EFP速度分析模型，但該模型忽略徑向稀疏波作用，簡化內(nèi)外罩微元碰撞作用，且EFP與射流藥型罩微元壓垮過程不盡相同，尚不能很好地描述帶殼裝藥雙層藥型罩形成的串聯(lián)EFP速度。

針對(duì)典型雙層藥型罩裝藥結(jié)構(gòu)，考慮徑向稀疏波對(duì)裝藥軸向驅(qū)動(dòng)藥型罩微元的影響，基于藥型罩微元徑向壓垮及碰撞理論建立了雙層藥型罩形成的串聯(lián)EFP速度計(jì)算模型，并分析裝藥結(jié)構(gòu)參數(shù)對(duì)雙層EFP速度影響。同時(shí)采用數(shù)值模擬及對(duì)計(jì)算模型進(jìn)行驗(yàn)證，計(jì)算結(jié)果與數(shù)值模擬結(jié)果吻合較好。本文的研究結(jié)果對(duì)雙層藥型罩EFP設(shè)計(jì)具有一定參考意義。

1 雙層藥型罩EFP速度理論模型

1.1 雙層藥型罩EFP裝藥結(jié)構(gòu)

圖1 雙層藥型罩裝藥結(jié)構(gòu)示意圖Fig.1 Structure of double liner shaped charge

根據(jù)文獻(xiàn)[9]，雙層藥型罩受裝藥爆轟驅(qū)動(dòng)初期內(nèi)外兩層罩始終貼合在一起，當(dāng)作單層罩來考慮，等效單罩微元同時(shí)受裝藥軸向驅(qū)動(dòng)及徑向壓縮作用(如圖2所示)。為推導(dǎo)雙層藥型罩形成的EFP速度計(jì)算模型，作如下假設(shè)：

1)假設(shè)瞬時(shí)爆轟，不考慮起爆位置影響；

2)藥型罩微元受爆轟波驅(qū)動(dòng)，軸向速度未達(dá)到最大前內(nèi)外雙罩不分離，按等體積原則當(dāng)?shù)刃握痔幚恚?/p>

3)藥型罩微元受爆轟波驅(qū)動(dòng)，軸向速度達(dá)到最大后，內(nèi)外罩受裝藥爆轟徑向作用徑向壓垮導(dǎo)致雙層藥型罩微元軸向發(fā)生碰撞和動(dòng)量交換，導(dǎo)致內(nèi)外雙罩分離；

4)藥型罩材料為剛塑性模型，動(dòng)態(tài)屈服強(qiáng)度為定值；

5)EFP成形后徑向速度為0 m/s，軸向速度即為最終速度。

圖2 等效罩受裝藥軸向裝藥徑向裝藥作用示意圖Fig.2 Effect of explosive driving on liner element

1.2 等效單罩微元速度模型

圖2為等效單罩微元N受軸向、徑向裝藥驅(qū)動(dòng)示意圖。罩微元徑向位置為r，距罩頂距離為l.CNx為罩微元所對(duì)應(yīng)軸向裝藥質(zhì)量，CNy為罩微元所對(duì)應(yīng)徑向裝藥質(zhì)量，mN為罩微元質(zhì)量，mNkx為罩微元所對(duì)應(yīng)軸向殼體微元質(zhì)量，mNky為罩微元所對(duì)應(yīng)徑向殼體微元質(zhì)量。A面為罩底與裝藥軸線垂直面。

1.2.1 等效單罩微元軸向速度

經(jīng)炸藥爆轟驅(qū)動(dòng)平板Gurney公式[10]可得等效單罩微元N的軸向極限速度為

(1)

等效單罩微元質(zhì)量[7]為

(2)

式中：ρ為等效單罩密度；ε為r處雙層藥型罩厚度；α為罩微元N切線處與x軸夾角。

根據(jù)有效裝藥理論，圓柱裝藥對(duì)大錐角、球缺罩的爆轟驅(qū)動(dòng)可簡化為等效平板驅(qū)動(dòng)模型[11]?？紤]裝藥徑向膨脹引起的軸向能量輸出損失，藥型罩微元N的軸向速度[12]為

vNx=vNx0cos (0.766 8r/R)·

(1-0.376 1e(-0.720 7δ/(2R))(H/(2R))0.466 4).

(3)

1.2.2 等效單罩微元徑向速度

由炸藥驅(qū)動(dòng)圓筒向內(nèi)壓垮計(jì)算公式[13]可得藥型罩微元N的徑向極限速度為

(4)

式中：r0是炸藥內(nèi)虛擬剛壁(即炸藥爆轟產(chǎn)物速度為0 m/s處)半徑。

對(duì)于球缺罩，罩微元被壓垮到軸線或A面處時(shí)已不受炸藥徑向壓縮作用，無法被驅(qū)動(dòng)到徑向最大速度，等效罩微元徑向速度[10]為

(5)

式中：T為爆轟波到達(dá)微元的時(shí)間；t為罩微元運(yùn)動(dòng)到軸線處時(shí)間或運(yùn)動(dòng)到A面處時(shí)間的最小值；τ為時(shí)間常數(shù)。

1.3 雙層藥型罩內(nèi)外罩微元速度模型

根據(jù)假設(shè)3，爆轟波驅(qū)動(dòng)罩微元后期，內(nèi)外藥型罩微元發(fā)生碰撞和動(dòng)量交換，導(dǎo)致內(nèi)外罩分離。

1.3.1 內(nèi)外罩微元徑向速度模型

藥型罩微元爆轟加載下向內(nèi)壓垮，罩微元內(nèi)表面壓合速度高于外表面壓合速度。等效單罩微元N內(nèi)外表面徑向壓垮速度[14]為

vNyh=2vNy/(k+1),

vNyf=vNyh(1+1/((rNy/rN0)2-1))0.5,

(6)

式中：vNyh為等效單罩微元N靠近炸藥表面徑向壓垮速度；vNyf為等效單罩遠(yuǎn)離炸藥表面徑向壓垮速度；rN0為等效單罩N壓垮到軸線后微元半徑。

等效單罩微元N對(duì)應(yīng)內(nèi)外罩微元徑向速度為

vNyo=(((1-n)vNyh+nvNyf+vNyh))/2，

(7)

vNyi=((1-n)vNyf+nvNyh+vNyf)/2，

(8)

式中：vNyo為外罩微元徑向壓垮速度；vNyi為內(nèi)罩微元徑向壓垮速度；n為外罩微元與內(nèi)罩微元厚度比。

1.3.2 內(nèi)外罩微元軸向速度模型

定義藥型罩微元靠近裝藥一側(cè)為微元后端，遠(yuǎn)離裝藥一側(cè)為微元前端。根據(jù)能量守恒定律，雙層藥型罩微元徑向動(dòng)能除用于徑向壓縮塑性變形消耗外，其余全部轉(zhuǎn)化成為軸向的動(dòng)能。由于罩微元徑向壓垮帶來的額外軸向速度，外罩微元前端速度大于內(nèi)罩微元后端速度，內(nèi)外罩微元發(fā)生碰撞而產(chǎn)生動(dòng)量交換，從而導(dǎo)致內(nèi)外罩所形成的EFP由于軸向速度差而發(fā)生分離。罩微元由于徑向壓縮引起的軸向速度變化[9]為

(9)

式中：i=i，o分別代表內(nèi)外罩；ti為i罩微元徑向速度為0 m/s的時(shí)刻；σydi為i罩微元材料屈服強(qiáng)度；Vi為i罩微元體積；mNi為i罩微元質(zhì)量；vNxi為i罩微元軸向速度。

罩微元前后端軸向速度為

(10)

式中：vNxih、vNxif分別為i罩微元內(nèi)前后端軸向速度。

假定碰撞后外罩前端軸向速度與內(nèi)罩后端軸向速度最終一致為

(11)

外罩、內(nèi)罩微元最終軸向速度為

(12)

內(nèi)外罩EFP最終速度為

(13)

由上述理論推導(dǎo)可知，外罩EFP與內(nèi)罩EFP之間存在軸向速度差，不考慮速度衰減，常規(guī)等口徑雙層藥型罩串聯(lián)EFP飛行一定距離后會(huì)發(fā)生分離。

基于球缺藥型罩所建立的雙層EFP速度模型可擴(kuò)展應(yīng)用于大錐角藥型罩或球缺- 大錐角藥型罩組合雙層EFP速度計(jì)算。

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

為驗(yàn)證所推導(dǎo)雙層EFP速度計(jì)算模型準(zhǔn)確性，采用LS-DYNA軟件對(duì)圖1所示裝藥結(jié)構(gòu)進(jìn)行雙層EFP成型數(shù)值模擬。EFP結(jié)構(gòu)模型為：鋼殼厚δ=0.09R，裝藥長徑比H/R=2，銅制內(nèi)外罩頂總厚為0.1R，藥型罩壁厚差ξo/do=ξi/di=0.3. 裝藥選用JH-2炸藥，密度1.72 g/cm3.

數(shù)值模型中的殼體及炸藥均采用Lagrange算法。藥型罩采用自適應(yīng)網(wǎng)格方法，并適時(shí)使用重啟動(dòng)方法刪除炸藥和殼體網(wǎng)格。

圖3為計(jì)算模型所得等效單罩徑向速度隨外罩曲率變化示圖，圖4為計(jì)算所得等效罩微元徑向速度隨位置變化曲線圖。由圖3和圖4可知，等效罩微元徑向壓垮速度隨徑向位置增加先增加、后減小，隨罩曲率半徑增大而減小。這主要是因?yàn)殡S罩微元徑向位置及罩曲率半徑增加，罩微元到罩底端距離減小，受徑向裝藥驅(qū)動(dòng)時(shí)間減小，因而徑向壓垮速度發(fā)生如圖3和圖4所示變化。

圖3 等效單罩徑向速度隨外罩曲率變化圖Fig.3 Radial velocity of equivalent liner element vs. roh/2R

圖4 等效罩微元徑向速度隨位置變化圖Fig.4 Radial velocity of equivalent liner element vs. r/R

圖5為雙層EFP軸向速度隨曲率變化計(jì)算模型與數(shù)值模擬結(jié)果對(duì)比圖。由理論分析可知，相同徑向位置處罩微元徑向壓垮帶來的額外軸向速度隨外罩曲率增加而減小。故由圖5可知，隨外罩曲率增加，外罩EFP軸向速度增加，內(nèi)罩EFP軸向速度衰減，內(nèi)外罩EFP軸向速度差減小，有利于形成頭尾相連大長徑比的串聯(lián)EFP. 計(jì)算模型所得結(jié)果與數(shù)值模擬結(jié)構(gòu)吻合較好，最大誤差不超過5%.

圖5 雙層EFP軸向速度隨曲率變化圖Fig.5 Axial velocities of tandem EFP vs. roh/R

圖6為雙層EFP軸向速度隨外罩厚度變化圖(roh=2R)。由理論分析可知，內(nèi)外罩微元碰撞后外罩微元前端軸向速度與內(nèi)罩微元后端軸向速度隨外罩厚度增加而增加。故由圖6可知，內(nèi)外罩所形成EFP軸向速度均隨外罩所占厚度比增加而增加。計(jì)算模型所得結(jié)果與數(shù)值模擬結(jié)構(gòu)吻合較好，最大誤差不超過6%.

圖6 雙層EFP軸向速度隨外罩厚度變化圖Fig.6 Velocities of tandem EFP vs. do/(do+di)

3 試驗(yàn)驗(yàn)證

為驗(yàn)證上述理論計(jì)算模型及數(shù)值模擬結(jié)果，采用圖1所示裝藥結(jié)構(gòu)開展雙層藥型罩EFP成型脈沖X光驗(yàn)證試驗(yàn)。表1為試驗(yàn)工況表。

表1 脈沖X光雙層藥型罩EFP成型試驗(yàn)方案表

表2和表3為起爆后不同時(shí)刻雙層EFP脈沖X光試驗(yàn)與數(shù)值模擬結(jié)果的對(duì)比表。表4和表5為脈沖X光雙層藥型罩EFP軸向速度對(duì)比表。雙層EFP脈沖X光成型結(jié)果與數(shù)值模擬結(jié)果吻合較好，計(jì)算模型、數(shù)值模擬及試驗(yàn)所得雙層EFP內(nèi)外罩速度基本一致，表明本文所采用的數(shù)值模擬方法及其所驗(yàn)證理論模型是可靠的。

表2 雙層EFP成型X光試驗(yàn)與數(shù)值模擬結(jié)果對(duì)比表(工況1)

Tab.2 Comparison of X-ray test and simulated double-layer liner EFP formation results in Condition 1

表3 雙層EFP成型X光試驗(yàn)與數(shù)值模擬結(jié)果對(duì)比表(工況2)

表4 雙層藥型外罩EFP速度對(duì)比表

工況外罩速度/(m·s-1)試驗(yàn)值仿真值理論計(jì)算值11755169216712218920722038

表5 雙層藥型內(nèi)罩EFP速度對(duì)比表

4 結(jié)論

基于藥型罩微元爆轟驅(qū)動(dòng)及碰撞理論建立了雙層藥型罩串聯(lián)EFP速度計(jì)算模型。通過計(jì)算及分析得到以下結(jié)論：

1)雙層藥型罩串聯(lián)EFP速度計(jì)算模型與數(shù)值模擬及試驗(yàn)結(jié)果吻合較好，分析模型與數(shù)值模擬結(jié)果最大誤差不超過6%，能夠較好描述裝藥結(jié)構(gòu)參數(shù)對(duì)串聯(lián)EFP速度影響規(guī)律。

2)藥型罩微元徑向壓垮引起的內(nèi)外罩EFP軸向速度差是串聯(lián)EFP分離的核心因素。徑向壓垮速度隨藥型罩曲率增加而減小，常規(guī)等口徑雙層藥型罩結(jié)構(gòu)內(nèi)罩EFP與外罩EFP之間存在軸向速度差。

3)隨外罩曲率增加，罩微元徑向速度減小，外罩EFP軸向速度增加，內(nèi)罩EFP軸向速度衰減，有利于形成頭尾相連大長徑比的串聯(lián)EFP。

4)裝藥長徑比、殼厚、曲率等因素一定的條件下, 隨外罩所占厚度比增加，內(nèi)罩EFP、外罩EFP軸向速度增加。

References)

[1] Weimann K, Blache A. An explosively formed projectile with tantalum penetration and steel stabilization base[C]∥Proceedings of the 18th International Symposium on Ballistics. San Antonio, TX, US: IBC,1999:22-25.

[2] Fong R, Weimann K. Testing and analysis of multi-liner EFP warheads[C]∥Proceedings of the 20th International Symposium on Ballistics. Orlando, FL, US:IBC, 2002:578-582.

[3] 鄭宇, 王曉鳴, 李文彬. 基于雙層藥型罩成型裝藥的串聯(lián)EFP[J]. 爆炸與沖擊, 2012, 32(1):29-33. ZHENG Yu, WANG Xiao-ming, LI Wen-bin. Formation of shaped charge with double layer liners into tandem EFP[J]. Shock and Waves, 2012, 32(1):29-33.(in Chinese)

[4] 鄭宇, 王曉鳴, 李文彬, 等. 曲率半徑對(duì)雙層球缺罩形成串聯(lián)爆炸成型彈丸的影響[J]. 高壓物理學(xué)報(bào), 2009, 23(3):33-38. ZHENG Yu, WANG Xiao-ming, LI Wen-bin, et al. Effects of liner curvature radius on formation of double-layered spherical segment charge liner into tandem explosively formed projectile (EFP)[J]. Chinese Journal of High Pressure Physics, 2009, 23(3):33-38. (in Chinese)

[5] 鄭宇, 王曉鳴, 李文彬. 材料對(duì)雙層藥型罩形成串聯(lián)EFP的影響[J]. 兵器材料科學(xué)與工程, 2009, 32(1):38-41. ZHENG Yu, WANG Xiao-ming, LI Wen-bin. Effects of materials on formation of double-layer shaped charges into tandem EFP[J]. Ordnance Material Science and Engineering, 2009, 32(1):38-41. (in Chinese)

[6] 劉健峰, 龍?jiān)? 紀(jì)沖, 等. 不同藥型罩壁厚組合影響同軸EFP成型規(guī)律研究[J]. 工程爆破, 2012, 18(4):9-13. LIU Jian-feng, LONG Yuan, JI Chong, et al. Influence of different wall-thickness combinations of arc-cone liner on coaxial EFP forming law[J]. Engineering Blasting, 2012, 18(4):9-13. (in Chinese)

[7] 劉健峰, 毛振興, 龍?jiān)? 等. 曲率半徑對(duì)雙層藥型罩EFP戰(zhàn)斗部成型及侵徹的影響[J]. 爆破器材, 2016,45(3):5-10. LIU Jian-feng. MAO Zhen-xing, LONG Yuan, et al. Influence of curvature radius on the formation and penetration of explosively formed penetrator warhead of the double layer liner[J]. Explosive Materials, 2016, 45(3):5-10. (in Chinese)

[8] 劉健峰, 龍?jiān)? 紀(jì)沖, 等. 雙層藥型罩EFP戰(zhàn)斗部性能參數(shù)的灰色系統(tǒng)理論分析及實(shí)驗(yàn)研究[J]. 含能材料, 2016, 24(8): 728-734. LIU Jian-feng, LONG Yuan, JI Chong, et al. Theoretical and experimental study on performance parameters of double layer liners EFP warhead based on grey system theory[J]. Chinese Journal of Energetic Materials, 2016, 24(8): 728-734. (in Chinese)

[9] 鄭宇. 雙層藥型罩毀傷元形成機(jī)理研究[D]. 南京:南京理工大學(xué), 2008. ZHENG Yu. Study on the formation mechanism of kill element from shaped charge with double layer liners[D]. Nanjing:Nanjing University of Science and Technology, 2008. (in Chinese)

[10] 隋樹元, 王樹山, 終點(diǎn)效應(yīng)學(xué)[M]. 北京: 國防工業(yè)出版社, 2000. SUI Shu-yuan, WANG Shu-shan. Terminal effects[M]. Beijing: National Defense Industry Press, 2000. (in Chinese)

[11] 林加劍. EFP成型及其終點(diǎn)效應(yīng)研究[D]. 合肥:中國科學(xué)技術(shù)大學(xué), 2009. LIN Jia-jian. Research on the formability of EFP and its terminal effect[D].Hefei: University of Science and Technology of China, 2009.(in Chinese)

[12] Jiang J W, Feng S S, Men J B. Simplified approach to velocity of axial forward projection fragments from the end of warhead[J]. Journal of Beijing Institute of Technology, 2004, 13(4):402-405.

[13] Chou P C, Carleone J, Hirsch E, et al. Improved formulas for velocity, acceleration and projection angle of explosively driven liners[C]∥Proceedings of the 6th International Symposium on Ballistics. Orlando, FL, US: IBC, 1981.

[14] 夏杰.雙層藥型罩射流形成的理論研究與數(shù)值模擬[D]. 武漢:武漢科技大學(xué), 2012. XIA Jie. The theoretical research and numerical simulation on jet formation of double-layer liner[D]. Wuhan:Wuhan University Science and Technology, 2012. (in Chinese)

A Calculation Model of Velocity of Tandem EFP with Double Layer Liners

WANG Zhe, JIANG Jian-wei, WANG Shu-you, REN Rei-chi

(State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China)

A velocity calculation model of tandem EFP is proposed by analyzing the motion of elements in double-layer liners. The analysis shows that the difference in EFP axial velocities caused by radial velocity plays an important role in the separation of the tandem EFP. The velocity of EFP formed by outside liner increases with the increase in the curvature of outside liner and the thickness ratio of outside liner to inside liner. On the other hand, the velocity of EFP formed by inside liner decreases with the increase in curvature of outside liner, but increases with the increase in thickness ratio of outside liner to inside liner. The calculated results agree well with experimental data and simulated results by LS-DYNA.

ordnance science and technology; explosion; tandem EFP; double layer liners; velocity

2016-11-11

爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室自主課題探索性項(xiàng)目(QNKT16-01)；國家國防科技工業(yè)局技術(shù)基礎(chǔ)科研項(xiàng)目(B2620110006)

王哲(1989—)，男，博士研究生。E-mail:994353449@qq.com

蔣建偉(1962—)，男，教授，博士生導(dǎo)師。E-mail:bitjjw@bit.edu.cn

TJ410.3+33

A

1000-1093(2017)07-1301-06

10.3969/j.issn.1000-1093.2017.07.007