基于戰斗部動爆破片特性分析的雷達回波模型

侯建強，韓壯志，彭　剛，李新欣，黃　巍

(1.解放軍軍械工程學院電子與光學工程系,河北 石家莊　050003；2.解放軍77618部隊,西藏

拉薩　850000；3.中國電子科技集團公司第五十四所,河北 石家莊　050003)

?

基于戰斗部動爆破片特性分析的雷達回波模型

侯建強1，韓壯志1，彭剛2，李新欣3，黃巍3

(1.解放軍軍械工程學院電子與光學工程系,河北 石家莊050003；2.解放軍77618部隊,西藏

拉薩850000；3.中國電子科技集團公司第五十四所,河北 石家莊050003)

0引言

彈丸爆炸后，破片以一定的方向飛散出去，其在空間的分布情況是影響破片殺傷作用場的重要因素[1]。戰斗部正常破片的初速以及初速分布是武器總體計算引戰配合和殺傷概率時的必需參數[2]。通常，武器系統戰斗部設計、試驗中的破片威力試驗，一般采用靜爆法，即將戰斗部置于一定高度的托彈架上，根據戰斗部的威力大小和測試項目布置靶板、靶網、傳感器及各種測試儀器[3]，得到所需參數。由于受到戰斗部爆炸前運動狀態的影響，靜爆測量的方法不能準確體現動爆破片的參數。動態破片參數測量是針對彈藥終點效應試驗中對彈丸(戰斗部)爆炸時的動爆破片速度、破片速度衰減系數、破片空間分布等參數的測試，為毀傷效能評估試驗提供數據支撐。

在戰斗部動爆測量方面，國外起步較早。國外靶場已有采用多普勒雷達測量破片速度的動爆測量法：如丹麥BS250雷達能測初速，還能得到阻力系數曲線、速度距離曲線、速度時間曲線；美國陸軍白沙導彈靶場裝備的靶場測試雷達，工作在X波段，能跟蹤速度為30~3 000 m/s的目標，測速精度為0.3 m/s。

國內關于破片的測量，主要是依靠靜爆測量的方法進行測試，文獻[4-7]等提到的都屬于國內發展較快的靜爆測量方法。國內動爆測量技術發展相對緩慢。文獻[8]雖然提到了采用雷達回波實現對彈片和沖擊波速度進行測量，但該文側重尋找加速階段與減速階段分界點的方法研究，并未介紹戰斗部破片動爆測量的具體實現方法。本文針對此問題，提出了基于戰斗部動爆破片特性分析的雷達回波模型。

1戰斗部動爆破片特性分析

一般情況下，破片分布是以彈藥中心軸對稱分布；部分呈為定向爆炸分布。在戰斗部特性分析時，根據本文所研究的戰斗部特點和文獻[9]研究成果，設定破片為預置鋼珠球形破片，選取戰斗部引爆前導彈飛行速度v0為900 m/s，偏航角為π/36，俯仰角為π/18，導彈坐標為(500，500，500)。同時，確定各破片靜態初速在1 900~2 000 m/s之間隨機產生。

計算破片速度衰減系數時，選取破片大氣阻力系數c=0.97，海平面空氣密度ρ0= 1.225 kg/m3，580 m處的相對空氣密度H(500) =0.952 9，單枚破片的實際質量取為m=0.002 kg，破片迎風面s=0.000 049 m2，則衰減系數K=0.013 87。

仿真結果如圖1，圖2所示(圖中坐標表示破片在立體空間的分散距離)。圖1所示為破片球形分散，分散角為0°~360°的仿真；圖2所示為破片分散運動軌跡，時間較長，包含破片下落過程，分散角為100°~120°。從這兩幅圖中，基本上可以了解破片運動的過程，有利于后續工作的進行。

圖1　球形分散0°~360°Fig.1　Spherical dispersion of 0°~360°

圖2　分散運動軌跡100°~120°Fig.2　Dispersiontrajectory of 100°-120°

表1給出了破片的主要特點。

表1　測量參數特點匯總

從表1中可以得出，對破片進行測量時可由對稱軸上某一個方向或幾個方向的速度、散布來表征彈藥軸對稱上爆炸特性。在非對稱方向上進行多方向測量，即可表征彈藥爆炸空間分布特性。采用測量雷達體制進行動爆測量時，還存在目標RCS小、測量目標數量多、目標有效測量時間短、目標速度變化范圍大、目標空間散布廣、測量設備距離目標區域遠等顯著特點。

