龍海斌，吳裕平

(中國直升機設計研究所，江西 景德鎮(zhèn) 333001)

0 引言

導彈、火箭彈等武器通常采用外掛的方式加裝在飛行器機翼下方，但近年來部分無人直升機采用在機身左右兩側外掛的方式來加裝武器，形成具備打擊敵方目標能力的武裝無人直升機。如美國的“火力偵查兵”采用掛梁+掛架的形式掛載了空地導彈和70mm火箭彈等武器。由武器與懸掛裝置等組成的組合外掛體在前飛時與機身、尾梁和平尾等部件產生氣動干擾，與其單獨在前飛流場中時的氣動特性差別比較大。目前國內外主要采用風洞試驗、經驗公式和數值計算等方法分別對飛行器機載武器和外掛物的氣動特性進行計算與分析。文獻[1]通過風洞試驗的方法研究了八棱錐頭導彈的減阻特性；文獻[2]在風洞中采用同時測力試驗技術對某型飛機進行了外掛測力試驗,獲得了該機翼下外掛物及其部件的氣動特性；文獻[3]則分別采用數值計算和試驗測量的方法計算和實測了直升機武器掛架的氣動載荷。經驗公式具有計算速度快的優(yōu)點，因此在型號初步設計階段應用得比較多。文獻[4]根據細長體理論干擾因子法得出了導彈部件升力的計算公式，編寫了可用于估算導彈氣動特性的軟件并應用于某型號導彈氣動特性計算；文獻[5]提出了一種適于在初步設計中計算導彈大攻角氣動特性的工程方法。隨著計算機技術和數值計算方法的進步，研究人員開始采用求解N-S方程的方法來計算導彈和外掛物的氣動特性。文獻[6]、[7] 采用求解N -S方程的紊流場計算軟件對某型空空導彈和有折翼的導彈進行了氣動特性計算。文獻[8]、[9]采用結構重疊網格技術對不同掛彈狀態(tài)的機翼-掛架-外掛物組合體進行了氣動特性計算與分析。文獻[10]采用格子Boltmann方法計算了稀薄氣體效應對導彈氣動特性的影響。從飛行試驗數據中也可以獲取導彈的氣動特性。文獻[11]采用基于K-最近鄰算法的數據分析方法從飛行試驗數據中獲取對導彈氣動特性的準確描述，以預測導彈的氣動特性。但是國內外對于無人直升機的組合外掛體的氣動特性研究比較少，本文采用求解N-S方程的CFD方法對某型無人直升機機載狀態(tài)的組合外掛體(掛梁、掛架、導發(fā)架和導彈四部分組成)的氣動特性進行了數值計算與分析研究。

1 數值計算概述

某型無人直升機的武器外掛方式與“火力偵查兵”類似，如圖1所示。組合外掛體由掛梁、掛架、導發(fā)架和導彈四部分組成，其中掛梁形狀與“火力偵查兵”的相似，掛架為“八”字形，導發(fā)架外形為前后有斜切面的長方體，導彈前緣為半球形，尾部有氣動控制面。

圖1 “火力偵查兵”武器示意圖

數值計算之前對無人直升機以及組合外掛體的幾何外形進行修理，之后將流場域進行網格劃分，對幾何外形曲率變化較大和流場比較復雜的區(qū)域進行了網格局部加密。數值計算過程中采用Roe格式的空間離散方法，湍流模型選擇一方程S-A模型。設置來流速度為55.57m/s(200km/h)以計算大速度前飛狀態(tài)下的氣動特性。

設定無人直升機抬頭時的攻角為正，向左側偏轉時側滑角為正。力矩系數參考點為計算狀態(tài)的重心位置。數值計算過程中的無人直升機模型包括機身、主槳轂、穩(wěn)瞄、尾梁、尾槳轂、垂尾、平尾等部件。由于大速度前飛時外掛體位于旋翼尾流之外，因此沒有考慮旋翼下洗流場的影響。該型無人直升機采用對稱的方式在機身兩側外掛武器，而且組合外掛體氣動外形相同。為了研究方便，只選取了位于機身右側的組合外掛體的CFD計算結果進行氣動特性分析，左側組合外掛體的氣動特性與右側基本相同。

