Wind power tower damping vibration attenuation device whole life period vibration attenuation income evaluation method

Abstract

The invention belongs to the technical field of wind power generation, and relates to a whole-life vibration reduction income evaluation method for a damping vibration reduction device of a wind power tower. According to the method, starting from the full-life-cycle and full-load working condition dynamic response of a wind power tower structure, the benefits of the vibration reduction device are quantitatively evaluated by quantitatively calculating the full-life-cycle failure risk and the loss cost value of the wind power tower loaded with the vibration control device based on the angles of wind power tower structure design and risk evaluation; a reliable basis is provided for selection and optimization of parameters of the damping vibration attenuation device; therefore, the overall safety and economy of the wind power structure are improved, and cost reduction and benefit increase of the structure are achieved.

Description

Technical Field

[0001] The invention belongs to the technical field of wind power generation, and relates to a method for evaluating the vibration reduction benefits of a wind power tower damping and vibration reduction device over its entire life cycle. Background Art

[0002] Wind power generation is one of the important components of the green and low-carbon new energy system. As the single-unit capacity and hub height of wind turbines continue to increase, and as wind turbines continue to expand into complex geographical environments such as mountainous areas, plateaus, and deep seas, their natural vibration frequencies continue to decrease and the loads they are subjected to continue to increase. The dynamic response under the coupling of multiple physical fields such as wind, waves, and earthquakes becomes more complex, and the risk of abnormal vibration is high.

[0003] In response to the increasingly serious problems of abnormal vibration of wind turbine towers, academia and industry have proposed a variety of wind turbine tower vibration reduction devices based on dampers, referring to existing technologies and experience in the field of seismic resistance of building structures. These devices are used to reduce the dynamic response of wind turbine towers, achieve structural load reduction and vibration suppression, and improve structural safety.

[0004] However, the current technology for wind turbine vibration reduction devices mainly focuses on the structural innovation of the device itself. As for how to evaluate the overall benefits of vibration reduction devices on wind turbine tower structures, there are no reports so far, and related research only discusses the load reduction and vibration suppression effects of a specific wind turbine under certain specific load conditions. In fact, to achieve the benefit evaluation of wind turbine vibration reduction devices, it is necessary to conduct load reduction and vibration suppression efficiency research on the full life cycle and full load conditions of wind turbine structures, and to refine vibration reduction benefit evaluation methods from the perspectives of wind turbine tower structure design and risk assessment. Summary of the invention

[0005] In view of the lack of a full-life benefit evaluation method for wind tower damper vibration reduction devices, the present invention starts from the full life cycle of the wind tower structure and the dynamic response of all load conditions, and proposes a full-life benefit evaluation method for vibration reduction devices based on the wind tower structure design and risk assessment perspectives, thereby filling the current gap and providing a methodological reference for the industry and academia.

[0006] The technical solution provided by the present invention is: a method for evaluating the vibration reduction benefits of a wind tower damper during its entire life cycle, comprising the following steps: S1. Determine the basic parameters of the wind tower structure to be evaluated, including tower material, hub height, tower cross-sectional shape and size, weight and eccentricity of the wind turbine impeller and nacelle, and calculate the dynamic response and load effect of the wind tower structure under various working conditions; S2. Determine the indicators and limits of the normal use limit state and the bearing capacity limit state of the wind power tower. For various limit state indicators, determine the corresponding limit values ​​one by one according to the design requirements; S3. Determine the corresponding cost loss after various limit state indicators are exceeded; S4. Obtaining the probability distribution of environmental parameters at the location of the wind tower, including the probability distribution of wind speed and the probability distribution of earthquake intensity; S5. Conduct a vulnerability analysis of the wind turbine tower, and calculate the exceedance probability of a certain limit state index of the wind turbine tower under wind load conditions and the exceedance probability of a certain limit state index of the wind turbine tower under earthquake conditions for a certain limit state index of the wind turbine tower; S6. Calculate the failure risk of wind towers, including the failure risk of the limit state corresponding to the limit state index and total failure risk ; S7, calculating the optimal values ​​of the parameters of the wind tower damping and vibration reduction device; S8. Calculate the average load reduction and vibration suppression efficiency of each limit state index of the wind tower support after installing the damping vibration reduction device based on the optimal value of the damper parameters. : ; In the formula, and are the characteristic values ​​of the limit state indicators of the wind turbine tower without dampers and after installing dampers, respectively. For the limit state indicator of normal use, its characteristic value is the root mean square of the time history calculation result, and for the limit state indicator of bearing capacity, its characteristic value is the maximum value; S9. Failure risk of limit state corresponding to limit state index Total failure risk The weight coefficient of each limit state index is calculated based on the ratio of ; S10. Calculate the total load reduction and vibration suppression efficiency of the wind tower after installing the damping vibration reduction device The formula is: .

