曹蘇雅拉圖，王雨時，王 強(qiáng) ，聞 泉 ，王乃耀

(1南京理工大學(xué)機(jī)械工程學(xué)院，南京 2 10094;2九江45信箱，江西九江 3 32008)

0 引言

引信慣性觸發(fā)裝置設(shè)計通常要兼顧鈍感度和靈敏度要求。文中對頭部不同形狀的彈丸以不同著角碰擊薄目標(biāo)過程進(jìn)行數(shù)值仿真，得到彈丸碰擊薄目標(biāo)時的最大前沖過載系數(shù)值，從而給引信慣性觸發(fā)裝置設(shè)計提供工作環(huán)境。

眾多文獻(xiàn)對頭部不同形狀的彈丸侵徹過載特性的研究主要側(cè)重于混凝土和鋼板等硬目標(biāo)，Jan Arild Teland等[1]采用LS-DYNA程序研究了截錐形彈丸侵徹混凝土目標(biāo)的過載特性;聶明飛等[2]利用LS-DYNA軟件，對彈丸侵徹多層混凝土靶板進(jìn)行了數(shù)值仿真，討論了混凝土靶板的破壞現(xiàn)象，研究了卵形頭部彈丸對多層混凝土目標(biāo)的侵徹規(guī)律以及目標(biāo)厚度對彈丸前沖過載特性的影響;徐文崢等[3]采用LS-DYNA軟件仿真了尖頭形彈丸正侵徹混凝土過程，推導(dǎo)出具有不同幾何特點(diǎn)的尖頭彈丸侵徹過載數(shù)學(xué)模型，并得出了錐形和卵形頭部彈丸侵徹混凝土靶板的沖擊過載特性。孫煒海等[4]對固支薄金屬靶板的低速穿透過程進(jìn)行了評估，給出了錐形頭部彈丸穿透金屬靶板的耗能公式和彈道極限公式。

1 數(shù)值模型

1.1 有限元模型

以某非旋轉(zhuǎn)彈丸質(zhì)量4.25 kg，碰靶速度為140 m/s為例。另假設(shè)彈體外殼厚度1 mm，材料7A04鋁合金。根據(jù)文獻(xiàn)[5]所做的研究，彈體外殼內(nèi)按土壤處理。為了保證簡化后彈丸重量一致，彈長略有不同，但這并不影響仿真結(jié)果[5]。彈丸頭部幾何形狀選取為圓柱形、截錐形、半球形和圓錐形，不考慮尾翼結(jié)構(gòu)。靶板為鋁合金板、水曲柳膠合板和杉木板(厚3 mm)。整個建模過程采用g·cm·μs單位制。本模型中彈體與靶板之間采用面－面侵蝕接觸算法。為節(jié)省計算量，建立四分之一模型。不同彈頭幾何形狀的彈丸如圖1所示。

圖1 不同頭部形狀的彈丸

1.2 材料模型及其狀態(tài)方程

根據(jù)文獻(xiàn)[6]，采用 JOHNSON-COOK材料模型GRUNEISEN狀態(tài)方程定義彈體外殼材料的模型和狀態(tài)方程。JOHNSON-COOK材料模型的本構(gòu)方程描述了應(yīng)變、應(yīng)變速率和溫度對流動應(yīng)力的影響，其表達(dá)式為[7]:

式中:σ為流動應(yīng)變;ε為流動等效塑性應(yīng)變;ε*m為無量綱塑性應(yīng)變率;A、B、n、C、m 為材料參數(shù);T*=( T －Tr)/( Tm－Tr)為無量綱溫度;Tm為熔點(diǎn);Tr為室溫。

GRUNEISEN狀態(tài)方程定義了體積與壓力之間的關(guān)系。其表達(dá)式為[7]:

式中:p為材料壓力;μ = ρ/ρ0－1;C為 νs－ νp曲線截距;S1、S2、S3為 νs－ νp曲線斜率系數(shù);γ0為GRUNEISEN常數(shù);α為γ0的一階體積修正;E為材料內(nèi)能。

