Wind turbine generator life estimation method and device, storage medium and electronic equipment

Abstract

The invention relates to a wind turbine generator life estimation method and device, a storage medium and electronic equipment, and relates to the technical field of wind power. And life estimation can be carried out on different types of wind turbine generators. The method comprises the steps that parameter information of the wind turbine generator is obtained, a finite element model corresponding to the wind turbine generator is built according to the parameter information, and the parameter information comprises the size and material of all parts in the wind turbine generator and environment data of the installation position of the wind turbine generator; multi-working-condition simulation is conducted through the finite element model, simulation damage of all key positions in the wind turbine generator set under different working conditions is determined, and the key positions represent the positions with weak strength in the wind turbine generator set; the working condition frequencies of the wind turbine generator set under different machine positions are determined, and according to the working condition frequencies and the simulated damage of the key positions under different working conditions, the estimated service life of each component of the wind turbine generator set under different machine positions is determined; and determining the minimum value in the estimated life of each component in the wind turbine generator set as the estimated life of the wind turbine generator set.

Description

Technical Field

[0001] The present disclosure relates to the field of wind power technology, and in particular to a method, device, storage medium and electronic equipment for estimating the life of a wind turbine generator set. Background Art

[0002] At present, in terms of wind turbine life estimation, most of the wind turbine life is determined based on the research results and field experience of each manufacturer. However, the research results and field experience of each manufacturer are closely related to the R&D personnel and production personnel, and are affected by the subjectivity of the R&D personnel or production personnel, and cannot be promoted. There is a lack of wind turbine life estimation methods that are uniformly accepted by the industry. Summary of the invention

[0003] In order to solve the deficiencies of the prior art, the present disclosure provides a method, device, storage medium and electronic device for estimating the life of a wind turbine generator set.

[0004] In order to achieve the above objectives, in a first aspect, the present disclosure provides a method for estimating the life of a wind turbine generator set, comprising:

[0005] Acquire parameter information of a wind turbine generator set, and construct a finite element model corresponding to the wind turbine generator set according to the parameter information, wherein the parameter information includes the size and material of each component in the wind turbine generator set and the environmental data of the installation location of the wind turbine generator set;

[0006] Performing multi-operating condition simulations using the finite element model to determine simulated damage at key locations in the wind turbine generator set under different operating conditions, wherein the key locations represent locations with weak strength in the wind turbine generator set;

[0007] Determine the operating condition frequency of the wind turbine generator set at different machine positions, and determine the estimated life of each component of the wind turbine generator set at different machine positions according to the operating condition frequency and the simulated damage of each key position under different operating conditions;

[0008] The minimum value of the estimated lifespans of the components in the wind turbine generator set is determined as the estimated lifespan of the wind turbine generator set.

[0009] In a second aspect, the present disclosure provides a device for estimating the life of a wind turbine generator set, comprising:

[0010] A modeling module is configured to obtain parameter information of a wind turbine generator set and construct a finite element model corresponding to the wind turbine generator set according to the parameter information, wherein the parameter information includes the size and material of each component in the wind turbine generator set and the environmental data of the installation location of the wind turbine generator set;

[0011] A simulation module is configured to perform multi-operating condition simulation through the finite element model to determine simulated damage of each key position in the wind turbine generator set under different operating conditions, wherein the key position represents a position with weak strength in the wind turbine generator set;

[0012] A first execution module is configured to determine the operating condition frequency of the wind turbine generator set at different machine positions, and determine the estimated life of each component of the wind turbine generator set at different machine positions according to each operating condition frequency and the simulated damage of each key position under different operating conditions;

[0013] The second execution module is configured to determine the minimum value of the estimated lifespans of the components in the wind turbine generator set as the estimated lifespan of the wind turbine generator set.

[0014] In a third aspect, the present disclosure provides a non-temporary computer-readable storage medium having a computer program stored thereon, which implements the method described in the first aspect when executed by a processor.

[0015] In a fourth aspect, the present disclosure provides an electronic device, including:

[0016] a memory having a computer program stored thereon;

[0017] A processor is used to execute the computer program in the memory to implement the method described in the first aspect.

[0018] In a fifth aspect, the present disclosure provides a computer program product, including a computer program, which implements the method described in the first aspect when executed by a processor.

[0019] Through the above technical solution, a finite element model is constructed according to the parameter information of the wind turbine, and the working condition simulation is performed to determine the simulated damage of each key position in the wind turbine under different working conditions. According to each simulated damage and the working condition frequency of the wind turbine at different machine positions, the estimated life of each component of the wind turbine at different machine positions is determined, and then the life of the wind turbine is obtained. The whole process is not affected by the subjectivity of personnel and can be used to estimate the life of different types of wind turbines.

[0020] Other features and advantages of the present disclosure will be described in detail in the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] The accompanying drawings are used to provide a further understanding of the present disclosure and constitute a part of the specification. Together with the following specific embodiments, they are used to explain the present disclosure but do not constitute a limitation of the present disclosure. In the accompanying drawings:

[0022] Figure 1 The figure is a flow chart of a method for estimating the life of a wind turbine generator system according to an exemplary embodiment of the present disclosure.

[0023] Figure 2a to Figure 2b is a cross-sectional view of a flange in a wind turbine generator set according to an exemplary embodiment of the present disclosure.

[0024] Figure 3 SN curve of a flange weld according to an exemplary embodiment of the present disclosure.

[0025] Figure 4 It is a schematic diagram showing the relationship between various dangerous points on a flange weld and the resultant bending moment according to an exemplary embodiment of the present disclosure.

[0026] Figure 5 is a curve showing the relationship between the load and the maximum principal stress at each dangerous point in the flange weld according to an exemplary embodiment of the present disclosure.

[0027] Figure 6 It is a schematic diagram of the load-stress conversion relationship corresponding to the flange weld according to an exemplary embodiment of the present disclosure.

[0028] Figure 7 The flowchart of determining the simulated damage of a flange weld under fatigue conditions is shown according to an exemplary embodiment of the present disclosure.

