萬仁毅， 李德元

(沈陽工業(yè)大學 材料科學與工程學院， 沈陽 110870)

?

預控破片形成及侵靶行為數(shù)值模擬*

萬仁毅， 李德元

(沈陽工業(yè)大學 材料科學與工程學院， 沈陽 110870)

為了優(yōu)化預控破片彈的彈體預制槽加工方案，增強彈丸殺傷威力，提高彈體材料利用率，運用有限元軟件ANSYS模擬了局部彈體材料爆轟形成預控破片的過程，研究了預控破片分別以不同初速度和不同侵徹角侵徹三種靶板的行為，并對模型應力分布和破片速度衰減變化進行了對比分析.結(jié)果表明，預制槽附近首先發(fā)生破碎，形成了形狀規(guī)則的預控破片，破片的初速度、靶板材料和侵徹角都將顯著影響預控破片的侵徹能力.模擬結(jié)果可為預控破片彈的設計和優(yōu)化提供參考.

ANSYS軟件； 有限元； 侵徹； 靶板； 預控破片； 彈體材料； 數(shù)值模擬； 模型

爆炸型彈藥形成的破片存在自然破片、預制破片和預控破片三種形式.多數(shù)爆炸型彈藥和通用炸彈會在爆炸時形成自然破片，彈丸殼體在爆炸作用下發(fā)生破碎并形成自然破片后，會向四周高速飛散，這些自然破片形狀各異、尺寸不一，初速度也并不相同，而破片的形狀、尺寸和初速度會對殺傷效果產(chǎn)生較大影響[1].增加破片的動能是增加榴彈殺傷威力的主要方式之一，可以考慮通過降低破片質(zhì)量的方式增加破片數(shù)量[2].預控破片是指在炸藥爆炸后形成的爆轟波作用下，戰(zhàn)斗部殼體以特定的破碎形式形成的形狀規(guī)則的破片.常用的預控破片形成技術(shù)包括殼體刻槽、裝藥內(nèi)表面刻槽和殼體內(nèi)嵌金屬罩等.隨著計算機技術(shù)與數(shù)值仿真模擬技術(shù)的不斷發(fā)展，計算機仿真模擬成為解決這類問題的重要手段.由于預控破片的形成、飛散及其終點毀傷效應等都是非常復雜的物理、力學過程，理論分析比較困難，而數(shù)值仿真方法能夠模擬介質(zhì)間的高速碰撞、相互作用、變形相變和斷裂等復雜過程，從而獲得清晰、連續(xù)變化的物理圖像與各種物理量的變化規(guī)律，進而可以準確了解破片的形成過程[3].采用計算機模擬方法分析破片的形成過程，具有重要的現(xiàn)實意義，可以為預控破片彈的設計和優(yōu)化提供理論參考.

國外學者針對預控破片戰(zhàn)斗部方面的研究起步得較早，尤其是在破片侵徹方面取得了很多研究成果[4].Finnegan[5]等研究了在不同侵徹狀態(tài)下，戰(zhàn)斗部空間破片的質(zhì)量和速度分布規(guī)律.隨著研究的進一步深入，Joseph[6]等研究了質(zhì)量不等的自然破片撞擊裝甲復合靶板的過程，建立了自然破片侵徹速度、破片質(zhì)量與侵徹深度之間的關(guān)系.目前，國內(nèi)學者針對民用和軍用領域內(nèi)戰(zhàn)斗部侵徹、毀傷作用數(shù)值模擬方面的研究較多，取得的數(shù)值模擬結(jié)果可為各種結(jié)構(gòu)和新型彈丸的設計與研發(fā)提供理論參考[7-8].

1　實體模型的建立

