# Comparative Analysis between PI and Linear-ADRC Control of a Grid Connected Variable Speed Wind Energy Conversion System Based on a Squirrel Cage Induction Generator

Authors: Hammadi Laghridat; Ahmed Essadki; Tamou Nasser
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

This paper aims at contributing to the modeling and control of a variable speed Wind Energy Conversion System (WEC-System) based on a Squirrel Cage Induction Generator (SCI-Generator). The connection between the SCI-Generator and the main utility grid is achieved by back-to-back three phase power converters (Generator and Grid Side Converters). A new control strategy named the Active Disturbance Rejection Control (ADRC) is proposed and utilized to control the Wind Energy Conversion (WEC) system based on the SCI-Generator. The objective is to control both the generator and the grid side converters in order to operate the system and to ensure the connection with the power grid. The first converter is used to control the SCI-Generator speed and field to extract the available maximum power from the wind turbine by using a Maximum Power Point Tracking (MPPT) technique and, also, to ensure that the extracted power does not exceed its rated value in case of strong wind speeds; in this case a pitch actuator system is used to control the blades pitch angle of the wind turbine. The second converter is used to control the active and reactive powers injected into the utility grid as well as to regulate the DC-Link Voltage. This control takes into account the rejection of internal disturbances as the variation of electrical parameters (the resistance, the inductance…) and the external disturbances as voltage dips and frequency droops in the main grid. To test and validate the performances of the proposed controller, a series of simulations were developed under MATLAB/Simulink environment, and the results have demonstrated the effectiveness of the proposed control under different case of simulations.


---

## Body

## 1. Introduction

Very early in the history of technology, wind was exploited to extract mechanical energy from it, where the conversion of wind energy into mechanical energy is indeed relatively easy; it is only necessary to have a satisfactory potential and to resist the whims of excessive winds. Commonly, mechanical devices called “wind turbines” are utilized in this purpose. And for the conversion of this power into electrical energy, electrical generators are used. Generally, there are three types of generators: Doubly Fed Induction Generators (DFIG), Squirrel Cage Induction Generators (SCIG), and Permanent Magnet Synchronous Generators (PMSG). For the high-power production, the SCI-Generators are the one to be used; their main interest lies in their simple structure where the connection with the utility grid is accomplished by back-to-back power converters.The proliferation of wind turbines has led electrical engineering researchers to conduct investigations to improve the efficiency of electromechanical conversion and the quality of the energy supplied. For better operation of these systems, Maximum Power Point Tracking (MPPT) algorithms are needed to improve the energy efficiency and to extract the maximum power available from the wind turbine. In literature, several techniques were used going from the utilization of classical extreme algorithms such as Perturb and Observes (P&O), Incremental Conductance (INC), and Hil Climbing (HCS) [1] to those based on artificial intelligence as Fuzzy Logic and Neural networks [2]. However, the wind turbine has an advantages that it can be reliably characterized by its specific function Cp(λ). This is why the literature commonly refers to MPPT algorithms that seek to bring the operating point to the maximum point that this function forms; we can find among them the tip speed ratio (TSR) and Optimal Torque Control (OTC) techniques; the last one is to be used in this paper.To ensure the accurate integration of the WEC-System with the utility grid, the regulation of active and reactive powers is needed as well as the synchronization of voltage and frequency. In this context numerous publications are found in literature as, in [3], the authors proposed an indirect vector control of a Squirrel Cage Induction Generator wind turbine where the objective is to control the DC-bus voltage and the reactive power by using the classical Proportional Integrator (PI) controller, authors in [4] used the Sliding Mode Control to drive the WEC-System based on a SCI-Generator, authors in [5] proposed a Fuzzy Logic Control of a variable speed cage machine wind generation system, and, in [6], the authors have used the Backstepping technique, Furthermore, additional types of controllers were also considered for squirrel cage-based WEC-Systems as cited in [3–16]. Nevertheless, another controller appeared in the last 10 years and it is used for the control applications in several area; this controller is known as the Active Disturbance Rejection Controller (ADRC); it is proposed by Jingqing Han [10]. However, to the authors’ best knowledge, very few publications are available in the literature that address the issue of controlling Wind Energy Conversion systems, where in [11] the authors proposed the control of a wind system based on a DFIG by the ADRC technique and in [12] the authors discussed the control of a direct driven wind turbine based on Permanent Magnetic Synchronous Generator by ADRC, but for the WEC SCI-Generator the application of this control cannot be found. The main advantage of this controller lies in the real time rejection of internal and external disturbances based on an extended state observer (ESO). There are different types of disturbances, whether internal disturbances as variation of the generator parameters (resistance or inductance) and external disturbances like the instability of the grid caused either by voltage droops or frequency variation. As a result, we propose in this paper the application of Linear-ADRC controller to Wind Energy Conversion System based on SCI-Generator.The contribution of this work lies in the development of a Linear Active Disturbance Rejection controller (Linear-ADRC) to control both the generator and the grid side converters. The first control is used to regulate the mechanical speed of the generator in order to guarantee the extraction of maximum power available from the wind at different speeds; also another attention is given to the pitch angle control in order to operate the wind turbine in high wind speed and to limit the extracted power (The generator side converter). The second control is used to regulate the DC-Link Voltage as well as the active reactive powers exchange with the grid (the grid side converter). To summarize, this paper seeks to evaluate the dynamical performances and sensitivity to SCIG parameter changes with comparing the classical Indirect Field Oriented Control based on the PI controller and the new Linear Active Disturbance Rejection Control.This paper is organized as follows, where Section2 presents the modeling of the Wind Energy Conversion System based SCI-Generator. The Optimal Torque Control and pitch angle control are given in Section 3. The Wind Energy Conversion System control by PI controller is detailed in Section 4, the application of the Linear Active Disturbance Rejection Control of the Wind Energy Conversion System is demonstrated in Section 5, and, finally, the simulation results and discussions are presented in Section 6.


