Adaptive Super-Twisting Sliding Mode Control For Wind Energy Conversion System
Adaptive Super-Twisting Sliding Mode Control For Wind Energy Conversion System
Adaptive Super-Twisting Sliding Mode Control For Wind Energy Conversion System
3524-3532
© Research India Publications. http://www.ripublication.com
3524
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
J J r ng2 J g
with K K r ng2 K g
Figure 1. WECS based on DFIG T n T
g g em
3525
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
P Vsd I sd Vsq I sq
M
Q V I V I , Tem n p
Ls
I rd sq I rqsd (8)
sq sd sd sq
and
1 M
I rd Vrd Rr I rd ss Lr I rq s
Lr Ls
(10)
1 V R I s L I s M
I rq L rq r rq s r rd s
Ls s
r
M2
where, 1 , n p is the pole pair number and s is the Figure 5. System to be controlled
Ls Lr
slip speed.
Hence, the electromagnetic torque becomes Adaptive Super-Twisting sliding mode control
3526
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
(ASTW) control strategy is proposed to improve the Vrd controls the reactive power ( I rd ) and Vrq controls the
performance and effectiveness of the WT-DFIG system.
active power ( Tem ).
The derivative of system (14) is
Assumption 1. We assume that G1 and G2 are bounded
1 M functions which satisfy the following conditions:
S1 Vrd Rr I rd ss Lr I rq s I rd ref
Lr Ls and G2 (t , S2 ) 2 S2
G1 (t , S1 ) 1 S1
12 12
, where 1 and 2
(15)
S n M V R I s L I s M T are some positive constants.
2 p s rq r rq s r rd s s ref
Ls Lr Ls
The theorem below which is inspired from [17], provides
sufficient conditions to ensure the stability and robustness of
Let’s consider functions G1 and G2 as follows:
current and electromagnetic torque control of the DFIG.
1 M Theorem 1: Consider the system described by equation (10)
G1 Rr I rd ss Lr I rq s I rd_ref
Lr Ls and the sliding surface given by (14). If we apply the control
(16) law (18) based on Second Order Sliding Mode Controller to the
G n M R I s L I s M T system. The closed-loop system (15) is ensured with finite-time
2 p s r rq s r rd s s ref
Ls Lr Ls convergence with the following adaptive gains
(ai ,bi for i 1, 2)
Assumption 1: We assume that functions G1 , G2 are bounded
G1 1 S1
12
, G2 2 S2
12
ki i Si , if Si i
where 1 , 2 are positive unknown constants. a
i (t , S i ) 2
0, if Si i
(20)
Then, equation (15) can be written as: b (t , S ) 2 a (t , S ) 4 2
i i i i i i i
1
S1 L Vrd G1 where i , ki ,i , i and i are arbitrary positive constants. The
r
(17) parameter i is introduced in order to get only positive values
S2 n p M sVrq G2
Ls Lr for ai (t , Si ) 0 .
y2 b2 ( S2 , t ) sgn S2
and
S qa S 1 2 sgn S + G
where, the adaptive gains are
1 1 1 1 * 1
* qb1 sgn S1
ai ( Si , Si , t ), (23)
, for i 1, 2 (19) 1
bi ( Si , Si , t )
q , * qy1
Lr
k1 1 S1 , if S1 1
a
1 (t , S1 ) 2
0, if S1 1
(24)
b (t , S ) 2 a (t , S ) 4 2
1 1 1 1 1 1 1
Figure 6. Block-scheme of the proposed WT-DFIG control We can easily deduce the derivative of (25)
strategy
3527
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
z1
1 qa1
1 1
z1 z2 G1 a1
Lr1 1 412 21 1 12 (35)
z1 2 2 2
(26) 1
z qb1 z
2 z1
1 From (33), we have
1 T 2
V0 z z Qz 1 zT z 1 z
2
Using assumption A1, we get (36)
2 z1 2 z1 z1
G1 ( x, t ) 1 ( x, t ) S1 sgn(S1 ) 1 ( x, t ) z1
12
(27)
We note min P and max P as the minimal and
where 1 ( x, t ) is a bounded function, 0 1 ( x, t ) 1 : maximal eigenvalues of matrix P , then we can write:
z A z1 , z2 z (28) min P z zT Pz max P z
2 2
T 1
1
Then, we can conclude
2 2
V z, a1 , b1 V0 a1 a1* b1 b1* (30)
2 1 2 1
V0 ( z ) rV01 2 (37)
where
min
12
( P)
where, r
V0 ( z ) 1 412 z12 z22 4 z1 z2 zT Pz (31) max ( P)
12
Q Q21 qb1 2q1a1 1 412 21 1. V1 V0
1
2 1
1
a1S1 2 b1S1 2
2 1
By choosing
From (39), we can rewrite (38) as
b1 21a1 1 412 (34)
3528
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
1 k
V z , a1 , b1 1V11 2 a1 a1 1 S1
1 2 1
(41)
1 k
b1 b1 1 S1
1 2 1
Thus, we obtain
V z, a1 , b1 1V1 z, a1 , b1
12
(42)
1
a1 k1 S1 (43)
2
From adaptation law (24) and (43), we have
1
b1 21a1 1 412 b1 21a1 k1 S1
2
k1 1
By selecting 1 , we ensure that 0 .
2k1 1
Hence, the derivative of system (30) is guaranteed definite Figure 8. Block diagram of the global system
negative and thus convergence of S1 0 which can be finally
written as follows: SIMULATION RESULTS AND DISCUSSION
V z, a1 , b1 1V11 2 (45) The proposed control strategy has been validated with
numerical validation using 1.5 MW 3-blade-WT FAST
with 1 0 . simulator, developed by National Renewable Energy
Laboratory (NREL), figure 8.
Similar development can be used to prove the stability
analysis in the case of Vrq . The wind inflow for the simulation generated by FAST
simulator is shown in figure 9. The wind speed profile
Simulation using FAST Code simulator represents a turbulent wind speed variation between 8 m/s and
14.5 m/s. Figure 10 shows the mechanical speed of the
The proposed Adaptive Super-Twisting algorithm has been generator until 1800 rpm. The desired DFIG reference are
tested for validation using FAST (Fatigue-Aerodynamics- directly tracked through an Adaptive Super-Twisting control.
Structures-Turbulence) software which is an effective aero- As can be clearly seen in figure 11 and figure 14, the
elastic simulator that can predict extreme and fatigue loads of a electromagnetic torque and rotor current of the DFIG show
spread of WT’s configurations (two and three-blade horizontal- good tracking performances under turbulent wind-speed
axis) [18]. FAST relies on advanced engineering models variation. The rotor current reference in figure 11 is determined
derived from basic laws with appropriate assumptions. This
as a function of stator voltage Vs . Figure 14 represents the
3529
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
3530
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
3531
International Journal of Applied Engineering Research ISSN 0973-4562 Volume 13, Number 6 (2018) pp. 3524-3532
© Research India Publications. http://www.ripublication.com
3532