羅少敏,徐誠(chéng),陳愛軍,張曉云,劉蘇蘇

(1.南京理工大學(xué)機(jī)械工程學(xué)院,江蘇南京 210094;2.南京理工大學(xué)理學(xué)院,江蘇南京 210094)

步槍彈侵徹帶軟硬復(fù)合防護(hù)明膠靶標(biāo)的數(shù)值模擬

羅少敏1,徐誠(chéng)1,陳愛軍2,張曉云1,劉蘇蘇1

(1.南京理工大學(xué)機(jī)械工程學(xué)院,江蘇南京 210094;2.南京理工大學(xué)理學(xué)院,江蘇南京 210094)

為研究步槍彈撞擊帶軟硬復(fù)合防護(hù)明膠靶標(biāo)的作用過程和作用機(jī)理,采用顯式有限元方法對(duì)7.62 mm步槍彈侵徹復(fù)合靶標(biāo)過程進(jìn)行數(shù)值模擬,分析侵徹過程中的典型現(xiàn)象及明膠靶標(biāo)動(dòng)態(tài)響應(yīng)。數(shù)值計(jì)算結(jié)果表明:陶瓷錐的形成是由壓縮應(yīng)力波和拉伸應(yīng)力波共同作用的結(jié)果;彈頭加速度變化存在明顯的分段與拐點(diǎn),侵徹陶瓷面板過程中,加速度達(dá)到最大,侵徹聚乙烯(PE)背板層時(shí),出現(xiàn)第二個(gè)拐點(diǎn);由于防護(hù)層存在多個(gè)界面,撞擊過程中PE背板界面存在速度多峰現(xiàn)象:當(dāng)彈頭運(yùn)動(dòng)加速度達(dá)到最大時(shí),PE背板界面出現(xiàn)第一個(gè)速度峰,明膠界面出現(xiàn)第一個(gè)壓力峰;當(dāng)彈頭開始侵徹PE背板時(shí),背板層出現(xiàn)第二個(gè)速度峰;在步槍彈撞擊過程中明膠內(nèi)壓力波傳遞呈現(xiàn)球形波基本形態(tài),壓力峰值隨距離增加呈指數(shù)衰減。

兵器科學(xué)與技術(shù);槍彈;數(shù)值模擬;軟硬復(fù)合防護(hù);明膠;侵徹

0 引言

現(xiàn)代戰(zhàn)場(chǎng)槍彈與破片對(duì)士兵的威脅最大,由硬質(zhì)防彈插板和軟質(zhì)防彈衣組成的軟硬復(fù)合防護(hù)成為最主要的單兵防護(hù)具。一般情況下,步槍彈的速度達(dá)到700～980 m/s,在該速度段內(nèi)針對(duì)步槍彈對(duì)帶軟硬復(fù)合防護(hù)的有生目標(biāo)作用效應(yīng)的研究具有重要的現(xiàn)實(shí)意義和工程應(yīng)用價(jià)值。

