Modeling and Simulation of the Variable Speed Wind Turbine Based on a Doubly Fed Induction Generator

Abstract

This chapter presents the modeling and simulation results of variable speed wind turbine driven by doubly fed induction generator (DFIG). The feeding of the generator is ensured through its stator directly connected to the electrical grid and by its rotor connected to the grid through two power converters, which are controlled by the pulse width modulation (PWM) technique. This configuration is the most used in the wind power generation systems. For the variable speed operation of the studied system, the maximum power point tracking strategy is applied for the turbine, and the stator flux-oriented vector control is used for the generator. The MATLAB/Simulink software is used for the system modeling and simulation. For the wind velocity model, a random wind profile is simulated, and the turbine and the generator parameters are extracted from an existing wind turbine system in the literature. The obtained results are addressed in this chapter.

1. Introduction

With the global warming issues and the climate changes, there is a serious need for the use of the renewable energy resources in the electricity generation industry. Currently, the wind represents one of the most important renewable energy resources, used for generating electrical energy in the world. In terms of the total installed wind capacity, it becomes up to 539 GW across the globe in 2017 [1]. The rapid rate of the wind energy industry growth is caused by the cost-effectiveness of electricity production from wind farms, compared to electricity production cost from fossil fuel energy [2], the stability of electricity cost [3], the short commissioning time of wind farms [4], and the ingenuity of skillful engineers.

According to a wind market survey, the doubly fed induction generator (DFIG) is the most popular generator used in the speed variable wind turbines (SVWT) [5]. It is a wound rotor asynchronous machine which has the stator windings directly connected to the electrical grid, and its rotor is linked to the constant frequency grid by means of two bidirectional power converters.

This configuration, known as “Scherbius structure” and shown in Figure 1, has several advantages of controlling independently active and reactive power [6, 7]. Moreover, the power converters used are sized to transfer only a fraction, equal to at most 30% of the turbine rated power [8, 9], which results in small-size, low-cost, less acoustic noise and reduced loss rate in the power converters [10]. Moreover, The DFIG-based wind turbine allows the rotor speed to be varied with the wind speed, and the speed variation range is around ±30% around the synchronism speed [11]. As a result, the wind generation system could operate in hyposynchronous and hypersynchronous mode, which would extract the maximum aerodynamic power for each wind speed value.

Figure 1.

Architecture Scherbius of DFIG-based wind turbine.

In order to design fault diagnosis and control approaches based on models for wind turbines, the development of a mathematical model, which represents as much details as technically possible and gives an accurate idea of the dynamic behavior of the system, seems to be an important step. For several purposes, different wind turbine models have been developed. In literature, we find for wind systems the aerodynamic model, which aims to verify and optimize the blade design, depending on predefined specifications, while the mechanical model is used by engineers for establishing a safe and economical dimensioning of the whole wind turbine system. Moreover, the economic model is used in the case of manufacturing and installing wind turbines with the purpose to evaluate the cost-effectiveness. In addition, there are models which predict the weather conditions and the power output of wind farms. Furthermore, there are models, which have the objective of evaluating the impact of wind turbines on the environment such as the evaluation of noise produced by the wind turbine operation. Finally, the general-purpose models concerned with the electrical properties of wind turbines are widely used.

The speed variable wind turbine (SVWT) model, developed and simulated in this work, is concerned with providing time simulation signals that can be exploited for designing fault diagnosis approaches based on models; the software tool used for simulation is the MATLAB/Simulink environment. In [12], a wind turbine model of a fixed-speed, stall-regulated system has been developed with the aim of measuring and evaluating the power quality impact of wind turbines on the grid. In addition, a model for a wind turbine generation system based on a DFIG, including the mechanical dynamics, the wind turbine electrical system, the converter, and the electrical grid has been presented in [13]. Luis et al. [14] presented the most commonly used wind turbine model meeting objectives as production energy, safety of turbine, grid connection, and others. Moreover, in [15], the detailed mechanical structural modeling of the wind turbine connected to the grid and based on DFIG has been developed, and it has been validated using NREL’s simulation tool, FAST v7; for analysis of the dynamic behavior of the wind power plant with DFIG under the grid fault conditions, modeling of the whole system has been established in [16]. Furthermore, with the use of an electromagnetic transient simulation software, the wind turbine driven by DFIG model is elaborated in [17].

