破片飛散角動靜爆關聯的試驗研究*

申景田，蘇健軍，姬建榮，陳 君

(西安近代化學研究所， 西安 710065)

0 引言

目前許多學者對彈藥靜爆威力場進行了大量的理論和試驗研究，為彈藥威力評估工作奠定了良好的基礎[1-3]。對于實際戰斗部來說，往往是在運動過程中爆炸，因此對彈藥動爆威力評估更有實用意義。目前已有學者對彈藥動爆威力場進行了一定的研究工作。

蔣海燕等[4]用仿真的方法分析了炸藥運動過程中爆炸時空氣沖擊波超壓的分布規律。魯衛紅等[5-8]對彈藥動爆破片的運動規律進行了大量研究：發現破片動態飛散區的多值問題并提出自己獨到的見解；提出一種破片動態分散密度的理論計算方法；建立了一種動爆破片運動規律的數學模型等。郭銳等[9]研究了破片飛行速度在彈體坐標系和地面固定坐標系之間的轉換關系。張龍杰、廖莎莎等[10-11]認為破片動態飛散速度可由靜態飛散速度和彈藥飛行速度矢量疊加得到。

由上述學者的研究工作來看，目前對彈藥動爆威力場的研究方法主要是理論分析和仿真試驗，沒有試驗數據的支撐。動爆試驗難度大、花費高，為了節省人力物力財力并且準確地評估彈藥的動爆威力，有必要研究動靜爆關聯問題。

1 試驗方法及結果

破片是彈藥的重要毀傷元之一，研究破片空間分布的動靜爆關聯對評估某些彈藥的動爆威力場有實際意義。為此，我們對同一批生產的某型彈藥分別進行動靜爆試驗。試驗采用直線排列的鋼靶板測量破片數量和空間分布；采用高速攝影裝置記錄彈藥的起爆點位置和飛行速度；采用梳狀靶測量破片著靶的平均速度。

試驗結果顯示破片的空間分布很廣泛，距離爆心較遠的位置破片分布密度較小，如圖1所示，橫向表示靶板位置，縱坐標表示著靶破片個數。

圖1 破片空間分布示意圖

由圖1可見，該型彈藥破片大部分集中分布在一個扇形殺傷帶中。經統計，動靜爆破片殺傷帶中破片數量分別占總破片數的85%以上，為簡化問題，至研究殺傷帶中破片飛散角的動靜爆關聯。破片動靜爆殺傷帶示意圖如圖2所示。

圖2 動靜爆殺傷帶示意圖

根據破片的著靶位置和彈藥起爆點位置可推算出破片的飛散角θs和θd(破片飛行方向和彈軸法線的夾角，偏向頭部為正，偏向尾部為負，下標s和d分別表示靜爆和動爆)，進而可計算出破片殺傷帶的飛散角和飛散方向角。

破片測速裝置測得的是破片從起爆點到靶板的平均速度。試驗中彈藥靜爆破片著靶時的平均速度分布范圍為1 130～1 695 m/s，破片飛散角絕對值越大，平均速度越小。根據高速攝影設備采集到的數據分析得出該型彈藥動爆試驗飛行速度為333 m/s。

2 試驗結果分析

有學者認為彈藥動爆破片飛散速度可由靜爆破片飛散速度和彈藥飛行速度進行矢量疊加得到[9-11]，如式(1)所示。

vd=vs+vm

(1)

式中：vd為破片動態飛散速度，vs為破片靜態飛散初速，vm為彈藥飛行速度。這種計算戰斗部動爆破片飛散角和飛行速度的方法稱為矢量疊加法。

只考慮式(1)中速度的方向，不考慮速度大小。矢量疊加后破片的動態飛散角計算公式如式(2)所示。

(2)

式中：θt為理論計算得到的破片動態飛散角。動爆破片飛散角理論與試驗結果如圖3。根據破片飛散角數據計算出的破片殺傷帶飛散角和飛散方向角數據如表1所示。由式(2)計算破片動態飛散角時用靜爆破片平均速度代替破片初速。

圖3 破片飛散角理論和試驗結果對比

