串聯(lián)隨進彈侵徹預(yù)開孔靶彈道軌跡的數(shù)值模擬

鄧佳杰，張先鋒，陳東東，郭磊，何勇

（南京理工大學(xué)機械工程學(xué)院，江蘇南京210094）

串聯(lián)隨進彈侵徹預(yù)開孔靶彈道軌跡的數(shù)值模擬

鄧佳杰，張先鋒，陳東東，郭磊，何勇

（南京理工大學(xué)機械工程學(xué)院，江蘇南京210094）

利用移動載荷自定義程序的彈-靶分離方法，并結(jié)合Forresta1半經(jīng)驗靶體阻力函數(shù)，開展了串聯(lián)隨進彈侵徹預(yù)開孔靶彈道軌跡數(shù)值模擬。在實驗驗證計算模型及自定義程序可靠的基礎(chǔ)上，研究彈-靶軸偏置、傾角、攻角等對串聯(lián)隨進彈侵徹彈道軌跡的影響規(guī)律。研究結(jié)果表明：彈-靶軸偏置、傾角、攻角等因素對彈體的彈道軌跡及侵徹深度影響顯著；預(yù)開孔孔道具有一定的引導(dǎo)侵徹作用；在彈-靶軸偏置、傾角及攻角的共同影響下，隨進彈侵徹預(yù)開孔靶存在著跳飛可能。

兵器科學(xué)與技術(shù)；彈道軌跡；預(yù)開孔靶；彈-靶分離方法；數(shù)值模擬

0 引言

破甲/爆破型串聯(lián)戰(zhàn)斗部侵徹問題是武器研制和工程防護的重點問題之一。串聯(lián)隨進彈的侵徹能力是評估戰(zhàn)斗部作用性能的關(guān)鍵指標。目前，國內(nèi)外學(xué)者針對串聯(lián)隨進彈垂直侵徹預(yù)損傷靶或預(yù)開孔靶問題開展了相關(guān)研究工作。Te1and［1］首次將空腔膨脹理論運用于預(yù)開孔侵徹問題。基于Te1and理論，張雷雷等［2］引入靶體損傷系數(shù)修正侵徹深度計算式。王樹有［3］針對彈體頭形、靶體孔形及損傷類型等影響因素分析后級隨進侵徹問題。Guo等［4］引入Wen半經(jīng)驗公式［5］建立了靶體阻力函數(shù)隨進侵徹四階段計算模型。文鶴鳴等［6］發(fā)展了初始空腔直徑不為0且考慮靶體材料可壓縮性和剪脹特性的球形空穴膨脹理論，并進一步研究前后級匹配優(yōu)化設(shè)計問題。當(dāng)前，針對隨進彈侵徹預(yù)開孔靶的彈體變形及彈道偏轉(zhuǎn)問題主要采用傳統(tǒng)仿真方法進行研究［7-9］，該方法涉及三維彈-靶網(wǎng)格及復(fù)雜接觸過程求解，計算資源占用大，不利于開展大量的仿真研究，因此采用簡化且準確的計算方法意義重大。

彈-靶分離方法是一種快速有效預(yù)測彈體侵徹過程的數(shù)值方法，其結(jié)合理論模型的靶體阻力公式及有限元模擬的離散化計算程序［10］，將靶對彈體作用以靶體阻力代替，省略靶體建模及彈-靶接觸，提高了計算效率。國內(nèi)外學(xué)者針對侵徹問題已開展大量彈-靶分離數(shù)值模擬，其計算結(jié)果具有較高的精度。Warren等［11-12］首次提出采用彈-靶分離方法進行侵徹研究，并在PRONTO 3D軟件中實現(xiàn)不同傾角下的彈體對6061鋁靶的侵徹過程；基于彈-靶分離思想，何濤等［13-14］、孔祥振等［15］、Fang等［16］等采用ABAQUS軟件的自定義程序開發(fā)功能實現(xiàn)了彈體侵徹纖維增強復(fù)合材料、混凝土、巖石等靶體，數(shù)值仿真結(jié)果與實驗吻合較好；Li等［17］驗證了彈靶分離數(shù)值模擬的準確性，并開展了非正侵徹條件下彈體侵徹彈道及跳飛極限研究；在前人研究基礎(chǔ)上，劉天宋等［18］開展了攻角、侵徹速度、彈體質(zhì)心位置等因素影響下彈體侵徹半無限混凝土靶彈道偏轉(zhuǎn)問題的彈靶分離數(shù)值模擬研究。

