Next Article in Journal
Surfactant-Assisted Perovskite Nanofillers Incorporated in Quaternized Poly (Vinyl Alcohol) Composite Membrane as an Effective Hydroxide-Conducting Electrolyte
Next Article in Special Issue
Wind Turbine Synchronous Reset Pitch Control
Previous Article in Journal
A Data-Driven Method for Energy Consumption Prediction and Energy-Efficient Routing of Electric Vehicles in Real-World Conditions
Previous Article in Special Issue
Automatic Tracking of the Modal Parameters of an Offshore Wind Turbine Drivetrain System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sliding Mode Control of a Variable- Speed Wind Energy Conversion System Using a Squirrel Cage Induction Generator

Department of Electrical Engineering, Kuwait University, P.O. Box 5969, Safat 13060, Kuwait
*
Author to whom correspondence should be addressed.
Energies 2017, 10(5), 604; https://doi.org/10.3390/en10050604
Submission received: 25 January 2017 / Revised: 19 April 2017 / Accepted: 19 April 2017 / Published: 1 May 2017
(This article belongs to the Collection Wind Turbines)

Abstract

:
This paper deals with the control of a variable-speed wind energy conversion (WEC) system using a squirrel cage induction generator (SCIG) connected to the grid through a back-to-back three phase (AC-DC-AC) power converter. The sliding mode control technique is used to control the WEC system. The objective of the controllers is to force the states of the system to track their desired states. One controller is used to regulate the generator speed and the flux so that maximum power is extracted from the wind. Another controller is used to control the grid side converter, which controls the DC bus voltage and the active and reactive powers injected into the grid. The performance of the controlled wind energy conversion system is verified through MATLAB simulations, which show that the controlled system performs well.

1. Introduction

Wind energy has been used to generate electricity for a long time. However, it is more prevalent nowadays because the cost of wind energy has continuously dropped, and it is approaching the competitive level of conventional energy [1]. Moreover, wind energy generation does not contribute to the pollution of the environment. The generators that are used to convert the mechanical power obtained from the wind turbine into electric power are generally either doubly-fed induction generators (DFIG), or squirrel cage induction generators (SCIG), or permanent magnet synchronous generators (PMSG). Therefore, variable-speed wind energy conversion systems can be broadly classified into three types: (i) DFIG-based wind energy conversion (WEC) systems; (ii) SCIG-based WEC systems; and (iii) PMSG-based WEC systems. For variable-speed operations, all three above-mentioned WEC systems need power electronic converters. The PMSG- and the SCIG-based systems need full-scale power electronic converters, whereas the DFIG-based systems need partial-scale power converters. Compared to the DFIG-based systems, the PMSG- and the SCIG-based systems could be more attractive due to the dropping cost of power electronics over time and due to the absence of brushes [2]. Hence, the main advantages of SCIG-based systems are their reliability, their ruggedness in design and their low operation and maintenance costs.
In an SCIG-based WEC system, the generator is coupled to the grid through back-to-back three-phase (AC-DC-AC) power converters. The AC/DC converter or the stator side converter controls the generator; the DC/AC converter or the grid side converter controls the DC bus voltage and the active and reactive powers injected into the grid. The controllers are designed so that maximum power is extracted from the wind at all wind speeds; thus, maximum system efficiency is achieved; the generated active power is transmitted through the DC-bus to the grid while ensuring the unity power factor. This paper proposes the use of the sliding mode control technique to control both the stator side and the grid side converters of an SCIG-based wind energy conversion system.
Extensive work was done on wind energy conversion systems. Mahela and Shaik [3] presented an overview of wind energy generation systems. Cheng and Zhu [4] surveyed the state of the art of wind energy conversion systems and their technologies. Chen et al. [5] reviewed the power electronics used for wind turbines systems. Papers by [1,6,7,8] discussed the state of the art of the control techniques applied to the different types of wind energy conversion systems.
Several control techniques were used to control the SCIG-based wind energy conversion system. Benchagra and his colleagues [9,10,11,12,13] published several papers on the control of SCIG-based WEC systems. Wang et al. [14] proposed a passivity-based robust controller for an SCIG-based WEC system. Hassan [15] used the backstepping control technique to control an SCIG-based WEC system. Zhao et al. [16] used an adaptive controller for maximum power point tracking (MPPT) for an SCIG-based WEC system. Heydari-doostabad et al. [17] proposed an optimal linear quadratic regulator (LQR) controller for a SCIG-based WEC system using a static compensator (STATCOM). Baloch et al. [18] employed the feedback linearization technique through field control concepts to control an SCIG-based WEC system. Domínguez-García [19] used indirect vector control for the squirrel cage induction generator WEC system. The sliding mode control (SMC) technique was used for induction generator-based WEC systems in [20,21,22,23,24]. Moreover, other types of controllers were also designed for squirrel cage-based WEC systems; for example, see [25,26,27,28,29,30,31].
The sliding mode technique is a very suitable approach to deal with the control of wind energy conversion systems. This is the case because this technique is robust with respect to variations in the parameters of the system and to bounded external disturbances. Therefore, the SMC technique is selected to control the squirrel cage induction generator-based WEC system.
The contribution of this paper involves the design of two controllers, which force the states of a squirrel cage induction generator-based wind energy conversion system to track the desired states of the system and the transmission of the the generated power to the grid. The first sliding mode controller enables maximum extraction of power from the wind at different speeds of the wind; the second sliding mode controller ensures that the generated active power is injected through the DC-bus to the grid with the unity power factor. The proposed controllers are validated through MATLAB simulation.
The rest of the paper is organized as follows. The model of an SCIG-based wind energy conversion system is presented in Section 2. The model of the system consists of two parts. The first part is composed of the model of the wind turbine and the induction generator; the second part consists of the model of the grid side converter. The computation of the desired states of the WEC system and the computation of the reference inputs are presented in Section 3. The design of sliding mode controllers for the WEC system is proposed in Section 4. Simulation results for the controlled system are presented and discussed in Section 5. Finally, the conclusion is given in Section 6.

2. Model of the Wind Energy Conversion System