[0029] Figure 8 The figure is a flow chart of a method for estimating the life of a wind turbine generator system according to an exemplary embodiment of the present disclosure.

[0030] Fig. 9 The present invention is a block diagram of a device for estimating the life of a wind turbine generator system according to an exemplary embodiment of the present invention.

[0031] Fig.10 is a block diagram of an electronic device according to an exemplary embodiment of the present disclosure. DETAILED DESCRIPTION

[0032] The specific implementation of the present disclosure is described in detail below in conjunction with the accompanying drawings. It should be understood that the specific implementation described herein is only used to illustrate and explain the present disclosure, and is not used to limit the present disclosure.

[0033] It should be noted that all actions of acquiring signals, information or data in the present disclosure are carried out in compliance with the relevant data protection laws and policies of the country where the device is located and with the authorization given by the owner of the corresponding device.

[0034] As mentioned in the background technology, at present, the life estimation of wind turbines is subject to the subjective influence of R&D personnel or production personnel and cannot be promoted.

[0035] In view of this, the present disclosure provides a method, device, storage medium and electronic device for estimating the life of a wind turbine generator set, which can estimate the life of different types of wind turbine generator sets without being affected by human subjectivity.

[0036] Figure 1 The following is a flow chart of a method for estimating the life of a wind turbine generator set according to an exemplary embodiment of the present disclosure. The method for estimating the life of a wind turbine generator set can be used in devices such as personal computers, laptops, and smart phones. The method for estimating the life of a wind turbine generator set may include the following steps:

[0037] In step S11, parameter information of the wind turbine is obtained, and a finite element model corresponding to the wind turbine is constructed according to the parameter information, wherein the parameter information includes the size and material of each component in the wind turbine and the environmental data of the installation location of the wind turbine.

[0038] It is worth noting that the components in a wind turbine include but are not limited to hubs, main shafts, frames, bearings, towers, and flanges. Each component may have multiple forms. For example, the flanges in a wind turbine include: Figure 2a The L-shaped flange shown in Figure 2b T-flange as shown. Figure 2a As shown in the figure, two L-shaped flanges are connected by bolts, where the first flange has a heavy wall of hw1, a bolt hole center circle diameter of bf1, and a thickness of t1, and the second flange has a heavy wall of hw2, a bolt hole center circle diameter of bf2, and a thickness of t2. Figure 2b As shown in the figure, the high neck of the T-type flange is hw, the height of the flange face is hf, and the thickness is t. In general, the flange is forged or rolled as a whole, and the outer end of the flange extends a section to be welded to the tower wall steel plate of the tower. Therefore, the key position on the flange is at the weld toe, that is, Figure 2a and Figure 2b The black triangle position in .

[0039] It is worth noting that the finite element model can be constructed according to internationally recognized wind power industry standards, such as the wind turbine certification guidelines issued by GL 2010 (Germanischer Lloyd).

[0040] In step S12, multi-operating condition simulation is performed using the finite element model to determine simulated damage of key positions in the wind turbine generator set under different operating conditions, wherein the key positions represent positions with weak strength in the wind turbine generator set.

[0041] It is worth noting that the key positions of the wind turbine with weak strength can be determined by performing load analysis on the structural characteristics of each component in the wind turbine. The load types involved in the load analysis process include but are not limited to wind load, gravity load and tower shadow effect, and the specific methods involved in the load analysis include but are not limited to static strength analysis and fatigue strength analysis.

[0042] In step S13, the operating condition frequencies of the wind turbine generator set at different machine positions are determined, and the estimated life of each component of the wind turbine generator set at different machine positions is determined according to the operating condition frequencies and the simulated damage of each key position under different operating conditions.

[0043] In step S14, the minimum value of the estimated lifespans of the components in the wind turbine generator set is determined as the estimated lifespan of the wind turbine generator set.

[0044] It is worth mentioning that the components in the wind turbine that need to be simulated for damage calculation, the corresponding simulation conditions of the finite element model, the input loads during the simulation process, and the number of dangerous points at key positions on each component can be found in Table 1 below.

[0045] Table 1

[0046]

[0047]

[0048] In the disclosed embodiment, a corresponding finite element model is constructed according to the parameter information of the wind turbine, and then the finite element model is used for simulation to determine the simulated damage of each key position in the wind turbine under different working conditions, and the estimated life of each component of the wind turbine under different machine positions is determined in combination with the working condition frequency of the wind turbine at different positions, and the estimated life of the wind turbine is obtained according to the estimated life of each component. The whole process does not require human intervention, and the life of the wind turbine can be estimated in a targeted manner, which is not limited by the experience of personnel and the type of wind turbine, and can realize the life estimation of different types of wind turbines, and the estimation result is reliable.

[0049] In order to facilitate those skilled in the art to better understand the wind turbine life estimation method provided by the present disclosure, the steps involved in the method are described below.

[0050] In a feasible implementation manner, in step S12, performing multi-operating condition simulations through the finite element model to obtain simulated damage of each key position in the wind turbine under different operating conditions may include:

[0051] Determine a plurality of fatigue working conditions corresponding to the wind turbine generator set, and simulate each of the fatigue working conditions for a preset time using the finite element model to obtain a simulated stress spectrum of each of the key positions in the wind turbine generator set, wherein the simulated stress spectrum represents the simulated stress amplitude of the key position under different fatigue working conditions and the number of simulation cycles corresponding to the simulated stress amplitude;

[0052] For each of the fatigue conditions, rain flow counting is performed on the simulated stress spectrum of the key position under the fatigue condition to obtain the simulated stress characteristics of the key position within the preset time period, wherein the simulated stress characteristics include a simulated stress range and a number of simulated stress changes;

[0053] According to the simulated stress characteristics and the material stress-life curve of the key position, linear cumulative damage processing is performed to obtain the simulated damage of the key position under the fatigue condition.

