姜海波 李艷茹 趙云鵬

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

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

姜海波 李艷茹 趙云鵬

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

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

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

0 引言

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

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

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

1 基本關(guān)系式

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

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

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

(1)

式中，B為葉片數(shù)；λt為尖速比；a為軸向速度誘導(dǎo)因子；x為相對(duì)半徑，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)力機(jī)運(yùn)行在穩(wěn)定的工況，由圖2可知，風(fēng)的入流速度為

(7)

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

(8)

(9)

2 考慮葉尖損失的弦長(zhǎng)曲線

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

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

(10)

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

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

(11)

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

(12)

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

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

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

(13)

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

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

(14)

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

(15)

(16)

圖4 阻力為0時(shí)葉尖損失修正后的轉(zhuǎn)矩系數(shù)變化趨勢(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時(shí)風(fēng)力機(jī)具有最高轉(zhuǎn)矩性能，文獻(xiàn)[9]推導(dǎo)出相同的結(jié)果，但未考慮葉尖損失的影響。將式(1)和式(2)代入考慮葉尖損失影響的式(15)進(jìn)行數(shù)值積分，圖4中的點(diǎn)陣所示為葉尖損失修正后的轉(zhuǎn)矩系數(shù)變化趨勢(shì)，點(diǎn)陣和式(15)表明阻力為0時(shí)葉尖損失對(duì)風(fēng)力機(jī)的最高轉(zhuǎn)矩性能有影響。

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

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

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

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

(17)

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

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

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

(18)

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

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

4 結(jié)論

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

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

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

(4)對(duì)于尖速比在4～10、升阻比在100以下的實(shí)用范圍內(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ī)的功率極限[J].機(jī)械工程學(xué)報(bào), 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)力機(jī)空氣動(dòng)力學(xué)[M].2版.北京: 中國(guó)電力出版社,2009. MARTIN O L H. Aerodynamics of Wind Turbines[M]. 2nd ed. Beijing: China Electric Power Press, 2009.

[8] 賀德馨. 風(fēng)工程與工業(yè)空氣動(dòng)力學(xué)[M]. 北京: 國(guó)防工業(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

國(guó)家自然科學(xué)基金資助項(xiàng)目(51375489)

姜海波，男，1962年生。海軍工程大學(xué)勤務(wù)學(xué)院教授、博士研究生導(dǎo)師。主要研究方向?yàn)榭稍偕茉撮_發(fā)利用。獲國(guó)家科技進(jìn)步一等獎(jiǎng)1項(xiàng)。出版專著1部，發(fā)表論文30余篇。E-mail：jianghaibo022@126.com。李艷茹，女，1974年生。海軍工程大學(xué)勤務(wù)學(xué)院副教授、博士。趙云鵬，男，1974年生。海軍工程大學(xué)勤務(wù)學(xué)院副教授。