劉 陽

(大唐新能源陜西公司,陜西 西安710065)

0 引言

近幾年,隨著可開發(fā)風(fēng)資源區(qū)域的多樣化,導(dǎo)致同一類型機(jī)組在不同環(huán)境下,運(yùn)行狀態(tài)存在較大差異。統(tǒng)計(jì)分析發(fā)現(xiàn),風(fēng)電機(jī)組運(yùn)行在較低空氣密度時(shí),機(jī)組葉片根部在額定風(fēng)速附近出現(xiàn)失速,導(dǎo)致機(jī)組滿發(fā)功率風(fēng)速向后推移,嚴(yán)重影響機(jī)組穩(wěn)定運(yùn)行[1]。

為了解決葉片在低空氣密度下的失速問題,科研人員進(jìn)行了一系列的研究。Butterfied[2]首次測量了水平軸風(fēng)力機(jī)的壓力分布,對動態(tài)失速進(jìn)行了定量研究;Leishman[3]開發(fā)了能夠描述動態(tài)失速翼型的升力,轉(zhuǎn)矩和阻力特性的模型;由于風(fēng)力機(jī)葉片失速發(fā)生在葉片根部,德國ENERCON公司對風(fēng)力機(jī)葉片的根部采用了翼型段代替葉根的過度段和圓柱段,通過對葉片的改型設(shè)計(jì),減小葉片根部的失速分離;邢作霞[4]、姚興佳[5]等采用獨(dú)立變槳技術(shù)以及智能控制算法等來降低機(jī)組載荷,減少葉片失速。劉庭瑞[6]針對風(fēng)力機(jī)葉片失速振顫斷裂失效問題,采用最優(yōu)模糊比例積分微分(proportional integral derivative,PID)控制方法和徑向吉函數(shù)神經(jīng)網(wǎng)絡(luò)PID控制方法進(jìn)行分析,為風(fēng)力機(jī)葉片失速非線性氣彈變槳提供了新思路。陳巖[7]通過對風(fēng)力機(jī)葉片翼型進(jìn)行穩(wěn)態(tài)氣動力分析,改善了失速工況下的負(fù)阻尼。俞國華[8]通過數(shù)值計(jì)算與理論分析研究了風(fēng)力機(jī)靜態(tài)失速和動態(tài)失速問題,對理解失速的復(fù)雜流動機(jī)理提供了參考。

本文在前人研究的基礎(chǔ)上,針對葉片在低空氣密度下,額定風(fēng)速附近失速問題,分別采用提前變槳策略以及葉片根部加裝根箱兩種解決方案,并采用三維數(shù)值模擬對這兩種方案進(jìn)行驗(yàn)證,通過對這兩種方案的對比分析,得出解決葉片失速的最優(yōu)解決方案。

1 風(fēng)力機(jī)葉片失速理論分析

1.1 翼型失速理論分析