用“*MAT_PLASTIC_KINEMATIC”定義彈體填充物。用“*MAT_PLASTI－C_KINEMATIC”定義靶板材料模型。塑性隨動模型通過在0和1間調(diào)整硬化參數(shù)β來選擇隨動硬化或各項(xiàng)同性，其應(yīng)變率用COWPER-SYMONDS 模型表示，屈服應(yīng)力為[8]:

式中:σ0為初始屈服應(yīng)力;C、P為 C OWPERSYMON DS應(yīng)變率參數(shù);εepff為有效塑性應(yīng)變;β為硬化參數(shù);ε˙為應(yīng)變率;Ep為塑性硬化模量。由參考文獻(xiàn)[9－11]，得7A04鋁合金、土壤、水曲柳膠合板和杉木板材料模型的仿真參數(shù)如表1和表2所列。按照表1和表2設(shè)置彈體和靶板材料模型和狀態(tài)方程參數(shù)。

表1 7A04鋁合金Johnson-Cook材料模型參數(shù)

表2 土壤、水曲柳膠合板和杉木板材料模型的主要參數(shù)

2 仿真結(jié)果

圖2 圓柱形彈丸以不同著角碰擊三種薄靶板的最大前沖過載系數(shù)

2.1 圓柱形彈丸

對圓柱形彈丸以不同著角碰擊三種薄目標(biāo)進(jìn)行了仿真。以不同著角碰擊鋁合金板、水曲柳膠合板和杉木板的最大前沖過載系數(shù)如圖 2所示。碰擊鋁合金板的最大前沖過載系數(shù)值比碰擊水曲柳膠合板和杉木板大很多，圓柱形彈丸以小著角和大著角碰擊薄靶板的最大前沖過載系數(shù)值相差較大，著角變化對最大前沖過載系數(shù)值影響很大，碰擊三種薄靶板時的最大前沖過載系數(shù)值均隨著角的增大而單調(diào)遞減，與著角呈反比關(guān)系。彈丸以大于10°著角碰擊木質(zhì)板時著角變化對前沖過載系數(shù)影響不大。

圓柱形彈丸以不同著角碰擊水曲柳膠合板的前沖過載系數(shù)曲線如圖3所示。

圓柱形頭部彈丸四種不同著角碰擊水曲柳膠合板時的最大前沖過載值從大到小依次為:10°、20°、30°、45°、60°和 75°。著角為10°時前沖過載曲線未見波動，隨著角的增大前沖過載系數(shù)曲線越來越抖動，著角較大時其它三種頭部形狀的前沖過載系數(shù)曲線也有抖動。

圖3 圓柱形彈丸以不同著角碰擊水曲柳膠合板的前沖過載系數(shù)曲線

2.2 半球形頭部彈丸

對半球形頭部彈丸以不同著角碰擊三種薄目標(biāo)進(jìn)行仿真。以不同著角碰擊鋁合金板、水曲柳膠合板和杉木板的最大前沖過載系數(shù)如圖4所示。半球形頭部彈丸以0°～10°著角碰擊三種薄靶板時的最大前沖過載系數(shù)值變化較小，在此范圍內(nèi)著角對彈丸前沖過載系數(shù)影響較小;當(dāng)以大于10°著角碰擊三種薄靶板時彈丸的最大前沖過載系數(shù)隨著角的增大而減小，與著角呈反比關(guān)系。

圖4 半球形頭部彈丸以不同著角碰擊三種不同目標(biāo)時的最大前沖過載系數(shù)

2.3 截錐形頭部彈丸

對截錐形頭部彈丸以不同著角碰擊三種薄目標(biāo)進(jìn)行仿真。以不同著角碰擊鋁合金板、水曲柳膠合板和杉木板的最大前沖過載系數(shù)如圖5所示。

圖5 截錐形頭部彈丸以不同著角碰擊三種靶板的最大前沖過載系數(shù)

