Comparison between the BEM and CFD methods.
With the negative impact of conventional energy resources that have been used worldwide, there is a demand for using other resources such as wind energy. Tons of researches have been applied around the globe on the process of designing and manufacturing wind energy conversion systems. In the present chapter, we are concentrating on wind turbine blades’ structural design process. The structural design of a wind turbine blade includes defining the wind turbine loads, selecting a suitable material, creating a structural model, and solving the model using the finite element method. This process will be repeated several times until a final design is achieved. The present chapter includes a discussion on the finite element method and wind turbine aeroelasticity.
- wind turbine blades
- composite structures
- finite element method
With the increase in clean energy demand, wind turbine power and subsequently its size increased dramatically. The power generated from wind turbine increased from approximately 30 GW in 2001 to about 600 GW around the globe as depicted in Figure 1.
While the total power capacity in Europe increased from about 70 GW to 180 GW in the period from 2008 to 2018, as shown in Figure 2, in the United States, the wind power generated from wind turbines increased from 0 to 100 GW in the period from 2000 to 2019 .
This increase in wind turbine size makes it important to efficiently design wind turbine structure. A blade structure must be stiff enough, so it does not fail due to wind turbine loads, and at the same moment it should be designed in lightweight and low cost. In the present chapter, we discuss the general structural design process for a wind turbine blade for maximum strength to weight ratio. This process starts by calculating the wind turbine loads, selecting a suitable material, creating a structural model, and applying the finite element method for predicting aeroelastic stability boundaries.
2. Blade aerodynamic loads
A wind turbine structure must be stiff enough to stand for the flow around. For this reason, the first step in the structure design process is to determine the aerodynamic loads applied to a wind turbine blade. There are mainly two methods for predicting a wind turbine aerodynamic load: blade element momentum (BEM) and computational fluid dynamics (CFD). The BEM can be classified as an analytical method that is fast and accurate. On the other hand, the CFD analysis is a numerical method which is based on numerical and empirical approximations. In general, the BEM is preferred for initial design and load estimation, while the CFD simulation is preferred for detailed design. Designers usually start using BEM in preliminary design and end up using the CFD simulation for detailed analysis and design. A comparison between the BEM and CFD analyses is listed in Table 1.
|Analytical method||Numerical method|
|Simple to be derived and used||Relatively complex in its formulation and application|
|Has relatively low computational time||Has relatively high computational time|
|In general, it does not consider 3D effects (just some corrections such as tip losses)||Depicts 3D effects in detail|
|Does not consider turbulence effects||Considers turbulence effects|
|Recommended for preliminary design||Recommended for detailed simulation|
Despite the simplicity of the BEM, it can obtain accurate results as the CFD analysis . However, the BEM is not preferred for detailed simulation. There are software available for both the BEM and CFD analyses, such as wind turbine design and simulation (Q-blade) that is based on the BEM  and ANSYS Fluent that is based on CFD simulation . For more details about ANSYS CFD simulation, we refer the reader to ANSYS Fluent Theory Guide . In the following section, we will summarize the BEM method.
2.1 Blade element momentum theory
The BEM is considered as a simple and fast method in calculating the applied aerodynamic loads on wind turbine blade. The blade is divided into
In Figure 3, denotes the upstream velocity, is the flow density,
In BEM, we calculate the forces on each blade element, and then the total forces over the blade can be calculated using summation. The lift force and blade angles are defined in Figure 4, which Figure 4 shows the triangle of in-plane velocity (no subscript), triangle of upstream velocity (subscript 1), and triangle of downstream velocity (subscript 3) , where, is the relative velocity angle.
To simplify the process of calculating blade loads using BEM, we describe it as follows:
where is the tip speed ratio.
Set the relative flow speed upper and lower limits:
Calculate the angle of attack from ; then, one can obtain the section lift and drag coefficient from the airfoil characteristics.
Calculate the parameter .
If , use the Glauert correction , or
if , use the Prandtl correction .
Solve the BEM nonlinear equation to obtain :
Repeat steps from 1 to 5 till convergence.
After convergence, calculate the relative velocity w:
Calculate the aerodynamic forces (lift and drag forces):
The thrust force can be calculated from the equation:
The blade tangential (circumferential force)
A schematic for a blade aerodynamic pressure is shown in Figure 5. It is seen that the maximum force occurs at the blade root and the minimum at the tip. For this reason, wind turbine blades are usually designed with taper in which the airfoil thickness increases toward the blade root. This property makes the blade structure stiffer at the root and lighter at the tip. For this reason, it is recommended to also design the wind turbine structure with taper. In other words, ribs, spars, and skins have thickness at the blade root higher than the thickness at the tip.
3. Material selection
Wind turbines have been made from different materials such as wood, aluminum, and composites. Modern wind turbines are usually made from composites such as carbon fibers and fiber glass. The wide use of composite materials turns to its relative high stiffness to weight ratio in addition to its ability to form complex shapes. Composites are found to be efficient with the large increase in wind turbines’ size and capacities. They also can be tailored to satisfy different stiffness and weight requirements. Composite fiber can be used in different orientation to improve the blade directional stiffness in addition to bending and torsional rigidity. Figure 6 shows the layup process for manufacturing a 6-m-diameter wind turbine that was designed and manufactured at Cairo University laboratories.
Two important things have to be considered when selecting composites: first, selecting the proper fiber direction for a blade structure, and second, insuring that the final product (blade structure) has the same material properties as it was desired in the design process. The former point can be overcome by applying a proper optimization process to select the best laminate configuration for maximizing the blade stiffness to weight ratio, while the latter issue can be resolved by testing the layup configuration after manufacturing some samples to make sure they have the same desired properties plus making experimental tests and measurements to compare between what was designed and what was manufactured.
There are an infinite number of orientations for composite structures. A composite layer can be oriented in any direction. However, it is important when selecting a blade laminate configuration to consider 0, 90, and ±45° angles for the blade skin. These directions are the most important fiber directions to increase the longitudinal, directional, torsional, and bending stiffnesses of a blade (Figure 7).
4. Blade construction and finite element analysis
The objective of wind turbine structure is to transfer and stand for wind turbine loads. Thus, it should be stiff enough to satisfy this objective. The structure weight is also important to be minimum as possible. A typical wind turbine structure consists of the skins, ribs, spar, and root or hub that connects between the blade and the wind turbine tower, as shown in Figure 8.
The ribs represent the aerodynamic profile shape for a blade. They distribute the aerodynamic loads and transform them to the main spar, in addition to maintaining the skin profile shape. The skins protect and cover the blade structure elements. Stringers stiffen the skins and connect between the structure of the skins and ribs. A spar represents the main structure element which carries the blade’s main loads and transforms them to its root which in turn connects between the blade and the hub. A spar is usually consist of upper and lower flanges (caps) in addition to a shear web (Figure 9). The shear web performs high resistance to shear force in which the bending moment over the blade is transformed into in-plane shear forces that are carried by the shear web.
The geometric model available in Figure 8 represents an equivalent model to the real wind turbine. This model usually is not suitable for finite element (FE) analysis. However, in finite element analysis, we use an equivalent model with some approximations that do not affect the accuracy of the analysis, but these approximations increase the speed and efficiency of the FE process. A geometric model that is physically a 3D model can be approximated into 2D or even 1D simple models. The blade skin and ribs, for instance, have a thickness dimension that is relatively smaller than the other blade dimensions; for this reason the skin and ribs can be approximated into 2D plates in finite element analysis. Geometrically, the blade skin and ribs are modeled as surfaces. Another example is the blade stringers. They have cross-section dimensions that are relatively smaller than the length dimension. Thus, stringers are usually modeled as 1D beams in the finite element analysis. These beams are geometrically represented by lines. These modeling concepts save significant time and effort in structural analysis and design process. An example of a surface approximation for a 6-m-diameter blade ribs and lower skin is shown in Figure 10.
In Figure 10 the blade ribs and skin are represented by surfaces with zero thickness in which we assume there is no change either in stress or strain through their thickness. This approximation is valid as long as the thicknesses of the ribs and skin are relatively smaller than the other surface dimensions. Based on this approximation, a blade is meshed as a 3D model, but it is modeled numerically in FE using quadrilateral element. It is found that higher-order elements such as nine-node element can obtain results more accurate than linear quadrilateral elements in elastic and aeroelastic analyses . So, the element selection is also important. After approximating the blade geometry, selecting the proper element, and defining the finite element model, it is important to perform a convergence test to select the best element size and density for the blade model . Figure 11 shows how the approximation from the 3D structure model to the 2D reference element is performed and the different reference coordinates that are used.
The finite element process starts by approximating the physical model as we discussed, then defining the material properties, and selecting a suitable element. After that start the mesh process which transforms the continuous geometric model into discrete elements and nodes. Finally, the model boundary conditions are defined, and the model is solved. The objective of the finite element analysis is to solve the general equilibrium equation 
in which is the Cauchy stress tensor,
4.1 Static structural analysis
We called this analysis one-way aeroelastic analysis in which the steady aerodynamic loads are defined and transformed to the structural model and then stresses, strains, and deformations are calculated. It is important in this analysis that the blade maximum deflection satisfy the minimum tower-to-tip clearance . The maximum stresses and strains must be lower than the model allowable stresses and strains. The allowable stresses and strains are defined based on the material yield stresses or strains divided by a factor of safety.
4.2 Modal analysis
In modal analysis, we aim to calculate blade natural modes. Natural modes define the blade natural frequencies and their corresponding mode shapes at which the blade can fail if it excited at any of these modes. Therefore, a blade should be designed with natural frequencies higher than any surrounding frequency.
The modal analysis can be achieved by solving an eigenvalue problem based on the general equilibrium equation that has the form
4.3 Dynamic analysis
In dynamic analysis, the unsteady aerodynamic loads are transformed into the wind turbine or blade structural model, and then the model is solved considering the time variation of both the loads and structural response. In dynamic analysis, the structure must be stable and safe from failure during its lifetime.
4.4 Aeroelastic analysis
A two-way aeroelastic analysis includes the calculation of blade aerodynamic loads and then transforms these loads into blade structure. After that the structure deformation is calculated and transformed back to the aerodynamic model to calculate new load distribution, and the process goes back and forth until a convergence point is determined. In the following section, we will discuss the wind turbine aeroelastic analysis process.
One of the concerns in aeroelastic analysis is the connection between the aerodynamic model and the structural model. The aerodynamic model is usually made using CFD analysis, while the structural model is usually made using FEM. Different methods have been used to connect between the aerodynamic and structural models such as spline methods that are used in ANSYS software. It is found that the use of the finite element shape functions is more efficient in connecting between the aerodynamic and structural models, than the spline methods  (Figure 13).
5. Aeroelastic analysis of wind turbine blades
With the large increase in wind turbine size, the determination of aeroelastic instability of wind turbine blades becomes a must. Nowadays, wind turbines are designed with high speed, large scale, and light structure. All these specifications make them vulnerable to aeroelastic instabilities. There are two important things in aeroelastic analysis: (1) calculating the exact loads of a wind turbine blade and (2) determining the static (divergence speed) and dynamic (flutter speed) stability limits for a wind turbine.
5.1 Aerodynamic loads
In relatively large-scale wind turbines, the structure elastic nature allows it to deform significantly due to the applied aerodynamic loads. This deformation in turn affects the wind turbine thrust and output power. An example is the wind turbine blade shown in Figure 14. It is designed for a 1550-kN thrust; however, the wind turbine actual thrust is 1500 due to the aeroelastic effect.
As the wind turbine size becomes larger and larger, these effects become more significant. Therefore, for large-size wind turbines, the aeroelastic effects become more important and can affect significantly the actual wind turbine power. In other words, the wind turbine aeroelastic deformation can significantly decrease the wind turbine power and at the same time decrease the tip-to-tower distance which can lead to wind turbine failure.
5.2 Static and dynamic aeroelastic instability
The first time for an aeroelastic instability accident to appear was in Washington, USA, in 1940, as a result of Tacoma Narrows Bridge collapse. Then, the aeroelastic instability analyses become important for aerospace structure such as aircraft wings  and wind turbines. Recently, researchers start to investigate different aeroelastic phenomena in wind turbines. Aeroelastic instability analysis includes two aeroelastic phenomena: divergence and flutter. The divergence is known as a static aeroelastic instability due to the blade’s large torsional displacement in response to static aerodynamic loads. It occurs at certain speed known as the divergence speed. On the other hand, flutter is a dynamic aeroelastic instability that occurs due to unstable vibrations that result from a coupling between twisting and bending modes at certain speed called the flutter speed.
No divergence or flutter failure is recorded for small-scale wind turbines (lower than 10 m diameter)  because small wind turbines have relatively stiff structure with low tip speed ratio. In such wind turbines, the flutter speed is found to be five times the wind turbine operating speed . For MW wind turbines, the critical flutter speed becomes lower and lower, and subsequently the aeroelastic analysis becomes more important.
A wind turbine divergence speed is calculated based on steady aerodynamic loads in the interaction with static elastic structure based on an iterative technique. Figure 15 shows a typical process for determining a wind turbine divergence speed based on the finite element method. On the other hand, a wind turbine flutter speed is calculated based on unsteady aerodynamic loads in the interaction with dynamic elastic structure based on an iterative technique. Figure 16 shows a typical process for determining a wind turbine flutter speed based on the finite element method.
Notice that the finite element equation given in Figures 15 and 16, , represents the general finite element equation derived from the structure equilibrium equation (Eq. (11)). The first term in the equation represents the structure inertia effects in which is the element mass matrix. In divergence analysis the inertia term vanishes, while it is an important and mandatory term in flutter analysis. The second term defines the structure damping effect in which is the damping term from both the structure and aerodynamic loads. Then structure damping can only calculate from the experiments. However, if no experiment data available for structure damping, the damping term is defined as a factor multiplied to the structure mass and stiffness matrices. The third term is the structure and aerodynamic stiffnesses. Both the structure and aerodynamic stiffnesses are calculated based on the structure stiffness matrix and the aerodynamic coefficient matrix, respectively. Once, the general finite element equation is defined based on the structure and aerodynamic analyses, it is transformed to an eigenvalue problem to determine the critical speeds for both divergence (Figure 15) and flutter (Figure 16).