張孝中, 趙太勇, 陳智剛, 郭光全, 李建飛

(1 中北大學(xué)地下目標(biāo)毀傷技術(shù)國(guó)防重點(diǎn)學(xué)科實(shí)驗(yàn)室, 太原 030051; 2 晉西工業(yè)集團(tuán)有限責(zé)任公司, 太原 030051)

0 引言

隨著現(xiàn)代科學(xué)技術(shù)的發(fā)展,尤其是軍事科技的發(fā)展,使得軍事目標(biāo)的抗打擊能力日益加強(qiáng)。應(yīng)用在末敏彈上的爆炸成型彈丸(EFP),可以用來攻擊坦克頂裝甲及艦船甲板。為了提高EFP遠(yuǎn)距離攻擊目標(biāo)的能力,可使用一種尾翼型EFP,使其具有良好的穩(wěn)定飛行能力。劉建青等[1]對(duì)三點(diǎn)起爆形成尾翼EFP進(jìn)行了實(shí)驗(yàn)研究及相關(guān)的數(shù)值模擬,路鳴等[2]對(duì)四點(diǎn)起爆形成帶尾翼EFP進(jìn)行了試驗(yàn)研究,于川等[3]采用了多點(diǎn)起爆方式,對(duì)帶尾翼翻轉(zhuǎn)爆炸成型彈丸進(jìn)行了實(shí)驗(yàn)研究,Bouet T H等[4]研究了多點(diǎn)起爆方式對(duì)尾翼EFP的影響,并得到具有良好尾翼結(jié)構(gòu)和穩(wěn)定飛行能力的尾翼型EFP。林加劍等人[5-6]使用在藥型罩上貼附惰性隔板的方法,對(duì)帶尾翼EFP形成進(jìn)行了相關(guān)實(shí)驗(yàn)研究。Weimann K等[7]就爆炸成型彈丸飛行穩(wěn)定性進(jìn)行了實(shí)驗(yàn)研究。唐平江等[8]主要研究了殼體對(duì)帶尾翼EFP成型過程的影響。門建兵等[9]進(jìn)行了帶尾翼EFP形成的三維數(shù)值模擬研究。而目前采用在藥型罩表面刻槽的方法形成尾翼型EFP的數(shù)值模擬研究較少。

文中提出一種在藥型罩表面周向?qū)ΨQ刻槽,并以45鋼填充槽體的結(jié)構(gòu),用該結(jié)構(gòu)來影響藥型罩的塑性變形過程,從而形成一種尾翼型的EFP。對(duì)比分析了3種藥型罩結(jié)構(gòu)方案,分別對(duì)3種方案進(jìn)行數(shù)值模擬,得到了最優(yōu)方案,繼而以該結(jié)構(gòu)方案對(duì)尾翼型EFP成型過程的影響進(jìn)行研究,得到刻槽參數(shù)及填充材料對(duì)尾翼型EFP侵徹性能和成型效果的影響規(guī)律,對(duì)今后尾翼型EFP研究工作具有指導(dǎo)意義。

1 計(jì)算模型與材料參數(shù)

1.1 計(jì)算模型

用LS-DYNA軟件[10]對(duì)刻槽藥型罩形成尾翼型EFP過程進(jìn)行數(shù)值模擬,采用Lagrange算法,單位制:cm-g-μs。因聚能裝藥為軸對(duì)稱結(jié)構(gòu),同時(shí)在藥型罩表面周向?qū)ΨQ刻4個(gè)槽,故建立1/4模型,以此為例進(jìn)行數(shù)值模擬。整個(gè)模型由炸藥、填充物和藥型罩三部分組成,模型如圖1所示,所有單元均為8節(jié)點(diǎn)實(shí)體單元solid164。計(jì)算中采用裝藥頂端中心點(diǎn)起爆方式。

