Open access peer-reviewed chapter

# Wind Energy Conversion System Modeling toward Different Approaches

By Ali El Yaakoubi, Lahcen Amhaimar and Adel Asselman

Submitted: May 28th 2019Reviewed: August 8th 2019Published: April 15th 2020

DOI: 10.5772/intechopen.89066

## Abstract

The main focus of this chapter is to modeling the different parts of the wind energy conversion system (WECS) and reviewing the different approaches used in this context. The chapter starts with the aerodynamic and the structural modeling of the wind turbines (WTs), and a description of the steps used to derive a linear time invariant (LTI) model. Thereafter, the chapter introduces models of the electrical actuators in the three phases (abc) and park phases (dq) reference frames, and recalls the assumptions considered. The chapter finishes by presenting the pulse width modulation (PWM) control strategy, the power converters and the pitch actuator models.

### Keywords

• wind turbines
• aerodynamics modeling
• structural modeling
• linear time invariant model
• wind turbines control

## 1. Introduction

Wind, which is the main driving source of the wind turbines (WTs), is characterized by speed and direction. Its origin is movements of air masses in the atmosphere due to temperature or pressure differences. In the lower layer of the atmosphere, the winds are delayed by frictional forces and obstacles that change the wind speed and direction. This is the origin of the turbulence flows. Turbulence increases with higher ground roughness and decreases with distance from the ground, which causes the wind speed to vary over a wide range of amplitudes and frequencies. The major characterization of the winds in the lower layer is that the kinetic energy is distributed in the frequency domain, which is known as the Van der Hoven spectrum [1]. According to this model, the kinetic energy is concentrated around two frequencies, which allows splitting the wind speed vin two components: a slowly variable deterministic component vm, and a rapid variable stochastic component vtas expressed by Eq.(1) [2]:

v=vm+vtE1

In the turbulence domain; the turbulence spectrum observed by a rotating blade element differs from that corresponding to at a fixed point, where a part of the kinetic energy moves toward higher frequencies and will concentrate around integers of rotation frequency [3]. In the stationary domain, the wind speed changes over the rotor disc due to wind shear and tower shadow effects. Wind shear is the changein the mean wind speed with the height due to the skin friction of the earth. Tower shadow effect is the reduction of the wind speed when the wind is felt by blades when passing the tower [4].

For an expectation of the realistic WT behavior as response to its control inputs, it is required to have an accurate model of the wind field that takes into consideration all of the above mentioned characteristics. The knowledge of the mean wind speed in a site, which represents the hour to annually changes, is necessary to predict if a wind energy installation is economically profitable. The turbulence component has low impacts on the annual energy produced, while it has significant impacts on the WT dynamics and thus on the power quality produced. Generally, two wind speed profiles including the two components are often employed in literature, an effective wind speed (i.e., a single point wind speed), and a fully field wind speed (i.e., multi points wind speed over the rotor area). The effective wind speed design is based on a mathematical model, while the fully field wind speed is derived from specialized software tools such as TurbSim [5]. The effective wind speed is sufficient in the view of maximum power capture point, while the full field wind profile is required to simulate the effects of the aerodynamic conditions on the WT structures.

## 2. WECS modeling

WECS includes various multidisciplinary subsystems, which can be classified as aerodynamic, structural and electrical. The aerodynamic subsystem represents the aerodynamic model of the WT. The structural subsystems include blades, tower and drive train models. The electrical subsystems include the generator, the back-to-back converter and the system control models. In order to model the dynamic behavior of the overall system, models of the different components need to be derived and coupled consecutively.

### 2.1 Aerodynamic modeling

The WT aerodynamic subsystem is often modeled by using blade element momentum (BEM) theory [6]. BEM is a quasi-steadymethod, i.e., the assumed local flow conditions and the resulting aerodynamic forces depend only on the current conditions, and would change immediately with any changes in the wind field or blade movement [4]. Under the wind effects, the WT is subjected to surface forces which vary over the rotor area. By applying this theory, the aerodynamic forces acting on each element of the blades can be calculated based on the section of that element. Drag and lift forces, as mentioned above, are the main aerodynamic forces acting on the blades structures, where the lift and the drag are the perpendicular and parallel forces to the incoming flow direction respectively, as illustrated in Figure 1. The lift FLand the drag FDforces acting on each blade section based on the local resultant air velocityv, and the lift cland drag cdcoefficients are given by Eqs. (2) and (3), respectively:

FL=12clαReρbcv2E2
FD=12cdαReρbcv2E3

where α, Re, b and c are local angle of attack, Reynolds number, Blade element length and chord respectively. The lift and drag coefficients are function of the local angle of attack and the Reynolds number. The angle of attack is function of the free stream wind speed, the rotational speed of the WT, and the local induction factor. The induction factor is the factor by which the free stream wind speed has been slowed down when it reaches the rotor plane. These coefficients are determined either by using wind tunnel tests or computational fluid dynamic simulations [7].

Due to the quadratic influence of the wind speed and the dependence of the lift and drag coefficients, the aerodynamic model of the WT is very complex and highly nonlinear. The lift and drag forces could be used to calculate aerodynamic torque and thrust force on the rotor, by integration from blade root to tip [7].

The aerodynamic torque and the thrust force are given by Eqs. (4) and (5), respectively.

Ta=12ρπR3Cqλβv2E4
F=12πR2Ctλβv2E5

where Cqand Ctare torque and thrust coefficients. λandβare tip speed ratio and pitch angle respectively. Cqand Ctcan be calculated by the BEM theory, approximated by a mathematical model for simplicity and computational cost or provided by the WT manufacturer in the form of a look-up table [8]. The power coefficient is a nonlinear function of λand β. Its value depends on the aerodynamic parameters of the turbine as well as on the metrological characteristic of the site. Theoretically this coefficient may take the value of 0.59, which is called the Betz limit, but in reality it is limited between 0.4 and 0.5 [9, 10].

For a detailed modeling, often the aerodynamic model is estimated by specialized aero-elastic simulation tools such as AeroDyn, based on the BEM theory. The outputs from the aerodynamic model are the aerodynamic thrust force and torque, and the distributed forces along each blade. These forces are considered as external forces acting on the structural components of the turbine.

### 2.2 Structural modeling

The structural model of the WT is a combination of flexible and rigid bodies. It includes tower, blades, drive train, nacelle, main bearing and hub. The choice of the modeling type depends strongly to the dynamics required to capture by the model. Generally, when interesting to the WT structural dynamics often tower, blades and drive train are modeled as flexible bodies whereas the others are modeled as rigid. In contrast, when interesting to the electrical dynamics the structural dynamics can be simplified and lumpedinto one rigid mass.

For computational cost simplicity and modal reduction, the flexible bodies motions are assumed to be a superposition of the first eigenmodes. Usually the first and second modes are considered in simulation. This is justified by the fact, that high modes being difficult to detect in actual measurements and having negligible energy content. Therefore, each blade will be modeled by two degree of freedoms(DOFs) for flap-wise and one for edge-wise deflection.

#### 2.2.2 Drive train

The drive train allows to convert the high torque with low speed Tlssto a low torque with a high speed Thss, by a set of differential equations. The dominant inertias in the WT drive train are: the rotor and generator inertias. The drive train is modeled in this chapterby two models: two mass and one mass. In the two mass model, the rotor and generator inertias are separated by a spring, defined by stiffness and damping coefficients, as illustrated by Figure 2. In the one mass model, the rotor and the generator inertias are lumped into one mass.

#### 2.2.2.1 Two mass model

In this model, the inertias of the low speed shaft (LSS), the gearbox and the high speed shaft (HSS) are neglected compared to the inertia of the generator and the rotor. The dominant masses in this model are connected by a flexible shaft and characterized by an equivalent torsional stiffness Klssand damping factor Dlss, which would be referred either on the LSS or the HSS. The model of the drive train in this chapteris referred to the LSS, which is described by the following differential equations:

Jrω̇m=TaTlssDmωmJgω̇g=TlssngThssDgωgTlss=Klssθmθgng+DlssωmωgngE6

where ωmand ωgare the rotor and the generator speed, θmand θgare the angular position of the LSS and the HSS side, Jris the rotor inertia, Jgis the generator inertia referred to the LSS. ngis the gearbox ratio, which is considered ideal and Thssis the high shaft speed torque. Dmand Dgare the rotor and generator friction coefficient, respectively. Tgis the generator torque.