截取彈丸圓柱部的中間段進行數(shù)值模擬分析，上下斷面均為軸向約束.由于圓柱形彈丸殼體具有對稱性，因此，采用自頂向下的方法建立四分之一實體模型.在實體模型中，圓柱形殼體的外徑為122 mm，壁厚為10 mm，高為92 mm.在圓柱形殼體外表面刻出方形槽，且槽寬為4 mm，槽深為2 mm，槽間隔為12 mm.為了提高數(shù)值模擬的精度，選用SOLID164三維實體單元進行建模，殼體材料選用50SiMnVB.由于戰(zhàn)斗部在爆炸過程中伴隨高溫、高壓和高應變率等現(xiàn)象，因此，可以釆用ANSYS/LS-DYNA軟件模擬金屬爆炸成型、彈道侵徹和沖擊過程.殼體和靶板材料的具體參數(shù)如表1所示.此外，試驗中選用HMX炸藥，采用高能炸藥材料模型和JWL狀態(tài)方程進行建模，且JWL狀態(tài)方程表達式為

(1)

式中：p為轟炸產(chǎn)物的壓力；V為轟炸產(chǎn)物的相對比容；ω、A、B、R1、R2為表征炸藥特性的相關(guān)參數(shù)；E為初始比內(nèi)能.JWL狀態(tài)方程各參數(shù)的具體數(shù)值如表2所示.

表1　模型殼體材料參數(shù)

表2　JWL狀態(tài)方程參數(shù)

靶板材料分別選用杉木板、均質(zhì)鋼板和LY12靶板，規(guī)格均為100 mm×100 mm×10 mm.彈塑性模型可以很好地描述材料的應力-應變曲線，因此，采用MAT-PLASTIC-KINEMATIC模型對靶板進行建模.建立破片侵徹靶板的實體模型后，為了使計算結(jié)果更加接近實際情況，劃分網(wǎng)格時破片和靶板的網(wǎng)格邊長均取為1 mm，并采用映射方式劃分網(wǎng)格.模型建立后，需要設定分析選項.在模型對稱面施加法向的對稱約束，靶板邊界處施加反射界面.采用單面侵蝕接觸算法分析破片和靶板之間的關(guān)系.建立實體模型后，采用拉格朗日算法對殼體材料和炸藥進行仿真.劃分網(wǎng)格時殼體材料和炸藥網(wǎng)格邊長也均取為1 mm.采用映射方式對殼體進行網(wǎng)格劃分，采用掃略方式對炸藥進行網(wǎng)格劃分.在建立的實體模型的對稱面上施加約束，殼體刻槽后的彈丸戰(zhàn)斗部與預控破片侵靶模型如圖1所示.

實體模型完成網(wǎng)格劃分后，需要設置分析選項.首先，需要設定能量選項，從而方便考察炸藥面與殼體表面在接觸碰撞過程中的能量傳遞和變化過程，以及炸藥爆炸過程中炸藥面與殼體表面接觸時的能量耗損過程.然后，需要設置求解時間和時間步，將時長同步比例因子設置為0.6，輸出結(jié)果的時間步大小設置為0.5 μs.同時，需要在模擬計算過程中，保證時間步的間隔大小與數(shù)值解的精度要求相匹配.分析選項設置完成后生成了k文件，需要修改k文件中的材料模型和狀態(tài)方程關(guān)鍵字，注意材料號和狀態(tài)方程號需要一一對應，并設置炸藥的起爆方式.當完成對k文件的參數(shù)設定后，調(diào)用LS-DYNA970求解器進行相應的數(shù)值求解.求解結(jié)束后，利用LS-PREPOST軟件調(diào)取計算結(jié)果，輸出應力、應變和速度曲線，同時對預控破片的速度隨時間的變化規(guī)律進行分析.利用通用后處理器POST1，觀察殼體和炸藥的變形、應力及應變狀態(tài).利用時間歷程后處理器POST2繪制時間歷程曲線.

圖1　彈丸戰(zhàn)斗部和預控破片侵靶模型Fig.1　Models for projectile warhead and penetration target of pre-formed fragments

2　分析與討論

利用LS-PREPOST軟件調(diào)用實體模型的計算結(jié)果后，可以動態(tài)演示預控破片的形成過程.圖2為在預控破片形成過程中，不同時刻下殼體所受的應力分布.由圖2可見，當炸藥采用中心起爆方式時，殼體由刻槽的位置開裂破碎，形成了形狀規(guī)則的破片.由于炸藥從中心開始起爆，爆炸后產(chǎn)生的爆轟波會向殼體兩端傳播，在殼體上下兩個端面處存在氣體大量溢出的現(xiàn)象，因此，殼體中間部分破碎后產(chǎn)生的破片速度明顯大于殼體兩端.