模型中藥型罩的形狀為球缺罩,其主要結(jié)構(gòu)參數(shù)裝藥直徑為Dk=50 mm,藥型罩外曲率半徑為R=1.64Dk,壁厚為δ=0.06Dk,裝藥高度為L(zhǎng)=1.2Dk。

圖1 計(jì)算模型

1.2 材料參數(shù)選擇

在數(shù)值模擬中,藥型罩材料為紫銅,本構(gòu)模型為Johnson-Cook,狀態(tài)方程為Gruneisen,屈服應(yīng)力σy表示為:

-(T*)m]

(1)

式中:T*=(T-Tr)/(Tm-Tr),Tm為常態(tài)下材料熔化溫度,Tr為參考溫度(一般為室溫)。式中參數(shù)值為A=0.09 GPa,B=0.29 GPa,n=0.31,C=0.025,m=1.08,Tm=1 356 K,Tr=293 K。

選擇8701炸藥,本構(gòu)模型為HIGH_EXHIGH_EXPL OSIVE_BURN,狀態(tài)方程為JWL,表達(dá)式為:

(2)

式中:P為壓力;E為內(nèi)能。所用參數(shù)值為E=10 GPa,A=564 GPa,B=6.801 GPa,R1=4.1,R2=1.3,ω=0.36。填充物材料為45鋼。

藥型罩和填充物材料參數(shù)如表1所示,其中ρ為密度(g/cm3),G為剪切模量(GPa)。炸藥材料參數(shù)如表2所示,其中D為爆速(m/s),Pcj為爆轟波波陣面壓力(GPa)。

表1 藥型罩和填充物材料參數(shù)[11]

表2 炸藥材料參數(shù)[12]

2 數(shù)值模擬結(jié)果及分析

2.1 尾翼型EFP成型原理

由爆轟理論知,炸藥爆炸時(shí)將產(chǎn)生高溫、高壓的爆轟產(chǎn)物,爆轟時(shí)這些產(chǎn)物將沿炸藥表面的法線方向向外飛散。端面具有球缺型凹坑的圓柱狀藥柱爆炸后,爆轟產(chǎn)物一方面以一定速度沿凹坑表面的法線方向向中心軸線匯聚,形成一股高速度、高密度的氣流,即聚能效應(yīng);另一方面由于稀疏波的作用,匯聚到中心軸線處的爆轟產(chǎn)物又迅速向周圍低壓區(qū)膨脹,使能量分散開。在成型裝藥的球缺型凹坑表面放置一個(gè)與之相切的金屬罩,則爆炸后的爆轟產(chǎn)物將推動(dòng)罩壁向中心軸線運(yùn)動(dòng),將能量傳遞給金屬罩形成EFP。為了形成一種飛行穩(wěn)定性較好的尾翼型EFP,需在罩表面周向?qū)ΨQ刻槽,并填充材料達(dá)到閉氣作用。當(dāng)經(jīng)過衰減后的沖擊波作用于罩刻槽處時(shí),在未刻槽處的罩微元壓力比其他位置處大,產(chǎn)生的塑性變形也大,將形成尾翼型EFP。

由于45鋼的壓縮性比銅小,因此在受到相等爆轟載荷作用時(shí),沖擊波壓力在鋼中衰減比銅快,所以選用45鋼填充槽體。

文中提出3種藥型罩結(jié)構(gòu)方案:1)無刻槽;2)有刻槽,無填充物;3)有刻槽,有填充物。下文將對(duì)EFP成型過程進(jìn)行數(shù)值模擬,將得到較好方案。

2.2 3種方案對(duì)EFP形成尾翼的影響

對(duì)以上3種方案分別進(jìn)行數(shù)值模擬,取140 μs時(shí)結(jié)果進(jìn)行分析。EFP形態(tài)見表3,計(jì)算結(jié)果見表4。

表3 3種方案EFP形態(tài)

表4 3種方案EFP速度和長(zhǎng)徑比

