王朋飛， 周前進

(1.西安航空職業(yè)技術學院， 西安 710089； 2.國營990廠， 合肥 230601)

鴨舵式火箭彈的氣動特性及控制能力分析

王朋飛1， 周前進2

(1.西安航空職業(yè)技術學院， 西安 710089； 2.國營990廠， 合肥 230601)

為了提高火箭彈的密集度，對鴨舵改變火箭彈氣動布局問題進行了分析?；邙喍媸交鸺龔椀氖芰Ψ治?，進行了不同攻角和舵偏角下火箭彈氣動特性數(shù)值仿真，并利用仿真數(shù)據(jù)進行了縱橫向控制能力計算。結(jié)果表明：鴨舵的存在使火箭彈的阻力系數(shù)和升力系數(shù)都有所增加，通過控制鴨舵的起控高度和舵偏角大小可得到火箭彈在縱橫向不同的控制量，從而為彈道修正彈的彈道設計提供依據(jù)。

鴨舵；氣動特性；舵偏角；修正能力

為了實現(xiàn)火箭彈完成精確打擊的能力，通過簡易控制技術對火箭彈加以改造，使之轉(zhuǎn)化成彈道修正火箭彈，達到提高射擊精度、降低成本的目的，已是火箭彈發(fā)展主要方向。實現(xiàn)二維彈道修正執(zhí)行機構主要有微型脈沖推力火箭、燃氣射流控制發(fā)動機和空氣動力鴨舵[1-2]。脈沖修正執(zhí)行機構的成本低、結(jié)構簡單，但只能提供有限次、推力恒定的離散式修正，且需對全彈進行全新設計；鴨舵式執(zhí)行機構能提供連續(xù)的控制力、控制能力大、控制精度高，且容易調(diào)整全彈質(zhì)心，因此可對散布較大的火箭彈進行改造，提高射擊精度[3-4]。

本文在分析了鴨舵式火箭彈受力情況的基礎上，利用多種網(wǎng)格技術建模的方式對無舵和有舵火箭彈進行了氣動仿真，分析了不同舵偏角下鴨舵上下面的流場分布情況及攻角、舵偏角對火箭彈氣動特性的影響；同時根據(jù)仿真的氣動數(shù)據(jù)對鴨舵式火箭彈進行彈道仿真，計算鴨舵起控高度和舵偏角大小對火箭彈控制能力的影響。

1 受力分析

與無控火箭彈相比，鴨舵式火箭彈的氣動力和氣動力矩增加了由舵面產(chǎn)生的控制力和控制力矩[5]，其受力如圖1所示。圖中δ為攻角，α為舵偏角；Rx、Ry和Mz為火箭彈攻角引起的氣動阻力、升力和俯仰力矩；DRx、DRy、DMs為由舵偏角引起的附加氣動阻力、升力、俯仰力矩。

圖1 舵面偏轉(zhuǎn)時的受力示意圖

2 氣動仿真分析

2.1 網(wǎng)格模型及仿真條件

為消除流場域邊界對仿真結(jié)果的影響，計算區(qū)域為長15倍彈長，直徑20倍彈徑的圓柱體。由于彈體結(jié)構比較復雜，所以采用分塊對接的方式將計算區(qū)域分成若干區(qū)域，采用多種網(wǎng)格技術，以保證網(wǎng)格的質(zhì)量[6-7]。彈體的網(wǎng)格模型如圖2、圖3所示。

圖2 彈頭部網(wǎng)格模型

圖3 彈尾部網(wǎng)格模型

2.2 外流場力學模型

流體流動需遵守的定律包括：質(zhì)量守恒定律、動量守恒定律、能量守恒定律。流體動力學控制方程就是這些守恒定律的數(shù)學描述[8-9]。

1) 質(zhì)量守恒方程

質(zhì)量方程又稱為連續(xù)方程，此方程是質(zhì)量守恒方程的一般形式，適用于可壓流動和不可壓流動。

2) 動量守恒方程

動量守恒方程又稱運動方程，式中:p是靜壓；τij為應力張量；ρgi和Fi分別為i方向的重力體積力和外部體積力。

3) 能量守恒方程

式中：keff為有效熱傳導系數(shù)；Jj′為組分j′的擴散流量。

2.3 氣動特性分析

仿真模型的彈體、舵面和尾翼表面采用無滑移固壁(wall)邊界條件，計算區(qū)域的入口、出口和外邊界全采用壓力遠場邊界條件，氣體模型為理想氣體模型，動力黏度計算采用薩蘭德定理。

2.3.1 不同攻角下的氣動特性分析

圖4是舵偏角α=0°時，阻力系數(shù)隨馬赫數(shù)、攻角的變化曲線。由圖4可知，在小攻角范圍內(nèi)，隨著攻角的增大，攻角引起的誘導阻力也隨之增大，使全彈阻力系數(shù)隨著增加，且全彈阻力系數(shù)隨馬赫數(shù)的變化規(guī)律基本一致[10]。

圖4 阻力系數(shù)變化曲線

圖5、圖6分別是攻角δ=4°，舵偏角α=0°時，無控火箭彈和鴨舵式火箭彈阻力系數(shù)和升力系數(shù)隨馬赫數(shù)的變化曲線。由其可知，在同一攻角下，鴨舵式火箭彈的阻力系數(shù)、升力系數(shù)比同馬赫數(shù)下的無控火箭彈要大。這說明有攻角時，鴨舵的存在增加了全彈的升力，同時全彈阻力也相應增加。

2.3.2 舵偏角對流場的影響

圖7～圖9為鴨舵式火箭靶彈在馬赫數(shù)Ma=1.5、攻角δ=0°時舵面截面處(平行于彈體縱向平面)壓力分布云圖。

