某運(yùn)輸機(jī)加裝失速條氣動(dòng)特性研究

劉毅，趙曉霞，歐陽紹修，袁志敏

(中航飛機(jī)股份有限公司研發(fā)中心，西安710089)

某運(yùn)輸機(jī)加裝失速條氣動(dòng)特性研究

劉毅*，趙曉霞，歐陽紹修，袁志敏

(中航飛機(jī)股份有限公司研發(fā)中心，西安710089)

為改善某運(yùn)輸機(jī)著陸襟翼構(gòu)型失速急劇滾轉(zhuǎn)問題，采用數(shù)值計(jì)算和風(fēng)洞實(shí)驗(yàn)方法優(yōu)選了機(jī)翼失速條的外形參數(shù)，并對(duì)氣動(dòng)力和流場(chǎng)特性進(jìn)行了研究分析。以失速條高度H和安裝位置距離前緣的長(zhǎng)度S為設(shè)計(jì)變量，采用求解RANS方程的方法研究了失速條對(duì)著陸構(gòu)型翼型二維特性的影響，表明S越小(即越靠近上翼面)失速迎角提前越多，H增大也能使失速迎角提前但敏感性小于S。失速條后方產(chǎn)生了分離氣泡且隨迎角增加而逐漸增大增長(zhǎng)，在破裂后導(dǎo)致翼型失速提前，使升力線出現(xiàn)圓弧形的失速特征。設(shè)計(jì)了4種失速條在機(jī)翼上的平面布局方案，通過縮比模型風(fēng)洞實(shí)驗(yàn)驗(yàn)證表明:40%半展長(zhǎng)處展向長(zhǎng)度2m，S=0的失速條使升力線由急劇失速變?yōu)槠巾斝褪俨⑾耸俸蟮牟粚?duì)稱滾轉(zhuǎn)力矩，將此失速條展長(zhǎng)縮小一半的2種方案也不同程度地改善了失速形態(tài)，15%半展長(zhǎng)處失速條對(duì)失速特性無明顯改善，主要原因是氣流分離從約40%半展長(zhǎng)處開始發(fā)生，失速條安裝在這一展向位置時(shí)才能發(fā)揮作用。

失速條;失速特性;急劇滾轉(zhuǎn);分離氣泡;計(jì)算流體力學(xué);風(fēng)洞實(shí)驗(yàn)

Keywords：stall strips；stall characteristics；violent roll motion；separation bubble；CFD；wind tunnel test

0 引 言

運(yùn)輸類飛機(jī)良好的失速特性是保證飛行安全的基礎(chǔ)，其中的重要參數(shù)之一是失速及改出過程中飛機(jī)出現(xiàn)的滾轉(zhuǎn)角度，根據(jù)中國民用航空規(guī)章第25部要求：“對(duì)于機(jī)翼水平失速，在失速和完成改出之間發(fā)生的滾轉(zhuǎn)大約不得超過20°左右”［1］。飛機(jī)的失速特性主要受機(jī)翼影響，緩和的失速特性通常比最大升力系數(shù)的量值更為重要，避免外翼失速是控制飛機(jī)失速滾轉(zhuǎn)并保持大迎角副翼操縱效率的有效手段，最佳的失速起始位置一般在40%半展長(zhǎng)處［2－4］。由于失速特性難以通過數(shù)值計(jì)算或風(fēng)洞實(shí)驗(yàn)準(zhǔn)確預(yù)測(cè)，往往在飛機(jī)進(jìn)入試飛階段后才發(fā)現(xiàn)失速特性問題，此時(shí)盡量采取局部改進(jìn)設(shè)計(jì)的方式來進(jìn)行調(diào)整，例如采用翼型前緣下垂、前緣縫翼、翼刀和失速條等。在外翼加裝前緣縫翼或翼型下垂使外翼氣流分離推遲能夠獲得更大的升力系數(shù)，但設(shè)計(jì)改動(dòng)量較大且仍然存在左右機(jī)翼失速不對(duì)稱的可能性。