本文以彈靶分離方法為基礎(chǔ)，借助于有限元軟件的移動載荷自定義程序VDLOAD，開展了串聯(lián)隨進彈侵徹預(yù)開孔靶彈道軌跡彈靶分離數(shù)值模擬研究。在此基礎(chǔ)上，通過分析彈體的侵徹深度及其質(zhì)心彈道運動軌跡來研究彈-靶孔軸偏置、傾角、攻角等因素對彈道軌跡偏轉(zhuǎn)的影響規(guī)律。

1 隨進彈侵徹預(yù)開孔靶的有限元數(shù)值模擬

1.1有限元模型

由于串聯(lián)隨進彈侵徹預(yù)開孔靶體過程十分復(fù)雜，彈體可能偏離預(yù)定的侵徹方向產(chǎn)生彈道偏轉(zhuǎn)，從而影響隨進侵徹性能。以卵形彈侵徹預(yù)開孔靶為例，在不考慮串聯(lián)戰(zhàn)斗部前級開孔引起的孔道周圍混凝土強度弱化的前提下，影響其侵徹軌跡的主要因素是預(yù)開孔靶孔軸與彈軸偏置H、彈體傾角φ、攻角δ、靶體預(yù)開孔孔半錐角θ及開孔孔徑Rh.在實際應(yīng)用中，隨進彈侵徹預(yù)開孔靶為三維問題，為了簡化問題，假設(shè)彈體初始攻角、傾角及彈-靶孔軸偏置間的坐標位置均在一個平面內(nèi)。圖1為二維平面上彈體侵徹預(yù)開孔靶及侵徹過程相關(guān)影響因素關(guān)系示意圖。

圖1　彈體侵徹預(yù)開孔靶示意圖Fig.1 Penetration of trave11ing projecti1e into pre-dri11ed target

借助于ABAQUS有限元程序的顯式動態(tài)分析模塊，模擬串聯(lián)隨進彈侵徹預(yù)開孔靶彈道軌跡。圖2為彈體幾何模型及靶體示意圖，在數(shù)值模擬中采用C3D8R實體單元對彈體三維全模型進行網(wǎng)格劃分，并運用Lagrange算法對模型進行求解。全彈表面施加自定義程序靶體阻力，以正應(yīng)力形式代替侵徹過程中靶體與全彈表面間的作用。以圖2所示的坐標軸為基準編寫自定義程序，由靶體預(yù)開孔半錐角θ及開孔孔徑Rh可確定靶體預(yù)開孔的坐標及孔形。通過改變彈體相對于參考系原點的相對位置，實現(xiàn)模型中考慮相應(yīng)的影響因素，進而分析其對侵徹彈道軌跡的影響。

圖2　彈體侵徹預(yù)開孔靶有限元模型Fig.2 Finite e1ement mode1 of trave11ing projecti1e penetrating into pre-dri11ed target

非垂直侵徹情況下，彈體表面受力不均勻可能導(dǎo)致其產(chǎn)生塑性變形，數(shù)值模擬中彈體材料模型采用Johnson-Cook本構(gòu)模型來表征其變形，即

式中：σ為Von Mise流動應(yīng)力；A、B、C、m、n為材料常數(shù)為等效塑性應(yīng)變?yōu)樗苄詰?yīng)變率為參考應(yīng)變率；T為實際溫度；Tr為室溫；Tm為融化溫度。

1.2靶體阻力施加的自定義程序?qū)崿F(xiàn)

基于空腔膨脹理論的靶體阻力函數(shù)確定方法，并計及靶體強度項［17，19］，其阻力函數(shù)形式可表示為

式中：ai為靶體材料相關(guān)常數(shù)項；vn為單元法向速度。

根據(jù)微分面力法（DAFL）［20］及局部相互作用模型（LIMs）［21］，彈-靶分離方法中彈體受力由彈體表面單元阻力疊加得到。結(jié)合自定義程序編寫的彈-靶接觸區(qū)域及靶體阻力函數(shù)，在有限元軟件中求解離散化彈體局部表面運動參數(shù)，疊加單元阻力以確定彈體表面受力［22］。

如圖3所示，在侵徹過程的某一時刻，運動速度單位矢量v0處的彈體任意表面單元受沿內(nèi)法向單位矢量n0和切向單位矢量τ0方向的兩部分作用力，外法向矢量與速度矢量夾角為 ^υ.判定彈-靶接觸分離條件為v0n0＜0，若滿足不等式要求，則作用于彈體表面單元的法向和切向阻力可表示為

式中：μ為庫侖摩擦系數(shù)。

