Icing fan blade wind energy utilization rate calculation method and system based on section segmentation

Abstract

The invention discloses a method and system for calculating the wind energy utilization rate of an ice-coated fan blade based on section segmentation, and relates to the field of wind energy and wind power generation, and the method comprises the steps: sectioning a wind driven generator blade, and obtaining two-dimensional section data; based on the wind speed and the operation parameters, an airflow field of the two-dimensional airfoil section is calculated, and then the lift-drag coefficient is obtained; establishing a wind energy utilization rate calculation model by utilizing the section data and the lift resistance coefficient, carrying out stress analysis on the blade element, and calculating lift resistance, thrust and torque; and the wind energy utilization rate under the normal operation condition is calculated according to the thrust and the torque. By focusing the key section, the calculation burden is remarkably reduced, and the influence of icing on the wind energy utilization rate can be accurately evaluated. Besides, the method is suitable for analysis under different icing forms and thicknesses, fan performance can be predicted in advance, performance reduction and equipment damage caused by icing are avoided, the maintenance cost is reduced, and the power generation efficiency is improved.

Description

Technical Field

[0001] The present invention relates to the field of wind energy and wind power generation, and in particular to a method and system for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation. Background Art

[0002] As a clean and renewable energy source, wind energy is gradually being accepted and widely used by more and more countries. my country's wind power industry is also growing rapidly. In recent years, the wind power industry in Sichuan, Chongqing, Yunnan and other regions has been increasing year by year. These regions are located in the southwest of my country, with a humid climate, prominent plateau and mountainous terrain features, and winter temperatures around 0°C, which provides favorable environmental conditions for icing. The problem of icing on wind turbine blades has gradually become prominent, and the average annual power generation loss caused by icing can be as high as 30%.

[0003] In order to further reveal the influence of icing on the power output characteristics of wind turbines, some scholars have experimentally studied the influence of icing environmental parameters on the power output characteristics of small wind turbines under artificial climate chambers and natural environmental conditions. The results show that the output power value of the wind turbine first decreases and then increases with the decrease of ambient temperature. The greater the liquid water content, the lower the output power of the wind turbine. As the thickness of the ice increases, the greater the wind speed, the greater the power loss of the icing wind turbine. In order to simplify the complex artificial icing process, light clay and putty tests were used to simulate the influence of different forms of icing on its power. The results show that compared with streamlined ice, angular ice will bring more aerodynamic load losses. Similarly, the prior art considers the influence of icing roughness on the aerodynamic performance of the wind turbine. The research results show that the decline in aerodynamic performance of iced blades is mainly caused by changes in the aerodynamic shape of the blades, and is less affected by the roughness of the icing. In fact, in addition to the thickness and shape of the icing, the icing range on the surface of the wind turbine blade under different environmental conditions is also quite different, especially under different median volume diameters of water droplets. As the icing develops, the icing tends to develop along the leading edge to the trailing edge of the wind turbine blade. When studying the influence of icing characteristics on the aerodynamic and power output characteristics of wind turbines, most researchers will establish a finite element calculation model of the three-dimensional flow field of wind turbine blades. However, wind turbines are large in size and have a large number of grids, which has certain requirements on computer performance and calculation time. This has to some extent restricted the application of this technology to the analysis of the power characteristics of icing wind turbines. Summary of the invention

[0004] In view of the above-mentioned problems, the present invention is proposed.

[0005] Therefore, the technical problem solved by the present invention is: how to quickly obtain the wind energy utilization rate under the condition of ice covering the blades of a wind turbine with a small amount of calculation.

[0006] In order to solve the above technical problems, the present invention provides the following technical solutions:

[0007] In a first aspect, an embodiment of the present invention provides a method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation, comprising:

[0008] Cut the wind turbine blades to obtain the two-dimensional cross-sectional data of the wind turbine blades at each position;

[0009] Based on the wind speed and wind turbine operating parameters, the airflow field of each wind turbine blade two-dimensional airfoil section is calculated to obtain the lift and drag coefficient of the two-dimensional airfoil section;

[0010] Based on the two-dimensional cross-sectional data and the lift-drag coefficient of the two-dimensional airfoil section, a calculation model for the wind energy utilization rate of the fan blade element is established, the force analysis of the two-dimensional cross-sectional blade element is carried out, and the lift-drag coefficient of the blade element is calculated;

