CN106777525B - 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法 - Google Patents

考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法 Download PDF

Info

Publication number
CN106777525B
CN106777525B CN201611050056.3A CN201611050056A CN106777525B CN 106777525 B CN106777525 B CN 106777525B CN 201611050056 A CN201611050056 A CN 201611050056A CN 106777525 B CN106777525 B CN 106777525B
Authority
CN
China
Prior art keywords
wind
blade
speed ratio
wind turbine
tip speed
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
Application number
CN201611050056.3A
Other languages
English (en)
Other versions
CN106777525A (zh
Inventor
殷明慧
陈小阳
杨志强
沈力
李政轩
李玉龙
卜京
谢云云
邹云
姚娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201611050056.3A priority Critical patent/CN106777525B/zh
Publication of CN106777525A publication Critical patent/CN106777525A/zh
Application granted granted Critical
Publication of CN106777525B publication Critical patent/CN106777525B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种考虑设计叶尖速比对风力机静、动态性能影响的风力机气动设计方法。该方法以风力机闭环性能指标——平均风能捕获效率作为目标函数，并将原先设计叶尖速比与弦长、扭角分离优化的流程改进为联合优化，使得原本忽略最大功率点跟踪（Maximum Power Point Tracking，MPPT）控制动态的风力机气动设计能够考虑到跟踪性能的不足，从而协调风力机的静态气动性能与跟踪动态过程。本发明所提方法能够有效提高变速风机在湍流风速下的风能捕获效率。

Description

考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法
技术领域
本发明属于风力机气动设计领域，特别是一种考虑设计叶尖速比对风力机静、动态性能影响的风力机气动设计方法。
背景技术
近年来，变速风机已成为大型风电机组的主流机型。变速风机大部分时间运行于低于额定风速的变速阶段。在该阶段，风机普遍采用最大功率点跟踪(Maximum PowerPoint Tracking，MPPT)控制策略。它会根据风速的变化实时调整风轮转速，使其始终运行于设计叶尖速比λdgn。因此，传统的针对变速风机的气动设计通常以λdgn处的功率系数极限Cp,max最大化为目标，如Glauert方法、Wilson方法以及其它基于叶素动量理论的逆设计方法。
但是，上述根据Cp,max确定λdgn的方法忽略了λdgn对MPPT动态过程的影响：
1、由于大转动惯量导致的慢动态性能，风轮的MPPT动态过程及其对风能捕获效率的影响是不容忽视的。变速风机常运行在跟踪λdgn的过程中而非λdgn处。且MPPT动态过程越长，风力机偏离λdgn运行的时间越长，这将导致较低的风能捕获效率。
2、除了上述的功率系数极限之外，λdgn还会影响MPPT动态过程。减小λdgn能够缩短风轮转速的跟踪范围，进而提升风能捕获效率；反之，则扩大跟踪范围，降低风能捕获效率。因此，如果仅为提高Cp,max而选择大的设计叶尖速比，则增大的MPPT过程反而使得风轮更长时间偏离设计叶尖速比，最终可能降低变速风机的闭环性能。
因此，针对风力机的气动逆设计，迫切需要一种能够综合考虑设计叶尖速比对静态气动性能和MPPT动态过程两方面影响的设计方法。但是现有技术中尚无相关描述。
发明内容
本发明的目的在于提供一种考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法。
实现本发明目的的技术解决方案为：一种考虑设计叶尖速比对风力机静、动态性能影响的风力机气动设计方法，包括以下步骤：
步骤1、进行初始化操作，具体为：
步骤1-1、确定初始设计参数，包括叶片数目B、叶片半径R、轮毂半径Rhub、翼型、弦长和扭角，并保持叶片数目B、叶片半径R、轮毂半径Rhub、翼型不变；
步骤1-2、根据初始风轮的叶片数目B、叶片半径R、轮毂半径Rhub、翼型、弦长和扭角，初始化PROPID程序；
步骤1-3、确定设计叶尖速比λdgn的初始搜索区间范围；
步骤2、将当前搜索区间等分为4个子区间；
步骤3、任选一个子区间的边界值λdgn；
步骤4、利用PROPID程序完成边界值为λdgn的气动逆设计，获得叶片气动外形参数——弦长与扭角，具体为：
步骤4-1、确定沿叶片展向的升力系数分布和诱导因子分布；
步骤4-2、执行PROPID程序，获得边界值λdgn对应的叶片弦长与扭角；
步骤5、计算目标函数，即平均风能捕获效率Pfavg；
步骤6、判断子区间边界值是否取完，若是，继续下一步；若否，返回步骤3；
步骤7、判断终止条件是否满足，如果不满足，生成新的搜索区间，然后返回步骤2；如果满足，则继续下一步；所述终止条件为子区间的长度小于或等于0.1；
步骤8、输出对应Pfavg最大值的优化风轮，流程结束。
本发明与现有技术相比，其显著优点为：1)本发明提出了一种考虑设计叶尖速比对风力机静、动态性能影响的风力机气动设计方法。该方法以风机闭环性能指标——平均风能捕获效率作为目标函数，并将原先设计叶尖速比与叶片气动外形参数(弦长和扭角分布)分离设计流程的改进为联合优化。2)本发明基于Bladed的仿真结果表明本发明所提方法能够有效提高变速风机在湍流风速下的风能捕获效率，验证了本发明方法的有效性。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为传统分离设计方法、本发明所提方法对应的优化风轮及初始风轮的叶片气动外形，其中图(a)为弦长对比图，图(b)为扭角对比图。
图2为优化风轮的功率系数曲线。
图3为优化风轮的运行叶尖速比的概率分布图。
具体实施方式
本发明的一种考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法，包括以下步骤：
步骤1、进行初始化操作，具体为：
步骤1-1、确定初始设计参数，包括叶片数目B、叶片半径R、轮毂半径Rhub、翼型、弦长和扭角，并保持叶片数目B、叶片半径R、轮毂半径Rhub、翼型不变；
步骤1-2、根据初始风轮的叶片数目B、叶片半径R、轮毂半径Rhub、翼型、弦长和扭角，初始化PROPID程序；
步骤1-3、确定设计叶尖速比λdgn的初始搜索区间范围；确定设计叶尖速比λdgn的初始搜索区间范围为5到9。
步骤2、将当前搜索区间等分为4个子区间；
步骤3、任选一个子区间的边界值λdgn；
步骤4、利用PROPID程序完成边界值为λdgn的气动逆设计，获得叶片气动外形参数——弦长与扭角，具体为：
步骤4-1、确定沿叶片展向的升力系数分布和诱导因子分布；沿叶片展向的升力系数分布为翼型最大升阻比附近对应的升力系数，诱导因子分布为0.3333。
步骤4-2、执行PROPID程序，获得边界值λdgn对应的叶片弦长与扭角；
步骤5、计算目标函数，即平均风能捕获效率Pfavg；所述目标函数为平均风能利用率Pfavg：
Figure BDA0001160570220000031
Figure BDA0001160570220000032
Pwy＝0.5ρπR2v3cos3ψ
n为一个统计时段内的采样次数，ψ为偏航误差角，Pcap为实际功率，Pwy为功率最优值。
步骤6、判断子区间边界值是否取完，若是，继续下一步；若否，返回步骤3；
步骤7、判断终止条件是否满足，如果不满足，生成新的搜索区间，然后返回步骤2；如果满足，则继续下一步；所述终止条件为子区间的长度小于或等于0.1；所述新的搜索区间为以当前Pfavg对应的λdgn为区间中点，左右区间长度为当前长度二分之一的区间。
步骤8、输出对应Pfavg最大值的优化风轮，流程结束。
本发明提出了一种考虑设计叶尖速比对风力机静、动态性能影响的风力机气动设计方法。该方法以风机闭环性能指标——平均风能捕获效率作为目标函数，并将原先设计叶尖速比与叶片气动外形参数(弦长和扭角分布)分离设计流程的改进为联合优化。
下面结合实施例对本发明做进一步详细的描述：
实施例
选取美国可再生能源实验室(National Renewable Energy Laboratory,NREL)研发的1.5MW风机作为基准，分别应用传统分离风轮优化设计方法和本发明提出优化方法进行气动优化设计。基准叶片参数见表1。
表1 1.5MW风力机叶片参数
Figure BDA0001160570220000041
计算Pfavg所用湍流风速参数如表2所示。
表2 湍流风速参数
Figure BDA0001160570220000042
Figure BDA0001160570220000051
按照传统分离设计原则(即以风轮的Cp,max最大化为目标去选择λdgn并优化风轮的气动外形)，通过调用PROPID程序完成不同λdgn的风轮气动设计，然后分别计算各风轮的Cp,max，从而获得Cp,max最大值所对应的λdgn及气动外形参数。
根据本发明提出的优化方法，搜索到Pfavg最大值所对应的优化风轮的气动参数。
将上述两个优化风轮分别命名为“分离优化风轮”、“联合优化风轮”。表3分别列出了两种风轮的静态气动性能Cp,max及其在同一湍流风速下的平均风能捕获效率Pfavg。
表3 不同风轮静态气动性能及闭环性能对比
分离优化风轮 联合优化风轮
λ<sub>dgn</sub> 8.5 6.2
C<sub>p,max</sub> 0.4790(+1.08％) 0.4739
P<sub>favg</sub> 0.4442 0.4547(+2.36％)
由表3可以看出：
1)分离优化风轮的静态气动性能优于联合优化风轮，前者的Cp,max较后者提高了1.08％；
2)联合优化风轮在湍流风速下的实际运行效率优于分离优化风轮，前者的Pfavg较后者提高了2.36％。
尽管分离优化风轮的Cp,max大于联合优化风轮，但是由于前者的运行叶尖速比分布在设计叶尖速比的概率较小，即前者以较大功率系数捕获风能的持续时间较短，最终导致分离优化风轮的Pfavg低于联合优化风轮。
综上，单点的静态气动性能最优并不能保证风力机在运行时的实际风能捕获效率最优。由于在设计时考虑到气动性能与MPPT动态过程的协调匹配，本发明公开的风轮气动设计方法可以进一步提高变速风机在湍流风速下的风能捕获效率。

