張立軍,馬東辰,胡闊亮,米玉霞,朱懷寶,顧嘉偉,劉靜,嚴(yán)強(qiáng)
(1.中國(guó)石油大學(xué)(華東)機(jī)電工程學(xué)院,山東青島,266580;2.上海麟風(fēng)風(fēng)能科技有限公司,上海,201508)
由于能源危機(jī)的日益嚴(yán)重,人們?cè)絹碓絻A向于對(duì)可再生能源進(jìn)行研究。其中,風(fēng)能以其綠色、無污染、可再生等特點(diǎn)已受到廣泛關(guān)注。我國(guó)風(fēng)能資源豐富,根據(jù)國(guó)家氣象局的資料,近海70 m高度的年平均風(fēng)功率密度可達(dá)300 W/m2以上,其中臺(tái)灣海峽和東海南部風(fēng)能資源甚為豐富,風(fēng)速大于6 m/s 時(shí)的累計(jì)小時(shí)數(shù)可達(dá)5 000 h[1]。目前對(duì)風(fēng)能利用的主要形式為風(fēng)力發(fā)電。根據(jù)其風(fēng)輪旋轉(zhuǎn)軸相對(duì)地面的安裝角度可以分為垂直軸風(fēng)力發(fā)電機(jī)(VAWT)和水平軸風(fēng)力發(fā)電機(jī)(HAWT),與水平軸風(fēng)力機(jī)相比,垂直軸風(fēng)力機(jī)有無需對(duì)風(fēng)裝置、易于安裝且便于維修等優(yōu)點(diǎn),尤其對(duì)于H 型垂直軸風(fēng)力機(jī)來說,其葉片采用等截面結(jié)構(gòu),更適用于采用計(jì)算流體力學(xué)方法得到更精確的結(jié)論;同時(shí),H型垂直軸風(fēng)力機(jī)噪聲小,比水平軸風(fēng)力機(jī)擁有更廣的應(yīng)用范圍,具有更高的商業(yè)使用價(jià)值。然而,現(xiàn)有的商業(yè)化H 型垂直軸風(fēng)力發(fā)電機(jī)的風(fēng)能利用率僅在30%~35%之間,遠(yuǎn)未達(dá)到其理論最大風(fēng)能利用率64%[2],且存在自啟動(dòng)能力弱等不足,在很大程度上限制了其規(guī)?;瘧?yīng)用。為解決垂直軸風(fēng)力機(jī)的低風(fēng)能利用率問題,人們研究了多種提高其性能的方法,這些方法主要?dú)w納為3種:一是風(fēng)力機(jī)的結(jié)構(gòu)優(yōu)化,包括葉片翼型和主軸結(jié)構(gòu)優(yōu)化、葉片傾角改變、升阻互補(bǔ)等;二是翼型區(qū)域流場(chǎng)控制,例如利用導(dǎo)流片產(chǎn)生干擾氣流,加裝葉頂端板、襟翼、圓臺(tái)型聚風(fēng)罩等;三是葉片攻角改變,即通過改變?nèi)~片槳距角來實(shí)現(xiàn)對(duì)攻角的改變。其中,改變?nèi)~片槳距角以其調(diào)節(jié)方式簡(jiǎn)便、高效、直接等優(yōu)點(diǎn)逐步受到人們的青睞。
現(xiàn)有的垂直軸風(fēng)力機(jī)旋轉(zhuǎn)過程中,葉片攻角隨方位角不斷變化,沒有維持在最佳攻角處,是其風(fēng)能利用率低的重要原因[3]。在槳距角調(diào)節(jié)過程中,主要角度是槳距角、攻角和方位角,這3個(gè)角度之間的關(guān)系如圖1所示。圖1中,W為誘導(dǎo)速度v和切向速度Rω的合成風(fēng)速,ω為風(fēng)輪的旋轉(zhuǎn)角速度;θ為葉片方位角,當(dāng)方位角θ位于0°~180°時(shí),該區(qū)域稱為風(fēng)輪的上風(fēng)區(qū);當(dāng)θ位于180°~360°時(shí),稱為風(fēng)輪的下風(fēng)區(qū);α為葉片攻角,是合成風(fēng)速方向與葉片弦線方向之間的夾角;β為葉片槳距角,是切向速度方向與葉片弦線方向之間的夾角。分析圖1中的速度關(guān)系,可以得出垂直軸風(fēng)力機(jī)葉片攻角的表達(dá)式為
式中:λ為風(fēng)輪的葉尖速比,
圖1 翼型角度關(guān)系圖Fig.1 Airfoil angle diagram
變槳技術(shù)的主要思路是通過調(diào)節(jié)槳距角使葉片攻角始終在重點(diǎn)旋轉(zhuǎn)區(qū)域內(nèi)處于最佳攻角處,即在一定葉尖速比下切向力系數(shù)達(dá)到最大時(shí)所對(duì)應(yīng)的攻角,切向力系數(shù)與攻角的關(guān)系如圖2所示。圖2中,CL和CD分別代表翼型的升力系數(shù)和阻力系數(shù),其值可由NACA 系列翼型空氣動(dòng)力學(xué)特性數(shù)據(jù)庫查得。CT為切向力系數(shù),是CL和CD在弦長(zhǎng)方向的合成值。由圖2中的幾何關(guān)系,得到CT的計(jì)算式為
圖2 葉片翼型氣動(dòng)力系數(shù)Fig.2 Aerodynamic coefficient of blade profile
隨著風(fēng)輪的旋轉(zhuǎn),葉片的合成風(fēng)速方向時(shí)刻變化,攻角也隨之變化,不易直接進(jìn)行調(diào)節(jié),而槳距角可通過調(diào)節(jié)葉片與旋轉(zhuǎn)切向方向之間的夾角來調(diào)節(jié)。同時(shí),由于槳距角與攻角之間有關(guān),槳距角的改變也會(huì)使攻角發(fā)生變化,因此,變槳技術(shù)的核心是通過合理調(diào)節(jié)槳距角來間接控制對(duì)應(yīng)攻角的變化,以實(shí)現(xiàn)對(duì)風(fēng)力機(jī)氣動(dòng)性能的改善。
垂直軸風(fēng)力機(jī)變槳距的研究重點(diǎn)在于變槳規(guī)律的可行性分析,目前國(guó)內(nèi)外學(xué)者運(yùn)用的主要研究方法有動(dòng)量模型法、數(shù)值計(jì)算法和測(cè)試試驗(yàn)法等。
動(dòng)量模型法利用流體力學(xué)(含空氣動(dòng)力學(xué))的相關(guān)理論如動(dòng)量理論、葉素理論、流管理論,借助水平軸風(fēng)力機(jī)的分析方法,通過提出假定,建立模型的方法來計(jì)算。隨著垂直軸風(fēng)力機(jī)氣動(dòng)特性研究及空氣動(dòng)力學(xué)的不斷完善,上述理論在風(fēng)力機(jī)的設(shè)計(jì)和結(jié)構(gòu)改進(jìn)過程中起到了越來越重要的作用。其中,對(duì)于變槳距垂直軸風(fēng)力機(jī)的理論研究,大部分學(xué)者運(yùn)用的主要是葉素理論和雙制動(dòng)盤多流管理論。
2.1.1 基于葉素理論分析
垂直軸風(fēng)力機(jī)在不同葉尖速比下運(yùn)行時(shí),葉片攻角的變化范圍不同。應(yīng)用葉素理論便于得到翼型在不同葉尖速比下的攻角變化范圍,因此,在研究特定葉尖速比下風(fēng)力機(jī)變槳規(guī)律時(shí),通常選用葉素理論來進(jìn)行分析。
顧華朋等[4]基于葉素理論,得到了NACA0012翼型在不同葉尖速比下攻角隨方位角的變化關(guān)系,如圖3所示。通過進(jìn)一步結(jié)合翼型的升阻力系數(shù)隨攻角變化關(guān)系,提出了一種以葉尖速比為分區(qū)的變槳方案,翼型的升阻力系數(shù)曲線如圖4所示,其變槳方案如式(4)所示。該變槳規(guī)律的合理性主要在于:在低葉尖速比、大攻角工況下,通過調(diào)節(jié)槳距角使葉片攻角在60°~140°之間,此時(shí)垂直軸風(fēng)力機(jī)主要靠葉片的阻力做功獲得較大啟動(dòng)力矩;在較高葉尖速比下,葉片攻角變化范圍較大,利用變槳方法控制葉片攻角在失速點(diǎn)±15°附近,使得葉片可以依靠升力做功從而獲得較大且穩(wěn)定的扭矩。通過建立風(fēng)輪扭矩系數(shù)模型,分析了該變槳運(yùn)行規(guī)律的有效性。
吳祥輝等[5]基于葉素理論,進(jìn)一步對(duì)葉尖速比λ<1.0 情況下的變槳規(guī)律進(jìn)行了討論,主要是通過求解不同的方位角和槳距角時(shí)所對(duì)應(yīng)的扭矩力,比較了各個(gè)方位角下最大扭矩力對(duì)應(yīng)的最佳槳距角,進(jìn)而得到1周槳距角的調(diào)節(jié)規(guī)律。得到的槳距角變化規(guī)律如圖5所示。
圖3 定槳葉片攻角變化曲線[4]Fig.3 Change curve of attack angle with fixed pitch angle[4]
圖4 NACA0012翼型大攻角升阻力系數(shù)曲線[4]Fig.4 NACA0012 airfoil lift drag coefficient curve of high attack angle[4]
2.1.2 基于雙制動(dòng)盤多流管理論分析
與葉素理論不同,雙制動(dòng)盤多流管理論不僅以單一葉片為分析對(duì)象,它還將風(fēng)輪旋轉(zhuǎn)分為上下風(fēng)區(qū),并將轉(zhuǎn)子作用盤面沿垂直于來流的方向細(xì)分成多個(gè)獨(dú)立微流管[6],理論上看對(duì)于垂直軸風(fēng)力機(jī)的氣動(dòng)性能分析更為合理。
圖5 槳距角變化規(guī)律[5]Fig.5 Changing rule of pitch angle[5]
PARASCHIVOIU 等[7]針對(duì)一種7 kW 的垂直軸風(fēng)力機(jī),將CARDAAV 代碼與基于遺傳算法的優(yōu)化器相結(jié)合,開發(fā)了一種新的優(yōu)化工具來定義葉片槳距角的變化。以基于雙制動(dòng)盤多流管理論求解的氣動(dòng)參數(shù)為評(píng)價(jià)指標(biāo),驗(yàn)證了優(yōu)化算法的可行性。
張立勛等[8]針對(duì)變槳距垂直軸風(fēng)力機(jī)憑某一尖速比下最優(yōu)變槳規(guī)律無法啟動(dòng)和變槳規(guī)律中的不連續(xù)性問題,重點(diǎn)對(duì)此不連續(xù)部分進(jìn)行了修正,如圖6所示,通過人工設(shè)計(jì)一條斜率最小的直線段,連接變槳規(guī)律中不連續(xù)段的兩端,并修改方位角n×180°附近的葉片槳距角為90°,實(shí)現(xiàn)規(guī)律連續(xù)性。采用雙制動(dòng)盤多流管理論對(duì)所設(shè)計(jì)的變槳規(guī)律進(jìn)行力矩特性分析,結(jié)果顯示,采用該變槳規(guī)律的風(fēng)力機(jī)能夠可靠自啟動(dòng)且高效發(fā)電。
張立軍[9]針對(duì)垂直軸風(fēng)力機(jī)下風(fēng)區(qū)風(fēng)況的復(fù)雜性,提出了局部葉尖速比的概念,基于雙制動(dòng)盤多流管理論,重新建立了下風(fēng)區(qū)葉片攻角與槳距角之間的關(guān)系式,對(duì)雙制動(dòng)盤多流管理論進(jìn)行了改進(jìn)和完善。
綜合上述研究發(fā)現(xiàn),采用動(dòng)量模型法對(duì)變槳距垂直軸風(fēng)力機(jī)進(jìn)行研究的目的主要有2個(gè):一是通過建立理論模型,得到槳距角調(diào)節(jié)規(guī)律;二是以風(fēng)力機(jī)的風(fēng)能利用率為評(píng)價(jià)指標(biāo),判定變槳前后風(fēng)力機(jī)的性能變化,驗(yàn)證變槳規(guī)律的可行性。這種研究方法簡(jiǎn)單易行,可以直接反映風(fēng)力機(jī)的出力效果,便于對(duì)比變槳前后風(fēng)力機(jī)的氣動(dòng)性能,但不能直觀觀察風(fēng)力機(jī)流場(chǎng)的變化,對(duì)于功率系數(shù)變化的原因也缺少探討。
圖6 修正前后變槳方案對(duì)比[8]Fig.6 Schemes comparison of before-and-after variable pitch modification[8]
近幾年來,隨著計(jì)算流體力學(xué)(CFD)技術(shù)的快速發(fā)展,數(shù)值計(jì)算方法以其精度高、計(jì)算速度快、成本低等優(yōu)點(diǎn),成為垂直軸風(fēng)力機(jī)研究的主要手段之一[10]。對(duì)變槳距垂直軸風(fēng)力機(jī)的研究通常采用二維計(jì)算方法,主要原因在于變槳技術(shù)并未改變翼型形狀,研究重點(diǎn)是變槳前后風(fēng)力機(jī)氣動(dòng)性能參數(shù)的改變,相對(duì)于三維計(jì)算,二維計(jì)算網(wǎng)格設(shè)置相對(duì)簡(jiǎn)單,工作量小,計(jì)算效率高[11]。一種經(jīng)典的垂直軸風(fēng)力機(jī)二維仿真模型如圖7所示,圖中d為風(fēng)輪直徑。
針對(duì)垂直軸風(fēng)力機(jī)的變槳問題,數(shù)值計(jì)算法的一個(gè)重要優(yōu)勢(shì)在于通過改變邊界條件,便于研究不同槳距角下風(fēng)力機(jī)的氣動(dòng)性能差異,模擬不同葉尖速比下風(fēng)力機(jī)所處工況條件。張立勛等[12]仿真研究了不同槳距角下風(fēng)力機(jī)的氣動(dòng)性能,通過分析不同槳距角下力矩系數(shù),如圖8所示,提出一種主動(dòng)式變槳規(guī)律:設(shè)置內(nèi)偏置槳距角或小幅值外偏置槳距角,使風(fēng)力機(jī)在小方位角下進(jìn)入有效力矩作用區(qū);設(shè)置外偏置槳距角使風(fēng)力機(jī)延緩?fù)顺鲇行Я刈饔脜^(qū)。研究表明,這種方法增加了有效切向力的作用范圍。
對(duì)于同一類問題,ABDALRAHMAN等[13]利用FLUENT 流體仿真軟件對(duì)不同葉尖速比下槳距角進(jìn)行研究,通過分析同一葉尖速比下槳距角在不同方位角下的出力情況,得到了風(fēng)輪旋轉(zhuǎn)1 周時(shí),各個(gè)葉片對(duì)應(yīng)的槳距角變化規(guī)律,如圖9所示。
圖7 垂直軸風(fēng)力機(jī)二維仿真示意圖Fig.7 Two-dimensional simulation diagram of VAWT
圖8 不同槳距角下力矩系數(shù)圖[12]Fig.8 Torque coefficients at different pitch angles[12]
除此之外,數(shù)值計(jì)算法在觀察流場(chǎng)分布上有其獨(dú)特優(yōu)勢(shì)。垂直軸風(fēng)力機(jī)風(fēng)能利用率較低的重要原因在于攻角的不斷變化,而攻角改變所導(dǎo)致的葉片周圍流場(chǎng)的變化是翼型升阻力系數(shù)變化的根本原因。借助CFD 仿真軟件,能夠直觀觀察葉片在各個(gè)方位角下的流場(chǎng)情況,對(duì)于研究變槳式風(fēng)力機(jī)有較好的輔助作用。左薇等[14]采用笛卡兒動(dòng)網(wǎng)格數(shù)值計(jì)算方法研究了動(dòng)態(tài)變槳式H 型風(fēng)力機(jī)的氣動(dòng)特性,重點(diǎn)分析了變槳前后不同時(shí)刻葉片周圍的渦量云圖,研究發(fā)現(xiàn)動(dòng)態(tài)變槳方式不會(huì)生成前緣分離渦,且可以避免葉片脫落尺度較大、強(qiáng)度較強(qiáng)的尾緣分離渦。
圖9 槳距角變化圖(葉尖速比λ=1.0)[13]Fig.9 Pitch angle variation diagram(tip speed ratioλ=1.0)[13]
彭偉等[15]提出一種新型隨動(dòng)變槳距垂直軸風(fēng)力機(jī)結(jié)構(gòu),如圖10所示,這種風(fēng)力機(jī)的主要特點(diǎn)在于:當(dāng)葉片處于迎風(fēng)面時(shí),葉片始終與來流風(fēng)垂直;當(dāng)葉片處于背風(fēng)面時(shí),葉片始終與來流風(fēng)平行,以此來減少阻力,增加其產(chǎn)生的力矩。通過對(duì)該變槳距風(fēng)力機(jī)進(jìn)行CFD壓力場(chǎng)仿真分析可知,該風(fēng)力機(jī)在轉(zhuǎn)動(dòng)過程中可以使葉片產(chǎn)生的壓力相對(duì)穩(wěn)定,有助于風(fēng)力機(jī)保持勻速轉(zhuǎn)動(dòng)。
圖10 新型隨動(dòng)變槳垂直軸風(fēng)力機(jī)示意圖[15]Fig.10 Schematic diagram of a new type of follow drive variable pitch vertical axis wind turbine[15]
REZAEIHA 等[16]指出,當(dāng)前對(duì)風(fēng)力機(jī)變槳規(guī)律的研究局限在于以穩(wěn)態(tài)值(如風(fēng)能利用系數(shù)等)為指標(biāo),而忽視了對(duì)瞬態(tài)值(如風(fēng)力機(jī)旋轉(zhuǎn)過程中的實(shí)時(shí)負(fù)載及力矩變化、邊界層分離現(xiàn)象及層流與湍流的轉(zhuǎn)化情況等)的研究。為此,通過建立風(fēng)輪CFD 仿真模型,對(duì)不同槳距角下風(fēng)力機(jī)的扭矩系數(shù)、葉片壓力分布、旋渦脫落等情況進(jìn)行了研究,綜合考慮瞬態(tài)值與穩(wěn)態(tài)值的變化,給出了最佳槳距角為β=-2°。其中,扭矩系數(shù)的變化趨勢(shì)如圖11所示[16]。
圖11 不同槳距角下風(fēng)力機(jī)扭矩系數(shù)圖[16]Fig.11 Torque coefficient diagram of wind turbine at different pitch angles[16]
上述研究均采用數(shù)值計(jì)算法,充分發(fā)揮了數(shù)值計(jì)算法精確度高、對(duì)仿真結(jié)果后處理能力強(qiáng)、流程顯示直觀等優(yōu)點(diǎn)。可以預(yù)見的是,隨著計(jì)算機(jī)算法的不斷完善,CFD 仿真計(jì)算的可靠性將越來越高,會(huì)為學(xué)者們的研究帶來更大便利。
數(shù)值計(jì)算法雖然便于求得風(fēng)力機(jī)的氣動(dòng)參數(shù),且易于觀察風(fēng)力機(jī)周圍流場(chǎng)的變化,但是由于所選計(jì)算模型、條件設(shè)置以及網(wǎng)格質(zhì)量等因素,計(jì)算仍存在一定的誤差,而測(cè)試試驗(yàn)法則更加直觀,并且試驗(yàn)得到的數(shù)據(jù)也更具有說服力。
目前,測(cè)試試驗(yàn)法主要分為2種。
一種是將風(fēng)力機(jī)樣機(jī)置于實(shí)際風(fēng)場(chǎng)中,直接利用傳感器等測(cè)試儀器和設(shè)備監(jiān)測(cè)風(fēng)力機(jī)的性能。例如LEE等[17]為驗(yàn)證基于數(shù)值計(jì)算得到的擺線式變槳規(guī)律[18]的正確性,搭建了風(fēng)力機(jī)裝置并進(jìn)行野外試驗(yàn),同時(shí)利用LabVIEW 軟件平臺(tái)對(duì)風(fēng)力機(jī)實(shí)時(shí)功率進(jìn)行采集,結(jié)果如圖12所示[17],對(duì)比發(fā)現(xiàn)試驗(yàn)結(jié)果低于CFD 仿真結(jié)果,但2 種結(jié)果的趨勢(shì)大致相同。這種研究方法的優(yōu)點(diǎn)在于風(fēng)力機(jī)位于實(shí)際流場(chǎng)中,對(duì)工況的還原度高,但存在的主要缺點(diǎn)是不能人為地改變風(fēng)速,對(duì)自然環(huán)境的依賴程度較強(qiáng)。同時(shí),由于需要搭建實(shí)際的風(fēng)力機(jī)測(cè)試系統(tǒng),研發(fā)成本較高,而且基于風(fēng)速的隨機(jī)性,受地勢(shì)、周圍環(huán)境等因素的影響,有時(shí)候測(cè)試結(jié)果未必真實(shí)反映風(fēng)力機(jī)的性能。
圖12 試驗(yàn)結(jié)果與仿真結(jié)果的對(duì)比[17]Fig.12 Comparison of experimental results and simulation results[17]
另一種測(cè)試試驗(yàn)方法是對(duì)變槳距風(fēng)力機(jī)在風(fēng)洞內(nèi)進(jìn)行性能測(cè)試,如圖13所示,通過調(diào)節(jié)軸流風(fēng)機(jī)的轉(zhuǎn)動(dòng)頻率,可以實(shí)現(xiàn)對(duì)風(fēng)速可控調(diào)節(jié),監(jiān)測(cè)不同工況下風(fēng)力機(jī)的性能[19]。但這種方法需要制造標(biāo)準(zhǔn)化的風(fēng)洞,研發(fā)成本也比較高。
圖13 垂直軸風(fēng)力機(jī)風(fēng)洞試驗(yàn)系統(tǒng)[19]Fig.13 Wind tunnel experimental system for VAWT[19]
張立軍等[20]研發(fā)了一種小型變槳距垂直軸風(fēng)力機(jī)樣機(jī),如圖14所示,利用自制風(fēng)洞,模擬了垂直軸風(fēng)力機(jī)的流場(chǎng)情況,對(duì)3~10 m/s 風(fēng)速下的風(fēng)力機(jī)的發(fā)電功率進(jìn)行了監(jiān)測(cè),試驗(yàn)結(jié)果如表1所示,結(jié)果顯示,變槳距風(fēng)力機(jī)的發(fā)電效率較定槳距風(fēng)力機(jī)的至少能提高7.86%。
在變槳規(guī)律研究中,測(cè)試試驗(yàn)法的主要優(yōu)勢(shì)是可信度大,但也存在制造成本高、研究周期長(zhǎng)、占地面積大等不足,成為阻礙其發(fā)展的主要原因。但隨著新能源技術(shù)越來越受到重視、國(guó)家政策的支持及相關(guān)產(chǎn)業(yè)鏈的完善,變槳距垂直軸風(fēng)力機(jī)的測(cè)試試驗(yàn)也將迎來更好的發(fā)展機(jī)遇。
圖14 變槳距垂直軸風(fēng)力機(jī)[20]Fig.14 Variable pitch vertical axis wind turbine[20]
表1 2種風(fēng)力機(jī)發(fā)電功率的對(duì)比Table 1 Comparison of power generated by two types of wind turbines
對(duì)于變槳距垂直軸風(fēng)力機(jī)而言,無論是采用動(dòng)量模型,還是采用數(shù)值計(jì)算以及測(cè)試試驗(yàn),最終目的在于實(shí)現(xiàn)其商業(yè)化應(yīng)用。而變槳距風(fēng)力機(jī)制作的關(guān)鍵在于使其葉片槳距角的變化與所求得的變槳規(guī)律相一致,為此,國(guó)內(nèi)外學(xué)者對(duì)此也進(jìn)行了相關(guān)探索和研究,目前實(shí)現(xiàn)變槳的方法主要分為2 種:一是機(jī)械傳動(dòng)變槳,二是液壓傳動(dòng)變槳。
采用機(jī)械傳動(dòng)機(jī)構(gòu)是實(shí)現(xiàn)葉片槳距角調(diào)節(jié)的有效途徑。2014年,楊凱等[21]利用伺服電機(jī)的精確性和蝸輪蝸桿減速機(jī)的自鎖功能,設(shè)計(jì)了一種自主變槳式風(fēng)力機(jī),如圖15所示。其主要工作原理是:利用光電編碼器檢測(cè)風(fēng)力機(jī)的轉(zhuǎn)速及葉片位置,通過變槳伺服電機(jī)調(diào)節(jié)葉片槳距角,并利用蝸輪蝸桿減速機(jī)的反向自鎖功能,保持槳距角在特定區(qū)域內(nèi)保持不變。
圖15 變槳風(fēng)力機(jī)剖面圖[21]Fig.15 Cross view of pitching wind turbine[21]
上述方案的局限性之一在于蝸輪蝸桿價(jià)格較為昂貴,制作成本較高。近幾年,許多學(xué)者提出了利用連桿機(jī)構(gòu)實(shí)現(xiàn)對(duì)槳距角的調(diào)節(jié)。YAMADA等[22]提出采用四連桿機(jī)構(gòu)來控制風(fēng)力機(jī)槳距角的垂直軸風(fēng)力機(jī)模型。JAIN 等[23]也提出類似的四連桿機(jī)構(gòu)模型,如圖16所示。
圖16 連桿式變槳距風(fēng)力機(jī)模型[23]Fig.16 Model of variable pitch VAWT with connecting rod[23]
廉正光等[24-25]提出一種采用雙曲柄調(diào)距機(jī)構(gòu)的垂直軸風(fēng)力機(jī),如圖17所示,該風(fēng)力機(jī)主要采用雙曲柄調(diào)距機(jī)構(gòu)及雙偏心軸機(jī)構(gòu)。
圖17 基于調(diào)距機(jī)構(gòu)的整體風(fēng)機(jī)模型[24]Fig.17 Overall structure of distance adjusting mechanism[24]
對(duì)于垂直軸風(fēng)力機(jī)變槳機(jī)構(gòu)中連桿機(jī)構(gòu)的設(shè)計(jì)重點(diǎn)在于各段桿長(zhǎng)的確定。張立軍等[20]采用隨機(jī)方向法對(duì)一種雙曲柄機(jī)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì),確定了各桿的桿長(zhǎng),得到了不同葉尖速比下變槳機(jī)構(gòu)的參數(shù)組合。雙曲柄機(jī)構(gòu)的幾何分析如圖18所示。
圖18 雙曲柄機(jī)構(gòu)幾何分析圖[20]Fig.18 Analysis diagram of double crank mechanism[20]
采用機(jī)械傳動(dòng)方式實(shí)現(xiàn)周期性槳距角調(diào)整的優(yōu)點(diǎn)在于結(jié)構(gòu)緊湊、傳動(dòng)平穩(wěn)、便于維護(hù)等,但是并不利于風(fēng)力機(jī)的大型化發(fā)展。
針對(duì)兆瓦級(jí)H 型垂直軸風(fēng)力機(jī),變槳機(jī)構(gòu)的設(shè)計(jì)需滿足驅(qū)動(dòng)力大、有足夠的強(qiáng)度和精度等要求[26],采用液壓傳動(dòng)變槳系統(tǒng),可以起到機(jī)構(gòu)響應(yīng)速度快、定位準(zhǔn)確、驅(qū)動(dòng)力大等優(yōu)點(diǎn)。如圖19所示,液壓變槳系統(tǒng)主要由變槳控制器、D/A 轉(zhuǎn)化器、電液比例閥和變槳執(zhí)行機(jī)構(gòu)等組成[27]。其工作原理是:通過電液比例閥控制液壓缸活塞桿驅(qū)動(dòng)曲柄連桿機(jī)構(gòu)實(shí)現(xiàn)控制變槳機(jī)構(gòu)的槳距角,通過位移傳感器,將活塞桿的位移反饋給PID槳距控制器,PID槳距控制器根據(jù)槳距角的給定值和實(shí)際響應(yīng)的測(cè)量值進(jìn)行比較形成偏差,進(jìn)而進(jìn)行參數(shù)放大,計(jì)算出電液比例閥的控制電壓。在此基礎(chǔ)上,再通過D/A 轉(zhuǎn)換器控制電液比例閥的輸出流量,從而驅(qū)動(dòng)液壓缸活塞桿的運(yùn)動(dòng),實(shí)現(xiàn)槳距角的大小調(diào)節(jié)[28]。
圖19 液壓變槳系統(tǒng)控制圖[27]Fig.19 Control diagram of hydraulic pitching system[27]
除可用于變槳外,液壓技術(shù)還可實(shí)現(xiàn)垂直軸風(fēng)力機(jī)的穩(wěn)壓恒頻發(fā)電。為了使垂直軸風(fēng)力機(jī)吸收更多的風(fēng)能,降低研發(fā)成本,張立軍等[29-30]提出了一種塔架式多層多軸垂直軸風(fēng)力機(jī),其結(jié)構(gòu)如圖20所示?!岸鄬印笨梢允癸L(fēng)力機(jī)充分吸收不同高度上的風(fēng)能;“多軸”可以降低風(fēng)剪效應(yīng)[31-32]對(duì)單一旋轉(zhuǎn)軸的制約,減小其旋轉(zhuǎn)軸在上下位置處的扭矩差。同時(shí),該風(fēng)力機(jī)采用液壓傳動(dòng),不但省去了機(jī)械式變速箱和整流逆變器,并用普通交流電動(dòng)機(jī)代替了永磁發(fā)電機(jī),制造成本降低了約30%,而且占地面積?。徊捎醚a(bǔ)償+穩(wěn)能技術(shù),系統(tǒng)能輸出(50±0.2) Hz 的恒頻率的交流電。目前該風(fēng)力發(fā)電機(jī)正處于初試階段,初步測(cè)試結(jié)果表明,基本能夠滿足并網(wǎng)要求。
圖20 多層多軸塔架式垂直軸風(fēng)力機(jī)[30]Fig.20 VAWT with multi-layer multi-axis tower[30]
現(xiàn)階段,雖然各國(guó)學(xué)者對(duì)于變槳距H 型垂直軸風(fēng)力機(jī)做了大量研究,且對(duì)于變槳規(guī)律的實(shí)現(xiàn)也提出了多種設(shè)計(jì)方案[33-35]。但目前除了用于風(fēng)速高于切出風(fēng)速時(shí)的“限速”外,垂直軸風(fēng)力機(jī)用于1周連續(xù)改變攻角進(jìn)而提高風(fēng)能利用率的變槳機(jī)構(gòu)并未得到商業(yè)化應(yīng)用,究其原因主要是在變槳機(jī)構(gòu)響應(yīng)時(shí)間、1周實(shí)時(shí)變槳規(guī)律計(jì)算、下風(fēng)區(qū)流場(chǎng)復(fù)雜性分析和動(dòng)態(tài)失速下可靠變槳等方面存在理論和技術(shù)瓶頸,機(jī)構(gòu)實(shí)現(xiàn)起來非常困難,尤其是對(duì)中小型垂直軸風(fēng)力機(jī)。
圖21 變槳系統(tǒng)控制流程圖Fig.21 Control flow chart of variable pitch system
目前的變槳系統(tǒng)控制流程圖如圖21所示,變槳結(jié)構(gòu)復(fù)雜,主要包括感應(yīng)機(jī)構(gòu)、調(diào)節(jié)機(jī)構(gòu)、機(jī)械與執(zhí)行機(jī)構(gòu)等多個(gè)環(huán)節(jié),而且存在閉環(huán)控制,很難實(shí)現(xiàn)機(jī)構(gòu)的快速響應(yīng)。特別是對(duì)于中小型垂直軸風(fēng)力機(jī)而言,風(fēng)力機(jī)的額定轉(zhuǎn)速高,對(duì)應(yīng)葉片攻角的變化頻率較快。以上海麟風(fēng)風(fēng)能科技有限公司所生產(chǎn)的小型垂直軸風(fēng)力機(jī)為例,風(fēng)力機(jī)的相關(guān)參數(shù)如表2所示。由表2可知:這些小型風(fēng)力機(jī)方位角的變化頻率較高,如額定功率為3 kW的風(fēng)力機(jī),在額定風(fēng)速和額定轉(zhuǎn)速下,方位角每變化10°,僅需要0.015 s[36]。由于風(fēng)速的隨機(jī)變化、葉片的慣性、連接機(jī)構(gòu)間的間隙、外部環(huán)境等因素的影響都會(huì)限制整個(gè)變槳機(jī)構(gòu)在實(shí)際運(yùn)行中的反應(yīng)速度,因此,變槳系統(tǒng)很難在較短時(shí)間內(nèi)適應(yīng)快速的風(fēng)速變化而完成葉片攻角的改變。
表2 小型垂直軸風(fēng)力機(jī)相關(guān)參數(shù)Table 2 Related parameters of small vertical axis wind turbine
在額定工況下,上述麟風(fēng)3 kW 垂直軸風(fēng)力機(jī)的葉尖速比λ=0.86,在該葉尖速比下,葉片攻角的1 周變化規(guī)律如圖22所示。由圖22可知:當(dāng)方位角變化10°時(shí),在大部分方位角下攻角變化5°左右,且在方位角θ=180°附近,攻角發(fā)生跳變。因此,在極短時(shí)間內(nèi),葉片攻角變化范圍較大。
而在實(shí)際工作狀態(tài)下,通常來流風(fēng)的風(fēng)速和風(fēng)向都在實(shí)時(shí)變化。風(fēng)速的改變帶來風(fēng)力機(jī)葉尖速比的變化,而在不同葉尖速比下,風(fēng)力機(jī)的攻角變化規(guī)律和槳距角調(diào)節(jié)規(guī)律不同;風(fēng)向的改變帶來葉片所處方位角的變化,而1周變槳規(guī)律強(qiáng)調(diào)在每個(gè)方位角處,葉片都對(duì)應(yīng)1個(gè)槳距角理論調(diào)節(jié)值,在這個(gè)理論調(diào)節(jié)值下對(duì)應(yīng)最佳攻角。因此,在實(shí)際工況下,為維持風(fēng)力機(jī)較高的輸出功率,每個(gè)時(shí)刻對(duì)葉片轉(zhuǎn)動(dòng)角度的調(diào)節(jié)值都不盡相同。而如何快速識(shí)別風(fēng)速隨機(jī)變化,確定葉片實(shí)時(shí)所處方位角,然后再計(jì)算槳距角調(diào)節(jié)值,已成為變槳機(jī)構(gòu)實(shí)現(xiàn)的難點(diǎn)。
圖22 葉尖速比λ=0.86時(shí)攻角變化規(guī)律Fig.22 Rule of attack angle with blade tip velocity ratioλ=0.86
與水平軸風(fēng)力機(jī)相比,垂直軸風(fēng)力機(jī)流場(chǎng)結(jié)構(gòu)更復(fù)雜,屬于典型的大分離非定常流動(dòng),且葉片在上風(fēng)區(qū)和下風(fēng)區(qū)要2次受到來流風(fēng)的作用。當(dāng)氣流經(jīng)過上風(fēng)區(qū)葉片到達(dá)下風(fēng)區(qū)時(shí),出現(xiàn)了擾流現(xiàn)象,此時(shí)下風(fēng)區(qū)各個(gè)方位誘導(dǎo)速度發(fā)生了改變,且方向也與來流風(fēng)方向不同。由于下風(fēng)區(qū)流場(chǎng)分布復(fù)雜,雙致動(dòng)盤多流管理論已經(jīng)很難完全反映誘導(dǎo)速度方向的變化,故在對(duì)垂直軸風(fēng)力機(jī)1周變槳規(guī)律制定時(shí),下風(fēng)區(qū)槳距角調(diào)節(jié)規(guī)律很難準(zhǔn)確得到。而目前關(guān)于垂直軸風(fēng)力機(jī)的分析,大都假設(shè)了垂直軸風(fēng)輪下風(fēng)區(qū)各個(gè)方位的誘導(dǎo)速度方向不變,這樣的假設(shè)雖簡(jiǎn)化了計(jì)算,但研究結(jié)果不能完全真實(shí)反映實(shí)際情況。
本文作者將動(dòng)量模型和數(shù)值計(jì)算技術(shù)相結(jié)合,在下風(fēng)區(qū)提出了局部葉尖速比的概念,并重新建立了下風(fēng)區(qū)葉片攻角的計(jì)算模型[9],如圖23所示。
定義下風(fēng)區(qū)局部葉尖速比X’為
圖23 下風(fēng)區(qū)葉片攻角分析模型Fig.23 Blade angle of attack analysis model in downwind area
根據(jù)圖23所示的幾何關(guān)系,得到下風(fēng)區(qū)葉片攻角的計(jì)算公式為
利用數(shù)值計(jì)算對(duì)推導(dǎo)出的垂直軸風(fēng)輪下風(fēng)區(qū)攻角的計(jì)算公式進(jìn)行了驗(yàn)證,在此基礎(chǔ)上制定了下風(fēng)區(qū)攻角調(diào)節(jié)策略。需要指出的是,這種方法是本文作者的一種合理嘗試,由于下風(fēng)區(qū)流場(chǎng)的復(fù)雜性,目前關(guān)于下風(fēng)區(qū)槳距角調(diào)節(jié)策略的研究還尚不成熟,或許需要進(jìn)一步建立針對(duì)下風(fēng)區(qū)流場(chǎng)精確計(jì)算的新方法。
4.4.1 風(fēng)力機(jī)運(yùn)行中動(dòng)態(tài)失速現(xiàn)象的發(fā)生情況
動(dòng)態(tài)失速是指在進(jìn)口來流攻角快速變化的過程中,風(fēng)輪葉片所表現(xiàn)出的與靜態(tài)風(fēng)洞實(shí)驗(yàn)完全不同的氣動(dòng)特性[37],是垂直軸風(fēng)力發(fā)電機(jī)低尖速比下運(yùn)行時(shí)的必然結(jié)果[38]。NACA0012翼型在動(dòng)態(tài)失速下上下風(fēng)區(qū)的失速攻角為17°和-17°[39]。圖24所示為在葉尖速比為0.5,1.0 和2.0 時(shí)該翼型攻角隨方位角的變化趨勢(shì)以及在這些葉尖速比下風(fēng)力機(jī)旋轉(zhuǎn)1周時(shí)失速區(qū)的位置和大小。可見葉片主要工作在失速區(qū),且葉尖速比越低,失速區(qū)越大,這與相關(guān)研究結(jié)論是一致的[37]。由圖24可知:λ=0.5時(shí)失速區(qū)占1周方位角的比例約為80%,此時(shí)葉片周圍出現(xiàn)分離流場(chǎng),受到的氣動(dòng)阻力急劇上升,攻角調(diào)節(jié)更加復(fù)雜。另外,現(xiàn)有的葉片動(dòng)態(tài)氣動(dòng)性能相關(guān)實(shí)驗(yàn)數(shù)據(jù)較少,這也極大地制約了垂直軸風(fēng)力機(jī)變槳規(guī)律的突破性研究。
圖24 3種葉尖速比下葉片攻角變化規(guī)律及失速區(qū)Fig.24 Variation of blade angle of attack and stall zone under three tip velocity ratios
4.4.2 考慮動(dòng)態(tài)失速下實(shí)時(shí)變槳規(guī)律
本文作者基于美國(guó)Sandia 國(guó)家實(shí)驗(yàn)室動(dòng)態(tài)失速下測(cè)得的風(fēng)力機(jī)實(shí)驗(yàn)數(shù)據(jù)與風(fēng)洞靜態(tài)實(shí)驗(yàn)結(jié)果[40-41],以獲得風(fēng)輪的最大切向力為目標(biāo),得到了垂直軸風(fēng)力機(jī)在上風(fēng)區(qū)和下風(fēng)區(qū)的最佳理論攻角。理論上,當(dāng)風(fēng)力機(jī)以最佳攻角運(yùn)行時(shí),將不再有動(dòng)態(tài)失速的現(xiàn)象產(chǎn)生?;陔p致動(dòng)盤多流管理論進(jìn)行Matlab編程計(jì)算,給出了垂直軸風(fēng)力機(jī)1周變槳距規(guī)律,變槳前后風(fēng)輪切向力對(duì)比如圖25所示。從圖25可知:利用該變槳距規(guī)律得到的風(fēng)能利用率可以由34.6%提高到42.8%。這說明這種變槳距規(guī)律能夠較好地提高垂直軸風(fēng)力機(jī)的切向力,進(jìn)而增加風(fēng)輪的轉(zhuǎn)矩[39]。
圖25 變槳前后風(fēng)輪所受切向力對(duì)比Fig.25 Contrast of tangential force before and after pitch angle adjustment
近幾年關(guān)于變槳距垂直軸風(fēng)力機(jī)的研究發(fā)展迅速,隨著流體理論的完善、計(jì)算流體技術(shù)的進(jìn)步和實(shí)驗(yàn)經(jīng)驗(yàn)的積累,垂直軸風(fēng)力機(jī)變槳技術(shù)將會(huì)越來越成熟。未來幾年,“大型化、高效率、低成本”仍將是垂直軸風(fēng)力機(jī)發(fā)展的主要趨勢(shì),盡管大型化垂直軸風(fēng)力機(jī)由于高度較高,所受風(fēng)剪效應(yīng)的影響較大[42-44],主軸偏振劇烈[45],但是其轉(zhuǎn)速下降,對(duì)變槳調(diào)節(jié)機(jī)構(gòu)的響應(yīng)時(shí)間可能要求變低,因此,新的流場(chǎng)模型的建立和空氣動(dòng)力學(xué)理論的提出,更準(zhǔn)確合理的槳距角調(diào)節(jié)規(guī)律和更可靠經(jīng)濟(jì)的變槳機(jī)構(gòu)以及智能控制系統(tǒng)的研發(fā),都將成為未來變槳距垂直軸風(fēng)力機(jī)的研究熱點(diǎn)和難點(diǎn)。
同時(shí),未來風(fēng)電發(fā)展也將會(huì)是大數(shù)據(jù)、人工智能、物聯(lián)網(wǎng)等新技術(shù)的時(shí)代,數(shù)據(jù)是一種潛力巨大、影響深遠(yuǎn)的能源,為盡快地利用這種能源,一些風(fēng)能公司創(chuàng)建了大數(shù)據(jù)監(jiān)控平臺(tái),利用多種傳感器融合技術(shù),廣泛采集機(jī)組運(yùn)行數(shù)據(jù),建立運(yùn)行數(shù)據(jù)庫和風(fēng)資源數(shù)據(jù)庫,預(yù)先感知運(yùn)行狀態(tài),為機(jī)組長(zhǎng)期穩(wěn)定高效運(yùn)行和優(yōu)化升級(jí)提供理論支撐。該技術(shù)在H 型變槳式垂直軸風(fēng)力機(jī)上應(yīng)用的可行性在于:通過大數(shù)據(jù)平臺(tái)的創(chuàng)建,可以快速根據(jù)來流風(fēng)風(fēng)速及風(fēng)向預(yù)先感知風(fēng)力機(jī)的流場(chǎng)情況,從而得到槳距角的調(diào)節(jié)值,調(diào)動(dòng)變槳機(jī)構(gòu)提前響應(yīng),以彌補(bǔ)響應(yīng)時(shí)間的遲滯性。
此外,為提高垂直軸風(fēng)力機(jī)風(fēng)能利用效率,除設(shè)計(jì)單一變槳機(jī)構(gòu)外,各國(guó)學(xué)者還致力于研究基于新材料的翼型智能柔性變形技術(shù)[46]、風(fēng)輪外圍大區(qū)域流場(chǎng)實(shí)時(shí)控制[47]及采用輔助設(shè)備[48]向高空收集高風(fēng)能等新思想,基于新材料和信息技術(shù),探索開發(fā)智能風(fēng)力機(jī),這些途徑都將為垂直軸風(fēng)力機(jī)的商業(yè)化應(yīng)用起到一定的推動(dòng)作用。
1)變槳技術(shù)的核心在于通過調(diào)節(jié)槳距角以使攻角維持在最佳角度,從而獲得較高的翼型升阻比以提高風(fēng)力機(jī)性能。目前對(duì)于垂直軸風(fēng)力機(jī)變槳規(guī)律的研究方法主要包括動(dòng)量模型法、數(shù)值計(jì)算法和測(cè)試試驗(yàn)法等,其中,數(shù)值計(jì)算法是H 型垂直軸風(fēng)力機(jī)變槳技術(shù)研究的主流方法。
2)無論是采用機(jī)械傳動(dòng)變槳,還是液壓傳動(dòng)變槳,局部方位角變槳或全方位角變槳未能應(yīng)用于工程實(shí)際中的主要難點(diǎn)在于:變槳控制系統(tǒng)的復(fù)雜性導(dǎo)致系統(tǒng)整體響應(yīng)時(shí)間較長(zhǎng),無法適應(yīng)風(fēng)速風(fēng)向的快速隨機(jī)變化;葉片攻角的不斷改變、方位角位置變化,以及實(shí)際工況下葉尖速比的波動(dòng),使得1周變槳規(guī)律實(shí)現(xiàn)很困難;葉片在上風(fēng)區(qū)和下風(fēng)區(qū)2次掃掠來流風(fēng),導(dǎo)致下風(fēng)區(qū)的誘導(dǎo)速度分布很復(fù)雜,給下風(fēng)區(qū)變槳規(guī)律的理論計(jì)算帶來了技術(shù)瓶頸;動(dòng)態(tài)失速下實(shí)時(shí)變槳控制規(guī)律的研究更復(fù)雜,目前僅限于理論研究,且受限于動(dòng)態(tài)氣動(dòng)性能相關(guān)實(shí)驗(yàn)數(shù)據(jù)的缺乏。
3)目前,商業(yè)化小型H 型垂直軸風(fēng)力機(jī)的風(fēng)能利用率為30%~35%,隨著變槳技術(shù)從原有的僅限于“限速”作用,逐漸地向局部方位角變槳,再到全方位角變槳方向演變,風(fēng)能利用率也將逐步提高。同時(shí),隨著人工智能、大數(shù)據(jù)、物聯(lián)網(wǎng)等新一代信息技術(shù)的發(fā)展,風(fēng)場(chǎng)外圍大區(qū)域的風(fēng)速、風(fēng)向等參數(shù)將會(huì)被提前預(yù)知;智能柔性變形機(jī)翼技術(shù)、智能材料與結(jié)構(gòu)和新型傳感控制技術(shù)的快速發(fā)展,將催生智能風(fēng)力機(jī)的誕生?;谶@些技術(shù)的融合,未來全方位角變槳技術(shù)真正應(yīng)用于大型垂直軸風(fēng)力機(jī)的可能性將會(huì)更大。