航空火箭彈立靶密集度仿真分析

黃若超,高宏超,毛 瑞,何慶華

(中國兵器工業第203研究所,西安 710065)

航空火箭彈立靶密集度仿真分析

黃若超,高宏超,毛 瑞,何慶華

(中國兵器工業第203研究所,西安 710065)

炮口速度及轉速對航空火箭彈的立靶密集度有非常重要的影響。為了提高航空火箭彈的立靶密集度,文中以常規火箭彈的散布理論和外彈道理論為基礎,對各隨機因素進行分析,建立符合其特征的數學模型,運用蒙特卡洛法進行數學仿真,研究不同炮口速度與轉速對航空火箭彈立靶密集度的影響,結果表明提高炮口速度有利于提高立靶密集度,而炮口轉速則需要控制在合適范圍內,此結果對今后工程運用提供了理論參考。

炮口速度;炮口轉速;立靶密集度;蒙特卡洛法

0 引言

航空火箭彈是武裝直升機配備的主要武器之一,因威力大、射程遠、戰斗部多樣化、成本低廉,廣泛的用于攻擊地面上的各類目標。但是航空火箭彈的落點散布大,密集度低,影響了其作戰效能。

文中對彈道偏差作了全面分析,建立了彈體氣動模型、發動機模型、陣風模型和大氣模型等,為了體現實際彈體的差異帶來的彈道偏差,建立了偏差的隨機項擾動模型[1],結合航空火箭彈結構參數與氣動計算結果,利用C++語言編程求解運動方程,運用積分、插值、散布統計、分布擬合檢驗等方法對偏差擾動過程進行了彈道仿真模擬。根據蒙特卡洛仿真與飛行試驗結果,分析了炮口速度和炮口轉速對航空火箭彈立靶密集度的影響。

1 立靶密集度簡介

文中提到的立靶密集度是描述航空火箭彈散布的常用指標,是指在立靶密集度試驗中,隨機彈著點相對火箭彈平均彈著點的偏離度[2],也稱為圓概率誤差,用E表示,其數值越小表示立靶密集度越好。

圖1 立靶密集度試驗示意圖

立靶密集度與高度散布標準差Sy和橫向散布標準差Sz密切相關[3]:

E=0.674 5×

2 仿真與結果分析

仿真試驗基于以下假設:

1)立靶處散布僅與上文提到的隨機項擾動有關,且其擾動規律不因炮口速度及轉速的改變而改變;

2)各項擾動模型參數根據試驗數據建模仿真擬合得到。

2.1 偏差項擾動模型

利用蒙特卡洛法仿真時,隨機誤差主要考慮推力偏心角、偏心距,結構參數差異,溫度氣象差異,以及發射時俯仰、偏航的初始擾動跳動量差異。主要擾動模型建立如下[4]:

1)俯仰、偏航初始擾動[5]

俯仰、偏航的初始擾動均服從正態分布,其標準差根據某型航空火箭彈飛行試驗多發彈落點,經數學仿真擬合確定,并且其標準差不因炮口速度改變而變化。

2)隨機陣風擾動

隨機陣風風速矢量W可表示為常值風速Wc和隨機風速ΔW的矢量和:

W=Wc+ΔW

其中ΔW是沿時間和空間隨機變化的。隨機風速ΔW的標準差為σΔW:

σΔW=

模型中隨機陣風的大小服從正態分布,標準差按照上式隨機生成,隨機陣風的方向為水平方向服從均勻分布(不考慮垂直風Wy)。

3)推力偏心角、偏心距擾動

推力偏心距d服從于正態分布,其分布特征用標準差σd表示。

初始偏心角δ是在[0,2π]內服從均勻分布的隨機變量。

4)靜不平衡擾動

質心偏離幾何軸線的線偏差h0服從正態分布,分布特性σh大小由某航空火箭彈彈體參數確定。

質心偏離方位是服從[0,2π]區間的均勻分布的隨機變量。

5)動不平衡擾動

慣性中心主軸在縱向對稱面的投影與彈體縱軸的夾角λ和慣性中心主軸與彈體縱向對稱平面之間的夾角μ服從正態分布,其標準差分別為σλ和σμ。

6)氣動偏差擾動[6]

阻力偏差是服從正態分布的隨機變量,其分布特性用標準差σX表示。

2.2 偏差項分析與仿真

建立火箭彈彈道偏差復數模型[7],對方程進行簡化求解,分別得出各項擾動所引起的彈道偏差復式方程,從而進一步分析影響原因:

航空火箭彈的初始角速度是初始擾動引起散布的主要原因,且當S?1(S為陀螺穩定因子)時,偏角引起的擾動基本等同非旋轉彈情況,因此初始擾動引起的偏角復式方程為:

橫風引起的偏角復式方程為:

推力偏心距引起的偏角復式方程為:

動不平衡引起的偏角復式方程為:

氣動偏心引起的偏角復式方程為:

建立偏差項仿真模型[8],圖2是航空火箭彈在不同炮口速度下的攻角值,可以看出隨著炮口速度的增加,航空火箭彈攻角值減小,即炮口速度大的航空火箭彈穩定性更好。

圖2 不同炮口速度下攻角值