This section presents the model of the wind energy conversion system, which is depicted in Figure 1. The model of the system consists of two parts. One part deals with the model of the turbine and the squirrel cage induction generator (SCIG). This part of the system will be controlled through the stator side AC/DC Converter. The second part of the model deals with the grid side DC/AC converter sub-system.

2.1. The Model of the Turbine and the Induction Generator

The mechanical power extracted from the wind can be expressed as follows:
P m = 1 2 C p ( λ , β ) ρ A V 3
where P m is the mechanical output power of the wind turbine; the power coefficient is represented by C p ( λ , β ) ; λ is the tip speed ratio ( T S R ) of the turbine blade; and β is the pitch angle. The constant ρ represents the air density; A is the area swept by the rotor blades; and V is the speed of the wind.
The gearbox transfers power from the rotating wind turbine shaft to the rotating generator shaft. The angular speed of the generator ω m is related to the angular speed of the wind turbine ω t through the following equation:
ω m = G ω t
where G is the gear ratio.
The tip speed ratio is related to the angular speed of the wind turbine ω t and the angular speed of the generator ω m through the following equation:
λ = R ω t V = R ω m G V
where R is the radius of the turbine blades.
The power coefficient C p is a nonlinear function of λ ; several nonlinear models of C p can be found in the literature. The model of the power coefficient adopted in this paper is as follows [32]:
C p ( λ , β ) = a 1 a 2 λ ¯ a 3 β a 4 e x p ( a 5 λ ¯ ) + a 6 λ
where λ ¯ = 1 / [ 1 / ( λ + a 7 β ) a 8 / ( β 3 + 1 ) ] with a 1 = 0.5109 , a 2 = 116 , a 3 = 0.4 , a 4 = 5 , a 5 = 21 , a 6 = 0.0068 , a 7 = 0.08 and a 8 = 0.035 . Figure 2 depicts the power coefficient C p versus the tip speed ratio λ for different pitch angles. The figure shows that the maximum value of the power coefficient C p is C p m a x = 0.47 ; this value is obtained when the pitch angle β equals zero degrees and the corresponding optimum tip speed ratio is λ o p t = 8.1 .
Note that the maximum mechanical power extracted from the wind is such that:
P m m a x = 1 2 C p m a x λ o p t ρ A V 3
Figure 3 shows the wind turbine mechanical (aerodynamic) power P m versus the generator speed ω m for different values of the wind velocity V. Furthermore, the graph of the maximum mechanical power P m m a x is shown in black in Figure 3.
The mechanical (or aerodynamic) torque is related to the the power P m captured by the wind turbine through the equation:
T m = P m ω t
The mechanical equation of the system is such that:
J d ω m d t = T e m T m B ω m
In the above equation, T e m is the electromagnetic torque of the induction generator, and the inertia J is the total inertia of the turbine and the induction generator. This inertia is calculated such that J = J t G 2 + J g where J t and J g represent the inertia of the turbine and the inertia of the generator respectively, and G is the gearbox ratio. Furthermore, B is the total external damping of the system.
In the direct-quadrature frame, the differential equations governing the flux linkages of the stator and the rotor of the generator are as follows [9]:
d ψ d s d t = R s i d s + ω s ψ q s + v d s
d ψ q s d t = R s i q s ω s ψ d s + v q s
d ψ d r d t = R r i d r + ω r ψ q r
d ψ q r d t = R r i q r ω r ψ d r
where ψ d s is the stator direct axis flux linkage and ψ q s is the stator quadrature axis flux linkage. The flux ψ d r is the rotor direct axis flux linkage, and ψ q r is the rotor quadrature axis flux linkage. The currents i d s and i q s are the stator currents in the direct-quadrature frame; and i d r and i q r are the rotor currents in the direct-quadrature frame. The voltages v d s and v q s represent the voltages of the windings of the stator in the direct-quadrature frame. The resistances R s and R r are the stator and rotor phase resistances. Furthermore, ω s is the synchronous speed, and ω r is the angular speed of the rotor.
The relationship between the stator and rotor flux linkages and the currents can be expressed as:
ψ d s = L s i d s + L m i d r
ψ q s = L s i q s + L m i q r
ψ d r = L r i d r + L m i d s
ψ q r = L r i q r + L m i q s
where L s is the stator inductance and L r is the rotor inductance; L m is the mutual inductance between the stator and the rotor.
The equation of the electromagnetic torque T e m in the d-q frame can be written as follows:
T e m = p L m L r ( ψ d r i q s ψ q r i d s )
where p is the number of pole pairs.
The active and reactive powers transmitted by the stator P s and Q s can be expressed as:
P s = 3 2 ( v d s i d s + v q s i q s )
Q s = 3 2 ( v d s i q s v q s i d s )
Combining Equations (7)–(16), we can write the following system of ordinary differential equations (o.d.e.’s):
d i d s d t = c 1 i d s + ω s i q s + c 2 ψ d r + c 3 ω m ψ q r + c 4 v d s d i q s d t = c 1 i q s ω s i d s + c 2 ψ q r c 3 ω m ψ d r + c 4 v q s d ψ d r d t = c 5 i d s c 6 ψ d r + ( ω s p ω m ) ψ q r d ψ q r d t = c 5 i q s c 6 ψ q r ( ω s p ω m ) ψ d r d ω m d t = c 7 ( ψ d r i q s ψ q r i d s ) c 8 ω m c 9 T m
where, c 1 = L r 2 R s + L m 2 R r σ L s L r 2 , c 2 = L m R r σ L s L r 2 , c 3 = p R r σ L s L r , c 4 = 1 σ L s , c 5 = L m R r L r , c 6 = R r L r , c 7 = p L m J L r , c 8 = B J , c 9 = 1 J .
By fixing the direct axis of the rotating d-q reference frame on the rotating flux vector, we obtain ψ q r = 0 . Therefore, the corresponding equation is removed from the model of the system, and the system of o.d.e’s is reduced accordingly. Furthermore, we will define the state variable vector x T comprised of the states x 1 x 4 , such that:
x T = [ x 1 x 2 x 3 x 4 ] T = [ i d s i q s ψ d r ω m ] T
Therefore, the model of the wind turbine and the induction generator in (20) can be written in compact form as follows:
x ˙ 1 = c 1 x 1 + ω s x 2 + c 2 x 3 + c 4 v d s x ˙ 2 = c 1 x 2 ω s x 1 c 3 x 3 x 4 + c 4 v q s x ˙ 3 = c 5 x 1 c 6 x 3 x ˙ 4 = c 7 x 2 x 3 c 8 x 4 c 9 T m
We will define the desired state vector x d such that, x d = [ x 1 d x 2 d x 3 d x 4 d ] T = [ i d s * i q s * ψ d r * ω m * ] T where x i d ( i = 1 , 2 , , 4 ) are the desired values of the states x i ( i = 1 , 2 , , 4 ) of the WEC system in (21). Since the desired states have to be an operating point of the WEC system, then these states must satisfy the model of the WEC system in (21). Therefore, the desired states of the model of the wind turbine and the induction generator system are governed by the following set of differential equations:
x ˙ 1 d = c 1 x 1 d + ω s x 2 d + c 2 x 3 d + c 4 v d s r
x ˙ 2 d = c 1 x 2 d ω s x 1 d c 3 x 3 d x 4 d + c 4 v q s r
x ˙ 3 d = c 5 x 1 d c 6 x 3 d
x ˙ 4 d = c 7 x 2 d x 3 d c 8 x 4 d c 9 T m r
Note that v d s r , v q s r and T m r in (22), (23) and (25) are the reference input voltages and the reference mechanical torque of the WEC system. The computation of desired states and the reference inputs will be discussed in Section 3.
Define the errors e i ( i = 1 , 2 , ..... , 4 ) such that:
e i = x i x i d
Using Equations (21)–(26), the dynamic model of the error of the wind turbine and the induction generator system can be written such that:
e ˙ 1 = c 1 e 1 + ω s e 2 + c 2 e 3 + c 4 u 1 e ˙ 2 = c 1 e 2 ω s e 1 c 3 ( e 3 e 4 + x 4 d e 3 + x 3 d e 4 ) + c 4 u 2 e ˙ 3 = c 5 e 1 c 6 e 3 e ˙ 4 = c 7 ( e 2 e 3 + x 3 d e 2 + x 2 d e 3 ) c 8 e 4 c 9 ( T m T m r )
where the inputs u 1 and u 2 are defined such that:
u 1 = v d s v d s r
u 2 = v q s v q s r
In Section 4, we will use the sliding mode control technique to design the controllers u 1 and u 2 to force the errors e 1 e 4 to converge to zero as t tends to infinity.