圖3、圖4給出了上述條件下，±0.6°雷達波束內破片的數量變化情況和距離-速度分布情況。

從仿真結果中可以得出，在20 ms之前，雷達波束內的破片數量急劇減少，且在8 ms左右破片數量下降90%；20 ms以后，波束內目標下降到20個以下。從圖4中可以看出，20 ms時，波束內不僅破片數量較少，而且速度、距離分布較為分散，且分為來向和去向兩個方向。如果對仿真圖進行放大，可以清晰地分辨出各個破片。

要想較完整測量爆炸分布形狀，就需要在多個方向上同時對爆炸點進行相關參數測量，然后進行數據綜合，最終得到整個爆炸范圍內的分布特性或某些重點關注方向上的爆炸分布特性。

測量時，按照所需測量方向將測試設備架設，采用統一測量標校手段進行測量單站定位定向；采用系統各測量分機的同步與啟動工作；各測量分機測量處理后的數據經時統、站址信息進行融合處理，即可得到綜合分布信息。測量示意圖如圖5所示。

圖3　±0.6°波束內破片數量隨時間的變化規律Fig.3　±0.6°The number variation of fragmentationwith time within beam

圖4　20 ms時破片距離-速度分布Fig.4　Distance-velocity distribution offragments by 20ms

圖5　爆炸點速度多角度測量示意圖Fig.5　Schematic on speed multi-anglemeasurement about explosion

2動爆破片雷達回波信號動態建模

2.1回波信號

要想完成測量，實現對破片目標的分辨，信號選擇對于雷達波束內多目標的分辨來說很重要。為完成對戰斗部動爆破片的測量，下面本文將從信號與回波建模角度進行分析。

相位編碼信號近似具有圖釘型的模糊圖[10]。理想圖釘型模糊函數不存在距離和多普勒耦合，能給出良好的鄰近目標的距離和速度分辨能力及測距、測速精度。在面對測量目標RCS小、測量目標數量多、目標有效測量時間短、目標速度變化范圍大、目標空間散布廣、測量設備距離目標區域遠等顯著特點時，本文決定采用二相編碼的信號形式進行分析，信號復數形式如下：

(1)

(2)

回波信號包含了目標距離信息、目標運動的動態信息(速度、加速度等，也是指多普勒頻率信息)、目標的RCS信息等內容。其中，距離信息和目標運動的動態信息主要包含在時間延遲τi中，目標的RCS信息主要包含在回波幅度系數Ai中。

2.2回波幅度系數的計算

記雷達到爆炸點的徑向距離為R0，對于散射截面積為σ的目標而言，其接收回波功率可表示為：

(3)

式(3)中，并未考慮電磁波能量在傳播過程中的衰減。在動爆破片測量中，對于任意一個破片(記到雷達徑向距離為Ri)，其接收回波功率為：

(4)

根據回波時間延遲τi的計算，可得回波測量位置的破片滿足：

(5)

為表示方便，對破片接收回波功率做歸一化。采用標準預置鋼珠作為參考目標，其散射截面積記為σ0(假設參考目標各向同性，運動過程中不存在翻轉和閃爍等變化)，以爆炸點作為參考點。則歸一化接收回波功率可記為：

(6)

(7)

(8)

針對預置鋼珠破片，可簡化記：

(9)

2.3動態模型的回波延遲和多普勒頻率計算

(10)

τi≈M/(c-v)=X1/(c-v)+X2/(c-v)

(11)

記A=X1/(c-v), B=X2/(c-v)；則：τi=A+B。其中，多普勒頻率主要由τi的A部分產生。則，多普勒頻率如下所示：

(12)

3仿真驗證

圖6—圖8為某一破片回波信號圖，以及對回波解調時有無多普勒對比圖。從仿真中可以看出，回波信號中比較全面的包含了破片信息，能夠詳細描述破片運動過程中回波的變化情況，從而描述出破片的運動狀態。

圖6　回波信號圖Fig.6　Echo signal diagram

圖7　解調后無多普勒回波頻譜圖Fig.7　Demodulation echo spectrumwithout doppler

圖8　解調后存在多普勒頻譜圖Fig.8　Doppler spectrum demodulation presence

4結論