#### 2.2.2.2 One mass model

As mentioned above, the turbine and generator inertias are lumped into one inertia in the one mass model. The one mass model referred to the LSS, is described by the following differential equation:

Jtω̇m=TaTgDtωmE7

where the inertia Jtand the friction Dtcoefficient of the lumped mass are expressed by:

Jt=Jr+ng2JgDt=Dm+ng2Dg

### 2.3 State space model of the WT

The aero-elastic equation of motion of a WT can be expressed as:

Mqutq¨+fqq̇uudt=0E8

where Mrepresents the mass matrix containing inertia and mass components, and fis the nonlinear forcing function vector that includes the stiffness and damping effects. q, q̇and q¨are the enabled DOFs displacements, velocities and accelerations. uis the control input which is the collective pitch angle, while udis the disturbance input. LTI (linear time invariant) model of the WT in the state space representation can be obtained by linearization of the nonlinear model by using aero-elastic software. Fatigue, aerodynamic, structural and turbulence software (FAST) linearizes numerically the equation of motion (Eq. (8)) by perturbing each variable around its respective operating point (OP) [11]. The linearization in FAST is realized in two steps. Firstly, a steady state OP of the enabled DOFs is calculated. Secondly, a numerical linearization about the resulting steady state OP is performed to form periodic matrices of the linear model. The periodic model of the WT in the state space representation, around an OP defined by the triplet β¯v¯ω¯gwhere β¯, β¯and ω¯gthe optimum value of pitch angle, wind speed and generator speed respectively, is given by:

ẋm=Amθmxm+Bmθmum+Bdmθmudmym=Cmθmxm+Dmθmum+DdmθmudmE9

where Amis the state matrix, Bmis the control matrix, Bdmis the disturbance matrix and Cmis the output matrix. The state vector is given by: xm=ΔqΔq̇T, where Δrepresents the perturbation of a parameter from its optimal value. It should be noted that the matrices of the periodic state space model depend on the rotor azimuth position θm. Therefore, it is important to mention that the dynamics of the model of Eq. (9) are expressed in a mixed reference coordinates, as denotes the subscript m, i.e., rotating and fixed references coordinates. In fact, the dynamics of the blades are expressed with respect to rotating reference coordinate, while those of the nacelle and tower are expressed with respect to a fixed reference coordinate. To derive a LTI model from the periodic linear model, it is required to apply a multi-blade coordinate (MBC) transformation [12]. The transformation provides a unifying reference frame for analysis and control design, where all DOFs are expressed with respect to the same reference coordinate. Employing this transformation, the model of Eq. (9) can be transformed to a periodic model expressed in a fixed reference coordinate system:

ẋFR=AFRθmxFR+BFRθmuFR+BdFRθmudFRyFR=CFRθmxFR+DFRθmuFR+DdFRθmudFRE10

The subscript FRsignifies that the model is expressed in the fixed reference coordinate. The states, control inputs and outputs measurement of the model in the rotating reference coordinate are transformed to those corresponding in the fixed reference coordinate using the matrix transformation (Tcθm), i.e., xm=TcθmxFR, um=TcθmuFRand ym=TcθmyFR. The disturbance input udmis not transformed because it is already expressed in the fixed reference coordinate.

Tcθm=1cosθmsinθm1cosθm+2π3sinθm+2π31cosθm+4π3sinθm+4π3E11

The WT averaged periodic model over the rotor rotational period (θm), which is the LTI model, is expressed as:

ẋ=Ax+Bu+Bdudy=Cx+Du+DdudE12

## 3. Actuators

The actuators aim to control WECS by receiving the desired signals from the control system and providing the real signals. They include generator, converters, pitch and yaw. The generator actuator receives a desired electromagnetic torque from the torque controller and provides a real electromagnetic torque to the drive train. The converters actuator aims to produce a current in the generator that produces a real torque that matches the desired value, by means of its control system. Generally, the electrical dynamics are much faster compared to the WT structural dynamics. Thus, they can be represented by a small communication delay between the torque demanded by the generator and the actual air-gap torque acting on the mechanical system when the structural dynamics are required. Such delay can be represented by a first order system with a small time constant. The pitch actuator receives the desired pitch angle for each blade from the pitch control, and provides the realistic pitch angle for the blade. The yaw mechanism rotates the nacelle around the vertical tower axis to place the turbine directly into the wind for a maximum power generation. Under high-speed winds, the yaw mechanism turns the blades 90° from the direction of the wind to reduce stress on internal components and avoid over-speed conditions. The yaw actuator requires the wind inflow direction and nacelle direction as inputs, and provides the yaw angle.

