一種帶副翼無傘末敏彈氣動特性仿真分析*

李 恒,郭 銳,劉榮忠

(南京理工大學機械工程學院,南京 210094)

一種帶副翼無傘末敏彈氣動特性仿真分析*

李 恒,郭 銳,劉榮忠

(南京理工大學機械工程學院,南京 210094)

為了探究不同尾翼末敏彈的氣動特性,設計了帶有副翼結構的無傘末敏彈模型,借助計算流體力學建立了無傘末敏彈的氣動力仿真模型。數值計算結果表明不同尾翼結構對末敏彈的穩態落角和阻力系數有較大影響。當翼展從150 mm增至250 mm時,末敏彈的穩態落角從43°減至16°,阻力系數從0.78增至1.61。當副翼斜置角從0°增至60°時,穩態落角、導轉力矩系數和阻力系數也在發生變化。這將為無傘末敏彈的設計提供有益參考。

無傘末敏彈;計算流體力學;氣動特性;副翼

0 引言

末敏彈穩態掃描技術主要包括有傘掃描和無傘掃描。大多數國家研制的末敏彈采用有傘掃描,如美國的“SADARM”末敏彈和德國的“SMART”末敏彈。末敏彈的無傘化是其一個重要發展方向。因此,國內外開展了一系列關于無傘末敏彈的理論研究和實驗研究。其中,顧建平和舒敬榮等[1-2]研究了單翼末敏彈的運動特性,胡志鵬、呂勝濤和郭銳等[3-7]研究了雙翼末敏彈的氣動特性,周志起等[8]研究了三翼末敏彈的氣動特性。

由于上述研究的無傘末敏彈的尾翼有效面積均較小,所以導致末敏彈的穩態落速較大,這就為末敏彈的掃描運動帶來極大的挑戰。基于此,文中提出一種新型可伸縮展開翼結構,以增大末敏彈的阻力,并

借助Fluent仿真軟件分析了該新型翼結構的氣動特性變化規律,文中主要研究內容包括兩部分:1)翼B展長變化時,末敏彈穩態落角和阻力系數的變化規律。2)不同副翼展長在不同斜置角下,末敏彈穩態落角、導轉力矩系數和阻力系數的變化規律。

1 氣動布局模型

圖1為兩種伸縮展開翼結構,末敏彈彈體為圓柱體,尺寸為Φ110 mm×135 mm;翼A弦長和展長分別為100 mm和300 mm;其中,圖1(a)所示為改變翼展結構,翼B弦長為100 mm,定義翼B展長為L;圖1(b)所示為增加副翼后結構,副翼弦長為100 mm,定義副翼展長為H,副翼斜置角為γ,翼B與副翼展向中心線共線;定義末敏彈的穩態落角為α,即彈軸與鉛垂方向的夾角。

2 數值仿真模型

計算域為圓柱形區域,流場軸向取22.2倍彈長,徑向取14.5倍彈徑,采用非結構網格,為了保證計算的精度,采用兩級分割,保證末敏彈附近網格致密。定義相應邊界條件。定義末敏彈參考面積S=0.069 5 m2,參考長度l=0.005 87 m(彈體質心至阻心距離)。

3 仿真結果及分析

獲得穩態落角α的方法:定義質心到彈體頂部的距離為0.6倍彈體長度,且質心位于彈軸上,在翼A中心線和彈軸所確定的平面內,對質心取矩,若矩為零,則所得攻角即為穩態落角。本例中通過不斷調試落角α的大小,采用逐漸逼近的方法,最終獲得穩態落角α。

3.1 翼B展長對氣動特性的影響

壓力分布:圖2為翼B展長取180 mm時壓力云圖。可以看出,高壓區主要分布在彈體頭部表面、彈體側壁迎風區和尾翼迎風區,低壓區主要分布在彈體尾部表面、彈體側壁背風區和尾翼背風區。大面積的高壓區有利于末敏彈減速。

圖2 翼B展長為180 mm時的壓力云圖

為了探究翼B展長L對末敏彈穩態落角α和阻力特性的影響,保持翼A弦長為100 mm,展長為300 mm,保持翼B弦長為100 mm,令翼B展長L分別為150、160、170、180、190、200、210、220、230、240和250 mm,進行仿真計算,通過不斷調試末敏彈穩態落角α的大小,最終獲得了不同翼B展長L下的穩態落角α的變化規律,如圖3(a)所示。另外,獲得了不同翼B展長的末敏彈在穩態落角工況下阻力系數的變化規律,如圖3(b)所示。

圖3 末敏彈翼B取不同展長時的氣動特性

由圖3(a)可以看出,隨著翼B展長的增大,末敏彈的穩態落角不斷減小,從43°減至16°,這是由于末敏彈質心右側的有效迎風面積隨著翼B展長的增大而增大,唯有將穩態落角減小才可以使得質心左右兩側的有效迎風面積相等,以形成穩態運動。由圖3(b)可以看出,隨著翼B展長的增大,末敏彈的阻力系數不斷增大,從0.78增至1.61,這是由于翼B展長的增大和穩態落角的減小均有利于末敏彈有效迎風面積的增大,從而使末敏彈阻力系數增加。

3.2 副翼展長和斜置角對氣動特性的影響