關(guān)于步槍彈撞擊軟硬復(fù)合防護(hù)的研究,國(guó)內(nèi)外主要集中在步槍彈直接侵徹陶瓷及侵徹陶瓷面板與背板組成的復(fù)合防護(hù)層,研究侵徹全過程動(dòng)態(tài)響應(yīng)和防護(hù)層防護(hù)能力。Krishnan等[1]研究了陶瓷/高聚乙烯(UHMWPE)復(fù)合板的抗彈性能,以數(shù)值模擬對(duì)比彈道穿深實(shí)驗(yàn)驗(yàn)證了有限元網(wǎng)格的精細(xì)程度對(duì)計(jì)算結(jié)果的影響,分析了彈頭對(duì)陶瓷/高聚乙烯復(fù)合板的侵徹過程,為單兵防彈衣的毀傷預(yù)測(cè)提供了參考。Feli等[2-3]研究了槍彈侵徹陶瓷—纖維復(fù)合材料過程,認(rèn)為隨著彈頭初始速度的減小陶瓷錐最大角度增大,同時(shí)隨著初始速度的增加,纖維凹面層數(shù)增多,頂部復(fù)合板的變形減小。Fawaz[4]等基于Hypermesh和LS-DYNA平臺(tái)研究了陶瓷復(fù)合防護(hù)的侵徹響應(yīng),認(rèn)為在有角度侵徹比垂直侵徹,破壞更為嚴(yán)重,為防護(hù)材料的優(yōu)化提供了合理的依據(jù)。但上述研究未直接涉及到步槍彈對(duì)軟硬復(fù)合防護(hù)后的有生目標(biāo)作用效應(yīng)。由于有生目標(biāo)的特殊性和復(fù)雜性,對(duì)其作用效應(yīng)的研究難以直接進(jìn)行,國(guó)內(nèi)外通常采用明膠靶標(biāo)來(lái)模擬。帶軟硬復(fù)合防護(hù)的明膠靶標(biāo)是彈道實(shí)驗(yàn)和評(píng)估步槍彈對(duì)帶防護(hù)有生目標(biāo)的作用效應(yīng)的主要靶標(biāo),目前國(guó)內(nèi)外未見采用數(shù)值模擬研究步槍彈對(duì)帶軟硬復(fù)合防護(hù)的明膠靶標(biāo)作用效應(yīng)相關(guān)文獻(xiàn)報(bào)道。

本文結(jié)合相關(guān)彈道實(shí)驗(yàn),建立了與之對(duì)應(yīng)的有限元數(shù)值計(jì)算模型,利用顯式有限元方法進(jìn)行數(shù)值模擬,再現(xiàn)了7.62 mm步槍彈對(duì)帶軟硬復(fù)合防護(hù)的明膠靶標(biāo)的侵徹過程。數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好,驗(yàn)證了所建立的有限元模型的正確性。在此基礎(chǔ)上,研究侵徹過程中典型物理現(xiàn)象和非貫穿侵徹時(shí)防護(hù)后明膠靶標(biāo)動(dòng)態(tài)響應(yīng),為非貫穿侵徹帶軟硬復(fù)合防護(hù)有生目標(biāo)的毀傷機(jī)理認(rèn)識(shí)、評(píng)估,槍彈優(yōu)化設(shè)計(jì)以及單兵防護(hù)裝備優(yōu)化設(shè)計(jì)提供了參考。

1 數(shù)值計(jì)算模型及實(shí)驗(yàn)驗(yàn)證

1.1 有限元模型

本文研究的帶軟硬復(fù)合防護(hù)明膠靶標(biāo)幾何結(jié)構(gòu)如圖1所示,它由陶瓷面板、聚乙烯(PE)背板和高聚乙烯纖維軟防護(hù)層及明膠組成。陶瓷輕質(zhì)高強(qiáng)防彈插板為PE材料陶瓷復(fù)合板,由高性能PE纖維“UD”無(wú)緯布和氧化鋁陶瓷材料構(gòu)成,符合美國(guó)NIJ標(biāo)準(zhǔn)防Ⅲ級(jí),規(guī)格尺寸為250 mm×300 mm.軟質(zhì)防護(hù)為超高分子量聚乙烯纖維防彈材料,符合GA141-2001警用標(biāo)準(zhǔn)防Ⅱ級(jí)。明膠材料為4℃、10%彈道明膠,作為仿生肌肉靶標(biāo),在制作過程中,將插板和軟防護(hù)與明膠材料一起脫模。本文主要研究7.62 mm步槍彈對(duì)該帶軟硬復(fù)合防護(hù)明膠靶標(biāo)的作用效應(yīng)。

