SlideShare a Scribd company logo
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

Analysis of Wind Turbine Driven Permanent Magnet
Synchronous Generator under Different Loading Conditions
Gaber El Saady* , El-Nobi A.Ibrahim, Hamdy Ziedan and Mohammed M.Soliman
Electric Engineering Department, Assiut University, Assiut, Egypt
* E-mail of the corresponding author: gaber1@yahoo.com
Abstract
This paper proposes the configuration of a Wind Turbine generating system equipped with Permanent Magnet
Synchronous Generator (PMSG). There are different types of synchronous generators, but the PMSG is chosen
which has better performance due to higher efficiency and less maintenance. Since it can be used without a gearbox
also implies a reduction of the weight of the nacelle and a reduction of costs. The model includes a wind turbine
model, drive train model and PMSG model. The equations that explain their behavior have been introduced. The
generator model is established in the d-q synchronous rotating reference frame. The proposed Wind Turbine
Generator System (WTGS) has been implemented in MATLAB/SIMULINK software. The PMSG is operating in
stand-alone which is loaded with different types of loads. The simulation results indicate the ability of wind driven
PMSG to operate over wide range of operating conditions at different loading conditions and show effect of different
load types in operation.
Keywords: Permanent Magnet Synchronous Generator (PMSG), Wind Turbine, Wind Energy and WTGS
MATLAB/SIMULINK.
1. Introduction
During the last few decades, wind energy became the most competitive form of clean, non-polluting and renewable
energy to provide a sustainable supply to the world development (Chen et al. 2012) which worldwide wind capacity
doubled approximately every three years. Currently, five countries (Germany, USA, Denmark, India and Spain)
concentrate more than 83% of worldwide wind energy capacity in their countries. The need for increased power
production from the wind and economic reasons, when the rated power of today’s wind turbines is still relatively
small (2MW units are now typical), makes it necessary to group wind turbines into so-called wind farms( Rolan' et al.
2009).Wind farms are built on land, but in recent years there has been(and will be in the future) a strong trend
towards locating them offshore this due to the lack of suitable wind turbine sites on land and the highest wind
speeds located near the sea(and consequently higher energy can be extracted from the wind).
Both induction and synchronous generators can be used for wind turbine systems (Slootweg et al. 2003). Mainly,
three types of induction generators are used in wind power conversion systems: cage rotor, wound rotor with slip
control and doubly fed induction rotors .The last one is the most utilized in wind speed generation because it
provides a wide range of speed variation. However, the variable speed directly driven multi-pole permanent magnet
synchronous generator (PMSG) wind architecture is chosen for this purpose and it is going to be modelled: it offers
better performance due to higher efficiency, simple structure, reliable operation, low noise and less maintenance
because it does not have rotor current. What is more, PMSG can be used without a gearbox, which implies a
reduction of the weight of the nacelle and reduction of costs and so on ( Rolan' et al. 2009).This investigation
presents the model of a PMSG WT able to work under low and fast wind speed conditions and during wind gusts
( López-Ortiz et al. 2012).
The general goal of this paper is to model the electromechanical energy conversion system of Standalone wind
turbine driven with PMSG. Optimum wind energy extraction is achieved by running the Wind Turbine Generator
(WTG) invariable speed because of the higher energy gain and the reduced stresses with using the (PMSG) the
design can be even more simplified.
Simulations have been implemented with the software MATLAB/ SIMULINK to validate the model.
2. System structure and analysis
A typical structure of wind energy conversion system with PMSG consists of a wind turbine, drive train and PMSG.

97
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

2.1. Wind Energy Conversation
The kinetic energy of the wind is given by the following equation (S.VIJAYALAKSHMI et al. 2011):
Ec =

1

(1)

2mv2

m = ρvA∆t

(2)

Where: m is the air mass, v is the wind speed, A is the covered surface of the turbine and 𝜌 is the air density.
The wind power Pw has the following expression:
𝑑

𝑃𝑤 =

𝐸𝑐 =

𝑑𝑡

1

𝜌𝐴𝑣 3

2

(3)

The power coefficient of the turbine Cp can be defined by following equation:
Cp =

Pm

;

Pw

CP < 1

(4)

So the extracted power is given by:
1

Pm = Cp ρπR2 v 3
2

(5)

2

Where: A is area swept by the rotor (A = πR ), R is radius of the turbine rotor and Pm is the mechanical power that
extracts from the wind.
The power coefficient Cp (betz coefficient) reaches maximum value =0.593. In practice, values of obtainable power
coefficient’s are in the range of 45 percent which depends of the tip speed ratio ʎ of the wind turbine and angle of the
blades ß.
Cp = Cp(ʎ, ß)

(6)

The amount of aerodynamic torque is obtained from the power:
𝑇𝑚 =
Substitute from ʎ =

Pm
w

=

1

Cp ρπR2

2

v3

(7)

𝑤

Rw

(8)

v

𝑇𝑚 =
Often the torque coefficient CT =

1

ρπR3 v 2

2

Cp

(9)

ʎ

CP

(10)

ʎ
1

Tm = CT ρπR3 v 2

So,

(11)

2

2.2 The power coefficient of the turbine
The power coefficient can be utilized in the form of look-up tables or in form of a function. The second approach is
presented below (Slootweg et al. 2003), where the power coefficient is defining as a function of the tip-speed ratio
ʎ and the blade pitch angle ß as
𝐶 𝑝 (ʎ, ß) = 0.5 (

116
𝑄

− 0.4ß − 5) 𝑒

21
𝑄

−

(12)

Where Q is represented as
Q=

1
1
0.035
−
ʎ+0.08ß 1+ß3

(13)

2.3 The variation of Cp
The simulation model is shown in Figure 1 which the variation of Cp with tip speed ratio ʎ at different values of
pitch angle ß, also the variation of Cp with pitch angle ß at different values of ʎ are obtained.

98
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

Figure 1. Power Coefficient cp
2.3.1 The variation of CP with the tip speed ratio ʎ for various values of the pitch angle β
The variation of CP with the tip speed ratio ʎ for various values of the pitch angle β is depicted in Figure 2. Thus, by
varying the pitch angle ß, the power coefficient can be changed and the power captured by the turbine can be
controlled.

power coefficent Cp & tip speed ratio curves

power coefficent Cp

0.5

ß
=0
ß
=5
ß
10
ß
=15
ß
=20
ß
=25
ß
=30

0.4
0.3
0.2
0.1
0
0

2

4

6

8

10

12

14

16

tip speed ratio(lambda)
Figure 2. Analytical Approximation of Cp (ʎ, ß) Curves

2.3.2

The variation of CP with pitch angle ß at different values of ʎ

The variation of CP with the pitch angle β for various values of the tip speed ratio ʎ is depicted in Figure 3. Thus, by
varying the tip speed ratio ʎ, the power coefficient can be changed and the power captured by the turbine can be
controlled.

99
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

power coefficent CP & pitch angle B curves

power coefficient CP

0.4

lambda=8
lambda=10

0.3
lambda=6

0.2
lambda=4

0.1

lambda=2

0
0

5

10

15
20
25
pitch angle B [degree]
Figure 3. Power Coefficient of Cp (ʎ, ß) with Pitch Angle ß Curves

30

3. Drive Train Model
The rotational speed is expressed by the following equation (Slootweg et al. 2003):
dwg
dt

=

(Te−Twg −Bm.wg)

(14)

Jeq

Where Sub-index g represents the parameters of the generator side, wg is the mechanical angular speed [rad/s]
of the generator, Bm is the damping coefficient [N.m/s], Te is the electromagnetic torque[N.m], Twg is the
aerodynamic torque that has been transferred to the generator side which ( Twg=Tw/ng )[N.m], Jeq is the equivalent
rotational inertia of the generator [Kg.m2 ].The equivalent rotational inertia of the generator is calculated from this
equation:
Jeq = Jg+

Jw

(15)

n2 g

Where Jg and Jw are the generator and rotor rotational inertia respectively , ng is the gear ratio ( ng=1 , no gear
box).
The parameters used in simulation models are shown in the table (1) :
Table 1: Drive Train Parameters
Parameter
Gear Ratio
Rotational Damping Coefficient

Symbol
Ng
Bm

Value
1
0

Equivalent Inertia

Jeq

0.3kg.m2

The block diagram of drive train model is shown in Figure 4.

100
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

Figure 4.

Drive Train Model

4. Generator Model
The PMSG is used to produce electricity from the mechanical energy obtained from the wind Turbine. In the PMSG,
the rotor magnetic flux is generated by permanent magnets which are placed on the rotor surface ( a non-salient-pole
PMSG)( Wu et al. 2011; Yin et al. 2007).
The main purpose of this case study is to modelling of PMSG from its equations. To simplify the analysis, The
PMSG is normally modeled in the rotor field (dq-axis) synchronous reference Frame, which the q-axis is 90o ahead
of the d-axis with respect to the direction of rotation. The rotor has two axes which the axis that is aligned with the
rotor and flux is called d-axis and the perpendicular axis to d-axis called q-axis ( Abedini 2008).
The flux caused by PM is in the direction of d-axis, the angle between stator axis and d-axis is called ϴe as shown in
Figure 5.

Figure 5. The Configuration of The Winding and PM in The PMSG
The synchronization between the d-q rotating reference frame and the abc-three phase frame is maintained by
utilizing a phase locked loop .To simplify the SG model of Figure 5, the following mathematical manipulations can
be performed.
The voltage equations for the synchronous generator are given by (16) and (17):
vds = −R s ids − wr ʎqs + pʎds
(16)
vqs = −R s iqs + wr ʎds + pʎqs
(17)
Where
ʎds = −Ld ids + ʎr
(18)
ʎqs = −Lq iqs
(19)
Where ʎr is the rotor flux which is constant in the PMSG so,
inductances .
Substitute from equations (18) and (19) in (16) and (17) yield

dʎr
dt

vds = −R s ids + wr Lq iqs − Ld pids

101

= 0 ; Ld and Lq are the stator dq-axis self-

(20)
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

vqs = −R s iqs − wr Ld ids + wr ʎr − Lq piqs

(21)

A simplified dq-axis model of the PMSG in the rotor-field synchronous reference frame is shown in figure 6.

Figure 6. Equivalent Circuit of The PMSG in The Synchronous Frame.
The electromagnetic torque produced by the PMSG can be calculated from the following equations:
3P
Te = (iqs ʎds − ids ʎqs )
(22)
2
3P

= ([iqs ʎr − (Ld − Lq )ids iqs ]
2
The rotor speed wr is governed by motion equation given by:
P
wr = (Te − Tm )

(23)
(24)

JS

To derive the PMSG model for dynamic simulation of synchronous generators equations are rearranged as:
1
ids = (−vds − R s ids + wr Lq iqs )/Ld
(25)
S
1

iqs = (−vqs − R s iqs − wr Ld ids + wr ʎr )/Lq
(26)
S
The active power delivered to the load is given by:
wr
Ps = Pm − Pcus = Tm ∗ wm = Te ∗ − 3Is 2 Rs
(27)
P
Where vds , vqs are the dq-axis stator voltages, R s is the stator resistance , ʎr is the rotor flux linkage, Tm is the
mechanical torque , ids , iqs are the dq-axis stator currents , wr is the rotor mechanical speed and Te is the
electromagnetic torque (Krause 2002 ; Boldea 2006) .
5. Simulation results
The main purpose of this case study is to investigate the operation of a stand-alone PMSG wind energy system
feeding different types of loads .The generator used in the study is a 2.45MW, 4000V, 53.33 HZ,400 rpm non salient
pole PMSG. The parameters of the PMSG used in this paper are shown in table 2.
Table2: The Parameters of The PMSG used in This Paper are listed below:
2.45 MW,4000V,53.33HZ, non-salient pole PMSG Parameters
Parameter

Values and units

Values in Pu

Rated power

2.4487MW

1.0 pu

Rated Apparent power
Rated Line voltage
Rated stator current
Rated rotor speed
Rated Stator Frequency
Number of pole pairs (p)

3.419 MVA
4000 V (rms)
490 A (rms)
400 rpm
53.33 Hz
8

1.0 pu
1.0 pu
1.0 pu
1.0 pu
1.0 pu

Rated mechanical torque

58.4585 KN.m

Rated Rotor Flux Linkage (ʎ )
Stator resistance(Rs)
Stator d-axis inductance(Ld)
Stator q-axis inductance(Lq)

4.971Wb (rms)
24.21 mΩ
9.816mH
9.816mH

0.7213 pu
0.00517 pu
0.7029 pu
0.7029 pu

102
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

5. 1

www.iiste.org

Loading of PMSG with Resistive Load

5.1.1 Loading PMSG with Fixed Resistive Load
The PMSG is loading with a three-phase resistive load R L= 5.5Ω, the block diagram and matlab simulation for the
model is shown in Figures. 7and 8 respectively.

Figure 7. Block Diagram for Simulation

Figure 8. Matlab Simulink Model
The dq-axis stator currents, ids and iqs in the synchronous frame rotating at the synchronous speed of wr are
calculated by the SG model. They are then transformed into the abc -axis stator currents ias, ibs, and ics in the
stationary frame through the dq/abc transformation. The calculated load voltages vas, vbs and vcs which are also
the stator voltages are transformed to the dq-axis voltages vds and vqs in the synchronous frame .These voltages
are then fed back to the SG model.
First , the PMSG is loaded with a three phase balanced resistive load Rl and operates at 320 rpm(.8 pu)at a given
wind speed ( the rotor speed is kept constant at 320 rpm due to assumption that the combined moments of inertia is
very large) .
The following Figure 9. Shows the currents and voltages at PMSG terminals.
ids,iqs,is (pu)

ids,iqs,is (pu)

0.5
0.4

iqs

0.3

is
ids

0.2
0.1
0
0

0.005

0.01

0.015

0.02

0.025

time(sec)
(a)

103

0.03

0.035

0.04

0.045

0.05
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

ias,ibs,ics (pu)

ias,ibs,ics (pu)

0.5

ics

ibs

ias

0

-0.5
0

0.005

0.01

0.015

0.02

0.025
time (sec)

0.03

0.035

0.04

0.045

0.05

0.035

0.04

0.045

0.05

0.035

0.04

0.045

0.05

(b)

vds,vqs,vs (pu)
0.6

vds,vqs,vs (pu)

vs
0.5

vqs

0.4

vds

0.3
0.2
0.1
0
0

0.005

0.01

0.015

0.02

0.025
time(sec)

0.03

(c)

vas,vbs,vcs (pu)
1

vas,vbs,vcs (pu)

vas

vbs vcs

0.5
0
-0.5
-1
0

0.005

0.01

0.015

0.02

0.025

0.03

time(sec)

(d)
Te,Ps (pu)
0.5

Te

Te,Ps (pu)

0.4
0.3

Ps

0.2
0.1
0
0

0.005

0.01

0.015

0.02

0.025
time(sec)

0.03

0.035

0.04

0.045

0.05

(e)

Figure 9. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d) abc voltages, (e)Torque and Power in case of
Resistive Load
From figure 9, the dq-axis stator currents, ids and iqs in the synchronous frame are DC variables, whereas the abcaxis stator currents, ias, ibs, and ics in the stationary frame are sinusoids in steady state. The magnitude of the stator
current ( is) represents the peak value of ias, ibs, and ics is given by equation (28):

104
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

is = √(iqs 2 + ids 2 )
A similar phenomenon can be observed for the stator voltages.
5.1.2

(28)

Decreasing The Load by Half

The generator initially operates in steady state with a resistive load of R l . The load resistance is reduced to R l /2 by
closing switch S at t = 0.0234 sec as shown in figure 10.

Figure 10. SG with a Three-Phase Resistive Load
After a short transient period, the system reaches a new steady-state operating point as shown in Figure 11. From
figures, a decrease in the load resistance results in an increase in the stator currents, but the stator voltages are
reduced mainly due to the voltage drop across the stator inductances. The electromagnetic torque Te and stator active
power Ps are increased accordingly when the system operates at the new operating point.
ids,iqs,is(pu)
0.8
0.7

is
iqs

ids,iqs,is(pu)

0.6
decrease load by half

0.5

ids

0.4
0.3
0.2
0.1
0

0.005

0.01

0.015

0.02

0.025
Time (sec)

0.03

0.035

0.04

0.045

0.05

0.035

0.04

0.045

0.05

(a)
ias,ibs,ics,is(pu)

ias,ibs,ics,is(pu)

1

ias

ibs

is

ics

0.5

0

-0.5

-1
0

0.005

0.01

0.015

0.02

0.025

time(sec)
(b)

105

0.03
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

vds,vqs,vs (pu)
0.6

vds,vqs,vs (pu)

vs
0.5
0.4

vqs

0.3

vds

0.2
0.1
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

0.05

time(sec)
(c)
vas,vbs,vcs,vs (pu)

vas,vbs,vcs,vs (pu)

1

vcs

vbs

vas

vs

0.5
0
-0.5
-1
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

0.05

time(sec)
(d)
Te,Ps (pu)
0.6

Te

Te,Ps (pu)

0.5

Ps

0.4
0.3
0.2
0.1
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

0.05

time(sec)

(e)

Figure 11. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d)abc voltages, (e)Torque and Power in Case of
Decreasing to the Half of The Resistive Load

