姜海波 李艷茹 趙云鵬

海軍工程大學(xué)勤務(wù)學(xué)院，天津，300450

水平軸風(fēng)力機的轉(zhuǎn)矩性能及葉尖損失對其的影響

姜海波 李艷茹 趙云鵬

海軍工程大學(xué)勤務(wù)學(xué)院，天津，300450

分析了葉尖損失機理，根據(jù)PRANDTL 和GLAUERT葉尖損失修正因子及葉素-動量理論，推導(dǎo)出考慮葉尖損失的葉片弦長公式，然后沿葉片展向積分，推導(dǎo)出與風(fēng)力機尖速比和翼型升阻比關(guān)聯(lián)的風(fēng)力機轉(zhuǎn)矩系數(shù)解析表達式，并得到風(fēng)力機在任何穩(wěn)定運行狀態(tài)轉(zhuǎn)矩系數(shù)的最高參數(shù)值，可作為風(fēng)力機轉(zhuǎn)矩系數(shù)的設(shè)計參考值。研究表明，葉尖損失對弦長的影響集中在葉尖部位，當(dāng)升阻比為無窮大且尖速比從常見范圍10變化到4時，葉尖損失導(dǎo)致轉(zhuǎn)矩性能損失約4%～10%，且損失的數(shù)值隨升阻比的變化極小。

風(fēng)能利用；風(fēng)力機；轉(zhuǎn)矩系數(shù)；葉尖損失；葉素-動量理論

0 引言

水平軸風(fēng)力機最重要的性能指標(biāo)是功率、轉(zhuǎn)矩和推力。功率系數(shù)的研究已有大量文獻報道，其中BETZ于1920年給出了風(fēng)力機功率系數(shù)的理論極限是16/27，約為0.593[1-3]。姜海波等[4]則給出了穩(wěn)定運行狀態(tài)下與尖速比和升阻比關(guān)聯(lián)的理想風(fēng)力機功率系數(shù)的具體表達式，以及阻力為0(即升阻比為無窮大)時功率的最佳性能公式，并證明只有當(dāng)尖速比和升阻比均為無窮大時功率系數(shù)才會趨近于Betz極限，因此人們就不必為追求超越功率系數(shù)理論計算值而花費太多的精力，設(shè)計的風(fēng)力機效率接近相應(yīng)條件下的理論值時，就可認(rèn)定是理想設(shè)計。

風(fēng)力機的另一個重要參數(shù)是轉(zhuǎn)矩系數(shù)，研究轉(zhuǎn)矩系數(shù)是否存在理論上的極限值具有十分重要的意義。本文將利用葉素-動量理論通過推導(dǎo)給出風(fēng)力機最高轉(zhuǎn)矩系數(shù)計算公式，對此問題給出肯定的回答。

動量理論假設(shè)風(fēng)力機的葉片數(shù)為無窮多，流過風(fēng)輪的每個空氣粒子都與葉片相互作用，但實際風(fēng)力機的葉片數(shù)為有限個，流過風(fēng)輪的一部分粒子與葉片相互作用，另一部分會穿越葉片之間的空隙。這兩種情況下空氣傳遞給風(fēng)輪的動量是有區(qū)別的，即誘導(dǎo)速度不同。如果葉片處的誘導(dǎo)速度很大，那么入流角就會很小，導(dǎo)致轉(zhuǎn)矩在切向的分量減小，從而減小轉(zhuǎn)矩，造成功率損失。本文將在理論葉片弦長的基礎(chǔ)上考慮有限葉片數(shù)對葉尖損失造成的影響，給出了實際風(fēng)力機最高轉(zhuǎn)矩性能的計算方法和具體參數(shù)值。

1 基本關(guān)系式

首先簡要考察風(fēng)力機穩(wěn)定運行狀態(tài)(既非加速又非減速的運行狀態(tài))，風(fēng)的作用情況和葉片受力狀態(tài)(圖1、圖2)[4]。

圖1 風(fēng)力機葉輪順風(fēng)向示意圖Fig.1 Schematic view of a wind turbine with downwind direction

圖2 圖1之A-A截面葉素風(fēng)速狀態(tài)與受力分析Fig.2 Wind speed state and force analysis of the blade element at A-A cross section in Fig. 1

用與轉(zhuǎn)軸同心的圓柱切割葉片，刨面位置見圖1的A-A位置，在半徑r處取一微段(亦稱葉素)dr進行受力分析，如圖2所示。設(shè)ω為風(fēng)力機自轉(zhuǎn)角速度，r為葉片展向微段dr到轉(zhuǎn)軸中心的距離，R為葉片的長度。微段dr沿旋轉(zhuǎn)平面向上運動，故存在逆風(fēng)相對速度W=ωr。W與切向誘導(dǎo)速度bW的合成速度為w=(1+b)W。無窮遠(yuǎn)處來流絕對速度U與軸向誘導(dǎo)速度aU的合成速度u=(1-a)U是通過葉片的軸向風(fēng)速，u與w的合成速度v對葉片產(chǎn)生升力和阻力，速度v的攻角為α。升力dL垂直于合成速度v，可分解為周向分量dLw和水平分量dLu；阻力dD與v方向一致，可分解為周向分量dDw和水平分量dDu。φ為v與旋轉(zhuǎn)平面所成的夾角，稱來流角或入流角；α為翼弦與v的夾角，稱為攻角；設(shè)λ為r處W與U的比值，稱為線速度比。

對于葉尖損失，PRANDTL等[5]給出了一種近似的修正方法，經(jīng)過GLAUERT改進后的修正因子為[6]

(1)

式中，B為葉片數(shù)；λt為尖速比；a為軸向速度誘導(dǎo)因子；x為相對半徑，x=r/R。

(2)

(3)

葉片局部誘導(dǎo)速度因子為

(4)

(5)

圖3 葉尖損失因子與速度誘導(dǎo)因子沿展向的分布曲線Fig.3 Distribution curve of the blade tip loss factor and velocity induction factor along the span

(6)

假定風(fēng)力機運行在穩(wěn)定的工況，由圖2可知，風(fēng)的入流速度為

(7)

由此還可以得到入流角的正弦和余弦表達式：

(8)

(9)

2 考慮葉尖損失的弦長曲線

根據(jù)葉素理論[7]，由圖2和式(7)～式(9)可得葉素推力為

dT1=dLu+dDu=dLcosφ+dDsinφ=

(10)

式中，ρ為空氣密度；C為r處翼弦的長度；CL為升力系數(shù)；CD為阻力系數(shù)。

根據(jù)動量理論[8]，來流對風(fēng)輪圓盤中半徑為r到r+dr的圓環(huán)的推力為

(11)

每個圓環(huán)徑跨過B個葉片的葉素，所以dT2=BdT1，可得理想的相對弦長公式為

(12)

理想弦長可以通過解析計算的方法得到。從圖形來看這種扭曲的形狀很難加工生產(chǎn)，強度等性能也不能滿足實際流體環(huán)境的需要，因此實際葉片必須進行改造(同時也會喪失部分性能)。但是作為理想弦長公式，可以指導(dǎo)設(shè)計，并能在理論上發(fā)揮重要作用，例如可以利用其推導(dǎo)最高功率、轉(zhuǎn)矩和推力性能公式。

3 最高轉(zhuǎn)矩性能計算

考慮葉尖損失后，根據(jù)圖2，由式(7)～式(9)可得，穩(wěn)定運行工況的葉素轉(zhuǎn)矩為

(13)

式中，F(xiàn)為升力的合力。

對于B個葉片組成的風(fēng)力機，將相對弦長式(12)代入上式積分，得轉(zhuǎn)矩系數(shù)為

(14)

令阻力系數(shù)CD=0(升阻比ζ=∞)，可得到與尖速比關(guān)聯(lián)的轉(zhuǎn)矩系數(shù)為

(15)

(16)