從表3和表4中可知,第一種方案形成的EFP速度及長(zhǎng)徑比最小,且無尾翼;第二種方案形成的EFP速度及長(zhǎng)徑比次之,且尾翼形態(tài)較差;第三種方案形成的EFP速度及長(zhǎng)徑比最大,且尾翼張的很開,形態(tài)良好,成形效果較好。因此有刻槽、有填充物是較好方案,繼續(xù)對(duì)該方案分析。

當(dāng)刻槽數(shù)目、長(zhǎng)度及填充物材料一定時(shí),尾翼型EFP主要受槽寬和槽深等參數(shù)影響。下文將對(duì)尾翼型EFP成型過程進(jìn)行數(shù)值模擬,并得到兩種參數(shù)的最優(yōu)結(jié)構(gòu)設(shè)計(jì)。

2.3 刻槽深度和寬度對(duì)EFP形成尾翼的影響

當(dāng)刻槽長(zhǎng)度L=1.0 cm、寬度W=0.21 cm時(shí),僅改變槽深H,依次取值為0.1 cm、0.13 cm、0.16 cm、0.19 cm、0.22 cm、0.25 cm進(jìn)行模擬計(jì)算。取140 μs時(shí)結(jié)果進(jìn)行分析,其計(jì)算結(jié)果見表5。不同槽深的EFP形態(tài)見表6,EFP速度及長(zhǎng)徑比隨槽深的變化曲線如圖2所示。

表5 不同槽深和槽寬的EFP速度及其長(zhǎng)徑比

表6 不同槽深的EFP形態(tài)

由計(jì)算結(jié)果可知,槽深由0.1 cm增加到0.25 cm,EFP的速度及長(zhǎng)徑比都呈先增大后減小的趨勢(shì)。但槽深過大時(shí),尾翼成形效果較差。當(dāng)槽深為0.19 cm時(shí),EFP速度及長(zhǎng)徑比都達(dá)到最大值,尾翼形態(tài)較好。

當(dāng)槽長(zhǎng)L=1.0 cm、槽深H=0.19 cm時(shí),槽寬W依次取值為0.15 cm、0.18 cm、0.21 cm、0.24 cm、0.27 cm、0.3 cm進(jìn)行計(jì)算,其結(jié)果見表5。不同槽寬的EFP形態(tài)見表7,EFP速度及長(zhǎng)徑比隨槽寬的變化曲線如圖3所示。

圖2 EFP速度及長(zhǎng)徑比隨槽深的變化曲線

表7 不同槽寬的EFP形態(tài)

圖3 EFP速度及長(zhǎng)徑比隨槽寬的變化曲線

由模擬結(jié)果可知,槽寬由0.15 cm增加到0.30 cm,EFP的速度及長(zhǎng)徑比都呈現(xiàn)先增大后減小的趨勢(shì)。當(dāng)槽寬為0.24 cm時(shí),EFP長(zhǎng)徑比達(dá)到最大值,雖其速度有所下降,但下降不明顯,且其尾翼形態(tài)良好。

2.4 填充物材料對(duì)EFP形成尾翼的影響

當(dāng)刻槽參數(shù)不變時(shí),分析不同填充物材料對(duì)尾翼EFP的影響。當(dāng)槽深為0.19 cm、寬為0.24 cm時(shí),填充材料為鋁合金、陶瓷、45鋼,分別進(jìn)行數(shù)值模擬。填充材料參數(shù)見表8,不同填充材料的EFP形態(tài)見表9,不同填充材料的EFP速度、長(zhǎng)徑比和尾翼厚度見表10。

表8 填充材料參數(shù)

表9 不同填充材料的EFP形態(tài)

表10 不同填充物材料的EFP速度和長(zhǎng)徑比

模擬結(jié)果表明,隨著填充材料密度的增大,EFP速度及長(zhǎng)徑比都增大。尾翼厚度是影響彈丸飛行穩(wěn)定性重要的因素[13],填充材料為45鋼時(shí),形成EFP尾翼厚度較大,尾翼成型效果較好。

3 總結(jié)