失速條是一種安裝在機(jī)翼前緣的固定式流動(dòng)控制裝置，一般對(duì)稱安裝在左右機(jī)翼靠?jī)?nèi)側(cè)，強(qiáng)制使左右內(nèi)側(cè)機(jī)翼提前對(duì)稱地分離失速，從而獲得滿意的失速特性［5－8］。失速條的發(fā)現(xiàn)源于結(jié)冰翼型失速提前并具有緩和形態(tài)的現(xiàn)象，其應(yīng)用可追溯至20世紀(jì)30年代，但其設(shè)計(jì)研究通常采用試湊法，沒有嚴(yán)格的設(shè)計(jì)準(zhǔn)則和依據(jù)［9－10］。20世紀(jì)50～60年代德哈維蘭公司的彗星4C和DH125飛機(jī)均加裝了失速條；70年代在S－3A反潛機(jī)上加裝失速條的研究表明，失速條保證了左右機(jī)翼大迎角時(shí)的流場(chǎng)對(duì)稱性，值得注意的是上述3型飛機(jī)的失速條均加裝在機(jī)翼40%半展長(zhǎng)附近［5］。20世紀(jì)90年代美國的F／A－18E／F戰(zhàn)機(jī)在跨聲速機(jī)動(dòng)時(shí)出現(xiàn)了急劇滾轉(zhuǎn)問題，為此NASA系統(tǒng)地研究了通過CFD或風(fēng)洞實(shí)驗(yàn)預(yù)測(cè)此類問題的可行性，表明傳統(tǒng)的靜態(tài)升力特性和CFD結(jié)果能夠與失速滾轉(zhuǎn)特性相關(guān)［11－15］。

某運(yùn)輸機(jī)在失速試飛過程中發(fā)現(xiàn)，隨著襟翼偏度增加，失速后出現(xiàn)的滾轉(zhuǎn)角度逐漸增大，在著陸襟翼構(gòu)型時(shí)滾轉(zhuǎn)角通常大于20°，不滿足適航標(biāo)準(zhǔn)要求。前期研究表明，縮比模型靜態(tài)測(cè)力試驗(yàn)中升力線失速的急劇程度以及失速區(qū)域滾轉(zhuǎn)力矩系數(shù)的偏離情況與試飛現(xiàn)象中滾轉(zhuǎn)的強(qiáng)弱有正比關(guān)系，若使升力線失速緩和、滾轉(zhuǎn)力矩在失速區(qū)域偏離減小，有可能解決試飛中的失速急劇滾轉(zhuǎn)問題［16］。為了以盡量小的代價(jià)改善飛機(jī)著陸構(gòu)型失速特性，提高飛機(jī)安全性，采用失速條方案是較為簡(jiǎn)單有效的方法。

1 二維特性數(shù)值計(jì)算研究

1．1 失速條二維參數(shù)設(shè)計(jì)

借鑒國內(nèi)外飛機(jī)失速條的設(shè)計(jì)經(jīng)驗(yàn)，失速條二維剖面外形通常為尖銳的三角形，安裝在翼型前緣附近使其失速提前。本文所研究的失速條設(shè)計(jì)工作點(diǎn)為著陸襟翼高升力狀態(tài)，前緣流場(chǎng)特性相對(duì)巡航構(gòu)型有較大變化，對(duì)失速條的高度和剖面安裝位置的設(shè)計(jì)提出了挑戰(zhàn)。

將失速條的2個(gè)關(guān)鍵幾何參數(shù)即高度H和安裝位置S（定義見圖1）在可能的取值范圍內(nèi)組成參數(shù)矩陣進(jìn)行了系統(tǒng)地分析研究，失速條前緣角固定為60°。對(duì)參數(shù)H范圍20～40mm，S范圍－60～60mm進(jìn)行組合形成了多個(gè)失速條方案，從中得到了可行的失速條二維參數(shù)，同時(shí)獲得了以上2個(gè)參數(shù)對(duì)著陸襟翼二維翼型失速特性的影響規(guī)律。

圖1 失速條參數(shù)H、S定義Fig．1 Definition of stall strip parameters H and S

1．2 數(shù)值計(jì)算方法