The organization of this chapter is as follows: the wind turbine structure is described in Section 2; the wind modeling turbine is presented in Section 3. In Section 4, the simulation results in MATLAB/Simulink environment are shown. Finally, the conclusion is presented in Section 5.

2. The wind turbine structure

The wind turbine is a complex system containing different components, which involved different domains: electrical, mechanical, and electronic areas and others. The complete model of the studied wind turbine system is represented by a set of blocks each representing a functional entity of the system. The general structure of a wind system is given by Figure 2.

Figure 2.

Wind turbine system structure.

As it is shown in the Figure 2, the wind velocity and the fixed-pitch angle represent the input of the wind turbine system. The aerodynamic conversion entity is formed of three blades which capture and convert the kinetic energy of the wind into mechanical energy, recovered on the slow rotating shaft. Then, the gearbox device increases the turbine low speed and makes it suitable with the generator rotational speed, which is about 1500 rpm. The generator receives mechanical energy and transforms it into electrical energy. The two power converters used are insulated gate bipolar transistor (IGBT) type and controlled by pulse width modulation (PWM) technique; they allow the independent control of the active and reactive powers and also the transfer of the slip power in two directions: from the generator to the network and from the network to the generator.

3. The wind turbine modeling

In this section, the mathematical model of each wind turbine block is presented.

3.1 Wind speed model

The wind resource is an important element in a wind energy system, and it represents a determining factor in the calculation of electricity production because, under optimal conditions, the power captured by the wind turbine is a cubic function of the wind speed. The wind is a moving air mass, and the wind kinetic energy is given by:

$$E=\frac{1}{2}·m·v^2$$ E1

where m is the moving air mass [g] and v is the air moving speed [m/s].

The wind power during is expressed as:

$$P_{\mathit{\text{wind}}}=\frac{E}{\mathit{\Delta t}}$$ E2

The wind speed v is generally represented by a scalar function evolving over time, given by V = f(t). It can also be divided into two components: a slowly varying part denoted as V₀ and a random varying part denoted as V_t; it represents the wind fluctuations. Therefore, the wind velocity can be written as:

$$V\left(t\right)=V_0+V_t\left(t\right)$$ E3

To mathematically model the wind speed profile, the literature offers three techniques:

• The first method is white noise filtering technique, in which the turbulence impact is corrected by the use of a low-pass filter having the following transfer function [18]:

$$F\left(s\right)=\frac{1}{1+\tau ·s}$$ E4

where τ is the filter time constant. It depends on the rotor diameter and the wind turbulence intensity and the average wind speed. Figure 3 shows the method of reconstruction of the wind profile using this technique.

• The second method of generating the wind speed profile is that which describes wind variations using the spectral density established by meteorologist I. Van der Hoven. In this model, the turbulence part is considered as a stationary random process, and therefore it does not depend on the variation of the mean wind speed [19, 20]. The variation of the wind speed v(t) is thus written in the form of the harmonic sum:

Figure 3.

Wind profile construction scheme by white noise filtering.

$$V_v\left(t\right)=A+\sum _{i=1}^n{a}_k·sin\left({\omega }_k·t\right)$$ E5

where A is the wind speed average value; a_k is the amplitude of k-order harmonic; ω _k is the pulsation of k-order harmonic; and i is the last harmonic rank retained in the wind profile calculation.

• The third method is the Weibull distribution in which a given site wind potential is obtained by measuring the average wind speed in regular time intervals. Then, the data obtained are then divided into numbers by wind speed classes using histogram [21]. The wind profile over a desired time period, respecting the Weibull distribution, is given by:

$$V_v\left(t\right)=\left(1+{\xi }_v\left(t\right)-{\xi }_v\right)·V_v$$ E6

where V_v is the wind speed average value and ξ_v is the disturbance mean value expressed by:

$${\xi }_v\left(t\right)={\left(-\frac{\mathit{ln}\left(\mathit{rand}\left(t\right)\right)}{c_v}\right)}^{\frac{1}{k_v}}$$ E7

where rand(t) is a function generating, in a uniform distribution, random numbers between 0 and 1 and (C_v, k_v) is a parameter pair, determined by analysis of the wind class histogram. C_v is a scale factor generally greater than 5. The shape factor k_v is greater than 3 if the histogram shape is like that of a normal distribution, characterized by an uniform distribution around a mean value.

