Steady State Analysis of PMSG PDF
Steady State Analysis of PMSG PDF
Steady State Analysis of PMSG PDF
www.iiste.org
97
www.iiste.org
(1)
2mv2
m = vAt
(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:
=
(3)
Pm
Pw
CP < 1
(4)
Pm = Cp R2 v 3
2
(5)
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
coefficients 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)
Pm
w
Cp R2
v3
(7)
Rw
(8)
=
Often the torque coefficient CT =
R3 v 2
Cp
(9)
CP
(10)
Tm = CT R3 v 2
So,
(11)
116
0.4 5)
21
(12)
Where Q is represented as
Q=
1
1
0.035
+0.08 1+3
(13)
98
www.iiste.org
power coefficent Cp
0.5
=0
=5
10
=15
=20
=25
=30
0.4
0.3
0.2
0.1
0
0
10
12
14
16
2.3.2
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
www.iiste.org
power coefficient CP
0.4
lambda=8
lambda=10
0.3
lambda=6
0.2
lambda=4
0.1
lambda=2
0
0
10
15
20
25
pitch angle B [degree]
Figure 3. Power Coefficient of Cp (, ) with Pitch Angle Curves
30
(TeTwg 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
100
www.iiste.org
Figure 4.
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.
dr
dt
101
(20)
www.iiste.org
(21)
A simplified dq-axis model of the PMSG in the rotor-field synchronous reference frame is shown in figure 6.
(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
Values in Pu
Rated power
2.4487MW
1.0 pu
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
58.4585 KN.m
4.971Wb (rms)
24.21 m
9.816mH
9.816mH
0.7213 pu
0.00517 pu
0.7029 pu
0.7029 pu
102
5. 1
www.iiste.org
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
www.iiste.org
ias,ibs,ics (pu)
ias,ibs,ics (pu)
0.5
ics
ibs
ias
-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
www.iiste.org
is = (iqs 2 + ids 2 )
A similar phenomenon can be observed for the stator voltages.
5.1.2
(28)
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.
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)
ias
ibs
is
ics
0.5
-0.5
-1
0
0.005
0.01
0.015
0.02
0.025
time(sec)
(b)
105
0.03
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)
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
106
(29)
(30)
www.iiste.org
(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.
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
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
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.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.5
0
0.005
0.01
0.015
0.02
0.025
time(sec)
(d)
109
0.03
0.045
0.05
www.iiste.org
vas,vbs,vcs(pu)
vas
vcs
vbs
vas,vbs,vcs(pu)
0.5
-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. Lpez-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
www.iiste.org
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 IISTEs 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. Theres 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