數(shù)值計(jì)算采用了基于有限體積法求解RANS方程的方法進(jìn)行。采用結(jié)構(gòu)化網(wǎng)格對(duì)著陸構(gòu)型襟翼二維剖面及失速條的繞流流域進(jìn)行了離散，物面第一層網(wǎng)格高度在10－5m量級(jí)以滿足y＋≈1的條件，對(duì)流場(chǎng)參數(shù)變化劇烈的區(qū)域進(jìn)行了網(wǎng)格加密，相鄰網(wǎng)格間的尺度比例小于1．2，總網(wǎng)格數(shù)在6×105左右，網(wǎng)格示意圖如圖2所示。不同失速條剖面方案僅對(duì)前緣局部網(wǎng)格進(jìn)行適應(yīng)性修改，使各方案網(wǎng)格條件盡可能對(duì)等。計(jì)算域在翼型的前方和上下為30倍弦長(zhǎng)，在后方為45倍弦長(zhǎng)，邊界條件均為速度入口。數(shù)值求解在FLUENT中完成，采用壓力與速度同時(shí)求解的耦合解法，對(duì)流項(xiàng)采用二階迎風(fēng)格式進(jìn)行離散，湍流模型為一方程的SA模型，渦粘性采用一階迎風(fēng)格式離散。

為驗(yàn)證建模和計(jì)算方法的可靠性，對(duì)相似構(gòu)型的某多段翼型二元特性進(jìn)行了計(jì)算并與實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比（見圖3），可見在升力失速以前計(jì)算值與實(shí)驗(yàn)值吻合良好，對(duì)失速迎角和最大升力系數(shù)2個(gè)關(guān)鍵參數(shù)也給出了較好的估計(jì)值，但對(duì)失速后的升力線形態(tài)未能準(zhǔn)確模擬。由于加裝失速條后誘發(fā)的前緣氣泡現(xiàn)象及氣動(dòng)力拐點(diǎn)均出現(xiàn)在失速迎角之前，因此所采用的計(jì)算方法是較為合理可信的，可以作為失速條初步選型分析的工具。

圖2 帶失速條翼型的二維網(wǎng)格Fig．2 Two dimensional grids of the airfoil section with stall strip

圖3 多段翼型升力特性計(jì)算值與實(shí)驗(yàn)值對(duì)比Fig．3 The comparison between calculation and test of the lift of a multi－element airfoil

1．3 失速條對(duì)翼型二維氣動(dòng)特性的影響

飛機(jī)著陸構(gòu)型加裝不同失速條方案的二元升力特性如圖4和5所示。從計(jì)算結(jié)果可見，各失速條方案使原始翼型的失速迎角有不同程度的提前，多數(shù)情況下失速形態(tài)從急劇失速改變?yōu)閳A弧形的緩和失速，失速迎角隨參數(shù)H、S有規(guī)律地變化。

失速條距離前緣位置參數(shù)S對(duì)失速特性的影響很大，隨著S減?。粗饾u向上翼面移動(dòng)）失速迎角顯著降低，S＝0mm時(shí)二元剖面失速迎角約為3°，而原始構(gòu)型失速迎角約為14°。失速條高度H對(duì)升力失速迎角也有影響，H越大失速提前越多，但高度20～40mm之間的失速條均能使失速提前。計(jì)算結(jié)果與理論分析的規(guī)律是吻合的，即失速條越靠近上翼面，高度越大，觸發(fā)氣流分離的迎角越早，且氣流分離發(fā)生前對(duì)升力特性的影響不大。計(jì)算獲得的定量結(jié)果還表明，失速特性對(duì)參數(shù)S的敏感性高于H，即失速條在二維剖面安裝位置的選擇對(duì)觸發(fā)翼型氣流分離失速更加重要。需要注意的是，失速條在三維機(jī)翼上的展向長(zhǎng)度通常不大，且受到相鄰的無失速條翼段壓力分布傳遞的影響，因此安裝到機(jī)翼后失速特性的變化量會(huì)顯著低于二維情況。

圖4 失速條參數(shù)S對(duì)升力特性的影響Fig．4 The impact of stall strip parameter S on the lift characteristics

圖5 失速條參數(shù)H對(duì)升力特性的影響Fig．5 The impact of stall strip parameter H on the lift characteristics

