基于雙橢圓弧型中弧線的系列翼型設(shè)計方法

張 森，席德科，劉治斌，陳寶峰

?

基于雙橢圓弧型中弧線的系列翼型設(shè)計方法

張 森1，席德科1※，劉治斌2，陳寶峰2

（1. 西北工業(yè)大學(xué)航空學(xué)院，西安 710072； 2. 山西安瑞風(fēng)機電氣股份有限公司，運城 044402）

為了能夠方便快捷的設(shè)計和修改翼型，采用兩段橢圓弧來構(gòu)造翼型的中弧線，并推導(dǎo)了描述中弧線的方程式。用該方法構(gòu)造的中弧線光滑連續(xù)，且不存在拐點。選用現(xiàn)有翼型的厚度分布，與中弧線分布函數(shù)進行疊加，并引入厚度比例因子來實現(xiàn)對厚度的調(diào)整，最終得到了一種基于雙橢圓弧型中弧線的翼型設(shè)計方法，稱之為DEA（double ellipse arcs）翼型。選用Clark-Y翼型作為基礎(chǔ)翼型，設(shè)計了多款DEA翼型，并利用X-foil軟件對翼型氣動性能進行求解，分別研究了最大相對彎度、最大彎度相對位置、最大相對厚度以及翼型中弧線的形狀因子對翼型氣動性能的影響。研究表明：增加最大相對彎度，可以提高翼型的升力系數(shù)，同時使翼型的升阻特性得到一定的改善；最大彎度位置前移，可以提高翼型在小攻角下的升力系數(shù)，同時增加翼型高效升阻比的攻角范圍；增加最大相對厚度可以提高翼型的最大升力系數(shù)，以及增大失速攻角，同時，高效升阻比的攻角范圍也隨著翼型最大相對厚度的增大而增加；中弧線前、后緣形狀因子對翼型氣動性能的影響相對較小。

設(shè)計；翼型；數(shù)值分析；中弧線；橢圓弧

0 引 言

翼型作為航空技術(shù)發(fā)展的產(chǎn)物，被廣泛的應(yīng)用于流體機械產(chǎn)品的設(shè)計中，如壓縮機、軸流風(fēng)機、軸流泵、風(fēng)力機等，翼型的氣動特性是決定流體機械性能優(yōu)劣的一個關(guān)鍵因素[1-5]。

當(dāng)選擇已有的翼型尚不能滿足工程要求時，需要重新設(shè)計或者對原有翼型進行修型。目前此項工作主要依靠計算空氣動力學(xué)的方法設(shè)計，以減少甚至取代翼型風(fēng)洞試驗，具體實施方法有直接法和逆設(shè)計[6-8]。直接法基本過程：首先確定目標，然后人工修改翼型型線或數(shù)據(jù)，計算氣動特性并與設(shè)計要求比較分析，重復(fù)進行修正，直到滿足要求，該方法要求設(shè)計者有較深厚的專業(yè)知識和豐富的設(shè)計經(jīng)驗。逆設(shè)計基本過程：給定壓力分布目標函數(shù)和約束條件，通過智能優(yōu)化方法自動修改翼型外形，經(jīng)過多次迭代，達到收斂條件為止，該方法可以進行多點/多目標的設(shè)計，但計算量大，特別是設(shè)計變量較多時尤為突出[9-15]。

翼型的型線決定了繞翼型流場特性參數(shù)的分布，即翼型的氣動性能，是翼型分析的基礎(chǔ)和關(guān)鍵。翼型的型線表達方法主要有：外形參數(shù)化方法[16]、形函數(shù)擾動法[17]、解析函數(shù)法[18-19]。Ray T等[20]采用翼型的特征參數(shù)對翼型外形參數(shù)化，設(shè)計變量幾何意義明確，但未能給出解析表達式。形函數(shù)擾動法是由原始翼型和擾動形函數(shù)的線性疊加實現(xiàn)[21-22]。解析函數(shù)法就是用一個解析函數(shù)直接表示翼型形狀，例如早期用多項式表達的NACA的4位數(shù)、5位數(shù)系列翼型，此外也有研究用級數(shù)表達翼型的方法[23]。

