Effects of wind field characteristics on pitch bearing reliability: a case study of 5 MW reference wind turbine at onshore and offshore sites

This paper presents a study on pitch bearing basic rating life affected by wind field characteristics at both onshore and offshore wind sites. The National Renewable Energy Laboratory 5 MW reference wind turbine is selected for the study. Wind field characteristics including reference hub height mean wind speed, wind speed distribution, wind shear, and vertical inflow are studied. A decoupled approach is employed where global analysis is performed first. Second, the load effects from the global analysis are applied on a reference pitch bearing designed based on best industrial practices. For the case study onshore site, it is found that the Kernel density estimation best fits the wind distribution, while the International Electrotechnical Commission proposed distribution appears to be not suitable. Moreover, it is shown that the seed number has high effect on the bearing life in turbulence wind and the wind speeds around rated have the highest contribution in both bearing fatigue damage and extreme load failure. The results contribute to better understanding of the wind field characteristics on the pitch bearing life.


Introduction
Wind turbines as a whole, and multimegawatt offshore wind turbines specifically, are considered as a solution for the large-scale realization of renewable power generation. Practically, a higher failure rate along with repair and maintenance requirements leads to a higher cost of energy [1]. Experiences show that faults or damages in the drivetrain, bearing, and gearbox contribute significantly to the wind turbine's downtime and nonavailability [2,3]. The maintenance and repair are particularly too costly, especially for offshore wind turbines where a special vessel and crane barge are needed to replace a faulty drivetrain. In addition, the harsh offshore weather reduces the repair period to only a few months, often in summer [4]. Pitch systems are among the largest contributor to overall downtime, representing more than 15% of wind turbine failures [5]. In another survey, 25% of total failures occurred in the pitch system in offshore wind turbines, and its failure rate in onshore wind is 10% [6].
From the wind field perspective, different characteristics are recorded during wind site measurement. The wind speed distribution is significant for the wind turbine design because it determines the frequency of occurrence of the individual load components. Different distributions are used in the wind industry. Carta et al. [7] reviewed different probability distribution functions (PDFs) with case studies in the Canary Islands. Even though Weibull distribution has advantages with respect to the PDFs analyzed, it cannot represent all the wind regimes encountered in nature such as, for example, those with high percentages of null wind speeds, and bimodal distributions. Many earlier works have focused on the distribution fitting, parameter estimations, and goodness of the fit metrics [8], but limited works are found on the effect of the wind speed on the pitch bearing load and life.
Each of the wind characteristics has its own importance. Wind site measurement sometimes does not include the vertical component, but its effect on the pitch bearing is valuable for the turbine designer. Seed number is a random number and different seeds affect bearing loads differently and due to the spatial coherence, different random seeds produce a Gaussian distribution of turbulence intensity in the longitudinal wind component [9]. Extreme and normal turbulence have a great role in the bearing's safety factor and life. Consequently, a comprehensive study on the pitch bearing in different aspects of the wind field in extreme load and fatigue life was performed. In addition, the contribution of wind speed on the bearing life and loads in the long term and short term is analyzed and is described in the methodology section. This paper is structured into five sections. In the next section, the reference wind turbine (RWT) is stated together with the nominated wind site. Additionally, the pitch bearing design process is presented. The methodology is described afterward, and the results and discussion present the findings with respect to the wind field, safety factors, and life. In the end, the conclusion is stated.

Reference wind turbine
The National Renewable Energy Laboratory's (NREL's) 5 MW reference wind turbine was considered as a case study [10]. The turbine was considered on two support structures: fixed-bottom and land-based. The turbine configuration is given in Table 1.

Wind sites
Two wind farms, onshore and offshore, have been studied. The offshore wind site is Anholt wind farm [11]    The measurements taken at 65.6 and 85.6 m for the offshore site and 30 and 40 m for the onshore site are considered. These measurements are projected to the hub height using the power law [13]. The position of the onshore and offshore sites taken from Google Maps are depicted in Fig. 1.