對(duì)圖1所示靶標(biāo)進(jìn)行有限元網(wǎng)格劃分,離散化后有限元網(wǎng)格如圖2(a)所示,為兼顧計(jì)算時(shí)長(zhǎng)與計(jì)算精度,對(duì)防護(hù)層和明膠塊中心彈著點(diǎn)區(qū)域進(jìn)行進(jìn)一步細(xì)化:面板、背板層和軟防護(hù)在彈著點(diǎn)中心5倍半徑區(qū)域加密,向外劃分逐漸稀疏,面板和背板中心區(qū)域劃分30段,如圖2(b)所示;彈頭由彈頭殼被甲、鉛套和鋼芯組成,彈頭殼錐部與圓柱部過渡部分由3個(gè)單元轉(zhuǎn)化為2個(gè)單元,鉛套前部精細(xì)劃分,鋼芯部分由頭部向尾部逐漸稀疏,如圖2(c)所示。彈頭、陶瓷面板層、PE背板層和明膠塊采用SOLID單元,其中彈頭中彈頭殼共 9 792個(gè)單元,鉛套共10 400個(gè)單元,鋼芯共9 216個(gè)單元;陶瓷面板共63 000個(gè)單元;PE背板共126 000個(gè)單元;明膠塊共180 000個(gè)單元。軟防護(hù)共46層,采用Shell單元建模,中心加密方式與面板和PE背板一致,共165 600個(gè)單元。

數(shù)值模擬中邊界條件設(shè)置與實(shí)驗(yàn)一致,彈頭殼、鉛套與鋼芯設(shè)置自動(dòng)面面侵蝕接觸,靶板與軟質(zhì)防護(hù)以及明膠塊之間設(shè)置為自動(dòng)面面接觸,陶瓷插板面板層與背板層之間設(shè)置固連失效接觸,實(shí)驗(yàn)中軟防護(hù)層與層之間無(wú)膠粘材料,因此軟質(zhì)防護(hù)內(nèi)部設(shè)置自動(dòng)單面接觸,彈頭與復(fù)合靶標(biāo)之間設(shè)置面面侵蝕接觸,整體采用單點(diǎn)積分,明膠塊設(shè)置沙漏控制。

圖2 彈頭侵徹靶板及明膠有限元網(wǎng)格Fig.2 Finite element mesh of bullet impacting armor and gelatin

1.2 材料模型與參數(shù)

彈頭的被甲、鉛套和鋼芯采用Johnson-Cook模型結(jié)合Gruneisen狀態(tài)方程來(lái)描述彈頭在侵徹過程中材料的力學(xué)行為,其材料參數(shù)與文獻(xiàn)[5]相同。陶瓷面板采用Johnson-Hoimquist模型來(lái)描述,根據(jù)文獻(xiàn)[6-7]來(lái)進(jìn)行主要參數(shù)的賦值。描述背板和軟防護(hù)的材料性能都采用 MAT_COMPOSITE_ FAILURE_OPTION_MODEL材料模型,對(duì)于PE背板采用SOLID單元,對(duì)于軟防護(hù)采用SHELL單元。該模型由Cheng等[8]嵌入LS-DYNA軟件中,是包含8種失效模式的正交異型材料模型,根據(jù)所做基礎(chǔ)材料實(shí)驗(yàn)及文獻(xiàn)[9]確定主要材料參數(shù)值。表1給出了PE背板采用的材料模型參數(shù),表2給出了PE纖維軟防護(hù)的材料模型參數(shù)。明膠靶標(biāo)采用流體彈塑性材料模型(MAT_ELASTIC_PLASTIC_HYDRO)結(jié)合LINEAR_POLYNOMIAL狀態(tài)方程來(lái)描述,材料參數(shù)參考文獻(xiàn)[10-11],與文獻(xiàn)[5,12]相同。

1.3 模型驗(yàn)證

為驗(yàn)證本文所建立有限元仿真模型的可信性和準(zhǔn)確性,進(jìn)行了7.62 mm步槍彈侵徹帶軟硬復(fù)合防護(hù)的明膠靶標(biāo)實(shí)驗(yàn)。實(shí)驗(yàn)中采用的軟硬復(fù)合防護(hù)靶標(biāo)結(jié)構(gòu)尺寸與圖1相同,入靶速度由光電測(cè)速靶測(cè)量,明膠內(nèi)壓力由壓電晶體傳感器測(cè)量,高速攝像機(jī)拍攝彈頭與防護(hù)及靶標(biāo)作用過程。下面從實(shí)驗(yàn)后彈頭結(jié)構(gòu)形態(tài)和穿深、背凸量、最大壓力、明膠靶標(biāo)瞬時(shí)凹入量等方面將計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行了比較。