通過圖3可以看出炮口轉速較小時,推力偏心矩周期長,引起的角偏差大,散布就大,當增加炮口轉速到一定值時,推力偏心矩周期趨于相同,引起的角偏差也就趨于相同,從而當炮口轉速增大到一定值時,立靶密集度幾乎不再變化。

圖3 不同炮口轉速下y軸方向推力偏心力矩

根據圖4可看出,隨著炮口轉速的增大,動不平衡擾動力矩增大,即動不平衡引起的角偏差增大,散布變大。因此當轉速增加到一定值時,立靶密集度反而隨著轉速的增加而變差。

圖4 不同炮口轉速下動不平衡擾動力矩

2.3 立靶密集度仿真

首先驗證炮口速度對立靶密集度的影響。在總沖不變的基礎上,建立5組推力模型,獲取不同炮口速度,利用蒙特卡洛打靶法,針對不同初速分別設置200條彈道仿真,仿真結果如表1所示。

表1 炮口速度變化對立靶密集度的影響

其次驗證不同炮口轉速對立靶密集度的影響。建立賦旋角影響炮口轉速數學模型。通過調整賦旋角角度來控制出炮口轉速,初始條件不變,利用蒙特卡洛打靶法,針對不同炮口轉速分別設置200條彈道仿真,仿真結果如表2所示。

表2 炮口轉速變化對立靶密集度的影響

根據表1、圖5,可看出立靶密集度隨著炮口速度的增加而提高。

根據表2、圖6可看出,在較低的轉速范圍內,立靶密集度隨著炮口轉速的提高而提高,當炮口轉速繼續增加時,對立靶密集度影響幾乎可以忽略。文中建立的航空火箭彈模型靜穩定度很大,動不平衡引起的密集度變化幾乎不能體現。

圖5 不同炮口速度靶面彈著點坐標

圖6 不同炮口轉速靶面彈著點坐標

3 立靶密集度試驗與分析

為了驗證偏差項仿真與蒙特卡洛仿真結果的正確性,以某型航空火箭彈為分析對象,進行飛行試驗驗證。

設計不同推力方案,通過立靶密集度試驗,驗證出炮口速度對密集度的影響,結果如表3所示。

表3 不同出炮口速度立靶測試結果

設計不同賦旋角方案,通過立靶密集度試驗,驗證了出炮口轉速對密集度的影響,測試結果如表4所示。

表4 不同出炮口轉速下立靶測試結果

根據仿真與試驗結果可知,立靶密集度隨炮口速度的增加而提高。炮口轉速需控制在合適范圍內,因為隨著轉速的增大立靶密集度變化不大,反而犧牲射程。航空火箭彈靜穩定度較高,所以轉速變化帶來的動不穩定影響不明顯。理論分析與試驗結果一致。

4 結論

文中分析了炮口速度與炮口轉速對航空火箭彈立靶密集度的影響。結果表明,提高炮口速度可以有效提高立靶密集度,而炮口轉速則需要合理設計,尤其對于靜穩定度不高的航空火箭彈,能使其同時兼顧立靶密集度與動不平衡擾動帶來的動不穩定影響。

[1] 高帥. 遠程火箭彈的散布研究 [D]. 南京: 南京理工大學, 2004.

[2] 鄭永杰. 減小航空火箭彈散布方法的研究及其優化 [D]. 南京: 南京理工大學, 2003.

[3] 王寶元, 衡剛, 周發明, 等. 基于坐標極差估算高炮立靶密集度的方法 [J]. 測試技術學報, 2014, 28(4): 356-357.

[4] 錢杏芳. 導彈飛行力學 [M]. 北京: 北京理工大學出版社, 2008: 137-148.

[5] 王志軍, 趙文宣. 火箭彈起始擾動數值仿真研究 [J]. 彈道學報, 1996, 8(3): 70-73.

[6] 湯江河. 影響直升機載航空火箭射擊精度的因素分析 [D]. 南京: 南京理工大學, 2007.

[7] 楊紹卿. 火箭外彈道偏差與修正理論 [M]. 北京: 國防工業出版社, 2011: 98-144.

[8] 毛保全, 鄧威, 馮帥, 等. 遠基于多學科協同仿真的立靶密集度計算方法 [J]. 計算機仿真, 2016, 33(1): 20-22.

Simulation and Analysis of Vertical Target Dispersion of Airborne Rocket

HUANG Ruochao,GAO Hongchao,MAO Rui,HE Qinghua

(No.203 Research Institute of China Ordnance Industries, Xi’an 710065, China)

The muzzle velocity and the muzzle rotating speed had very important influence on the vertica target dispersion of airborne rocket. In order to improve the vertical target dispersion of airborne rocket, base on the dispersion theory and exterior ballistic theory of conventional rocket, the random factors were analyzed to establish the mathematical model in line with its characteristics. By using the Monte Carlo method for mathematical simulation, the influence of different muzzle velocity and muzzle rotating speed on the vertica target dispersion of airborne rocket was studied. The results showed that increasing muzzle velocity was helpful to improve the vertical target density, while muzzle speed needed to be controlled within the appropriate range. This conclusion provided theoretical references for future engineering applications.

muzzle velocity；muzzle rotating speed；vertical target dispersion；Monte Carlo method

2016-11-03

黃若超(1991-),女,陜西藍田人,助理工程師,碩士研究生,研究方向:航空火箭彈總體技術。

V435；V211.48

A