5.2 PMSG Loaded with (R-L) Load:
The steady-state performance of a stand-alone salient-pole PMSG with an inductive load is analyzed using the dqaxis steady-state equivalent circuit (Wu et al. 2011 ; Krause 2002).The generator operates at the rotor speed of
320 rpm and supplies a three-phase RL load of 𝑅 𝐿 = 4.6797Ω and L = 13.966 mH, Since the q-axis leads the daxis by 90°, The generator dq-axis stator voltages which are also the load voltages can be calculated by following
equations:
vds + jvqs = (ids + jiqs )( R l + jwr Ll )
= (R l ids − wr Ll iqs ) + j(R l iqs + wr Ll ids )
Where wr is the rotor electrical speed which is also the speed of the dq synchronous reference frame.
Equation (29) can be arranged as
vds = R l ids − wr Ll iqs = R l ids − XL iqs

106

(29)

(30)
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

vqs = R l iqs + wr Ll ids = R l iqs + Xl ids

(31)

Where XL iqs = wr Ll iqs , Xl ids = wr Ll ids are referred to speed voltages due to the transformation of the load
inductance from abc-stationary frame to the dq-synchronous frame. A model for RL load is simulating which dq-axis
equivalent circuits of inductive load are shown in Figures. 12 and .The results are shown in Figure 13.

