基于Fluent的某小型無人機機翼布局對比研究

張東升 魏江鵬 文程祥

摘 要： 運用Solidworks對某規格的前掠翼、平直翼和后掠翼進行三維建模，利用CFD軟件，采用三維N?S方程及Spalart?Allmaras渦粘湍流模型對前掠翼、平直翼和后掠翼的空氣動力學特性進行研究，每間隔2°計算迎角從0°～38°時的不同情況，對比分析各自的優勢和不足，最后給出結論，為低速小型無人機的機翼布局選型提供了理論依據。計算結果表明低速飛行下平直翼布局升力系數較大，沒有氣流分離的情況，這種布局較為合理。

關鍵詞： 前掠翼； 平直翼； 后掠翼； 氣動分析； 失速； Solidworks

中圖分類號： TN876?34 文獻標識碼： A 文章編號： 1004?373X（2018）08?0147?04

Abstract： Solidworks is used to build 3D models for forward swept wing， straight wing and backward swept wing of a certain specification. The CFD software， three dimensional N?S equation and Spalart?Allmaras eddy viscosity turbulence model are adopted to study the aerodynamic characteristics of forward swept wing， straight wing and backward swept wing. Angles of attack from 0° to 38° are calculated at an interval of every 2°， whose advantages and disadvantages are contrasted and analyzed so that the conclusion can be finally drawn to provide a theoretical basis for wing layout selection of a low?speed small UAV. The calculation results show that the lift coefficient of straight wing layout is larger at a low speed flight as there exists no air separation， and the straight wing layout is more reasonable.

Keywords： forward swept wing； straight wing； backward swept wing； aerodynamic analysis； stall； Solidworks

無人機可完成某些載人飛機無法完成的特殊任務，使得無人機的發展受到空前的關注，各個國家紛紛投入大量資源研制新型多功能無人機 [1?3]。對于小型無人機等低速飛行的飛行器而言，機翼布局是一項值得研究的課題。現設計一款小型傾轉旋翼無人機，同時具有旋翼與固定翼，為了確定固定翼布局，利用Fluent對不同機翼進行空氣動力學仿真，進一步探討機翼布局氣動特性[4?6]，通過對比結果選出合適的機翼布局。NACA翼型有著良好的升力特性[7]，本文將對基于NACA4415翼型的前掠翼、平直翼與后掠翼進行升力特性對比，給出小型無人機機翼布局的設計參考，翼型如圖1所示。

1 模型的建立與網格劃分

機翼翼展均為1 600 mm，現取半個機翼進行流體計算。翼根部弦長300 mm，翼尖部弦長200 mm，前掠翼前緣前掠角為17°，后掠翼前緣后掠角為7°。計算區域采用Gambit劃分的非結構網格，劃分網格時進行了區域加密，對機翼周圍流場的網格進一步加密。圖2為網格放大示意圖。

計算域的邊界包括來流入口、來流出口、翼型固壁、Interior和固壁邊界。計算域前部外邊界長度為機翼最大弦長的10倍，后部外邊界長度為最大弦長的20倍，厚度邊界為機翼最大厚度的20倍，軸向邊界為翼長的10倍。計算采用的邊界條件：固壁表面為無滑移條件[8]、進口設置為pressure_far_field、出口設置為presssure_outlet、計算邊界由差值確定。

計算殘差收斂精度為10-3，來流馬赫數設置為0.08，雷諾數設置為4.997×105。計算基于RANS方法，采用二階迎風格式作為空間離散格式，采用的隱式解法為二維穩態分離解法，采用SIMPLE解法求解壓力?速度耦合[9?15]。

2 計算結果分析

氣動性能的計算結果隨著迎角的增加如圖3所示。

三種翼型的升力系數CL有著微小的的差距，整體趨勢相同，先增長再下降；三種翼型的阻力系數CD幾乎沒有差別，隨著迎角的增加持續增加。相同升力系數下，低于失速迎角，前掠翼阻力系數最大；高于失速迎角，前掠翼阻力系數最小，平直翼與后掠翼升阻比相似。