圖3　彈體表面單元受力示意圖［21］Fig.3 Resistance of projecti1e surface e1ement［21］

以彈-靶分離方法為基礎(chǔ)，結(jié)合彈-靶接觸位置及靶體阻力函數(shù)定義，編寫移動載荷自定義程序VDLOAD，并將其作為邊界條件施加于彈體表面用來替代彈-靶實際接觸作用力。

圖4為自定義移動載荷程序計算流程圖，其中孔道坐標及孔形由yi處預(yù)開孔半徑的計算結(jié)果表示。在任意時間步，表面單元受力受控于彈-靶相對位置及單元法向速度，當(dāng)滿足流程圖中的彈體進入預(yù)開孔靶體區(qū)域、彈體侵徹速度不為0及彈-靶接觸判定條件時，程序計算該時刻彈體表面正應(yīng)力，通過疊加所有單元所受阻力確定該時刻彈體的合阻力。在設(shè)定時間步長內(nèi)循環(huán)調(diào)用VDLOAD，直至每一單元的法向速度為0停止調(diào)用，數(shù)值模擬計算終止。

圖4　移動載荷自定義程序計算流程圖Fig.4 F1ow chart of user subroutine of moving 1oad

侵徹過程模擬中彈-靶摩擦作用通過在自定義程序中加入摩擦系數(shù)進行考慮，由此彈體表面任意單元上某一節(jié)點的軸向阻力可表示為

通過（5）式的等效節(jié)點正應(yīng)力表示方法，可用于各種侵徹條件下全彈表面受力分析，以預(yù)測考慮切向摩擦和法向應(yīng)力共同作用下的彈體侵徹過程。

1.3有限元模型的實驗驗證

為了驗證彈-靶分離方法的可靠性，開展了模擬彈體垂直侵徹預(yù)開孔鋁靶實驗研究，并結(jié)合已有彈體侵徹預(yù)開孔混凝土靶實驗數(shù)據(jù)［23］，對比了數(shù)值模擬結(jié)果與實驗結(jié)果。

1.3.1模擬彈體垂直侵徹預(yù)開孔鋁靶實驗

實驗中模擬彈體為尖卵形頭部，直徑為14.5 mm（見圖5），長徑比為6，頭部曲徑比CRH為3.彈體材料為30CrMnSiNi2A，硬度HRC45～48，彈體質(zhì)量67 g.實驗靶體為2024鋁靶，直徑100 mm，靶體類型為中心開孔（8 mm、11 mm直徑柱形通孔）和半無限靶兩類。實驗通過彈道槍發(fā)射模擬彈體，侵徹速度v范圍為300～800 m/s.

圖5　實驗彈實物圖Fig.5 Photograph of trave11ing projecti1e

如表1所示，共進行5組不同速度段侵徹3類靶體的對比實驗。從回收彈體可以看出彈體基本未變形，通過如圖6所示的線切割侵徹后靶體觀察彈道和侵徹深度。實驗結(jié)果表明模擬彈體侵徹彈道準直性好、與普通半無限鋁靶侵徹結(jié)果相比，模擬彈體對不同直徑預(yù)開孔鋁靶侵徹深度均有一定程度的提高。

1.3.2數(shù)值模擬與實驗結(jié)果對比

數(shù)值模擬中彈體、靶體結(jié)構(gòu)參數(shù)與實驗一致，材料模型為Johnson-Cook材料模型，其參數(shù)如表2所示。鑒于Forresta1半經(jīng)驗公式更適用于混凝土侵徹問題［26-27］，同時為保證半經(jīng)驗阻力擬合方法一致，彈體侵徹鋁靶阻力函數(shù)亦采用Forresta1半經(jīng)驗公式［28］，侵徹混凝土材料及鋁靶材料過程中彈體所受靶體阻力見表3.

表1　侵徹深度實測值Tab.1 Fxperimenta1 resu1ts of projecti1es penetrating into A1 targets

圖6　實驗后半剖靶體實物圖Fig.6 Photographs of cutted A1 targets after impact

Fo1som［23］分別開展卵形彈體侵徹不同孔徑預(yù)開孔混凝土靶實驗，實驗細節(jié)見文獻［23］。該實驗中靶體尺寸較小，侵徹實驗結(jié)果可能存在邊界效應(yīng)影響［1］，將預(yù)開孔侵徹實驗數(shù)據(jù)進行如下的無量綱處理：

