破片戰斗部有效破片數影響規律研究

謝 奕，尹建平，李旭東，楊 東，盧文杰,張兆杰

(1.中北大學 機電工程學院，太原 030051；2.河南北方紅陽機電有限公司，河南 南陽 473000；3.晉西工業集團有限責任公司，太原 030051)

1 引言

現代航空炸彈爆炸后會伴隨著破片的飛散，破片是戰場上最為常見的毀傷元之一[1]。相比較于常規破片戰斗部，預制破片戰斗部[2]可以通過控制破片的形狀、質量與數量產生各種種類的破片，以達到所預想的毀傷效果。破片戰斗部起爆方式可以影響破片的飛散效果，史志鑫等[3]通過數值模擬的方法分析了不同起爆方式對預制破片的初速和飛散方向角的影響。孫韜等[4]基于LS-DYNA軟件對比分析了預制破片戰斗部靜、動態起爆時的飛散規律。彈丸落角也會影響破片的飛散效果，洪豆等[5]結合理論和仿真計算得出了不同落角、落速和炸高情況下預制破片戰斗部的殺傷面積。劉武[6]，厲相寶[7]，龔柏林[8]，李翔宇[9]通過數值模擬研究了不同類型破片戰斗部破片飛散過程。

航空炸彈常配備“觸發+延時”復合引信[10]，一般碰撞目標后有一定的延期時間才會起爆彈丸，同時由于叢林戰場地面土質松軟，航彈空中轟炸地面目標時一般會出現侵徹地面一定深度后才發生爆炸的情況，這會導致侵入地面部分戰斗部飛散出的破片在土壤中侵徹，大大降低破片的毀傷效果。在這種情況下，除了破片初速、飛散方向角和殺傷面積等評判標準[11]，產生的有效破片數也是尤為重要一種指標，有效破片數可定義為除去侵徹地面的破片數，在地面上方能對目標造成殺傷的破片數目。關于現有對破片飛散的數值模擬研究，對于戰斗部侵入土質松軟地面條件下起爆的破片飛散情況的研究較少。在研究中運用AUTODYN軟件研究戰斗部侵徹地面深度對產生的有效破片數的影響；并分析同一侵徹深度下，落角和起爆方式對破片的飛散規律及有效破片數的影響。為破片戰斗部侵入地面條件下的威力評估提供參考依據。

2 仿真模型

2.1 模型建立

仿真模型由炸藥、預制破片、地面及空氣組成，預制破片戰斗部的仿真模型如圖1所示。

圖1 有限元模型示意圖

戰斗部裝藥為圓柱形，其直徑為122 mm，裝藥長140 mm；預制破片采用球形破片，其破片直徑為6 mm，球形破片環繞炸藥單層排列，每層56個破片，共20層，共計1 120個破片。為了方便觀察破片的飛散情況，將預制破片從下往上每層依次選取破片設置測點并編號，預制破片的排列方式與測點的選取如圖2所示。

圖2 預制破片排列方式和測點選取示意圖

2.2 仿真方案

針對無制導航空炸彈，為了研究深度對戰斗部產生的有效破片數影響及同一深度條件下起爆方式和落角對有效破片數的影響，建立不同侵徹深度、起爆方式和落角下的預制破片戰斗部有限元模型。無制導航空炸彈從空中自由落下時，航彈一般會在垂直或者大落角情況下侵徹地面，若在小落角情況下易發生跳彈，導致命中精度降低，也不會出現侵徹地面現象，故仿真選取60°、70°、80°及90°的大落角；同時在空氣阻力等因素下，戰斗部侵徹地面的深度也不恒定，但實際中要避免戰斗部侵徹過深，故選取侵徹深度分別為戰斗部徑長的1/10、1/5、3/10、2/5、1/2及3/5，以便研究侵徹深度對有效破片數的影響。具體仿真工況如表1。

表1 仿真工況

彈軸中點起爆、底部端面中心點起爆和底部端面多點環形起爆示意圖如圖3、圖4和圖5。因為當戰斗部有向下的速度時，破片也被賦予向下的速度，不利于戰斗部起爆后破片向空中飛散形成有效破片，所以戰斗部侵入地面后速度衰減為0時為最佳效果。現代航空炸彈中常用采用“觸發+延時復合”復合引信，航彈撞擊地面時引信被觸發，經過預先設置足夠長的延遲時間后航彈再被引爆[12]，以達到侵徹地面速度為0后的起爆的效果，故戰斗部都采用靜態起爆。數值模擬使用ALE流固耦合算法[13]，數值模擬不同工況下的仿真模型，分析數值模擬結果得出這些參數對破片飛散規律及有效破片數的影響。

圖3 彈軸中點起爆示意圖

圖4 底部端面中心點起爆示意圖

圖5 底部端面多點環形起爆示意圖

