兆瓦級風(fēng)力機葉片的改進型BEM理論研究

王會剛，楊增帥，劉曉雯

（1.唐山學(xué)院河北省智能裝備數(shù)字化設(shè)計及過程仿真重點實驗室，河北 唐山 063000；2.孚威吉自動化裝備河北有限公司，河北 涿州 072700）

0 引言

當(dāng)代社會，人們的生產(chǎn)生活對電能的依賴程度遠(yuǎn)遠(yuǎn)超過歷史上的任何階段。同時，由于傳統(tǒng)能源電力生產(chǎn)方式給人類賴以生存的環(huán)境造成了巨大的、不可逆轉(zhuǎn)的污染。風(fēng)能作為清潔可再生的新能源中的一員，具有容量大、技術(shù)相對成熟的優(yōu)勢，在各國政府的大力支持下得到了快速發(fā)展。

葉片是風(fēng)力機捕獲風(fēng)能的核心部件，其成本雖然只占風(fēng)力機總成本的1/10左右，但科學(xué)、合理的設(shè)計卻能使傳動部件降低高達(dá)1/5的載荷，進而大幅降低風(fēng)力機制造成本，最終降低能源價格[1]。葉片的形狀異常復(fù)雜，除了具有流線翼型的特點之外，空間扭曲、沿長度的截面突變等因素導(dǎo)致葉片建模十分復(fù)雜[2]。科研人員針對葉片氣動外形的設(shè)計，所使用的方法通常有圖解法、等升力系數(shù)法、等弦長法、Glauert法、Wilson法、簡化風(fēng)車設(shè)計法、經(jīng)典BEM法，由于種種原因，風(fēng)力機向低風(fēng)速、大型化發(fā)展趨勢日益明顯，部分理論方法不能滿足這一要求，因此對經(jīng)典理論進行一定程度的改進或者優(yōu)化顯得尤為必要。戴巨川等[3]針對2 MW風(fēng)力機模型通過建立葉片相關(guān)氣動外形參數(shù)控制方程，基于BEM法及CFD仿真分析，給出了參數(shù)化設(shè)計法；Raju Bharath-Koratagere-Srinivasa[4]應(yīng)用經(jīng)典BEM理論設(shè)計了1.6 MW風(fēng)力機葉片，但是因求解的弦長為9.5～4.9 m，不得不采用經(jīng)驗值代替；陳進等[5]結(jié)合一維BEM法和考慮軸向力和切向力的葉尖損失公式，給出了新的軸向誘導(dǎo)因子和切向誘導(dǎo)因子計算公式；Zhu Weijun等[6]利用自己改進并簡化的葉素動量理論設(shè)計了20 MW風(fēng)力機葉片氣動參數(shù)；汪泉等[7]采用了修正的BEM理論設(shè)計了全新的風(fēng)力機葉片并進行了仿真驗證。

大型風(fēng)力機乃至超大型風(fēng)力機必將占據(jù)未來風(fēng)電發(fā)展趨勢的主流，為有效提高葉片捕獲風(fēng)能的效率，針對大型風(fēng)力機葉片理論的深入研究顯得尤為重要。結(jié)合當(dāng)前兆瓦級風(fēng)力機領(lǐng)域葉片理論發(fā)展，研究了葉片理論計算方法并推導(dǎo)出一種軸向誘導(dǎo)因子函數(shù)表達(dá)式。基于經(jīng)典動量理論，引入葉根及葉尖損失因子的Prandtl修正和Glauert修正，結(jié)合R. Lanzafame等基于BEM理論和角動量理論推導(dǎo)出的推力及轉(zhuǎn)矩表達(dá)式，推導(dǎo)出軸向誘導(dǎo)因子及切向誘導(dǎo)因子表達(dá)式，再對軸向誘導(dǎo)因子進行Spera修正并進行必要簡化；考慮了制造等因素的影響，進一步對弦長和扭角進行擬合及線性化處理，重構(gòu)葉片幾何參數(shù)，從而實現(xiàn)了風(fēng)力機葉片幾何模型構(gòu)建。同時，也提出了適合大型風(fēng)力機葉片的翼型選擇原則。以1.5 MW風(fēng)力機葉片為實例研究，充分驗證了改進型BEM法的正確性，葉片各項特征參數(shù)較經(jīng)典Glauert法設(shè)計的葉片具有明顯技術(shù)優(yōu)勢，更適合在兆瓦級風(fēng)力機上應(yīng)用，也為葉片設(shè)計理論的進一步發(fā)展奠定了基礎(chǔ)。