2.2. The Model of the Grid Side Converter

The following equations relate the grid, converter voltages and the line currents [10]:
v g 1 = R t i g 1 + L t d i g 1 d t + v i 1 v g 2 = R t i g 2 + L t d i g 2 d t + v i 2 v g 3 = R t i g 3 + L t d i g 3 d t + v i 3
where, v g 1 , v g 2 and v g 3 are the three phase grid voltages; and i g 1 , i g 2 and i g 3 are the line currents. The voltages, v i 1 , v i 2 and v i 3 are the voltages of the voltage source inverter (VSI). The resistance and the inductance of the grid-side transmission line are R t and L t , respectively.
Using the d-q transformation, the equations in (31) can be written as:
v d g = v d i R t i d g L t d i d g d t + ω L t i q g v q g = v q i R t i q g L t d i q g d t ω L t i d g
where v d g is the direct component of the grid voltage, and v q g is the quadrature component of the grid voltage. The currents i d g and i q g are the direct and quadrature output line currents. The voltages v d i and v q i are the direct and quadrature VSI input voltages. Furthermore, ω = 2 π f , where f is the frequency of the grid.
Let the constants c 10 and c 11 be such that c 10 = R t L t and c 11 = 1 L t . Therefore, the above equations can be written as:
d i d g d t = c 10 i d g + ω i q g + c 11 ( v d i v d g ) d i q g d t = c 10 i q g ω i d g + c 11 v q i
The output active power P g and the output reactive power Q g of the wind energy conversion system can be written as:
P g = 3 2 ( v d g i d g + v q g i q g ) Q g = 3 2 ( v q g i d g v d g i q g )
Setting the initial angle of the direct and quadrature reference frame to π 2 and setting the initial angle of Phase 1 to 0o yield the voltage v q g = 0 and v d g = V g r i d , where V g r i d is the grid voltage. Hence, the output active and reactive powers of the system can be written as:
P g = 3 2 v d g i d g
Q g = 3 2 v d g i q g
The voltage equation for the DC-bus is such that:
C d U d c d t = i d c = i s i g
where U d c is the voltage of the DC-link, i s is the current from the induction generator side and i g is the current from the grid side; refer to Figure 4. The capacitor of the DC-link is denoted as C.
Remark 1.
The DC-bus voltage and the output active and reactive powers P g and Q g of the wind energy conversion system are controlled by the direct and quadrature components of the line current i g . Hence, the term i s in Equation (36) can be assumed as a disturbance in the control [15].
Motivated by the work in [9], we multiply both sides of (36) by U d c and use the fact that the SCIG active power can be expressed as P s = U d c i s and that P g = U d c i g , to obtain the following equation:
1 2 C d U d c 2 d t = P s P g
Using the expression P g = 3 2 v d g i d g , the above equation can be written as:
d U d c 2 d t = 3 C v d g i d g + 2 C P s
We will define the state variable vector x C , which contains the states of the converter x 5 x 7 , such that:
x C = [ x 5 x 6 x 7 ] T = [ i d g i q g U d c 2 ] T
Using Equations (32) and (38), the dynamic model of the grid side converter can be written as follows:
x ˙ 5 = c 10 x 5 + ω x 6 + c 11 ( v d i v d g ) x ˙ 6 = c 10 x 6 ω x 5 + c 11 v q i x ˙ 7 = 3 C x 5 v d g + 2 C P s
Let x 5 d = i d g * , x 6 d = i q g * and x 7 d = U d c 2 * represent the desired values of i d g , i q g and U d c 2 respectively. Furthermore, let c 12 be such that c 12 = 3 C . Define the errors e 5 , e 6 and e 7 such that:
e i = x i x i d ( i = 5 , 6 , 7 )
The model of the errors of the grid side converter can be written as:
e ˙ 5 = c 10 e 5 + ω e 6 + u 3 e ˙ 6 = c 10 e 6 ω e 5 + u 4 e ˙ 7 = c 12 e 5 v d g + d p
where:
u 3 = c 10 x 5 d + ω x 6 d x ˙ 5 d + c 11 ( v d i v d g )
u 4 = c 10 x 6 d ω x 5 d x ˙ 6 d + c 11 v q i
d p = c 12 x 5 d v d g x ˙ 7 d + 2 C P s
The objective of controlling the grid side converter is to design the controllers u 3 and u 4 so that the errors e 5 , e 6 and e 7 converge to zero even with the existence of the bounded “disturbance” d p . We will use the sliding mode control technique to accomplish this task.