由圖7～圖9可知，鴨舵的前緣出現(xiàn)了激波，迎風面氣流壓縮劇烈，形成高壓區(qū)；舵偏角的存在引起舵面上下區(qū)域流場不對稱，下區(qū)域流場的壓力明顯高于上區(qū)域，這個壓力差使鴨舵產(chǎn)生附加的氣動力和氣動力矩；隨著舵偏角的增大，舵面上下區(qū)域流場的壓力差越大。因此，可通過控制鴨舵起控高度和舵偏角，調(diào)整彈體飛行姿態(tài)，實現(xiàn)火箭彈在縱向和橫向上的彈道修正。

圖5 δ=4°、α=0°時阻力系數(shù)對比曲線

圖6 δ=4°、α=0°時升力系數(shù)對比曲線

圖7 舵偏角α=0°

圖8 舵偏角α=5°

圖9 舵偏角α=10°

3 控制能力分析

第2小節(jié)分析了不同舵偏角大小對應的氣動特性影響火箭彈縱橫向控制能力的情況。鴨舵在彈道上起控高度和舵偏角大小將影響火箭彈的控制能力[11-12]。為此基于前文的氣動力參數(shù)，進行了不同高度和舵偏角下火箭彈在縱橫向的控制能力計算，計算結(jié)果如表1所示。在縱向控制，10 km高度起控時，飛行高度隨時間和舵偏角變化曲線如圖10所示。由表1的結(jié)果可見：鴨舵在彈道上起控的高度越高，在縱向和橫向的控制量越大，反之起控的高度越?。欢嫫窃酱?，其控制量也越大。

表1 不同舵偏角下的控制能力 m

圖10 飛行高度-時間變化曲線

4 結(jié)論

分析了不同攻角和舵偏角對火箭彈氣動特性的影響規(guī)律，并針對鴨舵起控高度和舵偏角大小引起的控制能力問題進行了彈道仿真分析。隨著攻角的增大，全彈的阻力系數(shù)隨著增大；舵偏角的存在，引起了舵面上下表面的壓力差，可為火箭彈彈道修正提供控制力和力矩；在彈道上鴨舵起控高度越高，舵偏角越大，其縱橫向的控制能力越大。在10 km高度起控時，6°舵偏角下的縱向控制量將近是4°舵偏角下控制量的1.8倍。

[1] 胡金波，楊新民，何穎,等.二維彈道修正迫彈氣動特性仿真[J].四川兵工學報，2015，36(12)：43-47.

[2] 王朋飛，曹紅松.引信頭錐擺動角對火箭彈氣動特性及控制能力影響[J].彈箭與制導學報，2013，33(1)：152.

[3] 孟超.鴨舵式修正控制機構設計[D].北京：北京理工大學，2009：7-12.

[4] ERIC J A，MELISSA A.McDaniel.High angle of attack aerodynamic predictions using missile datcom[C]//23rd AIAA Applied Aerodynamics Conference.[S.l.]:[s.n.],2005：1-18.

[5] 馬寶成.基于旋轉(zhuǎn)舵的彈道控制技術研究[D].太原：中北大學，2012：40-42.

[6] SLOTNIK J P，KANDULA M，BUNING P G.Navier stokes simution of the space shuttle of the space launch vehicle flight trandonic flowfield using a large scale chimera crid system[J].AIAA 94-1860，1994.

[7] 伍星，盧永剛，宋瓊,等．基于Fluent的彈體氣動特性計算與分析[J].兵器裝備工程學報，2016(2)：22-25.

[8] ANDERSON J D.Computational Fluid Dynamics[M].北京：清華大學出版社，2002.

[9] KOOMULLIL R P，SONI B K，HUANG C T．Navier-stokes simulation on hybrid grids[J].AIAA 96-0768，1996.

[10] 張嘉易，王廣，郝永平．二維彈道修正鴨舵修正機構氣動特性研究[J].彈箭與制導學報，2013，33(2)：88-91.

[11] 汪亞利，邵偉平，郝永平,等．舵片修正彈丸氣動特性仿真研究[J].兵器裝備工程學報，2016，37(6)：65-68.

[12] 于濤．二維修正彈的彈道特性研究[D].南京：南京理工大學，2010：8-15.

(責任編輯周江川)

AnalysisofCanardRocketforAerodynamicCharacteristicsandControlAbility

WANG Pengfei1， ZHOU Qianjin2

(1 Xi’an Aeronautical Polytechnic Institute, Xi’an 710089, China； 2.State Factory No.990, Hefei 230601, China)

In order to improve the intensity of the rocket, the problem of changing the aerodynamic layout of the canard is analyzed. Based on the analysis of the canard rocket, the numerical simulation of the aerodynamic characteristics of the rocket under different angle of attack and elevator angle is carried out, and the vertical and horizontal control capability is calculated by using the simulation data. The results show that the existence of the canard has increased the resistance coefficient and the lift coefficient of the rocket. By controlling the height of the canard and the elevator angle, the rocket can get different control in the vertical and horizontal directions, the ballistic design of the bullet is provided.

canard; aerodynamic characteristics; elevator angle; correction capability

2017-07-04；

2017-07-31

王朋飛(1986—)，男，碩士，主要從事氣動仿真和飛行技術研究。

裝備理論與裝備技術

10.11809/scbgxb2017.11.009

本文引用格式：王朋飛，周前進.鴨舵式火箭彈的氣動特性及控制能力分析[J].兵器裝備工程學報，2017(11):41-44.

formatWANG Peng-fei，ZHOU Qian-jin.Analysis of Canard Rocket for Aerodynamic Characteristics and Control Ability[J].Journal of Ordnance Equipment Engineering,2017(11):41-44.

TJ714

A

2096-2304(2017)11-0041-04