1 改進型BEM理論設(shè)計法

BEM理論也稱葉素動量理論，是將葉素理論和動量理論結(jié)合而生成的一種理論方法，通常被應(yīng)用到風(fēng)力機葉片的設(shè)計當(dāng)中。但是經(jīng)典設(shè)計理論未考慮葉尖損失系數(shù)及葉根損失系數(shù)的影響，因為這兩個損失系數(shù)在大型風(fēng)力機設(shè)計中對風(fēng)能利用系數(shù)影響程度比較大，間接影響著風(fēng)輪直徑及葉片參數(shù)等。

1.1 經(jīng)典BEM理論

1.1.1 翼型扭角計算

假設(shè)風(fēng)輪以角速度Ω旋轉(zhuǎn)，來流風(fēng)速為v∞，取距離風(fēng)輪旋轉(zhuǎn)中心r處dr圓環(huán)為研究對象，如圖1所示。因渦流系的影響，風(fēng)輪附近流場流速發(fā)生了變化。根據(jù)渦流理論，經(jīng)過風(fēng)輪截面后的氣流軸向速度為v∞（1-a）；由角動量理論可知，此時的切向速度為Ωr·（1+a′），其中a為軸向誘導(dǎo)因子，a′為切向誘導(dǎo)因子，此處合成速度為W，由圖2可得到下式：

圖1 葉素掃略環(huán)

圖2 葉素速度合成

式中：Re為雷諾數(shù)；c為葉素弦長，m；μ為第一黏性系數(shù)，標(biāo)準(zhǔn)狀況下取1.698×10-5kg/（s·m）；ρ為風(fēng)場空氣密度，kg/m3。

葉素翼型的最佳局部攻角通常是根據(jù)式（3）計算的雷諾數(shù)由相關(guān)軟件輸出攻角和升阻比關(guān)系，按該翼型最大升阻比對應(yīng)的攻角（最佳攻角）進行計算。NACA4415升阻比曲線如圖3所示，可確定最佳攻角為6°。

圖3 NACA4415翼型特征

實際設(shè)計中，往往使攻角在最佳攻角附近，主要是考慮了制造難度，選擇最佳攻角時，應(yīng)盡量使各葉素扭角呈流線型布置。由式（1）入流角和半徑的關(guān)系可知，入流角隨半徑r的遞增（即由葉根向葉尖方向）而減小。如果按最大升阻比選擇的話，往往各葉素攻角不是單調(diào)遞減而是往復(fù)波動，扭角會出現(xiàn)局部遞增而出現(xiàn)波浪形扭曲的現(xiàn)象，顯然這樣的結(jié)果增加了制造難度。

1.1.3 角推力系數(shù)及切向牽引力系數(shù)的計算

利用式（1）計算出的入流角φ及選取的最佳攻角對應(yīng)的升力系數(shù)Cl及阻力系數(shù)Cd，根據(jù)葉素理論得到角推力系數(shù)及切向牽引力系數(shù)為：

式中：Cn為推力系數(shù)；Ct為切向牽引力系數(shù)；Cl為升力系數(shù)；Cd為阻力系數(shù)。

1.2 改進型BEM理論誘導(dǎo)因子的推導(dǎo)

經(jīng)典BEM理論并未考慮葉尖損失系數(shù)和根部損失系數(shù)，但是研究表明該項系數(shù)對風(fēng)力機尤其大型風(fēng)力機性能影響很大。因此，Prandtl等對經(jīng)典BEM理論進行了改進和優(yōu)化，引入了葉尖損失系數(shù)及葉根損失系數(shù)，如下式：

