Open access peer-reviewed chapter

Structural Optimization of Wind Turbine Blades for Improved Dynamic Performance

By Gerges Edwar Mehanny Beshay and Karam Yousef Maalawi

Submitted: October 30th 2018Reviewed: February 6th 2020Published: April 15th 2020

DOI: 10.5772/intechopen.91643

Downloaded: 45

Abstract

The design of the main structure of a wind turbine blade is optimized aiming at the improvement of the overall dynamic performance. Three optimization strategies are developed and tested. The first fundamental one is based on minimizing the total structural mass of the blade spar under frequency and strength constraints. The second and third strategies are concerned with the reduction of the overall vibration level by either minimizing a frequency-placement index or maximizing the natural frequencies and placing them at their target values to avoid large amplitudes and resonance occurrence. Design variables include cross-sectional dimensions and material properties along the spanwise direction of the blade spar. The optimization problem is formulated as a nonlinear constrained problem solved by sequential quadratic programming (SQP) technique. Two specific layup configurations, namely, circumferentially asymmetric stiffness ( CAS ) and circumferentially uniform stiffness ( CUS ), are analyzed. Exact analytical methods are applied to calculate the natural modes of vibration of a composite, thin-walled, tapered blade spar. The influence of coupling on the vibration modes is identified, and the functional behavior of the frequencies with the lamination parameters is thoroughly investigated and discussed. Finite element modeling using NX Nastran solver is performed in order to validate the analytical results. As a case study, optimized blade spar designs of a 750-kW horizontal axis wind turbine are given. The attained solutions show that the approach used in this study enhances the dynamic characteristics of the optimized spar structures as compared with a known baseline design of the wind turbine blade.

Keywords

  • wind turbine blades
  • structural optimization
  • natural frequencies
  • advanced composites
  • sequential quadratic programming
  • finite element method

1. Introduction

Among all renewable energies of different styles, wind energy is the most popularized and potentially applicable type of green energy. Because larger wind turbines have more power capture and economic advantages, the typical size of utility-scale wind turbines, as shown in Figure 1, has grown dramatically over the last three decades [1, 2]. Such large flexible configuration, operating in uncertain environments, gives rise to significant vibration problems and assesses the importance of analyzing structural dynamics in the design of successful wind turbine systems. The main supporting structures of the rotating blades are usually fabricated from thin-walled composite beams with a variety of cross-sectional types. These configurations are used extensively in the design of many aerodynamic structures because of their light weight-to-stiffness ratio and long fatigue life. To design these components, the dynamic characteristics, especially near-resonant conditions, need to be well examined to assure a safe operation.

Figure 1.

Size and power increase of commercial wind turbines.

The objective of this investigation is to optimize the structural dynamics of a thin-walled composite blade through the minimization of structural mass or reduction of the overall vibration level. The latter can be attained directly by maximizing the natural frequencies of the main blade structure under strength and mass constraints. In general conditions, however, we need a material that is as light as possible for a specified stiffness in order to satisfy the design criteria and to minimize the weight-induced fatigue loads.

The main advantages of fiber composite materials [3] are their high strength and stiffness combined with low density, their superior fatigue properties due to the prevention of crack propagation, and their ability to tailor the layup for optimum strength and stiffness. However, sharp transitions between component materials may cause stress and strain discontinuities that facilitate failure [4]. A solution that can be promising to enhance dynamics and aeroelastic stability of composite blades is the use of functionally graded materials (FGMs), in which the mechanical and physical properties vary spatially within the structure. The concept of functionally graded materials was originated in Japan in 1984 during a space project, in the form of proposed thermal barrier material capable of withstanding high-temperature gradients [5]. FGMs may also be developed using fiber-reinforced layers with a changing volume fraction of fibers, rather than constant, producing grading of the material with favorable properties [6].

Considering, next, optimization of wind turbine blades, Maalawi and Negm [7] presented an optimization model for the design of a typical blade structure of horizontal axis wind turbines. The main blade spar was represented by thin-walled tubular beam composed of uniform segments with the design variables chosen to be the cross-sectional area, radius of gyration, and length of each segment. The optimal design is pursued with respect to maximum frequency design criterion subject to mass and aeroelastic constraints. The optimization problem was solved by multidimensional search techniques, where the aeroelastic stability boundaries and steady-state response were calculated using Floquet’s transition matrix theory. Another work by Maalawi [8] developed an optimization model for placing the frequencies of a wind turbine tower/nacelle/rotor structure in free yawing motion. The mathematical formulation considered a single pole tower configuration having thin-walled circular cross section with constant taper along the tower height. The nacelle/rotor combination was modeled as a rigid mass elastically supported at the top of the tower by the torsional spring of the yawing mechanism. The resulting governing differential equation of motion was solved analytically by transforming it into a standard form of Bessel’s equation, which leads to the necessary exact solutions for the frequencies and mode shapes. Useful design charts were developed for placing the frequencies at their needed target values with no penalty of increasing the total structural weight of the system.

In the context of using the concept of material grading, Librescu and Maalawi [9] formulated an analytical approach for obtaining the optimal design of a class of solid nonuniform composite wings with improved aeroelastic stability. The objective function was measured by maximizing the wing divergence speed while maintaining the total structural mass at a value equal to that of a known baseline design. Exact solutions were obtained for different categories of unidirectionally reinforced composite wing structures: namely, the linear volume fraction (L-VF), the parabolic volume fraction (PR-VF), and the piecewise volume fraction (PW-VF) wing models. Results revealed that in general, the torsional stability of the wing can be substantially improved by using nonuniform, functionally graded composites instead of the traditional ones having uniform volume fractions of the constituent materials.

The optimization of large wind turbines was considered by Kun-Nan Chen and Pin-Yung Chen [10], who applied a two-step procedure for finding the optimum design of composite blades. The first step concerned the optimal aerodynamic shape of the blade as described by the chord and twist angle distributions in the spanwise direction. The second step yielded the optimal material distribution. A 3-MW wind turbine with blades having cross sections of NREL S818, S825, and S826 airfoil types is demonstrated as a case study. A parameterized finite element model of the aerodynamically optimized blade was created using the ANSYS software. The optimization results showed that the initial blade model is an infeasible design due to a high level of the maximum stress, exceeding the upper limit of the stress constraint, but eventually the process converges to a feasible solution with the expense of increased total mass of the blade. Another work by Maalawi and Badr [11] considered the excessive wind turbine blade vibrations induced by continuous pitching, which is necessary to limit the power output and protect the generator from damages in severe wind conditions. They utilized analytical Bessel’s functions of the first kind which yield to the exact solutions of the resulting governing differential equation. The associated optimization problem was formulated by considering two forms of the objective function. The first one was represented by a direct maximization of the fundamental frequency, while the second one considered minimization of the square of the difference between the fundamental frequency and its target or desired value. In both strategies, an equality constraint is imposed on the total structural mass in order not to violate other economic and performance requirements. Design variables encompass the blade tapering ratio, chord, and shear wall thickness distributions. Danny Sale et al. [12] developed a numerical methodology for the structural analysis and optimization of composite blades for wind and hydrokinetic turbines. They derived a structural mechanic model which is based upon a combination of classical lamination theory with a Euler-Bernoulli and shear flow theory applied to composite beams. The development of this simplified structural model was motivated by the need for an accurate and computationally efficient method that is suitable for parametric design and optimization studies of composite blades. An important characteristic of this structural model is its ability to handle complex geometric shapes and isotropic or anisotropic composite layups. For a specified design load, the objective of the structural optimization was to minimize the blade’s mass while satisfying constraints on maximum allowable stress, blade tip deflection, buckling, and placement of blade natural frequencies. Adam Chehouri et al. [13] presented an improved version of the preliminary optimization tool called CoBlade, which offers designers and engineers an accelerated design phase by providing the capabilities to rapidly evaluate alternative composite layups and study their effects on static failure and fatigue of wind turbine blades. In this study, the optimization formulations included nonlinear failure constraints, and a comparison between three formulations was made to show the importance of choosing the blade mass as the main objective function and the inclusion of failure constraints in the wind turbine blade design.