Pitch bearing design
Because the NREL 5 MW reference wind turbine does not include a pitch bearing design, a design is specified as described in this section. The main restriction for the bearing design is the blade's root dimensions. According to definition of 5 MW RWT [10], the blade root diameter is 3542 mm. It is assumed that the tension bolts in the blade side and hub side are M42 as a first guess. The bolt circle diameter of the blade root, D b,bcd , is derived in Eq. 1 from the geometrical design of the hole pattern, where e is the distance of the bolt axis from the edge of the clamped part which is equal to e = 1.5D h , and D b,out is the blade root's outer diameter [14]. According to ISO 273, D h , the hole diameter, is equal to 45 mm for bolt M42 [15].
The distance between the bolts is assumed 3D h [14], which leads to the determination of the number of bolts. Other assumptions on the bolt design are as follows: The friction of the surface of the bearing to the hub and blade is assumed to be 0.15. Hexagon head screw with shank according to DIN EN ISO 4014 is assumed [16]. Bolt length assumed to be 300 mm. This parameter can be modified after the bearing design is completed.
The bolt dimensions of M36, M39, M42, M45, and M48 with safety factors against exceeding the yield point, the surface pressure, and slipping with minimal preload are evaluated. The bolt design and calculation for safety factors are performed based on VDI-2230 [14]. The results of safety factors are given in Table 2.
As shown in the table, all bolt arrangements have safety factors above one, but according to VDI 2230 [14] the recommended minimum safety factor SG is 1.2. By this criterion, the arrangement with the M42 bolt satisfied the condition and is considered as the arrangement of the blade and pitch bearing connection. The material of the bolts is grade 10.9 according to ISO 898 [17].
The bearing ball diameter and pitch circle diameter sizing are based on the following conditions: It is assumed that around 80 percent of the raceway is filled with balls [18]. This assumption is consistent with DG03 [19]. It is considered that the distance between the raceway and bolt hole is 25% of the ball diameter [18]. The contact angle is assumed to be assumed 45° [18,20]. The raceway groove radius-to-diameter ratio is considered to be 0.52 [19]. Bearing geometrical dimensions adopted from Rothe Erde catalogue [21] in left and the general position of bearing in wind turbine in right [22] The surface hardness is assumed to be 58 hardness on Rockwell scale C (HRC) [19,20]. The Core hardness is assumed to be 250 Brinell hardness (HB) [19,20].
The sizing of the pitch bearing is performed in a way that all parameters are related to the ball diameter. The ball diameters of 50, 60, 75, 90, and 100 mm were studied. The designed pitch bearing main dimensions are stated in Table 3 and

Methodology
Extraction of the wind characteristics is the first stage in the analysis. It is described in the first subsection. The turbine response simulation, a global analysis, in which the loads and motions are extracted, is described afterward. In the end, the calculations of the bearing life and safety factors are presented.

Wind field characteristics
Several characteristics in the wind fields such as mean wind speed, wind shear, turbulence intensity, and wind speed distribution, which are important for the wind turbine, are described in Germanischer Lloyd (GL) guidelines [23] and standards [13]. Mean wind speed is considered as an annual average of 10-min means of the wind speed at hub height.
The vertical inflow wind is caused in different ways. Hills in the wind stream in front of the turbine and stable atmosphere can cause the inflow wind. In addition, other turbines in the wind farm can produce inflow wind. The International Electrotechnical Commission (IEC) 61400-1 recommends 8°for inclination of the mean flow with respect to a horizontal plane [13].
Generally, two mathematical models have been used to model the vertical profile of wind speed. The logarithmic profile (log law) and power law [24]. In this study, the power law that is recommended by IEC [13] is used. The power law is as follows: where˛is the power law exponent and wind speed (U ) at height (z) can be calculated from wind speed at reference height (z r ). The power law exponent according to GL 2010 guideline [23] and standard IEC 61400 [13] are 0.11 and 0.2, respectively, for extreme and normal onshore wind profile models, and it changes to 0.14 for the normal offshore wind profile model [25]. Turbulence intensity (TI) is defined as a fraction of the standard deviation ( u ) to the mean of wind speed, according to Eq. 3.
The wind speed distribution at the site is significant for the wind turbine design because it determines the frequency of occurrence of the individual load components. Different distributions are used in the wind industry. In this study the following distributions are nominated: Lognormal distribution Weibull distribution Rayleigh distribution Nakagami distribution Kernel density estimation The pdf and cumulative density functions (CDFs) of the nominated distribution are described in Appendix 1.
To evaluate the accuracy of each distribution, three different methods of root-mean-squared error (RMSE), coefficient of determination (R 2 ), and coefficient of efficiency (CE) as criteria for goodness-of-fit are considered. The criteria are calculated by Eqs. 4-6 [26], respectively: where x o;i and x p;i are the observed and predicted data of x i and the ones with the ( -) overbar are the mean of the parameter. Goodness levels for R 2 and CE are, respectively, 0.85 and 0.9 [26].