風(fēng)力機(jī)葉片翼型的攻角α為來流風(fēng)速方向與葉片翼型弦線之間的夾角。由圖1可知,葉片翼型截面的相對速度vrel為風(fēng)輪平面的軸向速度(1-α)v0和切向速度(1+α')ωr的合成,并且攻角α可由vrel和旋轉(zhuǎn)平面之間的夾角φ與槳距角θ表示,即：α=φ-θ。

圖1 葉片截面翼型圖Fig.1 Blade Section Airfoil Diagram

根據(jù)在翼型上升力阻力的定義,得到單位長度上的升力L和阻力D為：

式中：ρ為空氣密度,kg/m3;c為葉片弦長,m;Cl為升力系數(shù);Cd為阻力系數(shù)。

由公式(1)(2)可知,翼型的升力和阻力除與來流風(fēng)速以及空氣密度有關(guān)外,還與Cl和Cd有關(guān),且不同的翼型具有不同升阻力系數(shù),典型的Cl和Cd曲線如圖2所示。由圖2可知,當(dāng)攻角大于aM后,升力系數(shù)逐漸減小,阻力系數(shù)逐漸增加,從而導(dǎo)致升力迅速下降,阻力迅速上升。aM點(diǎn)稱為“失速點(diǎn)”,翼型在經(jīng)過失速點(diǎn)后,出現(xiàn)失速[4]。

1.2 低空氣密度下風(fēng)力機(jī)葉片失速

目前,主流的風(fēng)力機(jī)為變槳變速型風(fēng)力機(jī)。變槳變速型風(fēng)機(jī)在額定風(fēng)速以下時(shí),風(fēng)力機(jī)組槳距角始終維持在最佳槳距角(一般為0°),通過控制風(fēng)輪轉(zhuǎn)速來調(diào)節(jié)功率;額定風(fēng)速以上,通過不斷調(diào)整攻角避開其失速點(diǎn)以維持額定的功率輸出。并且同一臺風(fēng)力機(jī)在相同的最優(yōu)槳距角情況下,空氣密度不同,對應(yīng)的功率曲線不同。

圖2 Cl和Cd隨功角變化曲線Fig.2 Curves of Cl and Cd varying with work angle

風(fēng)力機(jī)組在實(shí)際運(yùn)行過程中,空氣密度越低,要達(dá)到額定功率,需要的來流風(fēng)速越大。而額定風(fēng)速前,風(fēng)輪轉(zhuǎn)速已達(dá)到額定轉(zhuǎn)速保持不變,隨著來流風(fēng)速的增加,φ角增加,但θ保持不變,從而α不斷增加,導(dǎo)致葉片攻角達(dá)到并超過失速點(diǎn),引起葉片失速。圖3為空氣密度為1.2 kg/m3變化至1 kg/m3時(shí)的功率曲線,由圖3可知,在空氣密度變化至1 kg/m3,風(fēng)力機(jī)功率曲線在13 m/s附近,明顯發(fā)生拐點(diǎn),此點(diǎn)為失速點(diǎn)。

圖3 不同空氣密度下的功率曲線Fig.3 power curve under different air density

2 風(fēng)力機(jī)葉片失速控制方案

2.1 葉片失速最佳槳距角調(diào)節(jié)方案

在低空氣密度下,為了避免風(fēng)力機(jī)在額定風(fēng)速前失速,可采取提前變槳控制策略,通過增大槳距角,達(dá)到減小葉片攻角,控制葉片失速的目的。同時(shí),為了避免葉片自身將造成較大的風(fēng)能損失,根據(jù)如圖4所示的在不同槳距角下的Cp-λ曲線可知,同一葉尖速比下,槳距角越大,風(fēng)能利用系數(shù)越低。因此綜合考慮,變槳角度最大不能超過2°。

圖4 在不同槳距角下的Cp-λ曲線Fig.4 Cp-λcurve at different pitch angles

為了對變槳效果進(jìn)行驗(yàn)證,利用Bladed軟件,對所研究的1.5 MW機(jī)組,在空氣密度1 kg/m3,變槳角度分別為0°、0.5°、1.0°、1.5°以及2.0°情形下的功率曲線進(jìn)行模擬,得到功率曲線如圖5所示。

圖5 不同槳距角下靜態(tài)功率曲線Fig.5 Static power curves at different pitch angles

由圖5可知,槳距角為0°時(shí),風(fēng)速為13 m/s附近時(shí),風(fēng)力機(jī)組功率出現(xiàn)拐點(diǎn),此點(diǎn)為葉片的失速點(diǎn)。隨著槳距角在0°到2°的變化過程中,風(fēng)力機(jī)功率曲線拐點(diǎn)明顯變平滑,其功率值提高。

2.2 葉片根部加裝根箱

根據(jù)經(jīng)驗(yàn),風(fēng)力機(jī)發(fā)生失速的區(qū)域?yàn)槿~片的葉根段,除改變槳距角外,風(fēng)力機(jī)葉片通過葉根改型,在物理上阻斷葉根的流動分離,可有效控制葉片的失速。學(xué)者們已經(jīng)通過將葉片根部圓柱度和過多段改為翼型段來增加功率捕獲減小失速分離,但其改善失速分離效果并不明顯,本文在其基礎(chǔ)上,將葉根圓柱段改為翼型段,同時(shí)加裝了翼刀結(jié)構(gòu),其結(jié)構(gòu)如圖6所示。

