張東升 魏江鵬 文程祥

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

關(guān)鍵詞： 前掠翼； 平直翼； 后掠翼； 氣動(dòng)分析； 失速； Solidworks

中圖分類(lèi)號(hào)： TN876?34 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 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

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

1 模型的建立與網(wǎng)格劃分

機(jī)翼翼展均為1 600 mm，現(xiàn)取半個(gè)機(jī)翼進(jìn)行流體計(jì)算。翼根部弦長(zhǎng)300 mm，翼尖部弦長(zhǎng)200 mm，前掠翼前緣前掠角為17°，后掠翼前緣后掠角為7°。計(jì)算區(qū)域采用Gambit劃分的非結(jié)構(gòu)網(wǎng)格，劃分網(wǎng)格時(shí)進(jìn)行了區(qū)域加密，對(duì)機(jī)翼周?chē)鲌?chǎng)的網(wǎng)格進(jìn)一步加密。圖2為網(wǎng)格放大示意圖。

計(jì)算域的邊界包括來(lái)流入口、來(lái)流出口、翼型固壁、Interior和固壁邊界。計(jì)算域前部外邊界長(zhǎng)度為機(jī)翼最大弦長(zhǎng)的10倍，后部外邊界長(zhǎng)度為最大弦長(zhǎng)的20倍，厚度邊界為機(jī)翼最大厚度的20倍，軸向邊界為翼長(zhǎng)的10倍。計(jì)算采用的邊界條件：固壁表面為無(wú)滑移條件[8]、進(jìn)口設(shè)置為pressure_far_field、出口設(shè)置為presssure_outlet、計(jì)算邊界由差值確定。

計(jì)算殘差收斂精度為10-3，來(lái)流馬赫數(shù)設(shè)置為0.08，雷諾數(shù)設(shè)置為4.997×105。計(jì)算基于RANS方法，采用二階迎風(fēng)格式作為空間離散格式，采用的隱式解法為二維穩(wěn)態(tài)分離解法，采用SIMPLE解法求解壓力?速度耦合[9?15]。

2 計(jì)算結(jié)果分析

氣動(dòng)性能的計(jì)算結(jié)果隨著迎角的增加如圖3所示。

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

2.1 壓力分析

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

2.2 切應(yīng)力分析

分析附面層氣流分離機(jī)理可知，越大的逆壓梯度，附面層分離越嚴(yán)重。此時(shí)壁面切應(yīng)力消失的點(diǎn)就是分離點(diǎn)，可以通過(guò)壁面切應(yīng)力的x分量是否為負(fù)來(lái)判斷。

圖5為三種翼型的三個(gè)z方向剖面壁面切應(yīng)力的x分量分布圖。由圖可以看出，平直翼沒(méi)有出現(xiàn)x分量為負(fù)值的情況，即沒(méi)有出現(xiàn)附面層分離現(xiàn)象，而前掠翼和后掠翼在翼身z=400剖面上出現(xiàn)了x分量為負(fù)的情況，說(shuō)明上吸力面出現(xiàn)輕微的氣流分離情況。

2.3 失速分析

由圖3a）可見(jiàn)，當(dāng)此款前掠翼的來(lái)流迎角處于0°～22°區(qū)間內(nèi)時(shí)，翼型的升力系數(shù)一直增大，這是隨著翼型迎角的增大，翼型下表面受力增大，而上表面受力變小導(dǎo)致的；當(dāng)翼型來(lái)流迎角在22°之后，翼型的升力系數(shù)減小，這是由于22°為此規(guī)格翼型的失速迎角，當(dāng)來(lái)流迎角大于22°時(shí)，翼型上表面氣流脫離形成分離渦，引起翼型升力減小。同時(shí)由圖3a）可以看出平直翼的失速迎角為22°，后掠翼的時(shí)速迎角也為22°。計(jì)算結(jié)果與三種翼型的空氣動(dòng)力學(xué)特性曲線(xiàn)吻合。

3 結(jié) 論

本文通過(guò)數(shù)值計(jì)算的方法，對(duì)機(jī)翼的三種布局，即前掠翼布局、平直翼布局和后掠翼布局的氣動(dòng)特性進(jìn)行計(jì)算。分析三種機(jī)翼的氣動(dòng)特性與流動(dòng)機(jī)理，對(duì)機(jī)翼流場(chǎng)進(jìn)行研究之后可知：

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

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

3） 平直翼升力系數(shù)與后掠翼近似，但平直翼沒(méi)有氣流分離情況，平直翼布局更適合文中小型無(wú)人機(jī)。

參考文獻(xiàn)

[1] 崔秀敏，王維軍，方振平.小型無(wú)人機(jī)發(fā)展現(xiàn)狀及其相關(guān)問(wèn)題分析[J].飛行力學(xué)，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] 劉文法，王旭，米康.一種新的變前掠翼無(wú)人機(jī)氣動(dòng)布局概念研究[J].航空學(xué)報(bào)，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] 郭劍東，宋彥國(guó).小型無(wú)人傾轉(zhuǎn)旋翼機(jī)氣動(dòng)與操縱特性試驗(yàn)研究[J].空氣動(dòng)力學(xué)學(xué)報(bào)，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ī)翼失速特性控制研究[J].空氣動(dòng)力學(xué)學(xué)報(bào)，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翼型失速特性的二維數(shù)值研究[J].內(nèi)蒙古工業(yè)大學(xué)學(xué)報(bào)，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] 王旭，黃萌，任智靜，等.前掠翼與平直翼布局氣動(dòng)特性的比較分析[J].空軍工程大學(xué)學(xué)報(bào)（自然科學(xué)版），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] 胡丹梅，杜朝輝，朱春建.水平軸風(fēng)力機(jī)靜態(tài)失速特性[J].太陽(yáng)能學(xué)報(bào)，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].空氣動(dòng)力學(xué)學(xué)報(bào)，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] 王福軍.計(jì)算流體動(dòng)力學(xué)分析CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社，2004.

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

[12] 袁新，江學(xué)忠.翼型大攻角低速分離流動(dòng)的數(shù)值模擬[J].工程熱物理學(xué)報(bào)，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] 方寶瑞，李天，余松濤，等.飛機(jī)氣動(dòng)布局設(shè)計(jì)[M].北京：航空工業(yè)出版社，1997.

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

[14] 劉文法，王旭，米康.前掠翼與后掠翼布局流動(dòng)機(jī)理的數(shù)值研究[J].空軍工程大學(xué)學(xué)報(bào)（自然科學(xué)版），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] 任智靜，王旭，劉文法.前掠翼布局中鴨翼氣動(dòng)影響的數(shù)值模擬[J].航空學(xué)報(bào)，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.