2 組合外掛體氣動特性

某型無人直升機右側的組合外掛體在各個計算狀態(tài)下的氣動力與力矩系數如圖2所示。分析圖2(a)可以發(fā)現掛梁對組合外掛體阻力的貢獻最大，在0°攻角時掛梁的阻力占組合外掛體總阻力的46.40%，同時組合外掛體的阻力占全機阻力的7.66%。隨著攻角的增大，掛梁、掛架和導發(fā)架的阻力基本上保持不變，而導彈的阻力卻不斷增大，因此大攻角時組合外掛體的阻力增大主要是由導彈引起的。表1給出了各部件在0°攻角時的迎風面積，對比表1中的迎風面積與圖2(a)中的阻力系數值可以看出掛梁與導發(fā)架的迎風面積基本相等，但是阻力卻相差比較大，說明部件的氣動外形對阻力影響比較大。

表1 各部件迎風面積(0°攻角)

由圖2(b)各曲線變化趨勢可以看出，攻角變化時掛梁、掛架和導發(fā)架的升力變化比較小，而導彈的升力呈線性增長。其中在2°攻角時組合外掛體的升力系數與導彈的比值為1.0045，而大部分攻角范圍內導彈與組合外掛體的升力變化趨勢基本一致。

而分析圖2(c)可以發(fā)現，隨著攻角的增大,組合外掛體的俯仰力矩不斷增大，說明其在縱向上是靜不穩(wěn)定的。其中掛架和導發(fā)架的俯仰力矩系數與組合外掛體變化趨勢一致，而掛梁和導彈的俯仰力矩基本上不受攻角變化的影響。

機身側滑角變化主要引起組合外掛體的側向力、滾轉力矩和偏航力矩的變化。由圖2(d)中的曲線變化趨勢可以發(fā)現組合外掛體和各部件的側向力隨側滑角呈線性變化，其中掛梁與掛架的側向力比較小，導彈的側向力比較大。分析圖2(e)可以看出組合外掛體和導彈的滾轉力矩隨著側滑角的增大而不斷減小，說明它們在橫向上是靜穩(wěn)定的。在-4°到2°攻角范圍內，導彈的滾轉力矩與組合外掛體滾轉力矩的比值在0.95～1.17之間，兩者的差別比較小。而導發(fā)架的滾轉力矩系數在大側滑角時不斷減小。從圖2(f)可以看出組合外掛體偏航力矩系數隨著側滑角的增大而不斷增長，說明組合外掛體在航向上是靜不穩(wěn)定的，其中組合外掛體的偏航力矩線斜率與導彈偏航力矩線斜率的比值為1.26，但兩者隨攻角的變化趨勢一致。

圖2 組合外掛體及各部件氣動特性圖

3 發(fā)射安裝角變化對氣動特性的影響

無人直升機載導彈的發(fā)射安裝角對導彈的發(fā)射距離和組合外掛體的氣動特性等有比較大的影響，導彈安裝角變化時組合外掛體的氣動特性如圖3所示。分析圖3(a)可以發(fā)現，導彈發(fā)射安裝角變化時組合外掛體的阻力變化趨勢基本上一致，在阻力系數變化圖上相當于偏移了一定的攻角值。而從圖3(b)可以看出在0°至9°范圍內，發(fā)射安裝角每增大3°，同一攻角的升力系數增大約0.006，組合外掛體在同一攻角時的升力增長與發(fā)射安裝角變化基本上呈線性變化關系。由圖3(c)可以看出，組合外掛體的俯仰力矩變化趨勢與升力基本一致。

