Open access peer-reviewed chapter

Modeling and Control of Wind Turbine to Damp the Power Oscillation

By Emad M. Elhaji

Submitted: November 20th 2017Reviewed: February 4th 2018Published: October 10th 2018

DOI: 10.5772/intechopen.74835

Downloaded: 62


Damping inter-area oscillation by using a permanent magnet synchronous generator (PMSG) wind turbine is considered. The PMSG wind turbine is connected to the IEEE-30 bus power system at different buses. H-infinity design controller is proposed to modulate the power where the input of the H-infinity control is the variation of the local grid generator speed and the output is feedback to activate the PMSG speed control, blade pitch angle control and dc voltage control. MATLAB/SIMULINK is used in this study. The IEEE-30 bus system is reduced to 7 buses based on the number of generators to simplify the stability study. The method is applied to a seven-area power system that exhibits undamped oscillations. Results presented in this study demonstrate the effectiveness of the wind generator in increasing system damping considerably.


  • permanent magnet synchronous generator
  • IEEE 30 bus system
  • reduced system
  • wind generator
  • H-infinity
  • inter-area oscillation
  • dynamic stability
  • small signal stability

1. Introduction

The reason to use the renewable energy is to reduce the pollution emissions such as carbon dioxide CO2 in the air. Renewable energy such as wind energy is one of the cost-effective forms of energy available today to generate electricity [1, 2]. The development of modern wind energy started almost 30 years ago but were not installed significantly for long period. Over the past decade, the power from wind turbines increased from 23,900 to 486,790 MW from 2001 to 2016, respectively [3, 4]. The blade size of modern wind turbines has increased to increase the amount of power from the wind [4, 5]. In general, two types of wind turbine can be installed to the power grid such as fixed-speed wind turbine and variable-speed wind turbine. The fixed speed wind turbines were installed in early 1990 to generate electricity regards of wind speed. This kind of wind turbines connected directly to the grid and its speed fixed at speed determined by the frequency of the grid [6, 7, 8]. After that in 1996, the variable speed wind turbine started connecting to the grid [9]. The classification of these turbines, fixed speed and variable speed wind turbines, is based on the principle of the generator and its operation [10]. The induction generator specifically squirrel cage induction generator is used with fixed speed wind turbine [11, 12]. In variable speed wind turbine two kind of generators can be used such as permanent magnet synchronous generator (PMSG) and doubly-fed induction generator (DFIG) [13, 14, 15]. Various of controls scheme can be applied in fixed and variable speed wind turbines such as stall control, blade pitch angle control, active stall control, dc voltage control, torque control and point of common coupling control (PCC) [16]. Nowadays, most of the wind turbines installed in power system are permanent magnet synchronous generator (PMSG) and doubly-fed induction generator (DFIG). PMSG connected directly to the power grid via back to back converter; therefore, the power from wind directly transferred to the grid. On the other hand, DFIG has two windings; stator winding connected directly to the power grid and rotor winding connected to the power grid via back to back converter [17, 18]. Therefore, the power from the wind transferees to the grid via stator and rotor. From the stator, the wind power directly transferees to the grid while from the rotor the power transferees to the grid or grid to the rotor depends on the operating conditions: over-synchronous condition if the power flows from rotor to grid and sub-synchronous condition if the power flows from grid to the rotor. DFIG required a slip rings to operate properly, but it has a high maintenance cost in long run for normal operation [13, 18].

The PMSG will be focused in this study. It is called also a brushless DC machine and just the rotor of the DC machine is replaced by a permanent magnet. Therefore, the PMSG will be a small in size compared to the DC machine, and the inertia of a PMSG will be higher. The result, PMSG has a good reliability and more higher power compared to the size [12, 14, 19]. In wind turbines application, the small generator in size with a higher power is preferred because the generator of wind turbine is placed hundreds of above the ground and then reduce the maintenance cost. The PMSG also has a back-to-back converter to control the flowing power from the wind, but it increases the cost of the PMSG in addition to the expensive permanent magnets [12, 14]. The maintenance of the PMSG is nearly free; therefore, it is highly desirable for wind turbine applications. Now, PMSG wind turbine is focused in area of research [12, 14].