Mark[24]于1989年開發(fā)了X-foil翼型分析與設(shè)計系統(tǒng)，該軟件對于黏性流體采用了面元法和邊界層理論，由于其計算簡單方便，適合于低速翼型的快速分析和設(shè)計，因此得到了廣泛應(yīng)用[25-29]。Ashok G等[30]對比了NLF(1)-0416和NLF(1)-0215F翼型的X-foil計算結(jié)果與試驗數(shù)據(jù)，指出X-foil程序可以用于優(yōu)化設(shè)計中計算不同外型的翼型氣動性能。鄧磊等[31]研究并發(fā)展了一套進行高升阻比自然層流翼型多設(shè)計點/多設(shè)計目標的優(yōu)化設(shè)計方法。為減少設(shè)計中的計算量，使用了X-foil程序進行流場計算，進一步驗證了該軟件的可靠性。

本文采用兩段橢圓弧來構(gòu)造翼型的中弧線，并推導(dǎo)了中弧線的方程式，最終得到了一種基于雙橢圓弧型中弧線的翼型設(shè)計方法，利用X-foil軟件分析了翼型幾何特征參數(shù)對該系列翼型氣動性能的影響規(guī)律，并給出翼型設(shè)計指導(dǎo)準則。

1 翼型幾何特征參數(shù)

圖1中給出了翼型的幾何特征參數(shù)。

注：c表示翼型的弦長，mm，將c設(shè)定為1；t表示翼型的最大相對厚度；f表示翼型的最大相對彎度；xt表示翼型最大厚度相對弦長的位置；xf表示翼型最大彎度相對弦長的位置。

2 雙橢圓弧型中弧線函數(shù)的構(gòu)造

2.1 雙橢圓弧型中弧線函數(shù)的求解

采用頂點重合的2個橢圓上截得的橢圓弧來構(gòu)造翼型中弧線，具體方法如圖2所示，為了方便描述，將翼型的弦長定義為1。

圖2 雙橢圓弧型中弧線示意圖

圖2中，點(0,0)與(1,0)分別為翼型的前緣點和后緣點，點(x,)為橢圓的頂點。中弧線型線方程y如式（1）所示。

式中1為翼型中弧線前緣形狀因子，?1<1≤0；2為翼型中弧線后緣形狀因子，?1<2≤0。

2.2 中弧線形狀因子k1和k2對中弧線形狀的影響

當(dāng)翼型中弧線的最大相對彎度和最大彎度相對位置x一定時，通過調(diào)整方程（1）中形狀因子1和2的大小可以局部調(diào)整中弧線的形狀。

為了研究形狀因子1和2對中弧線形狀的影響，在圖3中給出了最大相對彎度=3.43%，最大彎度相對位置x=40%時，中弧線前緣形狀因子1=?0.99，中弧線后緣形狀因子2分別為0、?0.2、?0.5、?0.7、?0.99，和中弧線后緣形狀因子2=?0.99，中弧線前緣形狀因子1分別為0、?0.2、?0.5、?0.7、?0.99時的型線圖。

圖3 不同k1和k2時的中弧線型線

由圖3可知，當(dāng)中弧線前緣形狀因子一定而中弧線后緣形狀因子變化時，前緣點至最大彎度處翼型中弧線型線保持不變，翼型后緣氣流出口角隨著中弧線后緣形狀因子的增大而增大，當(dāng)中弧線后緣形狀因子2=0時，達到最大值。同樣，當(dāng)中弧線后緣形狀因子一定而中弧線前緣形狀因子變化時，最大彎度至后緣點翼型中弧線型線保持不變，翼型前緣氣流出口角隨著中弧線前緣形狀因子的增大而增大，當(dāng)中弧線前緣形狀因子1=0時，達到最大值。結(jié)果表明，當(dāng)中弧線的最大相對彎度和最大彎度相對位置一定時，隨著形狀因子1與2增大，翼型中弧線型線變得飽滿，對應(yīng)的氣流角也逐漸增大，且形狀因子1與2對型線的影響以最大彎度位置為分界點，相互獨立。

3 系列翼型設(shè)計方法

翼型的型面函數(shù)可以表示為中弧線分布函數(shù)與厚度分布函數(shù)的疊加，如式（2）所示。

