張健東，武海軍，李 偉，李金柱，皮愛國

（北京理工大學(xué)爆炸科學(xué)與技術(shù)國家重點實驗室, 北京 100081）

在現(xiàn)代防御體系中，有限厚混凝土作為最常見的防御手段，在地面建筑和防御工事中應(yīng)用廣泛。彈體在對薄板類混凝土多層結(jié)構(gòu)進(jìn)行侵徹時，由于存在傾角、攻角等因素，彈體受力不對稱，彈道發(fā)生偏轉(zhuǎn)，導(dǎo)致彈體不能有效貫穿多層結(jié)構(gòu)，降低了彈體的有效侵徹能力。因此，研究彈體斜侵徹貫穿混凝土薄靶的作用機(jī)理和提高彈道穩(wěn)定性具有重要意義。

為了解決這些問題，研究人員對彈體侵徹混凝土靶的作用機(jī)理和彈道穩(wěn)定性等問題進(jìn)行了大量研究。Forrestal 等[1-2]基于空腔膨脹理論建立了卵形彈體侵徹混凝土的阻力模型，提出了Forrestal 半經(jīng)驗公式。Chen 等[3]基于空腔膨脹理論和沖塞破壞假定，針對剛性彈體斜侵徹有限厚混凝土靶的彈道偏轉(zhuǎn)問題，提出了包含開坑階段、隧道階段和剪切沖塞階段的三階段侵徹模型，但是該模型僅考慮了彈體在開坑階段的偏轉(zhuǎn)，而沒有考慮沖塞階段彈體姿態(tài)角減小的問題。馬兆芳等[4-6]根據(jù)剛性彈體斜貫穿混凝土靶的實驗結(jié)果，系統(tǒng)研究了彈體斜侵徹貫穿混凝土靶的偏轉(zhuǎn)過程，在Chen 等[3]的基礎(chǔ)上提出了沖塞階段彈體二次偏轉(zhuǎn)機(jī)制。另一方面，研究人員對不同彈形對侵徹能力的影響也進(jìn)行了研究。Kong 等[7]基于彈靶分離再接觸效應(yīng)、自由面效應(yīng)，研究了質(zhì)心位置、長徑比等對侵徹彈道的影響。鄧佳杰等[8]基于頭部非對稱刻槽彈體，建立了軸向壓-剪作用下的準(zhǔn)靜態(tài)柱形空腔膨脹理論并對此開展實驗研究，發(fā)現(xiàn)頭部非對稱刻槽彈體的自旋效應(yīng)能有效減小侵徹阻力，提高彈體的侵徹能力。Liu 等[9]基于空腔膨脹理論，分析了彈體頭部形狀系數(shù)與雙卵形特征參數(shù)的關(guān)系，得到了不同形狀系數(shù)對彈體侵深的影響規(guī)律。以色列軍事工業(yè)公司設(shè)計了多款頭部帶肋板的侵徹戰(zhàn)斗部[10]，并對不同肋板結(jié)構(gòu)對彈體侵徹能力和彈道穩(wěn)定性的影響進(jìn)行了深入研究，優(yōu)化后的彈體能在不明顯增加彈體質(zhì)量的前提下，提升彈體斜侵徹混凝土的能力，實驗中彈體能以10°傾角、375 m/s 初始速度無偏移貫穿3 層混凝土墻，如圖1 所示。綜上所述，目前對侵徹戰(zhàn)斗部的研究大多集中于彈體的侵徹能力以及彈體在半無限靶中的偏轉(zhuǎn)規(guī)律等方面，對新型彈體結(jié)構(gòu)設(shè)計和加強(qiáng)彈體斜貫穿混凝土薄靶時的彈道穩(wěn)定性研究相對較少。因此，為了提高彈體斜侵徹貫穿混凝土薄靶的彈道穩(wěn)定性，需要設(shè)計具有良好彈道性能的彈體結(jié)構(gòu)并結(jié)合實驗和數(shù)值模擬進(jìn)行研究。

