大型三槳葉上風(fēng)力機風(fēng)剪切塔影效應(yīng)的分析

大型三槳葉上風(fēng)力機風(fēng)剪切塔影效應(yīng)的分析

0風(fēng)輪掃掠面處不同葉片將風(fēng)中的動能轉(zhuǎn)化為機器能，這是風(fēng)動機的主要部件。其他部件的負荷主要是由葉片引起的。隨著風(fēng)力機向大型化的趨勢發(fā)展,風(fēng)力機塔架越來越高,葉片半徑越來越長,這種情況下風(fēng)剪切和塔影效應(yīng)對風(fēng)力機的影響也越發(fā)顯著。由于風(fēng)剪切和塔影效應(yīng)的存在,風(fēng)速值在整個風(fēng)輪掃掠面上處處不同,從而加劇了葉片在風(fēng)輪掃掠面上所受到的空氣動力載荷的周期性變化。作用在風(fēng)輪上的轉(zhuǎn)矩是風(fēng)輪軸功率的來源,它由葉片擺振方向載荷合成生成。由于作用于風(fēng)輪葉片上的周期性氣動載荷會引起葉片的動響應(yīng),而此響應(yīng)又反饋于外部氣動載荷,使得本就復(fù)雜的風(fēng)力機振動、疲勞、動力穩(wěn)定性等問題變得更加復(fù)雜且不容忽視。同時對并網(wǎng)型風(fēng)力發(fā)電機組而言,其輸出電壓和輸出功率也存在一定波動,輸出電能的品質(zhì)也會產(chǎn)生一定影響。目前,對于多兆瓦級風(fēng)力機來說,風(fēng)力機功率控制一般都是采用輪轂離地高度處風(fēng)速或平均風(fēng)速作為單一的計算風(fēng)速。這種控制方式隨風(fēng)力機容量的增大越發(fā)不合時宜。根據(jù)每個葉片所受載荷不同,對每個葉片進行獨立控制,這將是風(fēng)力機功率控制技術(shù)的一大發(fā)展方向。本文在考慮風(fēng)剪切和塔影效應(yīng)前提下,對大型三槳葉上風(fēng)向水平軸風(fēng)力機所受風(fēng)速進行了動態(tài)建模,并在此基礎(chǔ)上對一臺2MW風(fēng)力機空氣動力載荷和功率進行了詳細地計算和分析。1塔影效應(yīng)風(fēng)速模型大型三槳葉上風(fēng)向水平軸風(fēng)力機基本結(jié)構(gòu)如圖1。通常風(fēng)剪切是指風(fēng)速隨垂直高度的變化,塔影效應(yīng)是指風(fēng)遇到塔架堵塞改變大小和方向的變化,二者本質(zhì)上都是一個周期性變化過程。造成風(fēng)剪切現(xiàn)象的原因有動力因素和熱力因素,前者主要來源于地面的摩擦效應(yīng),即地面的粗糙度;后者主要表現(xiàn)為與近地層大氣垂直穩(wěn)定度的關(guān)系。風(fēng)剪切時風(fēng)速計算常用的是指數(shù)模型,如式(1):V(z)=Vh(zh)α(1)V(z)=Vh(zh)α(1)式中,V(z)——離地高度為z處的風(fēng)速;Vh——輪轂離地高度處風(fēng)速;h——輪轂離地高度;α——風(fēng)剪切系數(shù);z——離地垂直高度。式(1)還可寫成葉片微元半徑距離r和方位角θ的形式:W(r,θ)是風(fēng)剪切效應(yīng)對風(fēng)速施加的變化擾動,用三階泰勒函數(shù)展開有:W(r,θ)≈α(rΗ)cosθ+α(α-1)2(rΗ)2cos2θ+α(α-1)(α-2)6(rΗ)3cos3θ(3)W(r,θ)≈α(rH)cosθ+α(α?1)2(rH)2cos2θ+α(α?1)(α?2)6(rH)3cos3θ(3)由于塔架對氣流的堵塞,在塔架的上游和下游都將使來流速度減小,稱為塔影效應(yīng)。參考坐標(biāo)系如圖1,考慮塔影效應(yīng)時風(fēng)速的表達式為:V(y,x)=Vh+vt(y,x)(4)vt(y,x)=V0a2y2-x2(x2+y2)2(5)vt(y,x)=V0a2y2?x2(x2+y2)2(5)式中,V0——空間平均風(fēng)速;a——塔架半徑;y——槳葉微元到塔架軸線的y軸方向距離;x——槳葉微元到塔架軸線的x軸方向距離,即懸垂距離;vt(y,x)——塔影效應(yīng)對風(fēng)速施加的變化擾動。風(fēng)剪切和塔影效應(yīng)中使用的參考風(fēng)速不同,風(fēng)剪切使用的是輪轂離地高度處風(fēng)速Vh,塔影效應(yīng)中使用的是空間平均風(fēng)速V0,二者有如下關(guān)系式成立:V0=Vh[1+α(α-1)(R2)8h2]=mVh(6)V0=Vh[1+α(α?1)(R2)8h2]=mVh(6)式中,R——風(fēng)輪半徑;m=1+α(α-1)(R2)8h2m=1+α(α?1)(R2)8h2。實際中,Vh容易獲得,V0可通過式(6)計算獲得。由于y=rsinθ,式(5)可改寫成式(7):vt(r,θ,x)=mVha2r2sin2θ-x2(r2sin2θ+x2)2=Vhvtt(r,θ,x)(7)vt(r,θ,x)=mVha2r2sin2θ?x2(r2sin2θ+x2)2=Vhvtt(r,θ,x)(7)式中,vtt(r,θ,x)=ma2r2sin2θ-x2(r2sin2θ+x2)2vtt(r,θ,x)=ma2r2sin2θ?x2(r2sin2θ+x2)2,需要注意的是,塔影效應(yīng)只產(chǎn)生于下半風(fēng)輪掃掠面,即方位角0.5π≤θ≤1.5π。在上半風(fēng)輪掃掠面,即方位角為0≤θ<0.5π和1.5π<θ≤2π處,風(fēng)速計算模型為:V(r,θ)=Vh[1+W(r,θ)](8)在下半風(fēng)輪掃掠面,即方位角為0.5π≤θ≤1.5π處,風(fēng)速計算模型為:V(r,θ)=Vh[1+W(r,θ)][1+vtt(r,θ,x)](9)2機械轉(zhuǎn)速和風(fēng)輪軸功率的關(guān)系風(fēng)力機采用葉片坐標(biāo)系,根據(jù)動量-葉素理論,可求得作用在葉片上擺振方向(切向)的剪力fyb和彎矩Mxb、揮舞方向(法向)的剪力fxb和彎矩Myb分別為:fxb=12ρW2cCn(10)fyb=12ρW2cCt(11)Mxb=∫Rr0fybrdr(12)Myb=∫Rr0fxbrdr(13)葉片空氣動力扭轉(zhuǎn)力矩為:式中,R——風(fēng)輪半徑;下標(biāo)P——翼型壓力中心;下標(biāo)C——扭轉(zhuǎn)中心;c——葉素剖面弦長;r0——輪轂半徑;W——相對來流風(fēng)速。Cn=Clcosφ+Cdsinφ(15)Ct=Clsinφ-Cdcosφ(16)式中,Cn——法向力系數(shù);Ct——切向力系數(shù);Cl、Cd——翼型升力、阻力系數(shù);φ——入流角。作用在風(fēng)輪上的轉(zhuǎn)矩是風(fēng)輪軸功率的來源,它由葉片擺振力矩Mxb合成生成,與揮舞力矩Myb、扭轉(zhuǎn)力矩Mzb一樣隨葉片的方位角變化。在不考慮風(fēng)剪切和塔影效應(yīng)情況下,3個槳葉所受到的擺振力矩Mxb相同,此時風(fēng)力機機械轉(zhuǎn)矩如式(17):Tt=3Mxb(17)風(fēng)力機機械轉(zhuǎn)矩還可由式(18)近似表示:Τt=12ρπR3Cq(λ,β)V2(18)風(fēng)力機輸出功率與機械轉(zhuǎn)矩的關(guān)系為:Ρt=ΤtΩ=12ρπR2CΡ(λ,β)V3(19)式中,λ——風(fēng)力機葉尖速比,λ=ΩR/V;Ω——風(fēng)輪旋轉(zhuǎn)角速度;V——風(fēng)速;ρ——空氣密度;β——槳距角;CP(λ,β)、Cq(λ,β)——分別為風(fēng)能利用系數(shù)和轉(zhuǎn)矩系數(shù)。3風(fēng)速影響分析選取一2MW三槳葉上風(fēng)向水平軸風(fēng)力機,應(yīng)用Bladed和Matlab軟件進行功率控制和載荷分析。風(fēng)力機主要參數(shù)如表1。查翼型手冊得到作用在葉素上的翼型壓力中心、扭轉(zhuǎn)中心、升力系數(shù)Cl、阻力系數(shù)Cd。在考慮和不考慮風(fēng)剪切、塔影效應(yīng)兩種情況下,以槳葉A方位角為參量,分別對輪轂離地高度處風(fēng)速Vh取4個具有代表意義的數(shù)值,即6、12、18、25m/s,分別對2MW風(fēng)力機進行功率控制,額定風(fēng)速以下采用變速控制,額定風(fēng)速以上采用統(tǒng)一變槳控制,即3個槳葉接受同一槳距角改變命令。需要注意的是風(fēng)速取值6m/s時,采用變速功率控制方式,通過公式λ=ΩR/V計算,可知此時風(fēng)力機轉(zhuǎn)速為1.2rad/s,風(fēng)速取值12、18、25m/s時,采用變槳功率控制方式,風(fēng)力機轉(zhuǎn)速為額定轉(zhuǎn)速1.884rad/s?？紤]風(fēng)剪切、塔影效應(yīng)情況下,風(fēng)速變化曲線如圖2,兩種情況下風(fēng)力機擺振載荷、揮舞載荷和扭轉(zhuǎn)載荷響應(yīng)曲線如圖3～圖5。由圖2～圖5可得到:1)在上半風(fēng)輪掃掠面,槳葉可只考慮風(fēng)剪切的影響,在同一方位角的前提下,風(fēng)速隨風(fēng)輪半徑的增大而增大。方位角為0和2π時,即在槳葉葉尖處即掃掠面最高點處,風(fēng)速有最大值。在下半風(fēng)輪掃掠面,槳葉同時受到風(fēng)剪切和塔影效應(yīng)的雙重作用,風(fēng)速隨風(fēng)輪半徑的增大而減小,而且離風(fēng)輪中心越近,塔影效應(yīng)越嚴(yán)重,即風(fēng)速值趨向于變小,方位角為π時,即在槳葉葉尖處即掃掠面最低點處,風(fēng)速有最小值?？紤]風(fēng)輪掃掠面左右對稱,圖2只對槳葉在右半風(fēng)輪掃掠面,即方位角取值0≤θ≤π作了風(fēng)速變化分析。2)不考慮風(fēng)剪切及塔影效應(yīng)時,風(fēng)力機3個槳葉所受氣動載荷相同,當(dāng)考慮風(fēng)剪切及塔影效應(yīng)時,風(fēng)力機3個槳葉所受氣動載荷隨方位角呈周期性變化。如果風(fēng)輪旋轉(zhuǎn)頻率是1倍頻(1P),那么風(fēng)力機氣動載荷頻率就是3P。3)同一時刻3個槳葉,無論是擺振載荷還是揮舞載荷、扭轉(zhuǎn)載荷,數(shù)值相差較大。另外考慮和不考慮風(fēng)剪切、塔影效應(yīng)兩種情況下的載荷平均值也相差較大。這充分說明對于大型風(fēng)力機來說,風(fēng)剪切、塔影效應(yīng)對風(fēng)力機的影響不容忽視。要想對風(fēng)力機的輸出功率實現(xiàn)精確調(diào)節(jié),獨立變槳控制取代目前常用的統(tǒng)一變槳控制是必然趨勢。4)風(fēng)速在12m/s以上變化時,風(fēng)力機采用的是變槳功率控制方式,通過調(diào)整槳葉的槳距角,使得風(fēng)力機保持額定功率不變。由圖3b～圖3d可見,兩種情況下擺振載荷都基本維持在3.6×105Nm附近,由式(19)知風(fēng)力機基本保持額定功率2MW不變。由圖4b～圖4d可見,兩種情況下?lián)]舞載荷卻隨風(fēng)速的增大而大幅減小。5)由圖4可知,風(fēng)速在6m/s時的揮舞載荷值遠大于風(fēng)速在25m/s時的揮舞載荷值,可見即使在低風(fēng)速時,風(fēng)力機的氣動載荷對風(fēng)力機的振動、疲勞、動力穩(wěn)定性等也會產(chǎn)生很大影響。6)由圖5可知,由于存在變槳作用,風(fēng)力機扭轉(zhuǎn)載荷并不是線性變化的,在額定風(fēng)速處有負的最大值。4風(fēng)力機功率控制技術(shù)由于風(fēng)剪切和塔影效應(yīng)的存在,風(fēng)速在整個風(fēng)輪掃掠面上不是固定不變的。隨著風(fēng)力機向大型化的趨勢發(fā)展,塔架越來越高,風(fēng)輪直徑越來越大,風(fēng)剪切和塔影效應(yīng)對風(fēng)力機的影響也越來越大,使得風(fēng)力機每個葉片所受氣動載荷呈一個周期性的變化過程,同時葉片之間所受載荷相差也越來越大,對風(fēng)力機功率輸出、振動、疲勞、動力穩(wěn)定性等影響也越發(fā)不容忽視。顯然,對每個槳葉進行獨立控制取代對3個槳葉進行統(tǒng)一控制,將是風(fēng)力機功率控制