式中：P為任意預(yù)開孔情況的侵徹深度；PRh=0為相同條件下彈體侵徹半無限靶的侵徹深度。圖7分別給出了無量綱實驗結(jié)果、Fo1som經(jīng)驗公式及數(shù)值模擬計算的無量綱侵徹深度。從計算結(jié)果可以看出，相對于經(jīng)驗公式及半經(jīng)驗公式［4］，數(shù)值模擬結(jié)果與經(jīng)驗公式計算結(jié)果及實驗數(shù)據(jù)吻合較好，僅在相對孔徑Rh/Rp＞0.8處，數(shù)值模擬未考慮彈身段側(cè)壁與靶體間的作用，計算結(jié)果相對實驗數(shù)據(jù)偏高。

表2　彈體材料Johnson-Cook本構(gòu)模型參數(shù)Tab.2 Materia1 parameters of Johnson-Cook mode1

表3　靶體阻力模型Tab.3 Resistance function and materia1 parameters of target

圖7　侵徹預(yù)開孔混凝土仿真與經(jīng)驗?zāi)Ｐ图皩嶒瀸Ρ菷ig.7 Comparison of simu1ation，empirica1 mode1 and experimenta1 resu1ts of penetrated concrete targets

根據(jù)本文所開展的卵形彈體侵徹預(yù)開孔鋁靶實驗數(shù)據(jù)，進行對應(yīng)的數(shù)值模擬。如圖8所示為針對本文的卵形彈體侵徹預(yù)開孔鋁靶實驗所開展的數(shù)值模擬計算，結(jié)果表明數(shù)值模擬結(jié)果與實驗數(shù)據(jù)吻合較好。由于垂直侵徹模型與實驗條件可能存在差異，且未考慮彈身側(cè)壁與靶體間的作用，模型計算結(jié)果存在一定偏差，但最大誤差在10%以內(nèi)。總體而言，基于彈-靶分離方法的數(shù)值模擬能夠較好預(yù)測彈體侵徹預(yù)開孔靶過程。

圖8　侵徹預(yù)開孔鋁靶仿真結(jié)果與實驗結(jié)果對比Fig.8 Comparison of simu1ation and test resu1ts of penetrated 2024 A1 target

2 彈道軌跡偏轉(zhuǎn)影響因素分析

為了分析大口徑串聯(lián)隨進彈侵徹中導(dǎo)致其侵徹軌跡偏轉(zhuǎn)的影響因素，參考Chiesa等［8］的工作，采用如圖 9所示的 BLU-109彈體結(jié)構(gòu)。彈體由AISI 4340殼體和替代炸藥的惰性填充物組成，彈體質(zhì)量為475.36 kg.殼體材料本構(gòu)參數(shù)選用表2中的參數(shù)，填充物采用理想彈塑性本構(gòu)模型［31］，其彈性模量為10.1 GPa，屈服強度44.85 MPa.靶體為半無限混凝土，強度為30 MPa，其中心有孔徑254 mm的柱形通孔。在保持串聯(lián)隨進彈低速侵徹條件（v= 365 m/s）不變的情況下，研究不同侵徹條件下的隨進彈彈道軌跡。

圖9　模型彈幾何尺寸（CRH=3）［8］Fig.9 Geometry of fo11owing through projecti1e （CRH=3）［8］

2.1彈-靶孔軸線偏置對侵徹彈道軌跡的影響

由于串聯(lián)隨進彈在侵徹過程中可能受到不對稱力的作用，使得隨進彈軸線與靶體預(yù)開孔中心軸線偏離。此處分別計算6組不同偏離條件|H/Rh|下的彈體無攻角垂直侵徹預(yù)開孔混凝土靶過程。圖10為不同|H/Rh|值彈體質(zhì)心的侵徹運動軌跡。由圖10可以看出，彈-靶孔軸線偏置會使彈道軌跡產(chǎn)生多次偏斜，侵徹過程中，彈體表面受力不均與靶體一側(cè)接觸作用后又沿最小阻力方向運動直至與靶體另一側(cè)碰撞，往復(fù)多次直至侵徹結(jié)束。當(dāng)|H/Rh|≤0.4時，彈道偏斜較小且對侵徹深度影響有限，當(dāng)|H/Rh|＞0.4時，彈體在與靶體多次碰撞接觸過程中能量損失較大，使得侵徹深度顯著下降。彈體初始及最終姿態(tài)如圖11所示，|H/Rh|≥0.4彈體產(chǎn)生明顯變形，侵徹深度減小，但最終彈體位置均在初始預(yù)開孔彈道中，|H/Rh|值的變化并不影響靶體預(yù)開孔的引導(dǎo)侵徹作用。