When the wind turbines penetration to the power grid, the system stability become more significant. The PMSG is investigated for its capability to enhance the damping of dynamic oscillations in a multi-area system. The basic control applied to the PMSG is a field orientation control (FOC) [20, 21]. FOC is widely used in a research area especially in wind turbines. The FOC implemented in generator and grid sides of PMSG wind turbine to control the wind power that transferred to the grid [20, 21]. The synchronous frame of reference, a d-q frame of reference, is applied to the modeling of the generator side and grid side and a decoupled strategy is used to simplify the control design of the FOC. The proportional integral controls are used in the FOC in both sides of wind turbine as a result of a good steady state error. In generator side, the FOC is designed to control the torque of PMSG wind turbine by controlling the q-axis current of PM generator while in grid side, it is designed to control the dc voltage to guarantee the power transfers from the wind to the grid. The FOC cannot damp any oscillation if there is any disturbance happen on the system such as wind speed change, power load change or symmetrical and unsymmetrical fault. These disturbances will cause a low-frequency oscillation due to many synchronous generators and wind turbines located in many areas with long tie-line as in Figure 1 [22]. This oscillation with a range of frequency between 0.1 and 0.8 HZ is known as an inter-area oscillation [23, 24]. The system will be unstable if this oscillation increase in amplitude. Therefore, to make the system more stable, an auxiliary feedback control has to be added to damp the oscillation.

Figure 1.

Two area power system connected with many of synchronous generators and wind turbines.

In [25, 26, 27], numerous of controls are applied to the PMSG wind turbine concern with the stability under the presence of grid dynamics. In [28, 29], implemented some devices to improve the system stability such as power system stabilizer (PSS) and flexible alternative current transmission system (FACTS), but these devices may not be sufficient to damp inter-are oscillations at high penetrations levels. In some literature use the frequency feedback in DFIG to damp the oscillation [30, 31]. In this study, the dynamic of the PMSG wind turbine is discussed through a suitable small-signal model. The PMSG connected to the IEEE 30 bus system at different buses such as generator and load buses. The PMSG wind turbine model is linearized and reduced in order to simplify the analysis of the system. This model is applied to case of damping the inter-area oscillation. To damp the inter-area oscillation, three modulation feedback loops were applied: dc-ink voltage, PMSG speed and blade pitch angle modulation. The feedback design method used in the study is based on the H-infinity control concept [32, 33, 34, 35].

Section 2 describes a suitable dynamic model of the PM wind generator that can be used in dynamic studies. Section 3 demonstrates the small signal stability. Section 4 describes the H-infinity control method. Section 5 describes the simulation results and Section 5 provides the conclusion.

2. System modeling

The rotational of the wind turbine rotor blade converts the kinetic energy of the wind to the mechanical power. The mechanical power converts by the generator to the electrical power that is fed to the grid through a back-to-back converter. The wind turbine scheme as shown in Figure 2 combines the mechanical subsystem and electrical subsystem. The mechanical subsystem includes the two mass system and rotor shaft, and the electrical subsystem includes a PMSG, back-to-back converter, and grid. The interconnection system modeling of the generator and grid sides are described in an appropriate frame for developing a suitable control scheme which will be discussed in detail in the following section [35, 36, 37, 38, 39, 40, 41]. The modeling of the both sides of wind turbine as well as back to back converter will be in details in the following section.

Figure 2.

Configuration of the permanent magnet synchronous generator wind turbine.

2.1. Wind turbine model

The extracted power from the wind turbine can be expressed as in (Eq. (1)) [13, 14, 19]; where σis the air density, Ris the blade radius, vwis the wind speed, Cpλβis the turbine power coefficient and it is a nonlinear function as in (Eq. (2)). The value coefficients in (Eq. (2)) c1to c6c6are: c1=0.5167, c2=116, c3=0.4, c4=5, c5=21, and c6=0.0068[17]. Also, the torque of wind turbine can be presented in terms of the wind power and rotational speed of the shaft as in (Eq. (3)) [27, 37, 42, 43, 44].


Figure 3 shows the turbine power coefficient at different tip ratio and blade pitch angle. As seen from the same figure, the maximum value of power coefficient is 0.48 when the tip ratio is 8.1 and blade pitch angle is zero. The blade pitch angle control is zero that mean the wind speed is equal or less than the rated speed. The rated wind speed which extracted the maximum power is 12 m/s. The blade pitch angle control activates when the wind speed is more the rated value to decrease the turbine power coefficient so as to Cpλβachieve the maximum power extracted from the wind.

Figure 3.

Power coefficient and tip ratio.

2.2. Permanent magnet synchronous generator (PMSG) model

The dynamic model of the PMSG is derived based on the two-phase d-q synchronous frame of reference and Figure 4 shows its equivalent circuit [27, 43, 44]. The electrical model of the PMSG in the synchronous frame of reference is by


Figure 4.

Equivalent circuit of PMSG.

where Rsis the armature resistance, LlsLsd, Lsqare the inductance in d-q axes and they are equal in the round rotor, ψfis the permanent magnetic flux, ωe=pωgis the electrical rotating speed, ωgis related to the mechanical rotating speed of the generator, and pis the number of pole pairs.

The expressions of the electromagnetic torque Teand active power of the PMSG Peare given by [27, 44].


2.3. Rotor shaft model

The rotor shaft of PMSG wind turbine is presented as a two-mass system and the motion equations are given in (Eq. (6)) where Kis the stiffness of the shaft, θis the shaft angle, Jgis the inertia of PMSG and JTis the inertia of the wind turbine [44].


2.4. Grid model

The modeling of the grid side as shown in Figure 5a has two parts. The first part of modeling includes elements between the grid side convertor and the point of common coupling bus, which is known as a PCC bus. The first part dynamic model is based on the d-q frame of reference as shown in (Eq. (7)) [43, 44]; Where ωsis the angular frequency of the grid, Lgand Rgare the inductance and resistance of the connecting transformer and power lines between the inverter and the PCC bus, respectively, vdpcc, vqpcc, udivdi, and vqiuqirepresent the d- and q-axis components of the PCC voltage and grid side convertor, respectively, and igd,igqrepresent the d-q-axis components of the current flowing between the grid side convertor and the AC system. The expression of the active power is represented in (Eq. (8)) [43, 44].


Figure 5.

Grid side; (a) two parts of grid sides, (b) 30-bus system.

The second part of the modeling is the grid model as shown in Figure 5b which is presented the IEEE 30 bus system using here as a sample of a 7-area system, where Area 5 and 7 are equipped with the PMSG wind turbine [45]. To simplify study the system stability, the power system loads have to be converted to the impedance as in (Eq. (9)), Therefore, the generators and loads impedances of the power system have to be include in the Y-bus. Also, the number of the buses in the original IEEE 30-bus have to be reduced to the number of generators [46]. The IEEE-30 bus system has 6 generators and two PMSG wind turbines. One of the PMSG wind turbines connected to the load bus and the another wind turbine connected to the PV bus. Therefore, The IEEE 30-bus reduced to the 7 bus or 7 areas system. In the area 5, the load was added to provide the disturbance of the system.


The rotor variational dynamics of the synchronous generators, generator 1,2 to 7, in the reduced system are described in (Eq. (10)) [22], where Djis the machine damping coefficient, Ksij=ViVjYijsinθijδi+δjare the synchronizing coefficient, Pwis the wind power, Pmiis the power that is from grid generators, ωiis the speed of mass of the generators, and δiis the angular position of the mass. Area 7 connected to the load bus where ω7=KVdcand Kis a constant value. The unit of t is in s, ωin per-unit and ωo=314rs,.


