張照煌，高 迪，孫 佳，楊方超，胡德鵬

（華北電力大學 能源動力與機械工程學院，北京 102206）

1 研究背景

風能作為一種清潔的可再生能源，已經(jīng)成為世界各國的研究重點[1]。葉片是風力機將風能轉(zhuǎn)化為機械能的核心部件，其氣動性能的優(yōu)劣直接決定風力機對風能的利用效率。由于水平軸風力機運行環(huán)境復(fù)雜，所以對葉片不同展向位置的要求有所不同，風力機葉片內(nèi)側(cè)運行攻角較大，流動易分離且存在明顯的旋轉(zhuǎn)效應(yīng)，要求該部位的翼型具有較大的橫截面積和良好的厚度分布以獲得較高的截面慣性矩來傳遞足夠的扭矩，所以主要考慮結(jié)構(gòu)屬性方面；葉片外側(cè)表面易受污染，流動攻角和雷諾數(shù)變動范圍較大，要求葉片外側(cè)具有寬廣的工作范圍和平緩的失速特性，又由于葉片氣動性能隨著雷諾數(shù)、表面粗糙度等條件的變化具有一定的穩(wěn)定性，因此主要考慮氣動性能方面；葉片中部綜合考慮結(jié)構(gòu)和氣動等多方面的因素。在傳統(tǒng)的葉片設(shè)計中，根據(jù)內(nèi)側(cè)、外側(cè)、中部葉片的需求選定三種不同厚度的翼型，并計算葉片不同展向位置的弦長和扭角以此形成三維的葉片外形，其中段與內(nèi)外兩段之間的過渡要采用修型。

眾多學者對葉片設(shè)計中選取翼型的數(shù)量進行了相關(guān)研究。張思達等[2]選擇NACA和DU系列共5種翼型對葉片根部進行優(yōu)化設(shè)計；王曉靜等[3]分析S系列兩種翼型，得出選取S823翼型進行1.5MW風力機葉片設(shè)計的性能更加優(yōu)良的結(jié)論；Li等[4]選擇6種翼型NACA64-618（葉尖 ～75%）DU21、DU25、DU30、DU35（75%～25%）DU40和圓（25%～0）設(shè)計5MW風力機葉片幾何形狀；Yen等[5]基于S809和NACA 63215翼型設(shè)計一種新的混合翼型，從而組合設(shè)計風力機葉片；楊陽等[6]選擇厚度不等的5種DU翼型對風力機葉片進行全局優(yōu)化設(shè)計；劉小龍等[7]選擇三種翼型對風力機葉片和整機的氣動力性能進行分析；李文浩[8]沿葉片展向布置多翼型設(shè)計葉片，其中葉尖部位布置S832翼型，中部布置S831和S830翼型，葉根部位布置S818翼型研究葉片的氣動性能；楊涵[9]基于泛函集成理論針對3種不同厚度的DU系列翼型進行低風速風力機整機氣動與結(jié)構(gòu)一體化設(shè)計研究；王帆[10]選擇NACA系列6種翼型設(shè)計風力機葉片，對其氣動外形性能進行研究，此外研究葉素數(shù)目對風力機葉片氣動外形設(shè)計的影響，選取葉素數(shù)目為35左右可使其氣動性能良好；BHUPINDER Singh等[11]研究了兩種翼型NACA0008和NACA0012H對風力機葉片氣動特性的影響；葉濤等[12]選取3種風力機常用的NACA系列對稱翼型，分析3種對稱翼型對風力機葉片氣動性能的影響，結(jié)果表明葉片上布置的翼型越厚，風力機最佳功率越大，自啟動性能越好；于雷等[13]提出了一種新的混合翼型設(shè)計方法，采用多段翼的形式設(shè)計混合翼型，為風力機一體化設(shè)計提供了新的思路；P García等[14]采用DU91-W2-250（20%～45.6%）、RISOE A1-21（54.4%～65.6%）、NACA 63-418（74.4%～葉尖）三種翼型設(shè)計風力機葉片幾何形狀；Jie Zhu等[15]使用 4種翼型 DU400EU、DU300EU、DU91_W2_250、NACA_64_618設(shè)計 1.5 MW風力機葉片；李仁年等[16]選用7種翼型，以Wilson理論為基礎(chǔ)，采用Matlab對葉片的弦長扭角進行優(yōu)化設(shè)計，成功實現(xiàn)對1.5MW風力機葉片的優(yōu)化設(shè)計；孫雷等[17]選取7種WA系列翼型，根據(jù)翼型相對厚度對其進行順序排列，對1MW海上大型風力發(fā)電機組進行設(shè)計；Widad Yossri等[18]研究了4種翼型，即NACA 0012、NACA 4412、NACA 0015和NACA 4415以及50 cm、75 cm和100 cm三種轉(zhuǎn)子直徑尺寸的組合，并估計了每種配置產(chǎn)生的功率；Amer H Muheisen等[19]使用FX66-S-196 V、FX63-137 S和SG6043三種翼型沿葉片半徑分布，與使用具有相同尺寸的NACA4412單翼型風力機葉片來比較葉片的整體性能。結(jié)果表明，與單翼型的葉片相比，多種翼型水平軸風力機葉片表現(xiàn)出良好的性能，功率系數(shù)增加了約8%；廖書學等[20]選取4種對稱翼型探索了不同厚度翼型對風力機氣動效率和轉(zhuǎn)矩的影響；戰(zhàn)培國等[21]開展了對10MW大型風力機葉片的氣動布局研究，提出了葉片根部采用雙/多翼面氣動布局取代傳統(tǒng)的單一厚翼型，其氣動現(xiàn)象更為復(fù)雜，仍有待進一步深入研究。目前研究中少有分析風力機葉片翼型數(shù)量對其性能的影響，本文依據(jù)葉素動量理論設(shè)計葉片，并采用Qblade軟件分析不同翼型數(shù)量的葉片在額定工況下輸出功率的差異，以此表明葉片翼型數(shù)量對其性能的影響。