[0007] Preferably, in step S2, calculations are carried out according to various design load conditions specified in the IEC61400-1 specification; the normal use limit state of the structure is verified under normal operating conditions, and the normal use limit state indicators mainly include tower top cabin acceleration, tower top cabin displacement, structural mud surface angle, and mud surface displacement; the structural bearing capacity limit state is verified under extreme conditions, earthquake conditions, or fatigue conditions, and the bearing capacity limit state indicators mainly include structural member design stress ratio, node stress, and fatigue stress amplitude.

[0008] Preferably, in step S3, when the normal use limit state indicator is exceeded, the power generation efficiency will be reduced and the maintenance cost will be increased. The cost loss that may be incurred by the wind turbine is: Cost m = Single maintenance fee + temporary measures fee + downtime loss fee = a1 + b 1 + b 2 + b 3 + c 1 × c 2 %× c 3 ÷365× c 4 ; in: a 1 For labor and material costs; b 1 For transportation costs; b 2 For security assessment fees; b 3 For facility rental fees; c 1 is the annual power generation, kW; c 2 % is the power generation discount rate; c 3 Number of days of downtime for maintenance; c 4 is the revenue per kilowatt of electricity generated, ¥ / kW.

[0009] Preferably, in step S3, when the bearing capacity limit state index is exceeded, the wind tower collapse cost will be caused, and the possible cost loss of the wind turbine is: Cost m = Equipment fee + Cleaning fee + Downtime loss fee - Recycling fee = d 1 %× d 2 × d 3 × d 4 %× d 5 × d 6 +( e 1 + e 2 + e 3 + e 4 + e 5 + e 61 × e 62 %× e63 )+ f 1 × f 2 %× f 3 ÷365× f 4 - g 1 × g 2 %× g 3 ; in, d 1 % is the whole machine discount rate; d 2 is the actual power generation duration in one year, h; d 3 is the annual power generation, kW; d 4 % is the power generation discount rate; d 5 is the length of service, years; d 6 is the levelized cost of electricity, ¥ / kW·h; e 1 For processing fees; e 2 For transportation costs; e 3 for driver compensation; e 4 For lifting fee; e 5 Equipment rental fees; e 61 is the blade mass, t; e 62 % is the landfill reduction rate; e 63 % is the landfill fee per ton, ¥ / t; f 1 is the annual power generation, kW; f 2 % is the power generation discount rate; f 4 is the revenue per kilowatt of electricity generated, ¥ / kW; f 3 days for repair and reconstruction; g 1 is the total mass of the tower and generator, t; g 2 % is the recycling reduction rate; g 3 It is the recycling price per ton of cast iron, ¥ / t.

[0010] Preferably, in step S4, for the wind load, the wind speed is first calculated Corresponding wind pressure : (1) In the formula is the air density; Next, calculate the wind pressure Return period , that is, wind speed Return period : ; In the formula and is the basic wind pressure at 10- and 100-year return periods in different regions; Calculate the design life of the supporting structure Internal wind speed The probability of exceeding : ; In the formula is the strength index; Perform forward difference to obtain discrete wind speed The probability of wind load occurrence ,Right now: ; For earthquake action, first calculate the design life Internal earthquake intensity The exceedance probability is: ; In the formula is the upper limit of earthquake intensity, which is 12; For the moderate intensity, take the basic intensity minus 1.55 degrees; for 50 years; are the shape parameters of different regions, which can be calculated based on the 50-year probability of occurrence of the seismic fortification intensity of 10%; Perform forward difference to obtain discrete earthquake intensity The probability of earthquake action : .

[0011] Preferably, in step S5, when considering the wind turbine load, the probability of occurrence of each DLC is The calculation formula is: ; in, is the dimension of the design load condition, indicating the design load cases; The probability of exceeding a certain limit state indicator of a wind turbine tower under full load conditions The calculation formula is:

[0012] When considering earthquake action, the probability of exceeding a certain limit state indicator of the wind tower under earthquake action conditions The calculation formula is: .