2.5. DC-link model

The dc-link capacitor as shown in the Figure 2 is connected between the generator and grid converter sides. The energy stored in the capacitor Cis given by Edc=12CVdc2, where Vdcis a dc-link voltage, and it depends on the difference between the generator power Peand grid power Pgif the losses of both converters are considered very small as [8, 44].


To sending power from the wind turbine to the grid, the variation of the dc voltage has to be small or equal to zero which means the sending power Peand receiving power Pgare equal.

3. Small signal stability

In this study, the small signal stability of the PMSG wind turbine connected to the power grid is going to be discuss at a particular operating point. The small signal stability is ability of the system to maintain stable during any small disturbance occurrence. In the power system, disturbances with small magnitude are very common and can occur due to change in the generation power, power load. To study the small signal stability, the nonlinear system modeling such as wind turbine and dc-link modes has to be linearized to derive the state space matrix for the entire system by using the component connection techniques [22].

3.1. Linearization of the wind torque model

As in (Eq. (3)), the torque of the wind turbine is a function of the Cpλβ; therefore, it has to be linearized respect to the wind speed, rotational speed, and blade pitch angle at operating point as [44].




3.2. Linearization of a DC-link model

After substituting the sending and receiving power in to (Eq. (11)), the linearization of the d-link can be presented as


where Δvdcis a variation of a DC voltage.

After describing the system modeling and its linearization, the interconnection method between the system modeling with its control has to be apply to create the state space equation for the entire system to study the system stability. Figure 6(a–e) shows the generator and grid sides with their controls. The way how to connect all the system models and controls together will be presented by getting the state variable of the entire system modeling with their controls as follows [47].


Figure 6.

Connection of the PMSG models with its controls.

The system modeling of the PMSG with its controls as in Figure 6 has to be connected together to get the state space equation. The interconnection matrix which connect the PMSG wind turbine models together with their controls is as [47].


Matrix Eis the states input and control inputs, u, as described in Figure 6 while matrix Frepresents the reference input of the controllers. By substituting the output of (Eq. (14)) to (Eq. (15)), the input of the system will be as


By substituting the input Uin (Eq. (16)) into the (Eq. (14)), the (Eq. (14)) can be represented as [47].


where Xis a state variables of the whole system and its controls, Uis an input control, Yis an output of the system, Â=A+B1ED1ECis a state matrix for the whole system and its control, B̂=B1ED1Fis the input system matrix after connecting with its control, Ĉ=C+D1ED1ECis the output matrix of whole system, and D̂=D1ED1Fis the direct transition matrix of the entire system. By getting the state matrixes of the entire system the eigenvalue sensitivity method can be study. The original system of PMSG wind turbine with its control in Figure 2 is of 15th order which is presented as ΔisdΔisqΔωTΔωgΔθΔvdcΔigdΔigqΔx1Δx2Δx3Δx4Δx5Δx6Δθpllcorresponding to the d-current of PMSG, q-current of PMSG, mechanical speed of the wind turbine PM generator speed, rotor shaft angle, dc voltage, d-current of grid, q-current of grid, PMSG speed control, PMSG speed control, dc-voltage control, q-current of grid control, PCC voltage control, PLL control and PLL angle, respectively; where; this original system can be reduced by ignoring the fast dynamic mode into 5th order as in Figure 7. As a result, the slow dynamic modes of the system are the Δisq, ΔωT,Δωg,Δθand Δx. With reference to Figure 2, the PI speed controller determines the q-axis stator voltage. The state equation of this controller is given by (Eq. (18)) and its output by (Eq. (19)). According to the same figure, the vsddetermined such that isdis zero for unity power factor operation. Because of the response of this loop is considerably faster than the inter-area dynamics of interest, the d-axis of (Eq.(4)) can be eliminated from the motor state equation by assuming that Δisd=Δ Therefore, using (Eq. (19)) into q-axis of (Eq. (4)) and the previous condition, one obtains (Eq. (20)) that describes the simplified q-axis dynamics [44].

Figure 7.

Reduced system.


The reduced system equation is described in (Eq. (6)), (Eq. (18)), and (Eq. (20)), which can be written as [44].


4. Hcontrol method with wind generation

Figure 8 shows the H-infinity feedback control from grid to the PMSG speed, dc-voltage, and blade pitch angle controls. The wind generators in Area 5 and Area 7 act as a power source in the dynamic equation to increase the system damping by using the local frequency of machine 5 and wind turbine in Area 7, in order to achieve power modulation proportional to the frequency deviation of the area, e.g. ΔPw=DΔω5and ΔPw=DΔω7. From Figure 8, the input power to the grid ΔPwcan be presented as


Figure 8.

Power modulation control of the PMSG wind generator based on the H∞ controller.

The feedback signal is the locally derived system frequency deviation, Δω5and Δω7, and the control outputs are the variations of the blade pitch angle, Δβ, the rotor speed reference, Δωg, and the rate of change of the dc link voltage reference, ΔV̇dc. To design the H-infinity control, the plant consisting of the wind generation and the power system is considered controllable and observable. The state representation of the plant including the disturbance is given by (Eq. (23)) [48].


where is the state vector, z=Δω5Δω7Tis the fictitious output, y=Δω5Δω7is the measurement vector, u=ΔβΔω2ΔV̇dcis the control input vector, and d=ΔvwΔPLis the disturbance vector of the plant. A, Bd,Bu, C1, D11, D12, C, D21and D22are constant matrices of appropriate size depending on the power system order. For simplicity it is assumed that there is no direct feedback through term, i.e. D22=0. The state representation of the Hcontroller is given by (Eq. (24)) [48].


where K=ÂB̂ĈD̂is the controller matrix and x̂is the controller state. The state representation of the closed-loop system is obtained by combining (23, 24); therefore:


The closed-loop transfer function, Fs, from dto zwill be denoted as


Matrices Â, B̂, Ĉ, D̂are selected that the fsis stable and the infinite norm, fs<γwhere, γis a specified positive number.

5. Simulation studies

The IEEE 30-bus system as in Figure 5(b) was studied. The data for the lines, generators and loads are shown in [45]. Generator at bus 1 is the slack bus, in the simulations, and it is equipped with a governor for regulating frequency. The data of the wind generator on area 5 and area 7 are given in Table 1. Initially 50 wind generators were installed in each area, area 5 and area 7.

Pole pairs40
Rated power per generator2 MVA
Generator rated voltage575 V
Frequency50 Hz
Stator resistance0.0025 pu
Inductance in d-axis0.43 pu
Inductance in q-axis0.43 pu
Permanent magnet flux1pu
PMSG inertia2.5 × 105 kgm2
Wind turbine inertia3.7 × 106 kg m2
DC link voltage1500 V
DC-link capacitor per generator60,000 μF

Table 1.

Wind generator parameters [27].

At the steady state generators 1–7 generate respectively 38, 43, 59, 5, 38, 13, and 77 MW. The initial wind velocity at the wind generator site is 11 m/s, and the wind generator generates an additional 77 MW of power. The power system dynamics are of the 14th order according to (Eq. (10)). The combined system including the dynamics of the wind generator from (Eq. (21)) is of the 19th order. The system and controller above was simulated in the time domain in MATLAB/SIMULINK. Two case studies were considered and the results are presented below.

Case 1: Load change.

The load on Bus 5 changes from 50 to 60 MW at t = 1 s and the was operated at steady state. Figure 9 shows the responses of the system with and without the feedback control from the wind generator. The result shows a power unbalance in the area 5 initiating and the net power exported by each area causes a large oscillations. As the system is adjusting to the new steady state, the exported power from the slack Bus 1 is increasing and that of Area 5 is decreasing by the same amount equal to the load change of 10 MW. From the figure the changing of the load power is absorbed by the infinite bus. Another Areas start oscillation when the load change and goes to the same operating point.

Figure 9.

Case 1, load change: (a) net power of area 1; (b) net power of area 2; (c) net power of area 3, (d) net power of area 4; (e) net power of area 5; (f) net power of area 6.

Figure 10(a) shows the response of the wind power resulting from the disturbance. The wind generator response is immediate providing approximately a peak variation of 5 MW power. The modulation power is created by the deviation of the dc link voltage, blade pitch angle, and generator speed shown in Figure 10 (b–d) respectively. The deviation of the dc link voltage of only 100 V or 6.6% is sufficient, and variation of the blade pitch angle is less than 0.1 degree and PMSG speed is very small.

Figure 10.

Load change: (a) wind power; (b) DC voltage response; (c) pitch angle response; (d) PMSG speed response.

Case 2: Wind speed change.

The system operates at the same steady state described above. At t = 1 s, the wind generator power increases from 77 to 100 MW when the wind velocity changes from 11 to 12 m/s. Therefore, the additional power injected into areas 5 and 7 initiate acceleration of the machines at areas 5 and 7 and, subsequently, similar response at the other machines. Figure 11(a–f) shows the net power with and without power modulation feedback when the wind speed changed from 11 to 12 m/s at t = 1 s. The event causes a power unbalance in the area served by area 5 and 7 initiating exhibiting large oscillations. As the system is adjusting to the new steady state, the exported power from the slack Bus 1 is decreasing and that of area 5 and area 7 are decreasing by the same amount equal to the wind power change by of 46 MW. The effectiveness of the wind power modulation by using H-infinity feedback control can be seen in Figure 11(a–f) comparing to the responses without power modulation. The oscillation damping has increased significantly and the system arrives at the steady state approximately 20s.

Figure 11.

Case 1, load change: (a) net power of area 1; (b) net power of area 2; (c) net power of area 3. (d) Net power of area 4; (e) net power of area 5; (f) net power of area 6.

Figure 12 shows the responses of the system with and without the feedback control from the wind generator at t = 1 with changing of the wind speed from 11 to 12 m/s. The output of the wind generator power as in Figure 12(a) is smooth increasing to accelerate the PM-generator. The oscillations in the dc link voltage are greater in the case where modulation feedback is used resulting partly in the stabilizing power modulation. The deviation of the dc link voltage of only 20 V or 1.3% is sufficient. The blade pitch angle, in Figure 12(c), is activated only by the power modulation feedback and the variation of the blade pitch angle is less than 0.2°. The generator speed, in Figure 12(d), immediately accelerates as a result of the increasing wind power input and its variation is very small.

Figure 12.

Wind speed change: (a) wind power; (b) DC voltage response; (c) pitch angle response; (d) PMSG speed response.

6. Conclusion

The power modulation by using a PMSG wind turbine to damp the power oscillation in a multi-area power system is studied. H-infinity control is applied as a feedback control; where, the input is the variation of the local grid generator speed and the output feedback to the dc-voltage, blade pitch angle and PMSG speed controls. The high-order dynamic model of the wind generator system is reduced to simplify the analysis. The PMSG wind turbine connected to the IEEE 30 bus in different areas such as area 5 and 7. The results showed the effectiveness of the wind generator in significantly increasing damping. The dc link voltage modulation provides more damping to the grid power oscillations as it compares to the PMSG speed and rotor blade pitch angle modulations. The MATLAB/Simulink is used to design the controller and validate the results.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Emad M. Elhaji (October 10th 2018). Modeling and Control of Wind Turbine to Damp the Power Oscillation, Stability Control and Reliable Performance of Wind Turbines, Kenneth Eloghene Okedu, IntechOpen, DOI: 10.5772/intechopen.74835. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

62total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Power Quality and System Stability Impact of Large-Scale Distributed Generation on the Distribution Network: Case Study of 60 MW Derna Wind Farm

By Saad M. Saad, Naser El Naily, Jamal Wafi, Faisal A. Mohamed and Abdelsalam Elhaffar

Related Book

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us