A recent work by Maalawi [14] introduced a mathematical model for optimizing dynamic performance of thin-walled functionally graded box beams with closed cross sections. The objective function was measured by maximizing the natural frequencies and places them at their target values to avoid the occurrence of large amplitudes of vibration. Variables include fiber volume fraction, fiber orientation angle, and ply thickness distributions. Various power law expressions describing the distribution of the fiber volume fraction have been implemented, where the power exponent was taken as a main optimization variable. The mass of the structure is kept at a value equal to that of a known reference beam. Side constraints were also imposed on the design variables in order to avoid having unacceptable optimal solutions. A case study, including optimization of a cantilevered, a single-cell spar beam made of carbon/epoxy composite was considered. Conspicuous design charts were developed, showing the optimum design trends for the mathematical models implemented in the study. It was shown that the developed mathematical models are adequately satisfying the required global optimization of typical composite, functionally graded, thin-walled beam structures.

This chapter focuses on the optimization of the main structure of a wind turbine blade by either minimizing structural mass under frequency and strength constraints or maximizing the natural frequencies under mass, strength, and side constraints. This model is applied to tapered, anisotropic spar beam with thin-walled closed cross section made of laminated fibrous composites with variable thickness and stiffness. The study is focused on the spar structure that represents the main load-carrying component of the wind turbine blade. Material grading concept is utilized by changing the fiber content throughout the blade structure. Design variables include the volume fraction distribution of the constituent materials of construction and geometric and cross-sectional parameters of the blade spar. The blade is assumed to have a large span-to-chord ratio, which enabled us to model the main spar as an equivalent straight beam, positioned along the elastic axis. Structural analyses are performed using simplified mathematical expressions by implementing the conventional beam and classical lamination theories. The governing differential equations of motion are derived and solved by the transfer-matrix method for the coupled extensional-torsional and flexural-torsional modes of vibration. A case study is given considering thin-walled blade spar of a 750-KW horizontal axis wind turbine. Numerical results are presented and discussed showing the success of the developed mathematical model in producing efficient blade designs with improved dynamic performance. Finally, the relevant concluding remarks and recommendations for future studies are given and discussed.

2. Governing differential equations of motion

Among the dynamic characteristics of the blade main structures, determination of the natural frequencies and the associated mode shapes is of fundamental importance. An analytical model for the free vibration of anisotropic thin-walled beams with closed cross sections was developed by Armanios and Badir [15] using a variational asymptotic approach and Hamilton’s principles. This model was applied to arbitrary closed cross sections made of laminated fibrous composites with variable thickness and stiffness. The analysis was applied to two kinds of laminated composites: circumferentially uniform stiffness (CUS) and circumferentially asymmetric stiffness (CAS). The model was also implemented in Refs. [16, 17] to investigate the influence of coupling on the free vibration of thin-walled composite beams. Shadmehri et al. [18] studied the static and dynamic characteristics of composite thin-walled beams that are constructed from a single-cell cross section. The structural model considered incorporated a number of nonclassical effects, such as material anisotropy, transverse shear, warping inhibition, nonuniform torsion, and rotary inertia. The governing equations were derived using extended Hamilton’s principle and solved using extended Galerkin’s method. Phuong and Lee [19] presented a flexural-torsional analysis of thin-walled composite box beams. A general analytical model applicable to thin-walled composite box beams subjected to vertical and torsional loads was developed. Analytical solutions for the free vibration analysis of tapered thin-walled laminated composite beams with closed cross sections are given in Ref. [20]. The exact values of frequencies were obtained by means of power series schemes. A parametric analysis was performed for different taper ratios, stacking sequences, and materials.

Considering functionally graded constructions, Kargarnovin and Hashemi [21] investigated the free vibration of a fiber composite cylinder, in which the volume fraction of fibers varies longitudinally, using a semi-analytical method. The distribution of volume fraction of fiber in base matrix was based on power law model. Another study by Liu and Shu [22] developed an analytical solution to study the free vibration of exponential functionally graded beams with a single delamination. They showed that the natural frequencies increase as Young’s modulus ratio of the constituent materials becomes bigger.

Figure 2 shows the structural model of the blade spar, which is represented by a thin-walled cantilever beam which consists of Ns uniform segments. Each segment has different dimensions and material properties that satisfy the geometrical tapering and material grading distribution. Any segment k with length L k has a rectangular cross section with dimensions, width b k , depth a k , and wall thickness H k . Each segment is a uniform laminated fibrous composite beam which consists of N r layers, each of which has thickness h j , fiber volume fraction V fj , and fiber orientation angle θ j (j = 1,2,…., Nr).

Figure 2.

Blade spar structural model.

The constitutive relationships in terms of stress resultants and kinematic variables are [15, 16]:

TMxMyMz=C11C12C13C14C12C22C23C24C13C23C33C34C14C24C34C44U1φU3U2E1

where Tis the tensile force, Mxis the torsional moment, and Myand Mzare the bending moments about the y and z axes, respectively. Cmnare called the beam cross-sectional stiffness coefficients, and U1, U2, and U3are the average cross-sectional displacements along x, y, and z coordinates, respectively, and φxis the elastic twist about the x axis. The prime denotes differentiation with respect to x. Applying Hamilton’s principle, the equations of undamped free vibration are [15]:

C11U1+C12φ+C13U3+C14U2mU1´=0C12U1+C22φ+C23U3+C24U2Iφ´SzU3´+SyU2´=0C13U1+C23φ+C33U3+C34U2+Szφ´+mU3´=0C14U1+C24φ+C34U3+C44U2Syφ´+mU2´=0E2

where m, I, and Szand Syare the mass, polar, and first moments of inertia per unit length of the beam, respectively. The dot superscript denotes differentiation with respect to time.

A closed-form solution for the most general case of the equations of motion (Eq. 2) is not available. Two particular cases of fiber layup are considered in which some of the stiffness coefficients vanish. The first case is called circumferentially uniform stiffness (CUS) and the second circumferentially asymmetric stiffness (CAS). Figure 3 shows a rectangular cross-sectional beam segment with both CUS and CAS layup configurations. CUS layup configuration is manufactured by warping the composite layup using filament winding technique, θz=θz, while CAS layup configuration is manufactured such that the beam cross section is symmetric about the OXY plane. θz=θz.

Figure 3.

Spar segment with (a) CUS and (b) CAS layup configurations.

2.1 CUS layup configuration

For any segment k within the blade structure, a local coordinate system with X´-axis is introduced such that (see Figure 2).

0x´=XXkLkE3

In the special case of CUS layup, the equations of motion reduce to (for the kth segment):

C11kU1+C12kφmkU1´=0C12kU1+C22kφIkφ´=0C33kU3+mkU3´=0C44kU2+mkU2´=0E4