彈丸碰擊鋁合金板時著角變化對最大前沖過載系數(shù)的影響較小，最大前沖過載系數(shù)值變化范圍在250內(nèi)。碰擊水曲柳膠合板和杉木板時彈丸最大前沖過載系數(shù)隨著角的增大而先增大到最大值后逐漸減小，在著角75°左右時均出現(xiàn)峰值。

2.4 圓錐形頭部彈丸

對圓錐形頭部彈丸以不同著角碰擊三種薄目標(biāo)進(jìn)行仿真。以不同著角碰擊鋁合金板、水曲柳膠合板和杉木板的最大前沖過載系數(shù)如圖 6所示。頭部為圓錐形的彈丸以0°～10°著角碰擊水曲柳膠合板和杉木板時最大前沖過載系數(shù)值變化較小，此范圍內(nèi)著角變化對最大前沖過載系數(shù)影響不大;以大于10°著角碰擊時彈丸最大前沖過載系數(shù)隨著角的增大而增大，到峰值后逐漸變小，未呈現(xiàn)線性變化。彈丸碰擊鋁合金板時著角變化對最大前沖過載系數(shù)影響較小，最大前沖過載的變化范圍在200之內(nèi)，并且著角與彈丸最大前沖過載系數(shù)未呈現(xiàn)明顯規(guī)律，著角小于10°時最大前沖過載系數(shù)未呈線性變化，在某一特定角度處達(dá)到最大值。

圖6 圓錐形頭部彈丸以不同著角碰擊三種靶板的最大前沖過載系數(shù)

2.5 最大前沖過載系數(shù)仿真結(jié)果匯總

圓柱形、截錐形、半球形和圓錐形四種不同彈頭幾何形狀的彈丸以三種著速和0°～85°著角碰擊不同材質(zhì)薄目標(biāo)的最大前沖過載系數(shù)如表3所列，表中數(shù)值單位為g。由表3可看出，四種頭部形狀的彈丸以大于75°著角碰擊鋁合金板時均出現(xiàn)跳彈現(xiàn)象。圓柱形頭部彈丸以0°、1°、2°和 3°著角碰擊水曲柳膠合板時的最大前沖系數(shù)值為2190、1784、883和738;截錐形頭部彈丸以0°、1°、2°和 3°著角碰擊水曲柳膠合板時的最大前沖系數(shù)值為357、408、352和226，以上兩種情況均是著角對最大前沖過載系數(shù)影響大的例證。由圓柱形頭部彈丸碰擊鋁合金板可知，彈丸著角變化對最大前沖過載系數(shù)的影響很大，因此布靶引起的著角誤差對彈丸前沖過載系數(shù)影響很大，且無規(guī)律。

四種不同頭部形狀的彈丸以10°著角碰擊鋁合金板的前沖過載系數(shù)在初始階段有較大的波動，曲線呈現(xiàn)多個峰值，不像碰擊水曲柳膠合板和杉木板一樣過載系數(shù)曲線呈線性，四種不同頭部形狀彈丸之間的前沖過載系數(shù)峰值相差較大。

此外，按文獻(xiàn)[12]提供的彈丸穿透膠合板時的作用力公式計算截錐形彈丸以0°著角碰擊水曲柳膠合板的理論前沖過載系數(shù)為349.3，仿真結(jié)果為356.9，仿真結(jié)果與理論估算值相差不到3%，仿真結(jié)果與理論估算結(jié)果相符合，說明仿真結(jié)果是較為可信的。

表3 最大前沖過載系數(shù)

3 結(jié)論

圓柱形、截錐形、半球形和圓錐形頭部形狀彈丸小著角(小于10°)碰擊3 mm厚水曲柳膠合板的前沖過載系數(shù)依次分別為283～2190、103～408、210～229、11～24;碰擊3 mm厚杉木板的前沖過載系數(shù)依次分別為246～1407、99～249、113～125、10～14;碰擊3 mm厚鋁板的前沖過載系數(shù)依次分別為3917～6810、1854～1913、2707 ～2818、1011 ～1112。圓柱形頭部和截錐形頭部彈丸的前沖過載系數(shù)受著角變化影響較大，因此圓柱形和截錐形頭部彈丸應(yīng)考慮射擊試驗(yàn)時布靶誤差導(dǎo)致的著角變化以及彈丸攻角變化對引信靈敏度和鈍感度試驗(yàn)結(jié)論的影響。