[0011] According to the lift and resistance of the blade element, calculate the thrust and torque of the blade element;

[0012] Based on the thrust and torque acting on the blade element, the wind energy utilization rate of the blade element under normal operation is calculated.

[0013] As a preferred solution of the wind energy utilization rate calculation method of ice-covered wind turbine blades based on section segmentation, wherein:

[0014] The sectioning of the wind turbine blade comprises:

[0015] Assume that the radius of the wind turbine blade is R, and select the positions of 0.25R, 0.5R, 0.75R and R for cutting.

[0016] As a preferred solution of the wind energy utilization rate calculation method of ice-covered wind turbine blades based on section segmentation, wherein:

[0017] The calculation of blade lift resistance includes:

[0018] Based on the results of the two-dimensional cross-section blade element force analysis, the expression of the blade element lift resistance is obtained as follows:

[0019]

[0020] Where ρ is the air density in kg / m 3 ; c is the blade chord length, in m; r is the radius of the blade, in m; C L and C D are the lift coefficient of the two-dimensional airfoil section and the drag coefficient of the two-dimensional airfoil section; V rel It is the velocity of blade element relative to airflow, in m / s.

[0021] As a preferred solution of the wind energy utilization rate calculation method of ice-covered wind turbine blades based on section segmentation, wherein:

[0022] The calculation of the thrust and torque of the blade element includes:

[0023] The turning thrust dF received by the blade element t It is expressed as:

[0024]

[0025] As a preferred solution of the wind energy utilization rate calculation method of ice-covered wind turbine blades based on section segmentation, wherein:

[0026] The calculation of the thrust and torque of the blade element also includes:

[0027] Assume the number of blades is B; the blade angular velocity is Ω, in rad / s; the tip speed ratio of the blade element at r is λ r =Ω·r / V; lift-to-drag ratio e = C L / C D ; Then the torque of the micro-circular ring where the blade element is located is:

[0028]

[0029] As a preferred solution of the wind energy utilization rate calculation method of ice-covered wind turbine blades based on section segmentation, wherein:

[0030] The calculation of blade element wind energy utilization under normal operation includes:

[0031] The blade element position is selected to calculate the blade element wind energy utilization rate under normal operation. When calculating, it is assumed that the blade element relative airflow velocity and the components perpendicular and parallel to the blade rotation plane have the following relationship:

[0032]

[0033] Among them, V y and V x are the components perpendicular and parallel to the blade rotation plane respectively.

[0034] As a preferred solution of the wind energy utilization rate calculation method of ice-covered wind turbine blades based on section segmentation, wherein:

[0035] The calculation of blade element wind energy utilization under normal operation also includes:

[0036] Blade output power dP and wind energy utilization rate C p-r for:

[0037]

[0038] In a second aspect, an embodiment of the present invention provides a wind energy utilization rate calculation system for ice-covered wind turbine blades based on section segmentation, comprising:

[0039] A cross-sectional data acquisition module is used to section the wind turbine blades and obtain two-dimensional cross-sectional data of the two-dimensional wind turbine blades at each position;

[0040] A lift and drag coefficient acquisition module is used to calculate the airflow field of the two-dimensional airfoil section of each fan blade based on the wind speed and the fan operating parameters to obtain the lift and drag coefficient of the two-dimensional airfoil section;

[0041] The lift and drag calculation module is used to establish a wind energy utilization rate calculation model for the blade element of a wind turbine blade based on the two-dimensional cross-section data and the lift and drag coefficient of the two-dimensional airfoil section, perform a force analysis on the two-dimensional cross-section blade element, and calculate the lift and drag of the blade element;

[0042] A parameter calculation module, used to calculate the thrust and torque of the blade element according to the lift and drag of the blade element;

[0043] The utilization rate calculation module is used to calculate the wind energy utilization rate of the blade element under normal operation based on the thrust and torque received by the blade element.

[0044] In a third aspect, an embodiment of the present invention provides a computing device, including:

[0045] Memory and processor;

[0046] The memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions. When the one or more programs are executed by the one or more processors, the one or more processors implement the method for calculating wind energy utilization of ice-covered wind turbine blades based on profile segmentation as described in any embodiment of the present invention.