Turbine response simulation
In this section the global model of the 5 MW reference wind turbine in the OpenFAST as well as the nominated design load cases (DLCs) and wind condition are described.
The wind turbine simulation is performed in OpenFAST [27]. OpenFAST is a multiphysics, multifidelity tool for nonlinear aero-hydro-servo-elastic simulation of wind turbines. The wind turbine model is considered from the reference model in GitHub/OpenFAST [28].
Two types of DLCs are considered in this paper. One set of DLCs is considered in the ultimate strength analysis. The other set is considered for both ultimate and fatigue strength analysis. In all DLCs wind shear is considered according to site or standard conditions. The load results are according to the blade coordinate system, which is described in the FAST user manual [29]. The main load on the pitch bearing is the bending moment. According to Su et al. [30] and Technical University of Denmark (DTU) 10 MW RWT reports [31], DLCs 6.1 and 1.3 are the extreme value for the blade root bending moment. Subsequently, the mentioned DLCs are considered for ultimate strength analysis. According to IEC, DLCs 1.2, 2.4, 3.1, 4.1 and 6.4 are used for fatigue analysis. DLC 1.2 is considered in this study with higher resolution than proposed by IEC. DLC 2.4 is for fault condition does not cause an immediate shutdown of the turbine. In this study, it is considered that fault in electricity result in an immediate shutdown of the turbine. DLCs 3.1 and 4.1 are the startup and shutdown conditions of the turbine and do not relate to wind characteristics and they are not considered. DLC 6.4 is in park condition. Because the turbine is not operating (the turbine is in idle mode), the pitch system is not moving, and it is not possible to calculate the life according to method based on oscillation of the bearing. The list of DLCs and their characteristics are presented in Table 4.
Generally, the definition of DLCs and the wind models are according to IEC 61400-1 [13] and DNV 437 [32].
Two different post-analyses are carried out in the paper; short term and long term. The short-term analysis is referred to as simulation with a short-term wind condition i.e. 10min simulation at a given wind speed. In this case, the wind speed distribution is not considered and is assumed that all wind speeds occur with an equal probability. This is indeed not a real case, but when it is compared with the long term, the influence of the wind speed distribution -which is site specific -is better illustrated. In long-term analysis, the site specific wind speed distribution is considered for the loads and life calculation.

Pitch bearing life analysis
The radial, axial, and bending loads at each blade root obtained from the global analysis are applied on individual pitch bearings. To evaluate the design of oscillating pitch bearings, the classic Lundberg-Palmgren theory for estimating dynamic load capacity and fatigue life for continuously rotating rolling bearings [33] was extended and correlated to laboratory fatigue life data [34]. This included defining the critical amplitude of oscillation, the angle below which the raceway is stressed by only one roller and receives two stress repetitions per oscillation cycle because the pitch oscillation is relatively small. The calculation of the equivalent load and number of loading cycles for oscillatory applications was later presented [35]. These formulations were then simplified and applied on wind turbine pitch bearings [36]. A dither angle was defined for pitch motions even smaller than the critical amplitude, for which the stressed area is only partially uncovered and retraced. In this situation, fretting corrosion is often experienced and the developed fatigue life models are not applicable. A similar model was also applied to a four-point contact ball bearing in a wind turbine [37].
Based on these earlier works, in 2009 NREL published a pitch and yaw bearing design guide (DG03), which summarized the design criteria, calculation methods, and applicable standards recommended for use in performance and life analyses of these bearings [19]. However, other methods presented such as works done by Wöll et al. [38] or Breslau et al. [39], DG03 is very commonly used in the wind industry and guidelines for determining rating life of blade bearings [40]. The DG03 uses a comparison of the maximum contact stress, max , in the limit load condition to the maximum allowable stress of 4,200 MPa, to define the static safety factor as [41] where the maximum contact stress, max , is also expressed in megapascals. It recommends that the safety factor in this condition be greater than 1.5, which is the same as the ISO 76 recommendation for applications subjected to shock loads [42]. The DG03 modifies the basic dynamic axial load rating, C a , and defines a basic dynamic axial load rating for oscillating conditions, C a,osc . The formulation depends on a prorated angle of oscillation, Â , in comparison to the critical angle of oscillation, Â crit .
where the life exponent p = 3 for ball bearings, Z is the number of ball and the prorated angle of oscillation, Â is where K is the number of load cases. Parameters t k and N k denote the decimal fraction of time and speed of oscillation, respectively, for the bearing operated under the condition yielding Â k . A sample of pitch motion at 12 m/s in the normal turbulence model (NTM) and low turbulence model (LTM) in the IEC-onshore condition is depicted in Fig. 3. The dynamic equivalent load, P ea , is formulated as where F r , F a , and M denote the applied radial, axial, and moment loads, respectively. D pw denotes the pitch diameter of the bearing. This formulation takes into account the moment loading on the pitch bearing. To include various duty cycle loading, a prorated dynamic equivalent load N P ea is defined as The rating life is the predicted value of life before the first evidence of fatigue occurs on a raceway. The basic rating life, L 10 , is the rating life associated with 90% reliability for bearings manufactured with commonly used high-quality material, of good manufacturing quality, and operating under conventional operating conditions.
where the resulting value of L 10 is in millions of revolutions of the bearing rings relative to each other. The DG03 uses this basic dynamic axial load rating for oscillating conditions to estimate the basic rating life More information on the procedure and definition of the parameter are presented in reference [41] and DG03 [19].