2.1 壓力分析

圖4為22°迎角下三種翼型在z方向不同位置翼剖面的上下表面靜壓分布曲線。由圖可以看出：在翼根z=50處，后掠翼上翼面吸力比前掠翼和平直翼稍大；在翼尖z=750處，則是前掠翼的吸力比較大。這是由于三種翼型氣流展向速度方向不一致，導致附面層在不同的位置堆積，產生氣流分離引起的差異。

2.2 切應力分析

分析附面層氣流分離機理可知，越大的逆壓梯度，附面層分離越嚴重。此時壁面切應力消失的點就是分離點，可以通過壁面切應力的x分量是否為負來判斷。

圖5為三種翼型的三個z方向剖面壁面切應力的x分量分布圖。由圖可以看出，平直翼沒有出現x分量為負值的情況，即沒有出現附面層分離現象，而前掠翼和后掠翼在翼身z=400剖面上出現了x分量為負的情況，說明上吸力面出現輕微的氣流分離情況。

2.3 失速分析

由圖3a）可見，當此款前掠翼的來流迎角處于0°～22°區間內時，翼型的升力系數一直增大，這是隨著翼型迎角的增大，翼型下表面受力增大，而上表面受力變小導致的；當翼型來流迎角在22°之后，翼型的升力系數減小，這是由于22°為此規格翼型的失速迎角，當來流迎角大于22°時，翼型上表面氣流脫離形成分離渦，引起翼型升力減小。同時由圖3a）可以看出平直翼的失速迎角為22°，后掠翼的時速迎角也為22°。計算結果與三種翼型的空氣動力學特性曲線吻合。

3 結 論

本文通過數值計算的方法，對機翼的三種布局，即前掠翼布局、平直翼布局和后掠翼布局的氣動特性進行計算。分析三種機翼的氣動特性與流動機理，對機翼流場進行研究之后可知：

1） 在低速飛行狀態下，平直翼與后掠翼的升力系數相似，比前掠翼稍大，在0°迎角下三種翼型升力系數相似，隨著迎角增大，平直翼與后掠翼的升力系數曲線貼合，始終大于前掠翼。

2） 三種翼型的時速迎角相同，均為22°，符合機翼迎角范圍。

3） 平直翼升力系數與后掠翼近似，但平直翼沒有氣流分離情況，平直翼布局更適合文中小型無人機。

參考文獻

[1] 崔秀敏，王維軍，方振平.小型無人機發展現狀及其相關問題分析[J].飛行力學，2005，23（1）：14?18.

CUI Xiumin， WANG Weijun， FANG Zhenping. Present situation and some problems analysis of small?size unmanned air vehicles [J]. Flight dynamics， 2005， 23（1）： 14?18.

[2] JOHNSON W， YAMAUCHI G K， DERBY M R， et al. Wind tunnel measurements and calculations of aerodynamic interactions between tiltrotor aircraft [C]// Proceedings of 41st Aerospace Sciences Meeting and Exhibit. Reno： America Institute of Aeronautics and Astronautics， 2003： 1?10.

[3] ZHOU J， WEN C， YUN Y， et al. The history of the XV?15 tilt rotor research aircraft： from concept to flight [J]. Journal of aircraft， 2015， 8（1）： 1?12.

[4] 劉文法，王旭，米康.一種新的變前掠翼無人機氣動布局概念研究[J].航空學報，2009，30（5）：832?836.

LIU Wenfa， WANG Xu， MI Kang. A new aerodynamic configuration of UAV with variable forward?swept wing [J]. Acta aeronautica et astronautica sinica， 2009， 30（5）： 832?836.

[5] 郭劍東，宋彥國.小型無人傾轉旋翼機氣動與操縱特性試驗研究[J].空氣動力學學報，2015，33（1）：107?112.