2.3 模型材料參數

有限元仿真模型中的材料都是直接由AUTODYN軟件材料庫選出[14]，其材料參數也是默認給出的。空氣采用理想氣體狀態方程描述:

P=(γ-1)ρe

式中：P為空氣壓力；γ為理想氣體絕熱指數；ρ為空氣密度；e為空氣的初始比內能，在mm-mg-ms單位制下取2.068×105J/m3。空氣參數如表2所示。

表2 空氣參數

炸藥選用B炸藥，密度為1.717 g/cm3，爆速8 080 m/s；炸藥爆轟產物的狀態方程采用JWL(Jones-Wilkins-Lee)方程描述：

式中：P為炸藥爆炸產生的壓力；V是相對體積；E為單位體積內能；A、B、R1、R2、ω為材料常數，JWL狀態方程參數如表3所示。

表3 COMP B的JWL狀態方程參數

地面材料選用沙土SAND，密度為2.641 g/cm3，采用Compaction狀態方程、Mo Granual強度模型、Hydro失效模型。預制破片材料采用鎢合金TUNG.ALLOY，密度為 17 g/cm3，采用Shock狀態方程、Johnson Cook強度模型、Plastic Strain失效模型。每種材料的模型參數見表4。

表4 材料模型參數

3 數值模擬結果分析

3.1 侵徹地面深度對有效破片數的影響

數值模擬預制破片戰斗部在工況1～6下的破片飛散過程，以工況5為例，破片的飛散過程速度云圖如圖6，不同侵徹深度下戰斗部起爆后150 μs時破片的飛散效果云圖如圖7所示。

圖6 侵入地面42 mm時不同時刻下破片的飛散過程云圖

圖7 不同侵徹深度下t=150 μs時破片的飛散效果云圖

由數值模擬結果可以看出，炸藥爆轟所產生的能量賦予破片動能，驅動破片向外飛散，戰斗部中部產生的破片速度最快，兩端產生的破片速度較低。炸藥爆炸產生的沖擊波作用于土壤會產生反射波阻礙侵入地面部分破片的飛散，同時侵入地面部分破片飛散過程中也會受到土壤的阻礙作用，所以底部破片的速度最低，導致底部破片在空間中的飛散有著明顯的收斂，隨著侵徹深度的增大，侵徹地面土壤的底部破片數目增多，地面上能對目標造成殺傷的有效破片數目減少。由圖6可以看出，戰斗部裝藥爆炸產生的沖擊波會先于破片破壞土壤，導致原本戰斗部侵徹地面形成的坑進一步被破壞使其孔徑變得更大，隨后被驅動的破片才會開始侵徹土壤，破片起始的飛散并沒有出現侵徹土壤的現象，侵入地面部分破片的起始飛散沒有受到土壤的阻礙。隨著破片的持續飛散開始侵徹土壤，此時底部破片的速度衰減的加快，破片飛散方向也會發生改變。隨著戰斗部侵徹地面的深度越深，位于地下部分戰斗部的裝藥量越大，起爆后炸出的土塊越多。

有效破片數計算公式為：

ω=n+m-z

式中：ω為有效破片數；n為戰斗部起爆前完全位于地面上方的破片數，與侵徹深度及落角有關；m為起始位于土壤破片中起爆后從土壤飛出的破片數；z為起始位于地面上方破片中起爆后飛入土壤的破片數。通過數值模擬可得到破片觀測點的豎直方向速度及位移，觀測起爆后土壤中有豎直向上速度的破片及地面上方有向下速度的破片的豎直位移，判斷破片與地面的位置關系，可統計出m、z。

計算得出戰斗部不同侵徹深度下產生的有效破片數及其占總破片數的百分比如表5所示，結果表明隨著侵徹深度的增大，產生的有效破片的數量呈減小趨勢。戰斗部侵徹深度和有效破片數目總體成線性反比關系，故采用最小二乘法擬合該戰斗部侵徹地面深度與有效破片數占比的線性關系，擬合曲線如圖8，由擬合曲線得出該破片戰斗部侵徹地面深度與其產生的有效破片數占比的關系式為：

圖8 不同侵入深度下產生的有效破片占比及其擬合曲線

表5 不同深度下產生的有效破片數及其占總破片數的百分比

式中：y為有效破片數占比(%)；h為侵徹地面深度(mm)。

3.2 起爆方式對有效破片數影響

起爆方式對戰斗部有效破片數目的影響主要體現在對破片飛散方向角的影響[15]，破片的飛散方向角決定著地面附近破片是否侵入或飛出地面。破片飛散方向角定義為彈丸彈軸與破片飛散方向的夾角[16]，如圖9中θ所示。通過數值計算得到的測點破片速度和其在彈軸方向分量可以計算出破片的飛散方向角[17]，破片飛散方向角計算公式為：