圖1 戰(zhàn)斗部以10°傾角貫穿多層混凝土靶[10]Fig. 1 Projectile perforated the multilayer concrete targets at 10° inclination[10]

本研究設(shè)計并開展了頭部帶肋板的異形戰(zhàn)斗部在不同速度和傾角下貫穿混凝土薄靶實驗，通過高速攝影系統(tǒng)獲得彈體在貫穿靶板前后的速度和姿態(tài)變化，并根據(jù)實驗結(jié)果對侵徹過程和彈道穩(wěn)定性進(jìn)行分析。然后，利用動力有限元程序LSDYNA 對相同質(zhì)量、直徑的頭部帶肋板異形彈和尖卵形彈進(jìn)行數(shù)值模擬，分析兩種結(jié)構(gòu)彈體斜貫穿薄靶的作用過程，探究頭部帶肋板異形彈體的侵徹特性。

1 頭部帶肋板的異形彈體侵徹混凝土靶實驗

1.1 實驗方法

為了與常規(guī)尖卵形彈體進(jìn)行對比，以文獻(xiàn)[11]中的卵形彈為參考，在保證彈體質(zhì)量、長度和直徑相同的情況下，對彈體頭部進(jìn)行了重新設(shè)計，得到頭部帶肋板異形彈體，優(yōu)化后的彈體長度為317.0 mm、彈體直徑為64.0 mm。如圖2 所示，彈體頭部結(jié)構(gòu)由圓形凸臺和8 個肋板條構(gòu)成，其中凸臺頂面直徑為10.4 mm，凸臺底面直徑為21.6 mm，凸臺高度為18.0 mm；肋條外緣直徑為41.4 mm，肋板長度為63.0 mm，肋條最小寬度為4.5 mm，肋條高度為8.1 mm。彈體材料選用30CrMnSiNi2A，質(zhì)量為4.5 kg。彈體實物如圖3 所示，表面噴有白漆和黑色圓點，以便用高速攝影捕捉彈體運(yùn)動軌跡。實驗靶體為 ?1 200 mm×300 mm 的素混凝土靶，為了方便澆筑和保證實驗后靶體的完整性，外層用壁厚3 mm 的鋼箍固定，其澆筑和具體養(yǎng)護(hù)均按照標(biāo)準(zhǔn)規(guī)范執(zhí)行，實驗前測得靶體平均強(qiáng)度為38 MPa。

圖2 彈體頭部結(jié)構(gòu)示意圖（單位：mm）Fig. 2 Sketch of the projectile (Unit: mm)

圖3 彈體實物Fig. 3 Photograph of the projectile

該實驗在北京理工大學(xué)爆炸科學(xué)與技術(shù)國家重點實驗室的152 mm 一級輕氣炮發(fā)射平臺上進(jìn)行，實驗系統(tǒng)示意圖和實驗現(xiàn)場見圖4。炮口正對靶板中心，通過靶架和靶室軌道固定靶板，同時在靶板背面固定10 mm 厚的木板用以阻擋混凝土碎塊。為了回收彈體，在靶板后方放置一個較厚的混凝土塊。在靶室外側(cè)設(shè)置兩臺高速相機(jī)，分別記錄彈體入靶前和出靶后的姿態(tài)和速度，當(dāng)彈體經(jīng)過炮口時撞斷纏繞在炮口上的導(dǎo)線，觸發(fā)高速相機(jī)記錄彈靶作用過程。實驗結(jié)束后對靶體開坑尺寸、崩落區(qū)域大小等進(jìn)行測量。

圖4 實驗系統(tǒng)示意圖(a) 和實驗現(xiàn)場(b)Fig. 4 Sketch of the experimental system (a) and experimental site (b)

1.2 實驗結(jié)果及分析