式中y、y分別表示翼型的上表面型面函數(shù)和下表面型面函數(shù)；y表示翼型的厚度分布函數(shù)。

式（2）中的中弧線分布函數(shù)y使用本文推導(dǎo)的式（1），厚度分布函數(shù)則選用現(xiàn)有翼型的厚度分布。為了能夠調(diào)節(jié)翼型最大相對厚度的大小，引入了厚度比例因子，其表達式為

式中des表示設(shè)計翼型的最大相對厚度；ori表示原始翼型的最大相對厚度。

因此，可根據(jù)實際需要，通過調(diào)整厚度比例因子的大小得到理想的最大相對厚度。在式（2）的厚度項y前乘以系數(shù)，得到

式（4）給出了最終的翼型型面函數(shù)表達式。本文將該方法構(gòu)造的翼型稱之為DEA（double ellipse arcs）翼型。

4 結(jié)果與分析

選用Clark-Y翼型作為基礎(chǔ)翼型，運用上述DEA翼型設(shè)計方法進行翼型設(shè)計。Clark-Y翼型的幾何特征參數(shù)如下：最大相對厚度=11.71%，最大厚度相對位置x=28%，最大相對彎度=3.43%，最大彎度相對位置x=42%。在與Clark-Y翼型具有相同幾何特征參數(shù)的條件下，取中弧線形狀因子1=?0.8，2=?0.99，設(shè)計了一款DEA翼型，如圖4所示。

圖4 Clark-Y翼型與DEA翼型的型線對比

Clark-Y翼型與DEA翼型的型線對比結(jié)果顯示，兩者的型線有較好的重合度，采用DEA翼型設(shè)計方法能夠便捷的對現(xiàn)有翼型進行參數(shù)化。

為了研究幾何特征參數(shù)對DEA翼型氣動性能的影響，本文利用X-foil軟件對翼型氣動性能進行求解，取雷諾數(shù)=1.0×105。

4.1 最大相對彎度f對DEA翼型氣動性能的影響

在最大相對厚度=11.71%、最大厚度相對位置x=28%和最大彎度相對位置x=42%的情況下，取中弧線形狀因子1=?0.8，2=?0.99，分別設(shè)計了最大相對彎度為3.42%、4%、5%、6%時的翼型，氣動性能計算結(jié)果如圖5所示。

由圖5可看出，在最大相對厚度、最大厚度相對位置和最大彎度相對位置相同的情況下，基礎(chǔ)翼型Clark-Y與最大相對彎度=3.43%時的DEA翼型的性能曲線有較高的重合度，這里進一步驗證了用DEA翼型參數(shù)化設(shè)計方法設(shè)計翼型的可靠性。當(dāng)只改變最大相對彎度時，翼型的升力系數(shù)曲線在小攻角范圍內(nèi)基本保持平行，在同一攻角下，隨著最大相對彎度的增大，升力系數(shù)也隨之增大。此外，隨著最大相對彎度的增大，翼型的升阻比曲線也發(fā)生了變化，當(dāng)最大相對彎度=6%時，在0°、6°和12°攻角處，翼型的升阻比相對于Clark-Y翼型分別提高了6.88%、1.59%和24.49%。研究表明，在翼型厚度、最大厚度位置和最大彎度位置不變的情況下，增大相對彎度，可以提高DEA翼型的升力系數(shù)，同時使翼型的升阻特性得到一定的改善。

注：表示翼型攻角，(o)；Cl表示翼型的升力系數(shù)；Cd表示翼型的阻力系數(shù)；Cl/Cd表示翼型的升阻比。

4.2 最大彎度相對位置xf對DEA翼型氣動性能的影響

為了研究最大彎度相對位置x對DEA翼型氣動特性的影響，在最大相對厚度=11.71%、最大厚度相對位置x=28%和最大相對彎度=3.43%的條件下，取中弧線形狀因子1=?0.8，2=?0.99，分別設(shè)計了最大彎度相對位置x為25%、42%、50%、60%時的翼型，并對其氣動性能進行了計算，計算結(jié)果如圖6所示。

圖6 不同xf下的DEA翼型氣動性能曲線

4.3 最大相對厚度t對DEA翼型氣動性能的影響