式中：Ft為葉尖損失系數(shù)；N為葉片數(shù)；Fr為根部損失系數(shù)；R為輪轂與葉根連接面旋轉(zhuǎn)半徑（如圖4），m；rhub為葉根安裝面回轉(zhuǎn)半徑（如圖4），m。

圖4 風(fēng)輪半徑與葉片長度、輪轂半徑的關(guān)系

因此，考慮葉根及葉尖損失系數(shù)后的總損失系數(shù)為

由于葉根和葉尖損失系數(shù)的影響，計算推力和轉(zhuǎn)矩時需要考慮損失系數(shù)導(dǎo)致的實際輸出功率降低，故動量理論的推力及轉(zhuǎn)矩可表示為：

當(dāng)軸向誘導(dǎo)因子a＞0.4時，動量定理不再適用。因此，Spera進一步對BEM理論進行了修正[11]，從而實現(xiàn)了對軸向誘導(dǎo)因子求解理論做了補充。對于a≤0.2時軸向誘導(dǎo)因子無需修正，在此基礎(chǔ)上對軸向誘導(dǎo)因子求解公式進一步簡化：當(dāng)軸向誘導(dǎo)因子a＞0.2時，軸向誘導(dǎo)因子更新為：

根據(jù)上述理論推導(dǎo)，改進型BEM理論迭代求解軸向誘導(dǎo)因子和切向誘導(dǎo)因子流程圖如圖5所示。

圖5 改進型BEM理論迭代流程

迭代求解初始化軸向誘導(dǎo)因子和切向誘導(dǎo)因子，通常通過下式確定[12]：

式中，λ為葉尖速比。

此外，迭代容差選取時，如果容差過小則求解的軸向誘導(dǎo)因子過大；如果容差過大則求解的軸向誘導(dǎo)因子過小。

2 葉片翼型選擇原則

風(fēng)力機的功率特性，一定程度上取決于翼型的選擇，因此翼型的設(shè)計及其優(yōu)化設(shè)計一直是國內(nèi)外學(xué)者的研究重點。早期風(fēng)力機翼型多取自于航空領(lǐng)域內(nèi)的諸如滑翔翼等低速翼型，如WortmannFX-77、NASALS翼型等。后來各國相繼又研發(fā)很多翼型，如美國的NREL翼型，瑞典的FFA-W翼型，荷蘭的DU翼型，丹麥的Risφ-A1、Risφ-P、Risφ-B1翼型等。

翼型的選擇通常考慮以下幾點：

1）葉根承受了葉片全部載荷，諸如拉力、壓力、橫向剪切力和扭轉(zhuǎn)剪切力等，應(yīng)力狀態(tài)十分復(fù)雜，故而對強度和剛度要求極高，應(yīng)選用抗彎截面系數(shù)比較大的厚翼型，如WORTMANN FX 77-W-258翼型（最大厚度為26.15%）。

2）風(fēng)力機約有75%的電能來自于翼展方向上靠近葉尖的葉片長度的75%的范圍內(nèi)，因此該部分適宜選用較薄的翼型，便于高效地獲取風(fēng)能。如NACA4415翼型（最大厚度為15%）等。

3）由于現(xiàn)代大型風(fēng)力機葉片比較長，通常沿翼展方向分為若干功能截面（通常等分），不同截面根據(jù)所在位置需要完成的功能而選擇適宜的翼型，最大程度地提高葉片的風(fēng)能利用系數(shù)，同時保證在較低的固有頻率不發(fā)生失速顫振等異常現(xiàn)象。

3 1.5 MW風(fēng)力機葉片實例對比研究

分別應(yīng)用前述改進型BEM理論數(shù)學(xué)模型以及經(jīng)典Glauert理論法，設(shè)計1.5 MW風(fēng)力機葉片并進行對比研究，驗證改進型BEM理論在兆瓦級風(fēng)力機應(yīng)用領(lǐng)域的可行性及正確性。

3.1 設(shè)計參數(shù)