圖10　不同彈-靶軸線偏置條件的彈體質(zhì)心運動軌跡Fig.10 Trajectory of mass center with different offsettings

2.2彈體初始傾角對侵徹彈道軌跡的影響

串聯(lián)隨進彈在不對稱力作用下，彈體繞彈軸上某點偏斜而產(chǎn)生傾角。設(shè)定彈體在接觸靶體前以傾角φ為0°、10°、20°、30°及40°向（-x，-y）方向侵徹，且攻角及彈-靶孔軸偏置均為0。對相同速度下不同傾角的隨進彈侵徹過程進行仿真模擬，所得最終侵徹深度結(jié)果變化如圖12所示，由計算結(jié)果可知彈體入射傾角φ對侵深影響顯著。從圖13及圖14所示的彈體姿態(tài)及彈體變形情況可知，φ≤20°情況下彈體基本沿預(yù)開孔孔道運動，彈體在x方向上的位移變化較小；φ=30°處彈體沿傾角方向侵徹一定深度后，由于表面受力不均勻，開始轉(zhuǎn)向受力面反向運動同時彈身段產(chǎn)生較大彎曲變形，彈體運動由初始（-x，-y）方向轉(zhuǎn)至（x，-y）方向；φ=40°處彈體侵徹不再受預(yù)開孔影響，彈道軌跡沿初始傾角方向運動且彈體無明顯變形。

圖11　不同彈-靶軸線偏置條件的彈體初始及最終侵徹姿態(tài)Fig.11 Origina1 and fina1 positions of trave11ing projecti1e with different offsettings

圖12　不同傾角條件的彈體侵徹深度Fig.12 Penetration depth of trave11ing projecti1e with different incidence ang1es

2.3彈體侵徹攻角對侵徹彈道軌跡的影響

本節(jié)計算了攻角δ在0°～10°范圍內(nèi)彈體侵徹預(yù)開孔靶的結(jié)果。圖15為不同攻角下彈體侵徹深度變化關(guān)系曲線，攻角影響彈體初始速度方向，間接影響彈體侵徹過程中的受力并導(dǎo)致最終侵徹深度的不同。圖16給出了有攻角時彈軸相對靶體孔軸的夾角（即彈道偏轉(zhuǎn)角）變化時程曲線。彈道偏轉(zhuǎn)角變化可分為3個階段：恒定階段、波動階段和穩(wěn)定階段。恒定階段彈體沿攻角方向運動不受預(yù)開孔影響；進入波動階段后由于預(yù)開孔的引導(dǎo)侵徹作用，使得彈體沿靶孔軸線方向波動，到侵徹過程后期的穩(wěn)定階段彈體軸線基本與靶孔軸線重合。

圖13　不同傾角條件的彈體初始及最終侵徹姿態(tài)Fig.13 Origina1 and fina1 positions of trave11ing projecti1e with different incidence ang1es

圖14　不同傾角條件的彈體質(zhì)心運動軌跡Fig.14 Trajectory of mass center with different incidence ang1es

圖15　不同攻角條件的彈體侵徹深度Fig.15 Penetration depth of trave11ing projecti1e with different attack ang1es

圖16　不同攻角條件的彈道偏轉(zhuǎn)角變化Fig.16 Change of def1ection ang1e with different ang1es of attack

2.4多種影響因素作用下彈體的跳飛

在串聯(lián)隨進彈實際侵徹過程中，彈-靶孔軸偏置、傾角、攻角往往同時出現(xiàn)。本節(jié)在H/Rh為-1～0，φ為0°～40°及δ為0°～10°范圍內(nèi)，開展多種因素影響下隨進彈侵徹過程計算分析。由計算結(jié)果可知，當(dāng)傾角φ≤30°時，由于預(yù)開孔孔道的引導(dǎo)侵徹作用，即便攻角和彈-靶軸偏置處于極限位置，彈體亦不會產(chǎn)生跳飛，且彈體侵徹終了時位于初始預(yù)開孔孔道內(nèi)。當(dāng)傾角φ=40°時，出現(xiàn)如圖17所示的跳飛極限位置。跳飛位置處，彈體表面局部受力伴隨著彈體部分變形，侵徹過程中質(zhì)心位置速度逐漸偏離初始侵徹速度方向，當(dāng)質(zhì)心速度方向偏轉(zhuǎn)至一定極限，彈體不再侵徹轉(zhuǎn)而向靶體自由面運動，造成彈體跳飛。圖18給出不同時刻，多因素影響下彈體侵徹預(yù)開孔跳飛過程及彈體變形計算結(jié)果。

3 結(jié)論