## 2. Wind Energy Conversion (WEC) System Model

This section deals with the Wind Energy Conversion (WEC) system model. The configuration of the SCI-Generator and it’s illustrated in Figure1. It consists of two parts. The first part consists with the turbine model, the SCI-Generator, and the AC/DC converter (The stator side). The second part of the model consists of the DC-Link, the DC/AC converter, and a filter (The grid side).Figure 1
Structure of a Wind Energy Conversion (WEC) System Using SCI-Generator.The aim of the wind turbine is to capture the wind kinetic energy and convert it to a torque that rotates the rotors blades. Thereafter, the SCI-Generator converts the mechanical power (Pm) into an electrical power to be exchanged with the grid through the power converters (AC/DC/AC). Firstly, we began with the modeling the mechanical part and then the electrical part.


### 2.1. Turbine Model (Mechanical Part)

The mechanical power available from the mind turbine rotor is given by the following:(1)Pm=12Cpλ,βρAw3where Cp is the aerodynamic coefficient; β is the blades pitch angle, λ is the turbine tip speed ratio, and A is the turbine blades area. ρ and w are representing the air density and the wind speed, respectively.The tip speed ratio is given by(2)λ=RΩturvwhere Ωtur is the turbine speed and R is the turbine radius blades.The mechanical torqueTm is given by(3)Tm=PmΩturThe coefficient Cp is a nonlinear equation which depends on the tip speed λ, and it is expressed by(4)Cpλ,β=c1c21λi-c3β-c4exp-c5/λi+c6λwhere(5)1λi=1λ+c7β-c8β3+1where c1=0.5176, c2=166, c3=0.4, c4=5, c5=21, c6=0.0068, c7=0.08, and c8=0.035.Figure2 represents Cp according to λ for different values of pitch angle β. The figure shows that Cp is maximum when β=0deg and λ=8.1.Figure 2
The Aerodynamic coefficient for WEC-System.The gearbox is a device which adapts the (slow) turbine speedΩtur to the (fast) generator speed Ωmec. This gearbox can be mathematically modeled by the following:(6)Tg=TmG(7)Ωtur=ΩmecGwhere G is the gearbox ratio.The system mechanical equation is given by(8)JdΩmecdt=Tm-Tem-fΩmecwhere Tem is the induction generator electromagnetic torque and f is the friction coefficient. J represents the total inertia of the generator (JSCIG) and the turbine (Jtur) which can be calculated by(9)J=JturG2+JSCIG


### 2.2. SCI-Generator Model

The Squirrel Cage Induction (SCI) generator dynamic model in the direct-quadrature frame is presented by(10)dΨsddt=-Rsisd+ωsΨsq+vsd(11)dΨsqdt=-Rsisq-ωsΨsd+vsq(12)dΨrddt=-Rrird+ωrΨrq(13)dΨrqdt=-Rrirq-ωrΨrdwhere isd and isq are the direct and quadrature stator currents; ird and irq are the direct and quadrature rotor currents. Ψsd and Ψsq are the direct-quadrature frame stator fields; Ψrd and Ψrq are the direct-quadrature frame rotor fields. The stator voltages in the direct-quadrature frame are vsd and vsq. Also, Rs and Rr are the stator and the rotor phases resistances; and ωr is the rotor electrical speed and ωs is the stator electrical speed.The stator and rotor fields equations are given by(14)Ψsd=Lsisd+Mird(15)Ψsq=Lsisq+Mirq(16)Ψrd=Lrird+Misd(17)Ψrq=Lrirq+Misqwhere M is the mutual inductance and Ls and Lr are the stator inductance and the rotor inductance respectively.The electromagnetic torqueTem is expressed by(18)Tem=pMLrΨrdisq-Ψrqisdwhere p is the pole pairs number.The active and reactive powers(Ps,Qs) produced by the SCI-Generator can be written as follows:(19)Ps=32vsdisd+vsqisq(20)Qs=32vsqisd-vsdisqThe state model of the SCI-Generator can be written as follows:(21)disddt=-k1isd+ωsisq+k2Ψrd+k3ΩmecΨrq+k4vsd(22)disqdt=-k1isq-ωsisd+k2Ψrq-k3ΩmecΨrd+k4vsq(23)dΨrddt=k5isd-k6Ψrd+ωs-pΩmecΨrq(24)dΨrqdt=k5isq-k6Ψrq-ωs-pΩmecΨrdwhere(25)k1=Lr2Rs+M2RrσLsLr2;k2=RrMσLsLr2;k3=pRrσLsLr;k4=1σLs;k5=MRrLr;k6=RrLr.with σ=(1-M2/LsLr) being the dispersion coefficient.


### 2.3. Grid Side Converter (GSC) Model

The grid side equations can be expressed by(26)Vga=Rfiga+Lfdigadt+viaVgb=Rfigb+Lfdigbdt+vibVgc=Rfigc+Lfdigcdt+vicApplying the d-q transformation, (26) become as follows:(27)vid=Rfigd+Lfdigddt-ωgLfiiq+vgd(28)viq=Rfigq+Lfdigqdt+ωgLfiid+vgqwhere vgd and vgq are the grid voltage components (direct and quadrature) and igd and igq are the grid current components (direct and quadrature). The voltages vid and viq are the source inverter direct and quadrature components; Rf and Lf are the resistance and the inductance of the filter. ωg=2πfg; fg is the grid frequency, it is given by the PLL (Phase Locked Loop).The above equations can be written as(29)digddt=-k7igd+ωgigq+k8vid-vgd(30)digqdt=-k7igq-ωgigd+k8viq-vgqwhere k7=Rf/Lf and k8=1/Lf.The active and reactive powers(Pg,Qg) injected to grid are given by the following:(31)Pg=32vgdigd+vgqigq(32)Qg=32vgqigd-vgdigqThe voltage equation of the DC-Link can be written as(33)Cdvdcdt=idc=is-igwhere vdc and idc are the DC-Link Voltage and the DC-Link current. The currents is and ig are the SCI-Generator side current and the grid side current, respectively, and C is the DC-Link Capacitor.