共進(jìn)行5 發(fā)實驗，彈體速度在250～350 m/s 之間。由于輕氣炮靶室和觀測窗口大小有限，在高速攝影的視野內(nèi)無法安裝反光鏡對彈體在水平方向的姿態(tài)進(jìn)行觀察，所以只對彈體豎直方向偏轉(zhuǎn)的情況進(jìn)行討論。如圖5 所示，彈體以259.6 m/s 的速度、30°傾角斜侵徹時，觀察靶前和靶后高速攝影儀記錄的彈體運(yùn)動軌跡，彈體在豎直方向入靶姿態(tài)良好，沒有發(fā)生攻角侵徹，出靶之后彈體姿態(tài)清晰。實驗回收得到的彈體如圖6 所示，5 發(fā)彈體結(jié)構(gòu)均完整，僅有表面噴漆磨損。實驗彈道軌跡如圖7 所示，其中：β0為傾角（彈體軸線與靶體法線的夾角）， β為彈體姿態(tài)角（彈體運(yùn)動過程中彈體軸線與靶體法線的夾角）；為了明確彈體姿態(tài)偏轉(zhuǎn)程度，定義出靶角 β1為彈體出靶時彈體軸線與水平線的夾角（即姿態(tài)角的變化值），順時針為正、逆時針為負(fù)。具體實驗結(jié)果見表1。v0和v1分別為入靶速度和出靶速度，d1、d2分別為靶體正面開坑和背面崩落直徑。

圖5 彈體斜侵徹姿態(tài)Fig. 5 Oblique penetration attitude of the projectile

圖6 回收彈體Fig. 6 Residual projectiles

圖7 彈道軌跡示意圖Fig. 7 Attitude deflection of the projectile

表1 彈體貫穿混凝土靶實驗結(jié)果Table 1 Experimental results of projectiles penetrating concrete targets

由于無法保證每發(fā)彈體的發(fā)射速度完全一致，因此認(rèn)為發(fā)射速度近似的第1、2、4 發(fā)實驗為同一發(fā)射速度水平，第3、5 發(fā)實驗為同一發(fā)射速度水平。根據(jù)表1 的實驗數(shù)據(jù)，當(dāng)彈體發(fā)射速度在同一水平時，隨著傾角增大，彈體貫穿靶板后的剩余速度越低；當(dāng)傾角相同時，著靶速度降低，彈體偏轉(zhuǎn)越厲害。這是因為在斜侵徹過程中，彈體侵入靶體和沖出靶體時都會受到上下不對稱的偏轉(zhuǎn)力矩，導(dǎo)致彈道發(fā)生改變的同時降低了彈體的侵徹能力。

圖8 給出了5 發(fā)實驗靶體的開坑和崩落情況。從表1 和圖8 可以看到，與侵徹半無限靶不同，彈體對薄靶的貫穿不存在隧道區(qū)，靶體正面漏斗形的開坑區(qū)與靶體背部崩落區(qū)直接相連，在靶體中間形成一個約兩倍彈身直徑大小的孔洞。當(dāng)彈體侵入靶體時，靶體正面的混凝土在較高的沖擊壓力作用下破壞成碎塊，圍繞著靶點形成一個漏斗形的開坑，并產(chǎn)生數(shù)量不等的徑向裂紋。當(dāng)有傾角時，彈體開坑會造成靶體上下表面受力不均，使得開坑呈現(xiàn)扇形，同時裂紋也更明顯。當(dāng)彈體進(jìn)入靶體后，由于混凝土的抗拉和抗剪強(qiáng)度遠(yuǎn)低于抗壓強(qiáng)度，在臨近靶體背面，彈體撞擊靶板形成的球面壓縮波在離此點最近的自由面反射形成拉伸波，因而靶體內(nèi)部受到拉伸破壞的部位并不在彈體運(yùn)動路徑上[12]，造成了彈體出靶階段的二次偏轉(zhuǎn)。因此，彈體貫穿混凝土靶板造成的靶板正面開坑和背部崩落對彈體姿態(tài)的影響不可忽視。