Results and discussion
Results are categorized in four main subsections: wind field characteristics, extreme load, safety factor, and life.

Offshore wind site
The results of wind speed at hub height are projected using the power law. The power law exponent is equal to 0.1173. This value is derived from averaging the calculated power law in every 10-min measurement step. The wind speed with different distribution functions are depicted in Fig. 4. Different distribution parameters defining the behavior of wind speed are presented in Table 5. Summarized statistical errors of different criteria, goodness-of-fit, in different PDFs are presented in Table 6.
The Weibull distribution has the most favorable value in all three criteria, and the Nakagami and the Rayleigh po-    Table 7. The result shows the similarity and difference of the Anholt wind site with IEC class I offshore wind site.

Onshore wind site
The results of wind speed at hub height are projected by means of the power law. The power law exponent is equal to 0.0991. This value is derived from averaging the calculated power law in every 10-min measurement step. The wind speed distribution with different distribution functions are depicted in Fig. 5. Corresponding PDFs parameters at 90 m, are stated in Table 8.
The results of different criteria of goodness-of-fit in different PDFs are presented in Table 9.
The Kernel distribution has the most favorable value in all three criteria, followed by the Nakagami and Weibull distributions. The results of the lognormal and Rayleigh distribution fitting are below the other three; however, the later   one is the proposed distribution according to the IEC 61400 standard [13], guideline GL 2010 [23], and DNV 437 [32].
Other wind field features are presented in Table 10.
Comparing the result of the onshore and offshore wind distributions, it can be seen that the proposed Rayleigh distribution could not be fitted to the measurement in the onshore as precisely as the offshore wind site. One of the reasons for the lack of accuracy in onshore wind sites is the bimodal pattern of wind, in which the wind speed density curve has an extra hump [7].

Extreme loads at the pitch bearing
Nominated DLCs are simulated, and extreme loads in the blade coordinate system [29] are extracted. There are eight load components in each extreme case. Six of them are forces and moments in the blade coordinate system and the latter two are the in-plane vector resultant force and moment. Each of these eight components has negative and positive extremes that make a matrix of 8 16. The presented format of extreme loads is compatible with IEC 61400 [13]. The detailed results of the extreme loads in different sites are presented in Appendix 2. The equivalent load according to Eq. 9 is the indicator for the load on the pitch bearing. The results of P ea in different wind sites are illustrated in Fig. 6.
The results show that the load is dominant in the maximum of F x , M x , M y , F r , and M r , from which F x and F r are the radial load in the x-direction and the in-plane vector resultant force at the blade root plane, respectively; M x and M y are the bending moments in the xand y-directions, respectively; and M r is the in-plane vector resultant moment at the blade root plane. It should be noted that in the maximum of F x and F r , the bending moment is the rea-  Long-term and short-term contribution of extreme loads at pitch bearing with respect to wind speed at different wind sites K son for increase in the equivalent loads. In addition, all of these extremes resulted from DLC 1.3, ETM wind condition. This is consistent with the results of the DTU 10 MW design [31]. The GL guideline [23] does not consider this DLC, so wind turbine design based on this guideline may result in an underestimation of extreme loads in the pitch bearing.
The short-term and long-term contributions of wind speed on the pitch bearing equivalent loads are illustrated in Fig. 7 at different sites.
From both short-term and long-term results, it can be seen that a wind speed of 10 m/s is dominant in both IEC standards for onshore and offshore, however, the peaks will change in Anholt and Khaf to 11 and 12 m/s, respectively. These two wind speeds bracket the wind turbine rated speed of 11.4 m/s. In addition, the peaks in IEC cases are much sharper than those for Khaf. Khaf has a broader curve with more wind speeds contributing to extreme loads. Moreover, the maximum contribution of IEC cases and Anholt are double compared to Khaf's maximum contribution. Furthermore, the results also express that the peak of contribution increases from 24% to around 28% in the IEC onshore case in the long term while Khaf is still a little above 10%. The wider distribution for Khaf could be the result of higher extreme turbulence intensity with respect to other wind sites.