[0047] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the method for calculating wind energy utilization of ice-covered wind turbine blades based on profile segmentation.

[0048] Beneficial effects of the present invention: The present invention uses a method based on section segmentation to quickly obtain the wind energy utilization rate under ice-covered conditions of wind turbine blades with a small amount of calculation. This method avoids complex overall blade analysis and focuses on key sections, thereby significantly reducing the computational burden; through detailed two-dimensional airfoil section airflow field calculation, the lift coefficient and drag coefficient can be obtained more accurately, and then the lift and drag of the blade element can be accurately calculated. This is crucial for accurately evaluating the impact of icing on wind energy utilization; it is not only applicable to normally operating wind turbine blades, but can also be used to analyze the wind energy utilization rate of blade elements under different icing forms and thicknesses. This makes the method widely applicable and can cope with wind turbine performance evaluation under different working conditions; by predicting and analyzing the wind energy utilization rate of wind turbine blades under icing conditions in advance, performance degradation and equipment damage caused by icing can be effectively avoided, thereby saving maintenance costs and increasing power generation efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative labor.

[0050] Figure 1 It is an overall flow chart of the method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation according to the first embodiment of the present invention;

[0051] Figure 2 is a cross-sectional view of a fan blade of a simulation example of a method for calculating wind energy utilization rate of an ice-covered fan blade based on cross-sectional segmentation according to a second embodiment of the present invention;

[0052] Figure 3 The lift and drag coefficient C of the two-dimensional airfoil section of the wind turbine blade is a simulation example of the method for calculating the wind energy utilization rate of the ice-covered wind turbine blade based on section segmentation according to the second embodiment of the present invention. L and C D ;

[0053] Figure 4 It is a force analysis diagram of a blade element of a wind turbine blade of a simulation example of a method for calculating wind energy utilization rate of an ice-covered wind turbine blade based on section segmentation according to a second embodiment of the present invention;

[0054] Figure 5 Schematic diagram of wind energy utilization calculation results of wind turbine blade elements under different ice thicknesses, showing a simulation example of the method for calculating wind energy utilization of ice-covered wind turbine blades based on section segmentation according to the second embodiment of the present invention Figure 1 ;

[0055] Figure 6 Schematic diagram of wind energy utilization calculation results of wind turbine blade elements under different ice thicknesses, showing a simulation example of the method for calculating wind energy utilization of ice-covered wind turbine blades based on section segmentation according to the second embodiment of the present invention Figure 2 ;

[0056] Figure 7 Schematic diagram of wind energy utilization calculation results of wind turbine blade elements under different ice thicknesses, showing a simulation example of the method for calculating wind energy utilization of ice-covered wind turbine blades based on section segmentation according to the second embodiment of the present invention Figure 3 . DETAILED DESCRIPTION

[0057] In order to make the above-mentioned purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are described in detail below in conjunction with the drawings of the specification. Obviously, the described embodiments are part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in the art without creative work should fall within the scope of protection of the present invention.

[0058] In the following description, many specific details are set forth to facilitate a full understanding of the present invention, but the present invention may also be implemented in other ways different from those described herein, and those skilled in the art may make similar generalizations without violating the connotation of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.

[0059] Secondly, the term "one embodiment" or "embodiment" as used herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The term "in one embodiment" that appears in different places in this specification does not necessarily refer to the same embodiment, nor does it refer to a separate or selective embodiment that is mutually exclusive with other embodiments.

[0060] Example 1

[0061] Reference Figure 1 , which is the first embodiment of the present invention, provides a method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation, comprising:

[0062] S1: Cut the wind turbine blades to obtain the two-dimensional cross-sectional data of the wind turbine blades at each position.

[0063] In an embodiment of the present application, sectioning a wind turbine blade includes:

[0064] Assume that the radius of the wind turbine blade is R, and select the positions of 0.25R, 0.5R, 0.75R and R for cutting.

[0065] Specifically, the cross-section cutting method is used to cut the wind turbine blade to obtain two-dimensional cross-section data of the two-dimensional wind turbine blade at each position.

[0066] S2: Based on the wind speed and wind turbine operating parameters, the airflow field of each wind turbine blade two-dimensional airfoil section is calculated to obtain the lift and drag coefficient of the two-dimensional airfoil section;