[0013] Preferably, in step S6, the failure risk of a certain limit state indicator The calculation formula is:

[0014] Total failure risk of wind towers The calculation formula is: .

[0015] Preferably, in step S7, the method for calculating the optimal values ​​of the parameters of the wind tower damping device is: after installing different damper devices on the wind tower, repeat steps S1-S6 to calculate the total failure risk considering the damper; by adjusting the damper parameters, a fitting curve and a fitting formula of the total failure risk with respect to each damper parameter are obtained, and the lowest value of the total failure risk is obtained from the fitting curve and the fitting formula, at which time the corresponding parameters of the damping device are the optimal values.

[0016] The present invention provides a method for evaluating the vibration reduction benefits of a wind tower damper over its entire life cycle. The method accurately calculates the vibration risk of a wind tower under the coupling of multiple physical fields and the load reduction and vibration suppression efficiency of the damper to evaluate the vibration reduction benefits of the wind tower damper over its entire life cycle. The present invention reasonably adjusts the control strategy, number of dampers, damper parameters, and damper layout position of the wind tower damper to ensure that the wind tower structure meets the requirements of reducing vibration risks and improving load reduction and vibration suppression efficiency while increasing the annual power generation time of the wind tower. The method of the present invention can enable the wind tower to achieve the good effects of safety and reliability, load reduction and vibration suppression, and cost reduction and efficiency improvement at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 The present invention is a flowchart of a method for evaluating the vibration reduction benefits of a wind tower damper over its entire life cycle. DETAILED DESCRIPTION

[0018] The present invention is further explained below in conjunction with the accompanying drawings and specific embodiments.

[0019] Example 1 Taking a wind turbine tower as an example, the present invention provides a method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle. The process is as follows: Figure 1 As shown, the specific steps include: S1. Determine the basic parameters required for evaluating the wind tower structure, including tower material, hub height, tower cross-sectional shape and size, weight and eccentricity of wind turbine impeller and nacelle, etc., so as to realize multi-field coupled integrated dynamic modeling and finite element method modeling of the wind tower structure, and then provide means for calculating the dynamic response and load effect of the wind tower structure under various working conditions.

[0020] S2. Determine the indicators and limits of the normal use limit state and bearing capacity limit state of the wind tower. The working conditions of the wind tower are divided into four categories: normal conditions, extreme conditions, earthquake conditions and fatigue conditions. The normal operating load under normal conditions is the most unfavorable load during the normal operation of the wind turbine; the extreme load under extreme conditions is the most unfavorable load in all design load conditions (DLC) except transportation and installation; the fatigue load under fatigue conditions is the total load effect of the DLC corresponding to all fatigue limit states during the entire life cycle of the wind tower. The calculation is carried out specifically according to the DLC specified in the IEC61400-1 specification.

[0021] Under normal working conditions, the serviceability limit state of the structure is verified. The serviceability limit state index is mainly the acceleration of the tower top cabin. , Tower top nacelle displacement , Structural mud surface corner , mud surface displacement Etc. Under extreme working conditions, earthquake conditions or fatigue conditions, the ultimate bearing capacity of the structure is verified. The ultimate bearing capacity index is mainly the design stress ratio of the structural member. , node stress , fatigue stress amplitude For the above-mentioned limit state indicators, the corresponding limit values ​​are determined one by one according to the design requirements, that is, , , , , , , wait.

[0022] S3. Determine the corresponding cost loss after various limit state indicators are exceeded .

[0023] (1) When the normal use limit state indicators are exceeded, the power generation efficiency will be reduced and the maintenance cost will increase. The cost loss that the wind turbine may incur is: ; Including: Single maintenance fee A =a 1 ¥ (labor and material costs); Interim measures fee B = b 1 ¥ (transportation fee) + b 2 ¥(Safety Assessment Fee)+ b 3 ¥ (facility rental fee); Downtime loss fee C = c 1 kW (generated power per year) × c 2 % (power generation reduction rate) × c 3 days (maintenance downtime days) ÷ 365 × c 4 ¥ / kW (revenue per kilowatt of electricity generated).

