阅读说明：本技术 一种塔架极限载荷的估算方法、电子设备及存储介质 (Estimation method of tower ultimate load, electronic equipment and storage medium ) 是由 周颖 李俊 于 2021-09-15 设计创作，主要内容包括：本发明公开了一种塔架极限载荷的估算方法、电子设备及存储介质,所述方法包括：获取风机塔架的最大载荷,所述风机塔架的最大载荷包括所述风机塔架前后方向的最大载荷和/或左右方向的最大载荷；根据风机塔架所受最大载荷和风机塔架在预设期限内的极限载荷之间的关系估算出所述风机塔架的极限载荷值。本发明可以快速评估出风机塔架的极限载荷,从而便于校核材料的安全系数,保证风机安全运行。(The invention discloses a tower ultimate load estimation method, electronic equipment and a storage medium, wherein the method comprises the following steps: acquiring the maximum load of a fan tower, wherein the maximum load of the fan tower comprises the maximum load of the fan tower in the front-back direction and/or the maximum load of the fan tower in the left-right direction; and estimating the limit load value of the fan tower according to the relation between the maximum load borne by the fan tower and the limit load of the fan tower in a preset time limit. The method can quickly evaluate the limit load of the tower frame of the air outlet fan, thereby facilitating the check of the safety coefficient of materials and ensuring the safe operation of the fan.)

1. A method for estimating tower ultimate load, comprising:

acquiring the maximum load of a fan tower, wherein the maximum load of the fan tower comprises the maximum load of the fan tower in the front-back direction and/or the maximum load of the fan tower in the left-right direction;

and estimating the limit load value of the fan tower according to the relation between the maximum load borne by the fan tower and the limit load of the fan tower in a preset time limit.

2. The method for estimating tower limit load according to claim 1, wherein the maximum load in the fore-and-aft direction of the wind turbine tower is calculated based on the gravitational moment, the average moment and its fluctuation component applied to the wind turbine tower, and the moment and its fluctuation component applied to the wind turbine hub by wind shear.

3. Method for estimating tower limit load according to claim 2, wherein the maximum load M in the fore-aft direction of the wind turbine tower is the maximum load M_{fro_max}The following formula is used for calculation:

M_{fro_max}＝M_G+M_total+M_shear+M_{vol_total}+M_{vol_shear}

in the formula, M_GRepresenting moment of gravity action, M_totalRepresenting the average moment borne by the fan tower; m_shearRepresenting the acting moment of wind shear on the fan hub; m_{vol_total}Representing the fluctuation component of the average moment borne by the wind turbine tower; m_{vol_shear}Representing the fluctuating component of the wind shear's moment of action on the fan hub.

4. Method for estimating a tower limit load according to claim 3, wherein said gravitational moment M_GCalculated using the following formula:

M_G＝M_{G_nac}+M_{G_hub}+M_{G_b}

in the formula, M_{G_nac}Representing the initial applied moment of gravity, M, of the fan hub to the tower_{G_hub}Representing the initial applied moment of gravity, M, of the wind turbine blade against the tower_{G_b}The initial acting gravity moment of the wind turbine nacelle on the tower is represented.

5. Method for estimating the tower limit load according to claim 4, wherein said wind turbine tower is subjected to an average moment M_totalCalculated using the following formula:

M_total＝M_thr+M_nac+M_tow

in the formula, M_thrIndicating the thrust moment, M, of the fan wheel_nacRepresenting the moment of action of the nacelle of the wind turbine, M_towRepresenting the wind turbine tower moment.

6. Method for estimating tower limit load according to claim 5, wherein said wind shear exerts a moment M on the hub of the wind turbine_shearCalculated using the following formula:

M_shear＝f(u,α,D)

where u represents the annual average wind speed, α represents the wind shear, D represents the fan wheel diameter, and f represents M_shearFunctional relation with u, α, D.

7. Method for estimating the tower limit load according to claim 6, wherein the wind turbine tower is subjected to a fluctuating component M of the mean moment_{vol_total}Calculated using the following formula:

M_{vol_total}＝l_tσ_t

in the formula I_tRepresenting the average moment borne by the fan tower; sigma_tIndicating the fluctuation intensity; l_tAccording to the ratedCarrying out empirical formula value on the wind speed; sigma_tAccording to M_{RMS_total} ²＝M_total ²+σ_t ²Obtaining M_{RMS_total}Represents M_totalRoot mean square (rms).

8. The method for estimating tower ultimate load according to claim 7,

the fluctuation component M of the acting moment of the wind shear on the fan hub_{vol_shear}The following formula is used for calculation:

M_{vol_shear}＝l_sσ_s

in the formula I_sRepresenting a given crest factor; sigma_sThe standard deviation of the average moment borne by the fan tower is represented; wherein σ_sAccording to M_{RMS_shear} ²＝M_shear ²+σ_s ²Performing a calculation of M_{RMS_shear}F (u, α, I, D), I denotes the turbulence intensity, f (u, α, I, D) denotes M_{RMS_shear}Functional relation with u, alpha, I, D.

9. The method for estimating tower limit load of claim 8, wherein said obtaining of maximum load M in left and right direction of wind turbine tower_{side_max}The following formula is used for calculation:

M_{side_max}＝f(P)

wherein f (P) represents M_{side_max}And a functional relationship with P, wherein P represents the runner power.

10. The method for estimating tower limit load according to claim 9, wherein the relationship between the maximum load of the wind turbine tower and the limit load of the wind turbine tower in the preset time frame is as follows:

M_n＝kM_max

in the formula, M_nRepresenting the limit load value of the wind turbine tower in a preset time limit, wherein n is a positive integer and represents the preset time limit, and M_maxRepresenting the maximum load value borne by the fan tower; k represents outsideAnd (5) pushing the coefficient.

11. The method for estimating tower limit load of claim 10, wherein said wind turbine tower limit load comprises a wind turbine tower top limit load, and said wind turbine tower top limit load value M_{tow_top_max}Calculated using the following formula:

12. the method for estimating tower ultimate load of claim 11, wherein the wind turbine tower ultimate load comprises wind turbine tower bottom ultimate load,

ultimate load value M of bottom of fan tower_{tow_bot_max}Calculated using the following formula:

M_{tow_bot_max}＝M_{pro_max}。

13. an electronic device comprising a processor and a memory, the memory having stored thereon a computer program which, when executed by the processor, implements the method of any of claims 1 to 11.

14. A readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the method of any one of claims 1 to 11.

Technical Field

The invention relates to the technical field of wind power generation, in particular to a tower ultimate load estimation method, electronic equipment and a storage medium.

Background

In the design of a wind generating set (called a fan for short) and the construction stage of a wind power plant, load calculation of the wind generating set is a very important link. The load of the wind turbine refers to the forces or moments exerted on its components by the external environment and the internal environment, including aerodynamic loads, gravitational loads, inertial loads, and operational loads due to the action of the control system. The load of the wind generating set can be further divided into limit load and fatigue load according to the structural design requirement. The limit load refers to the maximum load which can be borne by the wind generating set, and the fatigue load refers to the alternating cycle load acting on the wind generating set. The purpose of the load analysis and calculation of the wind generating set is to calculate the stress and the strain of a structure under a specific working condition, and further carry out ultimate strength check and fatigue strength check.

The limit load affects the safety of the fan, and the load of the fan fluctuates continuously in each operating condition. The maximum value of the random load fluctuation which can be borne by the fan structure determines the safety of the fan in operation. The fan tower in the fan is loaded by itself during operation and also under the action of a rotating engine room and the like, and the stress is complex, so that the limit load of the fan tower needs to be determined, the safety coefficient of materials is checked, and the safe operation of the fan is ensured.

Disclosure of Invention

The invention aims to provide an estimation method of a tower frame limit load, electronic equipment and a storage medium, so as to realize the purpose of quickly estimating the limit load of a fan tower frame, thereby facilitating the check of the safety coefficient of materials and ensuring the safe operation of a fan.

In order to achieve the above purpose, the invention is realized by the following technical scheme:

a method of estimating tower ultimate loads, comprising: acquiring the maximum load of a fan tower; the maximum load of the fan tower comprises the maximum load of the fan tower in the front-back direction and/or the maximum load of the fan tower in the left-right direction; and estimating the limit load value of the fan tower according to the relation between the maximum load borne by the fan tower and the limit load of the fan tower in a preset time limit.

Optionally, the maximum load of the wind turbine tower in the front-rear direction is calculated according to the gravity acting moment, the average moment and the fluctuation component of the wind turbine tower, and the acting moment and the fluctuation component of the wind shear on the wind turbine hub.

Optionally, a maximum load M of the wind turbine tower in a fore-and-aft direction_{fro_max}The following formula is adopted for calculation：

M_{fro_max}＝M_G+M_total+M_shear+M_{vol_total}+M_{vol_shear}

In the formula, M_GRepresenting moment of gravity action, M_totalRepresenting the average moment borne by the fan tower; m_shearRepresenting the acting moment of wind shear on the fan hub; m_{vol_total}Representing the fluctuation component of the average moment borne by the wind turbine tower; m_{vol_shear}Representing the fluctuating component of the wind shear's moment of action on the fan hub.

Optionally, the gravitational moment M_GCalculated using the following formula:

M_G＝M_{G_nac}+M_{G_hub}+M_{G_b}

in the formula, M_{G_nac}Representing the initial applied moment of gravity, M, of the fan hub to the tower_{G_hub}Representing the initial applied moment of gravity, M, of the wind turbine blade against the tower_{G_b}The initial acting gravity moment of the wind turbine nacelle on the tower is represented.

Optionally, the average moment M experienced by the wind turbine tower_totalCalculated using the following formula:

M_total＝M_thr+M_nac+M_tow

in the formula, M_thrIndicating the thrust moment, M, of the fan wheel_nacRepresenting the moment of action of the nacelle of the wind turbine, M_towRepresenting the wind turbine tower moment.

Optionally, the wind shear acts on the wind turbine hub with a moment M_shearCalculated using the following formula:

M_shear＝f(u,α,D)

where u represents the annual average wind speed, α represents the wind shear, D represents the fan wheel diameter, and f represents M_shearFunctional with u, α and D.

Optionally, the wind turbine tower is subjected to a fluctuating component M of the mean moment_{vol_total}Is calculated by the following formula

Calculating:

M_{vol_total}＝l_tσ_t

in the formula I_tRepresenting a given crest factor; sigma_tRepresents M_totalStandard deviation of mean moment; l_tCarrying out value taking by an empirical formula according to rated wind speed; sigma_tAccording to M_{RMS_total} ²＝M_total ²+σ_t ²Obtaining M_{RMS_total}Represents M_shearRoot mean square (rms).

Optionally, the wind shear has a fluctuating component M of the moment of action on the fan hub_{vol_shear}The following formula is used for calculation:

M_{vol_shear}＝l_sσ_s

in the formula I_sRepresenting a given crest factor; sigma_sThe standard deviation of the average moment borne by the fan tower is represented; wherein σ_sAccording to M_{RMS_shear} ²＝M_shear ²+σ_s ²Performing a calculation of M_{RMS_shear}F (u, α, I, D), I denotes the turbulence intensity, f (u, α, I, D) denotes M_{RMS_shear}And the function relation of the wind shear alpha, the turbulence intensity I, the annual average wind speed u and the diameter D of the runner.

Optionally, the maximum load M in the left-right direction of the wind turbine tower is obtained_{side_max}The following formula is used for calculation:

M_{side_max}＝f(P)

wherein f (P) represents M_{side_max}And the functional relation exists with P, wherein P represents the power of the wind wheel.

Optionally, the relationship between the maximum load applied to the wind turbine tower and the limit load of the wind turbine tower in the preset time limit is as follows:

M_n＝kM_max

in the formula, M_nRepresenting the limit load value of the wind turbine tower in a preset time limit, wherein n is a positive integer and represents the preset time limit, and M_maxRepresenting the maximum load value borne by the fan tower; k represents the extrapolation coefficient.

Optionally, the extreme load of the wind turbine tower comprises the extreme load of the top of the wind turbine tower comprises the wind turbine towerUltimate load at the top, ultimate load value M at the top of the wind turbine tower_{tow_top_max}Calculated using the following formula:

optionally, the limit load of the wind turbine tower comprises a limit load of the bottom of the wind turbine tower, and a limit load value M of the bottom of the wind turbine tower_{tow_bot_max}Calculated using the following formula:

M_{tow_bot_max}＝M_{pro_max}。

in another aspect, the present invention also provides an electronic device comprising a processor and a memory, the memory having stored thereon a computer program which, when executed by the processor, implements the method as described above.

In yet another aspect, the present invention also provides a readable storage medium having stored therein a computer program which, when executed by a processor, implements a method as described above.

The invention has at least one of the following advantages:

the invention provides a new method for quickly estimating the limit load of a tower, which is characterized in that the limit load value is obtained by searching the relation between the limit load value and the maximum load value and extrapolating, and the limit load value is decomposed into the sum of an average component and a fluctuation component. And the fluctuation components of different extreme values adopt different peak value factors to display extreme value differences, the transmission values of all key parameters in the calculation process are represented by basic parameters, and estimated values of the extreme loads of the wind turbine tower at different heights are given.

The method reduces the time consumed by load time-sequence calculation, is convenient for checking the safety coefficient of materials, ensures the safe operation of the fan, and is favorable for quickly carrying out the initial design of the fan tower.

Drawings

FIG. 1 is a schematic flow chart of a method for estimating a tower ultimate load according to an embodiment of the present invention;

FIG. 2 is a schematic view of a tower coordinate system provided by an embodiment of the present invention;

FIG. 3 is a schematic view of the center of gravity and wind shear of various components of a wind turbine generator system according to an embodiment of the present invention;

FIG. 4 is a schematic diagram illustrating a comparison of tower top loads of a wind turbine tower according to an embodiment of the present invention;

fig. 5 is a schematic diagram illustrating a comparison of tower bottom loads of a wind turbine tower according to an embodiment of the present invention.

Detailed Description

The method for estimating tower ultimate load, the electronic device and the storage medium according to the present invention will be described in further detail with reference to the accompanying drawings and the following detailed description. The advantages and features of the present invention will become more apparent from the following description. It is to be noted that the drawings are in a very simplified form and are all used in a non-precise scale for the purpose of facilitating and distinctly aiding in the description of the embodiments of the present invention. To make the objects, features and advantages of the present invention comprehensible, reference is made to the accompanying drawings. It should be understood that the structures, ratios, sizes, and the like shown in the drawings and described in the specification are only used for matching with the disclosure of the specification, so as to be understood and read by those skilled in the art, and are not used to limit the implementation conditions of the present invention, so that the present invention has no technical significance, and any structural modification, ratio relationship change or size adjustment should still fall within the scope of the present invention without affecting the efficacy and the achievable purpose of the present invention.

As shown in fig. 1, the present embodiment provides a method for estimating a tower ultimate load, including:

step S100: the method comprises the steps of obtaining the maximum load of the fan tower, wherein the maximum load of the fan tower comprises the maximum load of the fan tower in the front-back direction and/or the maximum load of the fan tower in the left-right direction.

Step S200: and estimating the limit load value of the fan tower according to the relation between the maximum load borne by the fan tower and the limit load of the fan tower in a preset time limit.

In the embodiment, the factors influencing the tower top limit load of the fan tower and the tower bottom limit load of the fan tower are analyzed; the maximum load of the wind turbine tower is mainly divided into the maximum load in the front-back direction and the maximum load in the left-right direction, in this embodiment, the wind wheel faces the wind direction in the front-back direction, and the wind wheel faces the left-right direction in the vertical direction. As shown in fig. 2, the coordinates of the hub rotation coordinate system, the hub fixed coordinate system, the tower top coordinate system and the tower bottom coordinate system are respectively O1X1Y1Z1, O2X2Y2Z2, O3X3Y3Z3 and O4X4Y4Z4 coordinate systems, and the equivalent fatigue bending moment and force around each axis are determined by a right-hand rule.

As can be seen by analysis, the limit load at the top of the tower is mainly influenced by the cut-out wind speed, and the limit load at the bottom of the tower is mainly influenced by the rated wind speed.

According to the embodiment, the maximum load of the fan tower in the front-back direction can be calculated according to the action moment of the gravity, the average moment and the fluctuation component of the average moment borne by the fan tower, the action moment of wind shear on the fan hub and the fluctuation component of the action moment.

Specifically, the step S100 specifically includes: the maximum load M of the wind turbine tower in the front-back direction_{fro_max}The following formula is used for calculation:

M_{fro_max}＝M_G+M_total+M_shear+M_{vol_total}+M_{vol_shear} (1)

in the formula, M_GRepresenting moment of gravity action, M_totalRepresenting the average moment borne by the fan tower; m_shearRepresenting the acting moment of wind shear on the fan hub; m_{vol_total}Representing the fluctuation component of the average moment borne by the wind turbine tower; m_{vol_shear}Representing the fluctuating component of the wind shear's moment of action on the fan hub.

As shown in fig. 3, fig. 3 shows the initial gravitational moment mainly considering the nacelle of the wind turbine, the hub and the blades, and the gravity center of the blades is selected to be 1/3 from the blade root for the gravity center and the distance from the center of the tower top of each component (the blades, the hub and the nacelle) in the wind turbine generator set. In the calculation, the existence of the elevation angle cone angle of the wind wheel is considered, and the initial action gravity moment and the initial weight of the hub, the blade and the engine room on the tower are calculated according to the action of force and the arm of forceThe moment only affects the load mean value (average component) and has no influence of fluctuation component; therefore, the influence of the initial gravitational moment on the tower load is increased in the estimation of the tower top limit load. The gravity acting moment M_GCalculated using the following formula:

M_G＝M_{G_nac}+M_{G_hub}+M_{G_b} (2)

in the formula, M_{G_nac}Representing the initial applied moment of gravity, M, of the fan hub to the tower_{G_hub}Representing the initial applied moment of gravity, M, of the wind turbine blade against the tower_{G_b}The initial acting gravity moment of the wind turbine nacelle on the tower is represented.

Considering the influence of the wind wheel thrust, the nacelle moment and the tower moment on the tower load: according to different heights of the tower, the stress of the tower can be integrated and calculated to obtain the average moment applied to the tower, and the average moment applied to the tower fluctuates.

The average moment M borne by the fan tower_totalCalculated using the following formula:

M_total＝M_thr+M_nac+M_tow (3)

in the formula, M_thrIndicating the thrust moment, M, of the fan wheel_nacRepresenting the moment of action of the nacelle of the wind turbine, M_towRepresenting the wind turbine tower moment.

When the wind wheel of the fan works, the maximum fluctuation component exists in the torque response borne by the fan towers in the front and back directions. Mainly the influence of wind on the tower stress, i.e. the wind wheel thrust, the engine room moment, M brought by the tower moment_total(average moment M to which the wind turbine tower is subjected_total) The maximum load taking into account the fluctuating component M of this moment_{vol_total}Wave component M of_{vol_total}(ii) a The embodiment thus requires the maximum load on the wind turbine tower to take into account the fluctuating component of this moment.

The fluctuation component of the moment is related to gust, and the fluctuation component M of the average moment borne by the wind turbine tower can be obtained according to the derivation principle_{vol_total}Calculated using the following formula:

M_{vol_total}＝l_tσ_t (4)

in the formula I_tA peak coefficient representing a wind turbine thrust fluctuation; sigma_tRepresenting the deviation of the wind wheel thrust; l_tCarrying out value taking of an empirical formula according to a rated wind speed, specifically, carrying out calculation value taking of the empirical formula according to the fact that the wind speed at a hub is higher than the rated wind speed or lower than the rated wind speed, and determining the peak coefficient of wind wheel thrust fluctuation; sigma_tAccording to M_{RMS_total} ²＝M_total ²+σ_t ²Obtaining M_{RMS_total}The root mean square is indicated.

Above l_tThe values can be taken with reference to an empirical formula, based on the wind speed at the hub, and the following is a list of relevant variables involved in the empirical formula.

There are coefficients such that:

M_{max_total}＝A*M_total (5)

in the formula, M_{max_total}Representing the maximum value of the rotor thrust and a represents a parameter.

Fitting the coefficient, A and the cut-in wind speed by sample data (pre-acquired calculation data of a large number of fans in actual design and operation); cutting out wind speed; wind speed at the hub; there is a specific functional relationship for the rated wind speed, namely:

A＝f(v_in,v_out,v_h,v_r) (6)

wherein f (v)_in,v_out,v_h,v_r) Denotes A and v_inCutting in wind speed; v. of_outCutting out wind speed; wind speed at the hub; v. of_rA functional relationship of rated wind speed; v. of_inRepresenting a cut-in wind speed; v. of_outRepresenting the cut-out wind speed; v. of_hRepresenting the wind speed at the hub; v. of_rIndicating the rated wind speed.

Considering the wind shear effect on the maximum load of the tower: analyzing the moment of the tower in the front-back direction: when wind blows towards the wind wheel, the wind acts on the wind wheel, and the existence of wind shear generates moment on the hub of the wind wheel, so that the moment of the wind shear acting on the hub is increased in the estimation of the ultimate load of the tower, and the acting moment has fluctuation components.

Load limit M of rotor moment in runner hub_{max_shear}Is defined as:

M_{max_shear}＝M_shear+M_{vol_shear} (7)

in the formula, M_shearRepresenting the moment of wind shear acting on the fan hub (the average load component); m_{vol_shear}Representing the fluctuating component of the wind shear's moment of action on the fan hub.

Fitting to obtain the acting moment M of the wind shear on the fan hub through sample data (pre-acquired calculated data of a large number of fans in actual design operation)_shearThe functional relationship representation is as follows:

M_shear＝f(u,α,D) (8)

where u represents the annual average wind speed, α represents the wind shear, D represents the fan wheel diameter, and f (u, α, D) represents M_shearAs a function of u, α and D.

The fluctuation component M of the acting moment of the wind shear on the fan hub_{vol_shear}The following formula is used for calculation:

M_{vol_shear}＝l_sσ_s (9)

in the formula I_sRepresenting a given crest factor; sigma_sThe standard deviation of the average moment borne by the fan tower is represented; wherein σ_sAccording to M_{RMS_shear} ²＝M_shear ²+σ_s ²Performing calculation, wherein M is obtained by fitting sample data (calculated data obtained in advance when a large number of fans are actually designed and operated)_{RMS_shear}F (u, α, I, D), I denotes the turbulence intensity, f (u, α, I, D) denotes M_{RMS_shear}And the function relation of the wind shear alpha, the turbulence intensity I, the annual average wind speed u and the diameter D of the runner.

After the fluctuation influence caused by the wind wheel thrust and the wind shear thrust is uniformly considered, the maximum load of the wind turbine tower in the front-back direction can be converted into the following formula (1):

according to the statistics of a large amount of data, carrying out maximum load M on the tower in the left and right directions_{side_max}After normalization, M_{side_max}Has a certain functional relation with the wind wheel power P (the maximum load M of the tower in the left and right directions)_{side_max}After normalization, the linear relation exists between the normalized values and the power of the wind wheel), and the values are represented by f (P), so that the left and right loads of the tower can be estimated by fitting.

The method is characterized in that the load of the left side and the right side of a large number of units and the power P of a wind wheel (unit power) are counted, and the load of the left side and the right side of a tower can be estimated by performing formula fitting on two variables.

Acquiring the maximum load M of the left and right directions of the wind turbine tower_{side_max}The following formula is used for calculation:

M_{side_max}＝f(P) (11)

wherein f (P) represents M_{side_max}And P represents the unit power. The formula (11) can be obtained by counting sample data (pre-acquired calculation data of a large number of fans in actual design and operation).

Finding a relation coefficient between the maximum load borne by the fan tower and the limit load of the fan tower in a preset time limit to obtain the relation between the maximum load borne by the fan tower and the limit load of the fan tower in the preset time limit as follows:

M_n＝kM_max (12)

in the formula, M_nRepresenting the limit load value of the wind turbine tower in a preset time limit, wherein n is a positive integer and represents the preset time limit, and M_maxRepresenting the maximum load value borne by the fan tower; k represents the extrapolation coefficient. Thus, according to the above relation, only the maximum load M is needed to obtain the limit load_maxThat is, the maximum load value M_maxDecomposed into an average component and a fluctuationAnd (4) components.

In this embodiment, n is 50, i.e., the predetermined period is 50 years.

M_max＝M_ave+M_{max_vol}＝M_ave+σ (14)

In the formula, M_aveRepresenting the mean component of the maximum load value, M_{max_vol}Represents the fluctuation component of the maximum load value, and σ represents the deviation;

M_ave＝M_G+M_total+M_shear (15)

in the formula, M_GRepresenting moment of gravity action, M_totalThe average moment (acting moment of wind on a rotating wheel, a cabin and a tower barrel) borne by the wind turbine tower is represented; m_shearRepresenting the moment of wind shear on the fan hub.

σ＝M_{vol_shear}+M_{vol_total} (19)

In the formula, M_{vol_total}Representing the fluctuating component of the mean moment to which the tower of the fan is subjected, M_{vol_shear}Representing the fluctuating component of the wind shear's moment of action on the fan hub.

The difference of the influence of the left and right tower loads and the front and rear tower loads on the tower bottom load of the tower top is integrated: for tower bottom load, left and right moments in the front and rear directions are dominant; for the tower top limit load, the influence of left and right moments in the front-back direction and the left-right direction is large, and the tower top limit load is estimated and synthesized.

This embodiment thus further comprises: the ultimate load of the fan tower comprises an ultimate load value M of the top of the fan tower_{tow_top_max}Calculated using the following formula:

the ultimate load of the fan tower comprises an ultimate load value M of the ultimate load of the bottom of the fan tower_{tow_bot_max}Calculated using the following formula:

M_{tow_bot_max}＝M_{pro_max}。

in combination with fig. 4 and 5, tower limit loads of several units are selected for estimation, and the tower bottom and the tower top are mainly concerned. And (4) comparing the actually calculated limit load values, and normalizing all the limit load values. The obtained tower bottom limit load contrast is shown in figure 4, the tower top load contrast is shown in figure 5, the calculated result values are well matched, the error is within the acceptable range of the initial design load, and the estimation method has more accurate precision.

In another aspect, the present invention also provides an electronic device comprising a processor and a memory, the memory having stored thereon a computer program which, when executed by the processor, implements the method as described above.

In yet another aspect, the present invention also provides a readable storage medium having stored therein a computer program which, when executed by a processor, implements a method as described above.

The embodiment provides a new method for quickly estimating the limit load of the tower, which is characterized in that the limit load value is obtained by searching the relation between the limit load value and the maximum load value and extrapolating, and the limit load value is decomposed into the sum of an average component and a fluctuation component. And the fluctuation components of different extreme values adopt different peak value factors to display extreme value differences, the transmission values of all key parameters in the calculation process are represented by basic parameters, and estimated values of the extreme loads of the wind turbine tower at different heights are given.

The embodiment reduces the time consumed by load time sequence calculation, is convenient for checking the safety coefficient of materials, ensures the safe operation of the fan, and is favorable for rapidly carrying out the initial design of the fan tower frame

It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.

It should be noted that the apparatuses and methods disclosed in the embodiments herein can be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments herein. In this regard, each block in the flowchart or block diagrams may represent a module, a program, or a portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.

In addition, the functional modules in the embodiments herein may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

While the present invention has been described in detail with reference to the preferred embodiments, it should be understood that the above description should not be taken as limiting the invention. Various modifications and alterations to this invention will become apparent to those skilled in the art upon reading the foregoing description. Accordingly, the scope of the invention should be determined from the following claims.