Asynchronous Machine with Wind Turbine IRFO Control
A. ALAZRAG, L. SBITA
University of Gab´
es, National engineering school of Gab´
es (ENIG),
Process laboratory, Energetic, Environment and Electrical system, TUNISIA
Energy has become a fundamental element because of our
different demands in many domestic and industrial terms.
This prompts us to always think about new techniques for
producing this energy. Since the dawn of humanity, the pro-
duction of energy has been based largely on fuels such as
wood, fossils (coal, oil, gas. . . ), then uranium. But the
big problem accompanying the use of these materials is
the emission of gases and the massive releases of various
compounds; which unbalances our planet and pushes to-
wards an open non-renewable cycle. In just over a cen-
tury, energy with electricity as a modern form has taken
a prominent place. Its production covers a third of the
world’s energy consumption which is concentrated mainly
in thermo mechanical machines where combustion is on a
large scale with the direct emission of several million tons
of CO2 causing high degrees of pollution and temperature;
as well as the reduction of nature reserves. So, we need to
look for other alternatives to fossil fuels to produce elec-
tricity from renewable sources that are non-polluting and
more economical by making good use of the elements of
nature such as water, sun and wind. This is the goal of our
study, which focuses on one of the renewable energies in
development at the time, which is wind energy. The reli-
able and cost-effective wind turbine is the ideal source of
electricity for many applications. Wind turbines come in
many sizes, from small wind turbines of a few watts to
megawatt wind turbines feeding the electricity grid. The
largely dominant technology today is horizontal axis, three-
bladed and sometimes two-bladed windward rotor turbines.
These wind turbines have a nominal power between 5KW
and 5MW, they can operate at fixed speed or at variable
speed [12][13]. The types of generators associated with
wind turbines are asynchronous machines and synchronous
machines in their different variants. Among all the renew-
able energies contributing to the production of electricity,
wind energy currently holds the star of renewable for re-
gions and countries with enormous wind potential. It is
one of the most promising, in terms of ecology, competi-
tiveness, scope and creation of jobs and wealth. We will
focus on the technological advances that have allowed the
construction and proper functioning of wind turbines and
their integration into electricity production. In order to bet-
ter exploit wind resources for different wind conditions, this
study focuses on the Asynchronous Machine with indirect
orientation of the rotor flux which is the heart of a large part
of current wind turbines due to its advantages relative to
other electromagnetic actuators. Dynamic management of
interconnected grid-connected wind farm storage systems
would therefore be necessary to ensure the mitigation of
variable wind farm output variability and to maintain grid
power stability [29]. Based on the state of charge (SOC)
of the BESS authors proposed an algorithm for dynamic
programming in order to smoothly fluctuate wind power
[30]. Different application of strategies have been applied
for wind farm active power control, voltage control. SCIG
Abstract: The paper deals with a squirrel cage induction generator connected to the grid through a back-to-back converter
driven by vector control. The stator-side converter controls the generator torque by means of an indirect vector control
scheme. In order to reduce the system dependence from the mechanical system behavior, a torque loop is used in the current
reference calculations. The battery energy storage system (BESS) plays a fundamental role in controlling and improving
the efficiency of renewable energy sources. Stochasticity of wind speed and reliability of the main system components are
considered. The grid-side converter controls the DC bus voltage and the reactive power in order to accomplish the grid
codes. Speed control using flow directional control, indirect conventionally uses proportional integral (PI) type current
regulators, which achieve satisfactory objectives on torque and flow dynamics. The objective of this article is to present an
indirect vector control strategy with oriented rotor flux using current regulators of the proportional integral (PI) type,
applied to an asynchronous machine supplied by a voltage inverter, capable of supplying during restrictive stresses, more
satisfactory torque and flux responses. The obtained simulation results upon simulation tests of the global system are
developed under the MATLAB / Simulink environment and are satisfactory .
Keywords: Renewable Energy,Wind Turbine, Asynchronous machine (SCIG),Battery storage ,Indirect Rotor flux
orientation(IRFO),PI, electromagnetic torque, axis d,q.
Received: April 29, 2021. Revised: April 25, 2022. Accepted: May 29, 2022. Published: June 30, 2022.
,QWURGXFWLRQ
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
141
Volume 17, 2022
control can be implemented using different approaches :
scalar or vector control, direct or indirect field orientation,
rotor or stator field orientation [12][13]. Scalar control [15]
is simple to implement, but easily unstable. A better per-
formance is obtained with direct vector control, requiring
sensed flux values to define and control the field orientation
references. This, however means that it is necessary to use
hall-effect sensors, which, inpractise, is problematic, and
expensive,[14][16]. The indirect field-orientation method is
more sensitive to the machine parameters but removes the
necessity of direct flux sensing[15][16]. This paper pro-
poses an indirect vector control strategy less sensitive from
the machine parameters than the conventional scheme[8].
The wind turbines are also equipped with a control system
based on electronic converters to adapt to wind conditions.
This work is organized as follows: The Section II presents a
modeling of the wind turbine connected to the grid as well
as a modeling in a two-phase reference frame (the adapted
model is based on the Park transformation) linked to the ro-
tating field with a view to supply and control by static con-
verters. of the double feed asynchronous machine. The Sec-
tion III,study of storage system. The Section V implements
the indirect vector control of the squirrel cage asynchronous
generator based on a PI (Proportional Integral) regulator,
Section IV the control ensures the decoupling of the d and q
axes, the purpose of which is to improve the behavior static
and dynamic system. The Section VII, is reserved for sim-
ulation results with interpretation of these results.
Figure 1: Wind turbine with SCIG connected to electrical
grid
The system being analyzed can be seen in Figure 1. It is
important to point out that there are more elements, e.g. a
transformer and a wind farm grid, which are not included
in the present paper. The squirrel cage induction generator
(SCIG) is attached to the wind turbine by means of a gear-
box. The SCIG stator windings are connected to a back to
back full power converter.
A wind energy conversion system (WECS) transforms wind
kinetic energy to mechanical energy by using rotor blades.
This energy is then transformed into electric energy by a
generator, so the turbine is one of the most important el-
ements in wind turbine. In order to better understand the
process of wind energy conversion, descriptions of the ma-
jor parts of a wind turbine are given in this section. Accord-
ing to Newton’s law, the kinetic energy for the wind with
particular wind speed is described as:
Em=1
2m.V2
v(II.1)
Where m represents the mass of the wind, and its power can
be written as :
P
m=1
2.Cp(λ,β).ρ.S.V3
v(II.2)
Where : S : Blade swept area (m2), ρ: Specific density of
air(Kg/m3)
Vv: wind speed (m/s), Cp(λ,β): Power coefficient, λ:
speed ratio, β: pitch angle With :
λ=R.
Vv
(II.3)
:Mechanical turbine speed (rad/s), R : radius of the tur-
bine blade (m).
We can derive the formula of the torque which is :
Tt=1
2.t
.Cp(λ,β).ρ.π.R2.v3(II.4)
Where Cpis the power coefficient of the turbine which is
obtained from the following equation :
Cp(λ,β) = K1.(K2
λK3.βK4)e
K5
λ+K6.λ(II.5)
With :
1
λ=1
(λ+0.08.β)0.035
β3+1(II.6)
Where : k1=0.5176,k2=116,k3=0.4,k4=5,k5=
21andk6=0.006
This Figure shows the evolution of the power coefficient as
a function of λfor different values of β. The maximum is
obtained for an optimum pitch angle βopt =0 and a relative
speed λopt =0 .
0RGHOLQJRIDZLQGWXUELQHZLWK
LQGXFWLRQJHQHUDWRU
:LQGWXUELQH
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
142
Volume 17, 2022
Figure 2: Evolution of the power coefficient according to λ
and β
With the dependence on the λand β, maximum value of
Cp could be reached and maintained through controlling the
pitch angle and generator speed at particular wind speed. A
group of typical Cp λcurves for different βis shown in
Figure 2 and there is always a maximum value for Cp at
one particular wind speed. It is found that the power coef-
ficient is a function of λand β, and in order to reach the
maximum energy extracted from the wind turbine. The fol-
lowing equations show the values to be considered:
Topt =1
2.Cpmax(λ).ρ.S.R3
λ3
opt
.2(II.7)
re f =λopt
R.Vv(II.8)
The gear box provides speed and torque conversions from a
rotating power source to another device, using gear ratios.
A gearbox is often used in a wind turbine to increase the ro-
tational speed from a low-speed main shaft to a high-speed
shaft connecting with an electrical generator. A gear box is
mathematically modeled by the following equations:
(g=Kg.T
Tg=1
Kg.TT(II.9)
With : Kgis the speed gear box gain Kg0 For our study,
we chose a wind turbine with direct coupling
The differential equation which makes it possible to deter-
mine the evolution of the mechanical speed from the total
mechanical torque (tm)is given by:
J.dm
dt =Tm(II.10)
Where: J is the total inertia brought back to the generator
shaft, including the inertia of the turbine, the gearbox and
the generator. The mechanical torque deduced from this
coupling is the sum of all the torques applied to the rotor:
Tm=TgTem fg(II.11)
With:
Tm: is Electromagnetic torque developed by the generator,
Tg: gear box torque,
f: Total coefficient of friction of the mechanical coupling.
The conversion of mechanical energy to electric energy is
performed by the turbine and the generator. Different gen-
erator types have been used in wind energy systems. These
include the squirrel cage induction generator (SCIG), dou-
bly fed induction generator (DFIG), and synchronous gen-
erator (SG) [1]. The SCIG is simple and rugged in con-
struction. It is relatively inexpensive and requires minimum
maintenance. The SCIGs are also employed in variable-
speed wind energy systems. To date, the largest SCIG wind
energy systems are around 4 MW in offshore wind farms.
The voltage equations for the stator and rotor of the gener-
ator in the reference of Park with two axes (d,q), [3], [5].
Figure 3: Representation of fictitious windings of d-q axes
For the machine equations, assuming that the stator and
rotor windings are sinusoidal and symmetrical [14],[16], the
relation between voltage and currents on a synchronous ref-
erence dq can be written as:
vsd
vsq =Rs0
0Rs isd
isq +d
dt φsd
φsq +0ωs
ωs0 φsd
φsq
(II.12)
0
0=Rr0
0Rr ird
irq +d
dt φrd
φrq +0ωr
ωr0 φrd
φrq
(II.13)
2.1.1 Gearbox
2.1.2 Driveshaft dynamic equation
2.2 Asynchronous machine (SCIG)
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
143
Volume 17, 2022
Figure 4: Back to Back converter
φsd
φsq =LsM
MLr isd
isq (II.14)
And
φrd
φrq =LsM
MLr ird
irq (II.15)
With : Vs: Stator voltage, (is,ir): stator and rotor
current,(phis,phir): stator and rotor flux, (Rs,Rr):stator
and rotor resistance, (Ls,Lr):stator and rotor inductance , M
: mutual inductance , d : d-axis, q : q-axis. We become to
replace the flux with their expressions according to the cur-
rents; in using (10) and (11), the model of the three-phase
asynchronous machine in the axis coordinate system (d,q)
will be given in matrix form :
vsd
vsq
0
0
=
Rs+Ls(d/dt)LsωsM(d/dt)Mωs
LsωsRs+Ls(d/dt)MωsM(d/dt)
M(d/dt)MωrRr+Lr(d/dt)Lrωr
MωrM(d/dt)LrωrRr+Lr(d/dt)
isd
isq
ird
irq
(II.16)
This leads to the equivalent diagram coupled with an asyn-
chronous machine along axis d.
Figure 5: Equivalent diagrams coupled of an asynchronous
machine according to the axis d
Also, the equivalent diagram coupled of an asyn-
chronous machine according to the axis q.
Figure 6: Equivalent diagrams coupled of an asynchronous
machine according to the axis q
With : ls =LsM,lr =LrMExpression of elec-
tromagnetic torque In the general case, the instantaneous
electrical power Pe supplied to the windings stator and ro-
tor is expressed as a function of the axis sizes d, q:
P
e=vsd .isd +vsq.isq (II.17)
It breaks down into three series of terms: *.Power dissipated
in Joule losses:
Rs(i2
sd +i2
sq) + Rr(i2
rd +i2
rq)(II.18)
**. Power transmitted to the rotor in the form of a variation
of magnetic energy :
isd (dφsd /dt) + isq(dφsq/dt) + ird (dφrd /dt) + irq(dφrq/dt)
(II.19)
***. Mechanicanical power
P
m= (isq.φsd isd φsq)ωs+ (irqφrd ird φrq)ωr(II.20)
The electromagnetic torque developed by the machine is
given by the expression:
Cem =P
m
and=ω
p(II.21)
Then the scalar expression of torque :
Cem =p.(φsd isq φsqisd )(II.22)
It is possible to obtain other expressions of the instanta-
neous torque by using the stator flux expressions:
Cem =p.M
Lr
.(φrd isq φrqisd )(II.23)
Also :
Cem =p.M.(ird isq irqisd )(II.24)
This very important relationship highlights the fact that the
torque results from the interaction of components of quadra-
ture stator and rotor currents[9],[10]. The active and reac-
tive power yields :
P=3
2.(vsd isd +vsqisq)(II.25)
And :
Q=3
2.(vsqisd vsd isq)(II.26)
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
144
Volume 17, 2022
Active and reactive power, provided by the grid-side con-
verter, can be expressed as:
P
g=3
2.(vgd igd +vgqigq)(II.27)
/
Qg=3
2.(vgqigd vgd igq)(II.28)
For the proper functioning of the hybrid energy system, the
storage system plays a crucial role, it allows for continuity
of service and better quality of energy supplied. We recall
some electrical parameters used to characterize a battery,
these are:
Nominal capacitor (Qn): This is the maximum number of
ampere-hours (Ah) that can be extracted from the battery,
for given discharge conditions.
The state of charge “SOC” (State Of Charge): This is
the ratio between the capacity at time q(t) and the nominal
capacity Qn, is :
SOC(t) = q(t)
Qn
,with(0SOC 1)(III.1)
If SOC=1 the battery is totally charged and if SOC=0 the
battery is totally discharged.
The charging cycle (or discharging): This is the parameter
which reflects the relationship between the nominal capac-
ity of a battery and the current at which it is charged (or
discharged). It is expressed in hours.
•Cycle life : This is the number of charge/discharge cycles
that the battery can sustain before losing 20% of its nominal
capacity[5].
Figure 7: Battery discharge curve
By analyzing the figure above, we can see the presence
of three specific points on the characteristic (Q-V): these
three points are: the full load voltage (E0), the voltage cor-
responding to the end of the exponential zone ( Eexp) and
the corresponding voltage at the end of the nominal zone
(En). The charge and discharge equations are given as fol-
lows [7]:
*Discharge :
EB=E0Ri KQ
Qit
(it+i) + Eexp(t)(III.2)
//
*Charge :
EB=E0Ri KQ
Qit
iKQ
Qit
it+Eexp(t)(III.3)
With : Eexp(t) = B.|i(t)|.(Eexp(t) + A.sel(t)) Figure 8
shows the discharge characteristic of the storage system
used and the evolution of its voltage for different discharge
currents.
Figure 8: Characteristic VB = f(t) for different discharge
currents.
The bidirectional DC DC converter is a combination
of boost and buck converters. Such a converter is used to
charge and discharge the battery.
Figure 9: Circuit of the DC–DC bidirectional converter
The boost mode is applied for the discharging procedure of
the battery storage. Figure shows the circuit of the boost
mode operation of the converter, where the direction of the
inductor current is from the lower voltage side to the high-
erK1 voltage side . The averaged large signal inductor cur-
rent, ıL, and the DC-bus output voltage, Vdc, in a contin-
uous conduction mode (CCM) of operation can be found
using the equations below. is closed (K10=1) and is open
(K11=0).
LdiB
dt =VBVdc (III.4)
3. Battery storage
3.1 %oost 0ode
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
145
Volume 17, 2022
The buck mode is applied for the charging process of the
battery storage. Figure present the circuit of the buck mode
operation converter. In contrast to the buck mode operation,
the inductor current flows from the higher voltage side to
the lower voltage side . The averaged large signal inductor
current, ıL, and the output battery voltage, VB, are calcu-
lated by the equations below, and describe the buck-mode
operation in a CCM of the converter. is open (K10=0) and
is closed (K11=1)
LdiB
dt =VB(III.5)
By analyzing these two configurations, we can conclude
that the relationships between the input quantities (VB,iB)
and the output quantities (Vdc,iBm)of the converter are
given by the system of equations below
LdiB
dt =VBVdc(1u11)
iBm=IB(1u11)(III.6)
The three-phase voltage inverter allows the exchange of en-
ergy between a DC voltage source and a three-phase induc-
tive load [2].
Figure 10: Simplified model of the DC/AC converter
The three-phase system is considered balanced so we
can write:
Va+Vb+Vc=0 (IV.1)
The corresponding PWM gate signals for the converter. To
a 2 level three-phase voltage source inverter, there are six
switches of three legs in inverter controlling the phase volt-
age and thus the current of induction generator. By defining
the ON and OFF states of upper switch by 1 and 0, respec-
tively, for one leg, there exist up to eight different states for
inverter outputs. They are summarized in table 1 as well as
the resulted phase voltage in ABC and αβ frames. Eight
inverter output voltages can be considered as eight voltage
vectors [0, 0, 0] through [1, 1, 1].
Table 1: the phase-to-phase and phase to phase voltages of
the inverter according to the states of the switches
The composed voltages delivered by the two level three
phase inverter as a function of the state of the IGBTs are
given by the following equations:
UAB =Vdc(SASB)
UBC =Vdc(SBSC)
UCA =Vdc(SCSA)
(IV.2)
Knowing that composed voltages are expressed as a func-
tion of phase-to-neutral voltages by the following equa-
tions:
UAB =VAVB
UBC =VBVC
UCA =VCVA
(IV.3)
Then :
VA=Vdc
3(2SA(SB+SC))
VB=Vdc
3(2SB(SA+SC))
VC=Vdc
3(2SC(SB+SA))
(IV.4)
The integrated filter is of the low pass type, its purpose is to
improve the quality of the signals exchanged between the
converter and the grid parameter connection (Rf,Lf). Its
cut frequency (fc)is given by relation below:
fc=Rf
2.π.Lf
(IV.5)
Par application de la transform´
ee de Laplace, les courants
`
a la sortie du filtre (Rf,Łf)seront alors exprim´
es comme
suit:
Idgrid =Vdcr Vdgrid +ω.Lf.Iqgrid
Rf+p.Lf
Iqgrid =Vqcr Vqgrid ω.Lf.Idgrid
Rf+p.Lf
(IV.6)
Numerically , we write :
Idcr (n) = 1
2.π.fc.Te+1.Idcr (n1) + 1
Lf.Vdcr (n)+ω.Lf.Iqcr(n)Vdgrid (n)
2.π.fc+1
Te
Iqcr(n) = 1
2.π.fc.Te+1.Iqcr(n1) + 1
Lf.Vqcr(n)ω.Lf.Idcr (n)Vqgrid (n)
2.π.fc+1
Te
(IV.7)
Where Te is the sampling period.
3.2 Buck 0ode 2peration
4. Inverter modeling
4.1 Filter technology
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
146
Volume 17, 2022
Figure 11: General diagram of the oriented rotor flux vec-
tor control
The control vector orientation of the rotor flux is the most
used, because it eliminates the influence of the rotor and
stator leakage reactance and gives better results than the
methods based on the orientation of the stator flux or air
gap [1], [3], while ensuring the best torque behavior as a
function of the slip speed in steady state. According to the
general vector control scheme with oriented rotor flux is as
follows: The command is carried out by orienting the rotor
flux along the direct axis ”d” of the rotating frame is :
Isd :it is the excitation current generates and controls the
excitation flux φr.
Isq: it is the armature current, which at a given excitation
flux controls the torque. We impose that : Φrd =ΦretΦrq =
0 then :
Ce=pLm
Lr
Φrisq (V.1)
These expressions show that the flux depends only on the
direct component of the stator current ids, and that if the
latter is kept constant, the torque will only depend on the
quadrature component of the stator current iqs. However,
in the case of a voltage supply vsd and vsq influence both
ids and iqs, therefore on the flux and the torque, hence the
interest of adding compensation terms in order to make the
axes d and q completely independent.
The principle of this method consists in not using the mag-
nitude of the rotor flux but simply its position calculated
with reference sizes. In this method we do not need a sen-
sor, an estimator or an observer of flux [13], [15], [16], [17].
therefore we have no knowledge of the modulus and phase
of the rotor flux, this requires a measurement of the posi-
tion of the rotor. Figure represents the block diagram of an
indirect vector control with oriented rotor flux of an asyn-
chronous machine.
Figure 12: Block diagram of indirect rotor
The equation of a MAS used to define the transfer func-
tions required for the synthesis of regulators of speed and
stator currents. The different regulators used in the regula-
tion loops will be of the type proportional and integral (PI),
because the quantities to be regulated are continuous quan-
tities. We therefore have 3 regulators diagram :
: It takes the reference speed and the measured speed as
input. It acts on the torque (that is to say its output is the
reference torque) to regulate the speed.
Figure 13: functional diagram of speed control
(s)
re f (s)=(1
J.s+f)(Kp1+Ki1
s)
1+ ( 1
J.s+f)(Kp1+Ki1
s)=1+τ1.s
1+ (τ1+f
Ki1).s+J
Ki1.s2
(VI.1)
With : τ1=Kp1
Ki1This transfer function has a 2nd order
dynamic. Denominator in the canonical form, we have to
5. Oriented rotor flux control
5.1 Indirect rotor flux orientation
6. Grid side converter
6.1 Control strategy
6.1.1 The speed regulator
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
147
Volume 17, 2022
solve the following system of equations:
(J
Ki1=1
ω2
0
2ξ
ω0=τ1+f
Ki1
(VI.2)
For critical damping ξ=1, we obtained :
(Kp1=τ1.Ki1
Ki1=4.J
τ2
1
(VI.3)
: To ensure that the actual currents follow the set point cur-
rents, regulators currents acting on the control voltages are
essential, we are interested in the sizing of regulators. So
the two regulators are identical. The current regulation loop
respectively Isd and Isq can be represented by the figure :
Figure 14: functional diagram of Isd current regulation
Figure 15: functional diagram of Isq current regulation
With :
K2=1
Rs+ ( Lm
Lr)2.Rr
(VI.4)
τ=σ.L
Rs+ ( Lm
Lr)2.Rr
(VI.5)
The closed loop transfer function will therefore be :
Isd (s)
Isdre f (s)=(1
1+τ.s)(Kp2+Ki2
s)
1+ ( 1
1+τ.s)(Kp2+Ki2
s)=(Kp2.K2
τ).s+Ki2.K2
τ
s2+ ( Kp2.K2+1
τ) + Ki2.K2
τ
(VI.6)
The characteristic closed loop equation looks like this:
s2+2ξ ω0s+ω2=0 (VI.7)
For identification :
(Kp2=2ξ ω0τ1
K2
Ki2=ω2
0
K2
(VI.8)
The control of the electrical system of the grid is config-
ured in order to control the DC bus voltage and the reactive
power which is consumed or supplied by the converter of
the grid side. The DC bus voltage reference and the grid
voltage level are used to determine the current references
which determine the voltages to be applied on the grid side.
figure 16 depicts the block diagram of grid side control.
Figure 16: Block diagram of grid side converter control
To adjust the voltage Vdc, we can act on the active
power Pg, and its reference value is obtained by the ref-
erence value of the power dcwhich is given by:
P
dcre f =Vdc.Idcre f (VI.9)
By applying the Laplace transformation , we obtain:
Vdc
Idc
=1
C.s(VI.10)
Figure 17: functional diagram of DC bus voltage
The PLL is used to provide a unity power factor operation
which involves synchronization of the inverter output cur-
rent with the grid voltage and to give a clean sinusoidal
current reference [14]. The PI controller parameters of the
PLL structure are calculated in such a way that we can set
directly the settling time and the damping factor of this PLL
structure. The principle of the PLL is depicted in figure 18.
6.1.2 The current regulator
6.2 Grid-side converter
6.3 Phase locked loop
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
148
Volume 17, 2022
Figure 18: Block diagram of Phase locked loop
The PLL control strategy requires that the voltage be
equal to zero, which allows the decoupling of the active and
reactive power as presented by the equation below, [9],[10].
P
g=3
2.Vgd Igd (VI.11)
Qg=3
2.Vgd Igq (VI.12)
In addition, the relationship between the active power ex-
changed between the network and the DC bus is given by
Vdc.Idc =3
2.Vgd .Igd (VI.13)
The wind farm is based to asyncronous generator charac-
terised by :
Nominal power : Pn=2500000 W;
Nominal voltage : Vn=690V;
Table 2 : Characteristic of asynchronous machine .
The vector control performance with the orientation of
the indirect rotor flux of the asynchronous machine are eval-
uated using the software MATLAB / Simulink.
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 Tur-
bines), the speed ratio is equal to 1; for wind turbines with
a single blade, λis about 11. The three-bladed wind tur-
bines, as in our study, have a speed ratio of 6 to 7. The
speed ratio of Savonius wind turbines is less than 1 [8]. Cp
is the power coefficient or aerodynamic transfer efficiency
that varies with the wind speed.
Figure 19: Wind speed evolution
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 max-
imum limit of the power coefficient Cpmax= 16/27= 0.59.
Taking into account losses, wind turbines never operate at
this maximum limit, and the best-performing wind turbines
have a Cpbetween 0.35 and 0.45. Cpis specific to each
wind turbine, and its expression is given by the wind tur-
bine manufacturer or using nonlinear formulas. To calcu-
late the coefficient Cp, different numerical approximations
have been proposed in the literature.
Figure 20: Speed ratio evolution
7. Result of Simulation
7.1 Wind parameter
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
149
Volume 17, 2022
Figure 21: Power coefficient evolution
The wind speed varies between [6 m/s] and [11 m/s].
Figure 20 presents the aerodynamic power delivered by the
wind turbine and it reached 3kw when the wind speed is
up to 11 [m/s]. Figures 22 and 21 illustrate respectively the
mechanical speed of the generator shaft and the speed of the
turbine shaft. It can be noticed from during the simulation
time (150 seconds), the generator operates in both hypo and
hyper synchronous operating modes. The Figures 19and 20
show respectively the variation of the speed ratio λand the
variation of the power coefficient Cpand, which coincide
with the optimal speed ratio and with the maximum power
coefficient.
Figure 22: Pitch angle evolution
Figure 23: Mechanical speed evolution
The system has been exposed to a real wind speed profile, in
order to observe the behavior of the SCIG and its control. In
Figure 24 it is possible to observe how the system evolves
guided by the reference value. It is important to observe
how the electromagnetic torque instantly follows the refer-
ence value. That means, it is possible to choose the opera-
tional point or follow the wind variation. Thanks to this fast
control, the mechanical system, mainly the turbine, is able
to quickly reach a steady state. However, the turbine speed
always decreases due to the wind speed variation, despite
the wind speed variation appearing filtered on the turbine
by its own inertia. In Fig. 6(d), the Active Power value evo-
lution is represented, which presents almost the same shape
(just a little damped) as the wind profile. This damping is
due to the multiplication of the mechanical rotation speed
and the generator torque. As the generator torque almost
mimics the wind shape, the damping is due to the turbine
speed and its inertia, as explained before.
Figure 24: Electromagnetic Torque evolution
Normal PV production and average battery charge level.
This scenario is characterized by a high demand for grid
power to satisfy the load despite the batteries being dis-
charged. During PV production, the system supplies to the
grid and the batteries discharge.
7.2 speed variation
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
150
Volume 17, 2022
Figure 25: Power wind evolution
Figure 26: Axis d stator current
Figure 27: Axis q stator current
Figure 28: DC bus voltage
Figure 29: State of charge
Figure 30: current of battery
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
151
Volume 17, 2022
Figure 31: oltagev of battery
Figure 32: Power of battery
Figure 33: Axis d rotor current
Figure 34: Axis q rotor current
Figure 35: Reactive and active power
This paper presents a control technique in order to deal
with a SCIG connected to the grid through a full power
converter. Both stator side and grid side are taken into
consideration. The control schemes used in each converter
have been detailed. The vector control scheme enables
the system to control the stator side converter without a
flux sensor in side the machine. This fact assures fewer
mechanical problems during the life time of the machine.
Through this work, we have established the technique of
vector control by indirect rotor flux orientation. The model
flux oriented indirect control system is efficient both speed
and current regulators are used. The indirect vector control
proposed for the stator side converter is less dependent from
machine parameters than the conventional indirect vector
control, introducing to the system a faster response. . The
control strategy has been satisfactorily evaluated by means
of Matlab/ simulations.
8. Conclusion and future research
8.1 Conclusion
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
152
Volume 17, 2022
By closing this chapter by proposing solutions for some
problems of connection of wind system to the rid side. At
the end of this work, several direct perspectives are an-
nounced and the following points are quoted by way of il-
lustration:
* Experimental realization.
* optimizes the cost of renewable energy.
*Study of fault diagnosis and isolation algorithms.
* Study of fault-tolerant control algorithms.
[1] Tseligorov N. A., Ozersky A. I., Chubukin A. V.,
Tseligorova E. N., Development of a Robust Scalar
Control System for an Induction Squirrel-cage Motor
Based on a Linearized Vector Model, WSEAS Transac-
tions on Computers, ISSN / E-ISSN: 1109-2750 / 2224-
2872, Volume 21, 2022.
[2] Budhi Muliawan Suyitno, Reza Abdu Rahman, Ismail,
Erlanda Augupta Pane, Increasing the Energy Den-
sity and Power Ratio of a Staggered VAWT Wind
Farm by using The Rotor’s Diameter as a Reference ,
WSEAS Transactions on Fluid Mechanics, ISSN / E-
ISSN: 1790-5087 / 2224-347X, Volume 17, 2022.
[3] Salmi Hassan, Badri Abdelmajid, Zegrari Mourad,
Sahel Aicha, Bagudad Abdenaceur,PSO-Backstepping
Design for sharing Active and Reactive Power in grid
connected DFIG based wind turbine,WSEAS Transac-
tions on Circuits and Systems, ISSN / E-ISSN: 1109-
2734 / 2224-266X, Volume 20, 2021.
[4] B. Wu, Power Conversion and Control of Wind Energy
Systems, August 2011, WileyIEEE Press, ISBN: 978-
0-470-59365-3.
[5] Essamudin Ali Ebrahim, Abuelmaaty M. Ali, Per-
formance and Tracking Control of Three-Phase
Induction-Motor Drive Fed from a DC-Modified Nano-
grid,WSEAS Transactions on Power Systems, ISSN /
E-ISSN: 1790-5060 / 2224-350X, Volume 16, 2021.
[6] A. LOKRITI, Y. Z. (10 au 12 mai 2010). COMPARAI-
SON DES PERFORMANCES DES REGULATEURS
PI ET IP APPLIQUES POUR LA COMMANDE
VECTORIELLE A FLUX ROTORIQUE ORIENTE
D’UNE MACHINE ASYNCHRONE. Evaluation et
optimisation des syst`
emes innovants de production de
biens et de services . Hammamet - Tunisie .
[7] Eddine, C. B. (2015). Commande vectorielle indirecte
d’un moteur asynchrone aliment´
e par un onduleur sous
d´
efaut. Alg´
erie.
[8] Radhwane Sadouni, A. M. ( 2012). Indirect Rotor
Field-oriented Control (IRFOC) of a Dual Star Induc-
tion Machine (DSIM) Using a Fuzzy Controlle. Acta
Polytechnica Hungarica , Vol. 9, No. 4,.
[9] Tarek GALLAH, M. F. (2007). Commande vectorielle
d’un moteur asynchrone par orientation de flux ro-
torique.
[10] Z. S. WANG*, S. L. (2009). Indirect Rotor Field Ori-
entation Vector Control for Induction Motor Drives in
the Absence of Current Sensors.
[11] ZERBO, M. (2008). IDENTIFICATION DES
PARAM `
ETRES ET COMMANDE VECTORIELLE
ADAPTATIVE `
A ORIENTATION DU FLUX RO-
TORIQUE DE LA MACHINE ASYNCHRONE `
A
CAGE. Qu´
ebec.
[12] Report ReGrid: Basics of wind energy. Renewables
Academy (RENAC) AG, Sch¨
onhauser Allee 10-11,
10119 Berlin (Germany)
[13] J.Chatelain, Machines ´
electriques 1 , Presse Polytech-
niques Romandes, 1983.
[14] P.Brunet, Introduction `
a la commande vectorielle des
machines asynchrones , LTEG Henri BRISSON Vier-
zon, France.
[15] N.Caliao,Dynamic modelling and control of full
yratedc onverter wind turbines, Renewable Energy
36(8)(2011)2287–2297.
[16] M.Imecs, A survey of speed and flux control structures
ofs quirrel-cage induction motor drives, Tech.Rep.,
Acta Universitatis Sapientiae Electrical and Mechani-
cal Engineering ,2009.
[17] J.Holtz, Sensor less control of induction motor drives,
Proceedings of the IEEE 90(2002)1359–1394.
[18] M. Imecs, C. Szab, J.J. Incze, Stator-field-oriented
control of the variable-excited synchronous motor: nu-
merical simulation, in: 7th International Symposium of
Hungarian Researchers on Computational Intelligence,
Budapest,Hungary,2006.
[19] V. Vongmanee, Emulator of wind turbine gen-
erator using dual inverter controlled squirrel-cage
induction motor, in: The Eighth International
Conference on Power Electronics and Drive Sys-
tems,Taipei,Taiwan,2009.
[20] R.Leidhold ,G.Garca, M.I.Valla, Field-oriented con-
trolled induction generator with loss minimiza-
tion, IEEE Transactions on Industrial Electronics49
(2002)147–156.
8.2 Future research
References
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
153
Volume 17, 2022
[21] Sperstad, I.B.; Helseth, A.; Korpas, M. Valuation of
stored energy in dynamic optimal power flow of dis-
tribution systems with energy storage. In Proceedings
of the 2016 International Conference on Probabilistic
Methods Applied to Power Systems (PMAPS), Beijing,
China, 16–20 October 2016; pp. 1–8.
[22] Yang, D.; Wen, J.; Chan, K.W.; Cai, G. Dispatching
of Wind/Battery Energy Storage Hybrid Systems Us-
ing Inner Point MethodBased Model Predictive Con-
trol. Energies 2016, 9, 629. [CrossRef]
[23] Liu, W.; Liu, Y. Hierarchical model predictive control
of wind farm with energy storage system for frequency
regulation during black-start. Int. J. Electr. Power En-
ergy Syst. 2020, 119, 105893. [CrossRef]
Creative Commons Attribution License 4.0
(Attribution 4.0 International, CC BY 4.0)
This article is published under the terms of the Creative
Commons Attribution License 4.0
https://creativecommons.org/licenses/by/4.0/deed.en_US
WSEAS TRANSACTIONS on POWER SYSTEMS
DOI: 10.37394/232016.2022.17.15
A. Alazrag, L. Sbita
E-ISSN: 2224-350X
154
Volume 17, 2022