表1 破片殺傷帶飛散角及飛散方向角

獲取方式參 數飛散角/(°)飛散方向角/(°)靜爆試驗結果29.1-2.1矢量疊加法計算結果26.313.6動爆試驗結果26.916.1

從圖3可以看出，用矢量疊加理論得到的破片動態飛散角比試驗值小，并且破片殺傷帶中心線附近的破片飛散角理論值的誤差較大。用矢量疊加法計算的破片動態飛散角理論值誤差最大為7.6°，占動爆試驗破片殺傷帶飛散角的28.3%。

從表1可以看出，矢量疊加法計算的破片殺傷帶飛散角誤差為0.6°，占動爆試驗破片殺傷帶飛散角的2.2%。破片殺傷帶飛散方向角理論值誤差為2.5°，占動爆試驗破片殺傷帶飛散角的9.3%。

可見矢量疊加法計算破片殺傷帶飛散角時誤差較小，計算的破片殺傷帶飛散方向角和試驗結果有一定誤差，計算破片動態飛散角時與試驗結果相差較大。這3種矢量疊加法誤差共同作用會使理論計算的殺傷帶內破片更集中于分布在殺傷帶的后半區。如圖4所示。

除此之外，用式(2)計算破片動態飛散角時用破片著靶平均速度代替破片初速，會使計算的破片動態飛散角偏大，即矢量疊加法計算的動爆破片飛散角的誤差比上述更大。

圖4 殺傷帶內破片分布密度示意圖

3 修正矢量疊加法

從矢量疊加法的計算過程可以看出，該方法有幾條基本假設條件。首先，假設彈藥靜爆時破片立即以一定的初速度在未擾動空氣中運動，忽略破片的加速過程，這和破片的實際運動過程有偏差。王馬法、盧芳云等[12]研究表明，靜爆時爆轟波和彈藥殼體的相互作用會對殼體產生一定的加速作用，并且使殼體偏轉一個角度。隋樹元[13]、張龍杰[10]等的研究表明，靜爆彈藥爆炸后，殼體在爆轟波和爆轟產物的高壓作用下發生膨脹、變形、破裂，乃至破碎。殼體破碎后形成分散的破片，爆轟產物溢出并包圍破片。破片在爆轟產物作用下一直被加速，直到爆轟產物膨脹速度相對破片運動可以忽略為止。破片停止加速時，可視為戰斗部內彈道結束，但很難確切定出這個時間。之后，由于爆轟產物和沖擊波速度衰減迅速，而破片速度衰減較慢，破片將沖到爆轟產物和沖擊波的前面。可見破片加速過程很復雜，難以準確地定量描述，目前工程上多用古尼公式估算破片靜爆初速。

其次，矢量疊加法假設彈藥動爆后破片立即以一定的初速在未擾動空氣中飛行，忽略彈藥飛行速度對破片加速過程的影響。蔣海燕等[4]基于顯式動力學分析軟件AUTODYN，建立了運動裝藥空中爆炸的數值計算模型，認識到彈藥運動速度會使空氣沖擊波超壓周向分布不均勻。由沖擊波波后介質狀態的連續性可知，波后爆轟產物的壓力分布一定程度上也是不均勻的，而爆轟產物膨脹對破片有一定的加速作用，所以彈藥飛行速度會影響爆轟產物的膨脹規律進而對破片加速過程產生一定的影響。除此之外，矢量疊加法未考慮到彈藥運動速度對彈藥內爆轟波傳播及空氣沖擊波超壓分布等因素對破片加速過程和速度衰減過程的影響，而這些因素和破片的運動狀態息息相關。

根據試驗結果提出一種矢量疊加法的修正形式，使得該修正形式計算的破片動態飛散方向角更接近實際情況。對比表1動爆試驗和矢量疊加法計算得到的破片動態飛散角發現，矢量疊加法計算結果偏小。因此，用矢量疊加法計算破片殺傷帶動態飛散方向角時在彈藥飛行速度前乘以一個大于1的系數ξ，稱之為牽連速度修正系數，使修正后的破片飛散角計算結果增大，更接近動爆試驗結果，如式(3)所示。