3. Computation of the Desired States and the Reference Inputs of the System

This section deals with the computation of the desired states and the reference inputs of the wind energy conversion system. Assuming the d-axis of the rotating frame is aligned with the rotor flux vector, then the quadrature component of the flux of the rotor is equal to zero ( ψ q r = 0 ), and the direct component of the flux of the rotor is always at a maximum constant value. Therefore, we can write the following:
x 3 d = ψ d r * = ψ r m a x
where ψ d r * is the desired value of ψ d r , and:
ψ r m a x = L m i d s
Furthermore, the desired value of the current i d s can be calculated from (47), such that:
x 1 d = i d s * = 1 L m ψ d r * = 1 L m x 3 d
The wind turbine changes its speed by following the maximum of the power coefficient C p m a x , so that maximum power is extracted from the wind. Using Equation (3), the desired value of the rotor speed ω m is such that:
x 4 d = ω m * = G λ o p t V R
where λ o p t is the optimum value of tip speed ratio of the wind turbine.
The desired value of the current i q s can be calculated from the electromagnetic torque given in (16), such that:
x 2 d = i q s * = 2 3 L r p L m T m r x 3 d
Note that T m r is the input reference mechanical torque, and it is set to the maximum mechanical (aerodynamic) torque such that T m r = P m m a x x 4 d .
The desired values of the direct and quadrature line currents i d g and i q g can be found from Equations (17), (34) and (35) after setting P g = P s to transfer all active power from the side of the stator to the grid side while maintaining Q g = 0 to ensure the unity power factor. Hence:
x 5 d = i d g * = 2 3 P s V g r i d
and the desired value of the quadrature line current i q g is such that:
x 6 d = i q g * = 0
The DC-link voltage U d c must be maintained constant at a preset value. The desired value of the voltage of the DC-link U d c is such that:
U d c * = U N
where U N is the rated DC-bus voltage.
On the other hand, the input reference voltages v d s r and v q s r are computed using Equations (22) and (23) such that:
v d s r = 1 c 4 c 1 x 1 d ω s x 2 d c 2 x 3 d
v q s r = 1 c 4 c 1 x 2 d + ω s x 1 d + c 3 x 3 d x 4 d + 2 3 L r p L m T ˙ m r x 3 d

4. Control of the Wind Energy Power System

In this section, we will use the sliding mode control technique to control the wind energy conversion system. Figure 4 depicts the block diagram representation of the controlled SCIG-based WEC system.

4.1. Design of a Sliding Mode Controller for the Turbine and the Induction Generator