表1 PE背板材料主要參數(shù)Tab.1 Mechanics properties of materials of PE

表2 軟防護(hù)材料主要參數(shù)Tab.2 Mechanics properties of materials of so ft armor

數(shù)值仿真結(jié)果和實(shí)驗(yàn)結(jié)果皆顯示步槍彈均未穿透防護(hù):步槍彈穿過陶瓷面板,繼續(xù)侵徹PE背板層,最終停留在插板背板層中。如圖3所示,彈頭作為投射物在撞擊過程中受到防彈插板阻礙,不斷破碎及變形,彈芯頭部墩粗,呈明顯的“蘑菇”形,計(jì)算的彈頭形態(tài)與實(shí)驗(yàn)結(jié)果基本一致。

表3給出了計(jì)算獲得的典型侵徹特征量與實(shí)驗(yàn)結(jié)果的比較,可知計(jì)算獲得的步槍彈穿深、剩余長(zhǎng)度、明膠中壓力峰值以及明膠靶標(biāo)最大瞬時(shí)凹入量等與實(shí)驗(yàn)結(jié)果基本一致。綜合計(jì)算與實(shí)驗(yàn)對(duì)比情況可知,本文所建立的有限元模型具有較好的可信度和準(zhǔn)確性,能夠有效地模擬殺傷元對(duì)軟硬復(fù)合防護(hù)后的明膠靶標(biāo)的侵徹作用和動(dòng)態(tài)響應(yīng)。

圖3 彈頭形態(tài)對(duì)比Fig.3 Comparison of bullet configurations

2 計(jì)算結(jié)果分析

2.1 侵徹過程中防護(hù)層應(yīng)力變化特性

彈體撞擊陶瓷面板,產(chǎn)生一束壓力波,壓力波在橫向和縱向方向傳播,彈著點(diǎn)及其附近區(qū)域開始產(chǎn)生損傷,如圖4中10 μs時(shí)刻。縱向壓力波沿彈道方向傳播,到達(dá)陶瓷面板與PE背板界面處時(shí),由于陶瓷材料和高聚乙烯材料的密度和聲阻抗等性質(zhì)的差異,會(huì)反射形成一個(gè)拉伸波,被壓縮的固體材料會(huì)膨脹,拉伸波與彈體前進(jìn)方向反向傳播,不同材料界面開始出現(xiàn)損傷,當(dāng)時(shí)長(zhǎng)在30 μs時(shí),陶瓷錐基本形成。圖4描述了典型時(shí)刻段防護(hù)材料應(yīng)力分布情況。反射拉伸波在陶瓷面板內(nèi)不斷傳播,陶瓷錐不斷擴(kuò)展,至彈頭停止運(yùn)動(dòng)。