In this work, we adopted the second method to generate the random profile of the wind speed applied in the studied wind system input.

3.2 Aerodynamic conversion model

The aerodynamic conversion system is the wind turbine part, which is facing the wind; it generally comprises three blades of length R. Three-bladed wind turbines are much more common than two-bladed wind turbines. The turbine captures the kinetic energy of the wind and transforms it into mechanical energy recovered on the slow rotating shaft.

The kinetic power of the wind is given by:

$$P_{\mathit{\text{wind}}}=\frac{\rho ·S·{v_{\mathit{\text{wind}}}}^3}{2}$$ E8

The aerodynamic power is expressed as follows:

$$P_{\mathit{\text{aero}}}=C_p\left(\lambda ,\beta \right)·P_{\mathit{\text{wind}}}=C_p\left(\lambda ,\beta \right)·\frac{\rho ·S·{v_{\mathit{\text{wind}}}}^3}{2}$$ E9

The aerodynamic torque T_aer is given by the following expression:

$$T_{\mathit{aer}}=\frac{P_{\mathit{\text{aero}}}}{{\mathrm{\Omega }}_t}=\frac{1}{2·{\mathrm{\Omega }}_t}·C_p\left(\lambda ,\beta \right)·\rho ·S·{v_{\mathit{\text{wind}}}}^3$$ E10

where Ωt is the turbine speed [rad/s], ρ is the air density, ρ = 1.225 kg/m³, S = π R₂ is the rotor surface [m²], R is the blade length [m], and v_wind is the wind speed upstream of the wind turbine rotor [m/s]. λ is the speed ratio. It is a unitless parameter, related to the design of each wind turbine, and it represents the ratio between the speed of the blade’s end and that of the wind at the rotor axis or also called hub. λ is expressed as follows:

$$\lambda =\frac{{\mathrm{\Omega }}_t·R}{v}$$ E11

This parameter depends on the blade number of the wind turbine. If the blade number is reduced, the rotor speed is high, and a maximum of power is extracted from the wind. In the case of multiblade wind turbines (Western Wind Turbines), the speed ratio is equal to 1; for wind turbines with a single blade, λ is about 11. The three-bladed wind turbines, as in our study, have a speed ratio of 6 to 7. The speed ratio of Savonius wind turbines is less than 1 [22].

Cp is the power coefficient or aerodynamic transfer efficiency that varies with the wind speed. This coefficient has no unit, and it depends mainly on the blade aerodynamics, the speed ratio λ, and the blade orientation angle β. Betz has determined a theoretical maximum limit of the power coefficient C_pmax = 16/27 ∽ 0.59. Taking into account losses, wind turbines never operate at this maximum limit, and the best-performing wind turbines have a Cp between 0.35 and 0.45. Cp is specific to each wind turbine, and its expression is given by the wind turbine manufacturer or using nonlinear formulas. To calculate the coefficient Cp, different numerical approximations have been proposed in the literature. The Cp expressions frequently encountered in the literature are presented in Table 1.

Power coefficient type, Cp	Formula
Exponential	$0.22·\left[\frac{116}{{\lambda }_i}-0.4·\beta -5\right]·e^{\frac{12.5}{{\lambda }_i}}$ avec $\frac{1}{{\lambda }_i}=\frac{1}{\lambda +0.08·\beta }-\frac{0.035}{{\beta }^3+1}$ [23]
	$0.5·\left[\frac{116}{{\mathrm{\lambda }}_{\mathrm{i}}}-0,4·\mathrm{\beta }-5\right]·{\mathrm{e}}^{-\frac{21}{{\mathrm{\lambda }}_{\mathrm{i}}}}+0.068·\mathrm{\lambda }$ [8, 24, 25]; Avec: $\frac{1}{{\lambda }_i}=\frac{1}{\lambda +0.08·\beta }-\frac{0.035}{{\beta }^3+1}$
	$0.5176·\left[\frac{116}{{\lambda }_i}-0.4·\beta -5\right]·e^{-\frac{21}{{\lambda }_i}}+0.0068·\lambda$ [26]; Avec: $\frac{1}{{\lambda }_i}=\frac{1}{\lambda +0.08·\beta }-\frac{0.035}{{\beta }^3+1}$
	$0.5109·\left[\frac{116}{{\lambda }_i}-0.4·\beta -5\right]·e^{-\frac{21}{{\lambda }_i}}+0.0068·\lambda$ [27]; Avec: $\frac{1}{{\lambda }_i}=\frac{1}{\lambda +0.08·\beta }-\frac{0.035}{{\beta }^3+1}$
	$0.44·\left[\frac{125}{{\lambda }_i}-6.94\right]·e^{-\frac{16.5}{{\lambda }_i}}$ [21]; Avec: ${\lambda }_i=\frac{1}{\frac{1}{\lambda }+\mathrm{0,002}}$