## 2.1. Turbine Model (Mechanical Part)

The mechanical power available from the mind turbine rotor is given by the following:(1)Pm=12Cpλ,βρAw3where Cp is the aerodynamic coefficient; β is the blades pitch angle, λ is the turbine tip speed ratio, and A is the turbine blades area. ρ and w are representing the air density and the wind speed, respectively.The tip speed ratio is given by(2)λ=RΩturvwhere Ωtur is the turbine speed and R is the turbine radius blades.The mechanical torqueTm is given by(3)Tm=PmΩturThe coefficient Cp is a nonlinear equation which depends on the tip speed λ, and it is expressed by(4)Cpλ,β=c1c21λi-c3β-c4exp-c5/λi+c6λwhere(5)1λi=1λ+c7β-c8β3+1where c1=0.5176, c2=166, c3=0.4, c4=5, c5=21, c6=0.0068, c7=0.08, and c8=0.035.Figure2 represents Cp according to λ for different values of pitch angle β. The figure shows that Cp is maximum when β=0deg and λ=8.1.Figure 2
The Aerodynamic coefficient for WEC-System.The gearbox is a device which adapts the (slow) turbine speedΩtur to the (fast) generator speed Ωmec. This gearbox can be mathematically modeled by the following:(6)Tg=TmG(7)Ωtur=ΩmecGwhere G is the gearbox ratio.The system mechanical equation is given by(8)JdΩmecdt=Tm-Tem-fΩmecwhere Tem is the induction generator electromagnetic torque and f is the friction coefficient. J represents the total inertia of the generator (JSCIG) and the turbine (Jtur) which can be calculated by(9)J=JturG2+JSCIG


## 2.2. SCI-Generator Model

The Squirrel Cage Induction (SCI) generator dynamic model in the direct-quadrature frame is presented by(10)dΨsddt=-Rsisd+ωsΨsq+vsd(11)dΨsqdt=-Rsisq-ωsΨsd+vsq(12)dΨrddt=-Rrird+ωrΨrq(13)dΨrqdt=-Rrirq-ωrΨrdwhere isd and isq are the direct and quadrature stator currents; ird and irq are the direct and quadrature rotor currents. Ψsd and Ψsq are the direct-quadrature frame stator fields; Ψrd and Ψrq are the direct-quadrature frame rotor fields. The stator voltages in the direct-quadrature frame are vsd and vsq. Also, Rs and Rr are the stator and the rotor phases resistances; and ωr is the rotor electrical speed and ωs is the stator electrical speed.The stator and rotor fields equations are given by(14)Ψsd=Lsisd+Mird(15)Ψsq=Lsisq+Mirq(16)Ψrd=Lrird+Misd(17)Ψrq=Lrirq+Misqwhere M is the mutual inductance and Ls and Lr are the stator inductance and the rotor inductance respectively.The electromagnetic torqueTem is expressed by(18)Tem=pMLrΨrdisq-Ψrqisdwhere p is the pole pairs number.The active and reactive powers(Ps,Qs) produced by the SCI-Generator can be written as follows:(19)Ps=32vsdisd+vsqisq(20)Qs=32vsqisd-vsdisqThe state model of the SCI-Generator can be written as follows:(21)disddt=-k1isd+ωsisq+k2Ψrd+k3ΩmecΨrq+k4vsd(22)disqdt=-k1isq-ωsisd+k2Ψrq-k3ΩmecΨrd+k4vsq(23)dΨrddt=k5isd-k6Ψrd+ωs-pΩmecΨrq(24)dΨrqdt=k5isq-k6Ψrq-ωs-pΩmecΨrdwhere(25)k1=Lr2Rs+M2RrσLsLr2;k2=RrMσLsLr2;k3=pRrσLsLr;k4=1σLs;k5=MRrLr;k6=RrLr.with σ=(1-M2/LsLr) being the dispersion coefficient.


## 2.3. Grid Side Converter (GSC) Model

The grid side equations can be expressed by(26)Vga=Rfiga+Lfdigadt+viaVgb=Rfigb+Lfdigbdt+vibVgc=Rfigc+Lfdigcdt+vicApplying the d-q transformation, (26) become as follows:(27)vid=Rfigd+Lfdigddt-ωgLfiiq+vgd(28)viq=Rfigq+Lfdigqdt+ωgLfiid+vgqwhere vgd and vgq are the grid voltage components (direct and quadrature) and igd and igq are the grid current components (direct and quadrature). The voltages vid and viq are the source inverter direct and quadrature components; Rf and Lf are the resistance and the inductance of the filter. ωg=2πfg; fg is the grid frequency, it is given by the PLL (Phase Locked Loop).The above equations can be written as(29)digddt=-k7igd+ωgigq+k8vid-vgd(30)digqdt=-k7igq-ωgigd+k8viq-vgqwhere k7=Rf/Lf and k8=1/Lf.The active and reactive powers(Pg,Qg) injected to grid are given by the following:(31)Pg=32vgdigd+vgqigq(32)Qg=32vgqigd-vgdigqThe voltage equation of the DC-Link can be written as(33)Cdvdcdt=idc=is-igwhere vdc and idc are the DC-Link Voltage and the DC-Link current. The currents is and ig are the SCI-Generator side current and the grid side current, respectively, and C is the DC-Link Capacitor.


## 3. Optimal Torque Control and Pitch Angle Control

### 3.1. Optimal Torque Control (MPPT Technique)

The Maximum Power Point Tracking is a technique used to extract the maximum power from the turbine for different wind speed. It is necessary that the turbine speed must be adjusted relatively to the wind speed to obtain the maximum available power from the turbine. Therefore, the tip speed ratio must be maintained to an optimum value of the tip speed ratio (λ=λopt) and of the power coefficient (Cp=Cpmax).The maximum mechanical power can be expressed as a function of the rotational speed:(34)Pmmax=12ρπR5Cpmaxλopt3Ωmec3G3The optimal torque (desired torque) is such that(35)Temref=CoptΩmec2where the Optimal Torque Coefficient is(36)Copt=12ρπR5Cpmaxλopt31G3The mechanical speed reference is given by(37)Ωmecref=GλoptRw


### 3.2. Pitch Angle Control

In the previous part, it has been shown that the maximum power is generated by MPPT Technique. Nevertheless, it is necessary to limit the power of the WEC-System; As a pitch angle control is used it is operated by acting on the orientation of the turbine blades areaβ, to avoid the power not exceeding its rated value in case of strong wind speeds.The equation between the powerP and the pitch angle β is approached by a first-order system, and this approach allows controlling the SCI-Generator power by a closed loop PI controller [3].The equation betweenβ and βref can be written under the following form (closed loop):(38)β=11+Tβsβrefwhere s is a Laplace operator; Tβ is the time constant of the pitch servo. The pitch angle control scheme is depicted in the Figure 3.Figure 3
Pitch Angle Control.The PI controller is used to decrease the error signal between the rated power and the measured power. It is to produce at its output the reference pitch angleβref. In normal case, the pitch angle rates are at the maximum and minimum (β˙min,β˙max) allowable values between ±10degrees/s and the angle (βmin,βmax) allowable between (0degrees, 45degrees).


## 3.1. Optimal Torque Control (MPPT Technique)

The Maximum Power Point Tracking is a technique used to extract the maximum power from the turbine for different wind speed. It is necessary that the turbine speed must be adjusted relatively to the wind speed to obtain the maximum available power from the turbine. Therefore, the tip speed ratio must be maintained to an optimum value of the tip speed ratio (λ=λopt) and of the power coefficient (Cp=Cpmax).The maximum mechanical power can be expressed as a function of the rotational speed:(34)Pmmax=12ρπR5Cpmaxλopt3Ωmec3G3The optimal torque (desired torque) is such that(35)Temref=CoptΩmec2where the Optimal Torque Coefficient is(36)Copt=12ρπR5Cpmaxλopt31G3The mechanical speed reference is given by(37)Ωmecref=GλoptRw


## 3.2. Pitch Angle Control

In the previous part, it has been shown that the maximum power is generated by MPPT Technique. Nevertheless, it is necessary to limit the power of the WEC-System; As a pitch angle control is used it is operated by acting on the orientation of the turbine blades areaβ, to avoid the power not exceeding its rated value in case of strong wind speeds.The equation between the powerP and the pitch angle β is approached by a first-order system, and this approach allows controlling the SCI-Generator power by a closed loop PI controller [3].The equation betweenβ and βref can be written under the following form (closed loop):(38)β=11+Tβsβrefwhere s is a Laplace operator; Tβ is the time constant of the pitch servo. The pitch angle control scheme is depicted in the Figure 3.Figure 3
Pitch Angle Control.The PI controller is used to decrease the error signal between the rated power and the measured power. It is to produce at its output the reference pitch angleβref. In normal case, the pitch angle rates are at the maximum and minimum (β˙min,β˙max) allowable values between ±10degrees/s and the angle (βmin,βmax) allowable between (0degrees, 45degrees).


## 4. Wind Energy Conversion System Control by PI

This section focuses on the strategy adopted to control the WEC-System using a SCI-Generator by PI Controller. Two controllers are used. The first is used to control the stator side converter, where the objective is to extract maximum power form the wind turbines (regulates the electromagnetic torque and the field of the SCI-Generator). The second controller (grid side converter control) is used to regulate the DC-bus voltage and to control the grid powers (active and reactive), which is depicted in Figure4.Figure 4
Control Wind Energy Conversion (WEC) System Using SCI-Generator.


### 4.1. Stator Side Converter (SSC) Control by PI

The PI controller is considered as the most used to control the WEC-System based on a SCI-Generator, because of its simplest implementation. The objective in this part is to regulate the electromagnetic torque and the field of the SCI-Generator to extract the maximum mechanical powerPm available from the turbine for different wind speeds.To archive the decoupled control of electromagnetic torque and field of the generator, the Indirect Rotor Field Oriented Control (IRFOC) used let(39)Ψrq=0;Ψrd=ΨrThen the SCI-Generator Model(40)disddt=-k1isd+ωsisq+k2Ψrd+k4vsd(41)disqdt=-k1isq-ωsisd-k3ΩmecΨrd+k4vsq(42)dΨrddt=k5isd-k6ΨrdThe electromagnetic torque equation is given by(43)Tem=pMLrΨrdisqThe electromagnetic torque expression is a proportion with the quadrature stator current provided that the direct rotor field is constant, for this we will regulate field of the rotor.From (42) and (43), the reference quadrature stator current and the estimated direct component of the rotor field are given by(44)isqref=LrpMΨrdrefTemref(45)Ψrdest=k5k6+sisdwhere s is a Laplace operator and the torque Temref is the optimal torque (generated by MPPT).The rotor field angle is(46)θf=∫ωsl+pΩmecdtwith the slip frequency being obtained by(47)ωsl=RrMLsisqΨrdestThe compensation and the coupling terms are gathers in the following:(48)esd=1k4ωsisq+k2k4Ψrd(49)esq=1k4ωsisd+k3k4ωmecΨrdFigure 5 depicts the control of stator side converter scheme of the SCI-Generator by PI Controller.Figure 5
The Generator Side Converter Control by PI.


### 4.2. Grid Side Converter (GSC) Control by PI

The grid side converter is used to deliver the generated power into the electrical grid; the first controller is used to regulate the DC-Link Voltage and the other controller to control the power injected into the grid during wind variations to achieve a unit power factor. Therefore for independent control of the grid active and the grid reactive power, we use Voltage Oriented (VOC) Control. Hencevgd=0 and vgq=vg, and therefore the output Pg and Qg of the WEC-System can be written as(50)Pg=32vgigq(51)Qg=32vgigdAnd the grid currents(52)digddt=-k7igd+ωgigq+k8vid(53)digqdt=-k7igq-ωgigd+k8viq-vgFrom (51) the reference current igdref is calculated by the desired delivery of reactive power:(54)igdref=23vgQgrefThe compensation and the coupling terms are expressed in the following:(55)egd=vgd-1k8ωgigq(56)egq=vgq+1k8ωgigdwhere vgd=0 and vgq=vg. Figure 6 depicts the control of grid side converter scheme of the SCI-Generator by PI Controller.Figure 6
The Grid Side Converter Control by PI.


### 4.3. The Synthesis of PI Correctors

To tune the parameters of the PI controllers (the proportional gain controllerKp and the integral controller gain Ki), we use the compensation method for the closed loop system [6]; which is depicted in Figure 7.Figure 7
The Block diagram of the regulations.C s ( s ), F(s), Cg(s), and V(s) are the transfer functions of the stator currents, rotor filed, grid currents, and DC-Link Voltage, respectively, and are given by(57)Css=ks1+τss;Fs=kf1+τfs;Cgs=kg1+τgs;Vs=1Cswhere ks=1/Rs, τs=σLs/Rs; kf=M, τf=Lr/Rr and kg=1/Rf, and τf=Lf/Rf.The PI controllers parameters (Proportional gain and Integral gain) of the stator side converter are given by(i) For the stator currents(58)Kpis=3τstrcks;Kiis=3trcks(ii) For the rotor filed(59)KpΨr=3τftrfkf;KiΨr=3trfkfwhere trf and trc are the settling time for the rotor filed and the stator currents, respectively.The PI controllers parameters of the grid side converter can be written:(i) For the grid currents(60)Kpig=3τgtrgkg;Kiig=3trgkg(ii) For the DC-Link Voltage(61)Kpvdc=2ξCω0;Kivdc=Cω02;ω0=3trvξwhere trg and trv are the settling time for the grid currents and the DC-Link Voltage, respectively, and ξ is the damping coefficient generally taken 0.707.We notice that it is necessary to design a controller for the outer loop (filed and voltagevdc) with a slower response time than the dynamics of the inner loop (stator and grid currents).


## 4.1. Stator Side Converter (SSC) Control by PI

The PI controller is considered as the most used to control the WEC-System based on a SCI-Generator, because of its simplest implementation. The objective in this part is to regulate the electromagnetic torque and the field of the SCI-Generator to extract the maximum mechanical powerPm available from the turbine for different wind speeds.To archive the decoupled control of electromagnetic torque and field of the generator, the Indirect Rotor Field Oriented Control (IRFOC) used let(39)Ψrq=0;Ψrd=ΨrThen the SCI-Generator Model(40)disddt=-k1isd+ωsisq+k2Ψrd+k4vsd(41)disqdt=-k1isq-ωsisd-k3ΩmecΨrd+k4vsq(42)dΨrddt=k5isd-k6ΨrdThe electromagnetic torque equation is given by(43)Tem=pMLrΨrdisqThe electromagnetic torque expression is a proportion with the quadrature stator current provided that the direct rotor field is constant, for this we will regulate field of the rotor.From (42) and (43), the reference quadrature stator current and the estimated direct component of the rotor field are given by(44)isqref=LrpMΨrdrefTemref(45)Ψrdest=k5k6+sisdwhere s is a Laplace operator and the torque Temref is the optimal torque (generated by MPPT).The rotor field angle is(46)θf=∫ωsl+pΩmecdtwith the slip frequency being obtained by(47)ωsl=RrMLsisqΨrdestThe compensation and the coupling terms are gathers in the following:(48)esd=1k4ωsisq+k2k4Ψrd(49)esq=1k4ωsisd+k3k4ωmecΨrdFigure 5 depicts the control of stator side converter scheme of the SCI-Generator by PI Controller.Figure 5
The Generator Side Converter Control by PI.


## 4.2. Grid Side Converter (GSC) Control by PI

The grid side converter is used to deliver the generated power into the electrical grid; the first controller is used to regulate the DC-Link Voltage and the other controller to control the power injected into the grid during wind variations to achieve a unit power factor. Therefore for independent control of the grid active and the grid reactive power, we use Voltage Oriented (VOC) Control. Hencevgd=0 and vgq=vg, and therefore the output Pg and Qg of the WEC-System can be written as(50)Pg=32vgigq(51)Qg=32vgigdAnd the grid currents(52)digddt=-k7igd+ωgigq+k8vid(53)digqdt=-k7igq-ωgigd+k8viq-vgFrom (51) the reference current igdref is calculated by the desired delivery of reactive power:(54)igdref=23vgQgrefThe compensation and the coupling terms are expressed in the following:(55)egd=vgd-1k8ωgigq(56)egq=vgq+1k8ωgigdwhere vgd=0 and vgq=vg. Figure 6 depicts the control of grid side converter scheme of the SCI-Generator by PI Controller.Figure 6
The Grid Side Converter Control by PI.


## 4.3. The Synthesis of PI Correctors

To tune the parameters of the PI controllers (the proportional gain controllerKp and the integral controller gain Ki), we use the compensation method for the closed loop system [6]; which is depicted in Figure 7.Figure 7
The Block diagram of the regulations.C s ( s ), F(s), Cg(s), and V(s) are the transfer functions of the stator currents, rotor filed, grid currents, and DC-Link Voltage, respectively, and are given by(57)Css=ks1+τss;Fs=kf1+τfs;Cgs=kg1+τgs;Vs=1Cswhere ks=1/Rs, τs=σLs/Rs; kf=M, τf=Lr/Rr and kg=1/Rf, and τf=Lf/Rf.The PI controllers parameters (Proportional gain and Integral gain) of the stator side converter are given by(i) For the stator currents(58)Kpis=3τstrcks;Kiis=3trcks(ii) For the rotor filed(59)KpΨr=3τftrfkf;KiΨr=3trfkfwhere trf and trc are the settling time for the rotor filed and the stator currents, respectively.The PI controllers parameters of the grid side converter can be written:(i) For the grid currents(60)Kpig=3τgtrgkg;Kiig=3trgkg(ii) For the DC-Link Voltage(61)Kpvdc=2ξCω0;Kivdc=Cω02;ω0=3trvξwhere trg and trv are the settling time for the grid currents and the DC-Link Voltage, respectively, and ξ is the damping coefficient generally taken 0.707.We notice that it is necessary to design a controller for the outer loop (filed and voltagevdc) with a slower response time than the dynamics of the inner loop (stator and grid currents).


## 5. Wind Energy Conversion (WEC) System Control by Linear-ADRC

This section focuses on the strategy adopted to control the WEC-System based a SCI-Generator by Linear-ADRC. The same control scheme presented in Figure4 is used to control the stator side converter and the grid side converter while using the Linear-ADRC controller instead of the PI one. The general theory of the Linear-ADRC is detailed below.


### 5.1. Linear-ADRC Controller Design

The Active Disturbance Rejection Control (ADRC) is a robust control strategy, proposed by Han in 1995 [10]. The main advantage of this control is to compensate in real time the various external and internal disturbances. This control is based on the observer (ESO) which constitutes its kernel and which makes it possible to estimate the total disturbances. There are two types, the Linear-ADRC and NonLinear-ADRC. In order to reduce the model complexity and the controller computational, a Linear-ADRC design method is proposed and applied to control the SCI-Generator based WEC-System.The observer used to estimate the process state variables and total distributions is represented by the model:(62)x˙^=Ax^+Bu+Ly-y^y^=Cx^where x^=[x^1x^2]T the estimated variables, L=[β1β2]T is the observer gain vector, and A,B,C are gain matrices.The extended state observer gains are determined by [10](63)β1=2ωobs;β2=ωobs2ωobs is determined according to the closed loop system poles to ensure both a fast dynamic of the observer and a minimum sensitivity to noise.The canonical form of Linear-ARDC can be written:(64)y˙t=fy,d,u+b0utwhere f(y,d,u) represents the total (internal and external) disturbance, b0 is known part of the process (system), d(t) represents the external disturbance, and u(t) and y(t) are input and output of the process, respectively.The combination of the linear ES-Observer and the corrector represents the Linear-ADRC whose structure is given in Figure8.Figure 8
Linear-ADRC structure.In practice, we chosekp=ωc=4/tr with tr being the system desired time. Also, we chose ωobs in a range of ωobs=3~10ωc because the ES-Observer dynamic is faster than the system dynamic.


### 5.2. Stator Side Converter (SS-Converter) Control by Linear-ADRC

To control the stator side converter by ADRC, we elaborate three-Linear-ADRC controller; two Linear-ADRC controllers are used to regulate the stator currentsisd and isq and the other one to control the direct rotor field Ψrd to track their reference. Figure 9 depicts the Active Disturbance Rejection Control of the stator side converter (SSC).Figure 9
The Generator Side Converter Control by ADRC.(i) The Stator Currents Controllers. The stator currents regulations are achieved by two Linear-ADRC controllers, where (40) and (41) are adapted to the canonical form of Linear-ARDC:(65)i˙sdq=fsdqy,d,u+b0isdqusdqFor the d-axis stator current we have(66)fisd=-k1isd+ωsisq+k2Ψrd+k4-b0isdvsdb0isd=k4;uisd=vsdAnd for q-axis stator current we have(67)fisq=-k1isq-ωsisd-k3ΩmecΨrd+k4-b0isqvsqb0isq=k4;uisq=vsqor fisd and fisq are the total (external and internal) disturbances affecting the stator currents isd and isq, respectively. uisd and uisq are the control inputs of the isd and isq, respectively. b0isd and b0isq are the known parts of the system parameters.(ii) The Rotor Field Controller. We rewrite the field expression (42) under the canonical form (CF) of a Linear-ADRC regulator (65), and the expressions of the disturbance and the known part of the parameters are(68)fΨrd=-k6Ψrd+k5-b0Ψrdisdb0Ψrd=k5;uΨrd=isdwhere fΨrd is the total (external and internal) disturbances affecting the field Ψrd. uΨrd and b0Ψrd are the control inputs and the known parts of the field Ψrd, respectively.


### 5.3. Grid Side Converter (GSC) Control by Linear-ADRC

Similarly to the stator side converter control, we elaborate three-Linear-ADRC controller; but in this case the two Linear-ADRC controllers are used to regulates the active and reactive power delivery into the grid by controls the grid currentsigd and igq ((52) and (53)), and the other controller is used to maintain constant DC-Link Voltage.(i) The Grid Currents Controllers. The grid currents regulations are achieved by two Linear-ADRC controllers, where the equations of igd and igq are adapted to the canonical form of Linear-ADRC.(69)i˙gdq=fgdqy,d,u+b0igdqugdqSubstituting (52) and (53) into (69), we can write the grid currents equation as follows.For direct axis current,(70)figd=-k7igd+ωgigq+k8-b0igdvidb0igd=k8;uigd=vidAnd for quadrature axis current,(71)figq=-k7igq+ωgigd-k8vg+k8-b0igqviqb0igq=k8;uigq=viq(ii) The DC-Bus Voltage Regulation. The voltage expression of the DC-Link can be written as(72)Cdvdcdt=is-igIf we neglect all the losses in converters, the powers Ps and Pg can be expressed as Ps=isvdc and Pg=igvdc. Hence, to obtain the voltage equation for the DC-Link(73)C12dvdc2dt=Ps-Pg(74)dvdc2dt=2CPs-3Cvgqiqgwe put W=vdc2(75)dWdt=k9Ps-k10iqgwith k9=2/C and k10=(3/C)vg.The expressions of the disturbancefvdc and the known part of the parameters b0vdc are(76)fvdc=-k9Ps+-k10-b0vdcigqb0vdc=-k10;uvdc=igqFigure 10 depicts the Active Disturbance Rejection Control of the grid side converter (GSC).Figure 10
The Grid Side Converter Control by Linear-ADRC.


## 5.1. Linear-ADRC Controller Design

The Active Disturbance Rejection Control (ADRC) is a robust control strategy, proposed by Han in 1995 [10]. The main advantage of this control is to compensate in real time the various external and internal disturbances. This control is based on the observer (ESO) which constitutes its kernel and which makes it possible to estimate the total disturbances. There are two types, the Linear-ADRC and NonLinear-ADRC. In order to reduce the model complexity and the controller computational, a Linear-ADRC design method is proposed and applied to control the SCI-Generator based WEC-System.The observer used to estimate the process state variables and total distributions is represented by the model:(62)x˙^=Ax^+Bu+Ly-y^y^=Cx^where x^=[x^1x^2]T the estimated variables, L=[β1β2]T is the observer gain vector, and A,B,C are gain matrices.The extended state observer gains are determined by [10](63)β1=2ωobs;β2=ωobs2ωobs is determined according to the closed loop system poles to ensure both a fast dynamic of the observer and a minimum sensitivity to noise.The canonical form of Linear-ARDC can be written:(64)y˙t=fy,d,u+b0utwhere f(y,d,u) represents the total (internal and external) disturbance, b0 is known part of the process (system), d(t) represents the external disturbance, and u(t) and y(t) are input and output of the process, respectively.The combination of the linear ES-Observer and the corrector represents the Linear-ADRC whose structure is given in Figure8.Figure 8
Linear-ADRC structure.In practice, we chosekp=ωc=4/tr with tr being the system desired time. Also, we chose ωobs in a range of ωobs=3~10ωc because the ES-Observer dynamic is faster than the system dynamic.


## 5.2. Stator Side Converter (SS-Converter) Control by Linear-ADRC

To control the stator side converter by ADRC, we elaborate three-Linear-ADRC controller; two Linear-ADRC controllers are used to regulate the stator currentsisd and isq and the other one to control the direct rotor field Ψrd to track their reference. Figure 9 depicts the Active Disturbance Rejection Control of the stator side converter (SSC).Figure 9
The Generator Side Converter Control by ADRC.(i) The Stator Currents Controllers. The stator currents regulations are achieved by two Linear-ADRC controllers, where (40) and (41) are adapted to the canonical form of Linear-ARDC:(65)i˙sdq=fsdqy,d,u+b0isdqusdqFor the d-axis stator current we have(66)fisd=-k1isd+ωsisq+k2Ψrd+k4-b0isdvsdb0isd=k4;uisd=vsdAnd for q-axis stator current we have(67)fisq=-k1isq-ωsisd-k3ΩmecΨrd+k4-b0isqvsqb0isq=k4;uisq=vsqor fisd and fisq are the total (external and internal) disturbances affecting the stator currents isd and isq, respectively. uisd and uisq are the control inputs of the isd and isq, respectively. b0isd and b0isq are the known parts of the system parameters.(ii) The Rotor Field Controller. We rewrite the field expression (42) under the canonical form (CF) of a Linear-ADRC regulator (65), and the expressions of the disturbance and the known part of the parameters are(68)fΨrd=-k6Ψrd+k5-b0Ψrdisdb0Ψrd=k5;uΨrd=isdwhere fΨrd is the total (external and internal) disturbances affecting the field Ψrd. uΨrd and b0Ψrd are the control inputs and the known parts of the field Ψrd, respectively.


## 5.3. Grid Side Converter (GSC) Control by Linear-ADRC