由薄翼理論可知，對于理想不可壓縮流體的翼型繞流，如果氣流繞翼型的迎角、翼型厚度、翼型彎度都很小，則繞流場是一個小擾動的勢流場。這時，翼面上的邊界條件和壓強系數(shù)可以線性化為厚度、彎度、迎角三者影響的疊加。因此翼型的厚度對氣動性能也起到至關(guān)重要的影響。圖7中給出了最大彎度相對位置x=25%、最大厚度相對位置x=28%、最大相對彎度=3.43%，最大相對厚度分別為6%、8%、10%和11.71%時的DEA翼型氣動特性曲線，中弧線形狀因子取1=?0.8，2=?0.99。

圖7 不同t下的DEA翼型氣動性能曲線

4.4 中弧線形狀因子k1與k2對DEA翼型氣動性能的影響

為了研究中弧線形狀因子1和2對DEA翼型氣動特性的影響，分別設(shè)計了最大相對厚度=11.71%、最大彎度相對位置x=42%、最大厚度相對位置x=28%、最大相對彎度=3.43%時，中弧線形狀因子1=?0.2、?0.5、?0.7、?0.99，2=?0.99和1=?0.99，2=?0.2、?0.5、?0.7、?0.99的翼型，性能計算結(jié)果如圖8和圖9所示。

圖8 不同k1下的DEA翼型氣動性能曲線

圖9 不同k2下的DEA翼型氣動性能曲線

5 結(jié) 論

1）采用兩段橢圓弧來構(gòu)造翼型中弧線，并推導(dǎo)了中弧線方程式，該方程通過改變最大彎度和最大彎度相對位置對中弧線形狀進行控制，并由翼型中弧線形狀因子實現(xiàn)局部微調(diào)。用該方法構(gòu)造的中弧線光滑連續(xù)，且不存在拐點。

2）選用現(xiàn)有翼型的厚度分布，與中弧線分布函數(shù)進行疊加，并引入厚度比例因子來實現(xiàn)對厚度的控制，最終得到了一種基于雙橢圓弧中弧線的翼型設(shè)計方法，本文將該方法構(gòu)造的翼型稱之為DEA（double ellipse arcs）翼型。

3）選用Clark-Y翼型作為基礎(chǔ)翼型，設(shè)計了多款DEA翼型，并利用X-foil軟件對翼型氣動性能進行求解，分別研究了最大相對彎度、最大彎度相對位置、最大相對厚度，以及翼型中弧線的形狀因子對翼型氣動性能的影響。研究表明：增加最大相對彎度，可以提高翼型的升力系數(shù)，同時使翼型的升阻特性得到一定的改善；最大彎度位置前移，可以提高翼型在小攻角下的升力系數(shù)，同時增加翼型高效升阻比的攻角區(qū)間，但最優(yōu)升阻比會逐漸減小；增加最大相對厚度可以提高翼型的最大升力系數(shù)，以及增大失速攻角，同時，高效升阻比的攻角區(qū)間也隨著翼型厚度的增大而增加，但最優(yōu)升阻比會逐漸減小；中弧線形狀因子對翼型氣動性能的影響較小。

依據(jù)上述參數(shù)的調(diào)整變化原則，可設(shè)計出新的翼型或者對原翼型進行修型以滿足目標任務(wù)的需要。

[1] Mojtaba Tahani, Ghazale Kavari, Mehran Masdari, et al. Aerodynamic design of horizontal axis wind turbine with innovative local linearization of chord and twist distributions[J]. Energy, 2017, 131: 78－91.

[2] 項松，王吉，張利國，等. 一種高效率螺旋槳設(shè)計方法[J]. 航空動力學(xué)報，2015，30(1)：136－141.

Xiang Song, Wang Ji, Zhang Liguo, et al. A design method for high efficiency propeller[J]. Journal of Aerospace Power, 2015, 30(1): 136－141. (in Chinese with English abstract)