Sinusoidal	$0.5-0.167·\left(\beta -2\right)·sin\left[\frac{\pi ·\left(\lambda -3\right)}{18.9-0.3·\left(\beta -2\right)}\right]-0.00184·\left(\lambda -3\right)·\left(\beta -2\right)$ [28]
	$0.3-0.00167·\left(\beta -2\right)·\mathit{sin}\left[\frac{\pi ·\left(\lambda +0,1\right)}{10-0.3·\left(\beta -2\right)}\right]-0.00184·\left(\lambda -3\right)·\beta$ [29]
Polynomial	$6·{10}^{-7}·{\lambda }^5+{10}^{-5}·{\lambda }^4-65·{10}^{-5}·{\lambda }^3+2·{10}^{-5}·{\lambda }^2+76·{10}^{-3}·\lambda +0.007$ [30]
	$79.5633.{10}^{-5}·{\lambda }^{-5}-17.375·{10}^{-4}·{\lambda }^4+9.86·{10}^{-3}·{\lambda }^3-9.4·{10}^{-3}·{\lambda }^2+6.38·{10}^{-2}·\lambda +0.001$ [31]

Table 1.

Different numerical formulas of the power coefficient Cp.

Since we had as an objective the modeling and simulation of a three-bladed wind turbine with a nominal power of 3 kW; the parameters of both: the wind turbine and the generator have been used from [30]. For this reason, the analytical expression of the power coefficient C_p is given by:

$${\mathrm{C}}_{\mathrm{p}}=6·{10}^{-7}·{\lambda }^5+{10}^{-5}·{\lambda }^4-65·{10}^{-5}·{\lambda }^3+2·{10}^{-5}·{\lambda }^2+76·{10}^{-3}·\lambda +0.007$$ E12

This coefficient has a maximum value equal to 0,35 (C_pmax = 0,35 ) and an optimal value of relative speed equal to 7 (λ = 7).

The block diagram presenting the aerodynamic part is shown in Figure 4.

Figure 4.

Block diagram of the aerodynamic part.

3.3 Gearbox model

The mechanical part of the wind turbine consists of the turbine shaft rotating slowly at speed Ω_t, the gearbox having the multiplication gain G and driving the generator at a speed Ω_g, by means of a fast secondary shaft.

The gearbox is a device that allows to multiply the turbine speed of Ω_t by a multiplication gain G to make it adapt to the rapid speed of the generator Ω_g. This device is considered ideal, because the gearbox elasticity, friction, and energy losses are considered negligible. The two equations mathematically modeling the operation of this device are given as follows:

$$\left\{\begin{array}{c}{\mathrm{T}}_{\mathrm{g}}=\frac{{\mathrm{T}}_{\mathrm{aer}}}{\mathrm{G}}\\ {\mathrm{\Omega }}_{\mathrm{t}}=\frac{{\mathrm{\Omega }}_{\mathrm{g}}}{\mathrm{G}}\end{array}\right.$$ E13

where Tg is torque on the generator shaft (N·m), T_aer is the aerodynamic torque of the wind turbine (N · m), Ω_g is the speed generator shaft (rad · s⁻¹), Ω_t is the turbine speed shaft (rad · s⁻¹), and G is the multiplication gain; it is given by G = N1/N2.

Figure 5 shows the gearbox model for determining the multiplication gain G, and Figure 6 shows the gearbox block diagram.

Figure 5.

The gearbox model.

Figure 6.

Block diagram of the gearbox.

The total inertia J consists of the turbine inertia J_t and the generator inertia J_g; it can be written according to the following equation [28]:

$$J=\frac{J_t}{G^2}+J_g$$ E14

The total viscous friction coefficient f_v consists of the generator friction coefficient f_g and the turbine friction coefficient f_t. The coefficient f_v can be expressed as follows:

$$f_v=\frac{f_t}{G^2}+f_g$$ E15

Therefore, the mechanical part can be modeled according to the diagram shown in Figure 7.

Figure 7.

The several coefficients of the wind turbine mechanical part.

The generator speed Ω_g depends on the total mechanical torque T_mec. This torque is the result of the electromagnetic torque of the generator T_em, the viscous friction torque T_v, and the torque applied on the generator shaft T_g.

$$T_{\mathit{mec}}=J·\frac{d{\Omega }_g}{\mathit{dt}}$$ E16

$$T_{\mathit{mec}}=T_g-T_{\mathit{em}}-T_v$$ E17

$$T_v=f·{\Omega }_g$$ E18

Therefore, from these previously established equations, the differential equation of the mechanical system dynamics is expressed by:

$$J·\frac{d{\Omega }_g}{\mathit{dt}}=T_g-T_{\mathit{em}}-T_{\mathit{vis}}$$ E19

The block diagram of the wind turbine mechanical part is presented in Figure 8.

Figure 8.

Block diagram of the wind turbine mechanical part.

The diagram block of the whole wind turbine system is given in Figure 9.

Figure 9.

Block diagram of the whole wind turbine system.

In order to continuously reach the maximum power point provided by a wind turbine, operating over a wide range of wind speed, the maximum power point tracking (MPPT) control technique is used. In this chapter, the MPPT control without controlling the mechanical speed is presented [32]. This control strategy is based on the assumption that the wind speed little varies in steady state compared to the electrical constants of the wind turbine system. Therefore, at the maximum power point, the relative speed ʎ is equal to its optimum value λ_opt, and the power coefficient C_p is equal to its maximum value C_p-max, while the reference electromagnetic torque $C_{\mathit{em}}^{\ast }$ is given by:

$$C_{\mathit{em}}^{\ast }=\frac{C_{\mathit{aer}-\mathit{est}}}{G}$$ E20

$$C_{\mathit{em}}^{\ast }=C_{p-\mathit{max}}·\frac{\rho ·\pi ·R^5·{\mathrm{\Omega }}_g^2}{2·{\lambda }_{\mathit{opt}}^3·G^3}$$ E21

For simplification, the parameter K is expressed as:

$$K=C_{p-\mathit{max}}·\frac{\rho ·\pi ·R^5}{2·{\lambda }_{\mathit{opt}}^3·G^3}$$ E22

Therefore:

$$C_{\mathit{em}}^{\ast }=K·{\mathrm{\Omega }}_g^2$$ E23

The reference electromagnetic torque is proportional to the square of the generator speed Ω_g. The block diagram which presents the MPPT control strategy without the measurement of the generator speed is shown in Figure 10.

Figure 10.

Block diagram of the MPPT control without mechanical speed.

3.4 DFIG model

The doubly fed induction generator (DFIG) is a three-phase asynchronous machine, powered by two sources: by its stator and its rotor at the same time. Its main advantage is that it offers the possibility of controlling the power flows for the hypo- and hypersynchronous modes, either in the motor or generator operation. It also allows the variable speed operation of the system where it is integrated.

The DFIG model in the stationary reference frame, noted (α, β), is given in the state representation [33] as follows:

$$\frac{d}{\mathit{dt}}\left[\begin{array}{c}{i}_{\mathit{s\alpha }}\\ i_{\mathit{s\beta }}\\ i_{\mathit{r\alpha }}\\ i_{\mathit{r\beta }}\end{array}\right]=\left[A\right]·\left[\begin{array}{c}{i}_{\mathit{s\alpha }}\\ i_{\mathit{s\beta }}\\ i_{\mathit{r\alpha }}\\ i_{\mathit{r\beta }}\end{array}\right]+\left[B\right]·\left[\begin{array}{c}{v}_{\mathit{s\alpha }}\\ v_{\mathit{s\beta }}\\ v_{\mathit{r\alpha }}\\ v_{\mathit{r\beta }}\end{array}\right]$$ E24

where i_sα and i_sβ are the stator currents in the stationary reference frame (α, β); i_rα and i_rβ are the rotor currents in the reference frame (α, β); v_sα and v_sβ are the stator stresses in the stationary reference frame (α, β); v_rα and v_rβ are the rotor voltages in the stationary reference frame (α, β).

The matrices A ∈ ℝ^n×n, B ∈ ℝ^n×m, and C ∈ ℝ^p×n are, respectively, the state matrix, the input or control matrix, and the output or observation matrix. They are, respectively, given by:

$$\mathrm{A}=\left[\begin{array}{cccc}\frac{-{\mathrm{R}}_{\mathrm{s}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}}& \left(\frac{\left(1-\mathrm{\sigma }\right)}{\mathrm{\sigma }}·\mathrm{\omega }+{\mathrm{\omega }}_{\mathrm{s}}\right)& \frac{\mathrm{M}·{\mathrm{R}}_{\mathrm{r}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}& \frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}}·\mathrm{\omega }\\ -\left(\frac{\left(1-\mathrm{\sigma }\right)}{\mathrm{\sigma }}·\mathrm{\omega }+{\mathrm{\omega }}_{\mathrm{s}}\right)& \frac{-{\mathrm{R}}_{\mathrm{s}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}}& -\frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}}·\mathrm{\omega }& \frac{\mathrm{M}·{\mathrm{R}}_{\mathrm{r}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}\\ \frac{\mathrm{M}·{\mathrm{R}}_{\mathrm{s}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}& -\frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{r}}}·\mathrm{\omega }& \frac{-{\mathrm{R}}_{\mathrm{r}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{r}}}& \left({\mathrm{\omega }}_{\mathrm{s}}-\frac{\mathrm{\omega }}{\mathrm{\sigma }}\right)\\ \frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{r}}}·\mathrm{\omega }& \frac{\mathrm{M}·{\mathrm{R}}_{\mathrm{s}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}& -\left({\mathrm{\omega }}_{\mathrm{s}}-\frac{\mathrm{\omega }}{\mathrm{\sigma }}\right)& \frac{-{\mathrm{R}}_{\mathrm{r}}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{r}}}\end{array}\right]$$