在由硬質(zhì)插板和軟質(zhì)防護(hù)材料組成的防護(hù)體系中,插板防護(hù)對(duì)彈頭的侵徹阻滯和動(dòng)能消耗起主要作用。彈頭在界面擊潰過程中材料流動(dòng)發(fā)生質(zhì)量侵蝕和變形,損耗部分動(dòng)能,同時(shí)受到靶板阻力作用,侵徹速度下降,直至被攔截停止運(yùn)動(dòng)。彈頭與軟硬復(fù)合防護(hù)接觸全過程,伴隨了陶瓷擊潰、彈頭變形、背板背凸變形、軟質(zhì)防護(hù)變形和明膠靶標(biāo)動(dòng)態(tài)凹陷等主要物理現(xiàn)象。在侵徹過程中,陶瓷材料受到拉壓應(yīng)力波的反復(fù)作用發(fā)生裂紋擴(kuò)展、破碎及飛濺,如圖5所示。通過靶板z向應(yīng)變?cè)茍D與等值線可以看出,拉伸損傷和壓縮破壞對(duì)陶瓷錐的影響很大,陶瓷錐中心區(qū)域由于受到壓縮應(yīng)力波作用出現(xiàn)負(fù)應(yīng)變,在陶瓷錐邊沿區(qū)域由于受到反射拉伸波的作用出現(xiàn)正應(yīng)變,應(yīng)變?cè)茍D與陶瓷錐基本形態(tài)相似。背板層與彈頭直接接觸時(shí),阻擋彈頭繼續(xù)侵徹,發(fā)生變形,壓縮軟質(zhì)防護(hù)層變形,同時(shí)明膠靶標(biāo)受到撞擊產(chǎn)生凹陷,彈頭直接侵徹結(jié)束后明膠靶標(biāo)發(fā)生長(zhǎng)時(shí)間動(dòng)態(tài)卸載消耗傳遞至模擬靶標(biāo)的能量。

表3 數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比Tab.3 Comparison of numerical and experimental results

2.2 彈靶相互作用過程中速度、加速度變化特性

彈頭速度變化、加速度的變化與整個(gè)侵徹過程中陶瓷面板、PE背板的物理特性相關(guān),與陶瓷錐的成型、背板的變形直接相關(guān)。圖6給出了計(jì)算所獲得的彈頭侵徹帶軟硬復(fù)合防護(hù)的明膠靶標(biāo)過程中速度曲線和加速度曲線,在彈頭速度衰減曲線上標(biāo)出了3個(gè)特征點(diǎn):A0點(diǎn)為運(yùn)動(dòng)初始接觸點(diǎn),B0點(diǎn)為彈頭部分即將穿過陶瓷面板與背板接觸時(shí)刻,C0點(diǎn)為運(yùn)動(dòng)速度為0時(shí)刻。

在加速度曲線上,A1點(diǎn)為運(yùn)動(dòng)初始點(diǎn),B1點(diǎn)為加速度最大點(diǎn),在C1點(diǎn)時(shí)加速度曲線出現(xiàn)第二個(gè)拐點(diǎn),D1點(diǎn)加速度值減小為0.根據(jù)圖6和圖4,將彈頭侵徹過程和加速度變化分為3段:

1)A1B1段,約為0～20 μs,彈頭頭部與陶瓷面板初始接觸,加速度急劇增大,彈頭頭部材料向外流動(dòng)而發(fā)生質(zhì)量侵蝕和變形。在彈頭界面擊潰過程中,陶瓷面板受到高速?zèng)_擊作用,開始出現(xiàn)初始裂紋,在頭部“擠進(jìn)”陶瓷面板后,加速度達(dá)到最大,此階段視為彈頭侵徹過程的開坑階段。開坑部分直徑大于彈頭直徑。

2)B1C1段,約為20～50 μs,陶瓷面板由于短時(shí)間的強(qiáng)沖擊破碎,對(duì)彈頭的阻力減小,后期主要以摩擦阻力為主。在此階段內(nèi)陶瓷錐形成,彈頭與陶瓷錐一起運(yùn)動(dòng),同時(shí)受到背板阻力與面板摩擦力作用,加速度減小。此階段可視為彈頭在陶瓷面板內(nèi)運(yùn)動(dòng)的穩(wěn)定侵徹階段,速度曲線上的B0點(diǎn)與B1C1曲線段上的C1點(diǎn)很接近。

3)C1D1段,約為50～90 μs,在C1點(diǎn),加速度曲線出現(xiàn)第二個(gè)拐點(diǎn),彈頭頭部部分嵌入PE層壓板內(nèi)部,彈頭主要受到背板阻力作用,加速度逐漸衰減,直至減小為0.速度曲線上的C0點(diǎn)與C1D1段上的D1點(diǎn)不重合。隨著彈頭的運(yùn)動(dòng)PE背板的變形成為主要的吸能方式,背板背凸量增大直至彈頭停留在背板內(nèi),同時(shí)彈芯發(fā)生墩粗,整體長(zhǎng)度變短。

