International Journal of Computer and Electrical Engineering, Vol. 2, No.

5, October, 2010 1793-8163

Simulation of Wind-Turbine Speed Control by MATLAB

Furat Abdal Rassul Abbas and Mohammed Abdulla Abdulsada
environmental advantages of harnessing a clean and inexhaustible energy source and because of the economic incentives supplied by several governments [4].However, there are still many unsolved challenges in expanding wind power. The standard controls as well as recently developed advanced controls for variable-speed, horizontal axis wind turbines have been investigated [5]. Most of the wind turbines are equipped with self-excited induction generators (SEIG). They are simple and rugged in construction and offer impressive efficiency under varying operating conditions. Induction machines are relatively inexpensive and require minimum maintenance and care. Characteristics of these generators like the over speed capability make them suitable for wind turbine application [6].Recent advancements in Power Electronics have made it possible to regulate the SEIG in many ways, which has resulted in an increased interest in the use of SEIG for small scale power generation with wind power [7]. In the last years, Matlab-Simulink has become the most used software for modeling and simulation of dynamic systems. Wind turbine systems are an example of such dynamic systems, containing subsystems with different ranges of the time constants: wind, turbine, generator, power electronics, transformer and grid [8]. There are two principle-connections of wind energy conversion. The first one is connecting the wind-generator to grid at grid frequency. While connected to grid, grid supplies the reactive VAR required for the induction machines. Often, a DC-link is required to interface the wind-generator system with a certain control technique to the utility grid. The second is connecting the wind-generator system to isolated load in remote areas [9]. The terminal frequency of a SEIG connected to an isolated load depends mainly on the rotor speed. So, a constant rotor speed brings a constant terminal frequency. The rotor speed of the SEIG is adjusted by regulating the speed of the wind turbine. In this paper, a turbine blade pitch control has been assumed for this purpose. The PID controller which employed for turbine rotor speed control and hence the frequency regulation is proposed and simulated by Matlab-Simulink software package. II. DYNAMIC MODELING OF THE WIND TURBINE [10] The wind turbine is characterized by no dimensional curves of the power coefficient (Cp) as a function of both the tip speed ratio() and the blade pitch angle(). In order to fully utilize the available wind energy, the value of () should be maintained at its optimum value. Therefore, the power coefficient corresponding to that value will become maximum also.

AbstractIn this paper, the output frequency of a self excited induction generator (SEIG) driven by wind turbine and supplies static load are controlled. The principle connections of wind energy conversion are presented. The dynamic modeling of the wind turbine and its linearization are derived. The PID controller which employed for turbine rotor speed control and hence the frequency regulation is proposed. The block diagram of the proposed speed control system which consists of speed controller, actuator model and the turbine linearized model is simulated by Matlab-Simulink software package. Index TermsWind Turbine, Speed Control, PID Controller, Wind Generating System.

ABBREVIATIONS U: Output pitch angle (Degree). Jt: Wind turbine moment of inertia (kg.m2). B: Turbine blade angle (Degree). : Turbine tip speed ratio. Vw: Wind speed. : Air density (Kg/m2). CT: Torque coefficient. Cp: Power coefficient. R: Radius of area swept by the turbine bladed (m). A: Swept area of wind turbine (m2). Tm: Prime mover input torque (N.m.). Pm: Wind turbine output power (Watt). Km: Motor gain operator (volt/rad/sec). m: Motor time constant (sec).
t

: The angular shaft speed. : Linearization coefficients. I. INTRODUCTION

, ,

As energy demands around the world increase, the need for a renewable energy sources that will not harm the environment has been increased. Some projections indicate that the global energy demand will almost triple by 2050. Renewable energy sources currently supply somewhere between 15% and 20% of total world energy demand [1]. Wind energy conversion is a fast-growing interdisciplinary field that encompasses many different branches of engineering and science. According to the American Wind Energy Association, the installed capacity of wind grew at an average rate of 29% per year [2]. At the end of 2009, the worldwide installed capacity of wind energy was over 159 MW. The prediction capacity for 2010 is over 203 MW [3].The wind energy market has grown because of the

912

International Journal of Computer and Electrical Engineering, Vol. 2, No. 5, October, 2010 1793-8163

The tip speed ratio () can be defined as the ratio of the angular rotor speed of the wind turbine to the linear wind speed at the tip of the blades. It can be expressed as follows: (1) Where R is the wind turbine rotor radius, Vm is the wind speed and t is the mechanical angular rotor speed of the wind turbine. The output power of the wind turbine, can be calculated from the following equation
2 (2) Where () is the air density, and (A) is the swept area by the blades, and
C p = (0.44 0.0167 ) sin ( 3 ) 0.00184 ( 3 ) 15 0.3