[0054] The preset time can be preset according to the accuracy and cost requirements of the life estimation result. The longer the preset time is, the higher the accuracy of the life estimation result is and the higher the cost is. In the embodiment of the present disclosure, the preset time is 10 minutes.

[0055] It is worth noting that the load calculation can be performed by combining the wind resource data of the wind turbine at different positions, such as wind speed, wind direction, turbulence and wind shear, to obtain multiple fatigue conditions at key positions on the wind turbine. In the embodiment of the present disclosure, 1053 fatigue conditions are taken as an example.

[0056] It is worth noting that the rain flow counting processing of the stress spectrum can include: first, filtering out the small amplitude loads that do not cause damage in the stress spectrum, removing the intermediate transition points, and retaining only the peak-to-valley load values ​​to achieve peak-to-valley filtering; then, discretizing each simulated stress amplitude in the stress spectrum after peak-to-valley filtering, and dividing each simulated stress amplitude into multiple equivalent regions according to the peak-to-peak value, so as to simplify the stress spectrum into several load cycles, thereby determining the simulated stress range corresponding to several load cycles and the number of simulated stress transformations.

[0057] In a feasible implementation manner, performing linear cumulative damage processing according to the simulated stress characteristics and the material stress-life curve of the key position to obtain the simulated damage of the key position under the fatigue condition may include:

[0058] According to the material safety factor corresponding to the key position, the simulation stress range is corrected to obtain a simulation target stress range;

[0059] Determining a simulation target stress at the key position according to the simulation stress range and the stress-life curve at the key position;

[0060] According to the simulation target stress, the simulation target stress range and the number of simulation stress changes, the simulated damage of the key position under the fatigue condition is obtained.

[0061] It is worth noting that the material safety factor and the simulation stress range are substituted into the following calculation formula to obtain the corrected simulation target stress range:

[0062] `

[0063] Δσ i =γ m Δσ i ,

[0064] Among them, Δσ` i Characterize the simulation stress range, γ m Characterizes the material safety factor, generally taken as 1.15, Δσ i Characterize the target stress range for simulation.

[0065] The vertical axis of the stress life curve (hereinafter referred to as SN curve) represents the stress range (Stress rangeΔσ), and the horizontal axis represents the number of stress cycles (Number of stress cycles N). The DC value of the SN curve (representing the fatigue limit or endurance limit of the material) generally refers to the EN 1993-1-9:2005 standard. Since the number of cycles of Eurocode 3 in the EN 1993-1-9:2005 standard is greater than 1e8 times, it is determined to be infinite fatigue, and the number of cycles within the life cycle of the wind turbine will exceed 1e8 times. Therefore, the SN curve needs to be corrected according to GL2010.

[0066] For example, taking the SN curve corresponding to the flange as an example, for any fatigue condition, such as Figure 3 As shown in the figure, the DC value refers to the SN curve of Eurocode 3 as a solid line (i.e., S / N curve according to Eurocode 3), and the SN curve after the DC value of the SN curve is corrected according to GL2010 is a dotted line (i.e., S / N curve according to IIW). A The corresponding stress Δσ A =DC value, N A Characterizes the maximum stress level at which a material does not suffer fatigue failure under a specific number of cycles, also known as the conditional fatigue limit; turning point N B The corresponding stress m = 3. DCharacterizes the fatigue strength of a material under a specific number of cycles, also known as the edge load cycle number. When the DC value of the SN curve corresponding to Eurocode 3 is 90, N A =2.0e6, N D =5.0e6. According to the comparison results of the simulation stress range and the simulation target stress range, it is determined that the simulation target stress has the following two situations:

[0067] Case 1: When Δσ` i >Δσ D hour, m=3.

[0068] Case 2: When Δσ` i <Δσ D hour, m=5.

[0069] Set the simulation target stress N i , simulation target stress range n and simulation stress change times i, are substituted into the following calculation formula to obtain the simulated damage D of each dangerous point on the flange weld under this fatigue condition. j :

[0070]

[0071] In a feasible implementation manner, simulating each fatigue condition for a preset time period using the finite element model to obtain a simulated stress spectrum of each key position in the wind turbine generator set may include:

[0072] For each of the fatigue working conditions, the finite element model is used to simulate the fatigue working condition for a preset time period to obtain a simulated load time history of each of the key positions within the preset time period, wherein the simulated load time history represents a process in which the fatigue load at the key position changes with time within the preset time period;

[0073] Determine the simulated stress time history corresponding to the simulated load time history of each key position according to the load-stress conversion relationship corresponding to each key position;

[0074] According to the fatigue influence factors corresponding to the key positions, the stress at each moment in the stress time history corresponding to the key positions is corrected to obtain the simulated stress spectrum of the key positions.

[0075] It is worth noting that the fatigue influence factors corresponding to different key positions are different, which can be determined by looking up the table "Table 8.3 Transverse butt welds" on page 23 of the standard Eurocode3. The stress at each moment in the stress time history is corrected according to the fatigue influence factors corresponding to each key position. For example, "Table 8.3 Transverse butt welds" records that when the thickness of the L-type flange weld is greater than 25cm, its fatigue influence factor is the thickness influence factor, which can be obtained by the following calculation formula:

[0076] n=0.2.

[0077] Among them, f(t) represents the thickness influence factor, and t1 represents the thickness of the L-type flange.

[0078] For example, taking the L-type flange weld as an example, the finite element model is used to simulate each fatigue condition for a preset time to obtain the simulated load time history of the flange weld within the preset time; the flange weld is discretized into several dangerous points, and the simulated stress time history corresponding to the simulated load time history at each dangerous point is determined according to the load-stress conversion relationship corresponding to the flange weld; according to the thickness influence factor corresponding to the flange weld, the stress at each moment in the stress time history corresponding to each dangerous point is corrected to obtain the stress spectrum of the flange weld.

[0079] In a feasible implementation manner, determining the simulated stress time history corresponding to the simulated load time history of each key position according to the load-stress conversion relationship corresponding to each key position may include:

[0080] For each time point in the simulated load time history of each key position, determining the simulated combined bending moment and the simulated limit bending moment of the key position under the fatigue condition;

[0081] Determine the angle between the simulated combined bending moment and the key position in the horizontal plane coordinate system, and determine, based on the load-stress conversion relationship, the second simulated principal stress and the third simulated principal stress corresponding to the azimuth angle closest to the angle, wherein the second simulated principal stress represents the principal stress of the key position under the fatigue condition, and the third simulated principal stress represents the principal stress of the key position under all fatigue conditions;

[0082] According to the simulated combined bending moment, the simulated ultimate bending moment, the second simulated principal stress and the third simulated principal stress, a first simulated principal stress at the key position at the time point is obtained;

[0083] The simulated stress time history of the key position is determined according to the simulated principal stress at each time point in the simulated load time history of the key position.

[0084] For example, the L-type flange weld is divided into several dangerous points, such as Figure 4 As shown in the figure, taking the dangerous point corresponding to the azimuth angle α as an example, the stress spectrum of this point under the bolt preload displacement condition can be determined by the following method:

[0085] I) Determine the simulated external loads Mx and My at the dangerous point corresponding to the azimuth angle α under the bolt preload displacement condition, and calculate the simulated combined bending moment Mxy and the simulated ultimate bending moment M based on the simulated external loads Mx and My.

[0086] II) Determine the angle β between the simulated combined bending moment Mxy and the positive direction of the x-axis in the horizontal coordinate system XOY, and obtain the angle θ=α-β between the simulated combined bending moment Mxy and the dangerous point in the horizontal coordinate system XOY. And according to the load-stress conversion relationship, determine the second simulated principal stress σ corresponding to the azimuth closest to the angle θ under the action of the bolt preload displacement condition v , and the third simulated principal stress σ1 under the combined working condition of bolt preload displacement and simulated ultimate bending moment M, we can get Figure 5 The relationship curve between the load and the maximum principal stress at the dangerous point is shown.

[0087] III) The simulated bending moment Mxy, the simulated ultimate bending moment M, and the second simulated principal stress σ v And the second simulation principal stress σ1, substitute into the following calculation formula to obtain the first simulation principal stress σ of the danger point at this time point n :

[0088]

[0089] By solving the load at each time point in the simulated load time history in the above manner, the simulated stress time history of the dangerous point of the weld can be obtained.

[0090] In a feasible implementation manner, the load-stress conversion relationship corresponding to the key position can be determined in the following manner:

[0091] Performing static analysis on the structure of the wind turbine generator set by using the finite element model to obtain the first simulated principal stress of each key position in the wind turbine generator set under different working conditions;

[0092] The load-stress conversion relationship of the key position is determined according to the azimuth angle of each key position and the first simulated principal stress, wherein the azimuth angle represents the angle between the key position along the clockwise direction and the positive direction of the x-axis in the horizontal plane coordinate system.

[0093] It is worth noting that the simulated principal stresses of each key position under different working conditions can also be determined according to the requirements of the GL 2010 standard, and the present disclosure does not limit this.

[0094] For example, taking the L-type flange weld as an example, the L-type flange is subjected to static analysis through the finite element model. For example, the lower outer ring of the L-type flange weld position is taken, and the maximum principal stress of all nodes of the lower outer ring under the condition of only bolt preload displacement and the maximum principal stress under the condition of combined bolt preload displacement and simulated limit bending moment are calculated through the finite element model. Figure 6 The distribution relationship between the maximum principal stress and the node azimuth under the two working conditions shown is the load-stress conversion relationship of the flange weld. Figure 6 Where Time=1.0 represents the load-stress conversion relationship under the working condition of bolt preload displacement only, and Time=2.0 represents the load-stress conversion relationship under the working condition of combined bolt preload displacement and simulated ultimate bending moment M.

[0095] The following is an example of flange weld to fully explain the simulation damage of flange weld under fatigue conditions. Figure 7 As shown:

[0096] In step S71, the simulated load time history at the flange weld in the target wind turbine generator set and the load-stress conversion relationship corresponding to the flange weld are determined.

[0097] In step S72, the stress spectrum at the flange weld is obtained according to the inverse load-time history and the load-stress conversion relationship, and step S75 is executed.

[0098] In step S73, the SN curve of the flange weld material is determined, and step S75 is executed.

[0099] In step S74, the working condition frequency of the flange weld at different machine positions is determined, and step S75 is executed.

[0100] In step S75, the simulated damage of the flange weld under different working conditions is obtained according to the SN curve, stress spectrum and working condition frequency.

[0101] In step S76, a linear cumulative damage process is performed on the simulated damage under different working conditions to obtain the estimated life of the flange weld under different machine positions.

[0102] In a feasible implementation manner, in step S13, determining the operating condition frequency of the wind turbine generator set at different machine positions may include:

[0103] Determining wind resource data of the wind turbine generator set at different positions, wherein the wind resource data includes at least one of wind speed, wind direction, turbulence and wind shear;

[0104] According to the wind resource data, load simulation is performed on the wind turbine generator set to obtain the operating condition frequency of the wind turbine generator set at different machine positions.

[0105] It is worth noting that, based on wind resource data, the load conditions of wind turbines under different wind conditions are simulated to determine the equivalent fatigue load, and then the working condition frequencies under different machine positions are obtained. For example, by analyzing wind resource data, the inflow angle, turbulence and wind frequency of different sectors of wind turbines under different wind speeds are obtained, and working conditions are assigned to each wind speed segment of each sector. The frequency of each working condition is counted to obtain the working condition frequency of wind turbines under different machine positions.

[0106] For example, when the load simulation duration is 600 seconds, the operating condition frequency of the i-th sector in the wind turbine in the m-th wind speed range can be calculated by the following formula:

[0107]

[0108] Among them, n j Represents the frequency of the jth operating condition within 600 seconds, P im Represents the wind frequency corresponding to the i-th sector in the m-th wind speed range.