對(duì)加裝失速條后翼型的流場(chǎng)特性進(jìn)行了研究，結(jié)果如圖6所示?？梢钥闯?，由于失速條的尖銳外形使氣流強(qiáng)制分離，在其后方形成了分離氣泡。在小迎角時(shí)氣泡局限在失速條后方很小的區(qū)域內(nèi)，隨著迎角增加氣泡在長(zhǎng)度和高度方向逐漸增大，最終氣泡破裂并產(chǎn)生較大的分離區(qū)。氣泡的增長(zhǎng)過程類似于薄翼失速，但氣泡能夠發(fā)展的長(zhǎng)度相對(duì)較小，通常不超過剖面15%弦長(zhǎng)即出現(xiàn)升力失速。不同失速條位置均有類似的氣泡流動(dòng)結(jié)構(gòu)，但出現(xiàn)的迎角時(shí)機(jī)不同，S越大氣泡出現(xiàn)的迎角越遲。

根據(jù)計(jì)算結(jié)果確定了開展風(fēng)洞實(shí)驗(yàn)驗(yàn)證的主要參數(shù)為S，其取值范圍應(yīng)當(dāng)小于30mm，根據(jù)風(fēng)洞實(shí)驗(yàn)驗(yàn)證的結(jié)果進(jìn)行調(diào)整。H的變化沒有導(dǎo)致失速特性的顯著變化，可固定為30mm。另外三維失速條的展向長(zhǎng)度和布置位置也是主要因素，需要結(jié)合機(jī)翼流場(chǎng)特性進(jìn)行試湊來確定。

圖6 失速條后方分離氣泡的發(fā)展（H＝30mm，S＝0mm）Fig．6 The development of the separation bubble after the stall strip（H＝30mm，S＝0mm）

2 機(jī)翼加裝失速條實(shí)驗(yàn)研究

2．1 風(fēng)洞實(shí)驗(yàn)概況

飛機(jī)模型加裝失速條的風(fēng)洞實(shí)驗(yàn)驗(yàn)證在中國空氣動(dòng)力與發(fā)展中心的FL－12風(fēng)洞中完成，該風(fēng)洞為單回流閉口式風(fēng)洞，風(fēng)洞截面尺寸為4m×3m，空風(fēng)洞最大風(fēng)速為100m／s，湍流度小于0．12%。實(shí)驗(yàn)?zāi)Ｐ捅壤秊?．07，基于飛機(jī)平均氣動(dòng)弦長(zhǎng)的雷諾數(shù)為1×106。飛機(jī)模型材質(zhì)為金屬，失速條為木質(zhì)的縮比模型，粘接在機(jī)翼前緣指定位置，失速條安裝在飛機(jī)模型上的照片如圖7所示。模型機(jī)翼采用自由轉(zhuǎn)捩，考慮到機(jī)翼翼型為厚翼型不存在前緣氣泡，且前緣邊界層很薄接近無粘流動(dòng)狀態(tài)，因此模型機(jī)翼前緣附近的流態(tài)與真實(shí)飛機(jī)是接近的。

圖7 失速條風(fēng)洞試驗(yàn)?zāi)Ｐ虵ig．7 Test model of the stall strip

2．2 失速條在機(jī)翼上的平面布置

為了給失速條布局設(shè)計(jì)提供參考，采用絲線法進(jìn)行了飛機(jī)襟翼35°構(gòu)型上翼面流場(chǎng)測(cè)試，結(jié)果如圖8所示。從流譜照片可見，飛機(jī)失速發(fā)生和發(fā)展出現(xiàn)在外側(cè)發(fā)動(dòng)機(jī)以內(nèi)大約為40%半展長(zhǎng)處，而外翼后緣雖然有局部氣流分離但隨迎角增加發(fā)展緩慢。結(jié)合飛機(jī)機(jī)翼設(shè)計(jì)的一般經(jīng)驗(yàn)，在40%半展長(zhǎng)處開始失速較為理想，但著陸襟翼構(gòu)型由于內(nèi)側(cè)機(jī)翼剖面升力大，若左右機(jī)翼失速不對(duì)稱則產(chǎn)生的滾轉(zhuǎn)力矩仍然較高，因此在40%半展長(zhǎng)以內(nèi)對(duì)稱地加裝失速條使機(jī)翼在一定迎角強(qiáng)制對(duì)稱失速，有可能改善飛機(jī)失速滾轉(zhuǎn)問題。