由圖6可知,葉片結(jié)構(gòu)改型方式為：葉片根部改型為在葉片展向0.04R～0.25R(R為葉片長度)區(qū)域內(nèi),保持葉片前緣形狀不變,在后緣加裝增功組件,將原葉根的圓柱段改為翼型段,同時(shí)加裝翼刀。其翼型形狀為原葉片在最大弦線位置處的翼型,扭角大小與原葉片在同位置處的扭角相等。第1個(gè)翼刀距離葉根約2 m。

圖6 風(fēng)力機(jī)葉片根部改型結(jié)構(gòu)Fig.6 Wind turbine blade root modification structure

3 計(jì)算模型建立

3.1 模型建立

本文根據(jù)風(fēng)電機(jī)組的性能及尺寸建立計(jì)算模型,對兩種方案進(jìn)行驗(yàn)證。計(jì)算區(qū)域分為旋轉(zhuǎn)域和靜止域。旋轉(zhuǎn)域?yàn)榘~片在內(nèi)旋轉(zhuǎn)部件。靜止域?yàn)榘?機(jī)艙等在內(nèi)的靜止部件。計(jì)算模型圖如圖7所示,計(jì)算域長、寬、高分別為8倍風(fēng)輪直徑、3倍風(fēng)輪直徑、3倍風(fēng)輪直徑。

圖7 風(fēng)電機(jī)組計(jì)算模型Fig.7 Model of wind turbine calculation

風(fēng)電機(jī)組基本參數(shù)為：額定功率1.5 MW,葉片個(gè)數(shù)為3個(gè),葉片長度43.5 m標(biāo)準(zhǔn)空氣密度下的額定風(fēng)速12 m/s,額定轉(zhuǎn)速為17.6 r/min,風(fēng)輪仰角為5°,風(fēng)向?yàn)樯巷L(fēng)向。

3.2 網(wǎng)格劃分與無關(guān)性驗(yàn)證

運(yùn)用ICEM軟件對風(fēng)力機(jī)計(jì)算流域按區(qū)域進(jìn)行結(jié)構(gòu)網(wǎng)格劃分,如圖8所示為劃分好的風(fēng)力機(jī)網(wǎng)格。為了確保葉根區(qū)域的計(jì)算精度,對葉片根部進(jìn)行網(wǎng)格加密,改型后根部網(wǎng)格如圖9所示。

數(shù)值模擬中,為了消除網(wǎng)格數(shù)對計(jì)算結(jié)果的影響,分別對698萬、867萬、1 078萬、1 244萬和1 479萬網(wǎng)格數(shù)的模型進(jìn)行了模擬,如圖10所示。由圖10可知：隨著網(wǎng)格數(shù)的增加,計(jì)算功率逐漸穩(wěn)定。在網(wǎng)格超過1 124萬后,功率的變化很小。綜合考慮模擬精度和計(jì)算時(shí)長的影響,最終確定網(wǎng)格總數(shù)為1 124萬,根部改型后的網(wǎng)格為1 346萬。

圖8 風(fēng)力機(jī)輪轂周圍網(wǎng)格Fig.8 Grid around wind turbine hub

圖9 改型后葉片根部網(wǎng)格Fig.9 Modified blade root mesh

圖10 網(wǎng)格無關(guān)性驗(yàn)證Fig.10 Grid independence verification

3.3 計(jì)算方法和邊界條件設(shè)置

計(jì)算域入口處的來流風(fēng)速,出口的相對壓力設(shè)定為0,交界面采用動靜轉(zhuǎn)子邊界,地面、葉片表面、塔筒、機(jī)艙和導(dǎo)流罩均采用無滑移壁面。計(jì)算軟件為ANSYS CFX,湍流模型采用對壁面的流固分離具有很高的計(jì)算精度的SSTk-ε模型。對流向采用高階精度格式,隱式求解,同時(shí),對風(fēng)力機(jī)的功率和效率進(jìn)行監(jiān)測,當(dāng)功率和效率誤差小于5%且平均殘差小于10-4時(shí)認(rèn)為收斂。

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

4.1 功率分析

通過對機(jī)組運(yùn)行數(shù)據(jù)的統(tǒng)計(jì)分析得知：風(fēng)速在13 m/s時(shí),葉片的流動分離最大。同時(shí)綜合bladed分析結(jié)果,本次分析在13 m/s風(fēng)速,1 kg/m3空氣密下,分別對原葉片槳距角為0°、1.5°以及根部改型葉片在0°槳距角進(jìn)行模擬。

