朱東俊 王慶周 孫登成 趙 勇 龔小超

（哈爾濱工程大學 船舶工程學院 哈爾濱150001）

引 言

艦船在海上執(zhí)行任務時，容易受到穿甲彈的打擊威脅，穿甲彈利用其高速動能鉆入艦船內部，然后通過自身爆炸充分發(fā)揮其毀傷能效。由于軍事領域的發(fā)展需要，彈體侵徹問題上的研究投入也日益增多。侵徹問題可分為彈體對靶板的侵入、在靶板內運動和彈侵穿出等三個主要階段。

侵徹問題研究主要有實驗、理論分析和數(shù)值模擬三類基本方法［1-2］，初期多以實驗研究為基礎，然而實驗研究費用高昂，而且無法完全獲得實驗過程中所有侵徹現(xiàn)象與數(shù)據。理論分析方法對侵徹過程中的主要力學機理可以達到一定程度的反映，獲得相應的指導性結論；然而在理論分析中必然對材料性質和侵徹過程進行簡化，可能導致較大誤差甚至錯誤結果。隨著計算機性能的提高和研究方法的進步，彈侵研究主要采用數(shù)值模擬。數(shù)值模擬通過離散守恒方程，定義正確的本構模型和狀態(tài)方程以準確描述材料性能變化與力學過程，通過對不同條件下的侵徹問題計算分析，準確有效地獲得靶板破壞模式、侵徹過程與特性等的信息，且克服了實驗研究中的準備周期長、參數(shù)測量困難和成本高等不足，具備重復研究性。

數(shù)值模擬主要以有限元法為主。馬志濤等［3］運用AUTODYNA軟件建立彈丸高速侵徹厚鋁板的有限元模型，通過研究發(fā)現(xiàn)彈坑的深度和寬度隨著彈丸的撞擊角增加而減少，也隨著撞擊速度的增加而增加。周平等［4］采用LS-DYNA軟件準確模擬了鋁合金板在方形預制破片侵徹過程中的演變情況，設置相應的侵徹工況，研究不同速度和姿態(tài)下，破片的損傷模式、剩余速度以及彈道極限。梅志遠等［5］基于MSC/DYTRAN程序，通過考察立方體破片以 3種典型姿態(tài)侵徹厚度為4 mm的945鋼，得出了侵徹姿態(tài)不同且侵徹威力相差不大的結論。本文基于ABAQUS軟件建立相應的穿甲彈侵徹靶板的有限元數(shù)值模型，研究分析了侵徹過程中損傷機理以及靶板的損傷破壞模式，取得的數(shù)值結果合理準確，為日后穿甲彈侵徹與艦船防御的深入研究打下了良好的基礎。

1 材料失效理論

1.1 材料本構模型

彈侵過程中包含結構的大變形乃至損傷與破壞，材料的損傷與破壞問題是一個極其復雜與關鍵的問題。在有限元數(shù)值模擬中，正確定義材料的本構模型和狀態(tài)方程，是準確研究和描述材料力學行為的關鍵前提。在選擇描述沖擊問題的固體材料本構模型時，本文采用Johnson-cook材料屈服模型［6］計算金屬的屈服應力Y，其表達式為：

在穿甲彈侵徹目標靶板數(shù)值模型中，彈體與靶板分別采用鋼材料與鋁材料進行模擬，運用Johnson-cook材料模型定義穿甲彈和靶板的本構關系，彈體與靶板材料的本構關系參數(shù)如表1所示。

金屬材料的狀態(tài)方程采用固體力學中的Mie-Gruneisen狀態(tài)方程［7］進行描述：

表1 不同金屬材料本構關系參數(shù)

式中：Γ為Gruneisen常數(shù)，pH為沖擊Hugoniot曲線上的點的壓力，其表達式為：

表2 不同金屬材料Mie-Gruneisen狀態(tài)方程中參數(shù)

1.2 材料失效及單元刪除

此外，ABAQUS中侵徹模型要定義材料失效，實際采用與時間相關的Johnson-cook失效模型，引入失效參數(shù)D，其表達式如下：

表3 不同金屬材料中Johnson-cook失效模型參數(shù)

彈侵沖擊屬于大變形問題，在定義材料損傷失效后，侵徹過程中材料也將發(fā)生非線性失效。當材料沒有損傷時，D為0；而D為1時，則表示材料完全發(fā)生破壞，即認為該單元失效，此時ABAQUS中的標記狀態(tài)變量STATUS為1，該失效單元即被刪除，不再參加后續(xù)計算。

隨著彈體不斷侵入靶板，彈體與靶板相接觸的單元定會產生擠壓、扭曲等大變形問題，這將導致計算結果失真乃至計算終止。為了確保侵徹過程正常進行與滿足計算精度要求，在有限元中必然需要對嚴重畸變的網格進行網格重新劃分或刪除處理，本文將采用ABAQUS中的單元刪除技術處理網格畸變問題［9］。

2 有限元模型的建立

本文采用大型有限元軟件 ABAQUS 對穿甲彈侵徹靶板進行模擬分析。在圖1的三維彈侵穿甲數(shù)值模型中，將穿甲彈簡化為圓柱結構：長度H為30 mm，截面直徑D為10 mm。選用鋼材料進行模擬。靶板的長l為100 mm、寬b為100 mm、高h為10 mm，正方形中心區(qū)域邊長d為30 mm，材料為鋁合金。彈體以900 m/s的初始速度對靶板進行垂直侵徹，靶板四周邊界采用絞支約束，即只限制四邊的三方向位移，不約束旋轉自由度。

圖1 穿甲彈垂直侵徹靶板三維模型

在網格劃分階段，穿甲彈與靶板都采用八節(jié)點六面體單元進行劃分?？紤]到侵徹過程靶板中心區(qū)域網格變形較為劇烈，因此為提高計算精度，所以對靶板沖擊中心30 mm×30 mm區(qū)域網格進行細化處理，網格尺寸為1 mm×1 mm×1 mm，然后依次過渡網格尺寸到邊界區(qū)域，在保證網格質量的同時減少計算時間。在定義穿甲彈與靶板之間的接觸時，因為侵徹沖擊過程中接觸關系較為復雜，并且考慮到彈體侵入靶板過程中不可避免的存在失效單元刪除，這也涉及內部單元接觸問題。采用對接觸面類型限制較少的通用接觸算法［9］，所有實體的面都定義上自接觸，有效地改善單元接觸狀態(tài)，確保侵徹沖擊過程中穿甲彈與靶板的正確接觸。最終建立穿甲彈垂直侵徹靶板三維有限元模型如圖2所示。

在ABAQUS中穿甲彈侵徹靶板模型的初始條件設置與計算循環(huán)流程如圖3所示。

圖2 穿甲彈垂直侵徹靶板三維有限元模型

圖3 ABAQUS彈侵模型計算流程圖

3 計算結果

為了驗證本文穿甲彈侵徹靶板數(shù)值模型的有效性，圖4給出了本文的侵徹剩余速度變化曲線與參考文獻［10］的對比結果。從圖中可以看出，兩種結果的彈體沖擊速度時歷曲線的趨勢是一致的。在穿甲彈侵入靶板的過程中，伴隨著單元的失效刪除，導致沖擊過程中接觸面積變化，所以速度曲線在沖擊接觸階段振蕩比較明顯，然而數(shù)值結果依然準確可信。通過分析ABAQUS沖擊速度結果發(fā)現(xiàn)，在侵徹過程前期，彈體以較高的速度撞擊到靶板上，由于受到靶板材料較大的抗力，彈體的部分動能在較短時間內被消耗，導致其速度在較短時間內迅速衰減；在55 μs左右，隨著彈體逐漸穿透靶板及靶板材料的破壞與失效飛濺，靶板對彈體沖擊的抗力亦逐漸衰減，彈體速度逐漸趨于平穩(wěn)，整個穿甲過程基本完成。參考文獻中，穿甲彈最終剩余速度是829.02 m/s，而本文侵徹數(shù)值模型中的結果是826.49 m/s，兩者結果吻合良好，驗證了本文建立的穿甲彈侵徹靶板模型的準確有效性。

圖4 穿甲彈垂直侵徹靶板過程中彈體速度時歷曲線對比

圖5給出了側面視圖和正面視圖下穿甲彈侵徹靶板過程中典型時刻的Mises應力分布圖。從圖中可以較清晰地看到靶板從正表面被彈體碰撞的初期直到完全穿透的完整動態(tài)響應過程。

從圖5中可以看出，彈體垂直侵徹靶板過程中，靶板材料被不斷侵入的彈體擠壓而向四周（包括自由表面方向）高速運動并向外膨脹，因而在靶板正表面逐漸出現(xiàn)金屬材料的堆積，形成翻起的唇邊（靶前花瓣型變形）［11］，彈體與靶板之間的高速碰撞沖擊產生應力波且在結構中迅速傳遞開來。同時彈體頭部材料在靶板軸向抗力的作用下發(fā)生塑性破壞，從而導致彈體頭部迅速鈍化，最終呈“蘑菇狀”，如下頁圖6所示。隨著彈體的繼續(xù)侵徹，靶板正表面上，排開和飛濺的靶板材料逐漸增多并持續(xù)堆積，同時靶板背表面上逐漸形成“鼓包”；最終彈體以一定的剩余速度穿透“鼓包”，整個穿甲過程完成。