To facilitate the design of a controller for the wind turbine and induction generator system, we define the following nonlinear transformation:
z 1 = e 3 z 2 = c 5 e 1 c 6 e 3 z 3 = e 4 z 4 = c 7 ( e 2 e 3 + x 3 d e 2 + x 2 d e 3 ) c 8 e 4 c 9 ( T m T m r )
Using the above transformation, the dynamic model of wind turbine and induction generator error system in (27) can be written as follows:
z ˙ 1 = z 2 z ˙ 2 = f 1 + g 1 u 1 z ˙ 3 = z 4 z ˙ 4 = f 2 + g 2 u 2
where f 1 , f 2 , g 1 and g 2 are such that:
f 1 = c 5 c 1 e 1 + ω s e 2 + c 2 e 3 c 6 c 5 e 1 c 6 e 3
f 2 = c 7 ( e 3 + x 3 d ) c 1 e 2 ω s e 1 c 3 ( e 3 e 4 + x 4 d e 3 + x 3 d e 4 ) + c 7 ( e 2 + x 2 d ) c 5 e 1 c 6 e 3 + c 7 e 2 x ˙ 3 d + c 7 e 3 x ˙ 2 d c 9 ( T ˙ m T ˙ m r ) c 8 c 7 ( e 2 e 3 + x 3 d e 2 + x 2 d e 3 ) c 8 e 4 c 9 ( T m T m r )
g 1 = c 4 c 5
g 2 = c 4 c 7 ( e 3 + x 3 d ) .
Let the design constants β 1 , β 2 , K 1 , K 2 , W 1 and W 2 be positive scalars. Define the sliding surfaces s 1 and s 2 such that:
s 1 = z 2 + β 1 z 1
s 2 = z 4 + β 2 z 3
Proposition 1.
The following sliding mode control law:
u 1 = 1 g 1 f 1 β 1 z 2 K 1 s 1 W 1 s i g n ( s 1 )
u 2 = 1 g 2 f 2 β 2 z 4 K 2 s 2 W 2 s i g n ( s 2 )
when applied to the error system given by (27) guarantees that the errors e i ( i = 1 , , 4 ) converge to zero as t tends to infinity. Therefore, the states of the wind turbine and the induction generator given by (21) converge to their desired values as t tends to infinity.
Proof. 
Taking the time derivatives of s 1 and s 2 in (62) and (63), and using the dynamic model of the errors in (27) and the control laws given by (64) and (65), we obtain:
s ˙ 1 = z ˙ 2 + β 1 z ˙ 1 = f 1 + g 1 u 1 + β 1 z 2 = K 1 s 1 W 1 s i g n ( s 1 )
and:
K 2 s 2 W 2 s i g n ( s 2 ) = z ˙ 4 + β 2 z ˙ 3 = f 2 + g 2 u 2 + β 2 z 4 = K 2 s 2 W 2 s i g n ( s 2 )
From (66) and (67), we conclude that s ˙ i = K i s i W i s i g n ( s i ) for i = 1 , 2 . It is clear that s i s ˙ i = K i s i 2 W i | s i | for i = 1 , 2 . Hence, the equations given by (66) and (67) guarantee that s i s ˙ i < 0 when s i for i = 1 , 2 . Therefore, the trajectories of these discontinuous dynamic equations exhibit a finite time reachability to zero from any given initial conditions provided that the gains K 1 , K 2 , W 1 and W 2 are chosen to be sufficiently large, strictly positive scalars.
Since s 1 and s 2 are driven to zero in finite time, then after such a finite time, we have z 2 + β 1 z 1 = 0 and z 4 + β 2 z 3 = 0 . Thus, z 2 = β 1 z 1 and z 4 = β 2 z 3 . Using the equations in (57), we can write: z ˙ 1 = z 2 = β 1 z 1 , and z ˙ 3 = z 4 = β 2 z 3 . The solutions of these two o.d.e.’s are, z 1 ( t ) = e x p ( β 1 t ) z 1 ( 0 ) and z 3 ( t ) = e x p ( β 2 t ) z 3 ( 0 ) . Because β 1 and β 2 are positive scalars, it can be concluded that the states z 1 ( t ) and z 3 ( t ) tend to zero as t goes to infinity. Furthermore, since s 1 = 0 and s 2 = 0 on the sliding surfaces, it can be concluded that the variables z 2 ( t ) and z 4 ( t ) are driven to zero as t tends to infinity.
Using the transformation given by (56), it can be deduced that the errors e 1 ( t ) , e 2 ( t ) , e 3 ( t ) and e 4 ( t ) asymptotically converge to zero as t tends to infinity. ☐
Hence, it can be concluded that the application of the controller (64) and (65) to the model of the wind turbine and the induction generator system (21) guarantees the convergence of the states i d s , i q s , ψ d r and ω m to their desired values.
Remark 2.
The control laws v d s and v q s , which will be used to control the model of the wind turbine and the induction generator given by (21), are as follows:
v d s = u 1 + v d s r
v q s = u 2 + v q s r
where v d s r and v q s r are obtained from (54) and (55).

4.2. Design of a Sliding Mode controller for the Grid Side Converter

This subsection deals with the design of a sliding mode controller for the grid side converter.
Let β 3 , K 3 , K 4 , W 3 and W 4 be positive scalars. Furthermore, let B d be a positive scalar bounding the “disturbance” d p defined in (45) such that | d ˙ p + β 3 d p | B d .
Define the sliding surfaces s 3 and s 4 such that:
s 3 = e 6
s 4 = e ˙ 7 + β 3 e 7
Proposition 2.
The following sliding mode control law:
u 3 = 1 c 12 v d g c 12 v d g ( c 10 e 5 + ω e 6 ) + c 12 v ˙ d g e 5 + c 12 β 3 v d g e 5 K 4 s 4 ( B d + W 4 ) s i g n ( s 4 )
u 4 = c 10 e 6 + ω e 5 K 3 s 3 W 3 s i g n ( s 3 )
when applied to the model of the errors of the grid side converter given by (43) guarantees that the errors e 5 , e 6 and e 7 converge to zero as t tends to infinity. Therefore, the states of the grid side converter converge to their desired values as t tends to infinity.
Proof. 
Let the Lyapunov function candidate V x be such that:
V x = 1 2 s 3 2 + 1 2 s 4 2
Taking the derivative of V x with respect to time along the trajectories of (43) and using the control law given by (72) and (73), we obtain:
V ˙ x = s 3 s ˙ 3 + s 4 s ˙ 4 = s 3 e ˙ 6 + s 4 ( e ¨ 7 + β 3 e ˙ 7 ) = s 3 c 10 e 6 ω e 5 + u 4 + s 4 c 12 v d g ( c 10 e 5 + ω e 6 + u 3 ) c 12 e 5 v ˙ d g + d ˙ p + β 3 ( c 12 e 5 v d g + d p ) = K 3 s 3 2 W 3 s 3 s i g n ( s 3 ) + s 4 c 12 v d g ( c 10 e 5 + ω e 6 + u 3 ) c 12 e 5 v ˙ d g + d ˙ p + β 3 ( c 12 e 5 v d g + d p ) K 3 s 3 2 W 3 s 3 s i g n ( s 3 ) + s 4 c 12 v d g ( c 10 e 5 + ω e 6 + u 3 ) c 12 e 5 v ˙ d g β 3 c 12 e 5 v d g + | s 4 | B d = K 3 s 3 2 W 3 s 3 s i g n ( s 3 ) + s 4 c 12 v d g ( c 10 e 5 + ω e 6 + u 3 ) c 12 e 5 v ˙ d g c 12 β 3 e 5 v d g + B d s i g n ( s 4 ) = K 3 s 3 2 W 3 s 3 s i g n ( s 3 ) K 4 s 4 2 W 4 s 4 s i g n ( s 4 ) = K 3 s 3 2 W 3 | s 3 | K 4 s 4 2 W 4 | s 4 |
Therefore, V ˙ x is negative definite. Hence, it can be concluded that the sliding surfaces s 3 and s 4 exhibit a finite time reachability to zero from any given initial conditions provided that the gains K 3 , K 4 , W 3 and W 4 are chosen to be sufficiently large, strictly positive scalars. Therefore, the error e 6 ( t ) = s 3 ( t ) converges to zero in finite time. Furthermore, because β 3 is a positive scalar, we can conclude that the error e 7 ( t ) converge to zero as t tends to infinity.
Moreover, since s 3 and s 4 are driven to zero in finite time, then the controller u 3 after such time is: u 3 = c 10 e 5 v ˙ d g v d g e 5 β 3 e 5 . Hence, after such time, the dynamics of the error e 5 is such that e ˙ 5 = c 10 e 5 + u 3 = ( v ˙ d g v d g + β 3 ) e 5 . Therefore, the error e 5 ( t ) asymptotically converges to zero provided that the control gain β 3 is chosen to be sufficiently large, strictly positive scalar. ☐
Hence, it can be concluded that the application of the controller (72) and (73) to the model of the grid side converter (41) guarantees the convergence of the states i d g , i q g and U d c to their desired values.
Remark 3.
The control laws v d i and v q i , which will be used to control the model of the grid side converter given by (41), are as follows:
v d i = 1 c 11 u 3 + c 10 x 5 d ω x 6 d + x ˙ 5 d + v d g
v q i = 1 c 11 u 4 + c 10 x 6 d + ω x 5 d + x ˙ 6 d .