GUO Jiandong， SONG Yanguo. Testing study on aerodynamics and control characteristics of a small unmanned tilt rotor [J]. Acta aerodynamica sinica， 2015， 33（1）： 107?112.

[6] 白亞磊，李鵬.有限翼展機翼失速特性控制研究[J].空氣動力學學報，2013，31（1）：9?14.

BAI Yalei， LI Peng. Researches of a finite?span wing stall control [J]. Acta aerodynamica sinica， 2013， 31（1）： 9?14.

[7] 曲立群，汪建文，朱德臣，等.NACA4415翼型失速特性的二維數值研究[J].內蒙古工業大學學報，2009，28（1）：48?51.

QU Liqun， WANG Jianwen， ZHU Dechen， et al. A two?dimensional numerical study on stall of NACA4415 airfoil [J]. Journal of Inner Mongolia University of Technology， 2009， 28（1）： 48?51.

[8] 王旭，黃萌，任智靜，等.前掠翼與平直翼布局氣動特性的比較分析[J].空軍工程大學學報（自然科學版），2011，12（4）：1?4.

WANG Xu， HUANG Meng， REN Zhijing， et al. Comparison and analysis on the aerodynamic characteristics of forward swept and orthogonal wing configurations [J]. Journal of Air Force Engineering University （Natural science edition）， 2011， 12（4）： 1?4.

[9] 胡丹梅，杜朝輝，朱春建.水平軸風力機靜態失速特性[J].太陽能學報，2006，27（3）：217?222.

HU Danmei， DU Zhaohui， ZHU Chunjian. A study on the static stall for horizontal axis wind turbine [J]. Acta energiae solaris sinica， 2006， 27（3）： 217?222.

[10] 蘇繼川，黃勇，李永紅，等.小展弦比飛翼亞、跨、超聲速支撐干擾研究[J].空氣動力學學報，2015，33（3）：289?295.

SU Jichuan， HUANG Yong， LI Yonghong， et al. Support interference of low?aspect?ratio flying?wing from subsonic to supersonic speed [J]. Acta aerodynamica sinica， 2015， 33（3）： 289?295.

[11] 王福軍.計算流體動力學分析CFD軟件原理與應用[M].北京：清華大學出版社，2004.

WANG Fujun. Principle and application of computational fluid dynamics analysis CFD software [M]. Beijing： Tsinghua University Press， 2004.

[12] 袁新，江學忠.翼型大攻角低速分離流動的數值模擬[J].工程熱物理學報，1999，20（2）：161?165.

YUAN Xin， JIANG Xuezhong. Simulation of low?speed separated flows for the NREL airfoil at high angles of attack [J]. Journal of engineering thermophysics， 1999， 20（2）： 161?165.

[13] 方寶瑞，李天，余松濤，等.飛機氣動布局設計[M].北京：航空工業出版社，1997.

FANG Baorui， LI Tian， YU Songtao， et al. Aircraft aerodynamic layout design [M]. Beijing： Aviation Industry Press， 1997.

[14] 劉文法，王旭，米康.前掠翼與后掠翼布局流動機理的數值研究[J].空軍工程大學學報（自然科學版），2008，9（6）：11?15.

LIU Wenfa， WANG Xu， MI Kang. The numerical research on flow mechanism of forward?swept wing and backward?swept wing configurations [J]. Journal of Air Force Engineering University （Natural science edition）， 2008， 9（6）： 11?15.

[15] 任智靜，王旭，劉文法.前掠翼布局中鴨翼氣動影響的數值模擬[J].航空學報，2010，31（7）：1318?1323.

REN Zhijing， WANG Xu， LIU Wenfa. Numerical simulation on aerodynamic influence of canard on forward?swept configuration [J]. Acta aeronautica et astronautica sinica， 2010， 31（7）： 1318?1323.