Figure 12. dq-axis Equivalent Circuits

Te,Ps (pu)
0.6

Te

Te,Ps (pu)

0.5

Ps

0.4
0.3
0.2
0.1
0

0.005

0.01

0.015

0.02

ias,ibs,ics (pu)

0.4
0.2

0.025

0.03

0.035

0.04

0.045

0.05

0.035

0.04

0.045

0.05

0.035

0.04

0.045

0.05

time(sec)
(a)
ias,ibs,ics (pu)
ias

ibs

ics

0
-0.2
-0.4
0

0.005

0.01

0.015

0.02

0.025

0.03

time(sec)
(b)
vds,vqs,vs(pu)

vds,vqs,vs(pu)

0.6

vs

vqs

0.4
0.2

vds

0
-0.2
0

0.005

0.01

0.015

0.02

0.025

time(sec)

(c)

107

0.03
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

vas,vbs,vcs(pu)

vas,vbs,vcs(pu)

0.5

vcs

vbs

vas
0

-0.5
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

0.05

time(sec)

(d)
Te,Ps(pu)
0.3

Te

Te,Ps(pu)

Ps
0.2

0.1

0
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

0.05

time(sec)

(e)

Figure 13. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d)abc voltages, (e)Torque and Power in Case of
the Inductive Load
5.3 PMSG Loaded with RC Load

The steady state performance of a stand-alone salient-pole PMSG with a capacitive (RC) load is developed
as inductive (RL) load which is analyzed using the dq-axis steady-state equivalent circuit .The generator
operates at the rotor speed of 320 rpm and supplies a three-phase RC load of 𝑅 𝑙 = 5.5 Ω and C =
637.72 μF. Since the q-axis leads the d-axis by 90° so the generator dq-axis stator voltages which are also
the load voltages can be calculated by following equations:
𝑣 𝑑𝑠 + 𝑗𝑣 𝑞𝑠 = (𝑖 𝑑𝑠 + 𝑗𝑖 𝑞𝑠 )( 𝑅 𝑙 − 𝑗/𝑤 𝐶 𝑙 )
= (𝑅 𝑙 𝑖 𝑑𝑠 +
So,
𝑣 𝑑𝑠 = 𝑅 𝑙 𝑖 𝑑𝑠 −
𝑣 𝑞𝑠 = 𝑅 𝑙 𝑖 𝑞𝑠 −

𝑖 𝑑𝑠
𝑤 𝑟 𝐶𝑙

𝑖 𝑞𝑠
𝑤 𝑟 𝐶𝑙

𝑖 𝑞𝑠
𝑤 𝑟 𝐶𝑙

) + 𝑗(𝑅 𝑙 𝑖 𝑞𝑠 − 𝑖 𝑑𝑠 /𝑤 𝐶 𝑙 )

= 𝑅 𝑙 𝑖 𝑑𝑠 + 𝑋 𝑐 𝑖 𝑞𝑠

= 𝑅 𝑙 𝑖 𝑞𝑠 − 𝑋 𝑐 𝑖 𝑑𝑠

(32)
(33)
(34)

Where𝑋 𝑐 𝑖 𝑞𝑠 = 𝑖 𝑞𝑠 /𝑤 𝐶 𝑙 , 𝑋 𝑐 𝑖 𝑑𝑠 = 𝑖 𝑑𝑠 /𝑤 𝐶 𝑙 are referred to speed voltages due to the transformation of the load
inductance from abc-stationary frame to the dq-synchronous frame.The simulation results are shown in Figure 14.

108
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

ids,iqs,is(pu)

ids,iqs,is(pu)

0.4

is

ids

0.3
0.2

iqs

0.1
0
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

0.035

0.04

0.045

0.05

time(sec)
(a)
ias,ibs,ics(pu)

ias,ibs,ics(pu)

0.4
0.2

ibs

ias

ics

0

-0.2
-0.4
0

0.005

0.01

0.015

0.02

0.025

0.03

0.05

time(sec)
(b)

vds,vqs,vs(pu)
vs

vqs

0.4
0.2

vds

0
-0.2
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.035

0.04

0.045

0.05

time(sec)
(c)
vas,vbs,vcs(pu)
vas

vas,vbs,vcs(pu)

vds,vqs,vs(pu)

0.6

vcs

vbs

0.5

0

-0.5
0

0.005

0.01

0.015

0.02

0.025

time(sec)

(d)

109

0.03

0.045

0.05
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

vas,vbs,vcs(pu)
vas

vcs

vbs

vas,vbs,vcs(pu)

0.5

0

-0.5
0

0.005

0.01

0.015

0.02

0.025

0.03

0.035

0.04

0.045

0.05

time(sec)

(e)

Figure 14. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d) abc voltages, (e)Torque and Power in case of
capacitive load
6. Conclusion



This paper presents the model of stand-alone wind turbine driven permanent magnet synchronous generator
(PMSG) which the model has been implemented in MATLAB / SIMULINK.
The power coefficient ( Cp) curves have been developed and drawn with both tip speed ratio ʎ and pitch
angle ß. These curves show that with varying of the pitch angle ß or tip speed ratio ʎ, the power

coefficient can be changed and the power captured by the turbine also can be controlled.


The PMSG has been modeled in the d-q synchronous rotating reference frame.



The wind driven PMSG is operating stand alone at resistive load firstly where The values of
currents, voltages, torque and Power are obtained. The dq-axis stator currents, ids and iqs in the
synchronous frame are DC variables, whereas the abc-axis stator currents, ias, ibs, and ics in the
stationary frame are sinusoids in steady state.
Then decreasing the load by half, the values of currents, voltages, torque and Power obtained. Due
to a decrease in the load resistance results in an increase in the stator currents, but the stator
voltages are reduced mainly due to the voltage drop across the stator inductances. The
electromagnetic torque Te and stator active power Ps are increased accordingly when the system
operates at the new operating point.






Moreover, wind driven PMSG is loaded with inductive and capacitive loads also.
The results show that the ability of wind driven PMSG to operate over wide range of operating conditions
and at different load changing. Show effect of load type in operation.