5. Simulation Results of the Proposed Sliding Mode Control Scheme

The performance of the SCIG-based wind energy conversion system given by (21) and (41) using the sliding mode controllers (68), (69), (76) and (77) was simulated using the MATLAB software. The parameters of the system used for the simulation studies are given in Table 1.
The model of the wind speed used for simulation purposes is given by the following formula [33],
V = 10 + 0.55 [ sin ( 0.0625 w ) 0.875 sin ( 0.1875 w ) + 0.75 sin ( 0.3125 w ) 0.625 sin ( 0.625 w ) + 0.5 sin ( 1.875 w ) + 0.25 sin ( 3.125 w ) + 0.125 sin ( 6.25 w ) ]
with w = 2 π 10 t . The profile of the wind speed versus time is depicted in Figure 5.
The simulation results are presented in Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17. The currents i d s and i d s * versus time and the currents i q s and i q s * versus time are depicted in Figure 6 and Figure 7, respectively. It can be seen that the direct and quadrature currents are tracking their desired values. The fluxes ψ d r and ψ d r * versus time are depicted in Figure 8. Figure 9 shows the angular speeds ω m and ω m * of the generator versus time. It can be noticed that ω m is tracking the desired speed and following the maximum power coefficient C p m a x = 0.47 , as shown in Figure 10. The electromagnetic torque T e m versus time and the stator active power P s versus time are shown in Figure 11 and Figure 12. For the grid side, the direct currents i d g and i d g * versus time and the quadrature currents i q g and i q g * versus time are presented in Figure 13 and Figure 14, respectively. It is clear from these figures that the direct and quadrature currents track their desired values. The voltages U d c and U d c * versus time are depicted in Figure 15. It is clear that the DC link voltage is almost constant at all times. Finally, the output active power of the grid P g versus time and the output reactive power of the grid Q g versus time are shown in Figure 16 and Figure 17, respectively. Figure 16 shows that approximately all active power from the stator side has been transferred to the grid side. Furthermore, Figure 17 shows that the output reactive power of the grid Q g is approximately zero, which ensures that the power factor is approximately one.
The simulation results clearly indicate that the proposed sliding mode control scheme when applied to the a squirrel cage induction generator-based wind energy conversion system is able to force the states of the system to converge to their desired values. Hence, the controller enables maximum extraction of power from the wind at different speeds of the wind, and it ensures that the generated active power is transmitted through the DC-bus to the grid while maintaining the unity power factor.

6. Conclusions

This paper dealt with the control of a squirrel cage induction generator-based wind energy conversion system. Sliding mode controllers are proposed for the WEC system so that the actual states of the system track the desired states of the system. The control laws are designed to control the stator side converter, as well as the grid side converter. The stability of the controlled system is analyzed. Furthermore, simulation results are presented to verify the validity of the proposed control scheme. The simulation results indicate that the proposed sliding mode control scheme works well.
Future work will address the design of observer-based controllers for SCIG-based wind energy conversion systems.

Author Contributions