根據(jù)上述分析設(shè)計(jì)了4種失速條布局方案如圖9所示，實(shí)際尺寸失速條長(zhǎng)度有1和2m這2種方案，高度均為30mm，外側(cè)失速條中心點(diǎn)大致位于40%半展長(zhǎng)，內(nèi)側(cè)失速條中心點(diǎn)在15%半展長(zhǎng)。風(fēng)洞實(shí)驗(yàn)時(shí)失速條同樣采用了0．07的縮尺比。

圖8 著陸構(gòu)型機(jī)翼上翼面流場(chǎng)特性Fig．8 Flow characteristics of the upper wing at landing configuration

圖9 失速條平面布置方案Fig．9 Planform arrangement plans of stall strips

2．3 失速條布局方案1氣動(dòng)特性

以失速條布局方案1進(jìn)行了變參數(shù)S的氣動(dòng)特性影響實(shí)驗(yàn)，結(jié)果如圖10～12所示。從升力特性可見，隨著S降低，最大升力系數(shù)逐漸降低，失速迎角提前，S＝15mm和S＝0mm狀態(tài)的升力線在失速前出現(xiàn)了非線性區(qū)，失速形態(tài)變得緩和，S＝0mm狀態(tài)的升力特性呈平頂形。S＝0mm的失速條對(duì)升力失速特性改善明顯，但最大升力系數(shù)降低了0．15，降低了飛機(jī)的低速性能。

圖10 失速條參數(shù)S對(duì)飛機(jī)升力特性的影響（布局1）Fig．10 The impact of stall strip parameter S on the lift characteristics of aircraft（arrangement 1）

圖11 失速條參數(shù)S對(duì)飛機(jī)滾轉(zhuǎn)力矩特性的影響（布局1）Fig．11 The impact of stall strip parameter S on the roll moment of aircraft（arrangement 1）

圖12 失速條參數(shù)S對(duì)飛機(jī)阻力特性的影響（布局1）Fig．12 The impact of stall strip parameter S on the drag of aircraft（arrangement 1）

無側(cè)滑的滾轉(zhuǎn)力矩曲線表明，S＝30mm時(shí)Cl仍較大，S＝15mm時(shí)Cl明顯降低，S＝0mm時(shí)Cl隨迎角增大基本無波動(dòng)，表明飛機(jī)左右機(jī)翼失速對(duì)稱性顯著改善。原始構(gòu)型的實(shí)驗(yàn)?zāi)Ｐ陀捎谧笥覚C(jī)翼模型的微小偏差或流場(chǎng)不對(duì)稱性等因素出現(xiàn)不對(duì)稱滾轉(zhuǎn)力矩，真實(shí)飛機(jī)出現(xiàn)不對(duì)稱滾轉(zhuǎn)力矩的原因雖然不完全相同，但S＝0mm的失速條通過產(chǎn)生左右對(duì)稱的、顯著的局部氣流分離足以遮蔽外界擾動(dòng)的影響，將機(jī)翼的氣流分離形態(tài)固化，從而消除了隨機(jī)的不對(duì)稱滾轉(zhuǎn)力矩。隨著S量值增大，上述效應(yīng)減弱，不對(duì)稱滾轉(zhuǎn)力矩峰值也逐漸增大，不能起到預(yù)期的作用。

從阻力曲線可見，小迎角時(shí)各失速條沒有產(chǎn)生明顯的阻力增加，在接近失速點(diǎn)后阻力增大，S越小阻力增量越大。結(jié)合二維流場(chǎng)特性分析，這一現(xiàn)象是由于一方面失速條尺寸很小基本不增加飛機(jī)的摩擦阻力，另一方面小迎角時(shí)分離氣泡局限在較小的范圍內(nèi)，不會(huì)顯著改變機(jī)翼的壓力分布，產(chǎn)生的壓差阻力也較小，直到氣流分離點(diǎn)后才會(huì)出現(xiàn)阻力的激增。

2．4 失速條布局方案2、3、4氣動(dòng)特性

位于2個(gè)發(fā)動(dòng)機(jī)之間的失速條布局方案1、2、3的氣動(dòng)特性對(duì)比如圖13和14所示。從升力曲線對(duì)比可知，失速條縮短長(zhǎng)度后對(duì)失速特性的影響變?nèi)?，靠外的布?的升力失速優(yōu)于布局2，結(jié)合圖8分析，其原因是此處更加接近原始構(gòu)型氣流分離的嚴(yán)重區(qū)。從無側(cè)滑滾轉(zhuǎn)力矩曲線可見，布局3產(chǎn)生的滾轉(zhuǎn)力矩也小于布局2。

圖13 失速條不同布局方案的升力特性Fig．13 The lift characteristics of different stall strip arrangement plans

圖14 失速條不同布局方案的滾轉(zhuǎn)力矩特性Fig．14 The roll moment of different stall strip arrangement plans

布局方案4的升力曲線和滾轉(zhuǎn)力矩曲線如圖15和16所示，可見在S＝－15～30mm之間變化時(shí)僅使最大升力系數(shù)略有降低，但失速形態(tài)基本不變，滾轉(zhuǎn)力矩的偏離沒有明顯變化，沒有起到改善失速特性的作用，同時(shí)也說明失速條在機(jī)翼平面內(nèi)的布置位置是改善飛機(jī)失速特性的關(guān)鍵參數(shù)。

圖15 失速條布局方案4的升力特性Fig．15 The lift characteristics of stall strip arrangement plan 4

圖16 失速條布局方案4的滾轉(zhuǎn)力矩特性Fig．16 The roll moment of stall strip arrangement plan 4

3 結(jié) 論

采用CFD計(jì)算方法對(duì)失速條二元方案進(jìn)行了選型，獲取了能夠使失速迎角提前的失速條高度和剖面布置位置范圍，并通過流動(dòng)顯示發(fā)現(xiàn)了失速條后方的氣泡流動(dòng)結(jié)構(gòu)。通過風(fēng)洞實(shí)驗(yàn)對(duì)4種失速條布局方案進(jìn)行了驗(yàn)證，結(jié)果表明，失速條安裝在氣流分離起始位置并選擇合適的剖面參數(shù)時(shí)能夠有效改善失速特性，使升力曲線由急劇失速變?yōu)槠巾斝褪伲瑫r(shí)消除或削弱無側(cè)滑條件大迎角的滾轉(zhuǎn)力矩。在失速迎角之前，加裝失速條后的飛機(jī)構(gòu)型阻力沒有明顯增加。

［1］中國民用航空局．CCAR－25－R4運(yùn)輸類飛機(jī)適航標(biāo)準(zhǔn)［S］．2011．

［2］Torenbeek E．Synthesis of subsonic aircraft design［M］．Delft University Press，1982：227－262．

［3］Hoerner S F，Borst H V．Fluid－dynamic lift［M］．2nd ed．Published by the auther，1985．

［4］Mc Veigb M A，Kisielowski E．A design summary of stall characteristics of straight wing aircraft［R］．NASA CR－1646，1971：1－20．

［5］Obert E．Aerodynamic design of transport aircraft［M］．IOS Press，2009．

［6］Feistel T W，Anderson S B，Kroeger R A．Alleviation of spinentry tendencies through localization of wing－flow separation［J］．Journal of Aircraft，1981，18（2）：69－75．

［7］Newsom W A，Satran D R，Johnson J L．Effects of wing－leading－edge modifications on a full－scale，low－wing general aviation airplane［R］．NASA－TP－2011，1982．

［8］Chambers J R，Dicarlo D J，Johnson J L，et al．Exploratory study of the effects of wing－leading－edge modifications on the stall／spin behavior of a light general aviation airplane［R］．NASA－TP－1589，1979．

［9］Jacobs E N．Characteristics of two sharp－nosed airfoils having reduced spinning tendencies［R］．NACA－TN－416，1932．

［10］Weick F E，Scudder N F．The effect on lift，drag，and spinning characteristics of sharp leading edges on airplane wings［R］．NACA－TN－447，1933．

［11］Owens D B，Capone F J，Hall R M，et al．Free－to－Roll analysis of abrupt wing stall on military aircraft at transonic speeds［R］．AIAA－2003－0750，2003．

［12］Capone F J，Hall R M，Owens D B，et al．Recommended experimental procedures for evaluation of abrupt wing stall characteristics［R］．AIAA－2003－0922，2003．

［13］Parikh P，Chung J．A computational study of the Abrupt Wing Stall（AWS）characteristics for various fighter jets：part I，F(xiàn)／A－18E and F－16C［R］．AIAA－2003－0746，2003．

［14］Hall R M．Accomplishments of the Abrupt Wing Stall（AWS）program and future research requirements［R］．AIAA－2003－0927，2003．

［15］Woodson S H，Green B E，Chung J J，et al．Understanding abrupt wing stall with CFD［R］．AIAA－2003－0592，2003．

［16］Liu Y，Zhao X X，Ouyang S X，et al．A method to predict the magnitude of roll during stall for transport aircraft and its application［C］．The 55th Israel Annual Conference on Aerospace Sciences，2015．

［17］ANSYS．ANSYS FLUENT theory guide［M］．ANSYS Inc，Canonsburg，PA，2011．

Research on aerodynam ic characteristics of transport aircraft w ith stall strips

Liu Yi*，Zhao Xiaoxia，Ouyang Shaoxiu，Yuan Zhimin
(Research and Development Center，Avcation Industry Corporation of China Aircraft Co．Ltd．，Xi’an 710089，China)

In order to alleviate the violent rollmotion during stall of a transport aircraftwith landing flap configuration，the geometric parameters of stall strips are optimized and selected by numerical simulations and wind tunnel tests，and the aerodynamic force and flow field characteristics are studied．The height H and the install distance S from the leading edge are selected as design parameters for stall strips，and are evaluated by solving Reynolds Averaged Navior-Stokes(RANS)equations for the airfoil section of the landing flap configuration．The calculation indicates that smaller S value(installed closer to the upper surface)promotes earlier stall，while the larger H has similar butweaker effect．Separation bubble emerges after the stall stripswhen the angle of attack(AOA)of the airfoil becomes large，which grows larger and longer with increasing AOA．The bubble bursts eventually and causes the airfoil to stall earlier，leading to the rounded shape of the lift curve．The effect of the stall strips installed on thewing is studied by scaled model in wind tunnel tests，which shows that its spanwise length and arrangement have significant impact on the performance besides the cross section geometry．Four planform arrangements of stall strips are advanced and evaluated．Keeping the S parameter equal to 0，the stall trips installed at40%half span with spanwise length of 2m change the abrupt stall of the lift curve to a flat roof type one，and eliminate the asymmetric rollmoment after stall．The spanwise length of the stall strips is halved to form 2 new arrangements，which also ameliorate the stall of the lift curve and rollmoment to some extent．The stall trips installed at 15%half span have no obvious effects on the stall characteristics．The suggested explanation is that the flow separation starts at about40%half span of thewing at landing configuration，where the stall strips have the best performance．

V211．7

A

（編輯：李金勇）

1672-9897(2016)05-0036-06

10．11729/syltlx20160012

2016-01-12;

2016-04-20

*通信作者E-mail:evanliuyi@hotmail．com

Liu Y，Zhao X X，Ouyang S X，etal．Research on aerodynam ic charac teristics of transportairc raftw ith stallstrips．Journalof Experiments in Fluid Mechanics，2016，30(5):36-41．劉毅，趙曉霞，歐陽紹修，等．某運(yùn)輸機(jī)加裝失速條氣動(dòng)特性研究．實(shí)驗(yàn)流體力學(xué)，2016，30(5):36-41．

劉 毅（1982－），男，四川資陽人，碩士，高級(jí)工程師。研究方向：飛機(jī)氣動(dòng)力設(shè)計(jì)。通信地址：陜西省漢中市五一路陜飛大廈（723000）。E－mail：evanliuyi＠hotmail．com