圖8 5 發(fā)實驗靶體的正面(a) 和背面(b) 破壞情況Fig. 8 Damage condition of front (a) and back (b) of the experimental target

2 彈體侵徹過程數(shù)值模擬

2.1 有限元模型

為了研究頭部帶肋板異形彈體斜貫穿混凝土薄靶的能力和抗偏轉(zhuǎn)能力，在保證彈體的質(zhì)量、直徑、轉(zhuǎn)動慣量一致的前提下，利用TrueGrid 和LS-DYNA 建立頭部帶肋板異形結(jié)構(gòu)彈體和尖卵形頭部彈體的有限元模型，其中尖卵形頭部彈體參考文獻(xiàn)[11]，彈體結(jié)構(gòu)如圖9 所示。在建模過程中，為了消除邊界效應(yīng)對侵徹結(jié)果的影響，將靶體側(cè)面設(shè)置為無反射邊界，同時為了模擬靶架的約束，將靶體側(cè)面邊界設(shè)置為固定約束，靶體迎彈面和靶體背面保持為自由面。為了在保證計算結(jié)果可靠的同時縮短計算時長，對靶板中間10 倍彈徑的區(qū)域進(jìn)行網(wǎng)格加密。有限元模型（圖10）尺寸與實驗保持一致。

圖9 尖卵形彈體結(jié)構(gòu)示意圖（單位：mm）Fig. 9 Sketches of the ogive-nose projectile (Unit: mm)

圖10 有限元模型Fig. 10 Finite element model

2.2 材料模型及參數(shù)

實驗后回收的彈體結(jié)構(gòu)完整，沒有發(fā)生明顯侵蝕。同時根據(jù)Frew 等[13-14]的研究結(jié)果可知，中低速侵徹過程中彈體基本保持剛性特征。因此，數(shù)值模擬中彈體材料模型選用剛性模型，密度為7.85 g/cm3。由于彈體貫穿混凝土薄靶過程中靶體的開坑和崩落對彈體姿態(tài)的影響不可忽略，對混凝土分別選用HJC 模型和自定義TCK 模型[5]進(jìn)行計算。HJC 模型是主要針對沖擊加載下混凝土材料發(fā)生大應(yīng)變失效開發(fā)的材料模型，同時兼顧了應(yīng)變率效應(yīng)和塑性體積應(yīng)變效應(yīng)的壓縮損傷失效，但是該模型缺乏對混凝土脆性拉伸損傷的考慮；而TCK 模型則對材料的脆性拉伸損傷和裂紋擴(kuò)展等有很好的考慮，能夠較好地模擬裂紋的擴(kuò)展和開坑。兩種材料模型的基本參數(shù)分別如表2 和表3 所示。表2 中， ρ為材料密度，G為剪切模量，A、B分別為內(nèi)聚力強(qiáng)度和壓力強(qiáng)化系數(shù)，C為應(yīng)變率敏感系數(shù)，N為壓力硬化系數(shù)，fc為單軸抗壓強(qiáng)度，T為最大抗拉強(qiáng)度，?f,min為混凝土破碎的最小塑性應(yīng)變，Smax為混凝土歸一化強(qiáng)度，pc和μc分別為混凝土單軸受壓破碎時的壓力和體積應(yīng)變，pL和UL分別為壓實時的壓力和體積應(yīng)變，D1、D2為損傷參數(shù)，K1、K2、K3為壓力常數(shù)。表3 中，E為楊氏模量，K為體積模量，ν為泊松比，k、m為材料常數(shù)，KIC為材料的斷裂韌性。

表2 HJC 材料模型參數(shù)Table 2 Parameters of HJC material model

表3 TCK 材料模型參數(shù)Table 3 Parameters of TCK material model

3 結(jié)果分析

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