分析圖3(d)可以發(fā)現，發(fā)射安裝角變化時組合外掛體的側向力基本保持一致，由于發(fā)射安裝角變化時側向力迎風面積不變，因此側向力基本不變。由圖3(e)可以看出，發(fā)射安裝角每增大3°，組合外掛體的滾轉力矩系數減小約0.001，同一攻角時的滾轉力矩隨發(fā)射安裝角的變化而呈線性變化趨勢。根據圖3(f)可以看出在負側滑角范圍內組合外掛體的偏航力矩呈線性變化趨勢；而在正側滑角范圍內，偏航力矩的差別隨發(fā)射安裝角的增大逐漸變小。

4 機身側滑角變化對氣動特性的影響

無人直升機在轉彎飛行或遇到側向陣風(相當于機身側滑角發(fā)生變化)時組合外掛體的氣動特性如圖4所示。

分析圖4(a)可以看出，側滑角從-6°變化到0°時，組合外掛體的阻力變化比較小，而從0°變化到6°時阻力不斷增大。在0°攻角時，-6°側滑角時的阻力比0°時增大1.4%，而6°側滑角的阻力比0°時增大14.0%。這是由于選取的是位于機身右側的組合外掛體的氣動特性數據，因此負側滑角時機身對側風的干擾產生了阻礙作用。分析圖4(b)中的數據可以發(fā)現，在0°攻角時，-6°側滑角時的升力比0°時增大12.9%，而6°側滑角的升力比0°時減小20.7%，說明側滑角變化時組合外掛體的升力系數變化比較大，其中正側滑角范圍內影響更大。而由圖4(c)可以看出在0°攻角時，-6°側滑角時的俯仰力矩(絕對值)比0°時減小11.2%，而6°側滑角的俯仰力矩(絕對值)比0°時增大15.1%，組合外掛體的俯仰力矩系數的曲線斜率大小基本一致。由于組合外掛體的俯仰力矩為負值，因此正側滑角時組合外掛體的俯仰力矩變化更大。力和力矩隨機身側滑角的變化趨勢說明由于組合外掛體安裝在機身左右外側，側風或轉彎飛行時組合外掛體對無人直升機的升阻力和穩(wěn)定性影響比較大。

圖4 不同側滑角時組合外掛體氣動特性圖

5 掛載狀態(tài)變化對氣動特性的影響

武裝無人直升機在導彈發(fā)射之后需要返回基地進行加油、維修和重新裝彈等，因此無掛載飛行狀態(tài)的時間占無人直升機使用時間的比例比較大。組合外掛體無掛載狀態(tài)時的氣動特性如圖5所示。分析圖5(a)可以看出，在小攻角范圍內，無掛載狀態(tài)的阻力約為掛載狀態(tài)阻力的80%。而隨著攻角的增大，掛載狀態(tài)的阻力增長更快。由圖5(b)可以發(fā)現，無掛載與掛載狀態(tài)的升力系數隨攻角變化呈線性關系，其中掛載狀態(tài)的升力線斜率是無掛載狀態(tài)升力線斜率的2.37倍。

從圖5(c)可以看出，俯仰力矩系數的變化趨勢與升力系數變化趨勢一致，而掛載狀態(tài)的俯仰力矩線斜率是無掛載狀態(tài)俯仰力矩線斜率的1.66倍，說明掛載狀態(tài)的俯仰力矩變化更快，對機身的縱向穩(wěn)定性影響也越大。分析圖5(d)和圖5(e)可以發(fā)現，無掛載與掛載狀態(tài)的側向力系數和滾轉力矩系數與側滑角都呈線性變化關系，其中掛載狀態(tài)的側向力線斜率是無掛載狀態(tài)側向力線斜率的2.53倍，同時掛載狀態(tài)的滾轉力矩線斜率是無掛載狀態(tài)滾轉力矩線斜率的2.57倍。掛載狀態(tài)的側向力系數與滾轉力矩系數曲線斜率更小，說明掛載狀態(tài)組合外掛體對無人直升機橫向穩(wěn)定性影響更大。從圖5(f)可以看出，無掛載狀態(tài)的偏航力矩系數基本上不受側滑角變化的影響，而掛載狀態(tài)的偏航力矩系數隨側滑角的增大而呈線性增大的趨勢。由于導彈的側向迎風面積比較大，因此無掛載狀態(tài)的側向力、滾轉和偏航力矩與掛載導彈狀態(tài)差別比較大。