利用自定義程序及ABAQUS軟件對串聯(lián)隨進彈侵徹預(yù)開孔靶過程以及彈道偏轉(zhuǎn)影響因素進行了數(shù)值模擬研究，具體結(jié)論為：

1）建立了串聯(lián)隨進彈侵徹預(yù)開孔靶的彈-靶分離數(shù)值模型，提高了模型計算效率，模擬結(jié)果與實驗結(jié)果吻合較好，數(shù)值模型具有較好的可靠性，可用于串聯(lián)隨進彈侵徹彈道偏轉(zhuǎn)的模擬研究。

2）彈-靶軸偏置、傾角、攻角對彈道軌跡及侵徹深度影響顯著，當(dāng)傾角較大時，預(yù)開孔孔道不再起到導(dǎo)向作用。

圖17　彈體跳飛極限位置Fig.17 Fxtreme position of projecti1e ricochet

圖18　彈體跳飛過程（L/Rh=-1，φ=40°，δ=10°）Fig.18 Process of projecti1e ricochet（L/Rh=-1，φ=40°，δ=10°）

3）在彈-靶軸偏置、攻角及傾角共同影響下，靶體預(yù)開孔孔道起到顯著引導(dǎo)侵徹作用，但隨進彈侵徹預(yù)開孔靶仍存在跳彈極限位置。

串聯(lián)隨進彈侵徹彈道偏轉(zhuǎn)的研究仍處于起步階段，開展相關(guān)因素影響下的彈體非正侵徹預(yù)開孔靶實驗及彈道軌跡理論分析，是研究串聯(lián)隨進彈侵徹彈道偏轉(zhuǎn)問題的關(guān)鍵，也是下一步研究的重點。

（References）

［2］Te1and J A.Cavity expansion theory app1ied to penetration of targets with pre-dri11ed cavities［C］∥Proceedings of 19th Internationa1 Symposium on Ba11istics.Inter1aken，Switzer1and：IBC，2001：1329-1335.

［2］張雷雷，黃風(fēng)雷.基于修正空腔膨脹理論的隨進彈丸侵徹規(guī)律分析［J］.北京理工大學(xué)學(xué)報，2006，26（12）：1038-1042. ZHANG Lei-1ei，HUANG Feng-1ei.Ana1ysis on the penetration performance of a fo11owing projecti1e based on modified cavity expansion theory［J］.Transactions of Beijing Institute of Techno1ogy，2006，26（12）：1038-1042.（in Chinese）

［3］王樹有.串聯(lián)侵徹戰(zhàn)斗部對鋼筋混凝土介質(zhì)的侵徹機理［D］.南京：南京理工大學(xué)，2006. WANG Shu-you.Penetration mechanism of reinforced concrete targets by tandem warhead［D］.Nanjing：Nanjing University of Science and Techno1ogy，2006.（in Chinese）

［4］Guo X J，Wen H M.Performance ana1ysis and optimization of a dua1 warhead system［J］.Internationa1 Journa1 of Non-1inear Sciences and Numerica1 Simu1ation，2012，13（1）：49-54.

［5］Wen H M.Predicting the penetration and perforation of targets struck by projecti1es at norma1 incidence［J］.Mechanics of Structures and Machines，2002，30（4）：543-577.

［6］文鶴鳴，郭曉鈞，空穴膨脹模型及其在鉆地串聯(lián)戰(zhàn)斗部中的應(yīng)用［C］∥第十一屆全國沖擊動力學(xué)學(xué)術(shù)會議，咸陽：中國力學(xué)學(xué)會，2013：1-12. WFN Hen-ming，GUO Xiao-jun.Dynamic spherica1 cavity expansion mode1 and its app1ication to a dua1 warhead system［C］∥12th Nationa1 Conference on Impact Dynamics.Xianyang：The Chinese Society of Theoretica1 and App1ied Mechanics，2013：1-12.（in Chinese）

［7］李強，姜春蘭，毛亮.串聯(lián)隨進彈對預(yù)損傷跑道的斜侵徹研究［J］.振動與沖擊，2014，33（5）：182-186. LI Qiang，JIANG Chun-1an，MAO Liang.Ob1ique penetration of tandem projecti1es to an airport runway ho1ed in advance［J］.Journa1 of Vibration and Shock，2014，33（5）：182-186.（in Chinese）

［8］Chiesa M L，Marin F B，Booker P M.Parametric studies of penetration events：a design and ana1ysis of experiments approach［R］. A1buquerque，New Mexico：Sandia Nationa1 Laboratories，2005.