表4 給出了HJC 模型和TCK 模型的數(shù)值模擬和實驗結(jié)果?？梢钥闯觯?dāng)混凝土材料采用HJC 模型時，不同工況下彈體剩余速度的相對偏差在10% 以下，彈體出靶角度誤差不超過5.7°；當(dāng)采用TCK 模型時，彈體剩余速度的相對偏差保持在10%以下，且彈體出靶角度與實驗結(jié)果的誤差不超過3.33°。說明TCK 模型在預(yù)測彈體偏轉(zhuǎn)方面更有優(yōu)勢。但是，對于第3 發(fā)和第5 發(fā)實驗，兩種模型中的彈體都沒能貫穿靶板，主要原因是低速貫穿時彈體受靶體的阻力影響更明顯，而數(shù)值模擬并沒有完全模擬出混凝土的崩落，導(dǎo)致出靶過程中彈體受到的阻力仍然較大，彈體無法完全貫穿靶板。

表4 實驗與數(shù)值模擬結(jié)果對比Table 4 Comparison of experimental and numerical simulation results

圖11 和圖12 分別給出了335.9 m/s、30°傾角侵徹時HJC 模型和TCK 模型下靶板的損傷情況。因為HJC 模型考慮材料的壓縮損傷累積效應(yīng)、圍壓效應(yīng)和軟化效應(yīng)較完善，所以損傷主要集中在彈體周圍的高壓區(qū)域；而選用TCK 模型時，在靶體背面區(qū)域可以看到由于自由面反射壓縮波形成的拉伸損傷，損傷區(qū)域并不在彈體運(yùn)動路徑上，而是彈體下方的損傷區(qū)域更大。在真實情況下，靶板背部受到拉伸損傷的混凝土?xí)苯颖缆?，崩落后靶體對彈體產(chǎn)生的約束不對稱，從而導(dǎo)致彈體出靶過程中會發(fā)生二次偏轉(zhuǎn)，彈道偏轉(zhuǎn)如圖13 所示。根據(jù)模擬結(jié)果，由于TCK 模型能更好地預(yù)估彈體姿態(tài)，所以在頭部帶肋板異形彈體的侵徹能力和彈道穩(wěn)定性彈數(shù)值模擬研究中采用了自定義TCK 模型。

圖11 HJC 模型壓縮損傷云圖Fig. 11 Compression damage obtained by HJC model

圖12 TCK 模型拉伸損傷云圖Fig. 12 Tensile damage obtained by TCK model

圖13 彈道偏轉(zhuǎn)Fig. 13 Ballistic trajectory

3.2 頭部形狀對彈體侵徹能力的影響

圖14 為不同工況下兩種彈體貫穿靶板的入靶速度與剩余速度的關(guān)系。當(dāng)頭部帶肋板異形彈和卵形彈以30°傾角、不同速度貫穿混凝土靶時，隨著入靶速度增大，剩余速度也隨之增大，并且表現(xiàn)出近似線性關(guān)系，其中卵形彈的剩余速度高于異形彈25 m/s 左右，表明卵形彈具有更強(qiáng)的侵徹能力。圖15為兩種彈體以450 m/s、30°傾角貫穿靶板時的過載曲線，其中頭部帶肋板異形彈的最大過載為2.5×104g，卵形彈最大過載為2.1×104g。在侵徹初期彈頭侵入靶體階段，由于該異形彈頭部為凸臺加肋板結(jié)構(gòu)，其與混凝土靶板的接觸面積比卵形彈更大，所以軸向過載也更大；在彈頭侵入靶體后，彈體的阻力主要與彈體橫截面積和彈體速度有關(guān)，隨著速度降低，過載變小；當(dāng)彈頭出靶后，阻力只有彈洞的側(cè)壁摩擦力，過載迅速減小。

圖14 不同工況下彈體的剩余速度Fig. 14 Residual velocities of projectiles under different cases

圖15 侵徹過載-時程曲線Fig. 15 Over-load versus time curves of projectiles

3.3 頭部形狀對彈道偏轉(zhuǎn)的影響

