Wind turbine tower weld fatigue load-based operation optimization method and system

Abstract

The application discloses a kind of based on wind turbine tower cylinder weld fatigue load operation optimization method and system, this method includes: S01, obtain cabin running azimuth, cabin running azimuth is converted into absolute yaw position;Tower cylinder circumference is divided into corresponding sector, and the average cabin running azimuth in preset time is calculated;Obtain wind speed, calculate the average wind speed and effective turbulence intensity of preset time;S02, the above parameters are extended to full life cycle, the cumulative equivalent fatigue load of tower cylinder circumferential weld hot spot each sector is predicted;S03, according to the equivalent fatigue index of material, the equivalent fatigue load of tower cylinder circumferential point is calculated, and finally the maximum fatigue damage of tower cylinder circumferential point is obtained according to equivalent fatigue load-stress relationship;S04, if maximum fatigue damage is greater than the equivalent maximum fatigue damage of preset tower cylinder circumferential life cycle, then yaw sector management control is executed.The application has the advantages of reducing the fatigue damage of tower cylinder specific orientation and the like.

Description

Technical Field

[0001] This invention mainly relates to the field of wind power technology, specifically to an operation optimization method and system based on fatigue load of wind turbine tower welds. Background Technology

[0002] Wind turbine towers bear alternating fatigue loads over a 20-year lifespan. These fatigue loads primarily originate from turbulence, and the intensity of turbulence varies with wind direction due to surface roughness and wake effects. Traditional tower fatigue load calculations do not consider the influence of wind direction, assuming the wind always travels in one direction. However, at specific wind farm sites, wind turbines yaw according to the prevailing wind direction, ensuring they are always in the optimal windward orientation. Therefore, using a single-direction load for tower fatigue accumulation is conservative and cannot effectively control the turbine's operation. Summary of the Invention

[0003] The technical problem to be solved by the present invention is: in view of the problems existing in the prior art, the present invention provides an operation optimization method and system based on fatigue load of wind turbine tower weld seam to reduce fatigue damage in a specific position of the tower.

[0004] To solve the above-mentioned technical problems, the technical solution proposed by this invention is as follows:

[0005] An operation optimization method based on fatigue load of wind turbine tower welds includes the following steps:

[0006] S01. Obtain the nacelle operating azimuth angle and convert it into an absolute yaw position according to the calibrated relative position relationship; divide the circumference of the wind turbine tower into corresponding sectors and calculate the average nacelle operating azimuth angle over a preset time period.

[0007] Obtain wind speed and calculate the average wind speed Vave and effective turbulence intensity Ieff over a preset time.

[0008] S02. Extend the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff to the entire life cycle of the wind turbine to predict the hot spot P of the tower circumferential weld. j The cumulative equivalent fatigue load M in each sector pj_sector_prediction ;

[0009] S03. Based on the equivalent fatigue index m of the material, calculate the circumferential point P of the tower cylinder. j Equivalent fatigue load DEL prediction Based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles around the tower circumference were calculated, and finally, the point P on the tower circumference was obtained. j Maximum fatigue damage D max_prediction ;

[0010] S04, the maximum fatigue damage D max_prediction Compare with the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle; if D max_prediction If Dmax is reached, yaw sector management control will be executed to adjust the nacelle's yaw position and reduce the tower's circumferential point P. j Maximum fatigue damage; otherwise, continue operating according to the current operating strategy.

[0011] Preferably, the specific process of step S02 is as follows:

[0012] Based on the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff, the wind speed Weibull distribution and reference turbulence intensity Iref of the wind turbine within the operating time t are obtained.

[0013] Based on the Weibull distribution of wind speed and the reference turbulence intensity Iref within the operating time t, the equivalent fatigue load of the wind turbine with the calibrated load calculation simulation model is constructed.

[0014] Based on the principle of linear accumulation of damage, the equivalent fatigue load of the wind turbine with operating time t is extended to the cumulative equivalent fatigue load of the design life T.

[0015] Preferably, in step S04, yaw sector management control is performed based on the functional relationship between the equivalent fatigue damage of the hot spot Pj of the circumferential weld of the tower and the yaw azimuth angle of sector management.

[0016] Preferably, the process of yaw sector management and control is as follows: adjust the wind turbine to operate in the corresponding sector, that is, adjust the operating time in the wind direction θ to reduce the operating time of point Pj under high load, thereby reducing fatigue damage in the design life.

[0017] Preferably, in step S04, the specific process for obtaining the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle is as follows:

[0018] S41. Obtain the wind frequency distribution P(θ), Weibull distribution parameters (A, k) of each sector, and the effective turbulence intensity Ieff of each sector for each specific site location.

[0019] S42. Through load calculation, the nacelle azimuth angle θ and the hot spot P of the tower circumferential weld are obtained. j The corresponding equivalent fatigue load M pj (θ);

[0020] S43. Sort the wind frequency distribution P(θ) of each sector; the hot spot P of the circumferential weld of the tower. j The corresponding equivalent fatigue load M pj(θ) Sort; based on the fatigue curve index m of the material and the corresponding M after sorting. pj P(θ) and P(θ) are used to perform equivalent transformation on the equivalent fatigue load of all sectors to obtain the equivalent maximum value of the tower circumference, DEL_Max.

[0021] S44. Based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles of the tower circumference are calculated, and the equivalent maximum fatigue damage Dmax of the tower circumference life cycle is calculated.

[0022] Preferably, the specific process of step S42 is as follows: Divide the wind rose diagram into n equally divided sectors, assuming that the azimuth of the unit is at the center of each sector, then solve for the wind direction at point P. j P1, P2, ... P on the circumference of the tower corresponding to the position j …P n The equivalent fatigue load M at the point pj (θ).

[0023] Preferably, in step S43, the wind direction weights of each sector in the wind rose diagram are used to determine the values ​​of P1, P2, ... P on the circumference of the tower. j …P n The cumulative equivalent fatigue load at the point, and the hot spot P of the circumferential weld of the tower. j The equivalent transformation is performed according to the following formula:

[0024]

[0025] This invention also discloses an operation optimization system based on fatigue load of wind turbine tower welds, comprising:

[0026] The first module is used to obtain the nacelle operating azimuth angle and convert it into an absolute yaw position according to the calibrated relative position relationship; divide the circumference of the wind turbine tower into corresponding sectors and calculate the average nacelle operating azimuth angle over a preset time period.

[0027] Used to obtain wind speed and calculate the average wind speed Vave and effective turbulence intensity Ieff over a preset time;

[0028] The second module extends the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff to the entire life cycle of the wind turbine, predicting the hot spot P of the tower circumferential weld. j The cumulative equivalent fatigue load M in each sector pj_sector_prediction ;

[0029] The third module is used to calculate the circumferential point P of the tower cylinder based on the material's equivalent fatigue index m. j Equivalent fatigue load DEL predictionBased on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles around the tower circumference were calculated, and finally, the point P on the tower circumference was obtained. j Maximum fatigue damage D max_prediction ;

[0030] The fourth module is used to calculate the maximum fatigue damage D. max_prediction Compare with the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle; if D max_prediction If Dmax is reached, yaw sector management control will be executed to adjust the nacelle's yaw position and reduce the tower's circumferential point P. j Maximum fatigue damage; otherwise, continue operating according to the current operating strategy.

[0031] The present invention further discloses a computer-readable storage medium having a computer program stored thereon, characterized in that the computer program, when run by a processor, executes the steps of the operation optimization method based on the fatigue load of the weld seam of the wind turbine tower as described above.

[0032] The present invention also discloses a computer device, including a memory and a processor, wherein the memory stores a computer program, and the computer program, when run by the processor, executes the steps of the operation optimization method based on the fatigue load of the weld seam of the wind turbine tower as described above.