參照某公司1.5 MW風(fēng)力機組數(shù)據(jù)，整理為本研究對象的原始設(shè)計參數(shù)，如表1所示。

表1 1.5 MW風(fēng)力機組設(shè)計參數(shù)

根據(jù)以上技術(shù)參數(shù)，由一維動量理論計算得出，風(fēng)輪直徑D=77 m，葉片長度Lb=37.5 m。

依據(jù)前述翼型選擇原則，沿葉片翼展方向依次選擇WORTMANN FX 77-W-258翼型（最大厚度為26.15%，代號為AF1）、NACA4418翼型（最大厚度為18.02%，代號為AF2）、NACA4415翼型（最大厚度為15%，代號為AF3）以及NACA4412翼型（最大厚度為12.02%，代號為AF4）。

風(fēng)力機葉片外形非常復(fù)雜，除了具有流線外形，高度的空間扭曲和弦長的大幅度變化，加之葉片長度與最大弦長比值很大，葉片制造非常困難，因此需要在一定程度上進行優(yōu)化。如前所述，在選取最佳攻角時，對扭角進行了線性化處理。但弦長對制造難度的影響，比扭角要大得多。迭代法求解的弦長，往往變化十分劇烈，因此針對弦長進行了多項式擬合。

3.2 設(shè)計結(jié)果及分析

應(yīng)用上述理論迭代求解，并對弦長進行4次多項式擬合，完成了經(jīng)典BEM法葉片設(shè)計工作；同時基于前述機組參數(shù)，應(yīng)用經(jīng)典Glauert法得出了葉片參數(shù)。兩種方法計算的葉片參數(shù)如表2所示。

表2 1.5 MW風(fēng)力機葉片翼型參數(shù)表

通過表2可看出，在同樣設(shè)計參數(shù)下，改進型BEM理論法設(shè)計的葉片額定輸出轉(zhuǎn)矩為238.1 kN·m，滿足設(shè)計值；Glauert法設(shè)計的葉片額定輸出轉(zhuǎn)矩為89.2 kN·m，僅為前者輸出轉(zhuǎn)矩的37.46%；如圖6所示，兩種計算方法各段轉(zhuǎn)矩值隨半徑的變化趨勢基本相同，但在12.25～34.75 m之內(nèi)轉(zhuǎn)矩波動分別為11.39%和17.3%，表明改進型BEM法比經(jīng)典Glauert法有著更為穩(wěn)定的轉(zhuǎn)矩輸出，從上述分析得出，Glauert法在該種類型風(fēng)力機設(shè)計中的功率顯然遠(yuǎn)遠(yuǎn)達(dá)不到設(shè)計值，同時凸顯了改進型BEM 法在輸出轉(zhuǎn)矩上的高穩(wěn)定性。

圖6 轉(zhuǎn)矩分布圖

改進型BEM 法得出的扭角范圍為0.48° ～8.22° ，而經(jīng)典Glauert法計算的轉(zhuǎn)矩范圍為-3.38°～-1.33°，扭角小于零。由圖2所示速度合成原理可知，出現(xiàn)這種現(xiàn)象的原因之一是經(jīng)典Glauert法計算的切向誘導(dǎo)因子嚴(yán)重偏大等因素導(dǎo)致的；由圖7分析可知，經(jīng)典Glauert法計算的扭角值并未像改進型BEM法計算的扭角值一樣單調(diào)遞減，而是先遞減后遞增；由式（1）、式（2）可知，隨著半徑的增加，扭角應(yīng)單調(diào)遞減。上述現(xiàn)象也表明了改進型BEM法較經(jīng)典Glauert法更適合在兆瓦級風(fēng)力機領(lǐng)域應(yīng)用。

圖7 扭角分布圖