表1所示為各計(jì)算工況下的風(fēng)力機(jī)組功率,由表1可知,葉片在改變槳距角后,功率增加越1.84%,葉片根部改型后,功率增加2.14%。

表1 功率隨槳距角變化表Table 1 Variation of power with pitch angle

4.2 葉片根部表面流線分析

在風(fēng)力機(jī)旋轉(zhuǎn)葉輪的作用下,氣流在通過葉片前緣位置后,流動速度開始增大。由于空氣的粘性作用,當(dāng)氣流速度達(dá)到最大值后,葉片吸力面的氣流速度逐漸減小,同時(shí),逆壓梯度逐漸增加,當(dāng)吸力面上層流動無法克服逆壓梯度時(shí),將在葉片吸力面上發(fā)生倒流或逆流,使得氣體在吸力面上發(fā)生分離[9-11]。氣流的分離會造成大量的紊流和渦,使升力大量受損。同時(shí),氣流分離越靠近前緣位置,所產(chǎn)生的壓差阻力會越大,對氣動性能影響越明顯。風(fēng)力機(jī)在運(yùn)行過程中,葉片根部吸力面上產(chǎn)生的流動分離不可避免,因此,有效的減小葉根氣流分離,能夠顯著提升葉根的風(fēng)能捕獲。

圖11 葉片表面極限流線分布Fig.11 Distribution of limit streamline on blade surface

由圖11可知,來流風(fēng)速為13 m/s(槳距角β=0°)時(shí),從葉片根部起始位置開始便出現(xiàn)了氣體的流動分離。同時(shí)形成了一條清晰的分離線,沿著葉片展向,分離線逐漸向后緣位置靠近,在展向約0.25R位置處時(shí),分離線在葉片后緣位置處消失。來流風(fēng)速為13 m/s(槳距角β=1.5°)時(shí),對比β=0°,槳距角增加,葉片吸力面上流動分離線在葉根的起始位置遠(yuǎn)離前緣位置,且葉片根部的分離區(qū)域變小;葉片根部改型后,來流風(fēng)速為13 m/s(槳距角β=0°)時(shí),葉片根部流動分離的起始位置與原葉片相同。但葉片表面流動在接近第1個(gè)翼刀之后沿展向消失,有效控制了葉片。

由上分析可知,改變槳距角與葉根改型均能減小葉片根部的流動,但根部改型的效果更為明顯,改變槳距角僅減小了葉片表面流動分離。而葉片根部改型后,雖使根部出現(xiàn)了較小的附著渦,但顯著改善了葉根區(qū)域的流動分離,有效地將葉片根部的流動分離控制在第1個(gè)翼刀葉根側(cè)。

4.3 經(jīng)濟(jì)性分析

從分析結(jié)果可知,通過調(diào)節(jié)槳距角來控制葉片的在低空氣密度下的流動分離,僅需要對機(jī)組的主控制程序進(jìn)行修改,無需在機(jī)組上增加其他硬件設(shè)備。而針對葉片根部改型,需要在設(shè)計(jì)階段對葉片進(jìn)行優(yōu)化改型,對于已建成風(fēng)場需要很大的技改費(fèi)用,因此針對低空氣密度下的葉片失速問題,已建成的風(fēng)場易采用調(diào)節(jié)槳距角解決,而針對未建設(shè)風(fēng)場,采用根部改型后的效果更好。

5 結(jié)論

針對低空氣密度,風(fēng)電機(jī)組在接近額定風(fēng)速附近運(yùn)行,葉片根部出現(xiàn)失速問題,本文提出了提前變槳策略以及葉片根部加裝根箱兩種方案解決方案,同時(shí)采用三維數(shù)值模擬對兩種方案進(jìn)行驗(yàn)證。得到如下結(jié)論：

(1)采用兩種方案均能減小葉片表面失速分離,但加裝根箱的改善失速分離的效果更為明顯。

(2)采用提前變槳控制策略可提升機(jī)組功率約1.78%,葉片根部加裝根箱可提升機(jī)組功率約2.14%。

(3)對于已建成的風(fēng)場,從經(jīng)濟(jì)型考慮,采用提前變槳控制策略解決葉片失速問題效果更好。