(3)

式(3)即計算破片殺傷帶動態飛散方向角的修正矢量疊加法，θmt表示用修正矢量疊加法計算得到的破片動態飛散角，系數ξ通過試驗數據擬合得到。根據動爆試驗數據計算的破片牽連速度修正系數如表2所示。

從表2可以看出牽連速度修正系數ξ整體上會隨著破片靜態飛散角絕對值的增加而減小。破片殺傷帶邊緣處的破片牽連速度修正系數比殺傷帶中心的破片牽連速度修正系數小。牽連速度修正系數越大，表示矢量疊加法計算的破片動態飛散角與試驗結果相差越大。

表2 牽連速度修正系數

從矢量疊加法的兩個假設條件和式(3)可以看出，影響ξ大小的因素有4個：破片靜態飛散角及飛散速度、破片動態飛散角、彈藥動爆飛行速度。可以將ξ用θs,θd,vs,vm表示出來，如式(4)所示。

ξ=f(θs,θm,vs,vm)

(4)

對于同一型號同一批生產的彈藥來說，靜爆時破片飛散角和飛散速度相關，動爆時破片飛散角和破片靜態飛散角、彈藥飛行速度相關。式(4)中的4個自變量不獨立，可以只用θs和vm將ξ表示出來，如式(5)所示。

ξ=f(θs,vm)

(5)

將式(5)寫成式(6)的形式。

(6)

從動靜爆試驗數據中選取n組有代表性的數據(θs,vm,ξ)，代入式(6)中(截取等號右邊的前m項，m不大于n)計算出m個待定系數。

根據本次試驗情況選取式(6)等號右邊前5項(如果選式(6)前6項會因為系數矩陣第6列數值過大而接近于奇異矩陣，難以求解)和6組試驗數據(θs,vm,ξ)，組成五元二次方程組如式(7)所示。

AX=B

(7)

X=pinv(A)·B

(8)

pinv(A)表示系數矩陣A的最小范數最小二乘廣義逆矩陣。將計算結果X代入式(6)得式(9)。

(9)

將式(9)計算出的破片牽連速度修正系數代入式(3)計算出的破片動態飛散角與試驗結果對比如圖5所示。

圖5 破片飛散角修正理論和試驗結果對比

對比圖3和圖5發現修正矢量疊加法的誤差比矢量疊加法的小。經計算，矢量疊加法計算破片動態飛散角誤差的均值為5.58°，方差為5.28°；修正矢量疊加法誤差的均值為0.06°，方差為3.58°。可見修正矢量疊加法可極大地減少理論計算破片動態飛散角的誤差，計算結果雖然在試驗數據附近有一定的波動，但比矢量疊加法的波動幅度小。

修正矢量疊加法是在傳統矢量疊加法的基礎上以彈藥動靜爆試驗破片的主要空間分布范圍為基準來修正破片牽連速度，使破片飛散角的計算結果更接近試驗結果。但該方法還存在一些不足的地方，例如：未從基礎理論定量研究彈藥結構、爆轟波傳播、爆轟產物膨脹、破片速度衰減規律等因素在彈藥動爆情況下對破片飛散特性的影響，而是將這些因素對破片

動態飛散角的影響歸納為破片牽連速度的變化；該理論是否適用于不同類型、不同飛行速度的彈藥還需要進一步試驗研究來驗證；修正矢量疊加法只研究了破片飛散角的動靜爆關聯問題，未涉及到破片飛散速度的動靜爆關聯、破片和沖擊波耦合作用的動靜爆關聯等問題，因此，完善該理論還需要大量理論和試驗研究。

4 結論

通過分析動靜爆試驗數據和矢量疊加法得到以下結論。

1) 矢量疊加法忽略了部分影響破片飛散特性動靜爆關聯的因素，該理論不夠完善。用該方法描述破片飛散角的動靜爆關系有一定誤差。

2)在矢量疊加法的基礎上，提出一種修正矢量疊加法來描述動破片飛散角動靜爆關聯，使理論計算更接近試驗結果。