彈頭在侵徹過程中主要包括:界面擊潰、彈芯墩粗變形、磨損發(fā)熱等。彈頭變形主要發(fā)生在約前50 μs:在彈頭界面擊潰階段,彈頭頭部錐頭部分磨損侵蝕,彈頭形狀和質(zhì)量發(fā)生改變,與靶板接觸面積變大;在彈頭受到背板防護(hù)阻力作用過程中,彈頭的變形對(duì)侵徹的影響較為明顯,伴隨彈芯墩粗和質(zhì)量侵蝕,侵徹速度下降,同時(shí)會(huì)發(fā)生加速度的跳躍,如B1C1曲線段較為明顯的跳躍點(diǎn)。

圖4 侵徹過程典型時(shí)刻Fig.4 Penetration processes at some typical moments

圖5 靶板z向應(yīng)變?cè)茍D(t=30 μs)Fig.5 Strain contour of ceramic/PE composite target in direction z(t=30 μs)

圖6 彈頭速度和加速度變化曲線(v0=734 m/s)Fig.6 The variation of projectile velocity and acceleration(v0=734 m/s)

2.3 靶標(biāo)界面速度變化特性

在整個(gè)侵徹過程中靶標(biāo)存在3個(gè)界面,其中兩個(gè)明顯界面:一是插板防護(hù)與軟質(zhì)防護(hù)之間;二是軟質(zhì)防護(hù)與明膠靶標(biāo)之間。一個(gè)不明顯界面,即插板防護(hù)中陶瓷面板與PE背板材料之間。分析界面位置處單元及質(zhì)點(diǎn)的速度、加速度對(duì)于研究非貫穿侵徹作用下明膠靶標(biāo)的損傷機(jī)制尤為重要。

數(shù)值模擬中在彈頭彈道方向上選取PE背板與軟質(zhì)防護(hù)、軟質(zhì)防護(hù)與明膠靶標(biāo)界面處單元,分析明膠靶標(biāo)在侵徹過程中的動(dòng)態(tài)特性。與軟防護(hù)接觸的PE背板單元為71 986,與軟防護(hù)接觸的明膠靶標(biāo)單元為354 809,圖7為背板及明膠界面速度變化曲線,可以看出明膠界面單元的最大速度為66.8 m/s, PE背板單元運(yùn)動(dòng)的最大速度為87.3 m/s.

如圖7所示,背板界面存在速度多峰。對(duì)比圖6,背板界面單元第一個(gè)速度峰值點(diǎn)出現(xiàn)的時(shí)間約在20 μs左右,與彈頭加速度最大值出現(xiàn)的時(shí)刻基本一致。背板速度第一個(gè)峰值的出現(xiàn)與彈頭的強(qiáng)沖擊直接相關(guān):彈頭撞擊陶瓷面板層,在此過程中出現(xiàn)第一個(gè)速度上升沿并達(dá)到最大;在形成陶瓷錐后,彈頭與陶瓷錐一起運(yùn)動(dòng),繼續(xù)運(yùn)動(dòng)過程中伴隨陶瓷的破碎,對(duì)彈頭的阻力開始減小,出現(xiàn)“卸載”現(xiàn)象,但隨著彈頭到達(dá)PE背板,與PE背板接觸,PE背板的速度增大出現(xiàn)第二個(gè)峰值。軟防護(hù)最后一層與明膠界面單元第一個(gè)速度峰值時(shí)刻出現(xiàn)一定延遲,這與沖擊過程中壓力波的傳遞有關(guān),同時(shí)由于防護(hù)層變形,軟防護(hù)最后一層會(huì)擠壓明膠界面使得明膠界面速度出現(xiàn)第二個(gè)平緩峰值,隨后明膠界面與軟防護(hù)最后一層開始出現(xiàn)明顯的界面分離。

圖7 界面單元速度變化Fig.7 The variation of interface velocity

2.4 明膠內(nèi)壓力波傳遞特性

在撞擊過程中,當(dāng)彈頭接觸防護(hù)材料,壓力波形成,傳遞至界面,形成拉伸波。壓力波繼續(xù)經(jīng)由軟質(zhì)防護(hù)傳播至明膠模擬靶標(biāo),如圖8所示,在傳播初始階段,明膠靶標(biāo)內(nèi)典型時(shí)刻壓力波傳遞呈現(xiàn)球形波基本形態(tài)。沿彈頭侵徹方向每隔2 cm選取明膠上單元,共4個(gè)測(cè)點(diǎn),測(cè)點(diǎn)壓力隨時(shí)間變化曲線如圖9所示。

計(jì)算獲得的明膠界面壓力峰值達(dá)到94.1 MPa,界面壓力曲線會(huì)出現(xiàn)第二個(gè)壓力峰值,這是由于彈頭與PE背板撞擊所致;壓力波在明膠內(nèi)傳播速度為1 421.5 m/s,實(shí)驗(yàn)中測(cè)得傳播速度為1 501 m/s,二者基本一致。壓力峰值隨距離衰減情況如圖10所示,壓力波在傳播過程中,隨著傳播距離增加,壓力峰值呈指數(shù)下降。

3 結(jié)論

本文針對(duì)步槍彈侵徹帶軟硬復(fù)合防護(hù)的明膠靶標(biāo),進(jìn)行了有限元建模和數(shù)值模擬,通過與實(shí)驗(yàn)結(jié)果的對(duì)比驗(yàn)證了有限元模型的可信性和準(zhǔn)確性,計(jì)算結(jié)果表明:

1)彈頭在撞擊過程中,彈頭加速度存在明顯的分段與拐點(diǎn)。侵徹陶瓷面板過程中,加速度達(dá)到最大,侵徹PE背板層時(shí),出現(xiàn)第二個(gè)拐點(diǎn)。彈頭變形主要發(fā)生在侵徹陶瓷面板運(yùn)動(dòng)過程中,侵徹背板時(shí)主要吸能方式為PE背板的拉伸變形。

圖8 明膠內(nèi)壓力波傳遞Fig.8 Propagation of shock waves in gelatin

圖9 不同距離測(cè)點(diǎn)處壓力波形曲線Fig.9 Pressure of shock waves at different distances

2)由于靶標(biāo)存在多個(gè)防護(hù)層界面,撞擊過程中PE背板界面點(diǎn)存在速度多峰現(xiàn)象,出現(xiàn)時(shí)刻與彈頭運(yùn)動(dòng)存在一定的關(guān)聯(lián)。當(dāng)彈頭運(yùn)動(dòng)加速度達(dá)到最大時(shí),在彈頭強(qiáng)沖擊下,背板材料背面彈道方向單元速度達(dá)到第一個(gè)速度峰,隨后明膠界面出現(xiàn)第一個(gè)壓力峰,最大壓力達(dá)到94.1 MPa,防護(hù)層界面處明膠單元最大運(yùn)動(dòng)速度達(dá)到66.8 m/s.當(dāng)彈頭開始進(jìn)入背板層時(shí),PE背板第二個(gè)速度峰,導(dǎo)致明膠壓力曲線出現(xiàn)第二個(gè)峰值。

3)沖擊過程中,明膠內(nèi)壓力波傳遞呈球形波形態(tài),壓力峰值隨傳播距離呈指數(shù)規(guī)律衰減。

圖10 壓力峰值隨距離變化情況Fig.10 Peak pressure versus distance

References)

[1] Krishnan K,Sockalingam S,Bansal S,et al.Numerical simulation of ceramic composite armor subjected to ballistic impact[J]. Composites,2010,41:583-593.

[2] Feli S,Asgari M R.Finite element simulation of ceramic/composite armor under ballistic impact[J].Composite Structures,2011, 42(4):771-780.