2 翼型選擇及氣動性能

風力機一般運行于低雷諾數(shù)下，傳統(tǒng)航空翼型無法滿足風力機對翼型在低雷諾數(shù)下的氣動性能要求，為了提高風力機在某些特定的運行工況和工作環(huán)境下的整體性能，國外針對風力機的運行特點對風電葉片專用翼型開展了理論和實驗研究工作，設(shè)計并開發(fā)了一系列的適用于風電葉片的專用翼型族，其中主要有美國可再生能源實驗室的NREL-S系列翼型[22]、丹麥RIS?國家實驗室的RIS?系列翼型[23]、荷蘭代爾夫特大學、荷蘭航空實驗室和荷蘭國家能源研究中心聯(lián)合開發(fā)的DU系列翼型[24]。

每種系列翼型都會有其不同的氣動性能，有的是基于高升力、高升阻比設(shè)計的，有的是考慮翼型表面粗糙度的影響，有的基于厚度及失速特性設(shè)計翼型等。NACA翼型是現(xiàn)目前風力機葉片中被采用較多的翼型系列之一，因此，本文選取NACA24系列8種翼型為例，驗證模型的有效性，翼型型線示意如圖1所示，8種翼型的最大彎度及其相對位置相同，最大相對厚度逐漸遞增，其中c為翼型弦長。

圖1 翼型型線示意圖

翼型小攻角下的氣動性能可通過計算流體力學方法得到，大攻角下的氣動性能只能依靠試驗測量或半經(jīng)驗公式的失速模型外推得到。通過Qblade軟件計算翼型在小攻角范圍（-5°～15°）內(nèi)的升阻力系數(shù)，如圖2所示，翼型在不同的攻角和雷諾數(shù)下呈現(xiàn)不同的升力系數(shù)和阻力系數(shù)，選取最佳升阻比時對應(yīng)的攻角、升力系數(shù)以及阻力系數(shù)。

圖2 翼型氣動性能

3 空氣動力學模型

當不考慮風輪尾流旋轉(zhuǎn)時，致動盤模型[25]采用流管分析法來預(yù)測理想化的水平軸風力機軸向推力和輸出功率等參數(shù)。假設(shè)將風輪簡化為平面槳盤，沒有輪轂，葉片數(shù)無窮多，葉片旋轉(zhuǎn)時不受摩擦阻力影響，風輪前后氣流都是均勻的定常流，風輪前遠方和風輪后無窮遠處的氣流靜壓相等，作用在風輪上的推力是均勻分布，并且不考慮風輪后的尾流旋轉(zhuǎn)，當軸向誘導(dǎo)因子大于0.5時，致動盤理論將不再適用，需要對其進行經(jīng)驗修正。如圖3所示，下標符號∞為上游無窮遠處，d為致動盤處，w為下游無窮遠處，v為對應(yīng)位置處流速，p為對應(yīng)位置處壓強，A為某位置對應(yīng)截面面積?？紤]流動空氣流經(jīng)致動盤時沿軸向發(fā)生的速度變化，定義軸向誘導(dǎo)因子為a，則流動空氣在致動盤處的軸向速度為：

圖3 致動盤模型

將風輪葉片沿展向劃分為有限個葉素，分析每個葉素上所受的力和力矩，通過積分得到葉片所受的力和力矩。假設(shè)各葉素上氣流均視為二維流動，葉素之間互不影響，并且忽略風輪旋轉(zhuǎn)的三維效應(yīng)。如圖4所示，葉素上相對入流速度為：

圖4 葉素平面上速度和受力分析

4 實例計算與分析

風力機葉片的相關(guān)設(shè)計參數(shù)如表1所示，分別設(shè)計三種不同額定功率的風力機葉片，其中v是額定風速，Cp是風能利用系數(shù)，λ是葉尖速比，R是風輪半徑，R0是輪轂半徑，η1是主傳動系統(tǒng)的總效率，η2是發(fā)電系統(tǒng)的總效率，B是葉片數(shù)，ρ是空氣密度。