本文提出了基于戰斗部動爆破片特性分析的雷達回波模型。該模型結合戰斗部動爆破片特性，從歸一化幅度系數、時間延遲計算和多普勒頻率計算等方面全面地反映了破片運動情況，模型屬于動態模型，更接近于真實情況。仿真實驗表明，該模型能夠較全面地反應回波信息和動爆破片的運動狀態，對動爆測量有很大指導意義。當然模型經過簡化，還存在以下不足：

1)雖然理論分析中比較完善地分析了回波幅度的變化情況，但是由于參數不便獲取，仿真分析時認為破片散射截面積基本相等，且變化不大，這一點有待改進。

2)文中采用的是勻減速運動模型，雖然在實際中可以選取較短時間，認為破片做勻減速運動，但是理論上還可以根據空氣動力學和外彈道學原理，詳細分析破片的受力和加速度變化情況。由于此項內容不屬于文中重點，因此做了簡化處理。

3)為更完整地表征動爆破片運動情況，還可以增大仿真時間長度和目標數量等等。

總之，模型雖然進行了一定的簡化，但是文中模型是針對于任意時刻的，簡化也是針對某一時間段的，仿真驗證也表明了所建模型與理論分析的一致性，能夠為戰斗部動爆測試提供理論支撐。

參考文獻:

[1]張玉令, 羅興柏, 徐龍. 多發彈丸同時起爆下破片空間分布實驗研究[J]. 火工品， 2012(1):21-24.

[2]趙錦, 倪晉平. 光幕靶測量破片群初速的方法[J].測試技術學報, 2007, 21(3): 214-218.

[3]楊彥通， 龍洪明.導彈戰斗部破片殺傷實驗誤差分析[J].國防技術基礎，2003(3):28-31.

[4]王高, 尹國鑫, 李仰軍, 等. 電阻網靶破片群速度測量方法[J].探測與控制學報， 2011,33(3):47-55.

[5]龐秋紅, 田會, 倪晉平. 六幕光幕靶測量破片群飛行參數算法 [J]. 西安工業大學學報, 2008, 28(5): 417-421.

[6]韋宏強, 王勁松, 馮進良, 等. 基于激光靶的戰斗部破片群速度測量方法[J]. 儀器儀表學報, 2008, 29(10): 2225-2229.

[7]蔡榮立, 倪晉平, 田會. 一種抗光干擾光幕靶的設計與改進[J]. 光學技術, 2006, 32(5): 790-792.

[8]謝金鵬. 基于雷達回波信號的彈片與沖擊波速度檢測技術研究[D].鄭州：鄭州大學, 2011.

[9]王澤, 童幼堂, 李立緯. 防空導彈戰斗部破片飛散運動的仿真分析[J].戰術導彈技術， 2010(6)：93-96.

[10]朱曉華. 雷達信號分析與處理[M]. 此京：國防工業出版社， 2011:151-187.

摘要:針對戰斗部動爆測試研究現狀和靜爆研究模型的不足，提出了基于戰斗部動爆破片特性分析的雷達回波模型。該模型結合戰斗部動爆破片特性，從歸一化幅度系數、時間延遲計算和多普勒頻率計算等方面全面地反映了破片運動情況，模型屬于動態模型，更接近于真實情況。仿真實驗表明，該模型能夠較全面地反應回波信息和動爆破片的運動狀態，對動爆測量有很大指導意義。

關鍵詞:戰斗部；動爆測試；雷達；回波模型

Radar Echo Model of Warhead Fragments Dynamic AnalysisHOU Jianqiang1, HAN Zhuangzhi1,PENG Gang2, LI Xinxin3, HUANG Wei3

(Department of Electronics and Optics Engineering , Ordnance Engineering College , Shijiazhuang 050003, China;

2. No. 77618 Unit of PLA，Lasa 850000, China;

3. The Fifty-fourth Institute of China Electronics Technology Group Corporation, Shijiazhuang 050003, China)

Abstract:As a more effective way of weapon system evolution test, dynamic test can provide more exact data than static test. The characteristics of warhead dynamic fragments radar echo was presented in this paper. The model of radar echo based on dynamic test was given with the simulation. The acquired Doppler frequency and time delay could indicat the real movements of the fragments.

Key words:warhead; dynamic explosion; radar; echo model

中圖分類號:TN95；TJ160

文獻標志碼:A

文章編號:1008-1194(2015)06-0077-05

作者簡介:侯建強(1990—)，男，河北衡水人，碩士研究生，研究方向：信號處理和速度測量。E-mail：877784427@qq.com。

*收稿日期:2015-05-22