Claims (5)

1.一种考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法，其特征在于，包括以下步骤：
步骤1、进行初始化操作，具体为：
步骤1-1、确定初始设计参数，包括叶片数目B、叶片半径R、轮毂半径Rhub、翼型、弦长和扭角，并保持叶片数目B、叶片半径R、轮毂半径Rhub、翼型不变；
步骤1-2、根据初始风轮的叶片数目B、叶片半径R、轮毂半径Rhub、翼型、弦长和扭角，初始化PROPID程序；
步骤1-3、确定设计叶尖速比λdgn的初始搜索区间范围；
步骤2、将当前搜索区间等分为4个子区间；
步骤3、任选一个子区间的边界值λdgn；
步骤4、利用PROPID程序完成边界值为λdgn的气动逆设计，获得叶片气动外形参数——弦长与扭角，具体为：
步骤4-1、确定沿叶片展向的升力系数分布和诱导因子分布；
步骤4-2、执行PROPID程序，获得边界值λdgn对应的叶片弦长与扭角；
步骤5、计算目标函数，即平均风能捕获效率Pfavg；
步骤6、判断子区间边界值是否取完，若是，继续下一步；若否，返回步骤3；
步骤7、判断终止条件是否满足，如果不满足，生成新的搜索区间，然后返回步骤2；如果满足，则继续下一步；所述终止条件为子区间的长度小于或等于0.1；
步骤8、输出对应Pfavg最大值的优化风轮，流程结束。
2.根据权利要求1所述考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法，其特征在于，步骤1-3中，确定设计叶尖速比λdgn的初始搜索区间范围为5到9。
3.根据权利要求1所述考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法，其特征在于，步骤4-1中，沿叶片展向的升力系数分布为翼型最大升阻比附近对应的升力系数，诱导因子分布为0.3333。
4.根据权利要求1所述考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法，其特征在于，步骤5所述目标函数为平均风能捕获效率Pfavg：
Figure FDA0002387737740000021
Figure FDA0002387737740000022
Pwy＝0.5ρπR2v3 cos3ψ
n为一个统计时段内的采样次数，ψ为偏航误差角，Pcap为实际功率，Pwy为功率最优值。
5.根据权利要求1所述考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法，其特征在于，步骤7中新的搜索区间为以当前Pfavg对应的λdgn为区间中点，左右区间长度为当前长度二分之一的区间。
CN201611050056.3A 2016-11-24 2016-11-24 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法 Active CN106777525B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611050056.3A CN106777525B (zh) 2016-11-24 2016-11-24 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611050056.3A CN106777525B (zh) 2016-11-24 2016-11-24 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法