Similarly to the stator side converter control, we elaborate three-Linear-ADRC controller; but in this case the two Linear-ADRC controllers are used to regulates the active and reactive power delivery into the grid by controls the grid currentsigd and igq ((52) and (53)), and the other controller is used to maintain constant DC-Link Voltage.(i) The Grid Currents Controllers. The grid currents regulations are achieved by two Linear-ADRC controllers, where the equations of igd and igq are adapted to the canonical form of Linear-ADRC.(69)i˙gdq=fgdqy,d,u+b0igdqugdqSubstituting (52) and (53) into (69), we can write the grid currents equation as follows.For direct axis current,(70)figd=-k7igd+ωgigq+k8-b0igdvidb0igd=k8;uigd=vidAnd for quadrature axis current,(71)figq=-k7igq+ωgigd-k8vg+k8-b0igqviqb0igq=k8;uigq=viq(ii) The DC-Bus Voltage Regulation. The voltage expression of the DC-Link can be written as(72)Cdvdcdt=is-igIf we neglect all the losses in converters, the powers Ps and Pg can be expressed as Ps=isvdc and Pg=igvdc. Hence, to obtain the voltage equation for the DC-Link(73)C12dvdc2dt=Ps-Pg(74)dvdc2dt=2CPs-3Cvgqiqgwe put W=vdc2(75)dWdt=k9Ps-k10iqgwith k9=2/C and k10=(3/C)vg.The expressions of the disturbancefvdc and the known part of the parameters b0vdc are(76)fvdc=-k9Ps+-k10-b0vdcigqb0vdc=-k10;uvdc=igqFigure 10 depicts the Active Disturbance Rejection Control of the grid side converter (GSC).Figure 10
The Grid Side Converter Control by Linear-ADRC.


## 6. Simulation Results and Discussion

The simulations results of the presented control strategy are given and discussed in this section. The proposed WEC-System control has been simulated using the Matlab/Simulink. We have applied a variable wind speed profile as it is illustrated at Figure11; the parameters of the WEC-System are given in the Appendix.Figure 11
The Wind Speed Profile.As it can be noticed in Figures12 and 13, when the wind speed exceeds the rated speed, the mechanical speed of the turbine is limited by adjusting the pitch angle, and consequently the power coefficient and tip speed ratio of the turbine decrease.Figure 12
(a) The SCI-Generator Speed and (b) the blades pitch angel.
(a)
Ω m e c (b)
βFigure 13
(a) The aerodynamic coefficient and (b) the turbine tip speed ratio.
(a)
C p (b)
λFigures14 and 15 present the simulation results of the stator side converter, for the two control strategies (ADRC control and PI control).Figure 14
(a) The direct rotor fluxes of the SCI-Generator and (b) the electromagnetic torque of the SCI-Generator.
(a)
Ψ r d (b)
T e mFigure 15
The stator currents: (a) d-axisIsd and (b) q-axis Isq.
(a)
I s d (b)
I s qAs it can be seen in Figures14 and 15, the control by PI and the control by ADRC tracks their references, but the control by ADRC is faster (case of the electromagnetic torque and quadrature currents), more precise (case of the Rotor Field), and also without fluctuations (case of the direct current). In Figures 16–19 the simulation results of the grid side converter are presented.Figure 16
The DC-Link VoltageVdc.Figure 17
The grid currents: (a) d-axisIgd and (b) q-axis Igq.
(a)
I g d (b)
I g qFigure 18
(a) The output reactive powerQg of the grid and (b) the output active power Pg of the grid.
(a)
Q g (b)
P gFigure 19
The grid voltageVga and the grid current Iga.As illustrated, the Linear-ADRC approach offers a faster trackness characteristic while controlling the DC-Link Voltage compared to the conventional approach by PI, where the DC-Link Voltage is maintained at its reference without presenting an overshoot at the transitional regime. Also, it can be noticed that the Linear-ADRC regulates the grid currents to their references.The control of the active and reactive powers is also done. As can be seen in Figure18, the power injected into the grid tracks its reference (the extracted power from the wind turbine minus joules and DC-Link losses). And the reactive power was set to zero to ensure a unit power factor. Hence the voltage and current of the grid are in phase as seen in Figure 19.Finally, to test the robustness of the proposed control strategies, we have performed a change in the rotor resistance and the rotor inductance of the SCI-Generator by an increase of 50% and 10%, respectively, of their nominal value. The results obtained by LADRC controller and PI controller are highlighted in Figures20 and 21.Figure 20
(a) The direct rotor field for 50%Rr and (b) the direct rotor field for 10%Lr.
(a)
Φ r d (b)
Φ r dFigure 21
(a) The direct stator current for 50%Rr and (b) the direct stator current for 10%Lr.
(a)
I s d (b)
I s dFigures20–22 show that the Linear-ADRC control responses presents better results than the control based on PI one. The Linear-ADRC is more robust and stable during the variation of the generator parameters.Figure 22
(a) The quadrature stator current for 50%Rr and (b) the quadrature stator current for 10%Lr.
(a)
I s q (b)
I s qThe simulation results clearly indicate that the proposed Linear-ADRC strategy when applied to the WEC-System is able to track the desired values for both controls of the stator side converter and the grid side converter. Moreover, it is efficient in terms of stability and robustness regarding the variation of the generator parameters.


## 7. Conclusion

This paper deals with the modeling and control of a Wind Energy Conversion System using a Squirrel Cage Induction Generator. It seeks to evaluate the dynamical performances and sensitivity to SCI-Generator parameter changes with comparing the classical Indirect Field Oriented Control based on the PI controller and the new Linear Active Disturbance Rejection Control. The Linear Active Disturbance Rejection controller and the PI controller are proposed to operate the wind turbine so as to extract the maximum power from the wind energy and to deliver it to utility power grid. It consists, on one hand, in controlling the stator currents in order to adapt the SCI-Generator rotational speed with the wind speed by acting on the generator electromagnetic torque (MPPT technique) and also control the rotor field and, on the other hand, in controlling through the grid currents, the active and reactive powers injected into the grid.The simulation results show that the Linear-ADRC control strategy is efficient in terms of fast tracking and robustness to parameter sensitivity compared to the classical PI controller.



---

*Source: 8527183-2019-03-26.xml*