References
Junfei Chen, Hongbin Wu, Ming Sun, Weinan Jiang, Liang Cai and Caiyun Guo ,( 2012), "Modeling and Simulation
of Directly Driven Wind Turbine with Permanent Magnet Synchronous Generator," IEEE conference , Innovative
Smart Grid Technologies –Aisa.
Alejandro Rolan', Alvaro Luna, Gerardo Vazquez, Daniel Aguilar and Gustavo Azevedo, (2009)," Modeling of a
Variable Speed Wind Turbine with
Magnet Synchronous Generator," IEEE conference Industrial Electronics,
pp.734 – 739.
J. G. Slootweg, S. W. H. de Haan, H. Polinder and W. L. Kling,( 2003),"General Model for Representing Variable
Speed Wind Turbines in Power System Dynamics Simulations", IEEE Transactions on Power Systems, vol. 18, no.
1,pp. 144 - 151.
E. N. López-Ortiz, D. Campos-Gaona, E.L. Moreno-Goytia,( 2012)," Modelling of a Wind Turbine with Permanent
Magnet Synchronous Generator,"North American Power Symposium,2012,IEEE conference ,pp. 1 - 6,9-11 Sept..
S.VIJAYALAKSHMI,SAIKUMAR.S,SARAVANAN.S,R.V.SANDIP and VIJAY SRIDHAR,( 2011)," Modeling and
control of a Wind Turbine using Permanent Magnet Synchronous Generator, " International Journal of Engineering

110
Innovative Systems Design and Engineering
ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online)
Vol.4, No.14, 2013

www.iiste.org

Science and Technology ,Vol. 3 ,No. 3 ,pp. 2377 - 2384,March.
Bin Wu, Yongqiang Lang, Navid Zargan, and Samir Kouro,( 2011),”Power Conversion and Control of Wind Energy
Systems”, pp.73.
Ming Yin, Gengyin Li, Ming Zhou and Chengyong Zhao,( 2007)," Modeling of the Wind Turbine with Permanent
Magnet Synchronous Generator for Integration", IEEE Conference on power Engineering Society Meeting, pp. 1 - 6 .
Asghar Abedini.( 2008), Integration of permanent magnet synchronous generator wind turbines into power grid.
C. Krause. (2002), Analysis of electric machinery. 2nd Edition. United States of America: Willey.
I. Boldea. (2006), Synchronous Generators. United States of America: Taylor and Francis.

111
This academic article was published by The International Institute for Science,
Technology and Education (IISTE). The IISTE is a pioneer in the Open Access
Publishing service based in the U.S. and Europe. The aim of the institute is
Accelerating Global Knowledge Sharing.
More information about the publisher can be found in the IISTE’s homepage:
http://www.iiste.org
CALL FOR JOURNAL PAPERS
The IISTE is currently hosting more than 30 peer-reviewed academic journals and
collaborating with academic institutions around the world. There’s no deadline for
submission. Prospective authors of IISTE journals can find the submission
instruction on the following page: http://www.iiste.org/journals/
The IISTE
editorial team promises to the review and publish all the qualified submissions in a
fast manner. All the journals articles are available online to the readers all over the
world without financial, legal, or technical barriers other than those inseparable from
gaining access to the internet itself. Printed version of the journals is also available
upon request of readers and authors.
MORE RESOURCES
Book publication information: http://www.iiste.org/book/
Recent conferences: http://www.iiste.org/conference/
IISTE Knowledge Sharing Partners
EBSCO, Index Copernicus, Ulrich's Periodicals Directory, JournalTOCS, PKP Open
Archives Harvester, Bielefeld Academic Search Engine, Elektronische
Zeitschriftenbibliothek EZB, Open J-Gate, OCLC WorldCat, Universe Digtial
Library , NewJour, Google Scholar

More Related Content

PDF
Modeling and performance analysis of a small scale direct driven pmsg based w...
PDF
Analysis of PMSG in Wind Integration using T Source Inverter with Simple Boos...
PPTX
Open-End-Winding Permanent Magnet Synchonous Generator for Wind Energy Conver...
PDF
Wind Energy Conversion System Using PMSG with T-Source Three Phase Matrix Con...
PDF
Permanent magnet Synchronous machines
PDF
Modeling and simulation of pmsm
PDF
SIMULATION AND ANALYSIS OF PERMANENT MAGNET SYNCHRONOUS GENERATOR FOR RENEWAB...
PPTX
Magnitude and frequency control of grid connected doubly fed IG based synchro...
Modeling and performance analysis of a small scale direct driven pmsg based w...
Analysis of PMSG in Wind Integration using T Source Inverter with Simple Boos...
Open-End-Winding Permanent Magnet Synchonous Generator for Wind Energy Conver...
Wind Energy Conversion System Using PMSG with T-Source Three Phase Matrix Con...
Permanent magnet Synchronous machines
Modeling and simulation of pmsm
SIMULATION AND ANALYSIS OF PERMANENT MAGNET SYNCHRONOUS GENERATOR FOR RENEWAB...
Magnitude and frequency control of grid connected doubly fed IG based synchro...

What's hot (20)

PDF
DFIG_report
PDF
doubley fed induction motor
PPTX
Improved reactive power capability with grid connected doubly fed induction g...
PDF
DFIG control of WECS using indirect matrix converter
PPT
Lecture dc machines
PDF
Study of Vector Control Algorithm and Inverter design for BLDC Motor, V/f con...
PPTX
A STAND A
PPTX
Speed Control of Synchronous Motor
PDF
IRJET- Decoupled Control Technique of DFIG with Dual PWM Converters for Wind ...
PPTX
three phase induction motor drive
PDF
Independent Control Of Active And Reactive Powers From DFIG By Logic Fuzzy
PPTX
Presentation1
PDF
Active and Reactive Power control operation of DFIG forWind Power Generation ...
PPTX
Different type s of power converters fed relutance
PDF
Hr3513381342
DOCX
analysis of induction motor drive using slip power recovery scheme
PDF
Starting torque and torque ripple reduction using SVPWM based vector control ...
PDF
Study of Wind Turbine based Variable Reluctance Generator using Hybrid FEMM-M...
DFIG_report
doubley fed induction motor
Improved reactive power capability with grid connected doubly fed induction g...
DFIG control of WECS using indirect matrix converter
Lecture dc machines
Study of Vector Control Algorithm and Inverter design for BLDC Motor, V/f con...
A STAND A
Speed Control of Synchronous Motor
IRJET- Decoupled Control Technique of DFIG with Dual PWM Converters for Wind ...
three phase induction motor drive
Independent Control Of Active And Reactive Powers From DFIG By Logic Fuzzy
Presentation1
Active and Reactive Power control operation of DFIG forWind Power Generation ...
Different type s of power converters fed relutance
Hr3513381342
analysis of induction motor drive using slip power recovery scheme
Starting torque and torque ripple reduction using SVPWM based vector control ...
Study of Wind Turbine based Variable Reluctance Generator using Hybrid FEMM-M...
Ad

Similar to Analysis of wind turbine driven permanent magnet synchronous generator under different loading conditions (20)

PDF
F046013443
PDF
IRJET- Modelling of a PMSG Wind Turbine with Voltage Control
PDF
IRJET- Wind Turbine System Connected to a Permanent Magnet Synchronous Generator
PDF
Power Loss Minimization of Permanent Magnet Synchronous Generator Using Part...
PDF
Fl3610001006
PDF
Modeling of Wind Energy on Isolated Area
PDF
11.modeling and performance analysis of a small scale direct driven pmsg base...
PDF
Performance evolution of a PMSG based WECS using maximum power point tracking...
PDF
paper.pdf
PDF
Performance analysis of various parameters by comparison of conventional pitc...
PDF
Performance analysis of various parameters by comparison of conventional pitc...
PDF
PERMANENT MAGNET SYNCHRONOUS GENERATOR BASED WIND ENERGY CONVERSION SYSTEM
PDF
Paper id 42201607
PDF
Control Scheme for an IPM Synchronous Generator Based-Variable Speed Wind Tur...
PDF
DYNAMIC RESPONSES IMPROVEMENT OF GRID CONNECTED WPGS USING FLC IN HIGH WIND S...
PDF
Australia 1
PDF
Dynamic responses improvement of grid connected wpgs using flc in high wind s...
PDF
Indian journal 3
PPTX
Alternators in windmill
F046013443
IRJET- Modelling of a PMSG Wind Turbine with Voltage Control
IRJET- Wind Turbine System Connected to a Permanent Magnet Synchronous Generator
Power Loss Minimization of Permanent Magnet Synchronous Generator Using Part...
Fl3610001006
Modeling of Wind Energy on Isolated Area
11.modeling and performance analysis of a small scale direct driven pmsg base...
Performance evolution of a PMSG based WECS using maximum power point tracking...
paper.pdf
Performance analysis of various parameters by comparison of conventional pitc...
Performance analysis of various parameters by comparison of conventional pitc...
PERMANENT MAGNET SYNCHRONOUS GENERATOR BASED WIND ENERGY CONVERSION SYSTEM
Paper id 42201607
Control Scheme for an IPM Synchronous Generator Based-Variable Speed Wind Tur...
DYNAMIC RESPONSES IMPROVEMENT OF GRID CONNECTED WPGS USING FLC IN HIGH WIND S...
Australia 1
Dynamic responses improvement of grid connected wpgs using flc in high wind s...
Indian journal 3
Alternators in windmill
Ad

More from Alexander Decker (20)

PDF
Abnormalities of hormones and inflammatory cytokines in women affected with p...
PDF
A validation of the adverse childhood experiences scale in
PDF
A usability evaluation framework for b2 c e commerce websites
PDF
A universal model for managing the marketing executives in nigerian banks
PDF
A unique common fixed point theorems in generalized d
PDF
A trends of salmonella and antibiotic resistance
PDF
A transformational generative approach towards understanding al-istifham
PDF
A time series analysis of the determinants of savings in namibia
PDF
A therapy for physical and mental fitness of school children
PDF
A theory of efficiency for managing the marketing executives in nigerian banks
PDF
A systematic evaluation of link budget for
PDF
A synthetic review of contraceptive supplies in punjab
PDF
A synthesis of taylor’s and fayol’s management approaches for managing market...
PDF
A survey paper on sequence pattern mining with incremental
PDF
A survey on live virtual machine migrations and its techniques
PDF
A survey on data mining and analysis in hadoop and mongo db
PDF
A survey on challenges to the media cloud
PDF
A survey of provenance leveraged
PDF
A survey of private equity investments in kenya
PDF
A study to measures the financial health of
Abnormalities of hormones and inflammatory cytokines in women affected with p...
A validation of the adverse childhood experiences scale in
A usability evaluation framework for b2 c e commerce websites
A universal model for managing the marketing executives in nigerian banks
A unique common fixed point theorems in generalized d
A trends of salmonella and antibiotic resistance
A transformational generative approach towards understanding al-istifham
A time series analysis of the determinants of savings in namibia
A therapy for physical and mental fitness of school children
A theory of efficiency for managing the marketing executives in nigerian banks
A systematic evaluation of link budget for
A synthetic review of contraceptive supplies in punjab
A synthesis of taylor’s and fayol’s management approaches for managing market...
A survey paper on sequence pattern mining with incremental
A survey on live virtual machine migrations and its techniques
A survey on data mining and analysis in hadoop and mongo db
A survey on challenges to the media cloud
A survey of provenance leveraged
A survey of private equity investments in kenya
A study to measures the financial health of

Recently uploaded (20)

PDF
Reconciliation AND MEMORANDUM RECONCILATION
PDF
Business model innovation report 2022.pdf
PDF
Roadmap Map-digital Banking feature MB,IB,AB
PDF
A Brief Introduction About Julia Allison
PDF
MSPs in 10 Words - Created by US MSP Network
PDF
Unit 1 Cost Accounting - Cost sheet
PPT
340036916-American-Literature-Literary-Period-Overview.ppt
PDF
Katrina Stoneking: Shaking Up the Alcohol Beverage Industry
PDF
Nidhal Samdaie CV - International Business Consultant
PDF
WRN_Investor_Presentation_August 2025.pdf
PPTX
Principles of Marketing, Industrial, Consumers,
PPTX
Dragon_Fruit_Cultivation_in Nepal ppt.pptx
DOCX
Business Management - unit 1 and 2
PPTX
Probability Distribution, binomial distribution, poisson distribution
PPTX
Belch_12e_PPT_Ch18_Accessible_university.pptx
PPTX
Amazon (Business Studies) management studies
PDF
Power and position in leadershipDOC-20250808-WA0011..pdf
PPTX
5 Stages of group development guide.pptx
PDF
Laughter Yoga Basic Learning Workshop Manual
PPT
Data mining for business intelligence ch04 sharda
Reconciliation AND MEMORANDUM RECONCILATION
Business model innovation report 2022.pdf
Roadmap Map-digital Banking feature MB,IB,AB
A Brief Introduction About Julia Allison
MSPs in 10 Words - Created by US MSP Network
Unit 1 Cost Accounting - Cost sheet
340036916-American-Literature-Literary-Period-Overview.ppt
Katrina Stoneking: Shaking Up the Alcohol Beverage Industry
Nidhal Samdaie CV - International Business Consultant
WRN_Investor_Presentation_August 2025.pdf
Principles of Marketing, Industrial, Consumers,
Dragon_Fruit_Cultivation_in Nepal ppt.pptx
Business Management - unit 1 and 2
Probability Distribution, binomial distribution, poisson distribution
Belch_12e_PPT_Ch18_Accessible_university.pptx
Amazon (Business Studies) management studies
Power and position in leadershipDOC-20250808-WA0011..pdf
5 Stages of group development guide.pptx
Laughter Yoga Basic Learning Workshop Manual
Data mining for business intelligence ch04 sharda