圖5 不同掛載狀態(tài)時組合外掛體氣動特性圖

6 結論與探討

通過對某型無人直升機組合外掛體在不同導彈發(fā)射安裝角、機身側滑角以及不同掛載狀態(tài)時的氣動特性進行計算與分析，得出如下結論：

1)在負攻角與小攻角范圍內，掛梁對組合外掛體阻力的貢獻比較大。組合外掛體其他氣動參數的變化趨勢與導彈基本一致，說明導彈對整個組合外掛體的氣動特性影響很大。這為組合外掛體設計以及減阻等提供了研究方向。

2)無人直升機組合外掛體的升力、俯仰力矩和滾轉力矩隨導彈發(fā)射安裝角的改變呈線性變化關系。而阻力、側向力基本上不受導彈發(fā)射安裝角變化的影響。在負側滑角和小側滑角時，導彈發(fā)射安裝角對偏航力矩影響相對比較大。

3)由于機身對側風的阻礙作用，機身側滑角變化時，組合外掛體的氣動特性變化比較大。

4)導彈發(fā)射后組合外掛體的氣動特性變化趨勢與掛載狀態(tài)基本一致，掛載狀態(tài)的曲線斜率更大，說明掛載狀態(tài)時組合外掛體對無人直升機平衡性等影響更大。

5)通過對組合外掛體氣動特性的計算與分析，可根據分析結果來選擇類似外掛方式與氣動外形武器的無人直升機估算氣動特性。在已知無人直升機載武器的氣動特性之后，可以初步估算類似組合外掛體的氣動特性，以確定該外掛方式是否可行；若此方案可行，則可以初步確定武器發(fā)射安裝角。

參考文獻：

[1] 邢秀華. 八棱錐頭導彈減阻特性研究[J]. 航空兵器，1996(5):16-17.

[2] 李熙佩. 飛機外掛物及其部件同時測力試驗技術研究 [J]. 實驗流體力學，2005,19(3):7-9.

[3] 鄭甲宏，沈 靂，寇富軍. 直升機武器掛架飛行載荷研究[J]. 航空科學技術,2016,27(02): 60-64.

[4] 校金友,張 鐸. 導彈氣動特性與載荷計算[J]. 彈箭與制導學報，2004,24(2):62-64.

[5] 董月娟. 導彈大攻角氣動特性計算[J]. 戰(zhàn)術導彈技術，1995(4):10-17.

[6] 趙洪章,岳春國,李進賢. 基于Fluent的導彈氣動特性計算. 彈箭與制導學報,2007,27(2):203-205.

[7] 杜 廈,昂海松.一種改進折疊彈翼結構設計與氣動分析[J]. 航空兵器，2009(3):7-9.

[8] 王立強，董國國. 外掛物干擾流場特性數值仿真研究[J]. 航空計算技術，2012,42(1):57-60.

[9] 王立強，敬代勇.組合掛載外掛物氣動特性數值計算[J].戰(zhàn)術導彈技術，2012(6):16-21.

[10] 阮政委，何志強，周文雅, 等. 稀薄氣體效應對常規(guī)布局導彈氣動特性的影響[J]. 航空兵器，2016(5):3-7.

[11] 鄭鹍鵬，李海峰，華建林. 基于K-最近鄰的導彈靶試數據分析方法[J].航空兵器，2011(3):7-10.