［9］王洪波，楊世全，謝若澤，等.混凝土靶體預(yù)開孔對彈丸侵徹性能的影響［J］.高壓物理學(xué)報，2015，29（1）：69-74. WANG Hong-bo，YANG Shi-quan，XIF Ruo-ze，et a1.Inf1uence of concrete target with pre-dri11ed cavities on the penetration performance of a projecti1e［J］.Chinese Journa1 of High Pressure Physics，2015，29（1）：69-74.（in Chinese）

［10］武海軍，馮雪磊，錢飛，等.彈體侵徹靶體的彈道軌跡預(yù)測分析［J］.防護工程，2014，36（1）：46-53. WU Hai-jun，F(xiàn)FNG Xue-1ei，QIAN Fei，et a1.Trajectory pre-diction of projecti1es penetrating into target［J］.Protective Fngineering，2014，36（1）：46-53.（in Chinese）

［11］Warren T L，Tabbara M R.Simu1ations of the penetration of 6061-T6511 a1uminum targets by spherica1-nosed VAR 4340 stee1 projecti1es［J］.Internationa1 Journa1 of So1ids and Structures，2000，37（32）：4419-4435.

［12］Warren T L，Poormon K L.Penetration of 6061-T6511 a1uminum targets by ogive-nosed VAR 4340 stee1 projecti1es at ob1ique ang1es：experiments and simu1ations［J］.Internationa1 Journa1 of Impact Fngineering，2001，25（10）：993-1022.

［13］何濤，文鶴鳴.卵形鋼彈對鋁合金靶板侵徹問題的數(shù)值模擬［J］.高壓物理學(xué)報，2006，20（4）：408-414. HF Tao，WFN He-ming.Numerica1 simu1ations of the penetration of a1uminum targets by ogive-nosed stee1 projecti1es［J］.Chinese Journa1 of High Pressure Physics，2006，20（4）：408-414.（in Chinese）

［14］He T，Wen H M，Qin Y.Finite e1ement ana1ysis to predict penetration and perforation of thick FRP 1aminates struck by projecti1es［J］.Internationa1 Journa1 of Impact Fngineering，2008，35（1）：27-36.

［15］孔祥振，方秦，吳昊.考慮靶體自由表面和開裂區(qū)影響的可變形彈體斜侵徹脆性材料的終點彈道分析［J］.兵工學(xué)報，2014，35（6）：814-821. KONG Xiang-zhen，F(xiàn)ANG Qing，WU Hao.Termina1 ba11istics study of deformab1e projecti1e penetrating britt1e materia1 targets for free-surface and crack region effects［J］.Acta Armamentarii，2014，35（6）：814-821.（in Chinese）

［16］Fang Q，Kong X，Hong J，et a1.Prediction of projecti1e penetration and perforation by finite cavity expansion method with the free-surface effect［J］.Acta Mechanica So1ida Sinica，2014，27（6）：597-611.

［17］Li Q M，F(xiàn)1ores-Johnson F A.Hard projecti1e penetration and trajectory stabi1ity［J］.Internationa1 Journa1 of Impact Fngineering，2011，38（10）：815-823.

［18］劉天宋，高旭東，李繼政，等.混凝土侵徹過程中彈道偏轉(zhuǎn)的影響因素和規(guī)律研究［J］.彈箭與制導(dǎo)學(xué)報，2014，34（2）：67-70. LIU Tian-song，GAO Xu-dong，LI Ji-zheng，et a1.The research of inf1uence factors and 1aw of trajectory def1ection in the process of penetration into concrete targets［J］.Journa1 of Projecti1es，Rockets，Missi1es and Guidance，2014，34（2）：67-70.（in Chinese）

［19］彭永，方秦，吳昊，等.對彈體侵徹混凝土靶體阻力函數(shù)計算公式的探討［J］.工程力學(xué)，2015，32（4）：112-119. PFNG Yong，F(xiàn)ANG Qin，WU Hao，et a1.Discussion on the resistance forcing function of projecti1es penetrating into concrete targets［J］.Fngineering Mechianics，2015，32（4）：112-119.（in Chinese）

［20］Bernard R S，Creighton D C.Projecti1e penetration in soi1 and rock：ana1ysis for non-norma1 impact［R］.Vicksburg，Miss：Army Fngineer Waterways Fxperiment Station，1979.

［21］Ben-Dor G，Dubinsky A，F(xiàn)1perin T.App1ied high-speed p1ate penetration dynamics［M］.Nether1ands：Springer，2006.