[0109] In a feasible implementation manner, in step S13, determining the estimated life of each component of the wind turbine generator set at different machine positions according to the frequency of each operating condition and the simulated damage of each key position under different operating conditions may include:

[0110] For each key position, according to the simulated damage of the key position and the frequency of each working condition, the intermediate simulated damage of each position of the key position is obtained;

[0111] Adding the intermediate simulated damage of the key position at each position to obtain the total simulated damage of the component at the key position;

[0112] The estimated life of the components is determined based on the total simulated damage of each of the components.

[0113] For example, the simulated damage of the key position in the jth working condition and the working condition frequency are substituted into the following calculation formula to obtain the total simulated damage D (i.e., intermediate simulated damage) of the key position component in each position:

[0114]

[0115] Based on the total simulated damage of each component, the estimated life of the component is determined by searching the SN curve of the component.

[0116] The following is a complete example to illustrate the wind turbine life estimation method provided by the present disclosure. Figure 8 As shown:

[0117] In step S81, the size and material of each component in the target wind turbine generator set and the environmental data of the installation location of the target wind turbine generator set are obtained.

[0118] In step S82, a finite element model of the corresponding target wind turbine is constructed according to the environmental data and the size and material of each component.

[0119] In step S83, multi-operating condition simulation is performed through a finite element model to obtain the load-stress conversion relationship corresponding to each key position in the target wind turbine and the simulated load time history of each key position within a preset time length.

[0120] In step S84, the simulated stress time history corresponding to each key position in the target wind turbine is obtained according to the load-stress conversion relationship and the simulated load time history, and step S86 is executed.

[0121] In step S85, the SN curve of each component is determined, and step S86 is executed.

[0122] In step S86, the simulated damage of each key position in the target wind turbine under different working conditions is obtained according to the simulated stress time history of each key position and the SN curve of each component, and step S89 is executed.

[0123] In step S87, the wind resource data of the target wind turbine generator set at different machine positions are determined.

[0124] In step S88, the operating frequency of the wind turbine generator set at different machine positions is obtained according to the wind resource data, and step S89 is executed.

[0125] In step S89, the estimated life of each component of the target wind turbine generator set at different machine positions is obtained according to the simulated damage of each key position under different working conditions and the working condition frequency of the wind turbine generator set at different machine positions.

[0126] Taking 66 V42 wind turbines as an example, the estimated life of each component and the estimated life of the wind turbine obtained by the above method can be seen in Table 2 below.

[0127] Table 2

[0128]

[0129]

[0130]

[0131]

[0132] In the disclosed embodiment, a corresponding finite element model is constructed according to the life parameters of the wind turbine generator set, and multi-operating condition simulation is performed through the finite element model corresponding to the wind turbine generator set to determine the simulated damage of each key position in the wind turbine generator set under different operating conditions. There is no need to calculate the simulated damage of all positions in the wind turbine generator set, which reduces the amount of calculation and thus shortens the simulation time. According to the simulated damage of each key position under different operating conditions and the operating condition frequency of the wind turbine generator set at different machine positions, the estimated life of each component of the wind turbine generator set at different machine positions is determined, thereby determining the estimated life of the wind turbine generator set, thereby improving the accuracy and reliability of the life estimation results, and is not limited by the type of wind turbine generator set, and life estimation can be performed for wind turbine generator sets of different models.

[0133] Based on the same inventive concept, the present disclosure provides a device for predicting the life of a wind turbine generator set, such as Fig. 9 As shown, the wind turbine life estimation device includes a modeling module 901 , a simulation module 902 , a first execution module 903 and a second execution module 904 .

[0134] The modeling module 901 is configured to obtain parameter information of the wind turbine generator set, and construct a finite element model corresponding to the wind turbine generator set according to the parameter information, wherein the parameter information includes the size and material of each component in the wind turbine generator set, and the environmental data of the installation location of the wind turbine generator set;

[0135] A simulation module 902 is configured to perform multi-operating condition simulations through the finite element model to determine simulated damage of each key position in the wind turbine generator set under different operating conditions, wherein the key position represents a position with weak strength in the wind turbine generator set;

[0136] The first execution module 903 is configured to determine the operating condition frequency of the wind turbine generator set at different machine positions, and determine the estimated life of each component of the wind turbine generator set at different machine positions according to each operating condition frequency and the simulated damage of each key position under different operating conditions;

[0137] The second execution module 904 is configured to determine the minimum value of the estimated lifespans of the components in the wind turbine generator set as the estimated lifespan of the wind turbine generator set.

[0138] In the disclosed embodiment, a corresponding finite element model is constructed according to the parameter information of the wind turbine, and then the finite element model is used for simulation to determine the simulated damage of each key position in the wind turbine under different working conditions, and the estimated life of each component of the wind turbine under different machine positions is determined in combination with the working condition frequency of the wind turbine at different positions, and the estimated life of the wind turbine is obtained according to the estimated life of each component. The whole process does not require human intervention, and the life of the wind turbine can be estimated in a targeted manner, which is not limited by the experience of personnel and the type of wind turbine, and can realize the life estimation of different types of wind turbines, and the estimation result is reliable.

[0139] Optionally, the simulation module 902 is configured to determine a plurality of fatigue conditions corresponding to the wind turbine generator set, and simulate each of the fatigue conditions for a preset time using the finite element model to obtain a simulated stress spectrum of each of the key positions in the wind turbine generator set, wherein the simulated stress spectrum represents the simulated stress amplitude of the key position under different fatigue conditions and the number of simulation cycles corresponding to the simulated stress amplitude;

[0140] For each of the fatigue conditions, rain flow counting is performed on the simulated stress spectrum of the key position under the fatigue condition to obtain the simulated stress characteristics of the key position within the preset time period, wherein the simulated stress characteristics include a simulated stress range and a number of simulated stress changes;

[0141] According to the simulated stress characteristics and the material stress-life curve of the key position, linear cumulative damage processing is performed to obtain the simulated damage of the key position under the fatigue condition.