The first two equations express a coupled extension-twist vibration (ETV) mode, while the third and fourth equations express vertical bending vibration (VBV) and horizontal bending vibration (HBV) modes, respectively.

Applying the integral formulas given in Ref. [15], the stiffness coefficients can be determined from:

C11k=2KAkbk+ak2HkC12k=KBkbkHkakHkC22k=KCkbkHkakHk22bk+ak2HkC33k=KAkKBk2KCkakHk361+3bkHkakHkC44k=KAkKBk2KCkbkHk361+3akHkbkHkE5

where KA, KB, and KCare the reduced axial, coupled axial-shear, and shear stiffness coefficients, respectively, given by the expressions:

KA=A11A122A22KB=2A16A12A26A22KC=4A66A262A22E6

where Amn are called the membrane in-plane stiffnesses, which depend on the fiber orientation angle θ, volume fractions, and mechanical properties of the fiber and matrix materials [4].

The Young’s moduli in the longitudinal and lateral directions of the lamina E11 and E22 , the shear modulus G12 , and the major Poisson’s ratio v12 are calculated using the semiempirical methods by Halpin and Tsai [23]:

E11=Em1Vf+EfVfE22=Em1+ηEVf/1ηEVfG12=Gm1+ηGVf/1ηGVfν12=νm1Vf+νfVfηE=EfEm/Ef+EmηG=GfGm/Gf+GmE7

Subscripts “m” and “f” refer to the properties of matrix and fiber materials, respectively, and Vƒ is the volume fraction of fibers within each lamina. Considering the equations of coupled extension-twist vibration, the assumed solution is [24]:

U1xt=C´1eλx+iωtφxt=C´2eλx+iωtE8

where ωis the circular natural frequency of free vibration. Substituting from Eq. (8) into Eq. (4), the associated characteristic equation can be shown to be:

aλ4+bω2λ2+cω4=0a=C11kC22kC12k2b=C11kIk+C22kmkc=mkIkE9

The solution can be expressed as:

λ1,2=±iα1λ3,4=±iα2α1,2=ωbb24ac2aE10

If C11kIk>C22kmk, then the minus sign in Eq. (10) will generate natural frequencies with axial-mode dominated vibration, while the plus sign will generate natural frequencies with torsion-mode dominated vibration. If C11kIk<C22kmk, the domination is reversed. The general solution of the coupled extension-twist vibration problem (Eq. 4) can be expressed as:

U1xt=c1sinα1x+c2cosα1x+c3sinα2x+c4cosα2xeiωtφxt=c5sinα1x+c6cosα1x+c7sinα2x+c8cosα2xeiωtE11

In order to satisfy both the coupled equations of motion, there are specific relations between the constants such that:

c5=q1c1,c6=q1c2,c7=q2c3,c8=q2c4q1,2=2mkaC12kbb24acC11kC12kE12

2.2 CAS layup configuration

In this case, the equations of motion, applied to segment k , reduced to:

C11kU1mkU1´=0C22kφ+C23kU3Ikφ´=0C23kφ+C33kU3+mkU3´=0C44kU2+mkU2´=0E13

The second and third equations of motion express a coupled bending-twist vibration (BTV) mode, while the first and fourth equations of motion express extension vibration (EV) and horizontal bending vibration (HBV) modes, respectively. The non-zero stiffness coefficients are given by:

C11k=2KAkbk+ak2Hk2KBk2KCkbkHkC22k=KCkbkHkakHk22C23k=KBkbkHkakHk22C33k=KAkakHk361+3bkHkakHkKBk22KCkC44k=bkHk361+3akHkbkHkKAkKBk2KCkE14

Assuming harmonic solution similar to that given in Eq. (8), the characteristic equation is obtained by differentiating Eq. (13) to get:

aλ6+bω2λ4cω2λ2dω4=0a=C22kC33kC23k2b=C33kIkc=C22kmkd=mkIkE15

which has the solution:

λ1,2=±iγ1,λ3,4=±iγ2,λ5,6=±γ3E16

The general solution takes the form:

U3xt=c1sinγ1x+c2cosγ1x+c3sinγ2x+c4cosγ2x+c5sinhγ3x+c6coshγ3xeiωtφxt=c7sinγ1x+c8cosγ1x+c9sinγ2x+c10cosγ2x+c11sinhγ3x+c12coshγ3xeiωtE17

The relations between the constants are given by:

c7=k1c2,c9=k2c4,c11=k3c6c8=k1c1,c10=k2c3,c12=k3c5k1=C23kγ13C22kγ12Ikω2k2=C23kγ23C22kγ22Ikω2k3=C23kγ33C22kγ32+Ikω2E18

3. Natural frequencies and mode shapes

Applying the transfer-matrix method [25], the relation between the state vectors at one end of a segment Zkand that at the other end Zk+1is given by:

Zk+1=TkZkE19

where Tkis called the elementary transfer matrix associated with the kth segment. For a spar built up of Ns uniform segments, Eq. (19) can be applied at successive stations to obtain:

ZNs+1=ToZ1To=TNsTNs1TkT2T1E20

where the matrix Tois known as the overall transfer matrix, which relates the state vectors at the spar fixed end to the free end at which the boundary conditions are specified. Therefore, applying the boundary conditions at both ends and considering only the nontrivial solutions, the frequency equations can be readily obtained. Derivations of the elementary transfer matrices for different vibration modes are discussed in the following section.

3.1 CUS layup

In this case the elements of the state vector at both ends of the k th segment are related by the elementary transfer matrix for the coupled extension-twist vibration as follows:

U1LkφLkC11U1Lk+C12φLkC22φLk+C12U1Lk=T11T12T13T14T21T22T23T24T31T32T33T34T41T42T43T44kU10φ0C11U10+C12φ0C22φ0+C12U10E21

Imposing the cantilevered boundary conditions on the overall transfer matrix, the overall frequency equation for coupled extension-twist vibration of the blade spar is given by:

T33oT44oT34oT43o=0E22

3.2 CAS layup

Here the transfer matrix equation for the segment k is:

U3LkφLkU3LkC22kφLk+C23kU3LkC23kφLk+C33kU3LkC33kU3Lk=T11T12T13T21T22T23T31T32T33T14T15T16T24T25T26T34T35T36T41T42T43T51T52T53T61T62T63T44T45T46T54T55T56T64T65T66kU30φ0U30C22kφ0+C23kU30C23kφ0+C33kU30C33kU30E23

Imposing the cantilevered boundary conditions on the overall transfer matrix, the frequency determinant is given by:

T44T45T46T54T55T56T64T65T66o=0E24

The natural frequencies for coupled bending-twist vibration can be obtained numerically by solving the characteristic determinant of Eq. (24) for γ.

4. Optimization model formulation

4.1 Selection of design variables and pre-assigned parameters

In order to formulate a practical and inexpensive optimization model, the large number of design variables of a wind turbine blade has to be reduced to a reasonable number that can be easily dealt with. The design variables which are not subject to change in the optimization process are called the pre-assigned parameters. They are selected to be:

  1. Type of materials of construction

  2. Total blade length

  3. Chord distribution along the blade axis

  4. Twist angle distribution along the blade axis

  5. Airfoil type and dimensions

  6. Dimensions and spacing of the internal supporting ribs

  7. Covering skin thickness

  8. Blade-to-hub attachment, which is chosen to be of hingeless type

  9. Shape of the spar cross section, which is chosen to be rectangular