Tm V w

=
op

t ) (J t V w op

1 1 3 A C ( , ) * V w p top V w 2

op

= t R / Vw

T m V w

= K 21 + K 22 + K 23
op

(9)
=
Tm =
op

Pm =

3 AC p V w

t ) (J t op

3 V wop 1 A C p ( , ) top 2

]
op

Tm

= K 31 + K 32 + K 33
op

(10)

Where
op = Rtop Vwop

(3) Also, the torque available from the wind turbine can be expressed as: 1 Tm = A R CT V 2 w 2 (4) Where CT is the torque coefficient which is given by CT= Cp/ Then, the aerodynamic torque, Tm can be written as
t R V 3 3 R w 0.00184 t 3 V w Tm = 0.5 A (0.44 0.0167 ) sin V 15 0 . 3 w t

K=

1 AR 2

follows

(5) The fundamental dynamics of the variable-speed wind turbine are captured with the following simple mathematical model: (6) The moment of inertia of the turbine rotor is represented by Jt; ( t ) is the angular shaft speed; Tm is the mechanical torque necessary to turn the generator and is assumed constant value commanded by the generator. Because the generator moment of inertia of a direct-drive turbine is generally several orders of magnitude less than (Jt), it has been neglected. The aerodynamic torque has been described previously by equation (5). The wind turbine used in this investigation has the specification and parameters which given in appendix (A). III. DYNAMIC MODEL LINEARIZATION [11] A traditional approach to design commonly used linear controllers such as proportional-integral-derivative (PID) requires that the non-linear turbine dynamics be linearized about a specified operating point. Linearization of the turbine equation (6) would yield:
 t = t + V w + J t
Tm t =
op

 t = Tw Tm Jt

K11, K12, K13, K21, K22, K23, K31, K32, K33 are written in appendix (B). Here, , Vw, and represent deviations from the chosen operating point, tOP, VWOP, and OP. Selection of the operating point is critical to preserving aerodynamic stability in this system. The rotational speed operating point, tOP, was selected to be the desired constant speed of the turbine, 450 rpm (47.1 rad/sec). The blade-pitch and wind speed operating points were selected as (op = 90 and Vwop = 6.4 m/s). After Laplace Transformation equation (7) becomes:
J t s t = Vw (s ) + t + U (s )

(11)

D=

Jt

Let The turbine rotor shaft speed can be represented as

t = Vw (s ) + U (s ) Jt Jt 1 sD

(12)

Equation (12) describes the linearized model of the wind turbine. Such model is represented by the block diagram shown in Fig. (1).

(7)
Fig. 1 Block diagram of the linear turbine plant model.

Where the linearization coefficients are given by:

