Failure Analysis of Large-Scale Wind Power Structure under Simulated Typhoon

Recently, a number of wind power structures in tropical cyclone zones are damaged by typhoon. In order to study the failure mechanics and failure modes of wind power structure subjected to typhoon, the typhoon wind field in Dongtai wind farm is simulated based on the classical autoregressive (AR)model and a regional power-spectrum-density (PSD)model, and the simulated spectrum is verified to be in good agreement with the target spectrum. An integrated finite element (FE) model of wind power structure, composed of rotor, nacelle, tower, pile cap, andPHCpiles, is established.Modal analysis reveals that pile stiffness decreases the structure’s natural frequencies, especially for high order frequencies. Structural responses under the simulated typhoon are calculated by dynamic analysis. Results show that tower buckling is the most prone failure mode of the structure. The horizontal displacement of the hub and the axial force of the most unfavorable piles are both under the limit. This study provides a way to the antityphoon design of large-scale wind power structures.


Introduction
With the increase of investment of clean energy from governments, wind energy has grown enormously all over the world in the past decade and will stand to benefit from its role as both a source of energy security and a key to solving the problem of climate change in the future.Compared with the year 2011, wind power market grew by more than 10% in 2012, and the new global total installed wind power capacity at the end of 2012 was 282.5 GW, representing cumulative market growth of more than 19% [1].And China has already taken the leadership position of cumulative wind power capacity in 2011.
It should be noted that more and more wind farms are established in tropical cyclone zones, making the wind turbine and its support structure have to face the threat of typhoon.For instance, typhoon Maemi struck the Miyakojima Island with an average speed of 38.4 m/s and a maximum gust of 74.1 m/s on September 11, 2003.All of the wind turbines on the island were extensively damaged.Three of six turbines collapsed and the others suffered from destructive damage, whose blade were broken or the nacelle cover drooped.Based on FEM simulation and wind response analysis, Ishihara et al. [2] found that the overlarge maximum bending moment was the reason for the buckling and collapse of the tower.In September 2003, typhoon Dujuan attacked the Honghaiwan wind farm of Guangdong, and varying degrees of damage were caused to thirteen wind turbines.In August 2006, typhoon Saomai passed through the Hedingshan wind farm of Zhejiang and led to collapse of five wind turbines [3].
Although massive losses in coastal wind farms have been caused by typhoon, the failure mechanism and failure modes of wind power structures are not clear yet.Generally, modern wind turbines are designed mainly according to the international standard (IEC 61400-1) [4] that is based on European and North American conditions without the experience from areas attacked by typhoon.On the other hand, designers and engineers usually pay more attention to the wind turbine, including rotor, nacelle, and tower, but the soil-structure interaction (SSI) is always out of consideration.The two reasons above result in the ultimate wind load and the structural responses are inevitably underestimated; no need to mention that the characteristic of typhoon is essentially different from the normal wind.Recently, attention has been increasingly paid to the typhoon-induced damage of wind turbines.
Li et al. [3] analyzed the data of typhoon Saomai when it was passing through Hedingshan wind farm.It was pointed out that the violent wind, drastic turbulence, and sudden change of wind direction were major factors of wind turbine failures.From the view of economics, Clausen et al. [5] estimated that the cost of wind turbine in tropical cyclone zones is 20-30% more than that in normal areas.Chou et al. [6] examined the causes of blade damage of five wind turbines in Changhua Coastal Industrial Park of Taiwan.It was found that longterm effects such as blade fatigue could cause local crack damage and delamination, and the actual blade damage occurred mainly on the back edge of the blade near the wing section.
In offshore wind turbines, the structural responses are driven not only by wind but also by water wave.Consequently, the dynamic analysis of offshore wind turbines is more complex than that of land-based ones.Karimirad and Moan [7] employed an advanced blade element momentum theory, panel method, and the Morison formula to study the dynamic structural and motion responses of a spar-type floating offshore wind turbine caused by stochastic nonlinear wave and wind loads in harsh and operational environmental conditions.It is found that the wind turbulence affects power production rather than the dynamic motion and structural responses.Based on the Arbitrary Lagrangian-Eulerian (ALE) method and FEM, Zhang et al. [8] simulated a wind-induced aqueduct-water coupling system installed with isolated bearings.It is found that the resonance can be eliminated under high-order-frequency excitation and may happen when the excitation frequency is close to the firstorder water sloshing frequency.Karimirad and Moan [9] presented a stochastic dynamic response analysis of a tension leg spar-type wind turbine subjected to wind and wave actions.
To study the power performance and structural integrity of the system, negative damping, rotor configuration, and tower shadow effects are discussed.
This paper aims at analyzing the potential failure modes of wind power structure subjected to typhoon and is organized as follows.The location and wind condition of Dongtai wind farm is introduced in Section 2. In Section 3, the fluctuating wind speed of typhoon is simulated and verified based on a regional PSD mode and the AR method, followed by a presentation of wind load calculation in the same part.A FE model of a 1.5 MW wind turbine with pile foundation is established in Section 4, and structure's natural frequencies and windinduced responses are calculated by elastic dynamic analysis.The failure modes of the wind power structure are concluded finally.

Farm Location and Wind Condition
Dongtai wind farm is located in Dongtai City, Jiangsu Province, southeastern China, at 120 ∘ 54  east latitude and 32 ∘ 47  north longitude, as shown in Figure 1.Totally, one hundred and thirty-four 1.5 MW wind turbines are installed in a 60 km 2 site.Wind resource in this area is rich in east and poor in west, forming a long and narrow wind speed surging region.
Although Dongtai is not a typhoon landing site in general, it still has to face the threat of typhoon.For example, in 2009, 41.5-50.914-15 Super typhoon (SuperTY) ≥51.0 ≥16 typhoon Morakot passed through Dongtai with an average wind speed of 33 m/s and moved into East China Sea, as shown in Figure 2. In 2012, typhoon Damrey landed in Xiangshui of Jiangsu Province with an average wind speed of 35 m/s and brushed past Dongtai, as shown in Figure 3.It should be noted that, according to the extreme wind speed model (EWM) [4], the extreme wind speed with a recurrence period of 50 years of the wind farm is 24.37 m/s, which is far below the average wind speed of the two typhoons mentioned above.
According to [10], tropical cyclones are classified as 6 grades, as shown in Table 1.It can be seen that the lower limit of typhoon's average wind speed is 32.7 m/s, which is employed to simulate the fluctuating wind speed in this paper.latter represents turbulence and randomness.Consequently, the instantaneous velocity of wind can be described as

Typhoon Wind Load Simulation
where   and   () are the mean wind velocity and the fluctuating wind velocity at height  and time , respectively.  obeys the logarithmic law and can be described as [11]: where  ref is a reference wind velocity measured at the reference height  ref , and  0 is the roughness length.
According to the random vibration theory, the fluctuating wind velocity can be regarded as a zero-mean Gaussian process and described by the power-spectrum-density (PSD) model in the frequency domain [12].A lot of empirical spectrum models, mostly for the longitudinal wind velocity, have been published, such as von Karman spectrum [13], Davenport spectrum [14], and Simiu spectrum [15].Considering the features of typhoon, such as strong turbulence, high wind speed, and significant regional characteristics, the PSD model of typhoon is essentially different from that of normal wind.Shi et al. [16] proposed an empirical spectrum that does not change with altitude, according to the observed data in Shanghai since 1956 and taking the characteristics of the typhoon turbulence scale varying with altitude into account.Considering that the distance between Dongtai and Shanghai is acceptable, Shi's spectrum is employed here.The model is defined as where  V () is the PSD of the fluctuating wind speed,  is the ground roughness coefficient, V 10 is the average wind speed at the height of 10 m,  is the wind frequency, and  is the turbulence integral scale factor.

Fluctuating Wind Speed
Simulation.Discrete models of wind speed based on Box-Jenkins methods are used commonly in time-series analysis.These models, including autoregressive (AR) [17], moving average (MA) [18], autoregressive moving average (ARMA) [19], and autoregressive integrated moving average (ARIMA) models [20], are usually termed as linear filtering method and are able to reproduce the statistical properties of the fluctuating wind speed for a particular site with an acceptable computational cost.The classical AR model is employed to simulate the wind field of typhoon in this paper [21].To generate a family of  processes, the following equations are used: where   () and   ( − Δ),  = 1, . . ., , are the fluctuating wind speed at time  and time (−Δ), respectively, and   () is made of normally distributed random numbers with zero mean and unit variance.  are  ×  matrices,  = 1, . . ., , and  is regarded as the order of the model.The covariance between   () and   ( − Δ) can be denoted as [22]: Considering that   () and   ( − Δ) are stationary stochastic processes with zero mean, (10) can be rewritten as Postmultiplying ( 6) by u( − Δ)  and applying mathematical expectation derive The cospectrum and the covariance function satisfy the Wiener-Khintchine equation; that is, where    () denotes the cross-spectral density function of point (  ,   ,   ) and point (  ,   ,   ) and can be expressed as where    () represents the PSD of point  associated with the process   () and has the form of (3), and coh  () represents the coherence function of longitudinal fluctuations between point  and point .
Considering that wind power structure is a typical slender structure, a simplified expression of coh  () proposed by Shiotani and Avai [23] is adopted here: and   is proposed to be 60 [24].
Postmultiplying (6) by u() = [ 1 (), . . .,   ()] derives The vector N() of random series, having R  as their cross-correlation matrix, can be obtained by the linear combination where n() is a set of  independent random processes with zero mean and unit variance and L is a lower triangular matrix and can be calculated by a Cholesky factorization of R  and Substituting  (obtained from (10)) and N() (obtained from ( 17

Verification of the Simulated Wind Speed.
In order to verify the simulated wind speed history, the PSD of the simulated wind is compared with that of the target spectrum (see (3)) in Figure 5.It can be seen that the simulated spectrum agrees well with the target spectrum in the frequency domain, which implies that the simulated wind has a similar power distribution as the natural wind and can be used to calculate the structural typhoon-induced response.

Wind Load
Calculation.The wind rotor is a complex aerodynamic system that converts wind energy into mechanical power.The blade element theory is commonly used to calculate the wind load acting on the rotating rotor.Considering that the wind turbine should be shut down before the coming typhoon, the aerodynamic force acting on the rotor can be estimated by where   is the equivalent drag efficient of the rotor and 1.1 is used here [25],  is the air density, u  () is the total wind speed at the hub height, and  is the swept area.The relationship between wind speed and wind pressure under ambient conditions is described by the Bernoulli equation: where w() is the wind pressure,  is the unit weight of air, and  is the gravitational acceleration.
Assuming that the wind speed is parallel to the normal direction of the rotor, which is the most unfavorable condition for the support structure, the wind load acting on each element of tower can be calculated by the following expression: where   is the element area and   is the shape coefficient of the element.

Differential Equation of Motion.
In order to study the dynamic responses of the wind power structure subjected to typhoon load varying in time, a transient dynamic analysis should be conducted.Based on the D' Alemberts principle and due to the discretization process of a continuous structure with FEs, the following equation of motion can be derived: where M, C, and K denote the structural mass, damping, and stiffness matrices, respectively, ü  (), u  (), and u  () are the vectors of nodal accelerations, velocities, and displacements, respectively, and f() is the vector of applied forces.The Newmark integration method, which is an implicit time integration algorithm, is employed to solve (23).It should be pointed out that damping plays an important role in the dynamic response of the wind power structure.Ambient vibration tests are always employed to estimate the damping ratio because of its strong advantage of being practical and economical, using the freely available ambient wind wave excitation.Shirzadeh et al. [26] measured the damping ratio of a 3 MW offshore wind turbine under ambient vibration, and the measured values are 1.05% and 1.27% for, respectively, the first FA mode and the first SS mode.Similarly, Ma et al. [27] obtained the damping ratio of a 1.5 MW land-based wind turbine tower.Considering that the wind turbine and support structure of the latter case is the same as this study, the measured damping ratio 1.75% is used here.

Modal Analysis Algorithm.
A modal analysis is performed to calculate the natural frequencies of the wind power structure.Omitting the damping matrix and force vector in (23), the free vibration equation of the wind power structure for eigenvalue analysis is obtained Assume the general form of the solution is Equation ( 24) becomes The block Lanczos method is a very efficient and robust algorithm to perform a modal analysis for large models; thus it is employed on the platform of ANSYS [28].

FE Model.
As mentioned before, 1.5 MW wind turbines are installed in Dongtai wind farm.The tower is 62.75 m in height (the hub height is 65 m) and is fixed on the pile cap by an anchor ring.The tower consists of 3 segments of cylinder with dimensions shown in Table 2.The material parameters of the nacelle, rotor, tower, and pile cap are listed in Table 3.
There are thirty PHC piles arranged under the pile cap.Six of them are along the inner circle with a diameter of 4.1 m, and the others are along the outer circle with a diameter of 16.8 m, as shown in Figure 6.Based on pile tests, the horizontal stiffness and vertical stiffness of the pile are  ℎ = 1.8 × 10 7 N⋅m −2 and  V = 1.9 × 10 8 N⋅m −2 , respectively.
The FE model of the wind power structure is shown in Figure 7. Totally, 14898 elements are used with the types listed in Table 4. Considering that the failure mechanism of blade is beyond the scope of this study, solid element is employed to simulate the rotor with an equivalent mass.
Table 5 compares the natural frequencies of the structure with piles (case 1) and the structure fixed on a rigid foundation (case 2).It can be seen that the primary and secondary frequencies are around 0.41 Hz in both cases, which is close to the predominant frequency of typhoon.Consequently, a resonance is easy to happen.Furthermore, the natural frequencies of wind power structure in case 2 are higher than that in case 1, especially for 6th to 10th orders frequency.It reveals that piles decrease the structural stiffness and should be taken into account in the structural dynamic analysis.

Failure Mode Analysis
5.1.Horizontal Displacement of the Hub.As a typical slender structure, wind power structure is sensitive to horizontal displacement.Consequently, the horizontal displacement of the hub ( ℎ ) is usually regarded as a safety control index of the wind power structure.In Figure 8, it can be seen that the hub moves from −0.160 m to 1.125 m, and the maximum displacement is under the limit value of 1.3 m, which is 1/50 of the hub height [29].

Stress of the Tower.
Under the wind load, the vertical stress of the tower, which is predominant, is small at the top and large at the bottom.Accordingly, the shell thickness of the tower gradually increases from the top to the bottom.However, the tower bottom is still prone to buckling.Figures 9(a     the tower, respectively.Obviously, at quite a lot of time steps, the stress of the lowest elements, not only at the windward side but also at the leeward side, is over the steel's yield strength.As a result, buckling originates from the tower bottom.Considering that the maximum stress appears at time = 111.25 sec, the corresponding equivalent stress distributions are described in Figure 10.Due to the bending stress and the stress concentration, the stress reaches the maximum at the junction of the lower segment and the anchor ring.The yield range looks spindle-like with the height of 0.5 m approximately.

Stress of the Anchor Ring.
Owing to the large thickness of the anchor ring, the equivalent stress decreases rapidly and is under the steel's yield strength at most time steps.Assuming the stress concentration is eliminated by engineering means, the yield possibility of the anchor ring can be reduced effectively (Figure 11).The maximum equivalent stress of the anchor ring also appears at time = 111.25 sec and the stress distribution is described in Figure 12.It is shown that the ring protrudes at the windward side and indents at the leeward side, and it can be regarded as a typical buckling mode of the tower [2].

Axial Force of the Unfavorable Piles.
Each foundation pile bears different bending moment based on moment distribution method.Consequently, the piles at the windward side and the leeward side on the outer circle are the most unfavorable, and their axial force histories are shown in Figure 13.Pile tests that provided the ultimate compress and tensile bearing capacity of single pile are 1250 kN and 470 kN, respectively, and the ultimate tensile bearing capacity of the pile body is 2500 kN.Thus, the pile foundation has enough deign safety margin in this case.

Conclusion
Based on the classical AR model and a regional PSD model, the wind field of typhoon in a coastal wind farm is simulated and verified.A FE model of the wind power structure, including nacelle, rotor, tower, pile cap, and PHC piles, is established.Modal analysis shows that the natural frequencies of the wind power structure fixed on rigid foundation are higher than that of the structure with piles, which implies that piles decrease the structural stiffness and should be taken into account in the structural dynamic analysis.The primary and secondary frequencies of the structure are close to the predominant frequency of typhoon, leading to a resonance between the structure and the wind.Dynamic analysis reveals
)) into(6), u(), including  series of fluctuating wind speed, is obtained.Based on the AR model, the fluctuating wind speeds at the height of 5 m, 15 m, 25 m, 35 m, 45 m, 55 m, and 65 m in the wind farm are simulated.For simplicity, the wind speed histories at the height of 5 m, 35 m, and 65 m are shown in Figures 4(a)-4(c), respectively.The simulation parameters are defined as follows:  = 4, point number = 7, Δ = 0.1 sec, and Δ = 0.01 Hz.
) and 9(b) show the Von-Mises equivalent stress histories of the lowest element in the windward side and the leeward side of

Figure 8 :
Figure 8: The horizontal displacement history of the hub.
The element at the leeward

Figure 9 :
Figure 9: Von-Mises equivalent stress history at the bottom of tower.

Figure 10 :
Figure 10: Von-Mises equivalent stress at the bottom of tower while time = 111.25 sec.

Figure 11 :
Figure 11: Von-Mises equivalent stress histories at the bottom of anchor ring.

Figure 12 :
Figure 12: Von-Mises equivalent stress distribution of the anchor ring.

Table 2 :
Dimensions of the tower.

Table 4 :
Element types of structural components.

Table 5 :
Natural frequencies of wind power structures.