As mentioned above, the pitch angle is used as control input in the full load region, while the generator torque is used as control input in the partial and the full load regions. The yaw actuator is not activated where the WT is aligned with the wind direction, to allow a maximum power capture from the wind.

### 3.1 Generator modeling

The generator chosen in this manuscript is squirrel cageinduction generator (SCIG) because it is cheap, simple in construction, easy for maintenance, and easily replaceable. SCIG consists of a set of windings on the stator and rotor. The stator and rotor have three windings corresponding to the three phase system, which can be coupled either in star or delta. The stator phases of the machine are connected to the power grid through a back-to-back converter, whereas those of the rotor are short-circuited by a conductive ring at each winding side, as illustrates Figure 3.

#### 3.1.1 Energy generation theory

According to Ferraris theorem, a rotating magnetic field in the gap of the machine having a sliding speed with respect to stator (Ω=ωsp) is created [13]. The rotor turns, with respect to the stator, at the mechanical speed of ωr=dt,where θis the angle between the stator and the rotor windings. The relative speed between the rotating field and the rotor induces a variable voltage eindin each closed loop of the rotor conductors, where its magnitude linked to the stator flux Φsby Faraday’s law [5]:

dΦsdt=eindRciindE13

with iindis the induced current in the conductor, and Rcis the resistance of the conductor. The stator flux linkage is defined as the flux linkage of a single turn ϕmultiplied by the number of turns Nin the coils:

Φs=

#### 3.1.2 SCIG modeling assumptions

Thereafter some assumptions related to modeling of the SCIG are giving in the following of this manuscript [14]:

• The three phasesystem is symmetric and balanced with constant gap.

• The iron saturation is discarded.

• The magnetic flux density is assumed radial to the gap.

• All kind of losses in the iron are neglected.

• Both windings of the machine are sinusoidal distributed, whose axes are displaced by 120°, which produce a sinusoidal magnetic field distribution in the gap.

• The coefficients of the inductances are fixed and the coefficients of the mutual inductances depend to the windings position.

#### 3.1.3 abc and dq0 reference frames

The quadratic-direct-zero (dq0) (see Figure 4), proposed by Park, is a reference rotating frame, for modeling and analysis of the electrical systems [15]. It allows to remove the dependency of certain parameters like inductances, and to simplify the modeling process. Generally, the electrical systems are modeled in phases reference frame. For control and analysis purposes, it is required to have the dq0 model of the system. The passage from a reference frame to the other is performed by means of the Park transformation matrix. By application of this transformation, the three ac (a, b and c) quantities of the generator can be reduced to two DC (d and q) components and a homopolar component.

Assuming that the dqframe is oriented by an angle θfrom the abcreference frame, the Park transformation matrix to convert the abcquantities of the SCIG to the dq0 quantities is given by:

Tθ=23cosθcosθ2π3cosθ+2π3sinθsinθ2π3sinθ+2π3121212E14

Therefore, each quantity can be transformed from the abcto the dq0 reference frame by the following expression:

χdq0=TθχabcE15

where χdq0and χabcare the representation of a machine parameter in the dq0 and abcreference frames, respectively. χmay represent voltage, current or flux of the machine.

The inverse transformation of each quantity can be recovered by means of this expression:

χabc=T1θχdq0E16

where

T1θ=cosθsinθ1cosθ2π3sinθ2π31cosθ+2π3sinθ+2π31E17

#### 3.1.4 SCIG model in the abc reference frame

Applying Faraday law’s, under the above mentioned hypotheses and by using the induction motor’s illustration, the dynamic model of an induction machine in the abcreference frame is described by the following differential equation [14, 16, 17]:

vsabcvrabc=rsabc00rrabcisabcirabc+ddtΦsabcΦrabcE18

The voltage, current and flux in the abcreference frame are expressed as:

vsabc=vsavsbvsc,vrabc=vravrbvrc,isabc=isaisbisc,irabc=irairbirc,Φsabc=ΦsaΦsbΦscandΦrabc=ΦraΦrbΦrc

where Φsabcand Φrabcare the stator and rotor linkage fluxes in the abcreference frame, which are given by:

ΦsabcΦrabc=LssabcLsrabcLrsabcLrrabcisabcirabcE19

By which rsabc, Lssabc, rrabcand Lrrabcare the resistances and inductances of the rotor and stator in the abcreference frame, which are expressed as:

rsabc=rs000rs000rs,rrabc=rr000rr000rr,Lssabc=Ls+LmLm2Lm2Lm2Ls+LmLm2Lm2Lm2Ls+LmandLrrabc=Lr+LmLm2Lm2Lm2Lr+LmLm2Lm2Lm2Lr+Lm.

Lsand Lrare the inductance value of stator and rotor windings in each phase, respectively. Lmcan represent the stator magnetizing inductance, the stator-rotor mutual inductance, and the rotor magnetizing inductance all referred to the stator side.

Lsrabcand Lrsabcare the mutual inductances between the stator and rotor phases for the rotor and stator respectively, which are expressed as:

Lsrabc=Lmcosθrcosθr+2π3cosθr2π3cosθr2π3cosθrcosθr+2π3cosθr+2π3cosθr2π3cosθr,Lrsabc=Lmcosθrcosθr2π3cosθr+2π3cosθr+2π3cosθrcosθr2π3cosθr2π3cosθr+2π3cosθr.

rsand rrare resistance value of stator and rotor windings in each phase, respectively.

The induction machine coupling inductance matrix is defined by:

Lθr=LssabcLsrabcLrrabcLrsabc

The rotor quantities referred to the stator are given by: vrabc=Ns/Nrvrabc, irabc=Ns/Nrirabc, rrabc=Ns/Nr2rrabc, Lsrabc=Ns/Nr2Lsrabc, Lrsabc=Ns/Nr2Lrsabcand Lrrabc=Ns/Nr2Lrrabc. Where where Nsand Nrare the number of turns per phase in the stator and rotor windings, respectively.

The angular displacement of the rotor, with respect to the stator, and the stator are given by Eqs. (20) and (21), respectively:

θrt=ωrtdt+θr0E20
θst=ωstdt+θs0E21

where θr0and θs0are the initial position of the rotor and stator at t = 0, respectively. The rotor speed is linked to the mechanical speed of the generator ωgby:

ωr=pωgE22

The generator speed is linked to the turbine rotational speed (ωm)by:

ωg=ngωmE23

#### 3.1.5 SCIG model in dq0 reference frame

The SCIG model in the dq0 reference frame is obtained by applying the transformation matrices Tθsand Tθsθron the abcstator and rotor parameters, respectively. Applying these transformations on the Eqs. (14) and (17) we can write

vsdq0=Tθsrsdq0isdq0T1θs+TθsddtT1θsΦsdq0vrdq0=Tθsθrrrdq0irdq0T1θsθr+TθsθrddtT1θsθrΦrdq0E24

After simplifications, the stator and rotor voltages in the dq0 reference frame are given by:

vsdq0=rsdq0isdq0+ddtΦsdq0+ΩsΦsdq0vrdq0=rrdq0irdq0+ddtΦrdq0+ΩrΦrdq0E25

where the voltage, the current and the flux in the dq0 are expressed as:

vsdq0=vsdvsqvs0,vrdq0=vrdvrqvr0,isdq0=isdisqis0,irdq0=irdirqir0,Φsdq0=ΦsdΦsqΦs0andΦrdq0=ΦrdΦrqΦr0

By which Ωsand Ωrare the stator and rotor speed transformation matrices, which are given by: Ωs=0ωs0ωs00000and Ωr=0ωsωr0ωsωr00000, respectively.

The stator and rotor resistances in the abcand dq0 are equal, i.e., rsdq0=rsabcand rrdq0=rrabc.

Applying the previous transformations (Eqs. (14) and (17)), we obtain:

T1θsΦsdq0=Lssdq0T1θsisdq0+Lsrdq0T1θsθrirdq0T1θsθrΦrdq0=Lrsdq0T1θsθrisdq0+Lrrdq0T1θsθrirdq0E26

After simplifications, the stator and rotor fluxes of the SCIG in the dq0 reference frame are given as:

Φsdq0=Lssdq0isdq0+Lsrdq0irdq0Φrdq0=Lrsdq0isdq0+Lrrdq0irdq0E27
whereLssdq0=Ls+32Lm000Ls+32Lm000Lls,Lrsdq0=Lsrdq0=32Lm00032Lm00032LmandLrrdq0=Lr+32Lm000Lr+32Lm000Lr

Since the rotor windings are short circuited, the abcrotor voltage components are zero, vrabc=0. Moreover, as the electric system is supposed equilibrate, the electrical parameters of the model in the dq0 reference frame are reduced only to two components (d and q), i.e., the homopolar components equal zero (v0s=v0r=i0s=i0r=Φ0s=Φ0r). Therefore, the dynamic model of the SCIG machine in the dqreference frame is expressed as:

vsq=rsisq+dΦsqdt+ωsΦsdvsd=rsisd+dΦsddtωsΦsqvrq=rrirq+dΦrqdt+ωsωrΦrd=0vrd=rrird+dΦrddtωsωrΦrq=0E28

Figure 5 represents the electrical scheme of the SCIG in the dqreference frame based on the electrical model of Eq. (28).

After simplifications, the electromagnetic flux in the dqreference frame of the SCIG is given by the following equation:

Φsq=Lsisq+LMirqΦsd=Lsisd+LMirdΦrq=Lrirq+LMisqΦrd=Lrird+LMisdE29
withLs=Lls+LMLr=Llr+LMLM=32Lm.

#### 3.1.6 State space model of the SCIG machine in the dq0 reference frame

Combining Eqs. (28) and (29), the model of the SCIG in the state space representation can be derived, where xn=isqisdirqirdTis the state vector, un=vsqvsdTis the control input. The state (A) and control input (B) matrices of the model are given by:

A=rsσLsωs+pωgLm2σLsLrLmrsσLsLrpωgLmσLrωs+pωgLm2σLsLrrsσLspωgLmσLrLmLrσLsLrLmLrσLsLrpωgLmσLsrrσLrpωgσωspωgLmσLsLmLrσLsLrωspωgσrrσLr,B=1σLs00LmσLsLr01σLs0LmσLsLr

whereσ=1LM2/LsLr.

#### 3.1.7 Active power and electrical torque in the dq0 reference frame

The electromagnetic torque Tgdeveloped by the machine depends to the instantaneous current circulated in each of the six windings and the angle between the stator and rotor windings θ. Its expression in the abcreference frame is given as [17]:

Tg=p2isabcirabcTLθisabcirabcE30

Its model in the dq0 reference frame is given by [14]:

Tg=32pΦrdirqΦrqirdE31

The active power (Ps)yield by the SCIG is given as the product of the electromagnetic torque (Tg) and the generator speed (ωg):

Ps=TgωgE32

Its expression in the dq0 reference frame, function of voltages and currents stator, is given by:

Ps=32vsdisd+vsqisqE33

### 3.2 Converters modeling

The back to back converter interfaces the WT and the utility grid. It consists of a machine side converter (MSC), which acts as rectifier, and GSC, which acts as inverter, connected by a DC link as illustrated in Figure 6. The DC link allows an optimum injection of the generated energy to the grid, and a frequency decoupling between the machine and the grid. The structure of the converters is bi-directional power flow, i.e., the direction of the power flow can be reversed at any time. The converter outlined in this chapter is a voltage source converter (VSC), which is a switching IGBT-diode device. Its switching model describes the steady state and the dynamic behavior, including slow transients and high frequency components of voltages and currents. The high frequency of the components are not desirable for control and analysis purposes, thus it is common to use the average value of the variables rather than the instantaneous value. Therefore, only the fundamental ac dynamics which are considered for control design and analysis [5].

The main objective associated to the rectifier is to control the WT variable speed operation, i.e., to implement the MPPT (maximum power point tracking) strategy, whereas that associated to the inverter is to manage the power flow exchanged between the machine and the utility grid. As the modeling process of these converters are the same, we will here only to MSC. Figure 6 illustrates the connecting scheme of the MSC to the DC link. As can be shown, each phase of the converter is composed of a set of two IGBT-diode devices connected to the DC circuit. It has been proven that the efficiency of the converter is around 94–98%, thus it is common practice to represent the back to back converters as ideal switches. Therefore, each set of the IGBT-diode of the phase will be represented by a single switch, where its switches states are complementary. Moreover, for analysis purpose, it is assumed that there is a common point in the DC link that splits the voltage Vdcinto two voltages of amplitudes Vdc/2[17].

#### 3.2.1 MSC modeling assumptions

Thereafter are recalled some assumptions related to the converter modeling [14]:

• The switching of the components is instantaneous.

• The voltage drops at their terminals are negligible.

• The dead times are neglected.

• The load is three-phase balanced and coupled in star with isolated neutral point.

• It is assumed that the IGBTs switch at a frequency at least 10 times of the fundamental frequency of the ac grid, thus the action of this commutation can be represented by the average of the duty cycle [5]. Therefore, the currents and voltages of the VSC are expressed as function of the duty cycle of the IGBT.

#### 3.2.2 MSC model in abcreference frame

Applying Kirchhoff’s law to circuit given in Figure 6, the compound voltages between the three phases at the right side, are:

uAB=uAOuBOuBC=uBOuCOuCA=uCOuAOE34

At the load side, the compound voltages are function of the simple voltages, i.e., the voltages between the phases A, Band Cand the point O, are expressed as:

uAB=vAvBuBC=vBvCuCA=vCvAE35

The voltages equations of the VSC can be expressed also as [17]:

uABuCA=2vAvB+vCuBCuAB=2vBvC+vAuCAuBC=2vCvA+vBE36

Moreover, as the voltages vA, vB, vCform an equilibrium three phase system, the expression of Eq. (36) can be reduced to [5]:

uABuCA=3vAuBCuAB=3vBuCAuBC=3vCE37

Thus, the simple voltages function of the voltages between the phases and the effective point O, are expressed by:

vA=13uABuCAvB=13uBCuABvC=13uCAuBCE38
vAvBvC=13211121112uAOuBOuCOE39

As mentioned above, each phase is formed by two complementary ideal switches, i.e., if we associate to each phase a binary value of command Si(i = A, B, C) we have:

1. Si= +1, the upper switch is on and the lower switch is off;

2. Si= −1, the upper switch is off and the lower switch is on.

Therefore, the voltages between the phases and the effective point O, function of the switch functions and the DC link voltage are derived

uAOuBOuCO=Vdc2SASBSCE40

Combining Eqs. (39) and (40), the simple phase voltages, function of the switching functions and the DC link voltage, are obtained:

vAvBvC=13211121112Vdc2SASBSCE41

Finally, the simple voltages of the three phases at the output of the converter are given function of the switch functions. The switching functions are the commands of the converter, where their states will be provided from PWM control strategy.

Moreover, applying Kirchhoff’s law at the load side, the average model of the VSC, in the abcreference frame is given by:

where vSA, vSBand vSCare the induced voltages in the stator phases of the generator.

#### 3.2.3 MSC model in the dq0 reference frame

Expressing the voltages and currents of the VSC in the dqreference frame allows controlling independently the electromagnetic torque and rotor flux.

The model of the MSC in the dq0 reference frame can be obtained by application of the transformation matrix of Eq. (14) on Eq. (42) [8]:

Lsdisddt+rsisd=vsdSdVdc2Lsdisqdt+rsisq=vsqSqVdc2E43

where Sdand Sqare the dand qswitch functions, respectively.

This model will be used to control the SCIG machine by controlling the VSC based on PWM control strategy. Selective harmonic elimination, sinusoidal pulse width modulation (SPWM) and space vector modulation are the commonly found PWM strategies in literature. SPWM is the most used due to its simplicity and effectiveness [5].

#### 3.2.4 SPWM control strategy

SPWM is a multi-pulses-based modulation method that varies the pulse width of the converter output voltage in a sinusoidal manner, following a target reference voltage by comparing a low frequency reference signal fmwith a high frequency carrier signal fp, where fmfp[5]. The frequency of the reference signal should equal the frequency of the output voltage. For each phase leg, there are two complementary standing signals, one for the lower switch and the other for the upper switch. In fact, when the reference signal exceeds the carrier signal, the upper switch is turned on and the lower switch is off. Otherwise if the reference signal is less that the carrier signal, the lower switch is turned on and the upper switch is turned off.

For a three phase VSC, the three phase reference signals vmA, vmB, and vmCand a high frequency carrier (vp) are used

vmA=V̂msin2πfmtvmB=V̂msin2πfmt2π3vmC=V̂msin2πfmt+2π3E44
vp=V̂p1+4tTpV̂p34tTpift0Tp2tTp20

where V̂mand V̂pare the voltages amplitude, and Tpis the period of the carrier.

The states of the IGBT (Sj) (j = 1,2,3) are derived as follow:

Sj=+11ifvmivp0vmivp<0

where the states of the lower switches will be complementary with those of the upper switches. Thus the simple voltages (vi)at the terminal of the converter are:

vi=+Vdc2Vdc2ifSj=+1Sj=1

Therefore, the voltages at the output of the converter are a set of rectangular slots, whose width varies as function of the control law of the IGBT, the duration of switch and the time of functioning. Two parameters characterize this control strategy: the modulation index, which is the ratio of fmand fp; and the tuning coefficient, which is the ratio of V̂mand V̂p. For a proper control, these coefficients must be tuned carefully. Various works are interested to WTs control either for maximum capture, power flow exchange as well as power regulation and structural loads minimization. Readers can refer to [18, 19, 20, 21, 22, 23].

### 3.3 Pitch actuator

The pitch actuator system aims to drive the blades to the desired positions received from the pitch control. Due to the large inertia of the blades, there is a significant delay between a blade pitch command and the blade actually arriving in the desired position. The actuator dynamics are largely dominated by the blades and the motor inertias, and the capabilities of the pitch system; while the torque caused by the aerodynamic loading of the blades has only small impacts. Often the pitch actuator model is linear first order system with low pass filter, additional rate limiter acting on the demanded pitch angle βrefand angle limiter to set the maximum and the minimum pitch angles. The pitch system can be either hydraulic, electromechanical, or a hybrid of the both. Eq. (45) represents the model of the actuator, where τβis the time constant of the model.

ββref=1τβ.s+1E45

Figure 7 illustrates the pitch actuator scheme based on the mathematical model of Eq. (45).

## 4. Conclusions

This chapter presented the different models of the WECS such as wind turbines aerodynamics and structural dynamics, generator and converter actuators. The chapter has outlined the hypotheses and the approaches used in control and modeling such as state space representation, and the abcand the dqreference frames. The chapter finished by presentation of the SPWM control strategy and the pitch angle actuator.

## Conflict of interest

The authors do not have any conflicts of interest to declare.

## Nomenclature

abc

three phases system

dq

park phases

θs

position angle between the abcand the dqreference frames

θr

position angle between the rotor and the stator windings in the abcreference frame

ar, br, cr

rotor windings axes in the abcreference frame

as, bs, cs

stator windings axes in the abcreference frame

Ps

electrical power

Tg

electromagnetic torque

ωg

generator speed

is

stator current

Pg

electrical energy injected to grid

Lg

filter inductance

Rs

filter resistances

ig

grid current

Vdc

MSC

machine side converter

GSC

grid side converter

VSC

voltage source converter

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Ali El Yaakoubi, Lahcen Amhaimar and Adel Asselman (April 15th 2020). Wind Energy Conversion System Modeling toward Different Approaches, Design Optimization of Wind Energy Conversion Systems with Applications, Karam Y. Maalawi, IntechOpen, DOI: 10.5772/intechopen.89066. Available from:

### Related Content

#### Design Optimization of Wind Energy Conversion Systems with Applications

Edited by Karam Maalawi

Next chapter

#### Robust Nonlinear Control Strategy for Small Wind Turbines: A Case Study

By Ridha Cheikh and Hocine Belmili

#### Optimum Composite Structures

Edited by Karam Maalawi

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.