圖9 破片飛散方向角示意圖

式中:θ為破片飛散方向角(°)；vz為破片速度軸線方向分量(m/s)；v為破片速度(m/s)。

對比分析起爆方式為彈軸中點起爆、底部端面中心點起爆和底部端面多點環形起爆的3種工況下的數值模擬結果，結果表明3種起爆方式所產生的有效破片數目大致相同，起爆方式的微小影響主要體現在對破片的飛散方向角的影響；破片向空中或地下飛散的趨勢會影響產生的有效破片數目，破片飛散方向角為90°時表示破片水平飛散，大于90°時破片有偏向地面飛散的趨勢，小于90°時破片有偏向地面上方空中飛散的趨勢，3種起爆方式下每層破片的飛散方向角如圖10所示。采用彈軸中點起爆、底部端面中心點起爆和底部端面多點環形起爆所產生破片的平均飛散方向角分別為89.3°、87.7°和89.5°，采用底部端面中心點起爆后破片平均飛散方向角最小，表明破片平均向上飛散程度相比其他起爆方式更大。

圖10 不同起爆方式下破片的飛散方向角曲線

戰斗部起爆前在28 mm侵徹深度下編號為3、4、5、6，7和8的破片位于地面附近，3種起爆方式下處于地面附近戰斗部破片飛散方向角如表6。底部端面中心點起爆時編號為3、4、5、6和8的破片飛散方向角分別為86.3°、88.5°、93°、88.5°、87.8°、78.9°，其破片飛散方向角相比其他2種起爆方式較小，且只有編號為5的破片飛散方向角大于90°，而其他2種起爆方式下產生的飛散方向角大于90°的破片多于底部端面中心點起爆。表明采用底部端面中心點起爆時地面附近破片僅有一層有侵徹地面趨勢，而其他層數破片都有向上飛散的趨勢；整體向上飛散趨勢相對其他2種起爆方式較為明顯，更有利于起爆前處于地面土壤中及處于地表部分的破片向上飛散形成有效破片。

表6 3種起爆方式下起爆后地面附近破片的飛散方向角

3.3 落角對有效破片數的影響

數值模擬預制破片戰斗部在工況9、10和11下的破片飛散過程，對比分析不同落角下產生的有效破片數，圖11為不同落角下的破片飛散效果云圖。將戰斗部在同一侵徹深度、不同落角下侵徹形成的地面模型導入SOLIDWORKS軟件，通過軟件可計算出戰斗部侵入地面所形成坑的體積即戰斗部侵入地面部分的體積，60°、70°、80°和90°落角下體積分別為578.20 cm3、664.23 cm3、688.28 cm3、818.29 cm3；由此在同一戰斗部侵徹深度條件下，即地面水平面到戰斗部侵入地面最低點的豎直距離相同時，落角的改變導致侵入地面的戰斗部體積有所變化。由圖11可以看出，隨著落角的減小，炸出的土塊越少；由于落角的存在，向地面傾斜的破片有些會侵徹地面，地面阻礙了向地面傾斜部分破片的飛散，此處地面會聚集許多破片，其密集度高，而侵入地面的破片有些會飛出地面形成有效破片的現象。表7為不同落角下產生有效破片數及其占總破片數的百分比，由表可得當落角在60°～90°范圍內，同一侵徹深度情況下隨著落角的減小，戰斗部形成的有效破片數目越多。70°、80°和90°落角下有效破片數差別較小，60°落角下有效破片數差別較大。

圖11 不同落角下的破片飛散效果云圖

表7 不同落角下產生的有效破片數及其占總破片數的百分比

4 結論

1)戰斗部侵徹地面起爆后，先炸出更深更寬的坑，戰斗部侵入地面部分的破片會在空中飛散一段時間后，再出現侵徹地面土壤的現象，侵入地面部分破片的起始飛散沒有受到土壤的阻礙；有效破片數目受侵徹地面深度的影響程度很大，破片戰斗部侵徹地面深度越大，產生的有效破片數目越少。

2)起爆方式對戰斗部形成的有效破片數的影響很小，對于彈軸中點起爆、底部端面中心點起爆和底部端面多點環形起爆三種起爆方式，底部中心點起爆的破片平均飛散方向角最小，且地面附近的飛散方向角小于90°的破片最多，破片總體向空中飛散趨勢較大，相對容易形成有效破片。

3)落角對戰斗部形成的有效破片數影響較大，相同侵徹深度下，對于該口徑戰斗部60°、70°、80°及90°落角下產生的有效破片數占比依次為67.86%、58.21%、56.7%和52.05%，在60°落角下形成的有效破片多，采用底部端面中心點起爆且落角為60°時該口徑戰斗部的更容易形成有效破片。在同一侵徹深度下大落角條件下，相對小的落角所得到有效破片數多。