弦長分布曲線如圖8 所示，兩種理論設(shè)計法得出的葉尖部位的弦長值相差不大，沿半徑方向變化趨勢基本相同；但是對大型葉片來說，經(jīng)典Glauert法計算得出的弦長比改進型BEM法的計算結(jié)果偏小，葉素理論表明葉片輸出轉(zhuǎn)矩與弦長成正比，這也是導(dǎo)致經(jīng)典Glauert法計算出的轉(zhuǎn)矩嚴(yán)重偏低的原因。同時對于大型葉片來說，葉根處的強度和剛度也直接取決于弦長的大小；從制造角度分析，如表2所示，由截面2到截面3的如此大幅度的弦長漸變，增加了葉片的制造難度。從這個角度分析，改進型BEM法在大型風(fēng)力機葉片設(shè)計上表現(xiàn)出優(yōu)越的性能。

圖8 弦長分布圖

綜上所述，改進型BEM 理論法設(shè)計的1.5 MW風(fēng)力機葉片特征及各參數(shù)的分布規(guī)律如下：軸向誘導(dǎo)因子和切向誘導(dǎo)因子沿翼展方向分布特征如圖9所示，兩個誘導(dǎo)因子在葉根處比較大，隨半徑的增加大變化率下降，因子值減小，并逐漸趨于穩(wěn)定。對于軸向誘導(dǎo)因子而言，除葉根處比較大之外，其余部位都在0.2左右；切向誘導(dǎo)因子變化幅度比較大，葉根處很大，葉尖處接近0。

圖9 誘導(dǎo)因子分布圖

弦長及扭角（槳距角）翼展方向分布特征如圖10 所示，總體上均隨半徑的增加而減小，變化趨勢相同；其中，扭角在葉尖處接近0°，由式（1）的關(guān)系看出，此種現(xiàn)象充分考慮了半徑對入流角的影響；弦長在葉尖處達(dá)到了最小值，極大地降低葉片的轉(zhuǎn)動慣量，減小了葉根處所承受的彎矩。

圖10 弦長及扭角分布圖

4 結(jié)論

針對兆瓦級風(fēng)力機葉片設(shè)計理論，基于經(jīng)典動量理論，引入葉根及葉尖損失因子的Prandtl修正和Glauert修正，結(jié)合R. Lanzafame等基于BEM理論和角動量理論推導(dǎo)出的推力及轉(zhuǎn)矩表達(dá)式，推導(dǎo)出軸向誘導(dǎo)因子及切向誘導(dǎo)因子表達(dá)式，再對軸向誘導(dǎo)因子進行Spera修正并進行必要簡化；同時以1.5 MW風(fēng)力機葉片為實例研究，對比研究了改進型BEM法及經(jīng)典Glauert法在兆瓦級風(fēng)力機葉片設(shè)計領(lǐng)域的應(yīng)用特點，為風(fēng)力機葉片理論的進一步發(fā)展奠定了基礎(chǔ)。本文主要特點如下：

1）基于Prandtl、Glauert、Spera等的設(shè)計理論，結(jié)合當(dāng)前大型風(fēng)力機設(shè)計理論的發(fā)展，推導(dǎo)出軸向誘導(dǎo)因子及周向誘導(dǎo)因子函數(shù)表達(dá)式，并對Spera等提出的修正公式進一步簡化，完善了軸向誘導(dǎo)因子求解公式。

2）針對兆瓦級風(fēng)力機葉片特點，給出了大型風(fēng)力機葉片設(shè)計中翼型的選擇原則，為葉片設(shè)計提供了參考。

3）以1.5 MW風(fēng)力機葉片為設(shè)計實例，分別應(yīng)用改進型BEM法及經(jīng)典Glauert法完成該葉片的參數(shù)求解。分別從轉(zhuǎn)矩、扭角、弦長等關(guān)鍵參數(shù)進行對比研究，結(jié)果表明改進型BEM理論的正確性，且較經(jīng)典Glauert法更適合兆瓦級風(fēng)力機葉片的設(shè)計和應(yīng)用。

4）總結(jié)了葉片特征參數(shù)沿翼展方向的分布規(guī)律。系統(tǒng)地總結(jié)了軸向誘導(dǎo)因子、切向誘導(dǎo)因子、弦長、扭角沿翼展的分布規(guī)律，為深入研究兆瓦級風(fēng)力機葉片特征機理奠定了基礎(chǔ)。