利用圖2中的試驗數(shù)據(jù)，建立當具有規(guī)則形狀的預控破片以1 100 m/s的初速度垂直侵徹杉木板時的應力分布模型，并進行相應的數(shù)值模擬.侵徹杉木板時破片的應力分布和速度變化曲線如圖3所示.由圖3可見，在開坑過程中破片消耗的能量較少，破片剩余速度較大.當貫穿靶板后，破片剩余速度達到1 040 m/s，超過了其初速度的94.5%，即預控破片穿透10 mm厚的杉木靶板后仍具備殺傷威力.

當破片以500～1 100 m/s的初速度對LY12靶板進行侵徹時，侵徹鋁板時破片的速度變化與剩余速度-初速度關(guān)系曲線如圖4所示.由圖4可見，當破片侵徹10 mm厚的LY12靶板時，穿透靶板的破片初速度極限為550 m/s；當破片初速度大于極限穿靶速度后，破片的剩余速度隨著其初速度的增大近似呈線性關(guān)系增大.

圖2　不同時刻下殼體的應力分布Fig.2　Stress distribution of shell at different time

當模擬預控破片侵徹均質(zhì)鋼板時，預控破片的初速度為1 100 m/s，侵徹角分別為0°(垂直入射)、15°、30°和45°，侵徹鋼板時破片的速度變化與剩余速度-侵徹角關(guān)系曲線如圖5所示.由圖5可見，當侵徹角處于0°～30°范圍內(nèi)時，隨著侵徹角的增大，破片剩余速度逐漸減少，且侵徹角越大，破片的剩余速度衰減程度越大.這是由于當侵徹角增大時，破片在靶板水平方向的分速度增大，破片斜侵徹距離也隨之增大，因而消耗了大量動能，導致破片剩余速度降低.當侵徹角大于30°時，破片的剩余速度逐漸增大.這是由于此時破片并沒有穿透靶板，而是在斜侵徹一段距離后從靶板的右上方跳飛.

圖3　侵徹杉木板時破片的應力分布和速度變化曲線Fig.3　Stress distribution and velocity change curve of fragments when penetrating Chinese fir board

圖4侵徹鋁板時破片的速度變化與剩余速度-初速度關(guān)系曲線

Fig.4Curves for velocity change and residual velocity-initial velocity of fragments when penetrating aluminum board

圖5侵徹鋼板時破片的速度變化與剩余速度-侵徹角關(guān)系曲線

Fig.5Curves for velocity change and residual velocity-penetration angle of fragments when penetrating steel board

當預控破片以45°侵徹角對鋼板進行侵徹時，侵徹鋼板不同時刻的應力分布如圖6所示.由圖6可見，當侵徹68 μs后，破片與靶板脫離接觸，即發(fā)生了跳飛現(xiàn)象(見圖6d).當預控破片以1 100 m/s的初速度垂直侵徹不同靶板時，獲得的模擬結(jié)果如表3所示.由表3可見，當預控破片以相同初速度垂直侵徹不同靶板時，破片速度衰減程度顯著不同，靶板屈服應力是影響破片速度衰減的主要因素.

3　結(jié)　論

通過以上試驗分析，可以得出如下結(jié)論.

1) 當預控破片以相同初速度垂直侵徹不同靶板時，破片速度衰減程度顯著不同.

圖6　侵徹鋼板不同時刻的應力分布Fig.6　Stress distribution of penetrated steel board at different time

2) 當破片初速度低于550 m/s時，破片不能穿透LY12靶板；當初速度大于極限穿靶速度后，破片的剩余速度隨其初速度的增大近似呈線性關(guān)系增大.

3) 當預控破片垂直侵徹均質(zhì)鋼靶時，速度衰減很慢；當侵徹角增大至30°時，速度衰減較快；當侵徹角增大至45°時，侵徹過程中出現(xiàn)跳飛現(xiàn)象.

表3　預控破片垂直侵徹不同靶板的模擬結(jié)果

[1]吳廣，陳赟，馮順山，等.鋼靶運動對彈丸侵徹效應影響的仿真研究 [J].系統(tǒng)仿真學報，2012，24(2)：498-503.

(WU Guang，CHEN Bin，F(xiàn)ENG Shun-shan，et al.Simulation study of projectile normal penetration lateral movement target [J].Journal of System Simulation，2012，24(2)：498-503.)

[2]武敬博，茍瑞君，王聰，等.基于偏心起爆的戰(zhàn)斗部形狀對破片密度的增益研究 [J].兵器材料科學與工程，2016，39(2)：107-110.

(WU Jing-bo，GOU Rui-jun，WANG Cong，et al.Enhancement of warhead shape on fragments density based on eccentric initiation [J].Ordnance Materials Science and Engineering，2016，39(2)：107-110.)

[3]王瑤.預制破片毀傷效能數(shù)值仿真研究 [D].太原：中北大學，2012：7-11.

(WANG Yao.Study on numerical simulation of damage efficiency of prefabricated fragment [D].Taiyuan：North Central University，2012：7-11.)

[4]黃經(jīng)偉.破片式戰(zhàn)斗部破片規(guī)律及破片對步兵戰(zhàn)車的侵徹作用研究 [D].南京：南京理工大學，2014：15-20.

(HUANG Jing-wei.Fragment warhead department rules of fragment and fragments of infantry fighting vehicle intrusion penetration effect [D].Nanjing：Nanjing University of Science & Technology，2014：15-20.)

[5]Finnegan S A，Schulz J C，Heimdahl O E R.Spatial fragment mass and velocity distributions for ordnance and ultra-ordnance speed impacts [J].International Journal of Impact Engineering，1990，10(1/2/3/4)：159-170.

[6]Joseph B，Clay J.Calculating fragment impact velocity from penetration data [J].International Journal of Impact Engineering，2010，37：530-536.

[7]劉俊，李文彬.柱形破片侵徹鋼靶數(shù)值模擬研究 [J].信息技術(shù)，2016，45(2)：117-120.

(LIU Jun，LI Wen-bin.Numerical simulation and study of impact of cylindrical fragment against armor plates [J].Information Technology，2016，45(2)：117-120.)

[8]楊璐，石旭武.鋼筋混凝土梁塑性損傷模型的數(shù)值模擬 [J].沈陽工業(yè)大學學報，2016，38(1)：97-101.

(YANG Lu，SHI Xu-wu.Numerical simulation of plastic damage model for reinforced concrete beams [J].Journal of Shenyang University of Technology，2016，38(1)：97-101.)

(責任編輯：尹淑英英文審校：尹淑英)

Numerical simulation for formation and target penetrating behavior of pre-formed fragments

WAN Ren-yi, LI De-yuan

(School of Materials Science and Engineering, Shenyang University of Technology, Shenyang 110870, China)

In order to optimize the machining scheme for the pre-formed projectile body grooves of pre-formed fragment projectile, enhance the lethal power of projectile body and improve the utilization rate of projectile body material, the detonation forming process of pre-formed fragments was simulated with finite element software ANSYS. In addition, the penetration behavior of pre-formed fragments on three kinds of target boards at different initial velocities and penetration angles was researched. The model stress distribution and the velocity attenuation of fragments were compared and analyzed. The results indicate that the vicinity of pre-formed grooves is firstly broken, and the pre-formed fragments with regular shape form. Moreover, the initial velocity of fragments, target board material and penetration angle will remarkably affect the penetration capability of pre-formed fragments. The simulated results can provide the

for the design and optimization of pre-formed fragment projectile.

ANSYS software; finite element; penetration; target board; pre-formed fragment; projectile body material; numerical simulation; model

2016-06-17.

遼寧省教育廳重點實驗室基礎研究資助項目(LZ2014013).

萬仁毅(1979-)，男，重慶人，博士生，主要從事彈體新材料與新工藝等方面的研究.

10.7688/j.issn.1000-1646.2016.05.06

TG 156.33

A

1000-1646(2016)05-0508-05

*本文已于2016-09-07 16∶10在中國知網(wǎng)優(yōu)先數(shù)字出版. 網(wǎng)絡出版地址：http：∥www.cnki.net/kcms/detail/21.1189.T.20160907.1610.058.html