$$\mathrm{B}=\left[\begin{array}{cccc}\frac{1}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}}& 0& -\frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}& 0\\ 0& \frac{1}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}}& 0& -\frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}\\ -\frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}& 0& \frac{1}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{r}}}& 0\\ 0& -\frac{\mathrm{M}}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{s}}·{\mathrm{L}}_{\mathrm{r}}}& 0& \frac{1}{\mathrm{\sigma }·{\mathrm{L}}_{\mathrm{r}}}\end{array}\right]$$

$$\mathrm{C}=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right]$$

where R_s and L_s are, respectively, the single-phase resistance and the cyclic single-phase inductance of the stator winding; R_r and L_r are, respectively, the single-phase resistance and the cyclic single-phase inductance of the rotor winding; M is the mutual inductance between the stator phase and the rotor phase; $\sigma =1-\frac{M2}{L_s.L_r}$ is the leakage coefficient or the Blondel coefficient; ${\mathrm{\omega }}_{\mathrm{s}}$ is the synchronism angular speed [rad/s]; and $\mathrm{\omega }$ is the mechanical angular speed [rad/s].

In order to generate the reference rotor voltages which will be the input of the machine side converter, the stator flux-oriented vector control is applied to DFIG system [34].

3.5 Power converter models

The power electronic converters used consist of a rectifier made using semiconductors controlled at the opening and closing, and a three-phase voltage inverter consists of three reversible current switch arms, controlled at the opening and closing in the same time. Each arm consists of two switches, which contain each one insulated gate bipolar transistor (IGBT) and an antiparallel diode. The voltage capacitor DC allows the storage of the output rectifier energy. The passive filter type (L, R) is used to connect the inverter to the grid. Both converters used are controlled using pulse width modulation (PWM), and the power converter structure is given in Figure 11.

Figure 11.

Structure of the power converters (IGBT).

The input voltages of single phases of the rotor side converter (RSC) are described as follows:

$$\left\{\begin{array}{c}{V}_{\mathit{Sa}}=\frac{2·S_a-S_b-S_c}{3}·V_{\mathit{dc}}\\ V_{\mathit{Sb}}=\frac{2·S_b-S_a-S_c}{3}·V_{\mathit{dc}}\\ V_{\mathit{Sc}}=\frac{2·S_c-S_a-S_b}{3}·V_{\mathit{dc}}\end{array}\right.$$ E25

where S_i represents the switch states, supposedly ideal to facilitate the rectifier modeling, defined by:

$$S_i=\left\{\begin{array}{c}1,\overline{S_i}=0\\ 0,\overline{S_i}=1\end{array}\right.;\mathrm{i}=\mathrm{a},\mathrm{b},\mathrm{c}$$ E26

The rotor voltage equations and the DC capacitor equation are given, respectively:

$$\left[V_{\mathit{\text{rabc}}}\right]=R_r·\left[i_{\mathit{\text{rabc}}}\right]+L_r·\frac{d}{\mathit{dt}}\left[i_{\mathit{\text{rabc}}}\right]+\left[V_{\mathit{\text{Sabc}}}\right]$$ E27

$$C·\frac{d{V}_{\mathit{dc}}}{\mathit{dt}}=\left(S_a·i_{\mathit{Sa}}+S_b·i_{\mathit{Sb}}+S_c·i_{\mathit{Sc}}\right)-i_{\mathit{dc}}$$ E28

where V_rabc is the three-phase rotor voltage of DFIG [V]; i_rabc is the three-phase rotor current of DFIG [A]; C is the capacitor constant [F]; V_dc is the DC bus voltage [V]; i_dc is the DC output current [A].

The block diagrams of the rotor side converter (RSC) and the grid side converter (GSC) are given, respectively, in Figures 12 and 13.

Figure 12.

The block diagram of the RSC.

Figure 13.

The block diagram of the GSC.

4. Simulation results

The variable speed wind turbine model based on DFIG with a power of 3 Kw has been developed and simulated using MATLAB/Simulink software. The turbine and DFIG parameters are extracted from [30]. Some simulation results of the wind system modeled in this study are presented in the figures below. Figure 14 shows the random wind speed profile applied to the turbine. The wind speed varies between [6 m/s] and [11 m/s]. Figure 15 presents the aerodynamic power delivered by the wind turbine and it reached 3kw when the wind speed is up to 11 [m/s]. Figures 16 and 17 illustrate respectively the mechanical speed of the generator shaft and the speed of the turbine shaft. It can be noticed from the Figure 16 that during the simulation time (150 seconds), the generator operates in both hypo and hyper synchronous operating modes.

The Figures 18 and 19 show respectively the variation of the power coefficient Cp and the variation of the speed ratio λ, which coincide with the maximum power coefficient and with the optimal speed ratio.

Figure 14.

Wind speed profile [m/s].

Figure 15.

Aerodynamic power [W].

Figure 16.

Mechanical speed of DFIG [rpm].

Figure 17.

Turbine speed [rpm].

Figure 18.

Power coefficient Cp.

Figure 19.

Speed ratio ʎ.

The simulation results of the wind system electrical part, including the electrical characteristics of the DFIG, the power converters, and the capacitive bus, are presented in Figures 20–27.

Figure 20.

Three-phase stator current [A].

Figure 21.

Zoom on the stator currents [A].

Figure 22.

Three-phase rotor current of DFIG [A].

Figure 23.

Zoom on the rotor currents [A].

Figure 24.

The stator currents in the (α, β) reference [A].

Figure 25.

The rotor currents in the (α, β) reference [A].

Figure 26.

Three-phase stator voltage [V].

Figure 27.

Voltage of the DC capacitor [V].

By applying the first-order passive filter (R, L) to the square-wave signals, given in Figure 28, the rotor voltages in sinusoidal form are obtained and shown in Figure 29.

Figure 28.

The output voltages of the rotor side converter [V].

Figure 29.

Three-phase rotor voltage [V].

5. Conclusion

This chapter presents the modeling and simulation results of the most commonly used speed variable wind turbine driven by a doubly fed induction generator. In order to generate efficient and quick electrical power, the control techniques are applied, such as the MPPT control for wind turbine and the stator flux-oriented vector control are used for the generator. The wind turbine system and its control methods are established in the MATLAB/Simulink environment. The obtained signals are used for the design of fault diagnostic methods in future works.