四種頭部不同形狀的彈丸大攻角碰靶時前沖過載系數(shù)曲線均呈現(xiàn)較大的抖動，出現(xiàn)多次峰值，故有必要關(guān)注這種過載抖動對引信動態(tài)特性的影響。

截錐形頭部彈丸以不同著角碰擊木質(zhì)靶板時最大前沖過載曲線出現(xiàn)二次峰值，著角和最大前沖過載系數(shù)不呈線性關(guān)系。廣泛應(yīng)用的截錐形頭部彈射擊試驗(yàn)時應(yīng)考慮其前沖過載值不是在正碰時最大，應(yīng)關(guān)注小著角的微小變化對引信鈍感度和觸發(fā)靈敏度試驗(yàn)結(jié)果的影響。

半球形頭部彈丸以小著角(小于10°)碰擊鋁合金板的最大前沖過載系數(shù)比碰擊同厚度水曲柳膠合板大10倍左右，比碰擊杉木板大20倍左右。半球形頭部彈丸進(jìn)行射擊試驗(yàn)時，選用鋁合金靶板相比木質(zhì)靶板更能有效提高引信觸發(fā)靈敏度。

圓錐形頭部彈的最大前沖過載值比其它三種頭部彈丸小很多，說明圓錐形頭部彈相比其它三種頭部彈丸碰擊靶板時所受到的阻力較小。在強(qiáng)度條件得到滿足的情況下，設(shè)計鈍感度要求較高的引信時彈頭或彈頭引信形狀選為圓錐形，可有效提高引信鈍感度。

[1]Jan Arild Teland Henrik SjФl.Penetration into concrete by truncated projectiles[J].International Journal of Impact Engineering，2004，30(3):447－464.

[2]聶明飛，李玉龍.卵形頭部彈侵徹單多層混凝土靶板有限元仿真[J].探測與控制學(xué)報，2009，31(4):78－82.

[3]徐文崢，楊榕.彈藥幾何因素對過載特性影響的數(shù)值分析[J].彈箭與制導(dǎo)學(xué)報，2010，30(2):97－100.

[4]孫煒海，文鶴鳴.錐頭彈丸低速撞擊下薄金屬靶板的穿透[J].固體力學(xué)學(xué)報，2009，30(4):361－366.

[5]曹蘇雅拉圖，王雨時，王強(qiáng)，等.某尾翼火箭彈碰擊薄弱目標(biāo)前沖過載系數(shù)仿真研究[J].探測與控制學(xué)報，2012(2).

[6]LSTC，LS-DYNA KEYWORD USER’S MANUAL(VOLUME I)[M].California:Livermore Software Technology Corporation，2007.

[7]陳春黽，李伯龍，聶怍仁.穿甲彈侵徹鋁合金靶板數(shù)值模擬研究[J].科學(xué)技術(shù)與工程，2010(14):3442－3445.

[8]尚曉江，蘇建宇，王化鋒.ANSYS/LS-DYNA動力分析方法與工程實(shí)例[M].北京:中國水利水電出版社，2008

[9]王雨時，張德智，曹純柱.觸發(fā)引信防雨試驗(yàn)的等效靶板模擬[J].沈陽工業(yè)學(xué)院學(xué)報，1994(3):27－37.

[10]機(jī)械工程手冊、電機(jī)工程手冊編輯委員會.機(jī)械工程手冊(工程材料卷)[M].2版.北京:機(jī)械工業(yè)出版社，1996.

[11]樓劍鋒，王政，洪滔，等.鎢合金桿侵徹半無限厚鋁合金靶的數(shù)值研究[J].高壓物理學(xué)報，2009，23(1):65－70.

[12]高世橋，王寶興.截錐形彈低中速貫穿薄靶板時的動力分析與計算[J].應(yīng)用數(shù)學(xué)和力學(xué)，1986，7(11):1033－1038.