Mohamed Zribi and Muthana Alrifai contributed equally to this work. Mohamed Rayan contributed to the simulations of the performance of the system. The final manuscript was approved by all authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baloch, M.H.; Wang, J.; Kaloi, G.S. A Review of the State of the Art Control Techniques for Wind Energy Conversion System. Int. J. Renew. Energy Res. (IJRER) 2016, 6, 1276–1295. [Google Scholar]
  2. Li, H.; Chen, Z. Overview of Different wind generator systems and their comparisons. IET Renew. Power Gener. 2008, 2, 123–138. [Google Scholar] [CrossRef]
  3. Mahela, O.P.; Shaik, A.G. Comprehensive overview of grid interfaced wind energy generation systems. Renew. Sustain. Energy Rev. 2016, 57, 260–281. [Google Scholar] [CrossRef]
  4. Cheng, M.; Zhu, Y. The state of the art of wind energy conversion systems and technologies: A review. Energy Convers. Manag. 2014, 88, 332–347. [Google Scholar] [CrossRef]
  5. Chen, Z.; Guerrero, J.M.; Blaabjerg, F. A review of the state of the art of power electronics for wind turbines. IEEE Trans. Power Electron. 2009, 24, 1859–1875. [Google Scholar] [CrossRef]
  6. Nguyen, H.M.; Naidu, D.S. Evolution of wind turbine control systems. In Encyclopedia of Control Systems, Robotics, and Automation; EOLSS Publishers: Oxford, UK, 2010. [Google Scholar]
  7. Pao, L.Y.; Johnson, K.E. Control of wind turbines. IEEE Control Syst. 2011, 31, 44–62. [Google Scholar] [CrossRef]
  8. Babu, N.R.; Arulmozhivarman, P. Wind energy conversion systems-a technical review. J. Eng. Sci. Technol. 2013, 8, 493–507. [Google Scholar]
  9. Benchagra, M.; Errami, Y.; Hilal, M.; Maaroufi, M.; Cherkaoui, M.; Ouassaid, M. New control strategy for inductin generator-wind turbine connected grid. In Proceedings of the 2012 International Conference on Multimedia Computing and Systems (ICMCS), Tangiers, Morocco, 10–12 May 2012; pp. 1043–1048. [Google Scholar]
  10. Benchagra, M.; Hilal, M.; Errami, Y.; Ouassaid, M.; Maaroufi, M. Modeling and control of SCIG based variable-speed with power factor control. Int. Rev. Model. Simul. (IREMOS) 2011, 4, 1007–1014. [Google Scholar]
  11. Benchagra, M.; Maaroufi, M.; Ouassaid, M. A performance comparison of linear and nonlinear control o f a SCIG-wind farm connecting to a distribution network. Turk. J. Electr. Eng. Comput. Sci. 2014, 22, 1–11. [Google Scholar] [CrossRef]
  12. Benchagra, M.; Maaroufi, M. Control of wind farm connected distribution network. In Proceedings of the 10th IEEE International Conference on Networking, Sensing and Control (ICNSC), Evry, France, 10–12 April 2013; pp. 467–472. [Google Scholar]
  13. Benchagra, M.; Maaroufi, M.; Ouassaid, M. Nonlinear MPPT control of squirrel cage induction generator-wind turbine. J. Theor. Appl. Inf. Technol. 2012, 35, 26–33. [Google Scholar]
  14. Wang, P.; Wang, H.; Cai, X.; Han, Z. Passivity-based robust controller design for a variable speed wind energy conversion system. Turk. J. Electr. Eng. Comput. Sci. 2016, 24, 558–570. [Google Scholar] [CrossRef]
  15. Hassan, A. Nonlinear Control of Variable wind Speed Conversion System Based on a Squirrel Cage Induction Generator (SCIG). Am. J. Eng. Appl. Sci. 2015, 8, 275–284. [Google Scholar] [CrossRef]
  16. Zhao, H.; Wu, Q.; Rasmussen, C.N.; Blanke, M. calL1 Adaptive Speed Control of a Small Wind Energy Conversion System for Maximum Power Point Tracking. IEEE Trans. Energy Convers. 2014, 29, 576–584. [Google Scholar] [CrossRef]
  17. Heydari-doostabad, H.; Khalghani, M.R.; Khooban, M.H. A novel control system design to improve LVRT capability of fixed speed wind turbines using STATCOM in presence of voltage fault. Int. J. Electr. Power Energy Syst. 2016, 77, 280–286. [Google Scholar] [CrossRef]
  18. Baloch, M.H.; Wang, J.; Kaloi, G.S. Stability and nonlinear controller analysis of wind energy conversion system with random wind speed. Int. J. Electr. Power Energy Syst. 2016, 79, 75–83. [Google Scholar] [CrossRef]
  19. Domínguez-García, J.L.; Gomis-Bellmunt, O.; Trilla-Romero, L.; Junyent-Ferré, A. Indirect vector control of a squirrel cage induction generator wind turbine. Comput. Math. Appl. 2012, 64, 102–114. [Google Scholar] [CrossRef]
  20. Ouassaid, M.; Elyaalaoui, K.; Cherkaoui, M. Sliding Mode Control of Induction Generator Wind Turbine Connected to the Grid. In Advances and Applications in Nonlinear Control Systems; Springer International Publishing: New York, NY, USA, 2016; pp. 531–553. [Google Scholar]
  21. Amimeur, H.; Aouzellag, D.; Abdessemed, R.; Ghedamsi, K. Sliding mode control of a dual-stator induction generator for wind energy conversion systems. Int. J. Electr. Power Energy Syst. 2012, 42, 60–70. [Google Scholar] [CrossRef]
  22. De Battista, H.; Mantz, R.J.; Christiansen, C.F. Dynamical sliding mode power control of wind driven induction generators. IEEE Trans. Energy Convers. 2010, 15, 451–457. [Google Scholar] [CrossRef]
  23. Mi, Y.; Bao, X.; Yang, Y.; Zhang, H.; Wang, P. The sliding mode pitch angle controller design for squirrel-cage induction generator wind power generation system. In Proceedings of the 33rd Chinese Control Conference (CCC), Nanjing, China, 28–30 July 2014; pp. 8113–8117. [Google Scholar]
  24. Beltran, B.; Ahmed-Ali, T.; Benbouzid, M.E. High-order sliding-mode control of variable-speed wind turbines. IEEE Trans. Ind. Electron. 2009, 56, 3314–3321. [Google Scholar] [CrossRef]
  25. Leon, A.E.; Farias, M.F.; Battaiotto, P.E.; Solsona, J.A.; Valla, M.I. Control strategy of a DVR to improve stability in wind farms using squirrel-cage induction generators. IEEE Trans. Power Syst. 2011, 26, 1609–1617. [Google Scholar] [CrossRef]
  26. Senjyu, T.; Ochi, Y.; Muhando, E.; Urasaki, N.; Sekine, H. Speed and position sensor-less maximum power point tracking control for wind generation system with squirrel cage induction generator. In Proceedings of the 2006 IEEE PES Power Systems Conference and Exposition, Atlanta, GA, USA, 29 October–1 November 2006; pp. 2038–2043. [Google Scholar]
  27. Senjyu, T.; Ochi, Y.; Kikunaga, Y.; Tokudome, M.; Yona, A.; Muhando, E.B.; Urasaki, N.; Funabashi, T. Sensor-less maximum power point tracking control for wind generation system with squirrel cage induction generator. Renew. Energy 2009, 34, 994–999. [Google Scholar] [CrossRef]
  28. Faiz, J. Design and Implementation of a Solid-state Controller for Regulation of Output Voltage of a Wind-driven Self-excited three-phase Squirrel-cage Induction Generator. In Proceedings of the International Conference on Electrical Machines and Systems, Nanjing, China, 27–29 September 2005. [Google Scholar]
  29. Duong, M.Q.; Grimaccia, F.; Leva, S.; Mussetta, M.; Ogliari, E. Pitch angle control using hybrid controller for all operating regions of SCIG wind turbine system. Renew. Energy 2014, 70, 197–203. [Google Scholar] [CrossRef]
  30. Duong, M.Q.; Grimaccia, F.; Leva, S.; Mussetta, M.; Le, K.H. Improving transient stability in a grid-connected squirrel-cage induction generator wind turbine system using a fuzzy logic controller. Energies 2015, 8, 6328–6349. [Google Scholar] [CrossRef]
  31. Xu, P.; Shi, K.; Bu, F.; Zhao, D.; Fang, Z.; Liu, R.; Zhu, Y. A Vertical-Axis Off-Grid Squirrel-Cage Induction Generator Wind Power System. Energies 2016, 9, 822. [Google Scholar] [CrossRef]
  32. Pena, R.; Clare, J.C.; Asher, G.M. Doubly fed induction generator using back-to-back PWM converters and its application to variable-speed wind-energy generation. IEE Proc.-Electr. Power Appl. 1996, 143, 231–241. [Google Scholar] [CrossRef]
  33. Alrifai, M.; Zribi, M.; Rayan, M. Feedback Linearization Controller for a Wind Energy Power System. Energies 2016, 9, 771. [Google Scholar] [CrossRef]