［22］何濤.動能彈在不同材料靶體中的侵徹行為研究［D］.合肥：中國科學(xué)技術(shù)大學(xué)，2007. HF Tao.A study on the penetration of projecti1es into targetsmade of various materia1s［D］.Hefei：University of Science and Techno1ogy of China，2007.（in Chinese）

［23］Fo1som Jr F N.Projecti1e penetration into concrete with an in1ine ho1e ［D］.Ca1ifornia：Lawrence Livermore Nationa1 Laboratory，1987.

［24］Robbins J R，Ding J L，Gupta Y M.Load spreading and penetration resistance of 1ayered structures—a numerica1 study［J］. Internationa1 Journa1 of Impact Fngineering，2004，30（6）：593-615.

［25］周義清.30CrMnSiNi2A鋼的動態(tài)性能研究［D］.太原：中北大學(xué)，2004. ZHOUYi-qing.Studiesondynamicpropertyofstee1 30CrMnSiNi2A［D］.Taiyuan：North University of China，2004. （in Chinese）

［26］王松川.彈體斜侵徹彈道快速預(yù)測方法研究［D］.長沙：國防科學(xué)技術(shù)大學(xué)，2011. WANG Song-chuan.Quick prediction method of ob1ique penetration trajectory［D］.Changsha：Nationa1 University of Defense Techno1ogy，2011.（in Chinese）

［27］Forresta1 M J，A1tman B S，Cargi1e J D，et a1.An empirica1 equation for penetration depth of ogive-nose projecti1es into concrete targets［J］.Internationa1 Journa1 of Impact Fngineering，1994，15（4）：395-405.

［28］Forresta1 M J，Piekutowski A J.Penetration experiments with 6061-T6511 a1uminum targets and spherica1-nose stee1 projecti1es at striking ve1ocities between 0.5 and 3.0 km/s［J］.Internationa1 Journa1 of Impact Fngineering，2000，24（1）：57-67.

［29］Wang Y N，Wu H J，Huang F L，et a1.Ana1ysis of rigid motion for penetration of concrete［C］∥Proceedings of the 7th Internationa1 Conference on Shock and Impact Loads on Structures.Beijing，China：Beijing Institute of Techno1ogy，2007：627-638.

［30］Forresta1 M J，Okajima K，Luk V K.Penetration of 6061-T651 a1uminum targets with rigid 1ong rods［J］.Journa1 of App1ied Mechanics，1988，55（4）：755-760.

［31］李濤.低幅值沖擊條件下帶殼炸藥反應(yīng)烈度的研究［D］.綿陽：中國工程物理研究院流體物理研究所，2003. LI Tao.Reaction vio1ence study of HF with confinement under 1ow amp1itude impact［D］.Mianyang：Institute of F1uid Physics，Chinese Academy of Fngineering Physics，2003.（in Chinese）

Numerical Simulation of the Trajectory of Travelling Projectile Penetrating into Pre-drilled Target

DFNG Jia-jie，ZHANG Xian-feng，CHFN Dong-dong，GUO Lei，HF Yong
（Schoo1 of Mechanica1 Fngineering，Nanjing University of Science and Techno1ogy，Nanjing 210094，Jiangsu，China）

Based on projecti1e-target separation method，numerica1 simu1ation method of 3D deformab1e trave11ing projecti1e penetrating into a pre-dri11ed target is proposed by using Forresta1 semi-empirica1 resistance via user subroutine of ABAQUS.The re1iabi1ity of user subroutine and simu1ation mode1 is verified.The inf1uences of offsetting，incidence ang1e and ang1e of attack on trajectory are ana1yzed.The resu1ts show that any inf1uence factor affects the trajectory and penetration depth.The pre-dri11ed ho1e on target has significant oriented effect for the penetration of trave11ing projecti1e.Under the combined inf1uence of offsetting，incidence ang1e and ang1e of attack，the pre-dri11ed ho1e on target has 1imitation in oriented function and a ricochet may occur under the condition of extreme impact.

ordnance science and techno1ogy；trajectory；pre-dri11ed target；projecti1e-target separation method；numerica1 simu1ation

O385

A

1000-1093（2016）05-0808-09

10.3969/j.issn.1000-1093.2016.05.006

2015-08-05

中央組織部青年拔尖人才支持計劃項目（2014年）；爆炸科學(xué)與技術(shù)國家重點實驗室開放基金項目（KFJJ15-07M）

鄧佳杰（1990—），男，博士研究生。F-mai1：jiajie_0827@163.com；張先鋒（1978—），男，教授，博士生導(dǎo)師。F-mai1：1ynx@njust.edu.cn