Pitch bearing safety factor
The result of safety factors in different wind sites are presented in Fig. 8. The safety factor of the pitch bearing in Khaf is lowest, which is expected, due to the high equivalent load observed in Fig. 6. Normal turbulence in Khaf is lower than the IEC standard, but higher extreme turbulence results in higher load and lower safety factor. The normal Fig. 8 Pitch bearing safety factors in different wind sites turbulence in Anholt at 15 m/s is lower than the IEC standard, but the safety factors are comparatively the same. It should be noted that in this paper the safety factor is limited to the static load factor, and the safety factor of the fretting corrosion types of failure is not considered [19].

Pitch bearing life
The life of the bearing according to seed number in turbulence intensity, vertical inflow, and wind shear was studied and the results presented. Moreover, contribution of the wind speed at normal and low turbulence wind regimes in the life of the pitch bearing is stated.
The bearing life in the NTM condition at four different wind sites with six different seed numbers in long term and short term is depicted in Fig. 9. Different seeds in each wind site create different wind speed time series with the same turbulence intensity.
The results both in the long term and short term indicated that bearing life is not independent of the seed number but is affected by it. Anholt shows higher pitch bearing life compared to the other three cases. The calm wind regime in Anholt increases the life of the bearing. The IEC onshore case has the lowest overall and lowest short-term bearing life, while seed number 1 of Khaf has the lowest life in the long term. This change of lowest bearing life from short term to long term is the effect of wind speed distribution at Khaf. However, bearing life does not change considerably in the long term and short term, and the effect of seed numbers is dominant.
The effect of turbulence intensity is studied in all four sites in the short term and long term. The results based on NTM and LTM are depicted in Fig. 10. According to Fig. 10, the TI has a large effect on the life of the bearing; however, it does not happen at the Anholt wind site, especially in the long term. One of the reasons could be the lower TI of the Anholt site. Furthermore, the long-term condition decreases the life in LTM and does not considerably affect the life in NTM.
The results in the IEC offshore case are not comparable in LTM and NTM, because in LTM the wind shear power exponent intentionally changed to 0.01. The shortterm bearing life of four wind sites in Fig. 10 have the same parameters in wind pattern in LTM, and the only difference is the power law exponent that changes from 0.2 in the IEC onshore case to 0.01 in the IEC offshore case. The lowest and highest life belong to Anholt and Khaf, which have a power law exponent of around 0.1, and the bearing life in the IEC onshore and offshore cases is in between. In this regard, wind shear does not have an upward or downward effect on the bearing life.
The effect of the vertical component of wind speed is studied in LTM conditions. The bearing life results are il-  Fig. 11. In bars with a legend (W Vinflow), vertical inflow is considered and in bars with a legend (W/O Vinflow), vertical inflow is not considered.
The Khaf wind site did not have vertical component measurements, so it was omitted from this analysis. It is observed that considering the vertical inflow in the total wind speed increases the life of the bearing. Furthermore, longterm bearing life results, which considered wind speed distribution, have lower life compared to short-term bearing life results. Consequently, considering the vertical inflow increases the bearing life while the wind speed distribution decreases it.
In order to investigate the effect of each wind speed in the life of the pitch bearing, the short-term and long-term contributions of wind speed on the bearing life are presented in Fig. 12. The figure shows that the overall pattern in all four wind sites is the same, and the peaks increased by decreasing the turbulence content of the wind. The peaks are in over-rated wind speed in all wind sites in short term and shifted to under-rated wind speed in long term except Khaf. The largest peaks both in the long term and short term belong to Anholt. Life reduction has a narrower peak at the Anholt site compared to the three others. The contribution of high wind speed, practically above 20 m/s, is decreased drastically. This change has less effect on Khaf because of the high content of high wind speed. Fig. 13 shows the short-term and long-term contribution of LTM wind speed on the bearing life. In all cases, the wind model is LTM, therefore, the contribution of bearing life due to wind shear (short-term) and wind speed recurrence (long-term) are presented.  The short-term figure shows that the contributions in all wind conditions are to some extent the same; therefore, wind shear has a small effect on the contribution of life reduction in wind distribution. Life reduction has a narrower peak in the distribution in LTM compared to NTM. The Anholt NTM result has the same narrow peak as a result of a low TI value.
The result shows that wind speed contribution to bearing life reduction is increased by adding the effect of wind speed recurrence. In addition, compared to extreme cases, the contribution has a wider distribution in all four wind sites.