Publications (2)

Publication Number Publication Date
CN106777525A CN106777525A (zh) 2017-05-31
CN106777525B true CN106777525B (zh) 2020-07-31

Family

ID=58912918

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611050056.3A Active CN106777525B (zh) 2016-11-24 2016-11-24 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法

Country Status (1)

Country Link
CN (1) CN106777525B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107154648B (zh) * 2017-06-09 2019-10-18 南京理工大学 一种风电场双层有功分配控制方法
CN107194122B (zh) * 2017-06-23 2021-01-15 南京理工大学 一种变速风机叶片的多工况气动优化的改进方法
CN109359426A (zh) * 2018-11-26 2019-02-19 南京理工大学 一种风力机叶片气动参数与控制器参数的联合优化方法
CN111209638B (zh) * 2020-01-17 2022-09-20 南京理工大学 基于运行攻角偏差修正的低风速风力机叶片气动设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102242689A (zh) * 2011-06-24 2011-11-16 南京理工大学 基于风力发电的最大功率点跟踪控制的改进爬山算法
CN102322407A (zh) * 2011-08-30 2012-01-18 重庆大学 风力机叶片气动外形协同设计方法
CN104879273A (zh) * 2015-06-12 2015-09-02 浙江运达风电股份有限公司 一种提高风力发电机组风能捕获的控制方法
CN105134487A (zh) * 2015-08-24 2015-12-09 南京理工大学 一种考虑湍流频率因素的风力机最大功率点跟踪控制方法
CN105909469A (zh) * 2016-04-29 2016-08-31 南京理工大学 一种减小变桨动作的变速恒频风电机组的限功率控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102242689A (zh) * 2011-06-24 2011-11-16 南京理工大学 基于风力发电的最大功率点跟踪控制的改进爬山算法
CN102322407A (zh) * 2011-08-30 2012-01-18 重庆大学 风力机叶片气动外形协同设计方法
CN104879273A (zh) * 2015-06-12 2015-09-02 浙江运达风电股份有限公司 一种提高风力发电机组风能捕获的控制方法
CN105134487A (zh) * 2015-08-24 2015-12-09 南京理工大学 一种考虑湍流频率因素的风力机最大功率点跟踪控制方法
CN105909469A (zh) * 2016-04-29 2016-08-31 南京理工大学 一种减小变桨动作的变速恒频风电机组的限功率控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"Improved interpolation method based on singular spectrum analysis iteration and its application to missing data recovery";WANG Hui-zan等;《Applied Mathmatics and Mechanics(English Edition)》;20081015;第29卷(第10期);正文第1351-1361页摘要，第2节 *
"NREL Advanced Research Turbine (ART) Aerodynamic Design of ART-2B Rotor Blades";Dayton A.Griffin;《citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.841.451&rep=rep1&type=pdf》;20000831;正文第20页第4节 *
"考虑风力机一体化设计的叶片气动外形关键参数的探索";杨志强等;《第三十三届中国控制会议》;20140730;正文第7072-7076页摘要，第4节 *
"风电机组最大功率点跟踪控制的影响因素分析";张小莲等;《电力系统自动化》;20131125;第37卷(第22期);正文第15-21页摘要，第3-4节 *