[0142] Optionally, the simulation module 902 is configured to correct the simulation stress range according to the material safety factor corresponding to the key position to obtain a simulation target stress range;

[0143] Determining a simulation target stress at the key position according to the simulation stress range and the stress-life curve at the key position;

[0144] According to the simulation target stress, the simulation target stress range and the number of simulation stress changes, the simulated damage of the key position under the fatigue condition is obtained.

[0145] Optionally, the simulation module 902 is configured to simulate the fatigue condition for a preset time period through the finite element model for each fatigue condition, and obtain a simulated load time history of each key position within the preset time period, wherein the simulated load time history represents a process in which the fatigue load at the key position changes with time within the preset time period;

[0146] Determine the simulated stress time history corresponding to the simulated load time history of each key position according to the load-stress conversion relationship corresponding to each key position;

[0147] According to the fatigue influence factors corresponding to the key positions, the stress at each moment in the stress time history corresponding to the key positions is corrected to obtain the simulated stress spectrum of the key positions.

[0148] Optionally, the simulation module 902 is configured to perform a static analysis on the structure of the wind turbine generator set through the finite element model to obtain the first simulated principal stress of each of the key positions in the wind turbine generator set under different working conditions;

[0149] The load-stress conversion relationship of the key position is determined according to the azimuth angle of each key position and the first simulated principal stress, wherein the azimuth angle represents the angle between the key position along the clockwise direction and the positive direction of the x-axis in the horizontal plane coordinate system.

[0150] Optionally, the simulation module 902 is configured to determine the simulated combined bending moment and the simulated ultimate bending moment of each key position under the fatigue condition at each time point in the simulated load time history of each key position;

[0151] Determine the angle between the simulated combined bending moment and the key position in the horizontal plane coordinate system, and determine, based on the load-stress conversion relationship, the second simulated principal stress and the third simulated principal stress corresponding to the azimuth angle closest to the angle, wherein the second simulated principal stress represents the principal stress of the key position under the fatigue condition, and the third simulated principal stress represents the principal stress of the key position under all fatigue conditions;

[0152] According to the simulated combined bending moment, the simulated ultimate bending moment, the second simulated principal stress and the third simulated principal stress, a first simulated principal stress at the key position at the time point is obtained;

[0153] The simulated stress time history of the key position is determined according to the simulated principal stress at each time point in the simulated load time history of the key position.

[0154] Optionally, the first execution module 903 is configured to obtain, for each key position, an intermediate simulated damage at each position of the key position according to the simulated damage of the key position and the frequency of each working condition;

[0155] Adding the intermediate simulated damage of the key position at each position to obtain the total simulated damage of the component at the key position;

[0156] The estimated life of the components is determined based on the total simulated damage of each of the components.

[0157] Optionally, the first execution module 903 is configured to determine wind resource data of the wind turbine generator set at different positions, wherein the wind resource data includes at least one of wind speed, wind direction, turbulence and wind shear;

[0158] According to the wind resource data, load simulation is performed on the wind turbine generator set to obtain the operating condition frequency of the wind turbine generator set at different machine positions.

[0159] Regarding the wind turbine life prediction device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0160] Based on the same inventive concept, the present disclosure provides an electronic device, including:

[0161] a memory having a computer program stored thereon;

[0162] A processor is used to execute the computer program in the memory to implement the above-mentioned wind turbine life prediction method.

[0163] In the disclosed embodiment, a corresponding finite element model is constructed according to the parameter information of the wind turbine, and then the finite element model is used for simulation to determine the simulated damage of each key position in the wind turbine under different working conditions, and the estimated life of each component of the wind turbine under different machine positions is determined in combination with the working condition frequency of the wind turbine at different positions, and the estimated life of the wind turbine is obtained according to the estimated life of each component. The whole process does not require human intervention, and the life of the wind turbine can be estimated in a targeted manner, which is not limited by the experience of personnel and the type of wind turbine, and can realize the life estimation of different types of wind turbines, and the estimation result is reliable.

[0164] Fig.10 FIG. 1 is a block diagram of an electronic device 1000 according to an exemplary embodiment. Fig.10 As shown, the electronic device 1000 may include: a processor 1001 and a memory 1002. The electronic device 1000 may also include one or more of a multimedia component 1003, an input / output (I / O) interface 1004, and a communication component 1005.

[0165] The processor 1001 is used to control the overall operation of the electronic device 1000 to complete all or part of the steps in the above-mentioned wind turbine life estimation method. The memory 1002 is used to store various types of data to support the operation of the electronic device 1000, and these data may include, for example, instructions for any application or method used to operate on the electronic device 1000, and application-related data, such as the size and material of each component in the wind turbine, environmental information of the wind turbine installation location, etc. The memory 1002 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, disk or optical disk. The multimedia component 1003 may include a screen and an audio component. The screen may be, for example, a touch screen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signal may be further stored in the memory 1002 or sent via the communication component 1005. The audio component also includes at least one speaker for outputting audio signals. The I / O interface 1004 provides an interface between the processor 1001 and other interface modules, and the above-mentioned other interface modules can be a keyboard, a mouse, buttons, etc. These buttons can be virtual buttons or physical buttons. The communication component 1005 is used for wired or wireless communication between the electronic device 1000 and other devices. Wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, 4G, NB-IOT, eMTC, or other 5G, etc., or a combination of one or more of them, is not limited here. Therefore, the corresponding communication component 1005 may include: Wi-Fi module, Bluetooth module, NFC module, etc.

[0166] In an exemplary embodiment, the electronic device 1000 can be implemented by one or more application specific integrated circuits (ASIC), digital signal processors (DSP), digital signal processing devices (DSPD), programmable logic devices (PLD), field programmable gate arrays (FPGA), controllers, microcontrollers, microprocessors or other electronic components to execute the above-mentioned wind turbine life estimation method.

[0167] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided, and when the program instructions are executed by a processor, the steps of the above-mentioned wind turbine life estimation method are implemented. For example, the computer-readable storage medium can be the above-mentioned memory 1002 including program instructions, and the above-mentioned program instructions can be executed by the processor 1001 of the electronic device 1000 to complete the above-mentioned wind turbine life estimation method.