Conclusion
In this paper, the effect of wind field characteristics on the pitch bearing life was evaluated. The 5 MW reference turbine was employed in two onshore and offshore sites where different distributions were fitted to the wind measurements. While Rayleigh distribution is the proposed distribution by IEC, it was found that the Weibull is a better choice for offshore and Kernel performs best for onshore, at least for the case study sites.
It is also found that the bending moment has the highest contribution to the equivalent bearing load. The extreme turbulence model governs the extreme loads of the pitch bearing and it is crucial to consider it in pitch bearing design. The contribution of wind speeds in the extreme load in the short term and long term show that IEC sites have sharper peaks in wind speed contribution, and it is observed that wind sites similar to IEC can be designed on the extreme load conditions with around-rated wind speeds. For example, the offshore Anholt site pitch bearing can be designed based on extreme loads just by analyzing aroundrated wind speeds. On the contrary, the onshore wind site has a wide distribution, and it is not reasonable to design a bearing just by around-rated wind speeds in extreme conditions. The safety factor of the Khaf onshore wind site is lower than the IEC case, and if the bearing is designed and selected according to IEC results, it will be underestimated.
The bearing life contribution to wind speeds shows the importance of around-rated wind speed analysis. It is also realized that the power exponent in wind shear does not have a definite effect on the bearing life. Additionally, in equal total wind speed, the vertical component of the wind speed increases the bearing life in both short-term and longterm analysis. It is noticed that the bearing life is largely affected by the turbulence intensity of the wind speed. This change is more noticeable in the short term than in the long term of the bearing life. The effect of random seed number of turbulence on bearing life is considerable for both the short term and long term, and further investigation is needed besides studying more wind sites with different topography.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4. 0/.

Probability distributions
In this section, the definition and equations of the nominated distributions are presented. The following parts are extracted from references [43] and [44].

Lognormal distribution
The logarithmic normal or simply lognormal distribution is a probability that its PDF follows Eq. 14.
where and are mean and standard deviation of lnx respectively. The cumulative distribution function is stated in Eq. 15 where˚is CDF of standard normal distribution (N(0,1)).

Weibull distribution
The probability density function of the Weibull distribution is described in Eq. 16. The parameters k and C are the shape and scale parameters of the Weibull function.
The CDF of the Weibull distribution is presented in Eq. 17.
K There are different ways to estimate the scale and shape parameters of the distribution function. In this study, the maximum likelihood estimation (MLE) is used. Shape and scale parameters are estimated by following Eqs. 18 and 19, respectively [45].
In this method, numerical iterations are required to determine the shape parameter of the Weibull function.

Rayleigh distribution
Rayleigh distribution is a special case of the Weibull distribution with a shape parameter equal to 2.0. Subsequently, the PDF and CDF follow the Eqs. 20 and 21, respectively.

Nakagami distribution
The Nakagami distribution is commonly used in communication theory. The PDF of the distribution follows the Eq. 22.
The CDF of the Nakagami distribution is as Eq. 24 where P is the regularized (lower) incomplete gamma function.
The shape and speed parameters are calculated, respectively, as [46]: where x 2 is mean of x 2 .

Extreme load results
In this section the extreme loads at the blade root in different wind sites are presented. Figs. 14-17 relate to IEC-onshore, IEC-offshore, Anholt, and Khaf wind sites, respectively.