[0024] (2) When the bearing capacity limit state index is exceeded, it will cause the wind tower to collapse. The possible cost loss of the wind turbine is: ; Including: Machine fee D = d 1 % (whole machine discount rate) × d 2 h (Actual power generation time in one year) × d 3 kW (generated power per year) × d 4 % (power generation reduction rate) × d 5 years (length of service) × d 6 ¥ / kW h (levelized cost of electricity); Cleaning Fee E = e 1 ¥ (processing fee) + e 2 ¥ (transportation fee) + e 3 ¥(driver salary)+ e 4 ¥ (hoisting fee) + e 5 ¥ (equipment rental fee) + e 61 t (blade mass) × e 62 % (landfill reduction rate) × e 63¥ / t (landfill fee per ton); Downtime loss fee F = f 1 kW (generated power per year) × f 2 % (power generation reduction rate) × f 3 days (repair and reconstruction days) ÷ 365× f 4 ¥ / kW (revenue per kilowatt of electricity generated); Recycling Fee G = g 1 t (total mass of tower and generator) × g 2 % (recycling reduction rate) × g 3 ¥ / t (recycling price per ton of cast iron).

[0025] S4. Obtain the probability distribution of environmental parameters at the location of the wind tower , including wind speed probability distribution and earthquake intensity probability distribution; 1) For wind load, first use formula (1) to calculate the wind speed Corresponding wind pressure : (1); In the formula is the air density; Secondly, use formula (2) to calculate the wind pressure Return period , that is, wind speed Return period : (2); In the formula and It is the basic wind pressure under the 10-year and 100-year return periods in different regions of my country, which can be selected according to GB50009-2012; Then use formula (3) to calculate the design life of the support structure Wind speed (year) The probability of exceeding : (3); In the formula is the strength index; Finally, forward difference is performed on equation (3) to obtain the discrete wind speed The probability of wind load occurrence ,Right now: (4); in, is the dimension of load intensity, indicating the A load intensity.

[0026] (2) For earthquake action, first use equation (5) to calculate the design life Earthquake intensity within the year The exceedance probability of: (5); In the formula is the upper limit of earthquake intensity, which is 12; The most common earthquake intensity (frequent earthquake intensity) is the basic intensity (seismic fortification intensity) minus 1.55 degrees; for 50 years; are the shape parameters of different regions, which can be calculated based on the 50-year probability of occurrence of the seismic fortification intensity of 10%; Secondly, forward difference is performed on equation (6) to obtain the discrete earthquake intensity The probability of earthquake action : (6).

[0027] S5. Conduct a vulnerability analysis of the wind turbine tower. For a certain limit state indicator of the wind turbine tower:

[0028] (1) When considering wind turbine loads, first obtain the probability distribution of the turbulence model (e.g., NTM, ETM, EWM, etc. based on IEC61400-1 specification) at each wind speed. This will determine the probability distribution of different design load conditions (DLC) of the wind turbine at that wind speed. The probability distribution of wind speed obtained by S4 can be used to obtain the probability of occurrence of each wind speed. Multiply the wind speed occurrence probability by the DLC occurrence probability (turbulence model probability) of the wind turbine at that wind speed to obtain the occurrence probability of each DLC at each wind speed of the wind turbine. : (7); in, is the dimension of the design load condition, indicating the design load case.

[0029] On this basis, for each DLC at each wind speed, several random winds considering the rotor azimuth are selected according to the IEC61400-1 specification to perform dynamic time history analysis, and the probability of exceeding a certain limit state indicator of the wind tower under the DLC condition is calculated. ; Finally, by accumulating the product of the probability of occurrence of each DLC and the probability of exceeding the corresponding wind turbine tower index, the probability of exceeding a certain limit state index of the wind turbine tower under full load conditions is obtained. : (8); in, is a certain limit state indicator.

[0030] (2) When considering earthquake effects, first, the target response spectrum of the location of the wind turbine tower is obtained according to the GB50011-2010 standard, and several earthquake motion records are selected from the PEER database; secondly, the peak acceleration (PGA) amplitude modulation is performed on the intensity of the selected earthquake motion records; then, the dynamic time history analysis of the wind turbine tower based on the several earthquake motions after amplitude modulation is performed to obtain the IDA curve cluster; then, the exceedance probability of a certain limit state indicator of the wind turbine tower under a certain PGA is calculated. , obtain the probability of earthquake intensity at the location of the wind tower through S4 Finally, the probability of exceeding a certain limit state indicator of the wind tower under earthquake conditions is obtained by multiplying the probability of earthquake occurrence and the probability of exceeding a certain limit state indicator of the corresponding wind tower. : (9) .