[0168] In another exemplary embodiment, a computer program product is further provided. The computer program product includes a computer program executable by a programmable device. The computer program has a code portion for executing the above-mentioned wind turbine life estimation method when executed by the programmable device.

[0169] The preferred embodiments of the present disclosure are described in detail above in conjunction with the accompanying drawings; however, the present disclosure is not limited to the specific details in the above embodiments. Within the technical concept of the present disclosure, a variety of simple modifications can be made to the technical solution of the present disclosure, and these simple modifications all fall within the protection scope of the present disclosure.

[0170] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner without contradiction. In order to avoid unnecessary repetition, the present disclosure will not further describe various possible combinations.

[0171] In addition, various embodiments of the present disclosure may be arbitrarily combined, and as long as they do not violate the concept of the present disclosure, they should also be regarded as the contents disclosed by the present disclosure.

Claims

1. A method for estimating the life of a wind turbine generator set, characterized in that: include: Acquire parameter information of a wind turbine generator set, and construct a finite element model corresponding to the wind turbine generator set according to the parameter information, wherein the parameter information includes the size and material of each component in the wind turbine generator set and the environmental data of the installation location of the wind turbine generator set; Performing multi-operating condition simulations using the finite element model to determine simulated damage at key locations in the wind turbine generator set under different operating conditions, wherein the key locations represent locations with weak strength in the wind turbine generator set; Determine the operating condition frequency of the wind turbine generator set at different machine positions, and determine the estimated life of each component of the wind turbine generator set at different machine positions according to the operating condition frequency and the simulated damage of each key position under different operating conditions; The minimum value of the estimated lifespans of the components in the wind turbine generator set is determined as the estimated lifespan of the wind turbine generator set.

2. The method for estimating the life of a wind turbine generator set according to claim 1, characterized in that: The multi-operating condition simulation is performed by the finite element model to obtain simulated damage of each key position in the wind turbine under different operating conditions, including: Determine a plurality of fatigue working conditions corresponding to the wind turbine generator set, and simulate each of the fatigue working conditions for a preset time using the finite element model to obtain a simulated stress spectrum of each of the key positions in the wind turbine generator set, wherein the simulated stress spectrum represents the simulated stress amplitude of the key position under different fatigue working conditions and the number of simulation cycles corresponding to the simulated stress amplitude; For each of the fatigue conditions, rain flow counting is performed on the simulated stress spectrum of the key position under the fatigue condition to obtain the simulated stress characteristics of the key position within the preset time period, wherein the simulated stress characteristics include a simulated stress range and a number of simulated stress changes; According to the simulated stress characteristics and the material stress-life curve of the key position, linear cumulative damage processing is performed to obtain the simulated damage of the key position under the fatigue condition.

3. The method for estimating the life of a wind turbine generator set according to claim 2, characterized in that: The linear cumulative damage processing is performed according to the simulated stress characteristics and the material stress-life curve of the key position to obtain the simulated damage of the key position under the fatigue condition, including: According to the material safety factor corresponding to the key position, the simulation stress range is corrected to obtain a simulation target stress range; Determining a simulation target stress at the key position according to the simulation stress range and the stress-life curve at the key position; According to the simulation target stress, the simulation target stress range and the number of simulation stress changes, the simulated damage of the key position under the fatigue condition is obtained.

4. The method for estimating the life of a wind turbine generator set according to claim 2, characterized in that: The finite element model is used to simulate each fatigue condition for a preset time to obtain a simulated stress spectrum of each key position in the wind turbine generator set, including: For each of the fatigue working conditions, the finite element model is used to simulate the fatigue working condition for a preset time period to obtain a simulated load time history of each of the key positions within the preset time period, wherein the simulated load time history represents a process in which the fatigue load at the key position changes with time within the preset time period; Determine the simulated stress time history corresponding to the simulated load time history of each key position according to the load-stress conversion relationship corresponding to each key position; According to the fatigue influence factors corresponding to the key positions, the stress at each moment in the stress time history corresponding to the key positions is corrected to obtain the simulated stress spectrum of the key positions.

5. The method for estimating the life of a wind turbine generator set according to claim 4, characterized in that: The load-stress conversion relationship corresponding to the key position is determined by the following method: Performing static analysis on the structure of the wind turbine generator set by using the finite element model to obtain the first simulated principal stress of each key position in the wind turbine generator set under different working conditions; The load-stress conversion relationship of the key position is determined according to the azimuth angle of each key position and the first simulated principal stress, wherein the azimuth angle represents the angle between the key position along the clockwise direction and the positive direction of the x-axis in the horizontal plane coordinate system.

6. The method for estimating the life of a wind turbine generator set according to claim 4, characterized in that: The determining the simulated stress time history corresponding to the simulated load time history of each key position according to the load-stress conversion relationship corresponding to each key position respectively includes: For each time point in the simulated load time history of each key position, determining the simulated combined bending moment and the simulated ultimate bending moment of the key position under the fatigue condition; Determine the angle between the simulated combined bending moment and the key position in the horizontal plane coordinate system, and determine, based on the load-stress conversion relationship, the second simulated principal stress and the third simulated principal stress corresponding to the azimuth angle closest to the angle, wherein the second simulated principal stress represents the principal stress of the key position under the fatigue condition, and the third simulated principal stress represents the principal stress of the key position under all fatigue conditions; According to the simulated combined bending moment, the simulated ultimate bending moment, the second simulated principal stress and the third simulated principal stress, a first simulated principal stress at the key position at the time point is obtained; The simulated stress time history of the key position is determined according to the simulated principal stress at each time point in the simulated load time history of the key position.