[3] 沙毅，侯麗艷. 葉片厚度對軸流泵性能影響及內(nèi)部流場分析[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(18)：75－81.

Sha Yi, Hou Liyan. Effects of blade thickness on performance of axial flow pump and analysis of internal flow field[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(18): 75－81. (in Chinese with English abstract)

[4] 汪泉，陳進，王君，等. 基于連續(xù)攻角的風(fēng)力機翼型整體氣動性能提高的優(yōu)化設(shè)計[J]. 機械工程學(xué)報，2017，53(13)：143－149.

Wang Quan, Chen Jin, Wang Jun, et al. Wind turbine airfoil optimal design with high whole aerodynamic performance considering continuous angle of attack[J]. Journal of Mechanical Engineering, 2017, 53(13): 143—149. (in Chinese with English abstract)

[5] 徐浩然，楊華，劉超. 尾緣加厚的DU系列翼型氣動性能數(shù)值分析[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(17)：101－108.

Xu Haoran, Yang Hua, Liu Chao. Numerical value analysis on aerodynamic performance of DU series airfoils with thickened trailing edge[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(17): 101—108. (in Chinese with English abstract)

[6] Srilatha K R, Dwarakanath G S. Design of a natural laminar flow airfoil for a light aircraft[J]. Journal of Aircraft, 1990, 27(11): 966－968.

[7] 王驥月，叢茜，梁寧，等. 基于海鷗翼型的小型風(fēng)力機葉片仿生設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(10)：72－77.

Wang Jiyue, Cong Qian, Liang Ning, et al. Bionic design and test of small-sized wind turbine blade based on seagull airfoil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(10): 72—77. (in Chinese with English abstract)

[8] Green B E. An Approach to the Constrained Design of Natural Laminar flow Airfoils[M]. NASA Langley Technical Report Server, 1997.

[9] 嚴敬，劉小兵，周緒成，等. 基于奇點分布法的軸流泵葉片翼型設(shè)計與計算[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(7)：100－105.

Yan Jing, Liu Xiaobing, Zhou Xucheng, et al. Design and calculation of airfoil profile of blade in axial flow pump based on singularity approach[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(7): 100－105. (in Chinese with English abstract)

[10] 龍騰，劉莉，孟令濤，等. 基于代理模型的翼型多目標集成優(yōu)化設(shè)計[J]. 系統(tǒng)仿真學(xué)報，2010，22(7)：1604－1608.

Long Teng, Liu Li, Meng Lingtao, et al. Multi-objective integrated optimization of airfoil based on surrogate model[J]. Journal of System Simulation, 2010, 22(7): 1604－1608. (in Chinese with English abstract)

[11] Lee J D, Jameson A. Natural-laminar flow airfoil and wing design by adjoint method and automatic transition prediction [C]// 47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. Orlando, Florida, 2009.

[12] Youngren H. Multi-Point design and optimization of an natural laminar flow airfoil for a mission adaptive compliant wing[C]// AIAA Aerospace Sciences Meeting and Exhibit. 2008:485601-485605(5).

[13] Driver J, Zingg D W. Optimized natural-laminar flow airfoils, AIAA-2006-247[R]. Toronto: American Institute of Aeronautics and Astronautics, 2006.

[14] Ashok G, Michael S S. Low speed NLF airfoils: Case study in inverse airfoil design[J]. Journal of Aircraft, 2001, 38(2): 57－63.

[15] 朱國俊，馮建軍，郭鵬程，等. 基于徑向基神經(jīng)網(wǎng)絡(luò)-遺傳算法的海流能水輪機葉片翼型優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(8)：65－73.

Zhu Guojun, Feng Jianjun, Guo Pengcheng, et al. Optimization of hydrofoil for marine current turbine based on radial basis function neural network and genetic algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(8): 65－73. (in Chinese with English abstract)

[16] Liang Xiao, Meng Guanglei, Tong Shengxi, et al. Rapid design and optimization of airfoil based on improved genetic algorithm[J]. Acta Aerodynamica Sinica, 2016, 31(6): 803－812.

[17] 王迅，蔡晉生，屈崑，等. 基于改進CST參數(shù)化方法和轉(zhuǎn)捩模型的翼型優(yōu)化設(shè)計[J]. 航空學(xué)報，2015，36(2)：449－461.

Wang Xun, Cai Jinsheng, Qu Kun, et al. Airfoil optimization based on improved CST parametric method and transition model[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(2): 449－461. (in Chinese with English abstract)

[18] 姜海波，趙云鵬. 基于中弧線—厚度函數(shù)的翼型形狀解析構(gòu)造法[J]. 圖學(xué)學(xué)報，2013，34(1)：50－54.

Jiang Haibo, Zhao Yunpeng. Analytic expression method of airfoil profile shape based on a mean camber and thickness function[J]. Journal of Graphics, 2013, 34(1): 50－54. (in Chinese with English abstract)

[19] 姜海波，程忠慶，趙云鵬. 翼型型線的參數(shù)表達式探討[J]. 機械設(shè)計與制造，2015，3：69－72.

Jiang Haibo, Cheng Zhongqing, Zhao Yunpeng. Study on parameter expressions of airfoil contour[J]. Machinery Design and Manufacture, 2015, 3: 69－72. (in Chinese with English abstract)

[20] Ray T, Tsai H M. Swarm algorithm for single and multiobjective airfoil design optimization[J]. AIAA Journal, 2004, 42(2): 366－373.

[21] David W Z, Timothy M L, Laslo D, et al. Improvements to a newton-krylov adjoint algorithm for aerodynamic optimization, AIAA-2005-4857[R]. Toronto: American Institute of Aeronautics and Astronautics, 2005.

[22] 劉周，朱自強，付鴻雁，等. 高升阻比翼型的設(shè)計[J]. 空氣動力學(xué)學(xué)報，2004，22(4)：410－415.

Liu Zhou，Zhu Ziqiang，F(xiàn)u Hongyan，et al. Design of airfoil with high ratio of lift over drag[J]. Acta Aerodynamica Sinica, 2004, 22(4): 410－415. (in Chinese with English abstract)

[23] 陳進，張石強，Eecen P J，等. 風(fēng)力機翼型參數(shù)化表達及收斂特性[J]. 機械工程學(xué)報，2010，46(10)：132－138.

Chen Jin, Zhang Shiqiang, Eecen P J, et al. Parametric representation and convergence of wind turbine airfoils[J]. Journal of Mechanical Engineering, 2010, 46(10): 132－138. (in Chinese with English abstract)

[24] Mark D. XFOIL: An analysis and design system for low reynolds number airfoils[J]. Lecture Notes in Engineering, 1989, 54: 1－12.

[25] Bangga G, Hutomo G, Wiranegara R, et al. Numerical study on a single bladed vertical axis wind turbine under dynamic stall[J]. Journal of Mechanical Science and Technology, 2017, 31(1): 261－267.

[26] Ananda G K, Sukumar P P, Selig M S. Measured aerodynamic characteristics of wings at low Reynolds numbers[J]. Aerospace Science & Technology, 2015, 42: 392－406.

[27] Morgado J, Vizinho R, Silvestre M A R, et al. XFOIL vs CFD performance predictions for high lift low Reynolds number airfoils[J]. Aerospace Science and Technology, 2016, 52: 207－214.

[28] Ramanujam G, ?zdemir H, Hoeijmakers H W M. Improving airfoil drag prediction[C]//Wind Energy Symposium, AIAA Scitech. 2016.

[29] Paulo A S F S, Léo D S, Taygoara F O, et al. Analysis of cavitation for the optimized design of hydrokinetic turbines using BEM[J]. Applied Energy, 2017, 185: 1281－1291.

[30] Ashok G, Michael S S. Low-speed natural-laminar-flow airfoils: case study in inverse airfoil design[J]. Journal of Aircraft, 2001, 38(1): 57－63.

[31] 鄧磊，喬志德，楊旭東，等. 高升阻比自然層流翼型多點/多目標優(yōu)化設(shè)計[J]. 空氣動力學(xué)學(xué)報，2011，29(3)：330－335.

Deng Lei, Qiao Zhide, Yang Xudong, et al. Multi-point/ objective optimization design of high lift-to-drag ratio for NLF airfoils[J]. Acta Aerodynamica Sinica, 2011, 29(3): 330－335. (in Chinese with English abstract)

Design method for series airfoil based on mean camber line consisting of double ellipse arcs

Zhang Sen, Xi Deke, Liu Zhibin, Chen Baofeng

Institute of Aeronautics, Northwestern Polytechnical University, Xi’an, ChinaShanxi Anrui Fan Electric Co., Ltd, Yuncheng, China

Airfoil, as a product of aviation technology, has been widely used in the design of fluid machinery products. The aerodynamic characteristics of airfoils are a key factor in determining the performance of fluid machinery. When the existing airfoils are not able to meet the engineering requirements, it is necessary to redesign or trim the original airfoils. In this research, 2 ellipse arcs were used to form the mean camber line of the airfoil, and the corresponding equation was deduced. This equation controls the shape of the mean camber line by changing the maximum camber and the relative position of maximum camber, and adjusts the local shape by changing the 2 shape factors of the mean camber line. The mean camber line constructed by this method is smooth and continuous, and there is no knee point. Then the thickness distribution of the existing airfoil was superposed with the distribution function of mean camber line, and a thickness scale factor was introduced to adjust the thickness distribution. Ultimately, the design method for a series of airfoils based on mean camber line of double ellipse arcs is achieved, which is called DEA (double ellipse arcs) airfoil. The airfoil profile function constructed by this method has definite physical meaning, simple and reliable, and it is easy to realize serialization. In order to study the influence of airfoil characteristic parameters on aerodynamic performance of the DEA airfoil, the Clark-Y airfoil was taken as the basic airfoil, and a number of DEA airfoils were designed using the thickness distribution of the Clark-Y airfoil. Then the aerodynamic characteristics of the designed airfoils were solved by the X-foil software to study the influence of the maximum camber, the relative position of the maximum camber, the maximum thickness and the shape factors of the mean camber line on the DEA airfoil aerodynamic performance. There are 5 characteristic parameters in all that influence the shape of the DEA airfoil. We selected one of the 5 characteristic parameters as variable and fixed the other 4 characteristic parameters to design different DEA airfoils. And the aerodynamic characteristics were achieved at Reynolds number of 1.0×105. The calculation results of the 4 DEA airfoils with different values of maximum camber show that the increase of the maximum camber can improve the lift coefficient and ameliorate the characteristics of the lift-drag ratio. The calculation results of the 4 DEA airfoils with different values of relative position of the maximum camber show that as the relative position of the maximum camber moves forward, the lift coefficient under small angles of attack is improved, and the range of efficient lift-drag ratio gets broadened. The calculation results of the 4 DEA airfoils with different values of maximum thickness show that the increase of the maximum thickness can increase the maximum lift coefficient and the stall angle. At the same time, with the increase of thickness, the range of efficient lift-drag ratio also gets broadened. The calculation results of the DEA airfoils with different shape factors of the mean camber line also were achieved. At small attack angle, the change of the leading shape factors of the mean camber line has little influence on lift coefficient. With the decrease of the leading shape factors of the mean camber line, the interval of efficient lift-drag ratio has a tendency to move to high attack angle range. With the decrease of the trailing shape factors of the mean camber line, the lift coefficient and lift-drag ratio decrease gradually. Moreover, the interval of efficient lift-drag ratio also decreases and the decrease is mainly at the range of small attack angle. According to the adjustment principle of the above parameters, a new airfoil can be designed or modified to meet the needs of the target task.

design; airfoils; numerical analysis; mean camber line; double ellipse arcs

10.11975/j.issn.1002-6819.2018.02.006

TB126/TB21

A

1002-6819(2018)-02-0040-07

2017-07-07

2017-12-12

國家自然科學(xué)基金（11172243）；陜西省科技統(tǒng)籌創(chuàng)新工程計劃（2011KTCQ01-02）

張 森，河南新鄉(xiāng)人，博士生，研究方向為航空高科技軍轉(zhuǎn)民技術(shù)及流體機械設(shè)計。Email：sen96@mail.nwpu.edu.cn

席德科，教授，博士生導(dǎo)師，研究方向為風(fēng)洞設(shè)計及流體機械設(shè)計。Email：xideke@nwpu.edu.cn

張 森，席德科，劉治斌，陳寶峰. 基于雙橢圓弧型中弧線的系列翼型設(shè)計方法[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(2)：40－46. doi：10.11975/j.issn.1002-6819.2018.02.006 http://www.tcsae.org

Zhang Sen, Xi Deke, Liu Zhibin, Chen Baofeng. Design method for series airfoil based on mean camber line consisting of double ellipse arcs[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(2): 40－46. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.02.006 http://www.tcsae.org