Also Published As

Publication number Publication date
CN106777525A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106777525B (zh) 考虑设计叶尖速比风轮静、动态影响的风力机气动设计方法
CN107559143B (zh) 一种大型风力机尾缘襟翼结构参数寻优及多目标襟翼优化控制方法
CN104141591B (zh) 风力发电最大功率点跟踪的改进自适应转矩控制方法
CN106919749B (zh) 一种低噪声风力机叶片设计方法及低噪声风力机叶片
CN108488035B (zh) 永磁直驱风力发电机组失速和变桨混合控制方法
CN103244348A (zh) 变速变桨风力发电机组功率曲线优化方法
CN102562453B (zh) 变速恒频风力发电机在额定转速阶段的变桨控制方法
CN109268205B (zh) 一种基于智能风力机的风电场优化控制方法
CN102434391A (zh) 一种基于起始转速调整的改进最大功率点跟踪控制方法
CN104675629B (zh) 一种变速风力发电机组的最大风能捕获方法
CN111209638B (zh) 基于运行攻角偏差修正的低风速风力机叶片气动设计方法
CN108223266B (zh) 在额定风速以下控制风力发电机的方法和装置
CN107194122B (zh) 一种变速风机叶片的多工况气动优化的改进方法
CN107218175B (zh) 一种实现风力机最大化风能捕获效率的转速跟踪目标优化方法
Wang et al. Optimized chord and twist angle distributions of wind turbine blade considering Reynolds number effects
CN109783912B (zh) 一种升力型垂直轴风力发电机叶片翼型优化设计方法
CN109488525B (zh) 基于提高转速下限的转速跟踪目标优化方法
Anstock et al. A new approach for comparability of two-and three-bladed 20 MW offshore wind turbines
CN103883483B (zh) 一种100w风力发电机叶片
Pourrajabian et al. On the effect of altitude on the performance of a small wind turbine blade
CN109359426A (zh) 一种风力机叶片气动参数与控制器参数的联合优化方法
Soraghan et al. Influence of lift to drag ratio on optimal aerodynamic performance of straight blade vertical axis wind turbines
CN111310382A (zh) 考虑运行多工况的变速风机叶片设计攻角优化方法
CN209855955U (zh) 一种低风速高性能风力机叶片
Poole et al. Optimization of a mini HAWT blade to increase energy yield during short duration wind variations

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