CN108397344B - 风力发电叶片桨矩角的确定方法 - Google Patents
风力发电叶片桨矩角的确定方法 Download PDFInfo
- Publication number
- CN108397344B CN108397344B CN201810074567.1A CN201810074567A CN108397344B CN 108397344 B CN108397344 B CN 108397344B CN 201810074567 A CN201810074567 A CN 201810074567A CN 108397344 B CN108397344 B CN 108397344B
- Authority
- CN
- China
- Prior art keywords
- wind
- formula
- wind wheel
- footmark
- rotation center
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 11
- 238000010248 power generation Methods 0.000 title claims abstract description 10
- 230000005611 electricity Effects 0.000 claims abstract description 7
- 241001672694 Citrus reticulata Species 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 12
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 230000002093 peripheral effect Effects 0.000 claims description 2
- 238000011835 investigation Methods 0.000 claims 1
- 238000000354 decomposition reaction Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D1/00—Wind motors with rotation axis substantially parallel to the air flow entering the rotor
- F03D1/06—Rotors
- F03D1/0608—Rotors characterised by their aerodynamic shape
- F03D1/0633—Rotors characterised by their aerodynamic shape of the blades
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了属于风力发电技术领域的一种风力发电叶片桨矩角的确定方法。该方法基于风轮上考察叶素某微段速度分析分解的几何关系,根据风力机叶片不同区域对翼型的要求不同,选择合适的翼型,根据风能利用最大化为目标,直接建立了风电叶片桨矩角沿叶片半径变化的公式,确定翼型沿叶片展向的布置;再以每一叶素处取得最大风能利用系数为目标,推导出相应的计算公式,通过优化计算,获取每一叶素处翼型的桨矩角。本发明不仅提高了叶片桨矩角的精确性,还有效提高了风能利用率。一般估算,在相同风况条件下,此方法确定的叶片桨矩角,风能利用率可提高近3%。
Description
技术领域
本发明属风力发电技术领域,特别涉及一种风力发电叶片桨矩角的确定方法。
背景技术
目前,风力发电叶片桨矩角的确定一般是通过将叶片分段,并将每个截面作为一个叶素,根据风力机叶片不同区域对翼型的要求不同,选择合适的翼型,即确定了翼型沿叶片展向的布置;再以每一叶素处取得最大风能利用系数为目标,通过优化计算,获取每一叶素处翼型的桨矩角;最后通过修形使各翼型光滑联接。
发明内容
本发明的目的是提出一种风力发电叶片桨矩角的确定方法,其特征在于,基于风轮上考察叶素某微段速度分析分解的几何关系,根据风力机叶片不同区域对翼型的要求不同,选择合适的翼型,根据风能利用最大化为目标,直接建立了风电叶片桨矩角沿叶片半径变化的公式如下:
式中:Vrr为距离风轮旋转中心为风轮半径r处的空气相对于风轮叶片的流动速度,其中,说明;第一个角标r表示“相对”运动,第二个角标r表示“位置”为距离风轮旋转中心的距离=风轮半径r;
Ver为将距离风轮旋转中心为风轮半径r处的空气看做风轮部分所具有的速度,即此处空气的牵连运动速度,
Va为空气的绝对运动速度,即风速;
式中:为矢量间的夹角;
式中:ω为风轮转速;
亦即
因为
将式(4)代入式(5),得
亦即
在上述各式中,Vrr为入流风速;Va为空气绝对运动速度;Vrc为出流风速; Vrz为出入流风速变化量在风轮旋转轴线z轴上的投影;Vra为出入流风速变化量; Vrx为出入流风速变化量在该点牵连运动速度方向x轴上的投影;ω为风轮旋转角速度;为矢量间的夹角;
根据流量连续方程Qcr=Qjr,其中,Qcr为r处的出流流量,第一个角标c表示“出”流,第二个角标r表示“位置”为距离风轮旋转中心r处;Qjr为r处的入流流量,第一个角标j表示“进”或入流,第二个角标r表示“位置”为距离风轮旋转中心r处)或ScrVcr=SjrVrr;其中,Scr为r处的出流流量通流截面面积,第一个角标c表示“出”流,第二个角标r表示“位置”为距离风轮旋转中心r 处;Sjr为r处的入流流量通流截面面积,第一个角标j表示“进”或入流,第二个角标r表示“位置”为距离风轮旋转中心r处,考虑空气在远低于音速条件下不可压缩的特点,应有Scr=Sjr,所以
Vcr=Vrr (8)
根据风轮上考察叶素某微段速度分析分解的几何关系,得
式中:为两矢量的矢量差,箭线代表矢量,其中,第一个角标a表示“差”,第二个角标r表示“位置”为距离风轮旋转中心r处。其它符号意义同前。亦即
式中:为距离风轮旋转中心r处的风电叶片桨矩角,角标r代表距离风轮旋转中心r处。其它符号意义同前。
式中:ΔVxr为矢量在x轴上的投影,其中,第一个角标x代表坐标轴,第二个角标r表示“位置”为距离风轮旋转中心r处。其它符号意义同前。
根据冲量原理,此微段空气流速的变化,对叶片产生的周向力ΔFxr为
ΔFxr=ρdqΔVxr (12)
式中:ΔFxr为风速在x方向的变化引起的作用力在x方向上的变化量,其中,第一个角标x代表坐标轴,第二个角标r表示“位置”为距离风轮旋转中心r处;Ρ为空气密度;dq为空气流量;
而
式中:ds为叶片上垂直于弦长的空气微单元厚度;dl为叶片上沿弦长的空气微单元长度。其它符号意义同前。
将式(10)代入式(11)后,再代入式(12)并将式(13)也代入式(12),得
此力对风轮产生的转矩ΔMxr为
式中:ΔMxr为力ΔFxr对风轮产生的转矩,其中,第一个角标x代表坐标轴,第二个角标r表示“位置”为距离风轮旋转中心r处;
其它符号意义同前。
在全叶片上产生的转矩Mx为
式中:Mx为全叶片上的作用力对风轮产生的转矩;其它符号意义同前;定义函数F,
因为
所以式(17)亦可写为
根据欧拉—拉格朗日方程(Euler—Lagrange Equation),式(16)取极大值的必要条件是
式(19)可简化为
则
本发明的有益效果:本发明根据风能利用最大化为目标,直接建立了风电叶片桨矩角沿叶片半径变化的公式,不仅提高了叶片桨矩角的精确性,还有效提高了风能利用率。一般估算,在相同风况条件下,此方法确定的叶片桨矩角,风能利用率可提高近3%。
附图说明
图1为风轮示意图,图中:r为考察叶素所在风轮半径;dr为叶素径向长度。
图2为风轮上考察叶素某微段速度分析分解示意图;图中:Vrr为入流风速; Va为空气绝对运动速度(额定风速);Vrc为出流风速;Vrz为出入流风速变化量在风轮旋转轴线(z轴)上的投影;Vra为出入流风速变化量;Vrx为出入流风速变化量在该点牵连运动速度方向(x轴)上的投影;ω为风轮旋转角速度;为矢量间的夹角;
各量上方的箭线表示该量的矢量。
具体实施方式
本发明提出一种风力发电叶片桨矩角的确定方法,下面结合附图予以说明。
实施例:
已知某待建风电场额定风速Va为12m/s,风电机组风轮额定转速ne为15rpm (ω=2πne=1.5608rad/s),风轮半径31m。则
1)风电机组所在场地额定风速Va=12m/s;
2)风轮额定转速ω=1.5608rad/s;
3)根据式(21)可得额定工况时不同风轮半径处的翼型桨矩角
将不同的r(0.5m≤r≤31m)代入式(22)可的沿风轮半径而变化的翼型桨矩角,具体数值见附表所示。
Claims (1)
1.一种风力发电叶片桨矩角的确定方法,其特征在于,基于风轮上考察叶素某微段速度分析分解的几何关系,根据风力机叶片不同区域对翼型的要求不同,选择合适的翼型,根据风能利用最大化为目标,直接建立了风电叶片桨矩角,风电叶片桨矩角沿叶片半径变化的公式如下:
式中:Vrr为距离风轮旋转中心为风轮半径r处的空气相对于风轮叶片的流动速度,其中,说明;第一个角标r表示“相对”运动,第二个角标r表示“位置”为距离风轮旋转中心的距离=风轮半径r;
Ver为将距离风轮旋转中心为风轮半径r处的空气看做风轮部分所具有的速度,即此处空气的牵连运动速度,
Va为空气的绝对运动速度,即风速;
式中:为矢量间的夹角;
式中:ω为风轮转速;
亦即
因为
将式(4)代入式(5),得
亦即
在上述各式中,Vrr为入流风速;Va为空气绝对运动速度;Vrc为出流风速;Vrz为出入流风速变化量在风轮旋转轴线z轴上的投影;Vra为出入流风速变化量;Vrx为出入流风速变化量在该点牵连运动速度方向x轴上的投影;ω为风轮旋转角速度;为矢量间的夹角;
根据流量连续方程Qcr=Qjr,其中,Qcr为r处的出流流量,第一个角标c表示“出”流,第二个角标r表示“位置”为距离风轮旋转中心r处;Qjr为r处的入流流量,第一个角标j表示“进”或入流,第二个角标r表示“位置”为距离风轮旋转中心r处,或ScrVcr=SjrVrr;其中,Scr为r处的出流流量通流截面面积,第一个角标c表示“出”流,第二个角标r表示“位置”为距离风轮旋转中心r处;Sjr为r处的入流流量通流截面面积,第一个角标j表示“进”或入流,第二个角标r表示“位置”为距离风轮旋转中心r处,考虑空气在远低于音速条件下不可压缩的特点,应有Scr=Sjr,所以
vcr=Vrr (8)
根据风轮上考察叶素某微段速度分析分解的几何关系,得
式中:为两矢量的矢量差,箭线代表矢量,其中,第一个角标a表示“差”,第二个角标r表示“位置”为距离风轮旋转中心r处,其它符号意义同前;
亦即
式中:为距离风轮旋转中心r处的风电叶片桨矩角,角标r代表距离风轮旋转中心r处,其它符号意义同前,
式中:ΔVxr为矢量在x轴上的投影,其中,第一个角标x代表坐标轴,第二个角标r表示“位置”为距离风轮旋转中心r处,其它符号意义同前;
根据冲量原理,此微段空气流速的变化,对叶片产生的周向力ΔFxr为
ΔFxr=ρdqΔVxr (12)
式中:ΔFxr为风速在x方向的变化引起的作用力在x方向上的变化量,其中,第一个角标x代表坐标轴,第二个角标r表示“位置”为距离风轮旋转中心r处;Ρ为空气密度;dq为空气流量;
而
式中:ds为叶片上垂直于弦长的空气微单元厚度;dl为叶片上沿弦长的空气微单元长度;其它符号意义同前;
将式(10)代入式(11)后,再代入式(12)并将式(13)也代入式(12),得
此力对风轮产生的转矩ΔMxr为
式中:ΔMxr为力ΔFxr对风轮产生的转矩,其中,第一个角标x代表坐标轴,第二个角标r表示“位置”为距离风轮旋转中心r处;
其它符号意义同前;
在全叶片上产生的转矩Mx为
式中:Mx为全叶片上的作用力对风轮产生的转矩;其它符号意义同前;
定义函数F,
因为
所以式(17)亦可写为
根据欧拉—拉格朗日方程(Euler—Lagrange Equation),式(16)取极大值的必要条件是
式(19)可简化为
则
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810074567.1A CN108397344B (zh) | 2018-01-25 | 2018-01-25 | 风力发电叶片桨矩角的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810074567.1A CN108397344B (zh) | 2018-01-25 | 2018-01-25 | 风力发电叶片桨矩角的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108397344A CN108397344A (zh) | 2018-08-14 |
CN108397344B true CN108397344B (zh) | 2019-10-22 |
Family
ID=63094975
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810074567.1A Active CN108397344B (zh) | 2018-01-25 | 2018-01-25 | 风力发电叶片桨矩角的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108397344B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL8203016A (nl) * | 1982-07-28 | 1984-02-16 | Wiebe Draijer | Werkwijze voor het vervaardigen van wieken voor axiale ventilatoren en windturbines. |
CN102400847B (zh) * | 2011-11-29 | 2013-06-19 | 吉林大学 | 一种风力机叶片翼型 |
KR101296674B1 (ko) * | 2011-11-30 | 2013-08-14 | 현대로템 주식회사 | 풍력발전기용 블레이드의 루트 에어포일 |
CN105545583B (zh) * | 2016-03-16 | 2018-03-20 | 华北电力大学 | 风力发电叶片及背风面出流切线倾角的确定方法 |
-
2018
- 2018-01-25 CN CN201810074567.1A patent/CN108397344B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108397344A (zh) | 2018-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107559143B (zh) | 一种大型风力机尾缘襟翼结构参数寻优及多目标襟翼优化控制方法 | |
US20100054936A1 (en) | Vertical axis wind turbine | |
CN105545595A (zh) | 基于径向基神经网络的风力机反馈线性化功率控制方法 | |
CN108488035B (zh) | 永磁直驱风力发电机组失速和变桨混合控制方法 | |
Cardona | Flow curvature and dynamic stall simulated with an aerodynamic free-vortex model for VAWT | |
CN113268901B (zh) | 基于格子Boltzmann动压气体轴承间隙微流动仿真方法 | |
Danao | The influence of unsteady wind on the performance and aerodynamics of vertical axis wind turbines | |
CN108799199A (zh) | 高负荷宽裕度弯度可调串列叶栅的设计方法 | |
CN108397344B (zh) | 风力发电叶片桨矩角的确定方法 | |
Hamlaoui et al. | Improved bem method for hawt performance predictions | |
Pirrung et al. | Dynamic stall model modifications to improve the modeling of vertical axis wind turbines | |
CN104838134B (zh) | 风能设备以及用于运行风能设备的方法 | |
CN106246459B (zh) | 一种风力机导风筒 | |
Driss et al. | Evaluation of the Savonius Wind Rotor Performance for Different External Overlap Ratios | |
CN108691730A (zh) | 风力发电机系统及风能数据误差修正方法和偏航控制方法 | |
Zhao et al. | Study on approach of performance improvement of VAWT employing double multiple stream tubes model | |
Eboibi et al. | A numerical study of the influence of blade profile and solidity on the performances of vertical axis wind turbines | |
Khan et al. | Performance of a HAWT Rotor with a Modified Blade Configuration. | |
Brunner et al. | Unsteady effects on a pitching airfoil at conditions relevant for large vertical axis wind turbines | |
Chen et al. | Application of numercal simulation to obtain the optimization pitch angle for VAWT | |
CN106837683B (zh) | 迎风面的出流切线倾角的最优值确定方法 | |
Coiro et al. | Flow curvature effect on dynamic behaviour of a novel vertical axis tidal current turbine: Numerical and experimental analysis | |
CN109026519B (zh) | 风电叶片、风轮及风电叶片弦长的确定方法 | |
CN106687684B (zh) | 确定液压机操作点的方法与转化液压能的装置 | |
SahuP et al. | Determination of Torque Produced by Horizontal Axis Wind Turbine Blade Using FSI Analysis for Low Wind Speed Regime |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |