Introdução Variação
Introdução Variação
Introdução Variação
Mohammed. Ouassaid,
multivariable
with
many
uncertain
factors
and
to
achieve
Simulation
Maximum
results,
by
Power
Point
Matlab/Simulink,
Tracking
have
(MPPT).
shown
the
WEC;
MPPT;
I.
INTRODUCTION
MODELING O F WECS
II.
(6)
A.
1
2
PTurbine = - p1fR Cp(A,j3)v3
2
(1)
we = pnwm where wm
is the
mechanical angular velocity and Pn is the number of pole
pairs of the PMSG. For salient pole PMSG, Ld = Lq = L, and
(3)
1',; = 2L,
_ _Pn lv, Iv.f sineS)
(7)
where S and V, are the load angle torque and the stator flux,
respectively. Therefore, it's possible that the torque of the
generator is controllable by regulating the torque angle S if
the amplitude of stator flux linkage Vs is constant. So, if the
stator resistance is ignored then voltage vector applied to the
machine remains constant during one period of the control
algorithm: T,. As a result:
V," (k+I) = V," (k)+ vJ ,
(8)
where V," (k) and V, (k+l) are the stator flux vectors at the
sampling times k T, and (k+I) T, , respectively. Then, the
variation of the stator flux vector is [II]:
(9)
Since, according to (9), !1V, is achieved through the
Ai
0.035
A+ 0.0813
133 +I
(4)
Cp(A, 13) , that is CP,n" = 0.41, is achieved for Aopl = 8.1 and
(10)
where F is the viscous friction coefficient, J is the total inertia
and n is the mechanical torque developed by the turbine.
12000
10000
"-
8000
6000
4000
2000
0
0
50
100
Wind generator speed(rd/s)
150
III.
A.
vAoptG
OJm-opt = -- R-
(11)
ROJm-opt
Aopt
]3
(12)
by
(15)
So, to obtain commutation around the surface and good
dynamic performances, the control includes two terms [13]:
(16)
where ueq is valid only in the sliding surface. During the
sliding mode and in permanent regime, ueq is calculated from
(15), while un is used so as to guarantee the attractiveness of
the variable to be controlled towards the commutation surface.
Then
(17)
where kO)>- O. Moreover, nonlinear control Sliding Mode is a
discontinuous control. In order to reduce the chattering, the
continuous function as exposed in (18) where sgn(SO)) is a
sign function defined as [13]:
power.
B.
P_ra_te_d
> _.
Pg
Fig.2.
WECS
-------'
,gn(s)
-1
SO)>-t:
t: ISO)I
(18)
- t:>-SO)
(21)
where Ude is the DC-bus voltage and (SI ,S2' S3) the
switching states of the rectifier. FigA illustrates the voltage
vectors obtained this way [11]. Accordingly, the component of
the stator flux linkage, in the two phase stationary axis is given:
(22)
(23)
1f/1/3
(29)
Consequentely, reactive and active power control can be
achieved by controlling quadrature and direct current
components, respectively. The DC-side equation can be given
by:
(30)
The flux amplitude and its phase are derived respectively from:
(24)
/3
0, = arctg(If/s )
(25)
If/\a
and the position of stator flux linkage space vector (0, ). That's
illustrated in Fig.3. Therefore, the vector of voltage makes the
flux rotate and produce the desired torque. Also, with
comparators, the amplitude of the flux is kept in a pre-defmed
band.
C.
current
respectively.
The
reference
iq r-f for
the iq-f current is derived from the desired power factor. So:
(26)
(27)
didr-f
.
-_-
dt
(31)
(32)
SIMULATION RESULTS
fp
'w.
MPPT
Prated CONTROL
------+
(1Jm
roled
8,
Flux
&
Torque
calculation
C C1
Rj;Lj
SMC,
Q,.
3V
TABLE I .
C'I'
CT
C =O
CT=I
CT=O
CT= -I
CT=I
C'I'=I
CrO
CT= -I
81
SWITCHING VECTORS
82
V
Vo
V,
V2
V7
V6
83
V4
V7
V6
V3
Vo
VI
84
V6
V7
V2
V,
Vo
V3
V,
Vo
VI
V4
V7
V2
86
8,
VI
Vo
V3
V6
V7
V4
V2
V7
V4
VI
Vo
V,
0.5,----_-_--_-_
0.4
0.2
0.1
0:c------;--;2O---;--
is decreasing
l r
U
0.3-
2
Time(s)
Time(s)
1.5
1.5
0.5
0.5
-0.5
,O:c------:--;2O---::--:
Fig.6. Pitch angle fJ (in degree
Time(s)
0.5
1.5
2
2.5
Time(s)
3.5
120
.00
80
TABLE II.
03
02
- Generator speed
- Optimal speed
01
60
<I
40
20
0------2--7---
Time(s)
(2)
-0.2
.0.1
0.1
0.2
OJm
OJ
Voltage
-- Current
-200
_ 400 L--L-----L-----
0.42
0.44
0.46
0.48
0.5
Q)
0.635 m(H)
0.192(Wb)
pole pairs
[2]
M.A . Abdullah, AH.M. Yatim, Chee Wei Tan, " A study of maximum
power point tracking algorithms for wind energy system, " Proceedings
of Clean Energy and Technology (CET), IEEE First Conference on, pp.
321 - 326,201l.
[3]
[4]
F.O. Bennouna, S.N. Heraud, " Diagnosis & fault detection in wind
energy conversion system, " 10th IEEE International Conference on
Environment and Electrical Engineering (EEEIC),pp. 1-4,May 2011
[5]
[6]
[7]
[8]
[9]
CONCLUSIONS
0.05(
Time(s)
100 (rd/s)
[I]
V.
REFERENCES
10 (kW)
Pn
II
200
Value
Pr rated power
R stator resistance
of PMSG
400
Parameter