鎢桿超高速侵徹混凝土靶侵徹深度研究

王 鵬,付建平,崔 晉,袁 浩,昝守東,石曉山,楊 麗

(1.中北大學(xué) 機(jī)電工程學(xué)院, 太原 030051;2.中北大學(xué) 智能武器研究院, 太原 030051;3.國營803廠, 西安 710073 ; 4.國營524廠,吉林 132000;5.晉西工業(yè)集團(tuán)有限責(zé)任公司, 太原 030000)

0 引言

混凝土作為一種常見的復(fù)合材料被廣泛應(yīng)用于民用和國防工程,尤其在軍事防護(hù)工程中,混凝土建筑會(huì)受到各種彈體的毀傷,動(dòng)能侵徹是其中一種。動(dòng)能彈對混凝土的侵徹效應(yīng)一直是各國研究重點(diǎn)。桿式穿甲彈作為典型的穿甲武器被各國軍隊(duì)所裝備,隨著發(fā)射平臺技術(shù)的發(fā)展,給予桿式彈的速度逐漸提高。由于彈靶材料特性的影響,鎢桿在超高速與常規(guī)速度的情況下侵徹混凝土,毀傷效應(yīng)有所不同,毀傷機(jī)制也發(fā)生相應(yīng)的變化,因此,研究超高速情況下鎢桿侵徹混凝土是十分有必要的的。盧正操等[1]基于動(dòng)態(tài)空穴膨脹模型和A-T模型,建立不同侵徹模式下計(jì)算各種參數(shù)的理論公式,對長桿彈侵徹半無限混凝土靶進(jìn)行了系統(tǒng)的理論研究,結(jié)果表明:理論模型預(yù)測侵徹模式、侵徹深度、時(shí)程曲線和殘余質(zhì)量時(shí)與現(xiàn)有混凝土侵徹試驗(yàn)數(shù)據(jù)吻合得較好;周剛等[2]研究了鎢合金彈體超高速侵徹混凝土靶的相關(guān)理論,開展了鎢合金彈體超高速撞擊混凝土靶試驗(yàn),分析了成坑特性,研究了侵徹總深和殘余彈體長度隨撞擊速度的變化規(guī)律,結(jié)果表明:超高速撞擊條件下混凝土靶成坑為“彈坑+彈洞”形,成坑體積與彈體動(dòng)能近似成正比,侵徹深度隨彈速的提高呈現(xiàn)先增大后減小的現(xiàn)象,超高速侵深降低是彈體動(dòng)能經(jīng)歷銷蝕侵徹后“剛體侵徹階段”減少造成的;高飛等[3]研究砂漿靶目標(biāo)在動(dòng)能彈超高速撞擊下的破壞響應(yīng),利用2級輕氣炮開展卵形頭部鋼桿彈以1 200～2 400 m/s速度侵徹砂漿靶試驗(yàn),分析得到靶體開坑直徑和開坑深度與撞擊速度呈線性關(guān)系。侵深隨撞擊速度的增加呈現(xiàn)先線性增加、后逆減,在緩慢增加的趨勢,分別對應(yīng)剛體侵徹、半破碎侵徹和破碎侵徹。王明洋[4-5]運(yùn)用內(nèi)摩擦理論,結(jié)合巖石動(dòng)態(tài)可壓縮特性,建立了超高速動(dòng)能彈對地毀傷的侵徹深、成坑和安全厚度等設(shè)計(jì)計(jì)算方法。錢秉文等[6-8]研究彈體超高速撞擊混凝土靶,得到了侵深隨撞擊速度的增加呈先增加后減小的趨勢,在某一速度侵深存在最大值。程怡豪等[9]根據(jù)試驗(yàn)結(jié)果在4個(gè)沖擊速度區(qū)間內(nèi)圍繞混凝土靶的破壞特點(diǎn)及若干要素對靶體行為的影響規(guī)律展開論述,既低速區(qū)間內(nèi)(撞擊速度不超過40 m/s)簡要介紹梁、柱、板殼的破壞模式,中速區(qū)間內(nèi)(撞擊速度為40～1 000 m/s)論述彈靶相對尺寸、粗骨料、單軸抗壓強(qiáng)度和鋼筋對混凝土厚靶侵徹效應(yīng)的影響,高速區(qū)間內(nèi)(撞擊速度為1.0～2.0 km/s)描述半流體轉(zhuǎn)變速度和成坑效應(yīng),超高速區(qū)間內(nèi)(撞擊速度超過2 km/s)描述密度效應(yīng)占主導(dǎo)的流體侵徹行為。

針對上訴研究現(xiàn)狀,本文中開展0.2～3.0 km/s超高速鎢桿垂直侵徹混凝土研究,通過試驗(yàn)驗(yàn)證數(shù)值模型的合理性,通過LS-DYNA數(shù)值模擬結(jié)果分析超高速侵徹下無量綱侵深、彈體速度、質(zhì)量、殘余長度以及侵徹階段分析。

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

1.1 試驗(yàn)概況

圖1為試驗(yàn)現(xiàn)場布置示意圖,試驗(yàn)以125 mm火炮作為發(fā)射平臺,炮口距靶板125 m,通過調(diào)節(jié)發(fā)射藥量來控制彈體的著靶速度,前5 m放置章動(dòng)靶紙,靶用來觀察彈體著靶前飛行穩(wěn)定性。靶前41 m處放置天幕測速儀,靶板為兩塊圓柱形混凝土靶,前靶為靶Ⅰ,后靶為靶Ⅱ,靶Ⅰ和靶Ⅱ接觸放置。為了盡量避免邊界效應(yīng)的影響,對靶板進(jìn)行鋼箍加固。

圖1 試驗(yàn)現(xiàn)場布置示意圖

1.2 數(shù)值模擬及參數(shù)

仿真模型的彈體直徑24 mm,長徑比為11,質(zhì)量約為1.8 kg,彈體頭部為尖錐型,材料為鎢合金,彈體用MAT_JOHNSON_COOK本構(gòu)模型和EOS_GRUEISEN狀態(tài)方程描述,具體參數(shù)見表1[10]。其中Johnson-cook本構(gòu)關(guān)系為:

表1 鎢合金材料參數(shù)

Gruneisen狀態(tài)方程為:

式中:p為壓縮材料的壓力;E為材料的內(nèi)能;C0為Hugoniot關(guān)系式的截距;S1、S2、S3為Hugoniot關(guān)系式的系數(shù); ρ0為初始密度; γ0為Gruneisen系數(shù);a是對γ0的一維修正; μ=ρ/ρ0-1。

混凝土靶板用HJC本構(gòu)模型[11],該模型是一種專門針對混凝土受沖擊載荷作用而開發(fā)的動(dòng)態(tài)材料模型,模型考慮了大應(yīng)變、高應(yīng)變率和高壓情況,同時(shí)結(jié)合損傷,理論考慮了材料的拉伸斷裂行為,考慮了材料壓潰后的體積壓縮量與壓力的函數(shù)關(guān)系,并采用MAT_ADD_EROSION失效準(zhǔn)則,主應(yīng)變失效取0.4,剪應(yīng)變失效取0.2[12]。靶板為2塊?1.5m×1m混凝土靶,具體參數(shù)見表2。其中A、B、C為混凝土材料參數(shù),G為剪切模量, fc為單軸抗壓強(qiáng)度,彈體與靶板均采用Lagrange算法的solid164單元。彈體頭部網(wǎng)格尺寸0.15mm,靶板中心網(wǎng)格尺寸5mm,遠(yuǎn)離侵徹區(qū)域網(wǎng)格15mm。

表2 混凝土材料參數(shù)

由于靶板及彈丸具有對稱性,為節(jié)省計(jì)算時(shí)間,建立1/2模型,彈靶之間設(shè)置侵蝕接觸,在對稱面施加對稱邊界約束。有限元模型如圖2所示。

圖2 彈靶1/2模型

1.3 試驗(yàn)驗(yàn)證及數(shù)值模擬方法

圖3為章動(dòng)靶紙,圖4為彈體不同初速侵徹靶板的試驗(yàn)結(jié)果和模擬仿真結(jié)果。

圖3 章動(dòng)靶紙

圖4 彈體不同初速侵徹靶板的試驗(yàn)結(jié)果和模擬仿真結(jié)果

根據(jù)章動(dòng)靶紙(圖3)痕跡可知,彈體著靶前飛行姿態(tài)水平,著角可以忽略,視為彈體垂直侵徹混凝土靶板。

由圖4可知,由于彈體撞擊和沖擊波破壞,靶Ⅰ產(chǎn)生壓縮和剪切變形,使大塊混凝土破壞、碎裂,由于靶架約束,靶Ⅰ下半部分靶體處于原始位置,上半部分碎裂脹開鋼箍并飛散,靶板破壞嚴(yán)重,彈體大部分能量由靶Ⅰ吸收,彈體進(jìn)入靶Ⅱ,靶體表面形成對稱的局部破壞,破壞區(qū)分為中心粉碎區(qū)、漏斗形破壞區(qū)、徑向裂紋區(qū),靶Ⅱ破壞情況如圖5所示。由仿真結(jié)果可以看出,靶Ⅰ完全破壞,靶Ⅱ形成隧道,與試驗(yàn)結(jié)果基本一致。

圖5 靶Ⅱ破環(huán)情況

由圖5可知,彈體著速1 368 m/s時(shí),靶板入孔不規(guī)則,無徑向裂紋;著速提高至1 666 m/s時(shí),靶板入孔為圓孔,表面平整無崩落,在彈體侵徹后其彈著點(diǎn)呈放射型裂紋,裂紋一直延伸到靶體邊界;初速進(jìn)一步提高至1 724 m/s時(shí),靶板入孔明顯,靶板表面大面積崩落,裂紋至靶板邊緣,由此可見,隨著彈體初速的增加,靶Ⅱ破壞程度逐步加深。彈體在不同初速條件下侵徹混凝土靶板侵徹深度數(shù)值模擬和試驗(yàn)結(jié)果對比如表3所示。

表3 侵徹深度數(shù)值模擬與試驗(yàn)結(jié)果對比

由表3可知,侵徹深度數(shù)值模擬和試驗(yàn)結(jié)果誤差不大于7%。通過對比分析靶板破壞情況,極限侵徹深度,仿真數(shù)據(jù)和試驗(yàn)數(shù)據(jù)基本一致,說明該方法能夠較好地描述超高速侵徹物理過程,可再此基礎(chǔ)上進(jìn)行侵徹機(jī)理的探究。

2 超高速侵徹階段分析

第1節(jié)中因?yàn)樵囼?yàn)場地和設(shè)備等因素限制,靶板為2塊?1.5 m×1 m混凝土靶。為消除靶板尺寸效應(yīng),本節(jié)數(shù)值模擬的靶板尺寸為?1.5 m×2.5 m,彈體尺寸與第1節(jié)保持一致。

2.1 無量綱侵深

圖6為無量綱侵徹深度與速度的變化關(guān)系,由圖6可以看出,侵徹速度0.2～3.0 km/s區(qū)間內(nèi),無量綱侵徹深度隨速度先增加后減少,在1.2 km/s附近存在侵深最大值。彈體速度小于1.2 km/s時(shí),侵深隨速度的增加而增加,彈體速度在1.2 km/s附近侵徹深度達(dá)到最大值,約80倍彈徑,彈體速度大于1.2 km/s時(shí),無量綱侵深隨速度的增加而減小。因此,從侵徹深度方面討論,超高速彈體對于中低速彈體在侵徹能力上并無優(yōu)勢。

圖6 無量綱侵深與速度變化關(guān)系Fig.6 Penetration depth with velocity

2.2 彈體質(zhì)量侵蝕

圖7為不同速度下鎢桿侵徹混凝土彈體質(zhì)量時(shí)程圖,由圖7可知,彈體速度越大,質(zhì)量侵蝕速率越大,彈體速度0.2 km/s時(shí)彈體質(zhì)量基本無損失,隨著速度增加彈體質(zhì)量變化曲線下降段斜率逐漸增大,而3.0 km/s時(shí)彈體質(zhì)量幾乎損失殆盡。并且,彈體速度較低時(shí),彈體質(zhì)量下降段與平穩(wěn)段分界并不明顯,當(dāng)彈體速度大于1.4 km/s后,質(zhì)量曲線下降段與穩(wěn)定段轉(zhuǎn)折點(diǎn)清晰。

圖7 彈體質(zhì)量隨時(shí)間變化曲線

這是由于侵徹過程中,彈、靶相互作用產(chǎn)生高溫導(dǎo)致彈體表面材料軟化,尤其是作用劇烈的彈頭部,混凝土顆粒切削軟化的彈體表面材料,被切削下的彈體材料,一部分殘留在混凝土中而形成彈體質(zhì)量損失,另一部分與混凝土混合后隨彈體侵徹結(jié)束重新覆蓋于彈體表面,彈體速度越高材料軟化越多,被混凝土顆粒切削下來的材料越多。另一方面,隨著彈體速度的增加,侵徹機(jī)理發(fā)生改變,由剛體侵徹進(jìn)入流體侵徹,此時(shí)材料的強(qiáng)度在侵徹過程中作用逐漸減小,材料密度起主要作用。當(dāng)進(jìn)入流體侵徹,彈體材料侵蝕劇烈,彈體材料侵蝕,該階段稱為“流體侵徹階段”。當(dāng)彈體速度降低到彈靶界面壓力不足以侵蝕彈體時(shí),流體侵徹階段結(jié)束,此后彈體以剛體狀態(tài)繼續(xù)侵蝕,直到彈體速度降為0,該階段稱為“剛體侵徹階段”。不同速度下彈體質(zhì)量侵蝕率如表4所示。

表4 不同速度下彈體質(zhì)量侵蝕率

圖8給出了彈體質(zhì)量侵蝕率與初始動(dòng)能的關(guān)系,由圖8可知,0.5～1.6 km/s速度范圍彈體質(zhì)量侵蝕率與彈體初始動(dòng)能呈線性關(guān)系。根據(jù)Forrest et al,1996;Frew et al,1998的試驗(yàn)結(jié)果Silling and Forrest(2007)指出,在侵徹速度v≤1 km/s(或v2/2≤0.5(km2/s2))時(shí)彈體質(zhì)量的磨蝕與初始動(dòng)能存在線性關(guān)系;侵蝕速度v≥1 km/s后,質(zhì)量磨蝕維持在一常數(shù)水平。但是Forrest et al,1996;Frew et al,1998的試驗(yàn)彈體為不同標(biāo)號的鋼,而本文彈體為鎢合金,因此在彈體質(zhì)量的磨蝕與初始動(dòng)能存在線性關(guān)系的速度區(qū)間并不一致。所以不同彈體材料侵徹混凝土,彈體質(zhì)量磨蝕與初始動(dòng)能線性關(guān)系的速度區(qū)間要根據(jù)彈體材料而定。

圖8 彈體質(zhì)量侵蝕率和初始動(dòng)能呈線性關(guān)系

2.3 殘余彈長

圖9為殘余彈長隨速度變化曲線,由圖9可以看出,速度0.2 km/s時(shí)殘余彈長33.5 cm,與初始長度基本相同。1 km/s彈體殘余長度21.1 cm,侵蝕掉原長的40.9%,隨著速度的增加彈體被劇烈侵蝕,當(dāng)速度3.0 km/s彈體基本被完全侵蝕。侵徹開始時(shí)刻,靶板逐漸開坑,彈體0.2～1.5 km/s速度侵徹靶板時(shí),彈體材料并未全部侵蝕,說明彈體既經(jīng)歷了“流體侵徹”又經(jīng)歷了“剛體侵徹”。當(dāng)速度大于2.0 km/s時(shí),彈體已經(jīng)被完全侵蝕,因此不包括剛體侵徹階段。由于剛體侵徹深度在總侵徹深度中占很大比例,失去剛體侵徹階段會(huì)使總侵深降低,當(dāng)速度大于2.0 km/s時(shí)總侵徹深度隨速度的增加明顯降低。

圖9 殘余彈長隨速度變化曲線

2.4 侵徹階段分析

根據(jù)Eichelberger[14]和orphal[15]分別提出的撞擊4個(gè)階段和第3侵徹階段,將超高速侵徹可分為4個(gè)階段,即:① 瞬態(tài)高壓階段;② 流體侵徹階段;③ 剛體侵徹階段;④ 回彈階段。由于瞬態(tài)高壓階段和回彈階段侵徹深度占比很小,可以忽略不計(jì),因此本節(jié)重點(diǎn)關(guān)注流體侵徹階段和剛體侵徹階段。

彈體速度的增加會(huì)導(dǎo)致侵徹機(jī)制由剛體侵徹向流體侵徹轉(zhuǎn)變,在某一速度下2種侵徹機(jī)制均會(huì)出現(xiàn),因此研究剛體侵徹和半流體侵徹在侵徹過程中的占比,對研究超高速侵徹混凝土十分重要。