Analysis of wind turbine driven permanent magnet synchronous generator under different loading conditions

  • 1. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org Analysis of Wind Turbine Driven Permanent Magnet Synchronous Generator under Different Loading Conditions Gaber El Saady* , El-Nobi A.Ibrahim, Hamdy Ziedan and Mohammed M.Soliman Electric Engineering Department, Assiut University, Assiut, Egypt * E-mail of the corresponding author: gaber1@yahoo.com Abstract This paper proposes the configuration of a Wind Turbine generating system equipped with Permanent Magnet Synchronous Generator (PMSG). There are different types of synchronous generators, but the PMSG is chosen which has better performance due to higher efficiency and less maintenance. Since it can be used without a gearbox also implies a reduction of the weight of the nacelle and a reduction of costs. The model includes a wind turbine model, drive train model and PMSG model. The equations that explain their behavior have been introduced. The generator model is established in the d-q synchronous rotating reference frame. The proposed Wind Turbine Generator System (WTGS) has been implemented in MATLAB/SIMULINK software. The PMSG is operating in stand-alone which is loaded with different types of loads. The simulation results indicate the ability of wind driven PMSG to operate over wide range of operating conditions at different loading conditions and show effect of different load types in operation. Keywords: Permanent Magnet Synchronous Generator (PMSG), Wind Turbine, Wind Energy and WTGS MATLAB/SIMULINK. 1. Introduction During the last few decades, wind energy became the most competitive form of clean, non-polluting and renewable energy to provide a sustainable supply to the world development (Chen et al. 2012) which worldwide wind capacity doubled approximately every three years. Currently, five countries (Germany, USA, Denmark, India and Spain) concentrate more than 83% of worldwide wind energy capacity in their countries. The need for increased power production from the wind and economic reasons, when the rated power of today’s wind turbines is still relatively small (2MW units are now typical), makes it necessary to group wind turbines into so-called wind farms( Rolan' et al. 2009).Wind farms are built on land, but in recent years there has been(and will be in the future) a strong trend towards locating them offshore this due to the lack of suitable wind turbine sites on land and the highest wind speeds located near the sea(and consequently higher energy can be extracted from the wind). Both induction and synchronous generators can be used for wind turbine systems (Slootweg et al. 2003). Mainly, three types of induction generators are used in wind power conversion systems: cage rotor, wound rotor with slip control and doubly fed induction rotors .The last one is the most utilized in wind speed generation because it provides a wide range of speed variation. However, the variable speed directly driven multi-pole permanent magnet synchronous generator (PMSG) wind architecture is chosen for this purpose and it is going to be modelled: it offers better performance due to higher efficiency, simple structure, reliable operation, low noise and less maintenance because it does not have rotor current. What is more, PMSG can be used without a gearbox, which implies a reduction of the weight of the nacelle and reduction of costs and so on ( Rolan' et al. 2009).This investigation presents the model of a PMSG WT able to work under low and fast wind speed conditions and during wind gusts ( López-Ortiz et al. 2012). The general goal of this paper is to model the electromechanical energy conversion system of Standalone wind turbine driven with PMSG. Optimum wind energy extraction is achieved by running the Wind Turbine Generator (WTG) invariable speed because of the higher energy gain and the reduced stresses with using the (PMSG) the design can be even more simplified. Simulations have been implemented with the software MATLAB/ SIMULINK to validate the model. 2. System structure and analysis A typical structure of wind energy conversion system with PMSG consists of a wind turbine, drive train and PMSG. 97
  • 2. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org 2.1. Wind Energy Conversation The kinetic energy of the wind is given by the following equation (S.VIJAYALAKSHMI et al. 2011): Ec = 1 (1) 2mv2 m = ρvA∆t (2) Where: m is the air mass, v is the wind speed, A is the covered surface of the turbine and 𝜌 is the air density. The wind power Pw has the following expression: 𝑑 𝑃𝑤 = 𝐸𝑐 = 𝑑𝑡 1 𝜌𝐴𝑣 3 2 (3) The power coefficient of the turbine Cp can be defined by following equation: Cp = Pm ; Pw CP < 1 (4) So the extracted power is given by: 1 Pm = Cp ρπR2 v 3 2 (5) 2 Where: A is area swept by the rotor (A = πR ), R is radius of the turbine rotor and Pm is the mechanical power that extracts from the wind. The power coefficient Cp (betz coefficient) reaches maximum value =0.593. In practice, values of obtainable power coefficient’s are in the range of 45 percent which depends of the tip speed ratio ʎ of the wind turbine and angle of the blades ß. Cp = Cp(ʎ, ß) (6) The amount of aerodynamic torque is obtained from the power: 𝑇𝑚 = Substitute from ʎ = Pm w = 1 Cp ρπR2 2 v3 (7) 𝑤 Rw (8) v 𝑇𝑚 = Often the torque coefficient CT = 1 ρπR3 v 2 2 Cp (9) ʎ CP (10) ʎ 1 Tm = CT ρπR3 v 2 So, (11) 2 2.2 The power coefficient of the turbine The power coefficient can be utilized in the form of look-up tables or in form of a function. The second approach is presented below (Slootweg et al. 2003), where the power coefficient is defining as a function of the tip-speed ratio ʎ and the blade pitch angle ß as 𝐶 𝑝 (ʎ, ß) = 0.5 ( 116 𝑄 − 0.4ß − 5) 𝑒 21 𝑄 − (12) Where Q is represented as Q= 1 1 0.035 − ʎ+0.08ß 1+ß3 (13) 2.3 The variation of Cp The simulation model is shown in Figure 1 which the variation of Cp with tip speed ratio ʎ at different values of pitch angle ß, also the variation of Cp with pitch angle ß at different values of ʎ are obtained. 98
  • 3. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org Figure 1. Power Coefficient cp 2.3.1 The variation of CP with the tip speed ratio ʎ for various values of the pitch angle β The variation of CP with the tip speed ratio ʎ for various values of the pitch angle β is depicted in Figure 2. Thus, by varying the pitch angle ß, the power coefficient can be changed and the power captured by the turbine can be controlled. power coefficent Cp & tip speed ratio curves power coefficent Cp 0.5 ß =0 ß =5 ß 10 ß =15 ß =20 ß =25 ß =30 0.4 0.3 0.2 0.1 0 0 2 4 6 8 10 12 14 16 tip speed ratio(lambda) Figure 2. Analytical Approximation of Cp (ʎ, ß) Curves 2.3.2 The variation of CP with pitch angle ß at different values of ʎ The variation of CP with the pitch angle β for various values of the tip speed ratio ʎ is depicted in Figure 3. Thus, by varying the tip speed ratio ʎ, the power coefficient can be changed and the power captured by the turbine can be controlled. 99
  • 4. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org power coefficent CP & pitch angle B curves power coefficient CP 0.4 lambda=8 lambda=10 0.3 lambda=6 0.2 lambda=4 0.1 lambda=2 0 0 5 10 15 20 25 pitch angle B [degree] Figure 3. Power Coefficient of Cp (ʎ, ß) with Pitch Angle ß Curves 30 3. Drive Train Model The rotational speed is expressed by the following equation (Slootweg et al. 2003): dwg dt = (Te−Twg −Bm.wg) (14) Jeq Where Sub-index g represents the parameters of the generator side, wg is the mechanical angular speed [rad/s] of the generator, Bm is the damping coefficient [N.m/s], Te is the electromagnetic torque[N.m], Twg is the aerodynamic torque that has been transferred to the generator side which ( Twg=Tw/ng )[N.m], Jeq is the equivalent rotational inertia of the generator [Kg.m2 ].The equivalent rotational inertia of the generator is calculated from this equation: Jeq = Jg+ Jw (15) n2 g Where Jg and Jw are the generator and rotor rotational inertia respectively , ng is the gear ratio ( ng=1 , no gear box). The parameters used in simulation models are shown in the table (1) : Table 1: Drive Train Parameters Parameter Gear Ratio Rotational Damping Coefficient Symbol Ng Bm Value 1 0 Equivalent Inertia Jeq 0.3kg.m2 The block diagram of drive train model is shown in Figure 4. 100
  • 5. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org Figure 4. Drive Train Model 4. Generator Model The PMSG is used to produce electricity from the mechanical energy obtained from the wind Turbine. In the PMSG, the rotor magnetic flux is generated by permanent magnets which are placed on the rotor surface ( a non-salient-pole PMSG)( Wu et al. 2011; Yin et al. 2007). The main purpose of this case study is to modelling of PMSG from its equations. To simplify the analysis, The PMSG is normally modeled in the rotor field (dq-axis) synchronous reference Frame, which the q-axis is 90o ahead of the d-axis with respect to the direction of rotation. The rotor has two axes which the axis that is aligned with the rotor and flux is called d-axis and the perpendicular axis to d-axis called q-axis ( Abedini 2008). The flux caused by PM is in the direction of d-axis, the angle between stator axis and d-axis is called ϴe as shown in Figure 5. Figure 5. The Configuration of The Winding and PM in The PMSG The synchronization between the d-q rotating reference frame and the abc-three phase frame is maintained by utilizing a phase locked loop .To simplify the SG model of Figure 5, the following mathematical manipulations can be performed. The voltage equations for the synchronous generator are given by (16) and (17): vds = −R s ids − wr ʎqs + pʎds (16) vqs = −R s iqs + wr ʎds + pʎqs (17) Where ʎds = −Ld ids + ʎr (18) ʎqs = −Lq iqs (19) Where ʎr is the rotor flux which is constant in the PMSG so, inductances . Substitute from equations (18) and (19) in (16) and (17) yield dʎr dt vds = −R s ids + wr Lq iqs − Ld pids 101 = 0 ; Ld and Lq are the stator dq-axis self- (20)
  • 6. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org vqs = −R s iqs − wr Ld ids + wr ʎr − Lq piqs (21) A simplified dq-axis model of the PMSG in the rotor-field synchronous reference frame is shown in figure 6. Figure 6. Equivalent Circuit of The PMSG in The Synchronous Frame. The electromagnetic torque produced by the PMSG can be calculated from the following equations: 3P Te = (iqs ʎds − ids ʎqs ) (22) 2 3P = ([iqs ʎr − (Ld − Lq )ids iqs ] 2 The rotor speed wr is governed by motion equation given by: P wr = (Te − Tm ) (23) (24) JS To derive the PMSG model for dynamic simulation of synchronous generators equations are rearranged as: 1 ids = (−vds − R s ids + wr Lq iqs )/Ld (25) S 1 iqs = (−vqs − R s iqs − wr Ld ids + wr ʎr )/Lq (26) S The active power delivered to the load is given by: wr Ps = Pm − Pcus = Tm ∗ wm = Te ∗ − 3Is 2 Rs (27) P Where vds , vqs are the dq-axis stator voltages, R s is the stator resistance , ʎr is the rotor flux linkage, Tm is the mechanical torque , ids , iqs are the dq-axis stator currents , wr is the rotor mechanical speed and Te is the electromagnetic torque (Krause 2002 ; Boldea 2006) . 5. Simulation results The main purpose of this case study is to investigate the operation of a stand-alone PMSG wind energy system feeding different types of loads .The generator used in the study is a 2.45MW, 4000V, 53.33 HZ,400 rpm non salient pole PMSG. The parameters of the PMSG used in this paper are shown in table 2. Table2: The Parameters of The PMSG used in This Paper are listed below: 2.45 MW,4000V,53.33HZ, non-salient pole PMSG Parameters Parameter Values and units Values in Pu Rated power 2.4487MW 1.0 pu Rated Apparent power Rated Line voltage Rated stator current Rated rotor speed Rated Stator Frequency Number of pole pairs (p) 3.419 MVA 4000 V (rms) 490 A (rms) 400 rpm 53.33 Hz 8 1.0 pu 1.0 pu 1.0 pu 1.0 pu 1.0 pu Rated mechanical torque 58.4585 KN.m Rated Rotor Flux Linkage (ʎ ) Stator resistance(Rs) Stator d-axis inductance(Ld) Stator q-axis inductance(Lq) 4.971Wb (rms) 24.21 mΩ 9.816mH 9.816mH 0.7213 pu 0.00517 pu 0.7029 pu 0.7029 pu 102
  • 7. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 5. 1 www.iiste.org Loading of PMSG with Resistive Load 5.1.1 Loading PMSG with Fixed Resistive Load The PMSG is loading with a three-phase resistive load R L= 5.5Ω, the block diagram and matlab simulation for the model is shown in Figures. 7and 8 respectively. Figure 7. Block Diagram for Simulation Figure 8. Matlab Simulink Model The dq-axis stator currents, ids and iqs in the synchronous frame rotating at the synchronous speed of wr are calculated by the SG model. They are then transformed into the abc -axis stator currents ias, ibs, and ics in the stationary frame through the dq/abc transformation. The calculated load voltages vas, vbs and vcs which are also the stator voltages are transformed to the dq-axis voltages vds and vqs in the synchronous frame .These voltages are then fed back to the SG model. First , the PMSG is loaded with a three phase balanced resistive load Rl and operates at 320 rpm(.8 pu)at a given wind speed ( the rotor speed is kept constant at 320 rpm due to assumption that the combined moments of inertia is very large) . The following Figure 9. Shows the currents and voltages at PMSG terminals. ids,iqs,is (pu) ids,iqs,is (pu) 0.5 0.4 iqs 0.3 is ids 0.2 0.1 0 0 0.005 0.01 0.015 0.02 0.025 time(sec) (a) 103 0.03 0.035 0.04 0.045 0.05
  • 8. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org ias,ibs,ics (pu) ias,ibs,ics (pu) 0.5 ics ibs ias 0 -0.5 0 0.005 0.01 0.015 0.02 0.025 time (sec) 0.03 0.035 0.04 0.045 0.05 0.035 0.04 0.045 0.05 0.035 0.04 0.045 0.05 (b) vds,vqs,vs (pu) 0.6 vds,vqs,vs (pu) vs 0.5 vqs 0.4 vds 0.3 0.2 0.1 0 0 0.005 0.01 0.015 0.02 0.025 time(sec) 0.03 (c) vas,vbs,vcs (pu) 1 vas,vbs,vcs (pu) vas vbs vcs 0.5 0 -0.5 -1 0 0.005 0.01 0.015 0.02 0.025 0.03 time(sec) (d) Te,Ps (pu) 0.5 Te Te,Ps (pu) 0.4 0.3 Ps 0.2 0.1 0 0 0.005 0.01 0.015 0.02 0.025 time(sec) 0.03 0.035 0.04 0.045 0.05 (e) Figure 9. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d) abc voltages, (e)Torque and Power in case of Resistive Load From figure 9, the dq-axis stator currents, ids and iqs in the synchronous frame are DC variables, whereas the abcaxis stator currents, ias, ibs, and ics in the stationary frame are sinusoids in steady state. The magnitude of the stator current ( is) represents the peak value of ias, ibs, and ics is given by equation (28): 104
  • 9. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org is = √(iqs 2 + ids 2 ) A similar phenomenon can be observed for the stator voltages. 5.1.2 (28) Decreasing The Load by Half The generator initially operates in steady state with a resistive load of R l . The load resistance is reduced to R l /2 by closing switch S at t = 0.0234 sec as shown in figure 10. Figure 10. SG with a Three-Phase Resistive Load After a short transient period, the system reaches a new steady-state operating point as shown in Figure 11. From figures, a decrease in the load resistance results in an increase in the stator currents, but the stator voltages are reduced mainly due to the voltage drop across the stator inductances. The electromagnetic torque Te and stator active power Ps are increased accordingly when the system operates at the new operating point. ids,iqs,is(pu) 0.8 0.7 is iqs ids,iqs,is(pu) 0.6 decrease load by half 0.5 ids 0.4 0.3 0.2 0.1 0 0.005 0.01 0.015 0.02 0.025 Time (sec) 0.03 0.035 0.04 0.045 0.05 0.035 0.04 0.045 0.05 (a) ias,ibs,ics,is(pu) ias,ibs,ics,is(pu) 1 ias ibs is ics 0.5 0 -0.5 -1 0 0.005 0.01 0.015 0.02 0.025 time(sec) (b) 105 0.03
  • 10. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org vds,vqs,vs (pu) 0.6 vds,vqs,vs (pu) vs 0.5 0.4 vqs 0.3 vds 0.2 0.1 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 time(sec) (c) vas,vbs,vcs,vs (pu) vas,vbs,vcs,vs (pu) 1 vcs vbs vas vs 0.5 0 -0.5 -1 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 time(sec) (d) Te,Ps (pu) 0.6 Te Te,Ps (pu) 0.5 Ps 0.4 0.3 0.2 0.1 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 time(sec) (e) Figure 11. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d)abc voltages, (e)Torque and Power in Case of Decreasing to the Half of The Resistive Load 5.2 PMSG Loaded with (R-L) Load: The steady-state performance of a stand-alone salient-pole PMSG with an inductive load is analyzed using the dqaxis steady-state equivalent circuit (Wu et al. 2011 ; Krause 2002).The generator operates at the rotor speed of 320 rpm and supplies a three-phase RL load of 𝑅 𝐿 = 4.6797Ω and L = 13.966 mH, Since the q-axis leads the daxis by 90°, The generator dq-axis stator voltages which are also the load voltages can be calculated by following equations: vds + jvqs = (ids + jiqs )( R l + jwr Ll ) = (R l ids − wr Ll iqs ) + j(R l iqs + wr Ll ids ) Where wr is the rotor electrical speed which is also the speed of the dq synchronous reference frame. Equation (29) can be arranged as vds = R l ids − wr Ll iqs = R l ids − XL iqs 106 (29) (30)
  • 11. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org vqs = R l iqs + wr Ll ids = R l iqs + Xl ids (31) Where XL iqs = wr Ll iqs , Xl ids = wr Ll ids are referred to speed voltages due to the transformation of the load inductance from abc-stationary frame to the dq-synchronous frame. A model for RL load is simulating which dq-axis equivalent circuits of inductive load are shown in Figures. 12 and .The results are shown in Figure 13. Figure 12. dq-axis Equivalent Circuits Te,Ps (pu) 0.6 Te Te,Ps (pu) 0.5 Ps 0.4 0.3 0.2 0.1 0 0.005 0.01 0.015 0.02 ias,ibs,ics (pu) 0.4 0.2 0.025 0.03 0.035 0.04 0.045 0.05 0.035 0.04 0.045 0.05 0.035 0.04 0.045 0.05 time(sec) (a) ias,ibs,ics (pu) ias ibs ics 0 -0.2 -0.4 0 0.005 0.01 0.015 0.02 0.025 0.03 time(sec) (b) vds,vqs,vs(pu) vds,vqs,vs(pu) 0.6 vs vqs 0.4 0.2 vds 0 -0.2 0 0.005 0.01 0.015 0.02 0.025 time(sec) (c) 107 0.03
  • 12. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org vas,vbs,vcs(pu) vas,vbs,vcs(pu) 0.5 vcs vbs vas 0 -0.5 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 time(sec) (d) Te,Ps(pu) 0.3 Te Te,Ps(pu) Ps 0.2 0.1 0 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 time(sec) (e) Figure 13. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d)abc voltages, (e)Torque and Power in Case of the Inductive Load 5.3 PMSG Loaded with RC Load The steady state performance of a stand-alone salient-pole PMSG with a capacitive (RC) load is developed as inductive (RL) load which is analyzed using the dq-axis steady-state equivalent circuit .The generator operates at the rotor speed of 320 rpm and supplies a three-phase RC load of 𝑅 𝑙 = 5.5 Ω and C = 637.72 μF. Since the q-axis leads the d-axis by 90° so the generator dq-axis stator voltages which are also the load voltages can be calculated by following equations: 𝑣 𝑑𝑠 + 𝑗𝑣 𝑞𝑠 = (𝑖 𝑑𝑠 + 𝑗𝑖 𝑞𝑠 )( 𝑅 𝑙 − 𝑗/𝑤 𝐶 𝑙 ) = (𝑅 𝑙 𝑖 𝑑𝑠 + So, 𝑣 𝑑𝑠 = 𝑅 𝑙 𝑖 𝑑𝑠 − 𝑣 𝑞𝑠 = 𝑅 𝑙 𝑖 𝑞𝑠 − 𝑖 𝑑𝑠 𝑤 𝑟 𝐶𝑙 𝑖 𝑞𝑠 𝑤 𝑟 𝐶𝑙 𝑖 𝑞𝑠 𝑤 𝑟 𝐶𝑙 ) + 𝑗(𝑅 𝑙 𝑖 𝑞𝑠 − 𝑖 𝑑𝑠 /𝑤 𝐶 𝑙 ) = 𝑅 𝑙 𝑖 𝑑𝑠 + 𝑋 𝑐 𝑖 𝑞𝑠 = 𝑅 𝑙 𝑖 𝑞𝑠 − 𝑋 𝑐 𝑖 𝑑𝑠 (32) (33) (34) Where𝑋 𝑐 𝑖 𝑞𝑠 = 𝑖 𝑞𝑠 /𝑤 𝐶 𝑙 , 𝑋 𝑐 𝑖 𝑑𝑠 = 𝑖 𝑑𝑠 /𝑤 𝐶 𝑙 are referred to speed voltages due to the transformation of the load inductance from abc-stationary frame to the dq-synchronous frame.The simulation results are shown in Figure 14. 108
  • 13. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org ids,iqs,is(pu) ids,iqs,is(pu) 0.4 is ids 0.3 0.2 iqs 0.1 0 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.035 0.04 0.045 0.05 time(sec) (a) ias,ibs,ics(pu) ias,ibs,ics(pu) 0.4 0.2 ibs ias ics 0 -0.2 -0.4 0 0.005 0.01 0.015 0.02 0.025 0.03 0.05 time(sec) (b) vds,vqs,vs(pu) vs vqs 0.4 0.2 vds 0 -0.2 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.035 0.04 0.045 0.05 time(sec) (c) vas,vbs,vcs(pu) vas vas,vbs,vcs(pu) vds,vqs,vs(pu) 0.6 vcs vbs 0.5 0 -0.5 0 0.005 0.01 0.015 0.02 0.025 time(sec) (d) 109 0.03 0.045 0.05
  • 14. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org vas,vbs,vcs(pu) vas vcs vbs vas,vbs,vcs(pu) 0.5 0 -0.5 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 time(sec) (e) Figure 14. (a) dq-axis currents, (b)abc-currents,(c)dq-axis voltages,(d) abc voltages, (e)Torque and Power in case of capacitive load 6. Conclusion   This paper presents the model of stand-alone wind turbine driven permanent magnet synchronous generator (PMSG) which the model has been implemented in MATLAB / SIMULINK. The power coefficient ( Cp) curves have been developed and drawn with both tip speed ratio ʎ and pitch angle ß. These curves show that with varying of the pitch angle ß or tip speed ratio ʎ, the power coefficient can be changed and the power captured by the turbine also can be controlled.  The PMSG has been modeled in the d-q synchronous rotating reference frame.  The wind driven PMSG is operating stand alone at resistive load firstly where The values of currents, voltages, torque and Power are obtained. The dq-axis stator currents, ids and iqs in the synchronous frame are DC variables, whereas the abc-axis stator currents, ias, ibs, and ics in the stationary frame are sinusoids in steady state. Then decreasing the load by half, the values of currents, voltages, torque and Power obtained. Due to a decrease in the load resistance results in an increase in the stator currents, but the stator voltages are reduced mainly due to the voltage drop across the stator inductances. The electromagnetic torque Te and stator active power Ps are increased accordingly when the system operates at the new operating point.    Moreover, wind driven PMSG is loaded with inductive and capacitive loads also. The results show that the ability of wind driven PMSG to operate over wide range of operating conditions and at different load changing. Show effect of load type in operation. References Junfei Chen, Hongbin Wu, Ming Sun, Weinan Jiang, Liang Cai and Caiyun Guo ,( 2012), "Modeling and Simulation of Directly Driven Wind Turbine with Permanent Magnet Synchronous Generator," IEEE conference , Innovative Smart Grid Technologies –Aisa. Alejandro Rolan', Alvaro Luna, Gerardo Vazquez, Daniel Aguilar and Gustavo Azevedo, (2009)," Modeling of a Variable Speed Wind Turbine with Magnet Synchronous Generator," IEEE conference Industrial Electronics, pp.734 – 739. J. G. Slootweg, S. W. H. de Haan, H. Polinder and W. L. Kling,( 2003),"General Model for Representing Variable Speed Wind Turbines in Power System Dynamics Simulations", IEEE Transactions on Power Systems, vol. 18, no. 1,pp. 144 - 151. E. N. López-Ortiz, D. Campos-Gaona, E.L. Moreno-Goytia,( 2012)," Modelling of a Wind Turbine with Permanent Magnet Synchronous Generator,"North American Power Symposium,2012,IEEE conference ,pp. 1 - 6,9-11 Sept.. S.VIJAYALAKSHMI,SAIKUMAR.S,SARAVANAN.S,R.V.SANDIP and VIJAY SRIDHAR,( 2011)," Modeling and control of a Wind Turbine using Permanent Magnet Synchronous Generator, " International Journal of Engineering 110
  • 15. Innovative Systems Design and Engineering ISSN 2222-1727 (Paper) ISSN 2222-2871 (Online) Vol.4, No.14, 2013 www.iiste.org Science and Technology ,Vol. 3 ,No. 3 ,pp. 2377 - 2384,March. Bin Wu, Yongqiang Lang, Navid Zargan, and Samir Kouro,( 2011),”Power Conversion and Control of Wind Energy Systems”, pp.73. Ming Yin, Gengyin Li, Ming Zhou and Chengyong Zhao,( 2007)," Modeling of the Wind Turbine with Permanent Magnet Synchronous Generator for Integration", IEEE Conference on power Engineering Society Meeting, pp. 1 - 6 . Asghar Abedini.( 2008), Integration of permanent magnet synchronous generator wind turbines into power grid. C. Krause. (2002), Analysis of electric machinery. 2nd Edition. United States of America: Willey. I. Boldea. (2006), Synchronous Generators. United States of America: Taylor and Francis. 111
  • 16. This academic article was published by The International Institute for Science, Technology and Education (IISTE). The IISTE is a pioneer in the Open Access Publishing service based in the U.S. and Europe. The aim of the institute is Accelerating Global Knowledge Sharing. More information about the publisher can be found in the IISTE’s homepage: http://www.iiste.org CALL FOR JOURNAL PAPERS The IISTE is currently hosting more than 30 peer-reviewed academic journals and collaborating with academic institutions around the world. There’s no deadline for submission. Prospective authors of IISTE journals can find the submission instruction on the following page: http://www.iiste.org/journals/ The IISTE editorial team promises to the review and publish all the qualified submissions in a fast manner. All the journals articles are available online to the readers all over the world without financial, legal, or technical barriers other than those inseparable from gaining access to the internet itself. Printed version of the journals is also available upon request of readers and authors. MORE RESOURCES Book publication information: http://www.iiste.org/book/ Recent conferences: http://www.iiste.org/conference/ IISTE Knowledge Sharing Partners EBSCO, Index Copernicus, Ulrich's Periodicals Directory, JournalTOCS, PKP Open Archives Harvester, Bielefeld Academic Search Engine, Elektronische Zeitschriftenbibliothek EZB, Open J-Gate, OCLC WorldCat, Universe Digtial Library , NewJour, Google Scholar