CN101886619B - 风力发电机叶尖专用翼型 - Google Patents
风力发电机叶尖专用翼型 Download PDFInfo
- Publication number
- CN101886619B CN101886619B CN2010102188604A CN201010218860A CN101886619B CN 101886619 B CN101886619 B CN 101886619B CN 2010102188604 A CN2010102188604 A CN 2010102188604A CN 201010218860 A CN201010218860 A CN 201010218860A CN 101886619 B CN101886619 B CN 101886619B
- Authority
- CN
- China
- Prior art keywords
- airfoil
- coordinate
- aerofoil profile
- aerofoil
- blade tip
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- Y02E10/722—
Landscapes
- Wind Motors (AREA)
Abstract
本发明公开了一种风力发电机叶尖专用翼型,本基于翼型型线优化设计模型,以parsec翼型表达方法的9个控制参数为变量,设计出适用于叶片叶尖的翼型型线,在正常工作攻角范围内,具有很高的升力系数和升阻比,其失速性能非常平缓,并且具有很好的前缘粗糙度不敏感性,使得风力机即使在前缘结冰或被污染的情况下依然能够正常工作,很好的适应了风力机叶片叶尖部分的气动性能要求。
Description
技术领域
本发明涉及一种风力机部件,特别涉及一种风力发电机叶尖专用翼型。
背景技术
风能作为一种绿色能源,是解决能源问题、改善生态环境和减少二氧化碳排放的有效手段。风力机是将风能转化为电能等能源的中间设备,风力机依靠风轮叶轮汲取风能,叶轮的叶片翼型的气动性能直接影响着风力机风能利用的效率。
由于三维效应的影响,风力机叶片的叶尖区域往往产生复杂的流动现象,特别是叶尖涡的产生会使速度亏损值增加,风轮功率下降。在不增大风力机叶轮直径的情况下,对于风力机叶尖翼型性能的研究对提高风力机的发电功率有着非常重要的影响。
研究表明:对于大型风力机来说,具有较大的最大升力系数的翼型对减小风轮实度,增大启动力矩都有积极的作用。为使风力机的年输出功率得到提高,使单位发电量成本降低,在风力机叶片的叶尖部分应该使用拥有较高的最大升力系数的翼型。同时风力机翼型的阻力系数应该保持在一个较低的范围内,这体现为较高的升阻比系数。并且,由于叶尖翼型的失速特性对整个叶轮的动力学性能会产生巨大的影响,从而又影响到整个叶轮的结构设计,所以叶尖翼型的失速特性显得十分重要。为了减小叶片的气动力激励,如失速震颤等,应当使用具有较平缓失速性能的翼型。风力机风轮直径越大,高升力且拥有平缓失速性能的叶尖翼型就显得越为重要。另外,由于风力机叶片表面容易受到粉尘污染、昆虫污染、结冰和风蚀而破坏表面光洁度,使叶片表面变得粗糙,特别当翼型前缘粗糙时将导致翼面边界层转捩位置前移,转捩后边界层厚度增加,减少了翼型的弯度,从而减小最大升力系数,严重影响翼型的气动特性,所以还应该考虑翼型的前缘粗糙度敏感性。
然而,失速性能平缓的翼型的最大升力系数一般都比较小,而具有较大的最大升力系数的翼型又往往具有十分尖锐的失速性能。而即使两者都能够达到设计要求,翼型的气动性能仍有可能受到较大阻力系数或前缘粗糙度的影响。
因此,目前需要提出一种具有较大的最大升力系数,最大升阻比,失速性能平稳,而且前缘部分粗糙度增加后依然能够保持良好气动性能的叶尖翼型,以适用于风力发电机叶片。
发明内容
有鉴于此,本发明提供一种风力发电机叶尖专用翼型,具有较大的最大升力系数,最大升阻比,失速性能平稳,而且前缘部分粗糙度增加后依然能够保持良好气动性能的叶尖翼型,以适用于风力发电机叶片。
本发明的一种风力发电机叶尖专用翼型,翼型采用PARSEC翼型表达方法,方程为其中k=1时表示上翼面,k=2时表示下翼面,为翼面Z轴坐标,X为弦长方向的位置,为决定翼型形状的参数;将参数0.01<r<0.16、0.1<Xup<0.5、0.02<Zup<0.3、0<Zxxup<4、0.1<Xdw<0.5、0.05<Zdw<0.3、0<Zxxdw<5、0°<θ<30°和5°<β<30°代入可得以下2组6元方程组:
方程组一
方程组二
r-前缘半径;Xup--上翼面型线当Z坐标值取最大值时对应的X方向坐标值;Zup--上翼面型线的Z坐标最大值;Zxxup--上翼面型线当Z坐标取最大值时对应的曲率;Xdw--下翼面型线当Z坐标值取最小值时对应的X方向坐标值;Zdw--下翼面型线的Z坐标最小值;Zxxdw--下翼面型线当Z坐标取最小值时对应的曲率;θ--翼型尾缘处方向角;β--翼型尾缘角。
进一步,所述翼型优化模型采用基于非支配排序的NSGA Ⅱ遗传算法,设置初始种群大小c=50,最大进化代数g=30,交叉概率为Pc=1,变异概率为Pm=0.1。
本发明的有益效果:本发明的风力发电机叶尖专用翼型,基于翼型型线优化设计模型,以PARSEC翼型表达方法的9个控制参数为变量,设计出适用于叶片叶尖的翼型型线,在正常工作攻角范围内,具有很高的升力系数和升阻比,其失速性能非常平缓,并且具有很好的前缘粗糙度不敏感性,使得风力机即使在前缘结冰或被污染的情况下依然能够正常工作,很好的适应了风力机叶片叶尖部分的气动性能要求。
附图说明
下面结合附图和具体实施例对本发明作进一步描述。
图1为本发明的翼型图;
图2为本发明翼型的升力系数随攻角的变化曲线;
图3为本发明翼型的升阻比随攻角的变化曲线。
具体实施方式
图1为本发明的翼型图,图2为本发明翼型的升力系数随攻角的变化曲线,图3为本发明翼型的升阻比随攻角的变化曲线,如图所示:本实施例的风力发电机叶尖专用翼型,翼型采用PARSEC翼型表达方法,方程为其中k=1时表示上翼面,k=2时表示下翼面,为翼面Z轴坐标,X为弦长方向的位置,为决定翼型形状的参数;将参数0.01<r<0.16、0.1<Xup<0.5、0.02<Zup<0.3、0<Zxxup<4、0.1<Xdw<0.5、0.05<Zdw<0.3、0<Zxxdw<5、0°<θ<30°和5°<β<30°代入可得以下2组6元方程组:
方程组一
方程组二
r-前缘半径;Xup--上翼面型线当Z坐标值取最大值时对应的X方向坐标值;Zup--上翼面型线的Z坐标最大值;Zxxup--上翼面型线当Z坐标取最大值时对应的曲率;Xdw--下翼面型线当Z坐标值取最小值时对应的X方向坐标值;Zdw--下翼面型线的Z坐标最小值;Zxxdw--下翼面型线当Z坐标取最小值时对应的曲率;θ--翼型尾缘处方向角;β--翼型尾缘角。
本实施例中,翼型为工作在大中型风力机叶尖部分的翼型,选择设计工况为雷诺数Re=3.8×106,马赫数Ma=0.15。在翼型的设计过程中,选取PARSEC方法中的9个参数作为控制变量;
翼型的目标函数定义为:升力系数Cl、阻力系数Cd及失速特性Δ;
翼型在自由转捩条件下的最大升力系数ClFT、最大升阻比Cl/CdFT以及失速特性ΔFT;其中,失速特性通过最大升阻比附近的平方均差来表示:其中i为升阻比Cl/Cd取得最大值时对应的攻角a的取值。
考虑翼型的前缘粗糙度敏感性,将强制转捩条件下的最大升力系数ClCT、最大升阻比Cl/CdCT以及失速特性ΔCT也纳入优化模型的目标函数中去。
于是,优化模型的目标函数总共为6个:ClFT,Cl/CdFT,ΔFT,ClCT,Cl/CdCT,ΔCT。
优化模型采用基于非支配排序的NSGAII遗传算法,设置初始种群大小c=50,最大进化代数g=30,交叉与变异概率分别为Pc=1,Pm=0.1。
经过30代的优化得到一系列翼型。这里给出了在Z方向上的厚度与X方向的弦长的最大比值为0.1487的翼型,生成该翼型的PARSEC参数如下表:
r | Xup | Zup | Zxxup | Xdw | Zdw | Zxxdw | θ | β |
0.020664 | 0.35 | 0.085616 | 1.2605 | 0.24 | -0.066051 | 2.0055 | 11.85 | 7.4244 |
该翼型的几何参数如下表:
前缘半径 | 最大厚度 | 最大厚度发生位置(X方向) | 最大弯度 | 最大弯度发生位置(X方向) |
0.020664 | 0.14871 | 0.26 | 0.040884 | 0.66 |
使用XFOIL软件作为流体求解器,对本发明提出的翼型的气动性能进行了计算分析。雷诺数Re=3.8×106时,该翼型在自由转捩和强制转捩情况下的ClFT,Cl/CdFT,ΔFT,ClCT,Cl/CdCT,ΔCT见下表:
翼型的基本气动性能:
最大升力系数Cl | 最大升阻比Cl/Cd | 失速特性平方差 | |
自由转捩 | 2.0211 | 116.49 | 23.688 |
强制转捩 | 2.0211 | 116.49 | 122.29 |
自由转捩气动性能:
攻角 | 升力系数 | 阻力系数 | 升阻比 |
0 | 0.5139 | 0.00557 | 92.262 |
1 | 0.6241 | 0.00611 | 102.14 |
2 | 0.7336 | 0.00672 | 109.17 |
3 | 0.8421 | 0.00742 | 113.49 |
4 | 0.9484 | 0.00826 | 114.82 |
5 | 1.054 | 0.00914 | 115.32 |
6 | 1.1589 | 0.00998 | 116.12 |
7 | 1.2627 | 0.01084 | 116.49 |
8 | 1.3646 | 0.01172 | 116.43 |
9 | 1.4637 | 0.01271 | 115.16 |
10 | 1.5593 | 0.0138 | 112.99 |
11 | 1.6501 | 0.01502 | 109.86 |
12 | 1.7345 | 0.01638 | 105.89 |
13 | 1.8005 | 0.01788 | 100.7 |
14 | 1.8596 | 0.01964 | 94.684 |
15 | 1.9096 | 0.02211 | 86.368 |
16 | 1.9551 | 0.02513 | 77.799 |
17 | 1.9899 | 0.02927 | 67.984 |
18 | 2.0109 | 0.03521 | 57.112 |
19 | 2.0211 | 0.04331 | 46.666 |
20 | 2.0122 | 0.05538 | 36.334 |
强制转捩气动性能表:
攻角 | 升力系数 | 阻力系数 | 升阻比 |
0 | 0.4926 | 0.00762 | 64.646 |
1 | 0.6067 | 0.00781 | 77.682 |
2 | 0.7202 | 0.00805 | 89.466 |
3 | 0.8329 | 0.00833 | 99.988 |
4 | 0.9448 | 0.00863 | 109.48 |
5 | 1.054 | 0.00914 | 115.32 |
6 | 1.1589 | 0.00998 | 116.12 |
7 | 1.2627 | 0.01084 | 116.49 |
8 | 1.3646 | 0.01172 | 116.43 |
9 | 1.4637 | 0.01271 | 115.16 |
10 | 1.5593 | 0.0138 | 112.99 |
11 | 1.6501 | 0.01502 | 109.86 |
12 | 1.7345 | 0.01638 | 105.89 |
13 | 1.8005 | 0.01788 | 100.7 |
14 | 1.8596 | 0.01964 | 94.684 |
15 | 1.9096 | 0.02211 | 86.368 |
16 | 1.9551 | 0.02513 | 77.799 |
17 | 1.9899 | 0.02927 | 67.984 |
18 | 2.0109 | 0.03521 | 57.112 |
19 | 2.0211 | 0.04331 | 46.666 |
20 | 2.0122 | 0.05538 | 36.334 |
由以上三个表可见,本发明的翼型拥有较大的最高升力系数2.0122。当攻角为7°时,拥有最大升阻比116.49。由于在攻角大于5°以后,自由转捩的转捩发生点与强制转捩相同,导致其气动性能没有差别,所以该翼型具有十分良好的粗糙度不敏感性。自由转捩与强制转捩两种情况下的失速特性平方差都较小,分别为23.688和122.29;而且由图3也可直观的看出,该翼型的升阻比在最大升阻比附近的变化非常平稳。综上所述,本发明的翼型基于以上的翼型型线优化设计模型,以PARSEC翼型表达方法的9个控制参数为变量,以自由转捩和强制转捩条件下的最大升力系数、最大升阻比以及失速特性为目标函数,基于NSGA Ⅱ遗传算法,完全符合NREL所提出的风力发电机叶片叶尖翼型的要求。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种风力发电机叶尖专用翼型,其特征在于:翼型采用PARSEC翼型表达方法,方程为 其中k=1时表示上翼面,k=2时表示下翼面, 为翼面Z轴坐标,X为弦长方向的坐标, 为决定翼型形状的参数;将参数0.01<r<0.16、0.1<Xup<0.5、0.02<Zup<0.3、0<Zxxup<4、0.1<Xdw<0.5、0.05<Zdw<0.3、0<Zxxdw<5、0°<θ<30°和5°<β<30°代入 可得以下2组6元方程组:
方程组一
方程组二
r-前缘半径;Xup--上翼面型线当Z坐标值取最大值时对应的X方向坐标值;Zup--上翼面型线的Z坐标最大值;Zxxup--上翼面型线当Z坐标取最大值时对应的曲率;Xdw--下翼面型线当Z坐标值取最小值时对应的X方向坐标值;Zdw--下翼面型线的Z坐标最小值;Zxxdw--下翼面型线当Z坐标取最小值时对应的曲率;θ--翼型尾缘处方向角;β--翼型尾缘角;
3.根据权利要求2所述的风力发电机叶尖专用翼型,其特征在于:翼型优化模型采用基于非支配排序的NSGA II遗传算法,设置初始种群大小c=50,最 大进化代数g=30,交叉概率为Pc=1,变异概率为Pm=0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102188604A CN101886619B (zh) | 2010-07-07 | 2010-07-07 | 风力发电机叶尖专用翼型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102188604A CN101886619B (zh) | 2010-07-07 | 2010-07-07 | 风力发电机叶尖专用翼型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101886619A CN101886619A (zh) | 2010-11-17 |
CN101886619B true CN101886619B (zh) | 2012-01-18 |
Family
ID=43072667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102188604A Expired - Fee Related CN101886619B (zh) | 2010-07-07 | 2010-07-07 | 风力发电机叶尖专用翼型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101886619B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102235325B (zh) * | 2011-07-01 | 2012-10-31 | 重庆大学 | 基于翼型集成和混合尾缘改型的风力机叶尖翼型设计方法 |
CN104033331B (zh) * | 2014-06-06 | 2017-12-08 | 株洲时代新材料科技股份有限公司 | 一种风力发电叶片族系设计方法 |
CN104018999B (zh) * | 2014-06-18 | 2016-11-23 | 西北工业大学 | 一种用于兆瓦级风力机叶片的25%厚度主翼型 |
CN111859801B (zh) * | 2020-07-16 | 2022-07-19 | 湖北工业大学 | 失速型风力机翼型的设计方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100443720C (zh) * | 2005-11-09 | 2008-12-17 | 沈阳航空工业学院 | 大弯度风力机翼型 |
WO2007123552A1 (en) * | 2006-04-26 | 2007-11-01 | Midwest Research Institute | Adaptive pitch control for variable speed wind turbines |
EP2109713B1 (en) * | 2007-01-29 | 2013-07-24 | Bladena ApS | Wind turbine blade |
CN101615216B (zh) * | 2009-07-28 | 2011-01-05 | 重庆大学 | 对翼型尾缘进行流线形增厚的方法 |
-
2010
- 2010-07-07 CN CN2010102188604A patent/CN101886619B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101886619A (zh) | 2010-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101923584B (zh) | 风力机专用翼型设计方法及风力机专用翼型 | |
CN101532906A (zh) | 风力发电机叶片的流体动力学和结构力学分析方法 | |
CN101615216B (zh) | 对翼型尾缘进行流线形增厚的方法 | |
CN203374428U (zh) | 一族大厚度钝尾缘风力机翼型 | |
CN102094767B (zh) | 一族用于兆瓦级风力机叶片的翼型 | |
CN101458735A (zh) | 一种具有高升阻比的翼型 | |
CN101886619B (zh) | 风力发电机叶尖专用翼型 | |
CN106089569A (zh) | 一种适用于低雷诺数流动的小型风力机叶片翼型 | |
CN102003332B (zh) | 风力机叶片翼型族 | |
CN104863799A (zh) | 一种利用贝塞尔函数曲线的风力机翼型设计方法 | |
CN102444540B (zh) | 一种水平轴风力发电机的风轮叶片翼型 | |
CN103133272A (zh) | 一种大型风机的薄翼型叶片 | |
CN104018998B (zh) | 一种用于兆瓦级风力机叶片的21%厚度主翼型 | |
CN203383990U (zh) | 一种大型风机的大厚度翼型叶片 | |
CN103216381B (zh) | 一种风力发电机组叶片 | |
Sutikno et al. | Design and blade optimization of contra rotation double rotor wind turbine | |
CN103939283A (zh) | 垂直轴风力机专用高效叶片 | |
CN105781904A (zh) | 适用于兆瓦级风力机叶片的30%厚度翼型 | |
CN102278272B (zh) | 一种前突型水平轴风力机叶片 | |
CN104018999B (zh) | 一种用于兆瓦级风力机叶片的25%厚度主翼型 | |
CN201116510Y (zh) | 装有喷气襟翼翼型叶片的风力机叶轮 | |
CN203175769U (zh) | 一种大型风机的薄翼型叶片 | |
CN203321756U (zh) | 一种大型风机的中等厚度翼型叶片 | |
CN202348553U (zh) | 一种水平轴风力发电机的风轮叶片翼型 | |
CN106227985A (zh) | 潮流能水平轴水轮机用叶片翼型族设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120118 Termination date: 20140707 |
|
EXPY | Termination of patent right or utility model |