趙養(yǎng)正,呂鴻鷹,向玉偉

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

?

具有瓦片翼的旋轉(zhuǎn)子彈氣動(dòng)特性數(shù)值研究

趙養(yǎng)正,呂鴻鷹,向玉偉

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

瓦片翼;高阻力;旋轉(zhuǎn)子彈;數(shù)值模擬

0 引言

末修彈藥是常規(guī)彈藥制導(dǎo)化的重要部分。當(dāng)末修子彈拋出母彈彈體后,瓦片尾翼打開,不但對(duì)末修子彈進(jìn)行減速減旋,而且起到穩(wěn)定子彈姿態(tài)的作用。彈翼采用瓦片翼形狀是典型的阻力翼,可以增大子彈飛行阻力,同時(shí)迅速降低子彈的轉(zhuǎn)速,從而保證末修子彈在一定的速度和轉(zhuǎn)速下搜索目標(biāo)。目前對(duì)旋轉(zhuǎn)彈箭的氣動(dòng)特性的研究較多,文獻(xiàn)[1]通過試驗(yàn)的方法研究了平頭圓柱體帶四片升力翼的子彈的滾轉(zhuǎn)阻尼特性,文獻(xiàn)[2-5]通過準(zhǔn)靜態(tài)數(shù)值的方法研究了旋轉(zhuǎn)狀態(tài)下的動(dòng)導(dǎo)數(shù),文獻(xiàn)[6]研究了一種異形卷弧翼彈旋轉(zhuǎn)狀態(tài)下的氣動(dòng)特性,但研究這種具有瓦片翼子彈的氣動(dòng)特性的文獻(xiàn)較少,因此通過數(shù)值方法研究具有瓦片翼子彈的氣動(dòng)性能意義重大。

1 模型

模型結(jié)構(gòu)為臺(tái)階型彈體+3片瓦面型尾翼(見圖1)。此氣動(dòng)外形具有特殊的氣動(dòng)特性,即升力小、阻力大,可使彈體速度、轉(zhuǎn)速快速降低,增大彈道落角,短時(shí)間內(nèi)達(dá)到最終子彈穩(wěn)定工作狀態(tài)。

圖1 模型示意圖

2 數(shù)值方法

2.1 計(jì)算方法

文中使用CFD軟件cfx進(jìn)行數(shù)值計(jì)算。具體計(jì)算方法為:結(jié)合SST湍流模型,并采用有限體積中心差分格式對(duì)三維雷諾平均Navier-Stokes方程進(jìn)行求解。由于模型是旋轉(zhuǎn)的,因此采用了旋轉(zhuǎn)坐標(biāo)系[1],坐標(biāo)系固定在子彈體上,隨子彈一起以固定的角速度旋轉(zhuǎn)。入口邊界采用速度和溫度;出口邊界采用靜壓和溫度;遠(yuǎn)場邊界采用無粘滑移邊界墻;壁面邊界采用無滑移壁面條件。

2.2 計(jì)算網(wǎng)格

外場采用四面體網(wǎng)格,附面層采用三棱柱形網(wǎng)格。模型的表面網(wǎng)格和附面層網(wǎng)格見圖2和圖3。

圖2 表面網(wǎng)格

圖3 附面層網(wǎng)格

2.3 計(jì)算條件

1)靜態(tài)計(jì)算條件

馬赫數(shù):0.2～1.5;攻角:0°～30°;

2)動(dòng)態(tài)計(jì)算條件

馬赫數(shù):0.2～0.8;攻角:0°～12°;

轉(zhuǎn)速:5～15 r/s。

3 結(jié)果與討論

文中對(duì)靜態(tài)和旋轉(zhuǎn)狀態(tài)的氣動(dòng)特性進(jìn)行了研究,并對(duì)靜態(tài)的計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果進(jìn)行了比較分析。

3.1 靜態(tài)結(jié)果

由圖4可看出:正攻角時(shí),升力系數(shù)為負(fù)值,升力系數(shù)隨攻角變化為負(fù)線性的。此變化規(guī)律與常規(guī)彈剛好相反,因?yàn)槌Ｒ?guī)彈的法向力大,軸向力小,而此彈的法向力小(見圖5),軸向力大(見圖6),雖然法向力正攻角時(shí)為正值,但坐標(biāo)系轉(zhuǎn)化時(shí)要減去大的軸向力分量,導(dǎo)致升力在正攻角時(shí)為負(fù)值。

體軸系轉(zhuǎn)化為風(fēng)軸系的公式為:

CD=CAcosα+CNsinα

CL=CNcosα-CAsinα

其中:α為攻角;CA、CN為體軸系的軸向力系數(shù)和法向力系數(shù);CD、CL為風(fēng)軸系的阻力系數(shù)和升力系數(shù)。

從圖7的阻力系數(shù)曲線可以看出:在所研究的攻角范圍內(nèi),阻力系數(shù)比常規(guī)彈箭大得多,大約是標(biāo)準(zhǔn)彈箭模型(CBM)零升阻力的10倍,而且阻力隨攻角變化較小;同時(shí)阻力隨馬赫數(shù)的變化也小。因此這類彈的阻力受攻角的影響很小,阻力平穩(wěn)。

圖4 升力系數(shù)隨攻角的變化曲線

圖5 法向力系數(shù)隨攻角的變化曲線

圖6 軸向力系數(shù)隨攻角的變化曲線

圖7 阻力系數(shù)隨攻角的變化曲線

圖8 法向力系數(shù)導(dǎo)數(shù)隨馬赫數(shù)的變化曲線

圖9 俯仰力矩系數(shù)隨攻角的變化曲線

圖10 俯仰力矩系數(shù)隨法向力系數(shù)的變化曲線

3.2 旋轉(zhuǎn)狀態(tài)結(jié)果

在旋轉(zhuǎn)狀態(tài)結(jié)果中,將靜態(tài)結(jié)果作為轉(zhuǎn)速n=0 r/s處理,從圖11(a)可以看出隨著攻角的增大,在不同轉(zhuǎn)速下法向力系數(shù)出現(xiàn)了分化,在7 r/s以下法向力是隨著轉(zhuǎn)速的提高逐漸減小;然而,當(dāng)轉(zhuǎn)速處于10 r/s以上時(shí),法向力變?yōu)樨?fù)值,轉(zhuǎn)速越高,法向力絕對(duì)值越小,從圖11(b)可以看出法向力系數(shù)隨轉(zhuǎn)速的增大,經(jīng)歷了先減小后增大的過程,在8 r/s附近不同攻角下的法向力系數(shù)接近零,是變化的轉(zhuǎn)折點(diǎn)。比較圖12和圖13可以看出,轉(zhuǎn)速為5 r/s時(shí)(見圖12),下半表面的壓力大于上半表面的壓力,因此法向力為正;轉(zhuǎn)速為10 r/s時(shí)(見圖13),下半表面的壓力小于上半表面的壓力,因此法向力為負(fù),這是由于轉(zhuǎn)速增加,使近表面的流場壓力分布發(fā)生改變造成的。

圖14為俯仰力矩系數(shù)(對(duì)質(zhì)心)隨攻角及轉(zhuǎn)速的變化曲線,從圖14(a)可以看出不同轉(zhuǎn)速下的俯仰力矩系數(shù)較小,尤其是轉(zhuǎn)速為15 r/s時(shí)俯仰力矩系數(shù)最接近零,同一攻角下轉(zhuǎn)速為10 r/s時(shí)的俯仰力矩系數(shù)最大,而轉(zhuǎn)速為15 r/s時(shí)的俯仰力矩系數(shù)最小。從圖14(b)可以看出同一攻角下,攻角小于6°和轉(zhuǎn)速小于等于7 r/s時(shí)俯仰力矩的變化較小;在攻角大于6°和轉(zhuǎn)速大于7 r/s時(shí)俯仰力矩的變化較大。

圖11 法向力系數(shù)隨攻角及轉(zhuǎn)速的變化曲線(Ma=0.2)

圖12 模型法向力分布(r=5 r/s,Ma=0.2,α=6°)

圖13 模型法向力分布(n=10 r/s,Ma=0.2,α=6°)

圖14 俯仰力矩系數(shù)隨攻角及轉(zhuǎn)速的變化曲線(Ma=0.2)

圖15 俯仰力矩系數(shù)隨法向力系數(shù)的變化曲線(Ma=0.2)

3.3 靜穩(wěn)定性

3.4 動(dòng)穩(wěn)定性

下面通過火箭外彈道學(xué)方法用導(dǎo)彈穩(wěn)定性的判定準(zhǔn)則對(duì)旋轉(zhuǎn)子彈的靜、動(dòng)穩(wěn)定性加以判定。

1)靜穩(wěn)定條件

陀螺穩(wěn)定因子[7]:

2)動(dòng)穩(wěn)定條件

動(dòng)穩(wěn)定因子[8]:

慢圓運(yùn)動(dòng)穩(wěn)定條件:

通過對(duì)該彈Ma=0.2時(shí)動(dòng)穩(wěn)定性的計(jì)算表明,Sd<0,為慢圓運(yùn)動(dòng)。高速旋轉(zhuǎn)狀態(tài)下n=15 r/s時(shí),該彈滿足慢圓運(yùn)動(dòng)穩(wěn)定條件;當(dāng)n≤10 r/s時(shí),該彈為陀螺不穩(wěn)定,動(dòng)不穩(wěn)定,因此該彈在較高速旋轉(zhuǎn)狀態(tài)是動(dòng)穩(wěn)定的。

4 結(jié)論

文中通過數(shù)值計(jì)算和實(shí)驗(yàn)方法研究了具有瓦片翼子彈的氣動(dòng)特性,并通過對(duì)比分析得出如下結(jié)論:

1)采用文中的數(shù)值方法計(jì)算的氣動(dòng)參數(shù)實(shí)用、可靠,與試驗(yàn)結(jié)果一致性非常好。

2)與常規(guī)彈箭不同,此彈的法向力小,軸向力大,雖然法向力在正攻角時(shí)為正值,但正攻角時(shí),升力系數(shù)為負(fù)值,升力系數(shù)隨攻角變化為負(fù)線性關(guān)系。

3)該彈賦予不同的旋轉(zhuǎn)速度后出現(xiàn)了法向力正負(fù)分化的現(xiàn)象。

4)在小攻角和低轉(zhuǎn)速下旋轉(zhuǎn)運(yùn)動(dòng)對(duì)俯仰力矩的影響較小;在大攻角和高轉(zhuǎn)速下旋轉(zhuǎn)運(yùn)動(dòng)對(duì)俯仰力矩的影響較大。

5)在子彈低速旋轉(zhuǎn)時(shí)仍然可以應(yīng)用旋轉(zhuǎn)彈的靜態(tài)參數(shù)來判定子彈的穩(wěn)定性,但高速時(shí)就不能應(yīng)用。

6)Ma=0.2時(shí),該彈要通過高速旋轉(zhuǎn)的陀螺效應(yīng)才能實(shí)現(xiàn)動(dòng)穩(wěn)定。

[1] 趙養(yǎng)正, 蔣勝矩, 黨明利, 等. 旋轉(zhuǎn)彈體及減旋片滾轉(zhuǎn)阻尼數(shù)值模擬 [J]. 兵工學(xué)報(bào), 2015, 36(7): 1176-1180.

[2] 丁則勝, 劉亞飛, 徐琴, 等. 彈體及減旋片滾轉(zhuǎn)阻尼實(shí)驗(yàn)研究 [J]. 彈道學(xué)報(bào), 2001, 13(1): 62-65.

[3] 鄧帆, 陳少松, 陶鋼. 帶柵格翼導(dǎo)彈超聲速階段滾轉(zhuǎn)阻尼導(dǎo)數(shù)的數(shù)值研究 [J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2012, 30(2): 151-156.

[4] MACHAEL A P, LAWRENCE L G. steady-state computation of constant rotational rate dynamic stability derivatives: AIAA 2000-4321 [R]. 2000.

[5] DESPIRITO J, SILTON S I, WEINACHT P. Navier-Stokes predictions of dynamic stability derivatives-evaluation of steady-state methods: Maryland, USA: ARL-TR-4605 [R]. 2008.

[6] 金華, 戴金海, 吳蓓蓓. 一種異形卷弧翼彈氣動(dòng)特性的數(shù)值模擬 [J]. 彈箭與制導(dǎo)學(xué)報(bào), 2006, 26(3): 168-171.

[7] 吳甲生. 雷娟棉. 制導(dǎo)兵器氣動(dòng)布局與氣動(dòng)特性 [M]. 北京: 國防工業(yè)出版社, 2008: 299-306.

[8] MURPHY Charles H. Free flight motion of symmetric missiles: BRL Report No.1216 [R]. 1963.

The Aerodynamic Performance Simulation of a Spinning Projectile with Tile Wings

ZHAO Yangzheng,LYU Hongying,XIANG Yuwei

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

tile wing; larger drag; spinning projectile; numerical simulation

2015-11-03

趙養(yǎng)正(1965-),男,陜西藍(lán)田人,高級(jí)工程師,碩士,研究方向:彈箭空氣動(dòng)力學(xué)與實(shí)驗(yàn)流體力學(xué)。

V211.3

A