圖4 阻力為0時葉尖損失修正后的轉(zhuǎn)矩系數(shù)變化趨勢Fig.4 Change trend of torque coefficient after correction of blade tip loss when resistance is 0

圖4所示為阻力為0但不考慮葉尖損失的最大轉(zhuǎn)矩系數(shù)曲線，式(16)和圖4的曲線顯示阻力為0時風(fēng)力機具有最高轉(zhuǎn)矩性能，文獻[9]推導(dǎo)出相同的結(jié)果，但未考慮葉尖損失的影響。將式(1)和式(2)代入考慮葉尖損失影響的式(15)進行數(shù)值積分，圖4中的點陣所示為葉尖損失修正后的轉(zhuǎn)矩系數(shù)變化趨勢，點陣和式(15)表明阻力為0時葉尖損失對風(fēng)力機的最高轉(zhuǎn)矩性能有影響。

由圖4可以看出，尖速比較小時，葉尖損失的存在會明顯降低轉(zhuǎn)矩系數(shù)；升阻比為無窮大時，最高轉(zhuǎn)矩系數(shù)損失情況見表1。

表1 阻力為0時葉尖損失對最高轉(zhuǎn)矩系數(shù)的影響

由表1可以看出，不考慮阻力影響的情況下，尖速比從常見范圍10變化到4時，葉尖損失導(dǎo)致轉(zhuǎn)矩性能損失約4%～10%。

在實際流體環(huán)境，阻力不可能為0(升阻比ζ為有限值)，由式(14)得

(17)

式(17) 是依據(jù)理想弦長公式推導(dǎo)得到的，實際葉片會依據(jù)加工難度、結(jié)構(gòu)強度等因素的影響對弦長分布進行修改，被設(shè)計加工成簡化葉片，從而會降低部分性能，因此式(17)為任何實際風(fēng)力機能達到的最高轉(zhuǎn)矩性能。

式(17)給定升阻比ζ和尖速比λt一系列參數(shù)值后進行數(shù)值計算，可得到葉尖損失修正后的最高轉(zhuǎn)矩系數(shù)值如表2所示。這些參數(shù)值是任何處于穩(wěn)定運行狀態(tài)的實際風(fēng)力機都不可超越的極限。

表2 葉尖損失修正后的最高轉(zhuǎn)矩系數(shù)值

(18)

式(18)比式(16) 的計算值略小一點，但兩者的圖像十分相似，在尖速比大于4以后這兩個圖像幾乎重合。雖然推導(dǎo)出風(fēng)力機的最高轉(zhuǎn)矩性能的解析表達式(18)，但未考慮葉尖損失的影響，因此式(18)理論價值較高，實用價值較低。由式(18)直接計算出未考慮葉尖損失的最高轉(zhuǎn)矩系數(shù)的理論值，并與表2的數(shù)值計算結(jié)果進行比較，計算分析表明，對于尖速比在4～10、升阻比在100以下的實用范圍內(nèi)的葉片，葉尖損失導(dǎo)致轉(zhuǎn)矩性能損失均在10%以內(nèi)，且具體數(shù)值隨尖速比變化較大，而隨升阻比的變化極小。

實際上，式(17) 和表2已給出葉尖損失修正后與升阻比和尖速比關(guān)聯(lián)的最高轉(zhuǎn)矩性能和具體參數(shù)值，可作為設(shè)計風(fēng)力機轉(zhuǎn)矩性能最實際的追求目標(biāo)，因此兼具重要的理論與應(yīng)用價值。

4 結(jié)論

(1)由葉素-動量理論可導(dǎo)出考慮葉尖損失的弦長公式，葉尖損失對弦長的影響主要集中在葉尖部位。

(2)轉(zhuǎn)矩沿翼展積分可導(dǎo)出考慮葉尖損失的最高轉(zhuǎn)矩性能計算公式，這是實際風(fēng)力機轉(zhuǎn)矩系數(shù)的最高值。

(3)不考慮阻力影響的情況下，尖速比從10變化到4時，葉尖損失導(dǎo)致轉(zhuǎn)矩性能損失約4%～10%。

(4)對于尖速比在4～10、升阻比在100以下的實用范圍內(nèi)，葉尖損失導(dǎo)致轉(zhuǎn)矩性能損失在10%以內(nèi)，且隨尖速比變化較大，而隨升阻比的變化極小。

[1] BETZ A. Das Maximum der Theoretisch M?glichen Ausnützung des Windes durch Windmotoren[J]. Zeitschrift für das Gesamte Turbinewesen,1920, 20: 307-309.

[2] BERGEY K H. The Lanchester-Betz limit[J]. Journal of Energy, 1979,3: 382-384.

[3] van KUIK G A M. The Lanchester-Betz-Joukowsky Limit[J]. Wind Energy, 2007(10): 289-291.

[4] 姜海波,曹樹良,陽平.水平軸風(fēng)力機的功率極限[J].機械工程學(xué)報, 2011, 47 (10): 113-118. JIANG Haibo, CAO Shuliang, YANG Ping. Power Limit of Horizontal Axis Wind Turbine[J]. Chinese Journal of Mechanical Engineering, 2011, 47 (10): 113-118.

[5] PRANDTL L, TIETJENS O G. Applied Hydro and Aeromechanics[M]. New York: Dover Publications, 1957.

[6] TONY B, DAVID S, NICK J, et al. Wind Energy Handbook[M]. Hoboken: John Wiley & Sons Ltd.,2005.

[7] 馬丁·漢森. 風(fēng)力機空氣動力學(xué)[M].2版.北京: 中國電力出版社,2009. MARTIN O L H. Aerodynamics of Wind Turbines[M]. 2nd ed. Beijing: China Electric Power Press, 2009.

[8] 賀德馨. 風(fēng)工程與工業(yè)空氣動力學(xué)[M]. 北京: 國防工業(yè)出版社, 2006: 81-185. HE Dexin. Wind Engineering and Industrial Aerodynamics[M]. Beijing:National Defense Industry Press, 2006: 81-185.

[9] JIANG Haibo, LI Yanru, CHENG Zhongqing. Performances of Ideal Wind Turbine[J]. Renewable Energy, 2015, 83(11) : 658-662.

(編輯胡佳慧)

TorquePerformanceofHorizontalAxisWindTurbinesandItsImpactsofBladeTipLosses

JIANG Haibo LI Yanru ZHAO Yunpeng
Logistics College,Naval University of Engineering,Tianjin，300450

This paper analyzed the mechanism of blade tip losses, and derived the blade chord formula according to PRANDTL and GLAUERT blade tip loss correction factor and blade element-momentum theory. With the blade chord formula, an analytical expression of torque coefficient associated with tip speed ratio and airfoil lift drag ratio of wind turbines was deduced by integrating along the wingspan using the blade element-momentum theory, and highest torque coefficient value of wind turbines in any steady state was obtained, which might be used for pre-estimating torque coefficient of actual wind turbines in design. The research shows, the impacts of blade tip losses to chord concentrated in the blade tip area, blade tip loss will reduce the torque coefficient about 4 to 10 percent when tip speed ratio changes from 10 to 4, and the value of the loss has minimal changes with lift drag ratio.

wind energy utilization; wind turbine; torque coefficient; blade tip loss; blade element-momentum theory

TK89

10.3969/j.issn.1004-132X.2017.18.009

2016-11-14

國家自然科學(xué)基金資助項目(51375489)

姜海波，男，1962年生。海軍工程大學(xué)勤務(wù)學(xué)院教授、博士研究生導(dǎo)師。主要研究方向為可再生能源開發(fā)利用。獲國家科技進步一等獎1項。出版專著1部，發(fā)表論文30余篇。E-mail：jianghaibo022@126.com。李艷茹，女，1974年生。海軍工程大學(xué)勤務(wù)學(xué)院副教授、博士。趙云鵬，男，1974年生。海軍工程大學(xué)勤務(wù)學(xué)院副教授。