彈體尾部速度、彈體質(zhì)量、侵徹深度隨時(shí)間曲線如圖10所示。由圖10可知,彈體初速為0.6 km/s時(shí),彈體質(zhì)量和侵徹深度隨時(shí)間呈連續(xù)線性變化,質(zhì)量侵蝕率為15.32%,可以判斷侵徹全過程基本為剛體侵徹。彈體初速增加至1.0 km/s時(shí),彈體質(zhì)量變化曲線出現(xiàn)了明顯的轉(zhuǎn)折點(diǎn),侵度-時(shí)間曲線變?yōu)榉蔷€性,這是由于,侵徹開始瞬間彈體相對于靶體做減速運(yùn)動(dòng),產(chǎn)生一個(gè)強(qiáng)沖擊波向彈體和靶板傳遞,彈體承受極高的壓力,此壓力超過彈體材料所承受范圍,導(dǎo)致彈體材料發(fā)生塑性變形,侵徹機(jī)制變?yōu)榱黧w侵徹,彈體被侵蝕,質(zhì)量減少,阻力變大,彈體速度逐漸下降。當(dāng)彈體速度低于剛體臨界速度vcri,侵徹機(jī)制轉(zhuǎn)換為剛體侵徹,這個(gè)時(shí)間為1.5 ms(圖9(b)),此時(shí)彈體速度為572 m/s,在流體侵徹階段(0～1.5 ms)和剛體侵徹階段(1.5～3.6 ms)彈體質(zhì)量分別減少30.5%、12.1%,無量綱侵深分別為47.2、46.3,分別占總侵深50.4%、49.6%。這一現(xiàn)象在初速為1 400 m/s時(shí)更為顯著,與初速1.0 km/s不同的是,初速1.4 km/s情況下,在1.8 ms時(shí)彈體質(zhì)量變化曲線存在明顯拐點(diǎn),此時(shí)彈體速度為406 m/s。在流體侵徹階段(0～1.8 ms)和剛體侵徹階段(1.8～3.6 ms)內(nèi)彈體質(zhì)量分別減少66.2%、6.8%,侵深占比分別為70.4%、29.5%。

圖10 彈體尾部速度、彈體質(zhì)量、侵徹深度隨時(shí)間曲線

表5為不同初始速度下彈體流體/剛體侵徹深度數(shù)值模擬結(jié)果,圖11為不同初速下彈體流體/剛體侵徹變化曲線。

表5 不同初速下彈體流體/剛體侵徹深度數(shù)值模擬結(jié)果

圖11 不同初速下彈體流體/剛體侵徹變化曲線

通過表5和圖11可以看出,彈體初速較低時(shí)侵徹全過程為剛體侵徹,初速達(dá)到1.0 km/s時(shí)彈體材料開始侵蝕,出現(xiàn)流體侵徹階段,此時(shí)流體侵徹深度和剛體侵徹深度基本相等。隨著彈體初速的增加,流體階段侵徹深度呈現(xiàn)先增加后維持在一常數(shù)附近,剛體侵徹階段侵深逐漸減小,兩者疊加之后得到的總侵深呈現(xiàn)先增大后減小的趨勢。

3 結(jié)論

本文中通過試驗(yàn)驗(yàn)證方法驗(yàn)證了材料本構(gòu)模型的準(zhǔn)確性,又通過數(shù)值模擬方法進(jìn)一步研究了鎢合金桿彈超高速侵徹混凝土靶侵徹深度隨彈體初速的變化規(guī)律,得出以下結(jié)論:

1) 分析數(shù)值模擬所得數(shù)據(jù)可知,流體侵蝕階段,彈體質(zhì)量大幅度減少。0.2～1.8 km/s速度范圍內(nèi),彈體質(zhì)量侵蝕率與彈體初始動(dòng)能呈正比例線性關(guān)系,彈體初速大于2.0 km/s質(zhì)量侵蝕率高達(dá)95%。

2) 由數(shù)值模擬分析可得,侵徹深度隨彈體初速增加呈現(xiàn)先增加后減小現(xiàn)象,在初速1.2 km/s附近存在侵深最大值約80倍彈徑。超高速條件下侵徹深度對于中低速侵徹并無優(yōu)勢。

3) 彈體初速小于0.8 km/s侵徹機(jī)制為剛體侵徹,初速大于1.0 km/s時(shí)侵徹機(jī)制既包括剛體侵徹又包括流體侵徹。1.0～2.0 km/s流體侵深由47.2增加至68.6,3.0 km/s時(shí)降至56.1。在1.0～3.0 km/s剛體侵徹深度由46.3持續(xù)下降至1。流體侵深和剛體侵深的變化,導(dǎo)致體侵深曲線呈現(xiàn)先增大后減小的趨勢。