1)藥型罩結(jié)構(gòu)對(duì)形成尾翼型EFP影響明顯,藥型罩有刻槽且槽內(nèi)有填充物時(shí),EFP尾翼張的很開,形態(tài)較好。

2)槽寬W和槽深H對(duì)尾翼型EFP影響顯著。當(dāng)槽寬為0.048倍、槽深為0.038倍裝藥直徑Dk時(shí),形成的尾翼型EFP速度及長(zhǎng)徑比最大,尾翼形態(tài)較好。此結(jié)論只適用于中小口徑聚能裝藥。

3)填充材料對(duì)尾翼型EFP也有影響,填充材料為45鋼時(shí),尾翼型EFP速度及長(zhǎng)徑比較大,尾翼成型效果較好。這對(duì)尾翼型EFP研究工作有一定的指導(dǎo)意義,數(shù)值模擬的結(jié)果還有待進(jìn)行試驗(yàn)驗(yàn)證。

參考文獻(xiàn):

[1] 劉建青, 郭濤, 顧文彬, 等. 三點(diǎn)起爆形成尾翼EFP的數(shù)值模擬和實(shí)驗(yàn)研究 [J]. 爆炸與沖擊, 2013, 33(1): 38-46.

[2] 陸鳴, 井志勝, 張端平, 等. 四點(diǎn)起爆形成帶尾翼EFP的試驗(yàn)研究 [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2009, 29(1): 128-129.

[3] 于川, 董慶東, 孫承緯, 等. 帶尾翼翻轉(zhuǎn)型爆炸成形彈丸試驗(yàn)研究 [J]. 爆炸與沖擊, 2003, 23(6): 561-564.

[4] BOUET T H, TARAYRE P, GUILLON J P. Study of a multi-point ignition EFP [C]∥National Defence Research Establishment. 15th International Symposium on Ballistics. Israel: [s. n. ]. 1995: 159-166.

[5] 林加劍, 沈兆武, 任輝啟, 等. 貼隔板法形成尾翼型EFP的試驗(yàn)研究及數(shù)值模擬 [J]. 火炸藥學(xué)報(bào), 2009, 32(1): 74-78.

[6] 趙慧英, 沈兆武, 李成兵, 等. 帶尾翼爆炸成型彈丸的新技術(shù) [J]. 含能材料, 2006, 14(2): 102-104.

[7] WEIMANN K. Flight stability of EFP with star shaped tail [C]∥National Defence Research Establishment. 14th International Symposium on Ballistics. Queber: [s. n. ]. 1993: 755-763.

[8] 唐平江, 顧文彬, 羅智偉. 變壁厚殼體對(duì)形成帶尾翼EFP的影響 [J]. 四川兵工學(xué)報(bào), 2009, 30(1): 34-36.

[9] 門建兵, 蔣建偉, 萬(wàn)麗珍. 帶尾翼EFP形成的三維數(shù)值模擬研究 [J]. 北京理工大學(xué)學(xué)報(bào), 2002, 22(2): 166-168.

[10] 李裕春, 時(shí)黨勇, 趙遠(yuǎn). LS-DYNA基礎(chǔ)理論與工程實(shí)踐 [M]. 北京: 中國(guó)水利水電出版社, 2008: 264-282.

[11] 張先鋒, 陳惠武, 趙有守. EFP沖擊引爆帶殼炸藥數(shù)值模擬研究 [J]. 彈道學(xué)報(bào), 2006, 18(1): 90-92.

[12] 張鈞, 陳智剛, 李小軍, 等. 變壁厚球缺罩桿式射流的形成與侵徹性能研究 [J]. 爆破器材, 2016, 45(1): 39-42.

[13] MALEJKO G, KUHNLE W G. Aerodynamic predictions for explosively formed penetrators with thick fins using an inviscid cfd code [C]∥AIAA. 32nd Aerospace Sciences Meeting & Exhibit. [S.l.:s.n.], 1994: 1-8.