表1 風力機葉片的設(shè)計參數(shù)

沿展向?qū)⑷~片分成多段，各個截面處的翼型種類屬于離散取值，不同展向位置r/R處的翼型選取如表2所示，對每種不同額定功率的葉片選擇不同數(shù)量的翼型，并根據(jù)空氣動力學模型計算風力機葉片的軸向誘導(dǎo)因子和周向誘導(dǎo)因子沿展向的分布情況，如圖5所示，可以看出，不同展向位置處的軸向和周向誘導(dǎo)因子不同，從葉根到葉尖位置，軸向誘導(dǎo)因子逐漸增大，周向誘導(dǎo)因子逐漸減??；隨著額定功率的增大，軸向和周向誘導(dǎo)因子逐漸降低。再考慮葉尖損失和輪轂損失計算弦長和扭角，最后得到三維葉片外形。為驗證選取不同數(shù)量翼型葉片的性能，將葉片數(shù)據(jù)輸入仿真軟件Qblade，通過仿真得到葉片在額定工況下的輸出功率，如圖6所示。

表2 葉片各截面翼型選取

圖5 軸向誘導(dǎo)因子和周向誘導(dǎo)因子

圖6 風力機葉片翼型數(shù)量與輸出功率的關(guān)系

如圖6（a），風力機葉片的額定功率為300 W，翼型數(shù)量為3時，風力機葉片在額定工況下輸出功率為264.896 W；翼型數(shù)量分別為4、5、6、7、8時，風力機葉片的輸出功率分別為277.178 W、284.204 W、292.241 W、296.244 W、297.911 W；如圖6（b），風力機葉片的額定功率為1.5 MW，翼型數(shù)量為3時，風力機葉片在額定工況下輸出功率為1.42778 MW；翼型數(shù)量分別為4、5、6、7、8時，風力機葉片的輸出功率分別為1.43096 MW、1.44257 MW、1.4566 MW、1.46524 MW、1.47965 MW；如圖6（c），風力機葉片的額定功率為5 MW，翼型數(shù)量為3時，風力機葉片在額定工況下輸出功率為4.91771 MW；翼型數(shù)量分別為4、5、6、7、8時，風力機葉片的輸出功率分別為4.92405 MW、4.93314 MW、4.94878 MW、4.96373 MW、4.97894 MW?？梢钥闯?，隨著風力機葉片展向翼型數(shù)量的增加，輸出功率在額定功率以下并且數(shù)值逐漸增大。

其中，圖6（a）中翼型數(shù)量為4相較于翼型數(shù)量為3時，輸出功率高出12.282 W；翼型數(shù)量為8相較于翼型數(shù)量為7時高出1.667 W，相比3種翼型的葉片輸出功率高出33.015 W；圖6（b）中風力機葉片展向翼型數(shù)量為8相比7種翼型的葉片輸出功率高出14410 W，相比3種翼型的葉片輸出功率高出51870 W；圖6（c）中翼型數(shù)量為8相比7種翼型的葉片輸出功率高出15210 W，相比3種翼型的葉片輸出功率高出61230 W?？梢园l(fā)現(xiàn)，展向翼型數(shù)量分布逐漸增多，輸出功率的增量值逐漸減小。

3種不同額定功率下，風力機葉片展向翼型數(shù)量為4相較于展向翼型數(shù)量為3時，輸出功率的增量與額定功率的比值分別為：0.41%、0.212%、0.1268%。展向翼型分布數(shù)量為8相較于翼型數(shù)量為3時，增加的輸出功率占比分別為：11.005%、3.458%、1.2246%?？梢钥闯?，在相同展向翼型數(shù)量增加的情況下，隨著風力機額定功率的增大，增加的輸出功率占比逐漸減小。

5 結(jié)論

本文以NACA24系列翼型為例，基于葉素-動量理論，考慮葉尖和輪轂損失，分析三種不同額定功率下風力機葉片的展向翼型數(shù)量對其性能的影響，結(jié)果表明，風力機葉片的額定功率為300 W時，翼型數(shù)量為8的風力機葉片輸出功率相比翼型數(shù)量為3的葉片輸出功率高出0.55%；風力機葉片的額定功率為1.5 MW時，翼型數(shù)量為8的風力機葉片的輸出功率相比3種翼型的葉片輸出功率高出3.458%；風力機葉片的額定功率為5 MW時，翼型數(shù)量為8的風力機葉片的輸出功率相比3種翼型的葉片輸出功率高出1.2246%。對于水平軸風力機葉片設(shè)計時，選取展向翼型數(shù)量少，風力機葉片輸出功率降低，原因在于展向翼型數(shù)量少導(dǎo)致葉片表面不連續(xù)、不光滑，選取展向翼型數(shù)量越多，風力機葉片的輸出功率呈現(xiàn)遞增趨勢，展向翼型數(shù)量的增多保證風與葉片充分相互作用，從而提高風能利用率并且減少后期的修型。