Figure 1. A block diagram representation of the squirrel cage induction generator (SCIG)-based wind energy conversion (WEC) system.
Figure 1. A block diagram representation of the squirrel cage induction generator (SCIG)-based wind energy conversion (WEC) system.
Energies 10 00604 g001
Figure 2. The power coefficient C p versus the tip speed ratio λ for several values of the pitch angle β .
Figure 2. The power coefficient C p versus the tip speed ratio λ for several values of the pitch angle β .
Energies 10 00604 g002
Figure 3. The mechanical power P m versus the generator speed ω m for several values of the wind speed V.
Figure 3. The mechanical power P m versus the generator speed ω m for several values of the wind speed V.
Energies 10 00604 g003
Figure 4. A block diagram representation of the controlled SCIG-based WEC.
Figure 4. A block diagram representation of the controlled SCIG-based WEC.
Energies 10 00604 g004
Figure 5. The profile of the wind speed versus time.
Figure 5. The profile of the wind speed versus time.
Energies 10 00604 g005
Figure 6. The currents i d s and i d s * versus time.
Figure 6. The currents i d s and i d s * versus time.
Energies 10 00604 g006
Figure 7. The currents i q s and i q s * versus time.
Figure 7. The currents i q s and i q s * versus time.
Energies 10 00604 g007
Figure 8. The fluxes ψ d r and ψ d r * versus time.
Figure 8. The fluxes ψ d r and ψ d r * versus time.
Energies 10 00604 g008
Figure 9. The angular speeds ω m and ω m * versus time.
Figure 9. The angular speeds ω m and ω m * versus time.
Energies 10 00604 g009
Figure 10. The power coefficient C p versus time.
Figure 10. The power coefficient C p versus time.
Energies 10 00604 g010
Figure 11. The electromagnetic torque T e m versus time.
Figure 11. The electromagnetic torque T e m versus time.
Energies 10 00604 g011
Figure 12. The active power of the stator P s versus time.
Figure 12. The active power of the stator P s versus time.
Energies 10 00604 g012
Figure 13. The currents i d g and i d g * versus time.
Figure 13. The currents i d g and i d g * versus time.
Energies 10 00604 g013
Figure 14. The currents i q g and i q g * versus time.
Figure 14. The currents i q g and i q g * versus time.
Energies 10 00604 g014
Figure 15. The voltages U d c and U d c * versus time.
Figure 15. The voltages U d c and U d c * versus time.
Energies 10 00604 g015
Figure 16. The output active power of the grid P g versus time (negative sign for produced power).
Figure 16. The output active power of the grid P g versus time (negative sign for produced power).
Energies 10 00604 g016
Figure 17. The output reactive power of the grid Q g versus time.
Figure 17. The output reactive power of the grid Q g versus time.
Energies 10 00604 g017
Table 1. The parameters of the SCIG and the wind turbine.
Table 1. The parameters of the SCIG and the wind turbine.
ParameterValue
The rated power300 KW
The rated apparent power300 KVA
The rated voltage (line to line)575 v
The frequency/angular speed 2 π 50 rad/s
The nominal frequency of the system50 Hz
The resistance of the stator 0.0063 Ω
The leakage inductance of the stator11.8 mH
The resistance of the rotor 0.0048 Ω
The leakage inductance of the rotor11.6 mH
The magnetizing inductance11.6 mH
The capacitor of the DC bus20 mF
The rated voltage of the DC bus760 v
The grid resistance 0.1 Ω
The grid leakage inductance0.6 mH
The inertia of the generator10.0 kg m 2
The number of pole pairs2
The rated generator speed158.7 rad/s
The rotor diameter of the wind turbine14 m
The density of air1.22 kg/m 3
The rated wind speed12 m/s
The rated rotor speed19.7 rpm
The gearbox ratio23
The inertia of the turbine50.0 kg m 2

Share and Cite

MDPI and ACS Style

Zribi, M.; Alrifai, M.; Rayan, M. Sliding Mode Control of a Variable- Speed Wind Energy Conversion System Using a Squirrel Cage Induction Generator. Energies 2017, 10, 604. https://doi.org/10.3390/en10050604

AMA Style

Zribi M, Alrifai M, Rayan M. Sliding Mode Control of a Variable- Speed Wind Energy Conversion System Using a Squirrel Cage Induction Generator. Energies. 2017; 10(5):604. https://doi.org/10.3390/en10050604

Chicago/Turabian Style

Zribi, Mohamed, Muthana Alrifai, and Mohamed Rayan. 2017. "Sliding Mode Control of a Variable- Speed Wind Energy Conversion System Using a Squirrel Cage Induction Generator" Energies 10, no. 5: 604. https://doi.org/10.3390/en10050604

APA Style

Zribi, M., Alrifai, M., & Rayan, M. (2017). Sliding Mode Control of a Variable- Speed Wind Energy Conversion System Using a Squirrel Cage Induction Generator. Energies, 10(5), 604. https://doi.org/10.3390/en10050604

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop