Structural dynamic analysis method of offshore wind turbine based on typhoon whole process working condition

Abstract

This invention relates to a method for structural dynamic analysis of offshore wind turbines based on the entire typhoon operating conditions. It is applicable to the field of offshore wind power technology in new energy engineering. The technical solution adopted in this invention is as follows: Based on the user-selected design typhoon wind field parameters, the average wind speed at the turbine location during the period from entering to leaving the maximum influence radius of the typhoon is calculated. The average wind speed at each typical point on the turbine is determined based on the height of each typical point on the turbine and the average wind speed at the turbine reference point. The instantaneous wind speed V(z,t) at each typical point on the turbine is determined based on the average wind speed and the fluctuating wind speed v(z,t) at each typical point on the turbine. The wind direction angle θ(t) at the turbine location during the period from entering to leaving the maximum influence radius of the typhoon is calculated based on the design typhoon wind field parameters. A finite element model of the turbine structure is established, and dynamic time history calculations are performed by combining the instantaneous typhoon loads corresponding to each typical point on the turbine.

Description

TECHNICAL FIELD

[0001] The application relates to a method for structural dynamic analysis of an offshore wind turbine based on a whole-process typhoon working condition. BACKGROUND

[0002] With the exhaustion of traditional fossil energy such as oil and coal, the development and utilization of new energy have become the key driving force for economic development and social progress. Offshore wind power structure is a new, efficient and clean form of new energy utilization. Compared with onshore wind power, offshore wind power is closer to the southeast coastal areas of China where energy consumption is huge, which can significantly reduce the related costs and power line consumption of west-to-east power transmission; offshore wind speed is larger and turbulence intensity is smaller, so the wind resource quality is high; the valuable land resources are not occupied, and noise and visual interference are not caused.

[0003] China's coastal areas are one of the regions where typhoons land most frequently in the world. According to statistics of the China Meteorological Department, the number of typhoons and super typhoons landing in the southeast coastal areas of China is about 9 on average every year. Although considerable power generation can be obtained during the typhoon period, the extreme typhoon load also poses a serious threat to the development and construction of offshore wind power in the region.

[0004] The dynamic calculation of the wind turbine structure based on the typhoon working condition is an important aspect of the safety of the wind turbine, and the key lies in the simulation of the typhoon wind field load. The current analysis method is generally based on the design typhoon wind speed to simulate the flow field and calculate the structure under specific working conditions. The Chinese patent with the publication number CN107506521A, the publication date of December 22, 2017 and the invention name of a three-dimensional wind load simulation method for a power transmission tower discloses a numerical simulation method for simulating a spatially correlated three-dimensional wind field by using a double POD model and a Monte Carlo method; the disadvantage is that the time history characteristics of the wind field parameters cannot be obtained, and it is difficult to be directly used for structural dynamic calculation; the Chinese patent with the publication number CN108563825A, the publication date of September 21, 2018 and the invention name of a simulation method and system for a typhoon wind field discloses a method for determining the wind field profile characteristics and simulating the wind field according to the wind field characteristic parameters; the disadvantage is that additional wind tunnel simulation tests are needed, which is costly; the Chinese patent with the publication number CN11278255A, the publication date of May 4, 2021 and the invention name of a typhoon wind speed calculation method for a power transmission tower based on a typhoon model discloses a method for simulating the typhoon wind speed and wind direction of a power transmission line based on the monitored wind speed of a power transmission tower; the disadvantage is that detailed monitoring data are needed, and it is difficult to be popularized and applied in offshore environment.

[0005] Overall, the existing typhoon working condition wind field load simulation and dynamic calculation mainly exist the following problems: (1) the wind speed and wind direction change rapidly during the typhoon passing, and its fluctuating wind characteristics are significant, and the error of directly calculating the wind load according to the specification method is large; (2) due to the limitation of the bad weather conditions during the typhoon passing, it is difficult to obtain the field typhoon measured data, especially the typhoon measured data in the marine environment; (3) the numerical simulation method based on the flow field is large in work load and complex in process, and is high in requirements for calculation software and hardware; (4) the current dynamic calculation method of the wind turbine structure is generally for a specific wind speed working condition, and lacks the overall analysis of the wind turbine structure response during the whole process of the typhoon passing. SUMMARY

[0006] The technical problem to be solved by the present application is to provide a wind turbine structure dynamic analysis method based on a whole process working condition of a typhoon.

[0007] The technical scheme adopted by the present application is: a wind turbine structure dynamic analysis method based on a whole process working condition of a typhoon, characterized in that:

[0008] Based on the design typhoon wind field parameters selected by the user, the average wind speed of the wind turbine at the position of the wind turbine during the period from entering to leaving the maximum influence radius range of the typhoon is calculated ;

[0009] Based on the height of each typical point on the wind turbine and the average wind speed at the reference point of the wind turbine , the average wind speed of each typical point on the wind turbine is determined ;

[0010] Taking the von Karman power spectrum selected by the user as the standard spectrum, the fluctuating wind speed of each typical point on the wind turbine is calculated by using the harmonic superposition method ;

[0011] Based on the average wind speed and the fluctuating wind speed of each typical point on the wind turbine , the instantaneous wind speed of each typical point on the wind turbine is determined ;

[0012] Based on the design typhoon wind field parameters, the wind direction angle θ(t) of the wind turbine during the period from entering to leaving the maximum influence radius range of the typhoon is calculated

[0013] The instantaneous wind speed of each typical point on the wind turbine is vertically decomposed according to the wind direction angle θ(t), and the instantaneous wind speed of each typical point in the orthogonal direction is obtained ; ;

[0014] Based on the instantaneous wind speed of each typical point on the wind turbine in the orthogonal direction ​​and Calculate the instantaneous typhoon load at each typical point on the wind turbine;

[0015] A finite element model of the wind turbine structure was established, and dynamic time history calculations were carried out by combining the instantaneous typhoon loads of the corresponding parts of each typical point on the wind turbine.

[0016] Based on the user-selected design typhoon wind field parameters, the average wind speed at the wind turbine's location is calculated during the period from when the turbine enters to when it leaves the typhoon's maximum influence radius. ,include:

[0017] Obtain the wind field parameters of the user-selected design typhoon, including the maximum average wind speed of the design typhoon. and the corresponding maximum wind speed radius R1, the design typhoon's maximum influence radius R2 and the corresponding cut-out wind speed And the direction and speed of the typhoon's center path, V. w ;

[0018] Based on the typhoon's maximum influence radius R2, the vertical distance L from the wind turbine to the typhoon's center, and the typhoon's center's moving speed V. w The distance S(t) between the wind turbine and the center of the typhoon was calculated during the period from when the wind turbine entered the typhoon's maximum influence radius to when it left the typhoon's maximum influence radius.

[0019] Based on the distance S(t) between the wind turbine and the typhoon center, the maximum wind speed radius R1 of the designed typhoon, and the maximum average wind speed of the wind field Determine the average wind speed at the wind turbine location during the period from when the turbine enters to when it leaves the typhoon's maximum influence radius. .

[0020] The distance S(t) between the wind turbine and the typhoon center is...

[0021]

[0022] The average wind speed at the location of the fan for:

[0023] .

[0024] The height of typical points on the wind turbine and the average wind speed at the wind turbine reference point are used as the basis. Determine the average wind speed at typical points on the fan. ,include:

[0025]

[0026] In the formula, z is the height above sea level; This is the wind speed profile index.

[0027] The user-selected von Karman power spectrum is taken as a standard spectrum, and the fluctuating wind speed of each typical point on the fan is calculated , comprising:

[0028] (1) The von Karman power spectrum is selected as a standard spectrum, and the power spectrum density function is calculated according to the following formula:

[0029]

[0030] In the formula, is the circular frequency of the fluctuating wind, is the average wind speed of the first 10-minute period when the fan enters the typhoon influence range; is the friction velocity, is the longitudinal turbulence integral scale;

[0031]

[0032]

[0033]

[0034] In the formula, z0 is the surface roughness length;

[0035] (2) According to the self-power spectrum of each point , the cross-power spectrum is calculated according to the following formula :

[0036]

[0037] In the formula, is the coherence function, which can be obtained from the relative position coordinates of each point;

[0038] (3) Discretize the power spectrum, divide it into N parts in the frequency domain range 0~ fmax , and the frequency increment is

[0039]

[0040] Among them, is the upper limit cutoff frequency;

[0041] N discrete power spectrum matrices are obtained;

[0042] (4) Chlesky decomposition is performed on the power spectrum matrix to obtain a lower triangular matrix ;

[0043]

[0044] (5) Based on the lower triangular matrix Harmonic superposition was performed to obtain the 10-minute fluctuating wind speed time history at each typical point:

[0045]

[0046] In the formula: for phase angle, A random number between 0 and 2π;

[0047] (6) Repeat steps (1) to (5) K times to obtain the instantaneous wind speed time history of the wind field during the entire process of the typhoon;

[0048]

[0049] In the formula, T is the time required for the wind turbine to go from entering to leaving the typhoon's maximum influence radius; T0 is 10 minutes.

[0050] The average wind speed based on typical points on the wind turbine and pulsating wind speed Determine the instantaneous wind speed at typical points on the wind turbine. ,include:

[0051] .

[0052] The wind direction angle θ(t) of the wind turbine's position during the period from entering to leaving the maximum influence radius of the typhoon, calculated based on the design typhoon wind field parameters, includes:

[0053] .

[0054] A structural dynamic analysis device for offshore wind turbines based on the entire typhoon process, characterized in that it includes:

[0055] The location wind speed calculation module is used to calculate the average wind speed at the wind turbine location from the time the turbine enters to the time it leaves the maximum influence radius of the typhoon, based on the user-selected design typhoon wind field parameters. ;

[0056] The average wind speed calculation module is used to calculate the average wind speed based on the height of typical points on the wind turbine and the average wind speed at the wind turbine reference point. Determine the average wind speed at typical points on the fan. ;

[0057] The fluctuating wind speed calculation module is used to simulate and calculate the fluctuating wind speed at typical points on the wind turbine using a user-selected von Karman power spectrum as the standard spectrum and the harmonic superposition method. ;

[0058] The instantaneous wind speed determination module is used to determine the average wind speed at typical points on the wind turbine. and pulsating wind speed Determine the instantaneous wind speed at typical points on the wind turbine. ;

[0059] The wind direction angle calculation module is used to calculate the wind direction angle θ(t) of the wind turbine position during the period from entering to leaving the maximum influence radius of the typhoon, based on the design typhoon wind field parameters.

[0060] The wind speed decomposition module is used to decompose the instantaneous wind speed at typical points on the wind turbine. By performing vertical decomposition based on the wind direction angle θ(t), the instantaneous wind speed in the orthogonal direction at each typical point is obtained. and ;

[0061] The load calculation module is used to calculate the instantaneous wind speed in the orthogonal direction at typical points on the wind turbine. and Calculate the instantaneous typhoon load at each typical point on the wind turbine;

[0062] The model calculation module is used to establish a finite element model of the wind turbine structure and perform dynamic time history calculations by combining the instantaneous typhoon loads of the corresponding parts of each typical point on the wind turbine.

[0063] A storage medium storing a computer program executable by a processor, characterized in that: when the computer program is executed, it implements the steps of the method for structural dynamic analysis of offshore wind turbines based on the full-process operating conditions of a typhoon.

[0064] A dynamic analysis device for offshore wind turbine structure, comprising a memory and a processor, wherein the memory stores a computer program that can be executed by the processor, characterized in that: when the computer program is executed, it implements the steps of the dynamic analysis method for offshore wind turbine structure based on the full-process operating conditions of a typhoon.

[0065] The beneficial effects of this invention are as follows: This invention calculates the average wind speed at each point on the wind turbine based on the design typhoon wind field parameters, uses the selected von Karman power spectrum as the standard spectrum, and employs the harmonic superposition method for simulation to calculate the fluctuating wind speed at each typical point on the wind turbine. The instantaneous wind speed is determined by the average wind speed and the fluctuating wind speed. The wind direction angle is calculated based on the design typhoon wind field parameters. The instantaneous typhoon load on each point on the wind turbine is calculated based on the instantaneous wind speed and the wind direction angle, and model calculations are performed based on the instantaneous typhoon load.

[0066] This invention fully considers the characteristics of wind speed and direction changes in typhoon wind fields. The wind field time history simulation covers blade and tower structures above sea level, and the wind load calculation is accurate. The typhoon wind field time history simulation is specifically designed for the typhoon shutdown conditions of offshore wind turbines. The generated wind speed time history data can be directly used to calculate wind loads and then substituted into structural calculation software for calculation, which greatly improves calculation efficiency. Attached Figure Description

[0067] Figure 1 The flowchart is for an example.

[0068] Figure 2 This is a schematic diagram showing the positions of the typhoon and the wind turbine in the embodiment.

[0069] Figure 3 The variation law of typhoon with radius is designed for the example.

[0070] Figure 4 This is a schematic diagram of a typical point on the fan in the embodiment.

[0071] Figure 5 This is a schematic diagram illustrating the wind direction angle change process in the embodiment.

[0072] Figure 6 This is a schematic diagram illustrating the calculation of typhoon load on the wind turbine in the embodiment.

[0073] Figure 7 This is a specific embodiment of the example, showing the process of change in the average wind speed of a typhoon.

[0074] Figure 8 This is a specific embodiment of the example, showing the instantaneous wind speed time history of a typhoon.

[0075] Figure 9 This is a specific implementation of the embodiment, showing the process of typhoon wind direction angle change.

[0076] Figure 10 This is a three-dimensional finite element model of a wind turbine in one specific embodiment.

[0077] Figure 11 This is a cloud map showing the calculation results of the wind turbine dynamic response in a specific embodiment of the example.

[0078] 1. Wind turbine; 2. Typhoon center; 3. Path of typhoon center; 4. Vertical distance L from wind turbine to path of typhoon center; 5. Radius of maximum wind speed of typhoon R1; 6. Radius of influence of typhoon R2; 7. Sea level; 8. Wind turbine reference point 10m above sea level; 9. Wind turbine blades; 10. Wind turbine tower. Detailed Implementation

[0079] like Figure 1 As shown in the figure, this embodiment is a method for structural dynamic analysis of offshore wind turbines based on the entire typhoon process, characterized in that:

[0080] S1. Obtain the design typhoon wind field parameters selected by the user. Based on the location of the offshore wind turbine, query meteorological data to obtain historical typhoon data, and combine the wind turbine parameter requirements to determine: (1) the maximum average wind speed of the design typhoon wind field. (2) The maximum influence radius of the designed typhoon R2 and the corresponding cut-out wind speed. (3) The direction and speed of the typhoon's center path V w (4) The vertical distance L from the wind turbine to the typhoon center's path of movement, etc. (see Figure 2 ).

[0081] Due to the radial structure of the typhoon's wind field, consisting of the eye, cloud wall, and spiral rainbands, the wind speed exhibits a pattern of first increasing and then decreasing with the increase of the radius r, as shown in the attached figure. Figure 3 It can be simulated as a piecewise linear function using the following formula:

[0082] (1)

[0083] S2. Based on the user-selected design typhoon wind field parameters, calculate the average wind speed at the fan location during the period from entering to leaving the typhoon's maximum influence radius. .

[0084] During the movement of the typhoon's center, the average wind speed The distance S from the wind turbine to the typhoon center is mainly determined by the location of the wind turbine. Taking the moment when the wind turbine enters the typhoon's influence area as t=0, the distance S from the wind turbine to the typhoon center can be calculated using the following formula as the typhoon center moves:

[0085] (2)

[0086] In the formula, L is the vertical distance from the wind turbine to the typhoon center's path, which is less than R1; R2 is the typhoon's radius of influence; V w t represents the speed at which the typhoon center moves; t represents time, which begins when the wind turbine enters the typhoon's maximum influence radius.

[0087] The variation of distance S(t) and wind speed with radius The average wind speed can be obtained. The process of change is as follows:

[0088] (3)

[0089] S3. Based on the height of typical points on the wind turbine and the average wind speed at the wind turbine reference point. Determine the average wind speed at typical points on the fan. .

[0090] In this embodiment, the average wind speed in step S2 is the 10-minute average wind speed at a standard height of 10m. Based on the power function law of the average wind speed changing with height, the average wind speed time history of the typhoon wind field at M typical points within the offshore wind power structure range is obtained, including M1 points within the rotor range and M2 points within the tower range, as shown in the attached figure. Figure 4 As shown. The calculation formula is as follows:

[0091] (4)

[0092] In the formula, z is the height above sea level; The wind speed profile index can be taken as 0.12 for sea level.

[0093] S4. Using the user-selected von Karman power spectrum as the standard spectrum, the harmonic superposition method is employed to simulate and calculate the fluctuating wind speed at typical points on the wind turbine. .

[0094] (1) Based on the turbulence characteristics of typhoons, the von Karman power spectrum is selected as the standard spectrum, and the power spectral density function is calculated using the following formula:

[0095] (5)

[0096] In the formula, The angular frequency of the pulsating wind. The average wind speed during the first 10-minute period after the wind turbine enters the typhoon's influence area; For friction speed, The longitudinal turbulence integral scale can be calculated using the following formulas:

[0097] (6)

[0098] (7)

[0099] (8)

[0100] In the formula, z0 is the surface roughness length, which is taken as 0.005m for sea level.

[0101] (2) Based on the self-power spectrum of each simulation point Calculate the cross power spectrum using the following formula. :

[0102] (9)

[0103] In the formula, It is a coherence function, which can be obtained from the relative position coordinates of each simulation point.

[0104] (3) Discretize the power spectrum. In the frequency domain range 0~ If the interior is divided into N parts, then the frequency increment is:

[0105] (10)

[0106] in, Let be the upper cutoff frequency. Then, N discretized power spectrum matrices can be obtained. .

[0107] (4) By performing Chlesky decomposition on the power spectrum matrix, the lower triangular matrix can be obtained. .

[0108] (11)

[0109] (5) Based on the lower triangular matrix By performing harmonic superposition, as shown in the following formula, the 10-minute fluctuating wind speed time history for each simulation point can be obtained:

[0110] (12)

[0111] In the formula, for phase angle, A random number between 0 and 2π.

[0112] (6) By repeating steps (1) to (5) K times, the instantaneous wind speed time history of the wind field during the entire process of the typhoon can be obtained, where K is calculated by the following formula:

[0113] (13)

[0114] In the formula, T0 is 10 min; T is the time required for the wind turbine to go from entering to leaving the maximum influence radius of the typhoon.

[0115] S5, Average wind speed based on typical points on the wind turbine and pulsating wind speed Determine the instantaneous wind speed at typical points on the wind turbine. .

[0116] By superimposing the fluctuating wind speed time history at each point with the corresponding average wind speed, the corresponding instantaneous wind speed time history can be obtained, as shown in the following formula:

[0117] (14)

[0118] S6. Simulate the typhoon wind direction angle change process and decompose to obtain the orthogonal instantaneous wind field.

[0119] A. The wind direction angle θ(t) of the wind turbine position during the period from entering to leaving the maximum influence radius of the typhoon is calculated based on the design typhoon wind field parameters.

[0120] Influenced by the Coriolis force, the airflow of typhoons in the Northern Hemisphere exhibits a left-handed spiral pattern as it moves towards the low-pressure center, appearing as a counter-clockwise spiral when viewed from above. This results in significant changes in wind direction during the passage of a typhoon. The changes in wind direction at the location of the wind turbines as the typhoon center moves are shown in the attached figure. Figure 5 As shown.

[0121] At time t0, the typhoon center is located at position O0, and the wind turbine enters the typhoon's influence area with a wind direction angle of θ0. As the typhoon center gradually moves to O1, O2, O3, and O4, the wind direction angle at the wind turbine's location also gradually changes. At time t4, the wind turbine is about to leave the typhoon's influence area, and the wind direction angle at this time is θ4. During this process, the change in the wind direction angle at the wind turbine's location can be calculated using the following formula:

[0122] (15)

[0123] B. Calculate the instantaneous wind speed at typical points on the fan. By performing vertical decomposition based on the wind direction angle θ(t), the instantaneous wind speed in the orthogonal direction at each typical point is obtained. and As shown in the following formula:

[0124] (16)

[0125] (17)

[0126] S7. Instantaneous wind speed in the orthogonal direction at typical points on the wind turbine. and Calculate the instantaneous typhoon load at each typical point on the wind turbine.

[0127] This embodiment calculates the resistance experienced by the tower and blades. and the lift force on the blades As attached Figure 6 Based on the simulated wind speed points, the wind turbine is discretized into M parts. The formula for calculating the instantaneous typhoon load for each part is as follows:

[0128] (18)

[0129] (19)

[0130] In the formula, C D C represents the drag coefficient of the tower and blade sections. L The lift coefficient of the blade section; for the tower, CD The value is taken based on a smooth circular cross-section; for the blade, C... D and C L The value should be determined by referring to the airfoil parameters of the blade; V i V represents the instantaneous wind speed; A represents the wind speed perpendicular to the wind speed V. i The projected area of ​​the structure in the direction.

[0131] S8. Establish a finite element model of the wind turbine structure, and combine the instantaneous typhoon loads of the corresponding parts of each typical point on the wind turbine to carry out dynamic time history calculations to obtain the dynamic response of the wind turbine structure under the entire typhoon process.

[0132] To facilitate understanding, the following is a specific example:

[0133] (1) A certain offshore wind turbine unit 1 is located in the East China Sea and is in a shut-down, feathered state under typhoon conditions; the turbine hub height is 110m and the rotor diameter is 120m. Based on the historical typhoon data of the turbine's location, the maximum 10-minute average wind speed of the design typhoon is taken. =50m / s, reaching the level of a Category 15 strong typhoon, with a maximum wind speed radius R1=20km; cut-off wind speed =25m / s, corresponding to an influence radius R2=100km; the typhoon center's path direction is from east to west, and its moving speed is V. w =30km / h, the vertical distance L from the wind turbine to the typhoon center's moving path is 10km.

[0134] (2) Substituting the above parameters into equation (3), the change process of average wind speed over time is obtained as follows: Figure 7 As shown in the figure, the duration of the typhoon's impact is T=6.6h. During this period, as the typhoon center moves, the average wind speed exhibits an "M"-shaped double-peak trend, first increasing and then decreasing, then increasing again and decreasing once more.

[0135] (3) Select M1=73 simulation points in the impeller range and M2=4 simulation points in the tower range, and substitute the above wind speed time history into equation (4) to simulate and generate the multi-point average wind speed time history in the wind turbine range.

[0136] (4) Average wind speed Substituting parameters such as height z into equations (5) to (13), the time histories of fluctuating wind speeds at M=77 points within T=6.6h can be obtained.

[0137] (5) Average wind speed Substituting the fluctuating wind speed into equation (14), we can obtain the instantaneous wind speed time history at each point within time period T, such as... Figure 8 The figure shows the wind speed time history at the hub height, where the maximum instantaneous speed can reach 70 m / s.

[0138] (6) Calculate the change process of the wind direction angle. When the typhoon center 2 passes from west to east south of the wind turbine 1, substitute t=0~6.6h into equation (15) to calculate the change process of the wind direction angle θ(t) as follows: Figure 9 The wind direction gradually changes from northeast to southeast clockwise, with a wind direction angle change of 168° throughout the process; especially in the 2-4 hour period near the typhoon center, the wind direction angle change is about 140°. Then, the instantaneous wind speed V(t) is decomposed according to equations (16) and (17) to obtain the instantaneous wind speed time history in the two vertical directions of east-west and north-south.

[0139] (7) Calculate the typhoon instantaneous load time history of the blades and tower according to equations (18) and (19).

[0140] (8) Establish a three-dimensional finite element model of the wind turbine structure, and apply the instantaneous typhoon load to the corresponding locations in the model, as shown in the attached figure. Figure 10 Set up a dynamic time history analysis step to obtain the wind turbine's dynamic response, as shown in the attached figure. Figure 11 .

[0141] This embodiment also provides a dynamic analysis device for offshore wind turbine structures based on the entire typhoon process, including: a location wind speed calculation module, an average wind speed calculation module, a fluctuating wind speed calculation module, an instantaneous wind speed determination module, a wind direction angle calculation module, a wind speed decomposition module, a load calculation module, and a model calculation module.

[0142] In this example, the location wind speed calculation module is used to calculate the average wind speed at the reference point of the wind turbine's location during the period from when the wind turbine enters the typhoon's maximum influence radius to when it leaves the typhoon's maximum influence radius, based on the user-selected design typhoon wind field parameters. The average wind speed calculation module is used to calculate the average wind speed based on the height of typical points on the wind turbine and at a height of 10m. Determine the average wind speed at typical points on the fan. The fluctuating wind speed calculation module uses the user-selected von Karman power spectrum as the standard spectrum and employs the harmonic superposition method to simulate and calculate the fluctuating wind speed at typical points on the wind turbine. The instantaneous wind speed determination module is used to determine the average wind speed at typical points on the wind turbine. and pulsating wind speed Determine the instantaneous wind speed at typical points on the wind turbine. The wind direction angle calculation module is used to calculate the wind direction angle θ(t) of the wind turbine position during the period from entering to leaving the maximum influence radius of the typhoon, based on the design typhoon wind field parameters; the wind speed decomposition module is used to decompose the instantaneous wind speed at typical points on the wind turbine. By performing vertical decomposition based on the wind direction angle θ(t), the instantaneous wind speed in the orthogonal direction at each typical point is obtained. and The load calculation module is used to calculate the instantaneous wind speed in the orthogonal direction at typical points on the wind turbine. and The module calculates the instantaneous typhoon load at each typical point on the wind turbine. The model calculation module is used to establish a finite element model of the wind turbine structure and perform dynamic time history calculations based on the instantaneous typhoon load at each typical point on the wind turbine.

[0143] This embodiment also provides a storage medium storing a computer program that can be executed by a processor. When the computer program is executed, it implements the steps of the offshore wind turbine structure dynamic analysis method based on the entire typhoon process in this example.

[0144] This embodiment also provides a dynamic analysis device for offshore wind turbine structures, which has a memory and a processor. The memory stores a computer program that can be executed by the processor. When the computer program is executed, it implements the steps of the dynamic analysis method for offshore wind turbine structures based on the entire typhoon process in this example.

Claims

1. A method for structural dynamic analysis of offshore wind turbines based on the entire typhoon process, characterized in that: Based on the user-selected design typhoon wind field parameters, the average wind speed at the wind turbine's reference point was calculated during the period from when the turbine entered to when it left the typhoon's maximum influence radius. ; Based on the height of typical points on the wind turbine and the average wind speed at the wind turbine reference point Determine the average wind speed at typical points on the fan. ; Using the user-selected von Karman power spectrum as the standard spectrum, the simulation was performed using the harmonic superposition method to calculate the fluctuating wind speed at typical points on the wind turbine. ; Average wind speed at typical points on the wind turbine and pulsating wind speed Determine the instantaneous wind speed at typical points on the wind turbine. ; The wind direction angle θ(t) of the wind turbine position during the period from entering to leaving the maximum influence radius of the typhoon is calculated based on the design typhoon wind field parameters. Instantaneous wind speed at typical points on the wind turbine By performing vertical decomposition based on the wind direction angle θ(t), the instantaneous wind speed in the orthogonal direction at each typical point is obtained. and ; Instantaneous wind speeds in orthogonal directions at typical points on the wind turbine and Calculate the instantaneous typhoon load at each typical point on the wind turbine; A finite element model of the wind turbine structure was established, and dynamic time history calculations were carried out by combining the instantaneous typhoon loads of the corresponding parts of each typical point on the wind turbine. Based on the user-selected design typhoon wind field parameters, the average wind speed at the wind turbine's reference point is calculated during the period from when the wind turbine enters to when it leaves the typhoon's maximum influence radius. ,include: Obtain the wind field parameters of the user-selected design typhoon, including the maximum average wind speed of the design typhoon. and the corresponding maximum wind speed radius R1, the design typhoon's maximum influence radius R2 and the corresponding cut-out wind speed And the direction and speed of the typhoon's center path, V. w ; Based on the typhoon's maximum influence radius R2, the vertical distance L from the wind turbine to the typhoon's center, and the typhoon's center's moving speed V. w The distance S(t) between the wind turbine and the center of the typhoon was calculated during the period from when the wind turbine entered the typhoon's maximum influence radius to when it left the typhoon's maximum influence radius. Based on the distance S(t) between the wind turbine and the typhoon center, the maximum wind speed radius R1 of the designed typhoon, and the maximum average wind speed of the wind field Determine the average wind speed at the reference point of the wind turbine during the period from when the turbine enters to when it leaves the maximum influence radius of the typhoon. ; The distance S(t) between the wind turbine and the typhoon center is... ； The average wind speed at the wind turbine reference point for: ； z represents the height above sea level.

2. The method for structural dynamic analysis of offshore wind turbines based on the entire typhoon process as described in claim 1, characterized in that, The height of typical points on the wind turbine and the average wind speed at the wind turbine reference point are used as the basis. Determine the average wind speed at typical points on the fan. ,include: In the formula, This is the wind speed profile index.

3. The method for structural dynamic analysis of offshore wind turbines based on the entire typhoon process as described in claim 1, characterized in that, The method uses the user-selected von Karman power spectrum as the standard spectrum to calculate the fluctuating wind speed at typical points on the wind turbine. ,include: (1) Select the von Karman power spectrum as the standard spectrum, and calculate the power spectral density function according to the following formula: In the formula, The angular frequency of the pulsating wind. The average wind speed during the first 10-minute period after the wind turbine enters the typhoon's influence area; For friction speed, The longitudinal turbulence integral scale; In the formula, z0 is the surface roughness length; (2) Based on the power spectrum of each point Calculate the cross power spectrum using the following formula. : In the formula, It is a coherent function and can be obtained from the relative position coordinates of each point; (3) Discretize the power spectrum in the frequency domain range 0~ If the interior is divided into N parts, then the frequency increment is: in, The upper limit cutoff frequency; N discretized power spectrum matrices are obtained ; (4) Perform Chlesky decomposition on the power spectrum matrix to obtain the lower triangular matrix. ; (5) Based on the lower triangular matrix Harmonic superposition was performed to obtain the 10-minute fluctuating wind speed time history at each typical point: In the formula: for phase angle, A random number between 0 and 2π; (6) Repeat steps (1) to (5) K times to obtain the wind field pulsation wind speed time history during the entire process of the typhoon; In the formula, T is the time required for the wind turbine to go from entering to leaving the typhoon's maximum influence radius; T0 is 10 minutes.

4. The method for structural dynamic analysis of offshore wind turbines based on the entire typhoon process as described in claim 1, characterized in that, The average wind speed based on typical points on the wind turbine and pulsating wind speed Determine the instantaneous wind speed at typical points on the wind turbine. ,include: 。 5. The method for structural dynamic analysis of offshore wind turbines based on the entire typhoon process as described in claim 1, characterized in that, The wind direction angle θ(t) of the wind turbine's position during the period from entering to leaving the maximum influence radius of the typhoon, calculated based on the design typhoon wind field parameters, includes: 。 6. A structural dynamic analysis device for offshore wind turbines based on the entire typhoon process, characterized in that, include: The location wind speed calculation module is used to calculate the average wind speed at the reference point of the wind turbine during the period from when the turbine enters to when it leaves the maximum influence radius of the typhoon, based on the user-selected design typhoon wind field parameters. ; The average wind speed calculation module is used to calculate the average wind speed based on the height of typical points on the wind turbine and the average wind speed at the wind turbine reference point. Determine the average wind speed at typical points on the fan. ; The fluctuating wind speed calculation module is used to simulate and calculate the fluctuating wind speed at typical points on the wind turbine using a user-selected von Karman power spectrum as the standard spectrum and the harmonic superposition method. ; The instantaneous wind speed determination module is used to determine the average wind speed at typical points on the wind turbine. and pulsating wind speed Determine the instantaneous wind speed at typical points on the wind turbine. ; The wind direction angle calculation module is used to calculate the wind direction angle θ(t) of the wind turbine position during the period from entering to leaving the maximum influence radius of the typhoon, based on the design typhoon wind field parameters. The wind speed decomposition module is used to decompose the instantaneous wind speed at typical points on the wind turbine. By performing vertical decomposition based on the wind direction angle θ(t), the instantaneous wind speed in the orthogonal direction at each typical point is obtained. and ; The load calculation module is used to calculate the instantaneous wind speed in the orthogonal direction at typical points on the wind turbine. and Calculate the instantaneous typhoon load at each typical point on the wind turbine; The model calculation module is used to establish a finite element model of the wind turbine structure and perform dynamic time history calculations by combining the typhoon instantaneous loads of the corresponding parts of each typical point on the wind turbine. Based on the user-selected design typhoon wind field parameters, the average wind speed at the wind turbine's reference point is calculated during the period from when the wind turbine enters to when it leaves the typhoon's maximum influence radius. ,include: Obtain the wind field parameters of the user-selected design typhoon, including the maximum average wind speed of the design typhoon. and the corresponding maximum wind speed radius R1, the design typhoon's maximum influence radius R2 and the corresponding cut-out wind speed And the direction and speed of the typhoon's center path, V. w ; Based on the typhoon's maximum influence radius R2, the vertical distance L from the wind turbine to the typhoon's center, and the typhoon's center's moving speed V. w The distance S(t) between the wind turbine and the center of the typhoon was calculated during the period from when the wind turbine entered the typhoon's maximum influence radius to when it left the typhoon's maximum influence radius. Based on the distance S(t) between the wind turbine and the typhoon center, the maximum wind speed radius R1 of the designed typhoon, and the maximum average wind speed of the wind field Determine the average wind speed at the reference point of the wind turbine during the period from when the turbine enters to when it leaves the maximum influence radius of the typhoon. ; The distance S(t) between the wind turbine and the typhoon center is... ； The average wind speed at the wind turbine reference point for: ； z represents the height above sea level.

7. A storage medium storing a computer program executable by a processor, characterized in that: When the computer program is executed, it implements the steps of the offshore wind turbine structure dynamic analysis method based on the full-process typhoon working conditions as described in any one of claims 1 to 5.

8. A structural dynamic analysis device for offshore wind turbines, comprising a memory and a processor, wherein the memory stores a computer program executable by the processor, characterized in that: When the computer program is executed, it implements the steps of the offshore wind turbine structure dynamic analysis method based on the full-process typhoon working conditions as described in any one of claims 1 to 5.

Images