壓力分布:圖4為副翼斜置角取30°時壓力云圖。尾翼的形狀和安裝位置不同導致壓力分布也略有不同,當副翼有一定斜置角時,副翼下側迎風面為高壓區,上側迎風面為局部低壓區,這就產生了導轉力矩。副翼上側產生低壓區的原因類似于飛機機翼翼尖處產生的翼尖渦,在翼B的翼尖處產生翼尖渦,這導致副翼上側為局部低壓區。

圖4 副翼斜置角為30°時的壓力云圖

為探究副翼展長H和斜置角γ對末敏彈導轉特性、穩態落角α和阻力特性的影響,保持翼A弦長為100 mm,展長為300 mm,翼B弦長為100 mm,展長為180 mm。對弦長為100 mm,展長H分別為30、40和50 mm的副翼進行仿真,令副翼斜置角γ分別為0°、10°、20°、30°、35°、40°、45°、50°、和60°。通過不斷調試末敏彈落角α的大小,最終獲得不同展長副翼在不同斜置角γ下的穩態落角α的變化規律,如圖5(a)所示。另外,獲得了不同斜置角副翼末敏彈在穩態落角工況下導轉力矩系數和阻力系數的變化規律,如圖5(b)和圖5(c)所示。

由圖5(a)可以看出,當副翼斜置角一定時,副翼展長越大,則其穩態落角就越小,這是由于隨著副翼展長的增大,質心右側的有效迎風面積增大,唯有將穩態落角減小才可以使得質心左右兩側的有效迎風面積相等;當副翼展長一定時,隨著副翼斜置角的增大,穩態落角逐漸增大,這是由于隨著副翼斜置角的增大,質心右側的有效迎風面積不斷減小,唯有將穩態落角增大才可以使得質心兩側的有效迎風面積相等。由圖5(b)可以看出,當副翼斜置角一定時,副翼展長越大,則其導轉力矩系數就越大,這是由于產生導轉力矩的有效迎風面積在增大;當副翼展長一定時,隨著斜置角的增大,導轉力矩系數先增大后減小,并在45°附近取得最大值。由圖5(c)可以看出,當副翼斜置角一定時,副翼展長越大,則其阻力系數就越大,這是由于末敏彈的有效迎風面積在增大;當副翼展長一定時,隨著副翼斜置角的增大,其阻力系數不斷減小,這是由于末敏彈的有效迎風面積在減小。

圖5 末敏彈副翼取不同斜置角時的氣動特性

4 結論

1)隨著翼B展長的增大,末敏彈穩態落角不斷減小,若將末敏彈的穩態落角控制在30°附近,則可以將翼B展長取在190～210 mm。

2)可以通過增大副翼展長來增大導轉力矩系數,也可以通過將副翼斜置角設在45°附近取得最大的導轉力矩系數。

3)若副翼斜置角取45°,副翼展長取在30～50 mm,則可以實現末敏彈的穩態落角α在28°～32°范圍內變化。

[1] 顧建平, 韓子鵬. 單翼與有傘末敏子彈的掃描比較及分析 [J]. 彈箭與制導學報, 2011, 31(1): 123-124.

[2] 舒敬榮, 張邦楚, 韓子鵬, 等. 單翼末敏彈掃描運動研究 [J]. 兵工學報, 2004, 25(4): 415-420.

[3] 胡志鵬, 劉榮忠, 郭銳. 基于FLUENT的雙翼末敏彈氣動特性研究 [J]. 飛行力學, 2013, 31(1): 53-60.

[4] 胡志鵬, 劉榮忠, 郭銳. 兩種典型尾翼形狀對無傘末敏彈氣動特性的影響 [J]. 南京理工大學學報, 2012, 36(5): 739-744.

[5] 呂勝濤, 劉榮忠, 郭銳, 等. S-C型雙翼末敏彈尾翼結構方案設計 [J]. 彈道學報, 2014, 26(2): 6-12.

[6] 呂勝濤, 劉榮忠, 郭銳, 等. S-S雙翼末敏彈氣動外形優化設計 [J]. 兵工學報, 2013, 34(9): 1150-1155.

[7] 郭銳, 劉榮忠, 王本河, 等. 一種非對稱雙翼結構彈丸減速導旋特性試驗研究 [J]. 彈箭與制導學報, 2009, 29(5): 249-254.

[8] 周志起, 趙潤祥, 韓子鵬, 等. 末敏子彈氣動外形設計與氣動特性分析 [J]. 空氣動力學學報, 2013, 31(1): 15-21.

Simulation Analysis on Aerodynamic Characteristics of Non-fin TSP with Aileron

LI Heng,GUO Rui,LIU Rongzhong

(School of Mechanical Engineering. NUST, Nanjing 210094, China)

To study aerodynamic characteristics of terminal sensitive projectile (TSP) with different wings, non-fin TSP with aileron was designed. The simulation models of non-fin TSP were established by using computational fluid dynamics. The results show that TSP has different stable angles and drag coefficients when having different wings. When the span increases to 250 mm from 150 mm, the stable angle decreases to 16° from 43°, the drag coefficient increases to 1.61 from 0.78. The stable angle, torque coefficient and drag coefficient are different when the aileron oblique angle increases to 60° from 0°. This paper can provide helpful references for design of non-fin TSP.

non-fin TSP; computational fluid dynamics; aerodynamic characteristics; aileron

2014-09-19

國家自然科學基金-青年科學項目(11102088);國家自然科學基金(11372136)資助

李恒(1993-),男,山東煙臺人,學士,研究方向:流體力學、靈巧彈藥空氣動力學。

TJ414.5

A