為了進(jìn)一步研究兩種彈體的彈道特性，對不同傾角和著靶速度下兩種彈體侵徹混凝土的過程進(jìn)行了數(shù)值模擬。圖16 為彈體以不同速度侵徹時的姿態(tài)角-時間變化曲線，圖17 為400 m/s 初始速度下以不同傾角侵徹時的姿態(tài)角-時間變化曲線?？梢钥闯?，在貫穿混凝土靶時，兩種彈體都經(jīng)歷了姿態(tài)角先增大再減小的過程。姿態(tài)角增大發(fā)生在入靶階段，姿態(tài)角減小發(fā)生在沖塞階段。隨著入靶速度降低和傾角增大，彈體入靶階段姿態(tài)角變化更加明顯。

圖16 以不同初始速度侵徹時彈體的姿態(tài)角-時間變化曲線Fig. 16 Attitude angle versus time curves of projectile under different initial velocities

圖17 不同傾角侵徹時的姿態(tài)角-時間變化曲線Fig. 17 Attitude angle curves of projectile under different inclination angles

彈體的角加速度與繞質(zhì)心的轉(zhuǎn)動慣量相乘得到侵徹過程中彈體的偏轉(zhuǎn)力矩，如圖18 所示。當(dāng)彈體以450 m/s 進(jìn)行侵徹時，頭部帶肋板異形彈在入靶階段受到的最大偏轉(zhuǎn)力矩為792 N·m，而卵形彈則為1 462 N·m。因為彈體侵入靶板時靠近靶板一側(cè)的肋板會先接觸混凝土，產(chǎn)生一個繞彈體質(zhì)心，抑制彈體偏轉(zhuǎn)的力矩，提升了彈體姿態(tài)的穩(wěn)定性，所以頭部帶肋板異形彈姿態(tài)角變化幅度很小。在沖塞階段，由于靶板背面發(fā)生崩落，且崩落區(qū)域不在彈體運(yùn)動路徑上，此時彈體受到靶板的不對稱約束力將會使彈體姿態(tài)角減小。如圖18 所示，當(dāng)彈體以450 m/s 進(jìn)行侵徹時，0.75 ms 之后兩種彈體的偏轉(zhuǎn)力矩曲線基本重合，同樣說明沖塞階段彈體偏轉(zhuǎn)與彈體頭部形狀無關(guān)，只與彈身受到的不對稱約束有關(guān)。

圖18 彈體偏轉(zhuǎn)力矩變化曲線Fig. 18 Deflection moment curves of projectiles

4 結(jié) 論

開展了頭部帶肋板異形結(jié)構(gòu)彈體貫穿混凝土靶板的實驗研究，并根據(jù)實驗工況對頭部帶肋板異形彈體和尖卵形頭部彈體開展了不同傾角和速度下貫穿混凝土靶的數(shù)值模擬，得出以下結(jié)論。

(1) 彈體貫穿混凝土薄靶時，靶體存在正面開坑區(qū)和背面沖塞區(qū)而沒有隧道段。在250～360 m/s的實驗速度區(qū)間，傾角增大，會明顯降低彈體的貫穿能力和彈道穩(wěn)定性，傾角越大，剩余速度越低，彈體姿態(tài)偏轉(zhuǎn)越明顯。

(2) 彈體貫穿混凝土薄靶時，開坑階段的彈體受力不對稱會增大彈體的姿態(tài)角，而拉伸損傷引起的背部層裂崩落會減小彈體的姿態(tài)角，在研究貫穿機(jī)理時應(yīng)當(dāng)綜合考慮混凝土的拉伸和壓縮破壞對彈體姿態(tài)的影響。

(3) 相比尖卵形彈體，由于頭部帶肋板異形結(jié)構(gòu)彈體肋板和凸臺結(jié)構(gòu)的阻礙作用更強(qiáng)，侵徹過程中的過載更大，彈體貫穿靶板后的剩余速度更小。當(dāng)以相同速度斜貫穿靶板時，頭部帶肋板彈在侵入靶板時受到的偏轉(zhuǎn)力矩更小，彈道穩(wěn)定性更好。