On the other hand, the design variables, which are subject to change during the optimization process, are chosen to be the fiber orientation angle θj,k , fiber volume fraction Vfj,k , thickness of each lamina hj,k , cross-sectional dimensions (ak , bk , Hk ), and length Lk of each segment composing the blade spar (refer to Figure 2). The cross-sectional aspect ratio (ak /bk ) can be prescribed according to the airfoil cross-sectional dimensions, and only the variable ak is taken as a design variable. In addition, to formulate a normalized optimization model with scaled variables, it is necessary to start with a known reference beam to which all the design variables, constraints, and objective function are referred. It is chosen to be made of unidirectional laminated composites with fiber volume fraction Vfo and fiber orientation angle θo , and its cross-sectional dimensions are denoted by bo for width, ao for depth, Ho for wall thickness, and Lo for length. Therefore, the normalized variables are defined as follows:

L̂k=Lk/Lob̂k=bk/boâk=ak/aoĤk=Hk/Ho=j=1Nrĥjĥj=hj/HoE25

For a spar beam composed of Ns segments, the design variables are defined by the following matrix equation:

X=L̂1L̂2L̂kL̂Nsâ1θj,1â2θj,2âkθj,kâNsθj,NsVfj,1Vfj,2Vfj,kVfj,Nsĥj,1ĥj,2ĥj,kĥj,NsθNr,1θNr,2θNr,kθNr,NsVfNr,1VfNr,2VfNr,kVfNr,NsĥNr,1ĥNr,2ĥNr,kĥNr,NsE26

The total number of design variables in this case is equal to (3 Nr + 2)*Ns, which is a relatively large number of optimization variables.

4.2 Different optimization strategies

Several criteria are used in the dynamic optimization of wind turbine blades; some important ones are the minimal mass design, maximum frequency design, and frequency-placement criteria. In the first strategy, the mass of the blade is minimized while imposing constraints on the blade natural frequencies and strength. In the maximum frequency criterion, the reduction of vibration level is attained by maximizing the natural frequencies of the blade without regard to the complicated stiffness/mass ratio and the exciting frequencies constraints. Higher natural frequencies are favorable for reducing both of the steady-state and transient responses of any structure being excited. The last category of the objective functions to be considered is the placement of the blade natural frequencies, where the main goal is to separate the natural frequencies of the blade from the exciting frequencies to avoid large amplitudes from occurring near the resonant conditions.

4.2.1 Minimal mass design

The minimal mass design optimization problem can be stated as follows:

FindXWhichminimizesfX=m̂Subject to:Strength constraintResonance avoidanceSide constraintsαT101nω̂nω̂n1+n,n=1,2,XLXXUE27

where m̂=m/moand ω̂n=ωn/ωnoare the dimensionless mass and frequency and moand ωnoare the mass and frequency of the reference beam, respectively. XLXUare the lower and upper bounds imposed on the design variables.

Structural mass:m=2k=1nsj=1nrρjhjbk+ak2HkLkmo=2ρoHobo+ao2HoLoE28

The mass densities ρoand ρjare calculated according to the volume fractions Vfo and Vfj , respectively. The symbolαTis called the Tsai-Hill rupture coefficient [26], which depends on the ratios between the principle and rupture stresses (more details are given in appendix A). ω̂n, n = 1,2,… are called the target (or desired) frequencies of a known baseline design adjusted to be well separated from the exciting frequencies, and the increments nare the associated allowable tolerance of each frequency (e.g., n= 1%).

4.2.2 Maximum frequency optimization

Minimization of the overall vibration level is one of the most cost-effective solutions for a successful wind turbine design. It enhances other important design goals such as long fatigue life, high stability, and low noise level. Reduction of the overall vibration level can be attained by maximizing the stiffness-to-weight ratio of the wind turbine blade spar. It is well-known that natural frequency is a good indicator of structural stiffness-to-mass ratio. Thus, the optimization problem considered in this investigation will seek maximization of natural frequencies of the blade spar for different modes of vibrations while maintaining its total structural mass lower than or equal to that of a baseline design. Constraints are imposed on the optimization problem such that the blade spar has to be enclosed by the outer blade skin and has the sufficient strength to carry the applied loads without mass penalty.

The maximum frequency criterion, taking into considerations the different modes of vibration, may be cast in the following form:

FindXWhichminimizesfX=nWnω̂n,n=1,2,Subject to:Strength constraintMass constraintSide constraintsαT10m̂10XLXXUE29

The function f(X) represents a weighted sum of the nondimensional natural frequencies of successive modes of vibration with Wnthe weighting factor of the n th frequency ω̂n.

4.2.3 Frequency placement

A good design philosophy for vibration reduction is to separate the natural frequencies of the structure from the harmonics of other excitation sources. Thus, natural frequencies of the blade spar should be placed near a target value. Mathematically, the frequency-placement criterion may be expressed by minimizing an objective function constructed from a weighted sum of the squares of the differences between each important frequency ω̂nand its target value ω̂nas follows:

FindXWhichminimizesfX=nWnω̂nω̂n2Subject to:Strength constraintMass constraintSide constraintsαT10m̂10XLXXUE30

4.3 Side constraints

Side constraints are imposed on the design variables in order to verify various geometric, manufacturing, or logical reasons. In the present optimization model, these constraints are defined in the following.

4.3.1 Length of blade spar

The total length of the blade is kept equal to that of a known baseline design:

k=1NsL̂k1=0L̂k0E31

4.3.2 Spar cross-sectional dimensions

Since the blade spar is limited by a tapered configuration with a certain type of airfoil cross sections, the height of any segment k having its end at a specific position xk+1 (refer to Figure 2) must not exceed the outermost height (aU ) at this position. Figure 4 shows the outermost dimensions of a wind turbine blade spar of length L, where ar and at represent the spar heights at the root and tip locations, respectively. The dimensionless upper limiting height âUis defined by the relation:

Figure 4.

Outermost dimensions of a wind turbine blade spar.

âU=âr11atarx̂k+1x̂k+1=L̂1+L̂2++L̂kE32

Lower bounds are also imposed on the height of the k th segment as a reasonable percentage of the airfoil height near the tip (âL=0.5ât). Thus, the inequality constraints imposed on the dimensionless height of each segment of the blade spar are expressed as:

âLâkâU,k=1,2,.NsE33

Upper bounds on the segment height represent the main obstacle in blade spar optimization since it restricts the spanwise variation in spar height. Note that the width of each segment (bk ) depends on the height of that segment (ak ) according to the pre-assigned aspect ratio of the cross section.

In addition, the inequality constraints imposed on the walls thickness of each segment are:

ĤLĤkĤU,k=1,2,.NsE34

4.3.3 Fiber volume fraction and orientation angle within any segment (k)

0.2π/2Vfjθj0.8π/2,j=1,2,NrE35

The optimization models described by Eqs. (27–35) belong to the constrained nonlinear mathematical programming models since the objective functions and most of the constraints contain nonlinear algebraic expressions. Such problems can be solved by a variety of mathematical programming techniques [27, 28]. The sequential quadratic programming (SQP) is one of the most powerful developed and perhaps one of the best methods of optimization. The method has a theoretical basis that is related to the solution of a set of nonlinear equations using Newton’s method and the derivation of simultaneous nonlinear equations using Karush-Kuhn-Tucker (KKT) conditions to the Lagrangian of the constrained optimization problem [29, 30]. More details are given in Appendix B.

5. Behavior of the objective function

Focusing on CUS layup configurations, extensive studies have been carried out, using MATLAB, on the objective functions in order to be able to visualize the unconstrained behavior of natural frequencies and mass of specific thin-walled composite beams. Studies are performed in a two-dimensional design space such that only two design variables are allowed to change while the others are assigned to specific values. All variables are normalized with respect to a known reference beam parameters, as given in Eq. (25).

5.1 Case of two-segment cantilever beam

For a two-segment beam with each segment constructed only from one layer, the matrix of variables (Eq. 26) can be reduced to take the form:

X=L̂1â1θ1,1Vf1,1ĥ1,1L̂2â2θ1,2Vf1,2ĥ1,2E36

Note that in this case Ĥk=ĥj. The effects of changing two corresponding variables simultaneously in the two segments are studied in the following.

5.1.1 Effect of changing fiber orientation angles within the two segments (θ1, θ2)

Consider a two-segment beam with a height of the outboard segment half that of the inboard segment but with equal length. Fiber volume fraction and wall thickness of the two segments are the same as that of the reference beam. The following matrix of pre-assigned variables describes this case:

X=0.51.5θ10.510.50.75θ20.51E37

Figure 5 shows the developed contours of the normalized fundamental bending frequency in the (θ1, θ2)-design space. It is noticed that the maximum bending natural frequencies occurs at the point θ1θ2=00, which is an expected natural solution. Figure 6 shows the level curves of the normalized twist-dominated natural frequency. It indicates that the maximum occurs near the design point θ1θ2=4545. On the other hand, the design points (θ1, θ2) = (0, 0), (90, 90), (0, 90), and (90, 0) give the same values of the frequency representing a local minima solution, as shown in Figure 6. The angle sequence in this case affects the frequency level curves due to the change in heights of the two segments.

Figure 5.

Effect of fiber orientation angles on the normalized fundamental bending frequency (a two-segment cantilevered beam).

Figure 6.

Effect of fiber orientation angles on normalized twist-dominated natural frequency (a two-segment cantilevered beam).

5.1.2 Effect of changing length and height of the outboard segments (L̂2, â2)

Consider a two-segment beam with the following matrix of pre-assigned variables. The height and thickness of the inboard segment are equal to those of the reference beam, while thickness of the outboard segment is half that of the reference beam. Nondimensional length of the inboard segment is 1L̂2in order to keep the total length equal to that of the reference beam:

X=1L̂2100.51L̂2â200.50.5E38

Figure 7 shows the developed contours of the fundamental bending frequency in the (L̂2, â2)-design space. It is seen that the frequency domain is divided into three regions bonded by the reference frequency value. The extreme left and right regions contain two local minima. The central region contains the global maximum. The maximum bending natural frequencies occur near the point (L̂2, â2) ≈ (0.35, 0.3). This point corresponds to a total structural mass of about 70% of the reference beam, which represents a significant gain from the optimization process. This configuration of beam dimensions gives a tapered beam with a decreasing height and length. The same value of maximum natural frequency is also repeated in the design space at the uppermost right region at (L̂2, â2) ≈ (1, 2), but the total mass becomes higher and equal to that of the reference beam.

Figure 7.

Effect of length and height of the outboard segment on normalized fundamental bending frequency (a two-segment cantilevered beam).

From the previous investigations, it is observed that the frequency functions, even though implicit in the design variables, are well behaved and continuous in the selected design spaces. A weighted sum of successive natural frequencies is a good representative for the objective function of the optimization model. Another important observation is that the length and the height of the individual segments are so effective variables in the whole optimization process.

6. Case study

As a case study, the main spar of a medium scale composite blade of a 750-kW horizontal axis wind turbine is optimized. Full description and technical data can be found in Ref. [31]. The blades cross section having NACA 63-218 airfoil is made of E-glass/epoxy composites with properties given in Table 1. The effective length of the spar is 18.33 m extending from 20–98% of the blade length. The spar cross-sectional dimensions at root are 340 mm height and 680 mm width and at the tip 100 mm and 200 mm. The wall thickness at the root and tips are 15 mm and 5 mm, respectively. The chordwise location of the spar cross-sectional center coincides with the blade pitch axis.

6.1 Baseline design

An initial baseline design of the wind turbine blade spar with total mass = 580.55 kg is selected to be composed of 10 segments of equal length. Fiber orientation angle and volume fraction within each segment are pre-assigned to be 0° and 50%, respectively. Structural and dynamic properties of the spar baseline design are given in Tables 2 and 3, respectively.

PropertyEpoxy matrixE-glass fibers
Modulus of elasticity (GPa)Em = 4.5E11f = E22f = 74
Modulus of rigidity (GPa)Gm = 1.6G12f = 30
Poisson’s ratioνm = 0.4ν12f = 0.25
Density (kg/m3)ρm = 1200ρf = 2600
Tensile strength (MPa)σmr = 90σ11r ,σ22r (Eq. A-2)
Shear strength (MPa)τmr = 52

Table 1.

Material properties of the spar structure [4, 31].

Segment no.Length (mm)Height (mm)Width (mm)Thickness (mm)
1183331663214.4
229258413.3
326853612.2
424448811.1
522044010
61963928.9
71723447.8
81482966.7
91242485.6
101002004.5

Table 2.

Dimensions of the baseline design.

Vibration mode (dominated)Natural frequencies (Hz)
First modeSecond modeThird mode
Twist36.7364.8895.18
Extension94.9203.5319.8
Flap bending1.554.9011.15

Table 3.

Natural frequencies of the baseline design.

The calculated normal and tangential distributed forces (Pz , Py ) and the corresponding bending moments (Mz , My ) applying to the wind turbine blade are shown in Figures 8 and 9, respectively (refer to Appendix C). The wind speed of 60 m/s (survival loading condition) is taken in calculating the applied loads, while the blades were considered in the stationary horizontal position. These forces and moments are implemented in the strength constraints of the present optimization model.

Figure 8.

Shear forces applied to the blade at 60-m/s wind speed.

Figure 9.

Flapwise and chordwise bending moments applied to the blade at 60-m/s wind speed.

6.2 Definition of the reference beam

The reference beam, to which all the design variables and frequencies are normalized, is selected to be a uniform, thin-walled cantilevered spar with rectangular cross section constructed from single unidirectional E-glass/epoxy composite layer (θo = 0) and fiber volume fraction of Vfo = 50%. This fiber layup configuration gives maximum natural frequencies for bending and extension modes of vibration while maintaining moderate structural mass (mo = 580.5 kg) and strength. The dimensions of the reference beam are selected according to the mean cross-sectional dimensions of the spar baseline design (see Table 4).

Dimensions (mm)Natural frequencies: ωno/2(Hz)
ModeFirstSecondThird
Width bo = 440
Height ao = 220Twist17.5952.7687.93
Length Lo = 18,330Extension61.99185.97309.95
Thickness Ho = 10Bending0.654.0711.39

Table 4.

Dimensions and natural frequencies of the reference beam.

6.3 Optimization results

The blade spar described in Section 6.1 is optimized according to the three different optimization models under the specified constrains as given in Section 4.2. The method of sequential quadratic programming with MATLAB optimization toolbox is applied to obtain the needed optimal solutions (refer to Appendix B). The starting, lower, and upper bounds for each design variable are given in Table 5. Determination of these values affects the optimization results; thus, several trials were performed in order to select the appropriate values.

6.3.1 Minimum mass design

In this formulation the target frequencies are set equal to the first three flapping frequencies of the baseline design, which are well separated from resonant frequencies. The attained optimal solutions are given in Table 6 with the spar cross sections constructed from balanced and symmetric laminates. The achieved mass saving reached about 25%, where the natural frequencies are found to be too close to those of the baseline design, as given in Table 7. Considerable mass reduction can be observed with much savings for the inboard segments as compared to that for the outboard segments.

Design variableLower boundaryUpper boundaryStarting value
Dimensionless length of each segment L̂k0.011/Ns
Dimensionless height of each segment âk0.221
Fiber orientation angle θj,k−π/2π/20
Fiber volume fraction Vf j,k0.20.80.5
Dimensionless thickness of each layer ĥj,k0.010.20.1

Table 5.

Lower, upper, and starting boundaries of design variables for multiple-segment spar.

Segment no.Segment length (mm)Segment cross section (mm)Fiber volume fractionNo. of layersLayup
HeightThickness
1280.5323.89.90.3118(±45/03/±45/03/±45)s
2183.3316.79.90.3318(±45/03/±45/03/±45)s
3200.8313.88.90.3216(±45/03/±45/02/±45)s
4534.1301.48.90.3316(±45/02/±45/03/±45)s
5851.0286.38.90.3216(±45/03/±45/02/±45)s
61031.7263.87.90.3214(±45/02/±45/02/±45)s
71112.0235.27.80.3314(±45/02/±45/02/±45)s
81272.3207.47.80.3414(±45/05/±45)s
91482.5188.96.80.3312(±45/05)s
101589.4170.46.50.3312(±45/05)s
111273.1154.06.50.3312(±45/05)s
121132.9138.66.40.3212(±45/05)s
131100.9120.95.70.3310(±45/04)s
141051.8106.35.60.3310(±45/04)s
151177.694.15.20.338(±45/03)s
16880.582.44.80.338(±45/03)s
17693.270.74.60.338(±45/03)s
18540.059.04.50.348(±45/03)s
19456.848.04.00.336(±45/02)s
201485.744.03.60.206(±45/02)s

Table 6.

Minimum mass design optimization.

Vibration mode (dominated)DesignNatural frequencies (Hz)n=13Wnωn
First modeSecond modeThird mode
TwistBaseline36.7364.8894.9050.49
Min. mass36.1665.6289.1349.52
ExtensionBaseline94.90203.50319.80148.00
Min. mass89.13194.88299.66139.77
FlappingBaseline1.554.9011.152.48
Min. mass1.534.8811.172.45
Mass (kg)Baseline580.55
Min. mass437.76

Table 7.

Comparison of natural frequencies and mass between baseline and the minimum mass designs.

6.3.2 Maximum frequency design

In order to obtain a blade spar with higher natural frequencies, the maximum frequency optimization model given in Eq. 29 is applied under inequality mass constraint. Table 8 presents the attained optimal solutions for the blade spar with discrete grading model.

Segment no.Segment length (mm)Segment cross section (mm)Fiber volume fractionNo. of layersLayup
HeightThickness
184232918.70.834(±45/07/±45/07/±45)s
287331816.30.830(±45/06/±45/06/±45)s
390330614.40.826(±45/05/±45/05/±45)s
491929412.30.822(±45/04/±45/04/±45)s
5932281110.820(±45/03/±45/04/±45)s
69372698.80.816(±45/03/±45/02/±45)s
79302577.70.814(±45/02/±45/02/±45)s
89082456.60.812(±45/04/±45)s
98822345.50.810(±45/04)s
108612225.50.810(±45/04)s
118752115.50.810(±45/04)s
129361995.50.810(±45/04)s
1310351855.50.810(±45/04)s
1411281705.50.7810(±45/04)s
1510941565.50.6510(±45/04)s
16103114240.518(±45/03)s
1795613040.338(±45/03)s
187099130.336(±45/02)s
195436130.336(±45/02)s
2010374420.24(±45/0)s

Table 8.

Maximum frequency design with discrete spanwise material grading.

Figure 10 shows the optimization gains for the different natural frequencies of the blade spar with maximum frequency designs. It is noticed that gains of 48%, 23%, and 34% have been achieved in fundamental natural frequencies for twist, extension, and flapping modes of vibration, respectively.

Figure 10.

Maximum frequency optimization gains.

It can be concluded that a good blade design with maximum natural frequencies should have smaller wall thickness in the outboard portion of the blade spar. The maximum frequency optimization process recommends the segments located inboard to have higher wall thickness and higher fiber content. The spanwise variation in the height of each segment of the blade spar is always restricted by the airfoil envelope, which is the major obstacle in the optimization process.

When the wind turbine is operating, the rotating blades are the main source of vibration. The forcing frequencies are integer multiples of the rotational speed. A common way to represent natural frequencies and search for possible resonances is to plot the Campbell diagram as shown in Figure 11. The intersection of one of the radial lines with one of the system natural frequencies indicates a potential for resonance at the rotor speed corresponding to that point of intersection. The wind turbine considered in this investigation has two rotational speeds; they are 15 and 22 rpm. It is shown that for the baseline and maximum frequency designs, no resonance can occur at these operating speeds.

6.3.3 Frequency-placement design

Figure 11.

Campbell diagram for baseline and maximum frequency designs. (ـــــــBaseline, --- max. ∑ n ).

In the minimal mass design, the natural frequencies of the optimized blade spar are constrained to be equal to those of the baseline design, while in maximum frequency optimization, the natural frequencies are entirely maximized with constraints imposed on its structural mass. If it is desired to have a blade spar with specific target frequencies, the frequency-placement optimization model has to be applied. The desired target frequencies used in the present analysis are selected to be 1.9, 6.4, and 12.65 Hz for the first three modes of flapping vibration, respectively. These frequencies are placed away from resonant frequencies at the operating speeds according to the Campbell diagram shown in Figure 12. The attained optimal solutions of the blade spar are given in Table 9, with the flapping frequencies placed at their prescribed target frequencies.

Figure 12.

Campbell diagram for baseline and frequency-placement designs. (ــــــ baseline, ---- target).

Segment no.Segment length (mm)Segment cross section (mm)Fiber volume fractionNo. of layersLayup
HeightThickness
175333014.50.8026(±45/05/±45/05/±45)s
294831812.90.8024(±45/04/±45/05/±45)s
375430812.10.8022(±45/04/±45/04/±45)s
478729812.30.8022(±45/04/±45/04/±45)s
5100828411.90.8022(±45/04/±45/04/±45)s
654027511.10.7920(±45/03/±45/04/±45)s
746027110.80.7920(±45/03/±45/04/±45)s
82562579.90.7818(±45/07/±45)s
92162519.90.7818(±45/08)s
1018624110.90.7120(±45/09)s
1114312447.80.7414(±45/06)s
1225632106.40.8012(±45/05)s
1319301856.50.8012(±45/05)s
1415361656.40.8012(±45/05)s
1513741475.60.8010(±45/04)s
1613101256.40.4312(±45/05)s
17585775.30.7710(±45/04)s
181855110.90.7220(±45/09)s
19431715.50.3610(±45/04)s
201077445.50.6310(±45/04)s

Table 9.

Frequency-placement design with discrete spanwise material grading.

7. Finite element analysis

The finite element method (FEM) has been demonstrated as a powerful approach which can handle dynamic analysis of laminated composite structures. Nowadays, current commercial finite element software is capable of simulating nonlinearity of many engineering problems. Commercial software also come with advanced preprocessing and post-processing abilities [32] The preprocessing is just a way used for the data input, since the finite element method requires a large amount of data, while the post-processing is another way for presenting the results in the form of deformed shapes and contour maps. The core of the analysis is what occurs in between the two processes. In selecting a finite element analysis software, it is essential to regard the pre- and post-processor, which can directly affect the analysis speed and accuracy. One of the most appropriate finite element analysis software is Femap (finite element modeling and post-processing). Femap is an engineering analysis program developed by Siemens PLM Software that is used to build finite element models of complex engineering problems and view solution results [33]. Femap uses NX Nastran solver to analyze and solve finite element problems.

Femap is an advanced engineering simulation application for creating, editing, and importing/reusing mesh-centric finite element analysis models of complex products or systems. Femap provides powerful data-driven and graphical result visualization and evaluation.

The attained optimum design of the blade spar is modeled using FEM as a tapered beam with rectangular cross section. The root and tip cross sections are drawn using Femap, and then each line in the root is ruled to the corresponding line in the tip in order to generate the required surfaces as shown in Figure 13. The geometric model is divided into 20 segments with lengths according to the maximum frequency design given in Table 8.

Figure 13.

Finite element model of a 20-segment cantilevered spar meshed with 10,000 elements.

Material properties, element type, and thickness of each segment are to be defined. According to Table 8, a number of six different orthotropic materials with different fiber volume fractions and properties of E11 , E22 , G12 , ν12 , and ρ are to be defined and input to Femap. A total of 20 different layups are to be defined due to the change in fiber orientation angles and thickness of each segment within the spar design. Furthermore, 20 element types are to be defined according to the corresponding material and layup properties.

The normal mode eigenvalues analysis using Lanczos modal analysis is used by Femap in order to define the different mode shapes and the associated natural frequencies of the optimized blade spar design. Figure 14 shows the mode shapes with the associated natural frequencies for bending vibration of the optimized blade spar. Table 10 shows a comparison between analytical and Femap natural frequencies of the optimized blade spar design.

Figure 14.

Mode shapes and natural frequencies for vertical bending vibration of the blade spar with maximum frequency design.

AnalysisNatural frequencies (Hz)
1VBV2VBV3VBV
Analytical2.096.8713.30
FEM2.177.0414.49
Diff. %4.72.48.2

Table 10.

Comparison between analytical and FEM modeling for the maximum frequency design.

8. Conclusions and recommendations

This chapter presents an optimization model for enhancing the dynamic performance of the spar beam of a wind turbine blade. Design variables include the cross-sectional dimensions and material properties variation along the spanwise direction. Three optimization strategies are developed and tested, including the minimal mass design, maximum frequency design, and frequency-placement design, by placing the frequencies at their target values to avoid large amplitudes and resonance occurrence. Side constraints are imposed on the design variables in order to avoid abnormal-shaped optimized configurations. Based on the fact that an exact dynamic analysis of uniform thin-walled beam segment is available and well established, the dynamic analysis of tapered blade spar has been obtained by applying the transfer matrix method to calculate the natural mode of vibrations. The proposed model deals with dimensionless quantities in order to be applicable to thin-walled beams with arbitrary dimensions. Results indicated that the optimization process leads to significant increase of natural frequencies of the optimized spar when compared to the reference or baseline design without mass penalty. Finite element model showed a good agreement with the analytical model developed in this study with a variation of up to 10%. The main conclusions that can be revealed from the present work are:

  1. Tapered multiple-segment spar with spanwise material grading gives natural frequencies higher than that of the reference design. However, it is proved that maximization of the fundamental frequency alone does not guarantee maximization of the other higher frequencies. Higher frequencies have been found to have many local minima and maxima in the defined design space.

  2. There are optimum design variables of each segment such as the length, height, wall thickness, and fiber content at which the structural dynamic performance can be enhanced. Good designs favor minimum wall thickness and higher fiber volume fraction.

Finally, the analytical model formulated in this chapter can be extended and applied to study the forced dynamic response of a wind turbine blade. Other cross-sectional types of the blade spar, such as D-shape spar, can be considered, and the effects of blade twist, shear deformation, and rotary inertia are to be considered in future studies.

The fracture processes induced in fibrous composite materials depend upon the nature of constituents, the architecture of the laminate, and the type of mechanical loading imposed to the laminate. The rupture of fibrous composite materials is the result of one or combined effect of fiber fracture, transverse fracture in matrix, longitudinal fracture in matrix, fracture of fiber-matrix interface, and delaminations [34] (refer to Figure A.1).

Figure A.1.

Schematic of the fracture of a unidirectional fibrous composite at critical values of (a) axial, (b) transverse, and (c) shear stresses.

Figure A.2.

Lamina under loading condition of (a) principal axes and (b) general axes.

Failure theories for composites have been proposed by extending and adapting isotropic failure theories to account for the anisotropy in stiffness and strength of the composites. One of the first fracture criteria applied to an anisotropic materials was introduced by Tsai-Hill theory [26], which is applied to unidirectional lamina under principal axis in-plane loading condition as shown in Figure A.2 (a). The theory states that no fracture will occur in the lamina if αT1such that:

αT=σ11σ11r2+σ22σ22r2σ11σ22σ11r2+τ12τ12r2EA.1

where the subscript r refers to rupture strength of the material.

Rupture strengths for fibrous composite materials are given by the following approximate relations [34]:

σ11r=σ11frVf+VmEm/E11fσ22r=σmr14Vf/πτ12r=τmrEA.2

Subscripts f and m refer to fiber and matrix materials, respectively.

For a lamina under general loading condition as shown in Figure A.2 (b), the principal stresses are given by the following equation [4, 34]:

σ11σ22τ12=c2s22css2c22cscscsc2s2σxσyτxyc=cosθs=sinθEA.3

The Tsai-Hill failure theory is expressed in terms of a single criterion instead of the multiple subcriteria required in the maximum stress and maximum strain theories. The theory allows for considerable interaction among the stress components. One disadvantage, however, is that it does not distinguish directly between tensile and compressive strengths. The strength must be specified and used according to the given state of stress.

Sequential quadratic programming (SQP) is one of the mainly developed and perhaps one of the best techniques of optimization [29]. The method has a theoretical base that is related to the solution of a group of nonlinear equations using Newton’s method and the derivation of concurrent nonlinear equations using Karush-Kuhn-Tucker (KKT) conditions to the Lagrangian of the constrained optimization problem. The fundamental idea of SQP is to model the optimization problem at the present iterate xk by a quadratic programming subproblem and to employ the minimizer of this subproblem to identify a new iterate x k + 1 [30]. The challenge is to design the quadratic subproblem so that it yields a good step for the constrained optimization problem and so that the overall SQP algorithm has good convergence properties and good practical performance. Possibly the simplest derivation of SQP methods views them as an application of Newton’s method to the KKT optimality conditions for the optimization problem. If the problem is a so-called convex programming problem, that is, f(X) and Gi ( X ), i = 1,…,m, are convex functions, then the KKT equations are both necessary and sufficient for a global solution point. The Kuhn-Tucker equations can be stated as:

fx+i=1mλi.Gix=0λi.Gix=0,i=1,.,meλi0,i=me+1,,mEB.1

The first equation describes a canceling of the gradients between the objective function and the active constraints at the solution point. For the gradients to be canceled, Lagrange multipliers (λi, i = 1,…,m) are necessary to balance the deviations in magnitude of the objective function and constraint gradients. Because only active constraints are included in this canceling operation, constraints that are not active must not be included in this operation and so are given Lagrange multipliers equal to 0. This is stated implicitly in the last two Kuhn-Tucker equations.

B.1 Optimization using MATLAB

MATLAB is popular software that is used for the solution of a variety of scientific and engineering problems. The specific toolbox of interest for solving optimization problems is called the optimization toolbox [35]. MATLAB optimization toolbox contains a library of programs or m-files, which can be used for the solution of different optimization problems. A commonly implemented function denoted by fmincon is applied to most constrained objective functions. It is the most suitable function to the optimization problem of this investigation, since it uses sequential quadratic programming (SQP) method.

fmincon starts at “x0” and attempts to find a minimizer “x” of the objective function described in the m-file named “fun” subject to the linear inequalities “A×xb” and the linear equalities “Aeq×x=beq.” If there are no linear equalities or inequalities, A, b, Aeq, and beq are replaced with “[].” fmincon can define lower and upper bounds on the design variables in “x” so that the solution is always in the range “lb ≤ x ≤ ub.” fmincon can subject the minimization process to the nonlinear inequalities “c(x)” or equalities “ceq(x)” defined in the m-file constraint function named “nonlcon.” fmincon optimizes such that “c(x) ≤ 0” and “ceq(x) = 0.” fmincon minimizes with the optimization options involved in the structure “options.” Table B.1 defines all input and output arguments related to fmincon optimization function.

ArgumentDescription
funThe function to be minimized
x0, lb., and ubStarting, lower, and upper boundaries of design variables
nonlconThe function that computes the nonlinear inequality constraints “c(x) ≤ 0” and the nonlinear equality constraints “ceq(x) = 0”
fvalValue of the function “fun” at x
exitflagInteger identifying the reasons causing the algorithm terminated
gradGradient at x
hessianHessian at x
lambdaStructure containing the Lagrangian multipliers at the solution x
outputStructure containing information about the optimization process such as number of function evaluations and iterations taken and the used optimization technique

Table B.1.

Arguments related to fmincon function [35].

The distributed load vectors are expressed in the undeformed coordinates (see Figure C.1), as follows:

Figure C.1.

Deformed and undeformed configurations of a blade segment.

Distributed forces:P¯=P¯A+P¯I+P¯G+P¯DEC.1
Distributed moments:q¯=q¯A+q¯I+q¯G+q¯DEC.2

where the subscripts A, I, G, and D refer to the aerodynamic, inertial, gravitational, and damping contributions, respectively. The aerodynamic forcesPA and momentsqA can be obtained using the quasi-steady blade-element strip theory [1, 2].

Figure C.2 shows the velocity triangle and the coordinate system used for computing aerodynamic loads. Vyr and Vzr are the tangential and axial velocity components, respectively. The resultant velocity Vr can be calculated from:

Figure C.2.

Velocity triangle of an arbitrary airfoil section in the deformed state.

VrVyr2+Vzr2=Vyr1+VzrVyr2EC.3

The distributed lift and drag force vectors on an arbitrary airfoil section are given by the aerodynamic formulas (see Figure C-2):

LiftL¯=12ρaVr2CCLsinφiĵ+cosφik̂DragD¯=12ρaVr2CCDcosφiĵ+sinφik̂EC.4

where ρa = air density; C = local chord at spanwise location r; CL = CL (α) = CL αα, lift coefficient; CD = CD (α), drag coefficient; CL α= lift – curve slope; α= angle of attack; φi=inflow angle=θB+α.

Expressed in the undeformed system, the aerodynamic force vector per unit length of the blade is therefore given as:

P¯A=L¯+D¯=PxAî+PyAĵ+PzAk̂EC.5

The distributed inertia loads per unit blade span are obtained by applying D’Alembert’s principle as follows:

Inertial force vector:P¯IB=ρBr¨¯BdydzInertial moment vector:q¯IB=ρBr¯BrE.C×r¨¯BdydzEC.6

where ρB is the blade material mass density, r¨¯Bis the acceleration vector, and r¯Band r¯E.Care the position vectors of an arbitrary point and the elastic center of the blade cross section, respectively.

Next, considering the equilibrium of a differential element of the deformed blade, the equilibrium equations expressed in the rotating (xyz) axes are:

Force equilibrium:F¯x+P¯=0¯Moment equilibrium:M¯x+î×F¯+q¯=0¯EC.7
Internal force vector:F¯=Fxî+Fyĵ+Fzk̂Internal moment vector:M¯=Mxî+Myĵ+Mzk̂EC.8

A wind turbine blade cross section is composed of thin-walled, closed, single or multicellular composite beams. The periphery of each beam cross section is assumed to be constructed of flat composite laminates as shown in Figure C.5. The stress-strain relations of the composite laminates which discretize each cross section are computed using the classical laminate theory CLT. Although each laminate is actually an assembly of multiple layers with different constitutive properties, CLT is used to calculate a set of effective stiffness coefficients that allows a composite laminate to be treated as a single structural element [12]. Therefore the blade cross section may be considered to be composed of discrete sections of homogenous materials.

Generally, for a heterogeneous composite section shown in Figure C.3 and C.4, the modulus-weighted cross-sectional properties are defined as [12]:

Figure C.3.

Blade cross section with discretized composite laminated plates.

Figure C.4.

Cross section for a heterogeneous composite beam.

Ac=1Erefi=1nEiAiyc=1ErefAci=1nEiAiyizc=1ErefAci=1nEiAiziIy=1Erefi=1nEiIvo,i+Aizi2Iz=1Erefi=1nEiIwo,i+Aiyi2Iyz=1Erefi=1nEiIvwo,i+AiyiziEC.9

where Erefis a reference modulus of elasticity, (yi,zi) denotes the geometric centroid of each discrete element of the cross section, and (vo,i , wo,i ) denotes the principal axes of each discrete element.

The parallel axes theorem can be applied to compute the second moments of area about the cross-sectional principal axes as follows:

Iv=IyAczc2Iw=IzAcyc2Ivw=IyzAcyczcEC.10

Once the global cross-sectional properties are computed using the method of modulus-weighted properties, the effective axial stress applied to the blade cross section can be given by:

σxyz=FxAcMzIv+MyIvwIvIwIvw2yyc+MyIw+MzIvwIvIwIvw2zzcEC.11

Finally, by converting the distribution of the effective beam stresses into equivalent in-plane distributed loads on the flat laminates composing the cross section, as shown in Figure C.5, the lamina-level strains and stresses can be computed using the CLT.

Figure C.5.

Resultant loads and moments applying to a composite laminate.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Gerges Edwar Mehanny Beshay and Karam Yousef Maalawi (April 15th 2020). Structural Optimization of Wind Turbine Blades for Improved Dynamic Performance, Design Optimization of Wind Energy Conversion Systems with Applications, Karam Y. Maalawi, IntechOpen, DOI: 10.5772/intechopen.91643. Available from:

chapter statistics

45total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Design Optimization of Wind Energy Conversion Systems with Applications

Edited by Karam Maalawi

Next chapter

Aerodynamic, Structural and Aeroelastic Design of Wind Turbine Blades

By Mohamed Abdou Mahran Kasem

Related Book

First chapter

Introductory Chapter: An Introduction to the Optimization of Composite Structures

By Karam Maalawi

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us