=
t
t ) (J t

=
OP

1 C p ( , ) 3 A V wop 2 t t Tm t

OP

IV. MODEL OF ACTUATOR The permanent magnet DC motor is used as an actuator turbine blade adjustment, which may be represented by the block diagram shown in Fig.(2) where Ua(s) and Uo(s) are the Laplace transform of the pitch angle input and output

= K11 + K12 + K13

op

(8) 913

International Journal of Computer and Electrical Engineering, Vol. 2, No. 5, October, 2010 1793-8163

respectively, Km is the gain constant, and m is the time constant of the permanent magnet DC motor.
Ua(s)

Km m s +1

Uo(s)

Fig. 2 Block diagram of the actuator

V. SPEED CONTROLLER USING PID The PID controller is used for controlling the rotor speed. Such controller is illustrated in Fig. (3), where t(s) represents the input rotor speed (error signal), and Uc(s) represents the output pitch angle change. The transfer function between the input rotor speed and the output pitch angle change can be described as follows:
C( s ) = 2 U c ( s ) K p s + Ki + K d s = s t ( s )

Fig. 4 Block diagram of the wind turbine speed control system

(13)

T(s)

Kp +

Ki + Kd s s

Uc(s)

Fig. 3 Block diagram of the speed controller VI. COMPLETE SPEED CONTROL SYSTEM The proposed wind-turbine speed control, shown in Fig. (4), is simulated by Matlab-Simulink software package. The systematic approach to PID controller design provides a means of visually observing the effect of gain changes on the Root Mean Square (RMS) speed error [10].So in order to assess controller performance, the root mean square of the error between the actual rotational speed and the desired rotational speed indicates the capability of the controller to reject the wind speed fluctuations. The simulation was used repeatedly. Each of the gains was varied over a wide region. Observation of the system response inputs provides direction in choosing gain values. The PID controller parameters are tuned according this approach to give the optimal performance: KP = 15, Ki = 20 and Kd = 0.1. The wind speed is changed between 4.7 m/s to 8.1 m/s. In response, the control system will exhibits a corresponding change in blade angle between 5 deg. to 15 deg. in order to keep the rotor speed constant at the reference value 47.1 rad/s, as shown in Fig.(5).

Fig. 5 Time series traces of turbine Performance using PID-controller

VII. CONCLUSION In this paper, the output voltage frequency of SEIG driven by wind turbine and supplies static load is controlled. Frequency regulation has been presented in this paper via controlling the rotor speed. A PID controller is employed for turbine rotor speed control and hence the frequency regulation. The proposed wind-turbine speed control is simulated by Matlab-Simulink software package. The PID controller parameters are tuned to give the optimal performance: KP = 15, Ki = 20, and Kd = 0.1. The wind speed is changed between 4.7 m/s to 8.1 m/s. In response, the control system will exhibits a corresponding change in blade angle between 5 deg. to 15 deg. in order to keep the rotor speed constant at the reference value 47.1 rad/s. REFERENCES
[1] H. H. El-Tamaly and Adel A. Elbaset Mohammed, Modeling and Simulation of Photovoltaic/Wind Hybrid Electric Power System

914

International Journal of Computer and Electrical Engineering, Vol. 2, No. 5, October, 2010 1793-8163
Interconnected with Electrical Utility ", 12th International Middle-East Power System Conference, MEPCON 2008, pp 645-649, 2008. [2] L. Y. Pao and K. E. Johnson, A Tutorial on the Dynamics and Control of Wind Turbines and Wind Farms, Proc. Amer. Ctrl. Conf., June 2009. [3] World Wind Energy Association, World wind energy installed capacity. Available: http://www.wwindea.org, Accessed April / 1 / 2010. [4] Editor: K. Heinloth, " Energy Technologies Sub volume C: Renewable Energy", Springer-Verlag Berlin Heidelberg, 2006. [5] Jason H. Laks et. al., " Control of Wind Turbines: Past, Present, and Future", 2009 American Control Conference, June 10-12, 2009, St. Louis, Missouri, USA. Available: http://ecee.colorado.edu/~pao/anonftp/LaksPaoWright_ACC09.pdf [6] Sathyajith Mathew, "Wind Energy Fundamentals, Resource Analysis and Economics", Springer-Verlag Berlin Heidelberg, 2006. [7] Raja Singh Khela et. al., "Application of Artificial Neural Network for Analysis of Self-Excited Induction Generator", Journal of Computer Science and Technology (JCS&T) Vol. 6, No. 2, PP 73-79, October 2006. [8] Florin Iov et. al.," Wind Turbine Blockset in Matlab/Simulink", UNI.PRINT Aalborg University, Denmark, March 2004. [9] Ezzeldin. S. Abdin and X. Wilson, Control design and dynamic performance analysis of a wind turbine-induction generator unit IEEE Trans. On Energy Conversion, Vol. 15, No. 1, pp 91-96, March 2000. [10] Maureen Hand, Variable-Speed Wind Turbine Controller Systematic Design Methodology: A Comparison of Non-Linear and Linear Model-Based Designs, NREL Report No. TP-500-25540, National Renewable Energy Laboratory, July 1999. [11] R. M. Hilloowala, Control and interface of renewable energy systems", Ph. D. Thesis, The University of New Brunswick, Canada, 1992.
K 32 = 2 0.0167 KV wop op op 3 op 3 * 0.44 0.0167 op * 0.3 * cos 2 15 0.3 op 15 0.3 op

K 33 =

2 0.00184 K op 3 V wop op

APPENDIX A Wind Turbine Parameters: Rating : 1 kw , 450 rpm ( low speed side ) at Vw = 12 m/s . Size : Height = 4 m , Equator radius = 1 m , Swept area A= 4 m2 , = 1.25 kg/ m2. APPENDIX B Results of Linearization of the Wind Turbine Equation
KV 3 R wop 0.44 0.0167 op * * cos K11 = Vwop (15 0.3 op ) RTop

op 3 15 0.3 op

K12 =

3 KVwop 2 RTop

op 3 * 0.44 0.0167 op * sin 15 0.3 op

3 2 + 3 Vwop K 13 = 0.00184 * opVwop 2 RTop

3KV 2 wop K 21 = 0.44 0.0167 * op * RTop

op 3 * sin 15 0.3 op

KV 3 wop K 22 = 0.44 0.0167* op * * R Top

op 2 V wop 15 0.3 op

*cos

op 3 15 0.3 op

9 opVwop K 23 = 0.00184 K 2V wop op op

K 31 =

2 0.0167 KVwop op

3 op * sin 15 0.3op

915