在此過程中，受到彈體高速沖擊作用的靶板部分產生相對于靶板結構其他部分的運動，形成剪切帶，且剪切帶內的材料劇烈變形；此外由于剪切帶內的材料變形較快，因而變形生成的熱量來不及傳遞到周圍結構中，而在剪切帶內積累，從而使得該部分結構溫度大大升高，且材料軟化、強度降低，當這種材料降低的程度超過材料因變形硬化而提高的強度時，就會發(fā)生剪切帶材料的失穩(wěn)，即所謂的絕熱剪切過程。在彈體的沖擊作用下，靶板結構中剪切帶不斷發(fā)展，最終連接、交叉形成柱面形狀的塞塊，并被彈體的剩余能量推出靶板背面，完成整個沖塞作用過程。這一系列現(xiàn)象與國內外相關實驗［12］得到的結果較為一致，因而進一步證實本文建立的穿甲彈侵徹靶板模型是準確有效的。

圖5 穿甲彈垂直侵徹靶板過程典型時刻Mises應力分布

圖6 侵徹過程中穿甲彈典型時刻彈體頭部形狀

4 結 論

本文基于Johnson-cook本構方程定義描述靶板金屬材料在侵徹大變形中的行為特性，考慮失效參數(shù)影響，運用ABAQUS建立了穿甲彈侵徹靶板的有限元模型，對穿甲彈侵徹靶板的過程及動力響應進行分析，結論如下：

（1）ABAQUS中對于沖擊侵徹等大變形問題的處理，采用通用接觸算法，考慮靶板內部接觸關系，結合單元刪除技術，可以有效地消除單元畸變而導致計算結果錯誤甚至計算終止問題。

（2）通過與參考文獻對比，發(fā)現(xiàn)單元失效刪除的定義，導致侵徹沖擊接觸過程中接觸面積發(fā)生變化，速度變化曲線出現(xiàn)震蕩，然而數(shù)值模擬結果依然有效可信，充分驗證了本文建立的數(shù)值模型在預測侵徹過程中的動力響應及損傷的合理準確性。

（3）平頭穿甲彈垂直侵徹靶板的損傷模式主要是沖塞剪切破壞。篇幅所限，本文只討論了平頭彈的垂直侵徹，針對彈頭形狀、侵徹角度等影響研究，均可采用同樣方法開展進一步研究。

［1］ 尉志同，余文力，王濤，等.動能彈垂直侵徹混凝土靶板的三維數(shù)值模擬研究［J］.彈箭與制導學報，2007（4）：121-123.

［2］ 王政.彈靶侵徹動態(tài)響應的理論與數(shù)值分析［D］.復旦大學，2005.

［3］ 馬志濤，張偉，賈斌，等.彈丸超高速撞擊半無限厚鋁板數(shù)值模擬［J］.材料科學與工藝，2005（3）：283-286.

［4］ 周平，姚武文，侯日立.方形破片對鋁板的侵徹研究［J］.機電產品開發(fā)與創(chuàng)新，2011 （3）：21-22.

［5］ 梅志遠，朱錫.利用 MSC/DYTRAN 程序仿真分析導彈戰(zhàn)斗部立方體破片的侵徹威力［J］.海軍工程大學學報，2002（2）：39-42.

［6］ （美）邁耶斯著.張慶明，劉彥，黃風雷，等譯.材料的動力學行為［M］.北京：國防工業(yè)出版社.2006.

［7］ Shin Y S，Lee M，Lam K Y，et al.Modeling mitigation effects of water shield on shock waves［J］.Shock and Vibration，1998（5）：225-234.

［8］ 侯滿義，李曙林，孫旭，等.飛機壁板結構戰(zhàn)傷的動力有限元仿真［J］.空軍工程大學學報：自然科學版，2007（1）：1-3.

［9］ Dassault Systemes.ABAQUS 6.10-1 Documentation.Abaqus Analysis User’s Manual（21.3.1-1）-（21.3.2-4）［S］.2010.

［10］ 初文華.處理非連續(xù)問題的三維SPH算法及其在沖擊動力學問題中的應用［D］.哈爾濱：哈爾濱工程大學，2013.

［11］ 趙國志.穿甲工程力學［M］.北京：兵器工業(yè)出版社，1992.

［12］ 李金泉，黃德武，段占強，等.高速侵徹裝甲鋼絕熱剪切帶特性研究［J］.彈道學報，2003（3）：86-91.