[3] Feli S,Asgari M R.An analytical model for perforation of ceramic/multi-laryered planar woven fabric targets byblunt projectiles [J].Composite Structures,2011,93(2):548-560.

[4] Fawaz Z,Zheng W,Behdinan K.Numerical simulation of normal and oblique ballistic impact on ceramic[J].Composite Structures, 2004,63(2):387-395.

[5] 溫垚珂,徐誠(chéng),陳愛軍,等.步槍彈侵徹明膠靶標(biāo)的數(shù)值模擬[J]兵工學(xué)報(bào),2013,34(1):14—19.

WEN Yao-ke,XU Cheng,CHEN Ai-jun,et al.Numerical simulation of penetration of bullet on gelatin target[J].Acta Armamentarii,2013,34(1):14-19.(in Chinese)

[6] Simha C H M.Computational modeling of the penetration response of a high-purity ceramic[J].International Journal of Impact Engineering,2002,27(1):65-86.

[7] Cronin D S,Bui K,Kaufmann C,etal.Implementation and validation of the Johnson-Holmquist ceramic material model in LS-Dyna[C]//4th European LS-DYNA Users Conference.Ulm,Germany:DYNA more Gmh,2003:D-I-47-D-I-59.

[8] Cheng W,Hallquist J.Implementation of three-dimensional composite failure model into DYNA3D[M].US:Livermore Software Technology Corporation,2004.

[9] Menna C,Asprone D,Caprino G,et al.Numerical simulation of impact tests on GFRP composite laminates[J].International Journal of Impact Engineering,2011,38(8/9):677-685.

[10] Segletes S B.Modeling the penetration behavior ofrigid spheres into ballistic gelatin,ARL-TR-4393[R].Aberdeen Proving Ground,MD:Army Research Laboratory,2008.

[11] Kwon J,Subhash G.Compressive strain rate sensitivity of ballistic gelatin[J].Journal of Biomechanics,2010,43(3):420-425.

[12] Wen Y K,Xu C,Wang H S,et al.Impact of steel spheres on ballistic gelatin at moderate velocities[J].International Journal of Impact Engineering,2013,62:142-151.

Numerical Simulation of Bullets Penetrating into Gelatin Target with Hard/Soft Composite Armor

LUO Shao-min1,XU Cheng1,CHEN Ai-jun2,ZHANG Xiao-yun1,LIU Su-su1
(1.School of Mechanical Engineering,Nanjing University of Science and Technology,Nanjing 210094,Jiangsu,China; 2.School of Science,Nanjing University of Science and Technology,Nanjing 210094,Jiangsu,China)

In order to explore the interaction process and mechanism of bullet impacting a gelatin target with hard/soft composite armor,the penetration of 7.62 mm bullet into composite armor and gelatin is numerically simulated using an FEA method,and the typical phenomena of impacting process and the dynamic response of gelatin target are analyzed.The simulation results show that the formation of the ceramic cone results from the compressive and tensile stress waves.The obvious segmentation and inflection points display on the curve of projectile acceleration.The acceleration reaches its maximum when the projectile penetrates into a ceramic faceplate.The second inflection point appears when the projectile penetrates into PE rear plate.The multi-peaks appear on the velocity curve of PE rear plate due to multiple interfaces on protection layer:When the projectile reaches the maximum acceleration,the first peak appears,and the first pressure peak appears on gelatin interface;the second peak appears when the projectile penetrates into PE rear plate.The propagation of pressure wave in gelatin presents the basic form of spherical wave,and the pressure peak propagation complies with the law of exponential decay.

ordnance science and technology;bullet;numerical simulation;hard/soft composite armor; gelatin;penetration

TJ012.4

:A

:1000-1093(2014)08-1172-07

10.3969/j.issn.1000-1093.2014.08.006

2013-10-25

羅少敏(1987—),男,博士研究生。E-mail:1119lsm@163.com;徐誠(chéng)(1962—),男,教授,博士生導(dǎo)師。E-mail:xucheng62@mail.njust.edu.cn