[0031] S6. Calculate the failure risk of wind turbine towers. The exceedance probability of a certain limit state indicator obtained in S5 is The corresponding cost loss after the limit state index obtained by S2 is exceeded Multiply and calculate the failure risk of a certain limit state indicator ; Sum the failure risks of all limit state indicators to obtain the total failure risk of the wind tower : (10); (11).

[0032] S7. Select the optimal parameters of the wind tower damper. After installing different damper devices on the wind tower, repeat the calculation process of S1-S6 above to calculate the total failure risk of the damper. By adjusting the damper parameters (such as control strategy, damper type, damper parameters, number of dampers and damper layout, etc.), the total failure risk is obtained. The fitting curves and fitting formulas of the damper parameters are obtained from the fitting curves and fitting formulas. The lowest value, at which point the corresponding damper parameters are the optimal values.

[0033] S8. Calculate the load reduction and vibration suppression efficiency of the limit state index. Based on the optimal value of the damper parameters, use formula (12) to calculate the average load reduction and vibration suppression efficiency of each limit state index of the wind tower support after installing the damper: :

[0034] ; (12) In the formula and They are the characteristic values ​​of the limit state indicators of the wind turbine tower without dampers and with dampers installed, respectively. For the normal use limit state indicator, its characteristic value is the root mean square of the time history calculation result, and for the bearing capacity limit state indicator, its characteristic value is the maximum value.

[0035] S9. Calculate the importance coefficient of the limit state index. According to the risk of the limit state corresponding to the limit state index Total failure risk The weight coefficient of each limit state index is calculated based on the ratio of : (13);

[0036] S10. Calculate the total load reduction and vibration suppression efficiency of the wind tower after installing the damper : .

Claims

1. A method for evaluating the vibration reduction benefits of a wind tower damping and vibration reduction device over its entire life cycle, characterized in that: The following steps are involved: S1. Determine the basic parameters of the wind tower structure to be evaluated, including tower material, hub height, tower cross-sectional shape and size, weight and eccentricity of the wind turbine impeller and nacelle, and calculate the dynamic response and load effect of the wind tower structure under various working conditions; S2. Determine the indicators and limits of the normal use limit state and the bearing capacity limit state of the wind power tower. For various limit state indicators, determine the corresponding limit values ​​one by one according to the design requirements; S3. Determine the corresponding cost loss after various limit state indicators are exceeded; S4. Obtaining the probability distribution of environmental parameters at the location of the wind tower, including the probability distribution of wind speed and the probability distribution of earthquake intensity; S5. Conduct a vulnerability analysis of the wind turbine tower, and calculate the exceedance probability of a certain limit state index of the wind turbine tower under wind load conditions and the exceedance probability of a certain limit state index of the wind turbine tower under earthquake conditions for a certain limit state index of the wind turbine tower; S6. Calculate the failure risk of wind towers, including the failure risk of the limit state corresponding to the limit state index and total failure risk ; S7, calculating the optimal values ​​of the parameters of the wind tower damping and vibration reduction device; S8. Calculate the average load reduction and vibration suppression efficiency of each limit state index of the wind tower support after installing the damping vibration reduction device based on the optimal value of the damper parameters. : ； In the formula, and are the characteristic values ​​of the wind turbine tower limit state indicators before and after the damping and vibration reduction device is installed. For the normal use limit state indicator, its characteristic value is the root mean square of the time history calculation result, and for the bearing capacity limit state indicator, its characteristic value is the maximum value; is a certain limit state indicator; S9. Failure risk of limit state corresponding to limit state index Total failure risk The weight coefficient of each limit state index is calculated based on the ratio of ; S10. Calculate the total load reduction and vibration suppression efficiency of the wind tower after installing the damping vibration reduction device The formula is: 。 2. The method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle according to claim 1, characterized in that: In the step S2, calculations are performed according to the design load conditions specified in the IEC61400-1 specification; the normal use limit state of the structure is verified under normal working conditions, and the normal use limit state indicators mainly include the tower top cabin acceleration, the tower top cabin displacement, the structural mud surface angle, and the mud surface displacement; the structural bearing capacity limit state is verified under extreme working conditions, earthquake conditions, or fatigue conditions, and the bearing capacity limit state indicators mainly include the structural member design stress ratio, node stress, and fatigue stress amplitude.

3. The method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle according to claim 1, characterized in that: In step S3, when the normal use limit state indicator is exceeded, the power generation efficiency will be reduced and the maintenance cost will increase. The cost loss that may be incurred by the wind turbine is: Cost m = Single maintenance fee + temporary measures fee + downtime loss fee = a 1+ b 1+ b 2+ b 3+ c 1× c 2%× c 3÷365× c 4; in, a 1 is labor and material costs; b 1 is the transportation fee; b 2 is the safety assessment fee; b 3 is the facility rental fee; c 1 is the annual power generation, kW; c 2% is the power generation discount rate; c 3 is the number of days of downtime for maintenance; c 4 is the revenue per kilowatt of electricity generated, ¥ / kW.

4. The method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle according to claim 1, characterized in that: In step S3, when the bearing capacity limit state index is exceeded, the wind tower collapse cost will be caused, and the possible cost loss of the wind turbine is: Cost m = Equipment fee + Cleaning fee + Downtime loss fee - Recycling fee = d 1%× d 2× d 3× d 4%× d 5× d 6+( e 1+ e 2+ e 3+ e 4+ e 5+ e 61 × e 62 %× e 63 )+ f 1× f 2%× f 3÷365× f 4- g 1× g 2%× g 3; in, d 1% is the whole machine discount rate; d 2 is the actual power generation duration in one year, h; d 3 is the annual power generation, kW; d 4% is the power generation discount rate; d 5. Length of service, years; d 6 is the levelized cost of electricity, ¥ / kW·h; e 1 is the processing fee; e 2 is the transportation fee; e 3 for driver salary; e 4 is the lifting fee; e 5. Equipment rental fee; e 61 is the blade mass, t; e 62 % is the landfill reduction rate; e 63 % is the landfill fee per ton, ¥ / t; f 1 is the annual power generation, kW; f 2% is the power generation discount rate; f 4 is the revenue per kilowatt of electricity generated, ¥ / kW; f 3 is the number of days for repair and reconstruction; g 1 is the total mass of the tower and generator, t; g 2% is the recycling reduction rate; g 3 is the recycling price of cast iron per ton, ¥ / t.

5. The method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle according to claim 1, characterized in that: In step S4, for wind load, the wind speed is calculated. Corresponding wind pressure : ； In the formula is the air density; Calculate wind pressure Return period , that is, wind speed Return period : ； （2） In the formula and is the basic wind pressure at 10- and 100-year return periods in different regions; Calculate the design life of the supporting structure Internal wind speed The probability of exceeding : ； In the formula is the strength index; Perform forward difference to obtain discrete wind speed The probability of wind load occurrence ,Right now: ； in, is the dimension of load intensity, indicating the load intensity; (2) For earthquake action, calculate the design life Internal earthquake intensity The exceedance probability of: ； In the formula is the upper limit of earthquake intensity; For the moderate intensity, take the basic intensity minus 1.55 degrees; for 50 years; are the shape parameters of different regions, which are calculated based on the 50-year probability of occurrence of the seismic fortification intensity of 10%; Perform forward difference to obtain discrete earthquake intensity The probability of earthquake action : 。 6. The method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle according to claim 5, characterized in that: In step S5, when considering the wind turbine load, the probability of occurrence of each design load condition is The calculation formula is: ； in, is the dimension of the design load condition, indicating the design load cases; The probability of exceeding a certain limit state indicator of a wind turbine tower under full load conditions The calculation formula is: ； When considering earthquake action, the probability of exceeding a certain limit state indicator of the wind tower under earthquake action conditions The calculation formula is: 。 7. The method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle according to claim 6, characterized in that: In step S6, the failure risk of a certain limit state indicator The calculation formula is: ； Total failure risk of wind towers The calculation formula is: 。 8. The method for evaluating the vibration reduction benefits of a wind turbine tower damping and vibration reduction device over its entire life cycle according to any one of claims 1 to 7, characterized in that: In step S7, the method for calculating the optimal values ​​of the parameters of the wind tower damping device is as follows: after installing different damper devices on the wind tower, repeat steps S1-S6 to calculate the total failure risk of the damper; by adjusting the damper parameters, a fitting curve and a fitting formula of the total failure risk with respect to each damper parameter are obtained, and the lowest value of the total failure risk is obtained from the fitting curve and the fitting formula, at which time the corresponding parameters of the damping device are the optimal values.

Images