7. The method for estimating the life of a wind turbine generator set according to any one of claims 1 to 6, characterized in that: Determining the estimated life of each component of the wind turbine generator set at different machine positions according to the frequency of each operating condition and the simulated damage of each key position under different operating conditions includes: For each key position, according to the simulated damage of the key position and the frequency of each working condition, the intermediate simulated damage of each position of the key position is obtained; Adding the intermediate simulated damage of the key position at each position to obtain the total simulated damage of the component at the key position; The estimated life of the components is determined based on the total simulated damage of each of the components.

8. The method for estimating the life of a wind turbine generator set according to any one of claims 1 to 6, characterized in that: The determining the operating frequency of the wind turbine generator set at different machine positions includes: Determining wind resource data of the wind turbine generator set at different positions, wherein the wind resource data includes at least one of wind speed, wind direction, turbulence and wind shear; According to the wind resource data, load simulation is performed on the wind turbine generator set to obtain the operating condition frequency of the wind turbine generator set at different machine positions.

9. A wind turbine life prediction device, characterized in that: include: A modeling module is configured to obtain parameter information of a wind turbine generator set and construct a finite element model corresponding to the wind turbine generator set according to the parameter information, wherein the parameter information includes the size and material of each component in the wind turbine generator set and the environmental data of the installation location of the wind turbine generator set; A simulation module is configured to perform multi-operating condition simulation through the finite element model to determine simulated damage of each key position in the wind turbine generator set under different operating conditions, wherein the key position represents a position with weak strength in the wind turbine generator set; A first execution module is configured to determine the operating condition frequency of the wind turbine generator set at different machine positions, and determine the estimated life of each component of the wind turbine generator set at different machine positions according to each operating condition frequency and the simulated damage of each key position under different operating conditions; The second execution module is configured to determine the minimum value of the estimated lifespans of the components in the wind turbine generator set as the estimated lifespan of the wind turbine generator set.

10. The wind turbine life prediction device according to claim 9, characterized in that: The simulation module is configured to determine a plurality of fatigue working conditions corresponding to the wind turbine generator set, and simulate each of the fatigue working conditions for a preset time through the finite element model to obtain a simulated stress spectrum of each of the key positions in the wind turbine generator set, wherein the simulated stress spectrum represents the simulated stress amplitude of the key position under different fatigue working conditions and the number of simulation cycles corresponding to the simulated stress amplitude; For each of the fatigue conditions, rain flow counting is performed on the simulated stress spectrum of the key position under the fatigue condition to obtain the simulated stress characteristics of the key position within the preset time period, wherein the simulated stress characteristics include a simulated stress range and a number of simulated stress changes; According to the simulated stress characteristics and the material stress-life curve of the key position, linear cumulative damage processing is performed to obtain the simulated damage of the key position under the fatigue condition.

11. The wind turbine life prediction device according to claim 10, characterized in that: The simulation module is configured to correct the simulation stress range according to the material safety factor corresponding to the key position to obtain a simulation target stress range; Determining a simulation target stress at the key position according to the simulation stress range and the stress-life curve at the key position; According to the simulation target stress, the simulation target stress range and the number of simulation stress changes, the simulated damage of the key position under the fatigue condition is obtained.

12. The wind turbine life prediction device according to claim 10, characterized in that: The simulation module is configured to simulate the fatigue condition for a preset time period through the finite element model for each fatigue condition, and obtain a simulated load time history of each key position within the preset time period, wherein the simulated load time history represents a process in which the fatigue load at the key position changes with time within the preset time period; Determine the simulated stress time history corresponding to the simulated load time history of each key position according to the load-stress conversion relationship corresponding to each key position; According to the fatigue influence factors corresponding to the key positions, the stress at each moment in the stress time history corresponding to the key positions is corrected to obtain the simulated stress spectrum of the key positions.

13. The wind turbine life prediction device according to claim 12, characterized in that: The simulation module is configured to perform static analysis on the structure of the wind turbine generator set through the finite element model to obtain the first simulated principal stress of each key position in the wind turbine generator set under different working conditions; The load-stress conversion relationship of the key position is determined according to the azimuth angle of each key position and the first simulated principal stress, wherein the azimuth angle represents the angle between the key position along the clockwise direction and the positive direction of the x-axis in the horizontal plane coordinate system.

14. The wind turbine life prediction device according to claim 12, characterized in that: The simulation module is configured to determine the simulated combined bending moment and the simulated limit bending moment of each key position under the fatigue condition at each time point in the simulated load time history of each key position; Determine the angle between the simulated combined bending moment and the key position in the horizontal plane coordinate system, and determine, based on the load-stress conversion relationship, the second simulated principal stress and the third simulated principal stress corresponding to the azimuth angle closest to the angle, wherein the second simulated principal stress represents the principal stress of the key position under the fatigue condition, and the third simulated principal stress represents the principal stress of the key position under all fatigue conditions; According to the simulated combined bending moment, the simulated ultimate bending moment, the second simulated principal stress and the third simulated principal stress, a first simulated principal stress at the key position at the time point is obtained; The simulated stress time history of the key position is determined according to the simulated principal stress at each time point in the simulated load time history of the key position.

15. The wind turbine life prediction device according to any one of claims 9 to 14, characterized in that: The first execution module is configured to obtain, for each key position, the intermediate simulated damage of each key position according to the simulated damage of the key position and the frequency of each working condition; Adding the intermediate simulated damage of the key position at each position to obtain the total simulated damage of the component at the key position; The estimated life of the components is determined based on the total simulated damage of each of the components.

16. The wind turbine life prediction device according to any one of claims 9 to 14, characterized in that: The first execution module is configured to determine wind resource data of the wind turbine generator set at different positions, wherein the wind resource data includes at least one of wind speed, wind direction, turbulence and wind shear; According to the wind resource data, load simulation is performed on the wind turbine generator set to obtain the operating condition frequency of the wind turbine generator set at different machine positions.

17. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 8 is implemented.

18. An electronic device, characterized in that: include: a memory having a computer program stored thereon; A processor, configured to execute the computer program in the memory to implement the method according to any one of claims 1 to 8.

19. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the method according to any one of claims 1 to 8 is implemented.

Images