NONLINEAR MPPT CONTROL OF SQUIRREL CAGE INDUCTION GENERATOR-WIND TURBINE 1 M. BENCHAGRA, 2 M. MAAROUFI. 3 M. OUASSAID 1, 2 Department of Electrical Engineering, Mohammadia School of Engineering- BP 767-Rabat- Morocco 3 Department of Industrial Engineering, National School of Application Sciences, UCA, Safi- Morocco E-mail: m.benchagra@gmail.com, maaroufi@ac.ma, ouassaid@yahoo.fr ABSTRACT This paper presents a nonlinear control method to track the maximum power point of wind turbine equipped with Squirrel Cage Induction Generator (SCIG). Since the mechanical power of wind turbine is function of its shaft speed at a given wind velocity, the proposed controller provides the desired voltage at the output of Voltage Source Converter (VSC) so as to control the generator speed and then extracting maximum power from wind. The controllers exploit the backstepping scheme to calculate the required converter s control to eliminate the instantaneous errors of mechanical speed and flux. The proposed control laws are derived from the Lyapunov approach which is well suited for this nonlinear system. Finally, the simulation results verify the stability and effectiveness of the control strategy proposed. Keywords: Wind Power, SCIG, VSC, Nonlinear Control, Backstepping Controller, Lyapunov Approach, MPPT 1. INTRODUCTION Renewable energies are becoming increasingly important as alternative energy sources. Wind power is one of the most-effective systems available today to generate electricity from renewable sources. This energy is transformed into mechanical energy by wind turbine and then converted to electrical energy by electrical generator and a power converter [1]. To obtain maximum benefit from the wind energy, variable speed wind turbines are being used in general. Variable speed wind turbine systems produce variable voltage and frequency when no controller element is used. In order to extract maximum power in desirable values from variable speed wind generation systems, they must be operated together with the controller element [2]. The conventional back-to-back voltage source converter is usually used to connect the generator to the grid. The VSC can increase the system robustness and MPPT tracking control, many studies have been conducted to analyze, develop and improve the extraction of maximum power in literature [3-11]. Different types of generator can be used to realize a wind power generation (WPG). The common candidates for the generator are Doubly-fed Induction Generators (DFIG), Squirrel Cage Induction Generators (SCIG), and Permanent Magnet Synchronous Generators (PMSG). (WPGs) have so many uncertainties due to erratic nature of wind-based systems. Therefore, the controller should accommodate the effects of uncertainties and keep the system stable against a large variation of system parameters. The conventional PI-based controllers cannot fully satisfy stability and performance requirements. On the other hand, the system is highly nonlinear and has a large range of operating points. Thus, linearization around one operating point cannot be employed to design the controller. Nonlinear control methods can be used to effectively solve this problem [1] and [3]. In attempt to achieve high performances in the steady state as well as during the transients, a different nonlinear control structure must be applied. In the recent two decades, many modified nonlinear state feedback such as input-output feedback linearization, control based on sliding mode have been applied to more improve the induction machine control performances [12]. However, the behavior of the SCIG-wind turbine system using backstepping approach with virtual contron has not been tested. The objective of this paper is to implement the backstepping control to wind turbine equipped with SCIG, so as to regulate mechanical speed provides by MPPT algorithm and flux to their input references even during a wind speed conditions, and that a perfect conversion of wind power to electrical power. 26
Dynamic modeling and simulation of the SCIGwind turbine power generation system and backstepping controllers are performed by means of Matlab/Simulink. Behavior of SCIG-wind turbine depending on wind velocity variation is observed through the simulation study. 2. WIND TURBINE MODEL (3) ω t is the turbine speed, it is a function of mechanical speed ω m and gear ratio G. We consider a generic equation to model a power coefficient C p, based on the modeling turbine characteristics described by [13]: 2.1 System Configuration The schematic diagram of the studied wind Wind turbine SCIG PWM abc/dq Nonlinear MPPT Backstepping Controller VSC abc/dq turbine equipped with SCIG is shown in Fig.1. The rotor shaft of the studied SCIG is coupled to the hub of turbine by through a gearbox for transforming low rotational speeds of turbine to the required higher rotational speeds of the SCIG. The mechanical energy is transforming to electrical energy by the generator. The nonlinear backstepping controller provides proper switching signals for the VSC so as to extract maximum power from the wind. 2.2 Wind Turbine The wind turbine input power usually is 1 2 (1) The output mechanical power of wind turbine is: 1 2, (2) where δ is defined as the ratio of the tip speed of the turbine blades to wind speed: Figure 1. The overall diagram of the proposed control,.519 116.4 5 exp 21.68 1.35.8 1 (4) In order to make full use of wind power, in low wind speed β should be equal to zero [14]. Fig. 2 illustrates the wind turbine power curve when β is equal to zero. From the Fig.2 we can see there is one specific mechanical speed (or generator speed) at which the output power of wind turbine is optimal. Connected all the maximum power point of each power curve, the optimal power curve (P m- opt curve) is gotten. When the wind turbine is in the P m-opt curve, the turbine will get the maximum power P max. In this case, the maximum value of power coefficient C p is C Pmax =.47, is achieved for β= and δ opt =8.1, Fig. [3]. 27
Mechanical Power, P m (W) Power coefficient, C p 3.5 2.5 1.5.5 4 x 15 3 2 1.6.5.4.3.2.1 P m opt Figure. 2. Wind turbine power curve β = β =2 v =9 m/s v =7.1 m/s β =5 β =15 v =12 m/s v =11 m/s 5 1 15 2 25 3 Generator speed, ω m (rad/s) (8.1;.47) 5 1 15 Tip speed ratio, δ Figure. 3. Characteristics Cp vs δ, for various values of pitch angle β 2.3 Squirrel Cage Induction Generator Model The SCIG mathematical model can be expressed in a reference frame rotating with arbitrary angular speed by the following equations [14]: (5) where subscripts s and r refer to the stator and rotor side respectively, subscripts d and q refer to the d-axis and q-axis respectively. R s and R r are, respectively, the stator and rotor phase resistances, ω is the electrical speed, λ is the rotor-flux The stator and rotor flux can be expressed as [14] (6) where i ds, i qs, i dr, i qr are, respectively, the direct and quadrature stator and rotor currents and L m is the mutual inductance. The electromagnetic torque can be calculated as [14] (7) where p, is the pair-pole number The voltage and flux equation can be supplemented by the mechanical equation for the drive train (gearbox) by (8) to complete the model of generator used in this paper (8) where T m is the mechanical torque. J is the total inertia, it can be calculated as: (9) where J t and J g are respectively the turbine inertia and the generator inertia The above model (5) to (9) can be presented as differential equations for the stator currents and rotor flux vector components under the following form: (1) (11) (12) (13) (14) Where: ; 28
; 1 ; ; ; ; 1 3. NOLINEAR MPPT CONTROL STRATEGY FOR SCIG-WIND TURBINE The basic idea of the backstepping design is the use of the virtual control quantities to decompose systematically a complex non-linear control structure problem into simpler and smaller ones. Backstepping design is divided into two various design steps. In each one we deal with an easier and single-input-single-output design problem and each step provide a reference for the next design step [9] From the presented above model it seen that the squirrel cage induction generator is a nonlinear multivariable system with coupling between direct and quadrature axis. This model is dependent on the stator electrical speed ω s, this variable can be calculated if dq reference frame synchronized with rotor flux. Fixing the d-axis of the rotating reference frame on the rotor flux vector, we have λ dr =λ r and λ qr =. In this case, two state variables have been proposed for describing the SCIG model in order to extracting maximum power from turbine, mechanical speed and rotor flux, as follows: (15) (16) Therfore, the error is defined using the rotational speed and rotor flux and the first Lyapunov function is defined as 1 2 1 2 (17) The time derivative of V is given by (18) i ds and i qs are considered to be the virtual controls input. Thus, there references are obtained as: 1 (19) 1 (2) k ω and k λ are positive design constants that determine the closed loop dynamics. Equations (19) and (2) indicate that the virtual controls should be in order to satisfy the control objectives. So they provide the references for the next step and try to make the signals i ds and i qs behave as desired. So, we define again the error signals involving the desired variables in (19) and (2): (21) (22) Finally, we extended the Lyapunov function in (17) to include the states variables e ids, e iqs as: 1 2 (23) Control laws are derived by differentiating the Lyapunov function with respect to time. 1 1 (24) (25) In this case, the Lyapunov function derivative is given by (26) The controller of nolinear MPPT is implemented by (24) and (25) which is introduced in Fig. 4 and simulated by using Matlab/simulink. 29
λ ref =1 + - Equation (16) λ dr-est ω m-ref ω m + - Equation (24) v ds Equation (25) v qs i ds i qs - + - + Equation (2) Equation (19) i * qs i * ds Figure 4. Detailed diagram of the nonlinear backstepping controller scheme 15 2 v (m/s) 1 5 ω mref and ω m 15 1 5 2 4 6 8 1 Time [s] (a) 2 4 6 8 1 Time [s] (c).5 1.4 C p.4.3.2.1 λ ref and λ dr 1.2 1.8.6.4.2 2 4 6 8 1 Time [s] (b) 2 4 6 8 1 Time [s] Figure 5. SCIG wind power performances with nonlinear Backstepping control (d) 3
4. SIMULATIONS RESULTS The model of the SCIG based variable speed wind turbine system of figure 1 is built using MATLAB\SIMULINK dynamic system simulation software. The parameters of the turbine and SCIG used are given in appendix, respectively in Table I and Table II. The VSC and the nonlinear backstepping control algorithm are also implemented and included in the model. The PWM VSC is operated at 1 khz. Figure (5) shows the response of the system for a change of wind velocity, from 1 m/s to 12 m/s to 9 m/s and then comes back to 1 m/s. Fig (a) shows the wind velocity profile imposed. At t=.5 s the mechanical speed reference provides by MPPT is applied. During test, the power coefficient C p is kept at their optimal value (C P =.47), simulation result is shown in Fig (b). Fig (c) illustrates the mechanical speed tracking trajectory for a change in the wind velocity with the proposed nonlinear backstepping controller, according to the wind turbine characteristics, the optimal mechanical speed of SCIG is calculate as: 23 8.1 14 So, SCIG mechanical speed reference is 14. According to ω m-ref, when wind velocity v is 1 m/s, the optimal mechanical speed is 14 rad/s, and when wind velocity v is 12 m/s, the optimal mechanical speed is 168 rad/s, etc. It is seen that in MPPT control, the generator speed has to change to keep at any wind velocity maximum power coefficient (C p =.47) and so, extraction of maximum power. Fig (d) shows the rotor flux, as desired, the direct rotor flux component is perfectly tracked to their input references with the proposed nonlinear control strategy. As can be seen, the adopted control strategy is capable of providing accurate tracking for mechanical speed and flux references and, meanwhile, keeping good decoupling of electromagnetic torque and flux regulations, finally, an excellent MPPT control is achieved. 5. CONCLUSION This paper presented a nonlinear backstepping approach for maximum power point tracking MPPT of SCIG-wind turbine generation system. MATLAB/similink simulation results show the stability and effectiveness of this strategy of control. The validity of this proposed method is confirmed by theoretical analysis and simulation results. The nonlinear backstepping control theory can be easily exploited in the wind power generation systems equipped with SCIG-based variable speed. APPENDIX: TABLE I PARAMETERS OF THE TURBINE Item Symbol Value Density of air ρ 1.22 kg/m3 Area swept by A 615.8 m2 blades Speed-up gear G 23 ratio Base wind speed v n 12 m/s Turbine inertia J t 5 kg.m2 TABLE II PARAMETERS OF THE SCIG Item Symbol Value Rated power P n 3 kw No. of poles p 2 Rated speed ω n 158.7 rad/s Stator resistance R s.63 Ω Stator inductance L s.118 H Mutual L m.116 H inductance Rotor resistance R r.48 Ω Rotor inductance L r.116 H REFRENCES: [1] M. Pahlevaninezhad, A. Safaee, S. Eren, A. Bakhshai, P. Jain, Adaptive Nonlinear Maximum Power Point Tracker for a WECS Based on Permanent Magnet Synchronous Generator Fed by a Matrix Converter, Energy conversion congress and exposition, ECCE 29 IEEE, pp. 2578-2583. [2] Y. Oguz, I. Guney, Adaptive neuro-fuzzy inference system to improve the power quality of variable-speed wind power generation system, Turk. J. Elec. Eng & Com. Sci. vol. 18. No. 4, 21, pp. 625-646. [3] X. Zheng, L. Li, D. Xu, J. Platts, Sliding Mode MPPT Control of Variable Speed Wind Power System Power and Energy Engineering Conference, APPEEC 29. Asia-Pacific, pp. 1-4. [4] A. Meharrar, M. Tioursi, M. Hatti, A. B. Stambouli, A variable speed wind generator maximum power tracking based on adaptive neuro-fezzy inference system, Journal of 31
Expert Systems with Application, Vol. 38, No. 6, june 211. [5] W. M. Lin, C. M. Hong, Intelligent approach to maximum power point tracking control strategy for variable-speed wind turbine generation system, 7th International Conference on sustainable Energy Technologies, Vol. 35, No 6, June 21, pp. 2444-2447. [6] M. Kesraoui, N. Korichi, A. Belkadi, Maximum power point tracker of wind energy conversion system,. Journal of Renewable Energy, Vol. 25, No 1, 21, pp. 181-19. [7] L. G. Gonzalez, E. Figueres, G. Garcere, O. Carranza, Maximum-power-point tracking with reduced mechanical stress applied to windenergy-conversion-systems. Journal of Applied Energy, Vol. 87, No. 7, 21, pp. 234-2312. [8] T. Senjyu, Y. Ochi, Y. Kikunaga, M. Tokudome, A. Yona, E. B. Muhando, N. Urasaki, T. Funabashi, Sensor-less maximum power point tracking control for wind generation system with squirrel cage induction generator Journal of Renewable Energy, vol. 34, No. 1, 29, pp. 994-999. [9] M. Benchagra, M. Maaroufi, M. Ouassaid, Study and Analysis on the Control of SCIG and its responses to Grid Voltage Unbalances, International Conference on Multimedia Computing and Systems, IEEE ICMCS 11, 211 [1] Y. Errami, M. Ouassaid and M. Maaroufi, Modelling and Control Strategy of PMSG Based Variable Speed Wind Energy Conversion System International Conference on Multimedia Computing and Systems, IEEE ICMCS 11, 211 [11] M. Hilal, M. Maaroufi, M. Ouassaid, Doubly Fed Induction Generator Wind Turbine Control for a maximum Power Extraction International Conference on Multimedia Computing and Systems, IEEE ICMCS 11, 211 [12] A. L. Nemmour, F. Mehazzem, A. Khezzar, M, Hacil, L. Louze, R. Abdessemed, Advanced Backstepping controller for induction generator using multi-scalar machine model for wind power purposes, Journal of Renewable Energy, Vol. 35, No 1, 21, pp. 2375-238. [13] R. Pena, J.C. Clare and G.M. Asher, Doubly fed induction generator using back-to-back PWM converters and its application to variable speed wind-energy generation, Electric Power Applications, IEE Proceeding, pp. 231-241 [14] M. Benchagra, M. Hilal, Y. Errami, M. Ouassaid, M. Maaroufi Modeling and Control of SCIG based Variable-Speed with Power Factor Control, International Review of Modelling and Simulations, Vol 4, No. 3, 211 pp. 17-114. 32