[0033] Compared with the prior art, the advantages of the present invention are as follows:

[0034] 1. This invention corrects the wind direction frequency distribution design calculation value and the yaw azimuth angle of the unit operation by using the micro-site selection of a specific site. It determines whether the fatigue life of the unit tower is within the design range, or performs online fatigue damage accumulation calculation based on the principle of managing the remaining fatigue life and minimizing power generation loss, to obtain the operating sector with the optimal yaw azimuth angle. If it is not within the operating design range, the operating sector of the wind turbine is optimized, the optimal unit operating yaw azimuth angle is adjusted, fatigue damage in a specific position of the tower is reduced, and the system is optimized.

[0035] 2. This invention divides the fatigue hotspots of the tower circumferential welds according to the wind rose diagram, which can be divided into any number of equal parts (16 equal parts in the example), to obtain an equivalent fatigue load for the tower circumferential weld hotspots corresponding to different nacelle azimuth angles; this invention sorts the equivalent loads of the tower circumferential weld fatigue hotspots and the occurrence time probabilities of each sector according to their magnitude, to obtain an optimized sequence of equivalent fatigue loads and occurrence time probabilities for each sector, and thus obtains the maximum equivalent fatigue load of the tower circumference.

[0036] 3. This invention employs tower fatigue load accumulation based on wind direction frequency distribution, which can effectively reduce the equivalent fatigue load of tower welds, thereby reducing the weight of the wind turbine tower, saving steel consumption, and ultimately lowering the investment cost of the wind turbine. The prediction of tower fatigue weld life allows for early or real-time monitoring of the tower's fatigue life status. For the entire wind farm, the output status (generating hours) and remaining life of wind turbines at different turbine locations are not equal. Therefore, life prediction based on turbine location allows for output control of all turbines across the farm. For example, under power curtailment conditions, turbines with lower predicted lifespans can be prioritized for shutdown to optimize the lifespan and power generation of all wind turbines. Early prediction of tower lifespan, and if the predicted lifespan exceeds the design lifespan, allows for sector management of tower lifespan, reducing fatigue damage at circumferential danger points Pj.

[0037] 4. This invention introduces wind speed and yaw azimuth sensor data from the nacelle anemometer and synchronizes the data clock to obtain the average yaw azimuth angle, average wind speed, and turbulence intensity over 10 minutes. It then constructs a functional relationship expression between the fatigue damage of the hot spot of the tower circumferential weld and the yaw azimuth angle, thereby realizing the operation sector management of the unit and reducing the equivalent fatigue load at the location of the maximum damage to the tower weld. Attached Figure Description

[0038] Figure 1 M is the tower load coordinate system in this invention. x M y Schematic diagram.

[0039] Figure 2 This is a flowchart illustrating the cumulative fatigue load of the wind turbine tower weld in this invention.

[0040] Figure 3 This is a schematic diagram of wind direction distribution in this invention.

[0041] Figure 4 The tower circumference P of the present invention j Figure showing the cumulative results of point equivalent fatigue load.

[0042] Figure 5 This is a flowchart of the method of the present invention in a specific embodiment. Detailed Implementation

[0043] The present invention will be further described below with reference to the accompanying drawings and specific embodiments.

[0044] like Figure 1 As shown in the figure, the operation optimization method based on the fatigue load of the wind turbine tower weld in this embodiment includes the following steps:

[0045] S01. Obtain the nacelle operating azimuth angle and convert it into an absolute yaw position according to the calibrated relative position relationship; divide the circumference of the wind turbine tower into corresponding sectors and calculate the average nacelle operating azimuth angle over a preset time period.

[0046] Obtain wind speed and calculate the average wind speed Vave and effective turbulence intensity Ieff over a preset time.

[0047] S02. Extend the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff to the entire life cycle of the wind turbine to predict the hot spot P of the tower circumferential weld. j The cumulative equivalent fatigue load M in each sector pj_sector_prediction ;

[0048] S03. Based on the equivalent fatigue index m of the material, calculate the circumferential point P of the tower cylinder. j Equivalent fatigue load DEL prediction Based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles around the tower circumference were calculated, and finally, the point P on the tower circumference was obtained. j Maximum fatigue damage D max_prediction ;

[0049] S04, the maximum fatigue damage D max_prediction Compare with the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle; if D max_prediction If Dmax is reached, yaw sector management control will be executed to adjust the nacelle's yaw position and reduce the tower's circumferential point P. j Maximum fatigue damage; otherwise, continue operating according to the current operating strategy.

[0050] This invention uses the wind direction frequency distribution design calculation value and the yaw azimuth angle of the unit operation to correct the wind turbine's tower fatigue life within the design range, based on the micro-site selection of a specific location. Alternatively, based on the principles of managing the remaining fatigue life and minimizing power generation loss, it performs online fatigue damage accumulation calculation to obtain the optimal yaw azimuth angle operating sector. If it is not within the design range, the operating sector of the wind turbine is optimized, the optimal unit operating yaw azimuth angle is adjusted, fatigue damage in a specific azimuth of the tower is reduced, and the system is optimized.

[0051] In one specific embodiment, the process of obtaining the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle is as follows:

[0052] S41. Obtain the wind frequency distribution P(θ), Weibull distribution parameters (A, k) of each sector, and the effective turbulence intensity Ieff of each sector for each specific site location.

[0053] S42. Through load calculation, the nacelle azimuth angle θ and the hot spot P of the tower circumferential weld are obtained. j The corresponding equivalent fatigue load M pj (θ);

[0054] S43. Sort the wind frequency distribution P(θ) of each sector; the hot spot P of the circumferential weld of the tower. j The corresponding equivalent fatigue load M pj (θ) Sort; based on the fatigue curve index m of the material and the corresponding M after sorting. pj P(θ) and P(θ) are used to perform equivalent transformation on the equivalent fatigue load of all sectors to obtain the equivalent maximum value of the tower circumference, DEL_Max.

[0055] S44. Based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles of the tower circumference are calculated, and the equivalent maximum fatigue damage Dmax of the tower circumference life cycle is calculated.

[0056] Specifically, in step S41, the micro-site selection results of each turbine site in the wind farm are obtained based on fluid dynamics simulation, and the wind frequency distribution P(θ), Weibull distribution parameters (A, k) of each sector and the effective turbulence intensity Ieff of each sector are obtained for each turbine site at a specific site.

[0057] In step S42, the hot spot P of the circumferential weld of the tower is... j Point load conversion relationship and sector division method, such as Figure 1 As shown; the wind rose diagram is divided into 16 equal sectors. Assuming the unit's azimuth angle is at the center of each sector, the wind direction at point P is then calculated. j P1, P2, ... P on the circumference of the tower corresponding to the position j …P 16 The equivalent fatigue load M at the point pj (θ); Finally, based on the wind direction weighting of each sector in the wind rose diagram, P1, P2, ... P on the tower circumference are calculated. j …P 16 Accumulated equivalent fatigue load at points;

[0058] Among them, the tower circumference P j Point P3, in the prevailing wind direction, according to the IEC 61400-1 standard specification, the tower load coordinate system M x M y Show opinions Figure 1 As shown. Where P j Based on the wind rose diagram, the sectors are typically divided into 16 sectors, namely P1, P2, P3…P… j …P 16 M x : Tower bending moment (left and right); M y: Tower front and rear bending moments; M pj Point P in the circumferential direction j The projected bending moment; θ: the angle between the wind direction and the north direction of the tower; α: the angle at point P in the circumferential direction. j The angle with due north; according to Figure 1 M x The direction is consistent with the direction of the incoming wind speed, for point P on the circumference of the tower. j The bending moment is based on M x and M y Obtained through projection transformation. Coordinate system load M x Circumferential direction P j The included angle β between the points is obtained through angle conversion, as shown in the following formula:

[0059] β=π+θ-α (1)

[0060] β: Coordinate system load M x Circumferential direction P j The angle between the points.

[0061] For any point P on the circumference j The corresponding bending moment can be obtained through coordinate system projection transformation, and the transformation formula is as follows:

[0062] M pj (β)=M x cos(β)+M y sin(β) (2)

[0063] In step S43, the wind frequency distribution P(θ) of the 16 sectors is sorted from largest to smallest; the hot spot P of the circumferential weld of the tower is... j The corresponding equivalent fatigue load M pj (θ) Sort from largest to smallest; based on the fatigue curve index m of the material and the corresponding M after sorting. pj P(θ) and P(θ) are used to perform equivalent transformation on the equivalent fatigue load of all sectors according to the formula, to obtain the equivalent maximum value of the tower circumference DEL_Max; where P1, P2, ... P on the tower circumference are determined according to the wind direction weight of each sector in the wind rose diagram. j …P 16 The cumulative equivalent fatigue load at the point, and the hot spot P of the circumferential weld of the tower. j The equivalent transformation is performed according to the following formula:

[0064]

[0065] Finally, in step S44, based on the equivalent fatigue load-stress relationship, the maximum stress range and cycle number (Si, ni) of the tower circumference are calculated, and the equivalent maximum fatigue damage Dmax of the tower circumference life cycle is calculated.

[0066] This invention divides the fatigue hotspots of the tower circumferential welds according to the wind rose diagram. The division can be into any number of equal parts (16 parts in the example), resulting in an equivalent fatigue load for the tower circumferential weld hotspots corresponding to different nacelle azimuth angles. Based on the equivalent load of the tower circumferential weld fatigue hotspots and the occurrence time probability of each sector, this invention sorts them by size to obtain an optimized sequence of equivalent fatigue loads and occurrence time probabilities for each sector, thereby obtaining the maximum equivalent fatigue load of the tower circumference.

[0067] This invention employs tower fatigue load accumulation based on wind direction frequency distribution, which can effectively reduce the equivalent fatigue load on tower welds, thereby reducing the weight of the wind turbine tower, saving steel consumption, and ultimately lowering the investment cost of the wind turbine. The prediction of tower fatigue weld life allows for early or real-time monitoring of the tower's fatigue life status. Since the output status (generating hours) and remaining life of wind turbines at different locations within a wind farm are not equal, lifespan prediction based on turbine location allows for output control across the entire farm. For example, during power curtailment, turbines with lower predicted lifespans can be prioritized for shutdown to optimize the overall lifespan and power generation of the wind turbines. Early prediction of tower lifespan, and if the predicted lifespan exceeds the design lifespan, allows for sector management of tower lifespan, reducing fatigue damage at circumferential danger points Pj.

[0068] The present invention will be further described below with reference to a specific complete embodiment:

[0069] S01. Collect the azimuth angle of the nacelle operation using the nacelle yaw position sensor, and convert the nacelle azimuth angle collected by the calibrated relative position relationship into the absolute yaw position; divide the circumference of the wind turbine tower into corresponding sectors, and calculate the average nacelle operation azimuth angle within 10 minutes.

[0070] Wind speed is collected by the cabin wind speed sensor and corrected. The wind speed is averaged every 10 minutes, and the average wind speed Vave and effective turbulence intensity Ieff are calculated every 10 minutes. The time period is synchronized with the cabin yaw azimuth clock.

[0071] S02. Extend the data from the above-mentioned units' operation to the entire life cycle to predict the hot spot P of the tower's circumferential weld. j The cumulative equivalent fatigue load M in each sector pj_sector_prediction ;

[0072] Specifically, the wind speed Weibull distribution and reference turbulence intensity Iref of the wind turbine within the operating time t are obtained based on the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff.

[0073] Based on the Weibull distribution of wind speed and the reference turbulence intensity Iref within the operating time t, the equivalent fatigue load of the wind turbine with the calibrated load calculation simulation model is constructed.

[0074] Based on the principle of linear accumulation of damage, the equivalent fatigue load of the wind turbine that has been in operation for time t is extended to the cumulative equivalent fatigue load of the design life T (e.g., 20 years).

[0075] Table 1 P j Point equivalent fatigue load calculation results

[0076]

[0077] The wind rose diagram with running time t is shown below. Figure 3 As shown in Table 1, the equivalent fatigue load at point Pj on the circumference extended to the design life T (20 years) is also shown.

[0078] S03. Based on the equivalent fatigue index m of the material, the circumference P of the tower cylinder is calculated. j Equivalent fatigue load DEL prediction Based on the equivalent fatigue load-stress relationship, the maximum stress range and cycle number (Si, ni) of the tower circumference are calculated, and the point P on the tower circumference is obtained. j Maximum fatigue damage D max_prediction ;

[0079] S04. Compare the predicted maximum fatigue damage with the equivalent maximum fatigue damage Dmax over the tower's circumferential lifespan; if the predicted maximum fatigue damage Dmax... max_prediction If the equivalent maximum fatigue damage Dmax is greater than the tower circumferential lifespan, then based on the constructed functional relationship F(DEL_opt,θ) between the equivalent fatigue damage Pj of the tower circumferential weld hotspot and the yaw azimuth angle of sector management, yaw sector management control is executed to adjust the nacelle's yaw position, thereby reducing the tower circumferential point P. j The purpose of maximum fatigue damage; if the predicted maximum fatigue damage D max_prediction If the damage is less than or equal to the equivalent maximum fatigue damage Dmax of the tower circumferential lifespan, the unit will continue to operate according to the current operating strategy.

[0080] The results of the functional relationship F(DEL_opt,θ) are shown in Table 1. Since the equivalent fatigue load (fourth column of Table 1) at point Pj on the tower circumference is different at different wind directions θ, a functional relationship F(DEL_opt,θ) between the equivalent fatigue damage of Pj and the sector management yaw azimuth angle is constructed. That is to say, if the damage at point Pj exceeds the set value, the wind turbine will adjust the operation of one sector according to the calculation in Table 1 (that is, adjust the operation time at wind direction θ, thereby reducing the operation time of point Pj under high load, thereby reducing the fatigue damage of the design life).

[0081] This invention introduces wind speed and yaw azimuth sensor data from the nacelle anemometer and synchronizes the data clock to obtain the average yaw azimuth angle, average wind speed, and turbulence intensity over 10 minutes. It then constructs a functional relationship expression between the fatigue damage of the hot spot of the tower circumferential weld and the yaw azimuth angle, thereby realizing the operation sector management of the unit and reducing the equivalent fatigue load at the location of the maximum damage to the tower weld.

[0082] This invention employs site-specific tower fatigue load prediction to determine whether the tower fatigue life of the unit is operating within the design range. Based on the functional relationship between the equivalent fatigue load of the tower circumferential weld and the yaw azimuth angle, fatigue damage accumulation is calculated online. When the predicted equivalent fatigue load exceeds the design value, sector management control is executed to optimize the operating sector of the wind turbine, adjust the optimal yaw azimuth angle of the unit, reduce fatigue damage in specific tower locations, obtain system optimization, and improve the reliability of the wind turbine.

[0083] This invention also discloses an operation optimization system based on fatigue load of wind turbine tower welds, comprising:

[0084] The first module is used to obtain the nacelle operating azimuth angle and convert it into an absolute yaw position according to the calibrated relative position relationship; divide the circumference of the wind turbine tower into corresponding sectors and calculate the average nacelle operating azimuth angle over a preset time period.

[0085] Used to obtain wind speed and calculate the average wind speed Vave and effective turbulence intensity Ieff over a preset time;

[0086] The second module extends the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff to the entire life cycle of the wind turbine, predicting the hot spot P of the tower circumferential weld. j The cumulative equivalent fatigue load M in each sector pj_sector_prediction ;

[0087] The third module is used to calculate the circumferential point P of the tower cylinder based on the material's equivalent fatigue index m. jEquivalent fatigue load DEL prediction Based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles around the tower circumference were calculated, and finally, the point P on the tower circumference was obtained. j Maximum fatigue damage D max_prediction ;

[0088] The fourth module is used to calculate the maximum fatigue damage D. max_prediction Compare with the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle; if D max_prediction If Dmax is reached, yaw sector management control will be executed to adjust the nacelle's yaw position and reduce the tower's circumferential point P. j Maximum fatigue damage; otherwise, continue operating according to the current operating strategy.

[0089] The optimization system of the present invention is used to perform the optimization method described above, and also has the advantages described above.

[0090] This invention further discloses a computer-readable storage medium storing a computer program. When executed by a processor, the computer program performs the steps of the operation optimization method based on fatigue load of wind turbine tower welds as described above. This invention also discloses a computer device including a memory and a processor. The memory stores a computer program, which, when executed by a processor, performs the steps of the operation optimization method based on fatigue load of wind turbine tower welds as described above. This invention can implement all or part of the processes in the methods of the above embodiments, or it can be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include: any entity or device capable of carrying computer program code, recording media, USB flash drives, portable hard drives, magnetic disks, optical disks, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signals, telecommunication signals, and software distribution media, etc. Memory can be used to store computer programs and / or modules. The processor performs various functions by running or executing the computer programs and / or modules stored in the memory, and by accessing data stored in the memory. Memory may include high-speed random access memory, and may also include non-volatile memory, such as hard disks, RAM, plug-in hard disks, smart media cards (SMC), secure digital cards (SD cards), flash cards, at least one disk storage device, flash memory device, or other volatile solid-state storage devices.

[0091] The above are merely preferred embodiments of the present invention. The scope of protection of the present invention is not limited to the above embodiments. All technical solutions falling within the scope of the present invention's concept are within the scope of protection of the present invention. It should be noted that for those skilled in the art, any improvements and modifications made without departing from the principles of the present invention should be considered within the scope of protection of the present invention.

Claims

1. An operation optimization method based on fatigue load of wind turbine tower welds, characterized in that, Including the following steps: S01. Obtain the nacelle operating azimuth angle and convert it into an absolute yaw position according to the calibrated relative position relationship; divide the circumference of the wind turbine tower into corresponding sectors and calculate the average nacelle operating azimuth angle over a preset time period. Obtain wind speed and calculate the average wind speed Vave and effective turbulence intensity Ieff over a preset time. S02. The average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff are extended to the entire life cycle of the wind turbine to predict the hot spots of the tower circumferential welds. The cumulative equivalent fatigue load in each sector; S03. Based on the equivalent fatigue index m of the material, the circumferential points of the tower cylinder are calculated. The equivalent fatigue load is calculated, and based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles around the tower circumference are obtained, ultimately yielding the tower circumference point. Maximum fatigue damage ; S04, Maximum fatigue damage Compare with the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle; if If Dmax is reached, yaw sector management control will be executed to adjust the nacelle's yaw position in order to lower the tower's circumferential point. Maximum fatigue damage; otherwise, continue operating according to the current operating strategy.

2. The operation optimization method based on fatigue load of wind turbine tower weld seams according to claim 1, characterized in that, The specific process of step S02 is as follows: Based on the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff, the wind speed Weibull distribution and reference turbulence intensity Iref of the wind turbine within the operating time t are obtained. Based on the Weibull distribution of wind speed and the reference turbulence intensity Iref within the operating time t, the equivalent fatigue load of the wind turbine with the calibrated load calculation simulation model is constructed. Based on the principle of linear accumulation of damage, the equivalent fatigue load of the wind turbine with operating time t is extended to the cumulative equivalent fatigue load of the design life T.

3. The operation optimization method based on fatigue load of wind turbine tower weld seams according to claim 1, characterized in that, In step S04, yaw sector management control is performed based on the functional relationship between the equivalent fatigue damage of the hot spot Pj of the circumferential weld of the tower and the yaw azimuth angle of sector management.

4. The operation optimization method based on fatigue load of wind turbine tower weld seams according to claim 3, characterized in that, The process of yaw sector management and control involves adjusting the wind turbine to operate in the corresponding sector, that is, adjusting it to the wind direction. The running time is reduced to decrease the running time of point Pj under high load, thereby reducing fatigue damage during the design life.

5. The operation optimization method based on fatigue load of wind turbine tower weld seams according to any one of claims 1 to 4, characterized in that, In step S04, the specific process of obtaining the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle is as follows: S41. Obtain the wind frequency distribution of each sector at a specific site. Weibull distribution parameters (A, k) and effective turbulence intensity Ieff of each sector; S42. Obtain the cabin azimuth angle through load calculation. Hot spots of circumferential welds in tower Corresponding equivalent fatigue load ; S43. Distribute the wind frequency of each sector. Sort by value from largest to smallest; Identify hot spots of the tower's circumferential weld seams. Corresponding equivalent fatigue load Sort by numerical value from largest to smallest; based on the fatigue curve index m of the material and the corresponding sorted values. and The equivalent fatigue load of all sectors is equivalently transformed to obtain the equivalent maximum value of the tower circumference, DEL_Max. S44. Based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles of the tower circumference are calculated, and the equivalent maximum fatigue damage Dmax of the tower circumference life cycle is calculated.

6. The operation optimization method based on fatigue load of wind turbine tower weld seams according to claim 5, characterized in that, The specific process of step S42 is as follows: Divide the wind rose diagram into n equally divided sectors. Assuming that the azimuth of the unit is at the center of each sector, then solve for the wind direction in each sector. On the circumference of the tower corresponding to the position , , Point equivalent fatigue load .

7. The operation optimization method based on fatigue load of wind turbine tower weld seams according to claim 6, characterized in that, In step S43, the wind direction weight of each sector in the wind rose diagram is used to adjust the wind direction around the tower circumference. , , Accumulated equivalent fatigue load at points, hot spots of the circumferential weld of the tower. The equivalent transformation is performed according to the following formula: 。 8. An operation optimization system based on fatigue load of wind turbine tower welds, characterized in that, include: The first module is used to obtain the cabin operating azimuth angle and convert the cabin operating azimuth angle into an absolute yaw position based on the calibrated relative position relationship. Divide the circumference of the wind turbine tower into corresponding sectors and calculate the average nacelle azimuth angle over a preset time period; Used to obtain wind speed and calculate the average wind speed Vave and effective turbulence intensity Ieff over a preset time; The second module extends the average nacelle operating azimuth angle, average wind speed Vave, and effective turbulence intensity Ieff to the entire life cycle of the wind turbine, predicting hot spots at the circumferential welds of the tower. Cumulative equivalent fatigue load in each sector ; The third module is used to calculate the tower circumference point based on the material's equivalent fatigue index m. Equivalent fatigue load Based on the equivalent fatigue load-stress relationship, the maximum stress range and number of cycles around the tower circumference were calculated, ultimately yielding the tower circumference point. Maximum fatigue damage ; The fourth module is used to determine the maximum fatigue damage. Compare with the equivalent maximum fatigue damage Dmax of the preset tower circumferential life cycle; if If Dmax is reached, yaw sector management control will be executed to adjust the nacelle's yaw position in order to lower the tower's circumferential point. Maximum fatigue damage; otherwise, continue operating according to the current operating strategy.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is run by the processor, it executes the steps of the operation optimization method based on fatigue load of wind turbine tower weld as described in any one of claims 1 to 7.

10. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the computer program is run by the processor, it executes the steps of the operation optimization method based on fatigue load of wind turbine tower weld as described in any one of claims 1 to 7.

Images