[0067] Specifically, the finite element method is used to calculate the airflow field of the two-dimensional airfoil section of each fan blade.

[0068] S3.: Based on the two-dimensional cross-section data and the lift-drag coefficient of the two-dimensional airfoil section, a calculation model for the wind energy utilization rate of the wind turbine blade element is established, the force analysis of the two-dimensional cross-section blade element is performed, and the lift-drag coefficient of the blade element is calculated;

[0069] In the embodiment of the present application, calculating the blade lift resistance includes:

[0070] Based on the results of the two-dimensional cross-section blade element force analysis, the expression of the blade element lift resistance is obtained as follows:

[0071]

[0072] Where ρ is the air density in kg / m 3 ; c is the blade chord length, in m; r is the radius of the blade, in m; C L and C D are the lift coefficient of the two-dimensional airfoil section and the drag coefficient of the two-dimensional airfoil section; V rel It is the velocity of blade element relative to airflow, in m / s.

[0073] S4: Calculate the thrust and torque on the blade element according to the lift and resistance of the blade element;

[0074] In the embodiment of the present application, calculating the thrust and torque of the blade element includes:

[0075] The turning thrust dF received by the blade element t It is expressed as:

[0076]

[0077] Assume the number of blades is B; the blade angular velocity is Ω, in rad / s; the tip speed ratio of the blade element at r is λ r =Ω·r / V; lift-to-drag ratio e = C L / C D ; Then the torque of the micro-circular ring where the blade element is located is:

[0078]

[0079] S5: Based on the thrust and torque of the blade element, calculate the wind energy utilization rate of the blade element under normal operation;

[0080] In the embodiment of the present application, calculating the blade element wind energy utilization rate under normal operation includes:

[0081] The blade element position is selected to calculate the blade element wind energy utilization rate under normal operation. When calculating, it is assumed that the blade element relative airflow velocity and the components perpendicular and parallel to the blade rotation plane have the following relationship:

[0082]

[0083] Among them, V y and V x are the components perpendicular and parallel to the blade rotation plane respectively.

[0084] Specifically, four typical blade element positions of 0.25R, 0.5R, 0.75R and R are selected to calculate the blade element wind energy utilization rate under normal operation.

[0085] It should be noted that when calculating, V rel and the components V perpendicular and parallel to the plane of rotation of the blade y and V x The relationship is derived under the premise that the fan blades work in an ideal stable state.

[0086] Blade output power dP and wind energy utilization rate C p-r for:

[0087]

[0088] It should be noted that the blade element output power dP and wind energy utilization rate are calculated under the ideal chord length.

[0089] In another possible implementation, a two-dimensional cross-section of an ice-covered blade element is drawn according to the ice-covering shape and thickness of the wind turbine blade, and the wind energy utilization of the blade element under different ice-covering shapes and thicknesses is calculated using the methods provided in S1-S5; the changes in the wind energy utilization of the blade element at typical positions of the wind turbine blades before and after ice covering are compared, and the overall operating efficiency of the wind turbine is inferred.

[0090] The above is a schematic scheme of the wind energy utilization rate calculation method for ice-covered wind turbine blades based on section segmentation in this embodiment. It should be noted that the technical scheme of the wind energy utilization rate calculation system for ice-covered wind turbine blades based on section segmentation and the technical scheme of the wind energy utilization rate calculation method for ice-covered wind turbine blades based on section segmentation belong to the same concept. For details not described in detail in the technical scheme of the wind energy utilization rate calculation system for ice-covered wind turbine blades based on section segmentation in this embodiment, please refer to the description of the technical scheme of the wind energy utilization rate calculation method for ice-covered wind turbine blades based on section segmentation.

[0091] The wind energy utilization rate calculation system of ice-covered wind turbine blades based on section segmentation in this embodiment includes:

[0092] A cross-sectional data acquisition module is used to section the wind turbine blades and obtain two-dimensional cross-sectional data of the two-dimensional wind turbine blades at each position;

[0093] A lift and drag coefficient acquisition module is used to calculate the airflow field of the two-dimensional airfoil section of each fan blade based on the wind speed and the fan operating parameters to obtain the lift and drag coefficient of the two-dimensional airfoil section;

[0094] The lift and drag calculation module is used to establish a wind energy utilization rate calculation model for the blade element of a wind turbine blade based on the two-dimensional cross-section data and the lift and drag coefficient of the two-dimensional airfoil section, perform a force analysis on the two-dimensional cross-section blade element, and calculate the lift and drag of the blade element;

[0095] A parameter calculation module, used to calculate the thrust and torque of the blade element according to the lift and drag of the blade element;

[0096] The utilization rate calculation module is used to calculate the wind energy utilization rate of the blade element under normal operation based on the thrust and torque received by the blade element.

[0097] This embodiment further provides a computing device, which is applicable to the case of a method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation, and includes:

[0098] A memory and a processor; the memory is used to store computer executable instructions, and the processor is used to execute computer executable instructions to implement the method for calculating wind energy utilization of ice-covered wind turbine blades based on profile segmentation as proposed in the above embodiment.

[0099] This embodiment further provides a storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation as proposed in the above embodiment.

[0100] The storage medium proposed in this embodiment and the method for calculating wind energy utilization rate of ice-covered wind turbine blades based on cross-section segmentation proposed in the above embodiment belong to the same inventive concept. Technical details not fully described in this embodiment can be referred to the above embodiment, and this embodiment has the same beneficial effects as the above embodiment.

[0101] Example 2

[0102] Reference Figure 2-Figure 7 , which is an embodiment of the present invention, provides a method for calculating the wind energy utilization rate of ice-covered wind turbine blades based on profile segmentation. In order to verify the beneficial effects of the present invention, scientific demonstration is carried out through simulation experiments.

[0103] Step 1: Use the cross-section cutting method to cut the wind turbine blades to obtain the two-dimensional cross-section data of the wind turbine blades at each position. Figure 2 As shown, assuming that the radius of the wind turbine blade is R, the wind turbine blade positions selected are 0.25R, 0.5R, 0.75R and R for sectioning.

[0104] Step 2: Based on the operating parameters such as wind speed and fan speed, the finite element method is used to calculate the airflow field of the two-dimensional airfoil section of each fan blade to obtain the lift and drag coefficient C of the two-dimensional airfoil section L and C D ;like Figure 3 As shown, the wind energy utilization rate of the blade elements of the fan blades at different positions is different and changes with the change of the rotation speed.

[0105] Step 3: Establish a calculation model for wind energy utilization of fan blades, such as Figure 4 As shown, the force analysis of the two-dimensional cross-section blade element is carried out, and the calculation formula of the blade element lift and resistance is derived;

[0106]

[0107] Where ρ is the air density, kg / m 3 ; c is the blade chord length, m; r is the radius of the blade, m; C L and C D are lift and drag coefficients respectively; V rel is the blade element speed relative to airflow, m / s.

[0108] Step 4: Further derive the thrust, torque and wind energy efficiency calculation formulas on the blade element based on the blade element lift resistance; the turning thrust dF on the blade element t It can be expressed as:

[0109]

[0110] Assume the number of blades is B; the blade angular velocity is Ω, rad / s; the tip speed ratio of the blade element at r is λ r =Ω·r / V; lift-to-drag ratio e = C L / C D ; Then the torque of the micro-circular ring where the blade element is located is:

[0111]

[0112] The blade output power dP and wind energy utilization rate C under the ideal chord length p-r for:

[0113]

[0114] Step 5: Select four typical blade positions of 0.25R, 0.5R, 0.75R and R to calculate the blade wind energy utilization rate under normal operation; when calculating, consider that the fan blade works in an ideal stable state, V rel and the components V perpendicular and parallel to the plane of rotation of the blade y and V x With the following relationship:

[0115]

[0116] Step 6: Draw the two-dimensional cross section of the blade element covered with ice according to the ice morphology and thickness of the wind turbine blade, and use the steps from step 2 to step 5 to calculate the wind energy utilization rate of the blade element under different ice morphology and thickness. Figure 5 As shown in the figure, the blade cross-section under different ice thickness and shape is drawn, and the lift and drag are calculated, as shown in the figure. Figure 6 As shown; Finally, the wind energy utilization rate C of blades with different thickness and shape can be obtained p-r Change rules, such as Figure 7 shown.

[0117] Step 7: Compare the changes in wind energy utilization of blade elements at typical positions of wind turbine blades before and after icing, and infer the overall operating efficiency of the wind turbine.

[0118] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

Claims

1. A method for calculating wind energy utilization rate of iced wind turbine blades based on section segmentation, characterized in that: include: Cut the wind turbine blades to obtain the two-dimensional cross-sectional data of the wind turbine blades at each position; Based on the wind speed and wind turbine operating parameters, the airflow field of each wind turbine blade two-dimensional airfoil section is calculated to obtain the lift and drag coefficient of the two-dimensional airfoil section; Based on the two-dimensional cross-sectional data and the lift-drag coefficient of the two-dimensional airfoil section, a calculation model for the wind energy utilization rate of the fan blade element is established, the force analysis of the two-dimensional cross-sectional blade element is carried out, and the lift-drag coefficient of the blade element is calculated; According to the lift and resistance of the blade element, calculate the thrust and torque of the blade element; Based on the thrust and torque acting on the blade element, the wind energy utilization rate of the blade element under normal operation is calculated.

2. The method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation according to claim 1, characterized in that: The sectioning of the wind turbine blade comprises: Assume that the radius of the wind turbine blade is R, and select the positions of 0.25R, 0.5R, 0.75R and R for cutting.

3. The method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation according to claim 2, characterized in that: The calculation of blade lift resistance includes: Based on the results of the two-dimensional cross-section blade element force analysis, the expression of the blade element lift resistance is obtained as follows: Where ρ is the air density in kg / m 3 ; c is the blade chord length, in m; r is the radius of the blade, in m; C L and C D are the lift coefficient of the two-dimensional airfoil section and the drag coefficient of the two-dimensional airfoil section; V rel It is the velocity of blade element relative to airflow, in m / s.

4. The method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation according to claim 3, characterized in that: The calculation of the thrust and torque of the blade element includes: The turning thrust dF received by the blade element t It is expressed as:

5. The method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation according to claim 4, characterized in that: The calculation of the thrust and torque of the blade element also includes: Assume the number of blades is B; the blade angular velocity is Ω, in rad / s; the tip speed ratio of the blade element at r is λ r =Ω·r / V; lift-to-drag ratio e = C L / C D ; Then the torque of the micro-circular ring where the blade element is located is:

6. The method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation according to claim 5, characterized in that: The calculation of blade element wind energy utilization under normal operation includes: The blade element position is selected to calculate the blade element wind energy utilization rate under normal operation. When calculating, it is assumed that the blade element relative airflow velocity and the components perpendicular and parallel to the blade rotation plane have the following relationship: Among them, V y and V x are the components perpendicular and parallel to the blade rotation plane respectively.

7. The method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation according to claim 6, characterized in that: The calculation of blade element wind energy utilization under normal operation also includes: Blade output power dP and wind energy utilization rate C p-r for:

8. A system using the method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation as claimed in any one of claims 1 to 7, characterized in that: include: A cross-sectional data acquisition module is used to section the wind turbine blades and obtain two-dimensional cross-sectional data of the two-dimensional wind turbine blades at each position; A lift and drag coefficient acquisition module is used to calculate the airflow field of the two-dimensional airfoil section of each fan blade based on the wind speed and the fan operating parameters to obtain the lift and drag coefficient of the two-dimensional airfoil section; The lift and drag calculation module is used to establish a wind energy utilization rate calculation model for the blade element of a wind turbine blade based on the two-dimensional cross-section data and the lift and drag coefficient of the two-dimensional airfoil section, perform force analysis on the two-dimensional cross-section blade element, and calculate the lift and drag of the blade element; A parameter calculation module is used to calculate the thrust and torque of the blade element according to the lift and drag of the blade element; The utilization rate calculation module is used to calculate the wind energy utilization rate of the blade element under normal operation based on the thrust and torque received by the blade element.

9. A computing device comprising: Memory and processor; The memory is used to store computer executable instructions, and the processor is used to execute the computer executable instructions. When the computer executable instructions are executed by the processor, the steps of the method for calculating wind energy utilization of ice-covered wind turbine blades based on profile segmentation as described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the steps of the method for calculating wind energy utilization rate of ice-covered wind turbine blades based on section segmentation as described in any one of claims 1 to 7.

Images