CN112855464A - 一种基于特定场址湍流分布优化风电机组疲劳载荷的方法 - Google Patents
一种基于特定场址湍流分布优化风电机组疲劳载荷的方法 Download PDFInfo
- Publication number
- CN112855464A CN112855464A CN202110069725.6A CN202110069725A CN112855464A CN 112855464 A CN112855464 A CN 112855464A CN 202110069725 A CN202110069725 A CN 202110069725A CN 112855464 A CN112855464 A CN 112855464A
- Authority
- CN
- China
- Prior art keywords
- turbulence
- dlc
- turbine generator
- wind turbine
- fatigue
- 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.)
- Pending
Links
Images
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
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/84—Modelling or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Combustion & Propulsion (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Sustainable Energy (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Sustainable Development (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,包括以下步骤:S1.采集特定场址湍流数据;S2.结合IEC‑4标准计算若干分位数dlc1.2工况每个风速段对应的湍流;S3.计算包括dlc1.2在内的所有疲劳工况;S4.按比例分配各个分位数湍流下的时长权重;S5.进行后处理评估。上述技术方案根据特定场址湍流数据,结合IEC 61400‑1(Edition 4.0)标准,按照多种分位数计算出dlc1.2工况每个风速段对应的湍流,实现包括dlc1.2在内的所有疲劳工况的测量,从而针对性地分配各个分位数湍流下的时长权重,实现降低疲劳载荷的效果,对风电机组的使用条件和使用寿命的意义非常重大,可以大大拓展风电机组的适用情况和发电能力。
Description
技术领域
本发明涉及风力发电技术领域,尤其涉及一种基于特定场址湍流分布优化风电机组疲劳载荷的方法。
背景技术
有资料显示,据统计机械零件破坏的80%以上为疲劳破坏,特别是大型复杂机械结构的疲劳问题更为突出。风电机组是具有高动态载荷的动力***,所承受的载荷主要是随时间变化的动态载荷的动态随机载荷,结构件都会因此产生动应力,引起疲劳损伤。根据行业标准要求,风电机组的设计寿命至少要20年,因此在设计时需要重点考虑其疲劳强度,以确保整机结构的安全可靠性。
在风电机组设计领域,通常用等效疲劳载荷(Equivalent Fatigue Loads)来衡量疲劳载荷大小和进行结构件疲劳强度评估。目前在风电机组设计行业被应用最为广泛的等效疲劳载荷计算,是基于IEC 61400-1(Edition 4.0)标准,借助商业仿真软件GH Bladed中对计算后的疲劳载荷时序进行后处理计算。
IEC 61400-1(Edition 4.0)标准中计算的dlc1.2工况使用的是NTM正常湍流风模型,使用的是90%分位数的湍流。在实际评估中显得较为保守,不能充分发挥风力发电机组的使用功效。
中国专利文献CN101839806B公开了一种“风力发电机组及其疲劳载荷监测***”。采用了第一检测装置(211),用于检测塔筒根部的疲劳状态,得到第一检测信号,第二检测装置(212),用于检测除塔筒根部以外至少一个部位的疲劳状态,得到第二检测信号;控制装置(22),用于接收第一检测信号和第二检测信号,并根据预定策略处理第一检测信号和第二检测信号,得到实际损伤值,比较实际损伤值与预定损伤值,判断风力发电机组的被测部位的损伤状态;输出装置(23),用于输出被测部位的实际损伤值。上述技术方案难以准确有效的评估风电机组在特定厂址下的疲劳载荷,并且无法针对性地进行优化。
发明内容
本发明主要解决原有的技术方案无法准确有效的评估风电机组在特定厂址下的疲劳载荷的技术问题,提供一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,根据特定场址湍流数据,结合IEC 61400-1(Edition 4.0)标准,按照多种分位数计算出dlc1.2工况每个风速段对应的湍流,实现包括dlc1.2在内的所有疲劳工况的测量,从而针对性地分配各个分位数湍流下的时长权重,实现降低疲劳载荷的效果,对风电机组的使用条件和使用寿命的意义非常重大,可以大大拓展风电机组的适用情况和发电能力。
本发明的上述技术问题主要是通过下述技术方案得以解决的:本发明包括以下步骤:
S1.采集特定场址湍流数据;
S2.结合IEC-4标准计算若干分位数dlc1.2工况每个风速段对应的湍流;
S3.计算包括dlc1.2在内的所有疲劳工况;
S4.按比例分配各个分位数湍流下的时长权重;
S5.进行后处理评估。
现行的IEC 61400-1(Edition 4.0)标准中,疲劳载荷主要包括dlc1.2工况,dlc2.4工况,dlc4.1工况,dlc6.4工况,其中又以dlc1.2工况所占时长比重最长,对风电机组疲劳影响最大,而dlc1.2工况的设置直接取决于特定厂址的湍流大小,因此风电机组的疲劳载荷与湍流大小呈高度正相关。IEC61400-1(Edition 4.0)标准中计算的dlc1.2工况使用的是NTM正常湍流风模型,使用的是90%分位数的湍流。IEC 61400-1(Edition 4.0)标准中对NTM正常湍流风模型给出了以下公式:
σ1=Iref*(0.75*Vhub+b);b=5.6m/s
但是在实际评估疲劳时,90%分位数的湍流要远远大于实际的湍流风情况。IEC-4标准认为同一风速下湍流的分布服从威布尔分布,而以90%分位数的湍流对全年的疲劳载荷进行评估就过于保守了。本技术方案使用类似割圆法的方式,对同一风速下不同的湍流分布分别进行疲劳载荷的计算,然后根据各个不同湍流的分布比重加权获得更为准确的风电机组疲劳载荷结果。
作为优选,所述的步骤S2中的若干分位数包括5%,10%,15%,20%,25%,30%,35%,40%,45%,50%,55%,60%,65%,70%,75%,80%,85,90%,95%,99%分位数。
作为优选,所述的步骤S2根据Edition 4.0标准,对湍流在不同风速下的概率分布的定义有如下公式:
k=0.27*Vhub+1.4 (2)
C=Iref*(0.75*Vhub+3.3) (3)
将(2),(3)式代入(1)式,推导出以下公式:
同时PR依次代入5%,10%,15%,20%,25%,30%,35%,40%,45%,50%,55%,60%,65%,70%,75%,80%,85,90%,95%,99%分位数各个分位数,计算出dlc1.2工况在各个分位数、风速下的湍流值。
作为优选,所述的步骤S3中疲劳工况包括dlc1.2工况,dlc2.4工况,dlc4.1工况,dlc6.4工况。采用动力学计算软件Bladed计算出所有疲劳工况时序,除了dlc1.2工况以外,还有dlc2.4,dlc4.1,dlc6.4工况。
作为优选,所述的dlc1.2工况所占时长比重最长,对风电机组疲劳影响最大,设置好包括dlc1.2工况在内的所有工况,以此推导出不同分位数和风速下的湍流。由于该方法对dlc1.2工况的影响最大,本技术方案只改变dlc1.2工况的湍流设置。
以某公司的WD156-3300机组为例进行风电机组疲劳载荷优化对比。dlc1.2工况的风速依次是3,4,6,8,10,12,14,16,18,20m/s,本次计算的湍流分位数依次是5%,10%,15%,20%,25%,30%~95%,99%(由于100%分位数无穷大无法计算,所有以99%分位数替代)。以5%湍流计算0~5%湍流范围的疲劳载荷;以10%湍流计算5~10%湍流范围的疲劳载荷……以95%湍流计算90~95%湍流范围的疲劳载荷,以99%湍流计算95~100%湍流范围的疲劳载荷。
作为优选,所述的步骤S4通过雨流算法得到风电机组的整机疲劳载荷,然后乘以风速威布尔分布的时间权重,得出各个风速段的时长。
作为优选,所述的dlc1.2工况的每个风速段下的湍流分成20组,将对应风速段的时长分布均分成20份。
作为优选,所述的步骤S5后处理评估包括风电机组等效疲劳计算,风电机组等效疲劳计算通过动力学计算软件Bladed进行后处理计算后,得到20年载荷的马可夫矩阵,将累加次数等效值为1e7得到的载荷范围作为等效载荷。风电机组主要部件的疲劳载荷在原方法的基础上有10%左右的优化效果。
本发明的有益效果是:根据特定场址湍流数据,结合IEC 61400-1(Edition4.0)标准,按照多种分位数计算出dlc1.2工况每个风速段对应的湍流,实现包括dlc1.2在内的所有疲劳工况的测量,从而针对性地分配各个分位数湍流下的时长权重,实现降低疲劳载荷的效果,对风电机组的使用条件和使用寿命的意义非常重大,可以大大拓展风电机组的适用情况和发电能力。
附图说明
图1是本发明的一种流程图。
图2是现有IEC-4标准对正常湍流风模型的定义图。
图3是本发明的一种对dlc1.2工况的正常湍流风的重新定义后的主轴疲劳效果示意图。
图4是本发明的一种风电机组疲劳载荷优化对比图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:本实施例的一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,如图1所示,包括以下步骤:
S1.采集特定场址湍流数据;
S2.结合IEC-4标准计算若干分位数dlc1.2工况每个风速段对应的湍流;
S3.计算包括dlc1.2在内的所有疲劳工况;
S4.按比例分配各个分位数湍流下的时长权重;
S5.进行包括风电机组等效疲劳计算在内的后处理评估。
步骤S2中的若干分位数包括5%,10%,15%,20%,25%,30%,35%,40%,45%,50%,55%,60%,65%,70%,75%,80%,85,90%,95%,99%分位数,计算dlc1.2工况每个风速段对应的湍流过程如下:
现行的IEC 61400-1(Edition 4.0)标准中,疲劳载荷主要包括dlc1.2工况,dlc2.4工况,dlc4.1工况,dlc6.4工况,其中又以dlc1.2工况所占时长比重最长,对风电机组疲劳影响最大,而dlc1.2工况的设置直接取决于特定厂址的湍流大小,因此风电机组的疲劳载荷与湍流大小呈高度正相关。IEC61400-1(Edition 4.0)标准中计算的dlc1.2工况使用的是NTM正常湍流风模型,使用的是90%分位数的湍流。IEC 61400-1(Edition 4.0)标准中对NTM正常湍流风模型给出了以下公式:
σ1=Iref*(0.75*Vhub+b);b=5.6m/s
但是在实际评估疲劳时,90%分位数的湍流要远远大于实际的湍流风情况。IEC-4标准认为同一风速下湍流的分布服从威布尔分布,而以90%分位数的湍流对全年的疲劳载荷进行评估就未免显得过于保守了。见图2中标记处:
那么就可以使用类似割圆法的方式,对同一风速下不同的湍流分布分别进行疲劳载荷的计算,然后根据各个不同湍流的分布比重加权获得更为准确的风电机组疲劳载荷结果。
根据IEC 61400-1(Edition 4.0)标准,对湍流在不同风速下的概率分布的定义有如下公式:
k=0.27*Vhub+1.4 (2)
C=Iref*(0.75*Vhub+3.3) (3)
将(2),(3)式带入(1)式,通过上述公式可以推导出以下公式:
同时PR依次代入5%,10%,15%,20%,25%,30%~95%,99%各个分位数,计算出dlc1.2工况在各个分位数、风速下的湍流值,设置好包括dlc1.2工况在内的所有工况。以此就可以推导出不同分位数和风速下的湍流,能更精确的计算风电机组疲劳载荷。
步骤S3所述的计算包括dlc1.2在内的所有疲劳工况过程如下:
采用动力学计算软件Bladed计算出所有疲劳工况时序,除了dlc1.2工况以外,还有dlc2.4,dlc4.1,dlc6.4工况。
以某公司的WD156-3300机组为例进行风电机组疲劳载荷优化对比。dlc1.2工况的风速依次是3,4,6,8,10,12,14,16,18,20m/s,本次计算的湍流分位数依次是5%,10%,15%,20%,25%,30%~95%,99%(由于100%分位数无穷大无法计算,所有以99%分位数替代)。以5%湍流计算0~5%湍流范围的疲劳载荷;以10%湍流计算5~10%湍流范围的疲劳载荷……以95%湍流计算90~95%湍流范围的疲劳载荷,以99%湍流计算95~100%湍流范围的疲劳载荷。所有dlc1.2工况对主轴My的等效疲劳作用如图3所示。
步骤S4中按比例分配各个分位数湍流下的时长权重过程如下:
通过雨流算法得到风电机组的整机疲劳载荷,然后乘以风速威布尔分布的时间权重,因为dlc1.2工况的每个风速段下的湍流分成了20组,因此对应风速段的时长分布也需要均分成20份。
步骤S5中进行包括风电机组等效疲劳计算在内的后处理评估过程如下:通过动力学计算软件Bladed进行后处理计算后,得到20年载荷的马可夫矩阵,将累加次数等效值为1e7得到的载荷范围作为等效载荷。
从图3中可以明显看到,不同分位数下的湍流对风电机组的疲劳差异是相当大的,而按90%分位数的湍流计算的风电机组疲劳处在图中线框位置,相比之下,基于特定场址湍流分布计算风电机组疲劳载荷的方法明显更为准确,且具有降低疲劳载荷的效果,详见图4。
从图4可见,基于特定场址湍流分布计算风电机组疲劳载荷的方法,风电机组主要部件的疲劳载荷在原方法的基础上有10%左右的优化效果。这一降载效果对风电机组的使用条件和使用寿命的意义非常重大,可以大大拓展风电机组的适用情况和发电能力。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
尽管本文较多地使用了湍流数据、工况等术语,但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质;把它们解释成任何一种附加的限制都是与本发明精神相违背的。
Claims (8)
1.一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,其特征在于,包括以下步骤:
S1.采集特定场址湍流数据;
S2.结合IEC-4标准计算若干分位数dlc1.2工况每个风速段对应的湍流;
S3.计算包括dlc1.2在内的所有疲劳工况;
S4.按比例分配各个分位数湍流下的时长权重;
S5.进行后处理评估。
2.根据权利要求1所述的一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,其特征在于,所述步骤S2中的若干分位数包括5%,10%,15%,20%,25%,30%,35%,40%,45%,50%,55%,60%,65%,70%,75%,80%,85,90%,95%,99%分位数。
4.根据权利要求1所述的一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,其特征在于,所述步骤S3中疲劳工况包括dlc1.2工况,dlc2.4工况,dlc4.1工况,dlc6.4工况。
5.根据权利要求4所述的一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,其特征在于,所述dlc1.2工况所占时长比重最长,对风电机组疲劳影响最大,设置好包括dlc1.2工况在内的所有工况,以此推导出不同分位数和风速下的湍流。
6.根据权利要求1所述的一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,其特征在于,所述步骤S4通过雨流算法得到风电机组的整机疲劳载荷,然后乘以风速威布尔分布的时间权重,得出各个风速段的时长。
7.根据权利要求6所述的一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,其特征在于,所述dlc1.2工况的每个风速段下的湍流分成20组,将对应风速段的时长分布均分成20份。
8.根据权利要求1所述的一种基于特定场址湍流分布优化风电机组疲劳载荷的方法,其特征在于,所述步骤S5后处理评估包括风电机组等效疲劳计算,风电机组等效疲劳计算通过动力学计算软件Bladed进行后处理计算后,得到20年载荷的马可夫矩阵,将累加次数等效值为1e7得到的载荷范围作为等效载荷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110069725.6A CN112855464A (zh) | 2021-01-19 | 2021-01-19 | 一种基于特定场址湍流分布优化风电机组疲劳载荷的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110069725.6A CN112855464A (zh) | 2021-01-19 | 2021-01-19 | 一种基于特定场址湍流分布优化风电机组疲劳载荷的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112855464A true CN112855464A (zh) | 2021-05-28 |
Family
ID=76007371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110069725.6A Pending CN112855464A (zh) | 2021-01-19 | 2021-01-19 | 一种基于特定场址湍流分布优化风电机组疲劳载荷的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112855464A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110555543A (zh) * | 2018-05-31 | 2019-12-10 | 中国电力科学研究院有限公司 | 一种风电场预测的不确定度分析方法及*** |
CN111291514A (zh) * | 2020-02-03 | 2020-06-16 | 华北电力大学 | 一种风电机组疲劳寿命的快速预测的方法 |
CN111648916A (zh) * | 2019-03-04 | 2020-09-11 | 通用电气再生能源德国有限公司 | 对于风力涡轮的推力极限 |
-
2021
- 2021-01-19 CN CN202110069725.6A patent/CN112855464A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110555543A (zh) * | 2018-05-31 | 2019-12-10 | 中国电力科学研究院有限公司 | 一种风电场预测的不确定度分析方法及*** |
CN111648916A (zh) * | 2019-03-04 | 2020-09-11 | 通用电气再生能源德国有限公司 | 对于风力涡轮的推力极限 |
CN111291514A (zh) * | 2020-02-03 | 2020-06-16 | 华北电力大学 | 一种风电机组疲劳寿命的快速预测的方法 |
Non-Patent Citations (3)
Title |
---|
蒋红武 等: "风电场湍流强度对机组结构疲劳安全性的研究", 《应用能源技术》 * |
蔡继峰 等: "基于湍流强度分布进行机组疲劳累积的方法探讨", 《风能》 * |
高俊云: "对IEC 61400-1(第四版)中风电机", 《风能》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106815771B (zh) | 一种风电场载荷的长期评估方法 | |
JP6906354B2 (ja) | 風車発電機の疲労損傷量算出装置、風力発電システム、及び風車発電機の疲労損傷量算出方法 | |
CN112219028B (zh) | 风力涡轮机控制方法 | |
CN104915747A (zh) | 一种发电机组的发电性能评估方法及设备 | |
CN111291514B (zh) | 一种风电机组疲劳寿命的快速预测的方法 | |
CN106979126B (zh) | 基于svr的风力发电机组高风速段有效风速估计方法 | |
CN105654207A (zh) | 基于风速、风向信息的风电功率预测方法 | |
CN103745024B (zh) | 基于三维尾流模型修正风电机组尾部风速功率特性评估法 | |
CN108072524B (zh) | 一种风电机组齿轮箱轴承故障预警方法 | |
CN112459965A (zh) | 考虑风电场尾流的偏航优化控制方法、装置、设备及介质 | |
CN110414102A (zh) | 海上风力发电结构疲劳分析方法 | |
CN105701337A (zh) | 风电机组的疲劳寿命预测方法和装置 | |
CN110365053B (zh) | 基于延迟优化策略的短期风电功率预测方法 | |
CN104074687B (zh) | 一种用于兆瓦级风电机组的载荷及性能测试方法及装置 | |
EP3771822A1 (en) | A method for computer-implemented determination of a vertical speed wind profile of a wind field | |
CN111287911A (zh) | 一种风电机组疲劳载荷的预警方法和*** | |
CN112065668A (zh) | 一种风电机组状态异常评估方法及*** | |
CN114564885A (zh) | 一种基于bp神经网络的水轮发电机组轴瓦温度预警方法 | |
CN115221785A (zh) | 一种神经网络预测模型的构建方法、预测方法、*** | |
US11674500B2 (en) | Method for controlling a wind energy farm taking wake effects into account | |
CN105023099A (zh) | 一种考虑湍流强度的风力发电机出力评估方法 | |
CN112131753A (zh) | 风机疲劳寿命的评估方法、***、设备及可读介质 | |
CN112855464A (zh) | 一种基于特定场址湍流分布优化风电机组疲劳载荷的方法 | |
CN106771370B (zh) | 一种风机测风仪检测方法及装置 | |
CN115898787A (zh) | 一种风电机组静态偏航误差动态识别方法及装置 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210528 |
|
RJ01 | Rejection of invention patent application after publication |