3.3 POWER TRANSFORMER MODELLING
Power transformers are essential plant components of the power system. In general, they
provide the interface between sections of the network with different rating voltages, for
example a generating plant and the transmission network, a static VAR compensation (SVC)
and the transmission network. Transformers consist of two or three copper windings per
phase and one or more iron cores. They are normally contained in metallic enclosures (i.e.
tanks), and are immersed in highgrade oil for insulation purposes (Grainger and Stevenson,
1994).
From the modelling point of view, it is convenient to separate the electric circuit, formed
by the copper windings, from the magnetic circuit, formed by the iron core. The reactances
of the windings can be found from shortcircuit tests, and the ironcore reactances can be
found from opencircuit tests. The threephase windings of power transformers may be
connected in a number of ways, but in highvoltage transmission the most popular
connections are: (1) star“star, (2) delta“delta, and (3) star“delta. Furthermore, the star
point can be either solidly grounded, grounded through an earthing impedance, or it may be
¬‚oating.
In power transformers the magnetising current usually represents only a small percentage
of the load current. However, this current is rich in harmonics, and a detailed representation
of the magnetic circuit is mandatory in studies involving harmonic frequencies. In
fundamental frequency studies, such as power ¬‚ows, this requirement is not as severe and it
is waved in most cases, unless the study is aimed at conducting an accurate assessment of
power system losses (Acha and Madrigal, 2001; Arrillaga et al., 1997).
3.3.1 Singlephase Transformers
The starting point for developing comprehensive power ¬‚ow transformer models is the
schematic representation of the basic twowinding transformer shown in Figure 3.10. The
windings contain resistance but it is assumed that the core does not saturate and exhibits no
hysteresis.
The two transformer windings, termed primary (p) and secondary (s) windings contain
Np and Ns turns, respectively. The voltages and currents existing in both windings are related
by a matrix of shortcircuit (sc) impedance parameters, as given by the following
71
POWER TRANSFORMER MODELLING
Is
Ip
1 3
+
+
Vs
Vp


2 4
Twowinding transformer
Figure 3.10
expression:
! ! !
Vp Zsc p Zsc m Ip
¼ °3:68Þ
Vs Zsc m Zsc s Is
where
9
Zsc p ¼ Rp þ jXsc p ; >
=
Zsc s ¼ Rs þ jXsc s ; °3:69Þ
>
;
Zsc m ¼ jXsc m :
These parameters are obtained by measurements of the actual transformer, where Rp and Rs
are the resistances of the primary and secondary windings, respectively. The reactances
Xsc p ; Xsc s , and Xsc m are shortcircuit reactances obtained by exciting two terminals of the
transformer shown in Figure 3.10, at reduced voltage, and short circuiting the other two. The
ratio of excitation voltage to shortcircuit current gives the relevant shortcircuit reactance:
E12 E34 E13
Xsc p ¼ ; Xsc s ¼ ; Xsc m ¼ : °3:70Þ
I34 I12 I24
From the point of view of system analysis there are advantages in expressing the short
circuit impedance matrix of Equation (3.68) in admittance form:
! ! !
ÀYsc m Vp
Ip Ysc p
¼ ; °3:71Þ
ÀYsc m
Is Ysc s Vs
where
9
Zsc s
;>>
¼ >
Ysc p
Zsc p Zsc s À Zsc m > >
2
>
>
=
Zsc p
¼ ; °3:72Þ
Ysc s
Zsc p Zsc s À Zsc m >
2
>
>
>
>
>
Zsc m
:>
¼ ;
Ysc m
Zsc p Zsc s À Zsc m
2
It is observed that up to three shortcircuit tests may be required to characterise the matrices
of shortcircuit parameters. However, the primary and secondary shortcircuit admittances
are almost the same when expressed in perunit values, say Ysc. Owing to the strong
magnetic coupling afforded by iron cores, the mutual admittance between primary and
72 MODELLING OF CONVENTIONAL POWER PLANT
secondary windings can also be taken to have a value of Ysc . Hence, the transformer short
circuit admittance matrix, in per units, is:
! ! !
Ysc ÀYsc Vp
Ip
¼ : °3:73Þ
ÀYsc Ysc Vs
Is
3.3.2 Simple Tapchanging Transformer
The effect of expressing the transformer parameters in the perunit system is to transform
the original voltage ratio Np : Ns into a unity voltage ratio 1 : 1. This enables a simple
equivalent circuit consisting of the shortcircuit admittance Ysc connected between the
primary bus (p) and the secondary bus (s) to describe adequately the system performance of
the twowinding transformer.
However, power transformers are often ¬tted with a tapchanging mechanism to enable a
degree of voltage magnitude regulation at one of the transformer terminals. This is achieved
by injecting a small variable voltage magnitude in phase (added or subtracted) with the
voltage magnitude at the output winding. Such transformers are termed load tapchanging
(ltc) transformers and play an important role in power ¬‚ow studies. The representation of an
ltc transformer may be achieved by the series connection of the shortcircuit admittance
representing a perunit transformer and an ideal transformer with taps ratio T : 1 (Laughton,
1968). This arrangement is shown in Figure 3.11.
Vp Vs
Ysc
Ip Is
I′
V
I
T :1
Simple tapchanging transformer
Figure 3.11
The following relationships exist in the ideal transformer,
T I0
V T
¼; ¼: °3:74Þ
and
Vs 1 1 I
The current across the admittance Ysc is:
À Á À Á
I ¼ Ysc Vp À V ¼ Ysc Vp À TVs ¼ Ip : °3:75Þ
Also,
À Á
I 0 ¼ TI ¼ Ysc TVp À T 2 Vs ¼ ÀIs : °3:76Þ
Combining Equations (3.75) and (3.76) in matrix form gives:
! ! !
ÀTYsc Vp
Ip Ysc
¼ : °3:77Þ
ÀTYsc T 2 Ysc
Is Vs
73
POWER TRANSFORMER MODELLING
3.3.3 Advanced Tapchanging Transformer
Following the same line of reasoning, a comprehensive power system transformer model is
´
derived for a singlephase threewinding transformer (Acha, AmbrizPerez, and Fuerte
Esquivel, 2000). Each winding is represented as the series combination of a shortcircuit
admittance and an ideal transformer. Furthermore, each winding is provided with a complex
tapchanging mechanism to allow for tapchanging and phaseshifting facilities. Moreover,
the magnetising branch of the transformer is included to account for the core losses.
Figure 3.12 shows the equivalent circuit of the threewinding transformer.
1:Uv
Vs
1:Ui
I2′ Ysc s
V2 I2 Is
Vp
Ysc p I ′
Ip I1 V0
V1
1
V3 I3 Ysc t
I3′
I0 It
Tv:1
Ti:1 1:Wv
Vt
1:Wi
B0
G0
Figure 3.12 Comprehensive tapchanging transformer. Reproduced, with permission, from E. Acha,
´
H. AmbrizPerez, and C.R. FuerteEsquivel, ˜Advanced Transformer Control Modelling in an Optimal
Power Flow using Newton™s Method™, IEEE Trans. Power Systems 15(1) 290“298, # 2000 IEEE
The primary winding is represented as an ideal transformer having complex tap ratios
Tv : 1 and Ti : 1 in series with the admittance Ysc p , where Tv ¼ TiÃ ¼ t þ j ¼ T¬t . The
symbol Ã denotes the conjugate operation. The secondary winding is represented as an ideal
transformer having complex tap ratios Uv : 1 and Ui : 1 in series with the admittance Ysc s ,
where Uv ¼ UiÃ ¼ u þ j ¼ U¬u . Similarly, the ideal transformer in the tertiary winding
has complex tap ratios Wv : 1 and Wi : 1 in series with an admittance Ysc t , where
Wv ¼ WiÃ ¼ w þ j
¼ W¬w . It is assumed here that Ysc p , Ysc s and Ysc t are the shortcircuit
admittances of the primary, secondary, and tertiary windings, respectively. The magnetising
branch of the transformer is represented by the admittance Y0 ¼ G0 þ jB0 .
The resistive path of the magnetising branch is directly related to the iron losses, and its
conductance G0 draws a current that varies linearly with the voltage across the magnetising
branch. However, in the inductive path the relationship between the current and the voltage
is dictated by the rms V“I characteristic, which under saturating conditions becomes non
linear.
74 MODELLING OF CONVENTIONAL POWER PLANT
The following relationships exist in the ideal primary, secondary, and tertiary
transformers:
0
V1 Tv Ti I1
¼ ; and ¼; °3:78Þ
V0 1 1 I1
0
V2 Uv Ui I2
¼ ; ¼; °3:79Þ
and
V0 1 1 I2
0
V3 Wv Wi I3
¼ ; ¼: °3:80Þ
and
V0 1 1 I3
The currents across the admittances Ysc p, Ysc s, and Ysc t are, respectively,
À Á À Á
I1 ¼ Ysc p Vp À V1 ¼ Ysc p Vp À Tv V0 ¼ Ip ; °3:81Þ
I2 ¼ Ysc s °Vs À V2 Þ ¼ Ysc s °Vs À Uv V0 Þ ¼ Is ; °3:82Þ
I3 ¼ Ysc t °Vt À V3 Þ ¼ Ysc t °Vt À Wv V0 Þ ¼ It ; °3:83Þ
and at the centre of the transformer the following relationship holds:
0 0 0
0 ¼ I1 þ I2 þ I3 À I0 ¼ Ti I1 þ Ui I2 þ Wi I3 À I0 : °3:84Þ
Substituting Equations (3.81)“(3.83) into Equation (3.84) gives:
À2 Á
Ã Ã Ã
0 ¼ ÀTv Ysc p Vp À Uv Ysc s Vs À Wv Ysc t Vt þ Tv Ysc p þ Uv Ysc s þ Wv Ysc t þ Y0 V0 : °3:85Þ
2 2
Putting Equations (3.81)“(3.83) and (3.85) in matrix form gives:
2 32 32 3
ÀTv Ysc p
Ysc p 0 0
Ip Vp
6 Is 7 6 76 Vs 7
ÀUv Ysc s
0 Ysc s 0
6 7¼6 76 7:
4 It 5 4 54 Vt 5
ÀWv Ysc t
0 0 Ysc t
Ã Ã Ã
ÀTv Ysc p ÀUv Ysc s ÀWv Ysc t Tv Ysc p þ Uv Ysc s þ Wv Ysc t þ Y0
2 2 2
0 V0
°3:86Þ
Equation (3.86) represents the transformer shown in Figure 3.12. However, it is possible
to ¬nd a reduced equivalent matrix that still models the transformer correctly while retaining
only the external buses p, s, and t. This is done by means of Gaussian elimination:
2 3 2 Ã
Uv Ysc p Ysc s þ Wv Ysc p Ysc t þ Ysc p Y0 ÀTv Uv Ysc p Ysc s
2 2
Ip
6 7 16 Ã
4 Is 5 ¼ 4 ÀTv Uv Ysc p Ysc s Tv Ysc s Ysc p þ Wv Ysc s Ysc t þ Ysc s Y0
2 2
D
Ã Ã
ÀTv Wv Ysc p Ysc t ÀUv Wv Ysc s Ysc t
It
32 3
Ã
ÀTv Wv Ysc p Ysc t Vp
76 7
Ã
ÀUv Wv Ysc s Ysc t 54 V s 5; °3:87Þ
Tv Ysc t Ysc p þ Uv Ysc t Ysc s þ Ysc t Y0
2 2
Vt
where
D ¼ Tv Ysc p þ Uv Ysc s þ Wv Ysc t þ Y0 :
2 2 2
75
POWER TRANSFORMER MODELLING
The nodal admittance representation of a twowinding transformer can be easily obtained
by introducing simplifying assumptions in Equation (3.87). For instance, when the tertiary
winding does not exist, the row and column corresponding to this bus become redundant and
they are removed from Equation (3.87). Moreover, the tap ratios Wv and Wi become zero.
Hence, the nodal admittance matrix equation representing the twowinding transformer is
arrived at:
! ! !
Ã
1 Uv Ysc p Ysc s þ Ysc p Y0 ÀTv Uv Ysc p Ysc s
2
Ip Vp
¼ ; °3:88Þ
Ã
ÀTv Uv Ysc p Ysc s Tv Ysc s Ysc p þ Ysc s Y0 Vs
D 2
Is
where
D ¼ Tv Ysc p þ Uv Ysc s þ Y0 :
2 2
It must be noted that owing to the ¬‚exibility of the twowinding transformer model in
Equation (3.88), it is possible to assemble a transformer model that represents the
transformer circuit shown in Figure 3.12 by using three of these twowinding transformer
models. An example of how this can be achieved is shown elsewhere (see Acha, Ambriz
´
Perez, and FuerteEsquivel, 2000).
Transformer models with more constrained tapping arrangements can also be derived
from Equation (3.88). For instance, take the case of the tapchanging transformer shown in
Figure 3.11, represented by Equation (3.77). Such a representation can be derived from
Equation (3.88) by including no magnetising branch, Y0 ¼ 0, and a nominal tapping position
for the secondary winding, Uv ¼ 1. Moreover, the tapping position of the primary winding
is real as opposed to complex, Tv ¼ T, and the shortcircuit admittance is assumed to be all
on the primary side, Ysc s ¼ 0 and Ysc p ¼ Ysc . The latter consideration requires application of
L™Hopital differentiation rule with respect to Ysc s .
A further strength of the transformer model in Equation (3.88) is that, owing to the
complex nature of their taps, it represents rather well the system behaviour of a phase
shifting (PS) twowinding transformer. This is more easily appreciated if it is assumed that
in Equation (3.88) both complex taps have unit magnitudes:
)
Tv ¼ 1¬t ¼ cos t þ j sin t ;
°3:89Þ
Uv ¼ 1¬v ¼ cos v þ j sin v :
Hence,
"#
Ysc p °Ysc s þ Y0 Þ
Ip 1
¼
Ysc p þ Ysc s þ Y0 ÀYsc p °cos t À j sin t Þ ‚ Ysc s °cos u þ j sin u Þ
Is
! !
ÀYsc p °cos t þ j sin t Þ ‚ Ysc s °cos u À j sin u Þ Vp
À Á : °3:90Þ
Ysc s Ysc p þ Y0 Vs
This is a comprehensive model of a PS transformer that yields very ¬‚exible power ¬‚ow and
optimal power ¬‚ow PS models, as will be shown in Chapters 4 and 7, respectively.
3.3.4 Threephase Transformers
Based on nodal analysis, quite general models for multiwindings, multiphase transformers
can be derived (Chen and Dillon, 1974; Laughton, 1968). The essence of the method is to
transform the shortcircuit parameters of the transformer windings, suitably arranged in a
76 MODELLING OF CONVENTIONAL POWER PLANT
matrix of primitive parameters Y , into nodal parameters Y . This is done with the help
of appropriate connectivity matrices, namely, C and C . The connectivity matrices relate
the voltages and currents in the unconnected transformer windings to the phase voltages and
currents when the threephase transformer is actually connected.
The primitive and nodal parameters are related by the following matrix expression:
Y ¼ C Y C : °3:91Þ
The primitive parameters of three identical singlephase transformers, for which the
terminals between transformers are not connected in any way but contain offnominal
tapping facilities on the primary winding, have the following arrangement:
! ! ! 9
ÀTv Ysc
Ysc
I1 V1 >
>
¼ ; >
>
Ã >
ÀTv Ysc 2
Tv Ysc >
I2 V2
>
2 3 2 32 3 >
>
ÀTv Ysc >
Ysc 0 0 0 0
I1 V1 >
>
>
6 I 7 6 ÀT Ã Ysc 76 V 7 >
>
2
6 27 6 v 76 2 7 >
Tv Ysc 0 0 0 0
>
>
676 76 7
! ! ! =
6 I3 7 6 0 76 V3 7
ÀTv Ysc ÀTv Ysc
Ysc 0 Ysc 0 0
I3 V3
) 6 7¼6 76 7 ;
¼ 6I 7 6 0 0 76 V4 7 >
Ã Ã
ÀTv Ysc ÀTv Ysc
2 2
76 7 >
6 47 6
Tv Ysc 0 Tv Ysc 0
I4 V4
76 7 > >
676
ÀTv Ysc 54 V5 5 >
>
4 I5 5 4 0 >
0 0 0 Ysc >
>
V6 > >
Ã
ÀTv Ysc >
2
0 0 0 0 Tv Ysc
I6 >
>
! ! ! >
>
ÀTv Ysc >
Ysc
I5 V5 >
;
¼ :
Ã
ÀTv Ysc 2
Tv Ysc
I6 V6
°3:92Þ
In general, these matrix equations may be expressed in compact form:
I ¼Y V : °3:93Þ
The three singlephase transformers, when suitably connected, electrically speaking, may
serve the purpose of transforming threephase voltages and currents. The assembly is termed
a ˜threephase bank™. Each singlephase unit in the bank is closely associated with one phase
of the threephase system. Depending on the electrical connection and operating conditions,
there may be currents from more than one phase circulating in one singlephase unit at any
one time, but there are not ¬‚ux interactions between windings of different units.
Quite a different situation prevails in multilimb transformers, where all windings of the
threephase unit are magnetically coupled. The primitive admittance matrix equation of
the twowinding, threephase transformers is a full matrix, and up to 21 shortcircuit tests
may be required to de¬ne fully this primitive admittance matrix. In the remainder of this
chapter only the threephase bank will be addressed.
The three most popular threephase transformer connections found in highvoltage
transmission are addressed below, namely the star“star, delta“delta, and star“delta. To
determine their nodal admittance matrix models, one requires information of the matrix of
primitive parameters, Y , and the relevant connectivity matrices, C and C .
3.3.4.1 Star“star connection
The threephase connection is shown in Figure 3.13 when the windings are connected in star“
star con¬guration, with both star points grounded through admittances, YN and Yn, respectively.
77
POWER TRANSFORMER MODELLING
Ia
IA
VA Va
Ysc
Ysc
Ib
Ysc
IB Ysc
VB Vb
Ysc Ysc
Ic
IC
VC Vc
YN Yn
Star“star connection
Figure 3.13
The transformation matrix, which relates the voltages existing in the unconnected
transformer to the voltages in the connected threephase transformer shown in Figure 3.13,
is given explicitly in Equation (3.94):
2 32 32 3
0 À1 VA
V1 1 0 0 0 0 0
6 V2 7 6 0 À1 76 VB 7
0 0 1 0 0 0
676 76 7
6 V3 7 6 0 0 76 V C 7
0 À1
1 0 0 0
676 76 7
6 V4 7 6 0 À1 76 Va 7
0 0 0 1 0 0
6 7¼6 76 7: °3:94Þ
6 V5 7 6 0 0 76 V b 7
0 À1
0 1 0 0
676 76 7
6 V6 7 6 0 À1 76 Vc 7
0 0 0 0 1 0
676 76 7
4 V7 5 4 0 0 54 V N 5
0 0 0 0 0 1
V8 0 0 0 0 0 0 0 1 Vn
In compact form, we have,
V ¼ C V : °3:95Þ
The nodal matrix representation of this transformer connection is obtained by substituting
Equations (3.92) and (3.94) into Equation (3.91):
2 32 32 3
ÀTv Ysc ÀYsc
Ysc 0 0 0 0 Tv Ysc
IA VA
6 IB 7 6 0 76 VB 7
ÀTv Ysc ÀYsc
Ysc 0 0 0 Tv Ysc
676 76 7
6 IC 7 6 0 76 VC 7
ÀTv Ysc ÀYsc
0 Ysc 0 0 Tv Ysc
676 76 7
6 Ia 7 6 ÀTv Ysc ÀTv Ysc 76 Va 7
Ã Ã
2 2
0 0 Tv Ysc 0 0 Tv Ysc
6 7¼6 76 7
6 Ib 7 6 0 ÀTv Ysc 76 Vb 7
Ã Ã
ÀTv Ysc 2 2
0 0 Tv Ysc 0 Tv Ysc
676 76 7
6 Ic 7 6 0 ÀTv Ysc 76 Vc 7
Ã Ã
ÀTv Ysc 2 2
0 0 0 Tv Ysc Tv Ysc
676 76 7
4 IN 5 4 ÀYsc À3Tv Ysc 54 VN 5
ÀYsc ÀYsc Tv Ysc 3Ysc þ YN
Tv Ysc Tv Ysc
Ã Ã Ã Ã
Tv Ysc ÀTv Ysc ÀTv Ysc ÀTv Ysc À3Tv Ysc 3Tv Ysc þ Yn
2 2 2 2
Tv Ysc Tv Ysc
In Vn
°3:96Þ
If both star points N and n are solidly grounded then the nodal voltages VN and Vn become
zero. Hence, the rows and columns corresponding to bus N and bus n become redundant and
78 MODELLING OF CONVENTIONAL POWER PLANT
are deleted from matrix Equation (3.96):
232 32 3
ÀTv Ysc
Ysc 0 0 0 0
IA VA
6 IB 7 6 0 7 6 VB 7
ÀTv Ysc
Ysc 0 0 0 76 7
676
6 IC 7 6 0 ÀTv Ysc 76 VC 7
0 Ysc 0 0
6 7¼6 76 7: °3:97Þ
6 Ia 7 6 ÀT Ysc 0 7 6 Va 7
Ã 2
0 0 Tv Ysc 0
676v 76 7
4 Ib 5 4 0 0 5 4 Vb 5
Ã
ÀTv Ysc 2
0 0 Tv Ysc
Ã
ÀTv Ysc 2
0 0 0 0 Tv Ysc
Ic Vc
3.3.4.2 Delta“delta connection
This transformer connection is shown in Figure 3.14. In the delta connection the following
relationships exist between the voltages and currents in the connected and unconnected
circuits:
1
V ¼ p¬¬¬ C V ; °3:98Þ
3
p¬¬¬
3I ¼ C I : °3:99Þ
Ia
IA
VA Va
Ib
IB Ysc Ysc Ysc Ysc
VB Vb
Ic
IC Ysc
Ysc
VC Vc
Delta“delta connection
Figure 3.14
The relevant connectivity matrices for this transformer connection are set up and, upon
substitution in Equation (3.91), the following nodal admittance matrix is arrived at:
23 2 32 3
ÀYsc ÀYsc À2Tv Ysc
2Ysc Tv Ysc Tv Ysc
IA VA
6 IB 7 6 ÀYsc Tv Ysc 76 VB 7
ÀYsc À2Tv Ysc
2Ysc Tv Ysc
67 6 76 7
6 I 7 1 6 ÀY À2Tv Ysc 76 VC 7
ÀYsc
6 C7 6 76 7
2Ysc Tv Ysc Tv Ysc
sc
6 7¼ 6 76 7:
Ã Ã Ã
6 Ia 7 3 6 À2Tv Ysc ÀTv Ysc ÀTv Ysc 76 Va 7
2 2 2
Tv Ysc Tv Ysc 2Tv Ysc
67 6Ã 76 7
4 Ib 5 4 T Ysc ÀTv Ysc 54 Vb 5
Ã Ã
À2Tv Ysc ÀTv Ysc
2 2 2
Tv Ysc 2Tv Ysc
v
Ã Ã Ã
À2Tv Ysc ÀTv Ysc ÀTv Ysc
2 2 2
Tv Ysc Tv Ysc 2Tv Ysc
Ic Vc
°3:100Þ
3.3.4.3 Star“delta connection
This transformer connection is shown in Figure 3.15 for the case when the star point is
solidly grounded. Following a similar procedure to that used to derive the nodal admittance
79
POWER TRANSFORMER MODELLING
Ia
IA
Va
VA
Ysc
Ib
Ysc
IB Ysc Ysc
Vb
VB Ysc
Ic
Ysc
IC
Vc
YN
VC
Star“delta connection
Figure 3.15
matrices of the star“star and delta“delta connections, the nodal matrix representation of this
transformer connection is:
p¬¬¬ p¬¬¬
2 32 32 3
ÀTv Ysc
Ysc 0 0 3 Tv Ysc 3 0
IA VA
p¬¬¬ p¬¬¬ 76 7
676
ÀTv Ysc 3
6 IB 7 6 Tv Ysc 3 76 VB 7
0 Ysc 0 0
676 p¬¬¬ 76 7
p¬¬¬
6I 7 6 ÀTv Ysc 3 76 VC 7
6 C7 6 76 7
0 0 Ysc Tv Ysc 3 0
6 7¼6 p¬¬¬ p¬¬¬ 76 7:
6 Ia 7 6 ÀTv Ysc 3 ÀTv Ysc 3 76 Va 7
Ã Ã
ÀTv Ysc 3
2 2 2
0 Tv Ysc 3 2Tv Ysc 3
676 76 7
6 I 7 6 Ã p¬¬¬ p¬¬¬ 76 7
Ã
4 b 5 4 Tv Ysc 3 ÀTv Ysc 3 ÀTv Ysc 3 ÀTv Ysc 3 54 Vb 5
2 2 2
0 2Tv Ysc 3
p¬¬¬ p¬¬¬
Ã Ã
Tv Ysc 3 ÀTv Ysc 3 ÀTv Ysc 3 ÀTv Ysc 3
2 2 2
Ic Vc
0 2Tv Ysc 3
°3:101Þ
3.3.5 Sequence Domain Parameters
Transformer parameters are also amenable to representation in the frame of reference of the
sequences (Chen and Dillon, 1974). The matrix of symmetrical components and its inverse,
given in Equations (3.57), are used to such effect. This requires that the order of all matrices
involved in the exercise be a multiple of three. This characteristic is met by matrices
representing the star“star connected transformer with both star points solidly grounded, the
delta“delta transformer, and the star“delta transformer with the star point solidly grounded.
It should be noted that the symmetrical components transform given in Equations (3.57)
cannot directly be applied to cases of starconnected windings, where one or two star points
are not grounded or are grounded through earthing impedances. In such cases, Kron™s
reductions are applied ¬rst to ¬nd out reduced equivalent representations which are a
function only of phase terminals. This follows the spirit of the procedure presented in
Section 3.2.2 for the elimination of transmission line ground wires.
A generic, compact representation of Equations (3.97), (3.100), and (3.101) correspond
ing to the star“star, delta“delta, and star“delta connections may be expressed as:
! ! !
ÀYII VABC
IABC YI
¼ ; °3:102Þ
ÀYIII
Iabc YIV Vabc
80 MODELLING OF CONVENTIONAL POWER PLANT
where the order of matrices YI , YII , YIII , YIV is 3 ‚ 3 and suitable for direct treatment by
the matrix of symmetrical components, to enable representation in the frame of reference of
the sequences. This is achieved by applying the following symmetrical component
operations:
9
Yi ¼ TÀ1 YI TS ; > >
>
S
>
Yii ¼ TS YII TS ; =
À1
°3:103Þ
Yiii ¼ TÀ1 YIII TS ; >
>
>
>
S
;
À1
Yiv ¼ TS YIV TS :
Table 3.1 shows matrices Yi , Yii , Yiii , and Yiv in explicit form, for the star“star, delta“delta,
and star“delta transformer connections.
Table 3.1 Transformer sequence domain admittances
Matrix type Star“star Delta“delta Star“delta
2 3 2 3 2 3
Ysc 0 0 00 0 Ysc 0 0
40 05 4 0 Ysc 05 40 05
Ysc Ysc
Yi
0 0 Ysc 00 Ysc 0 0 Ysc
2 3 2 3 2 3
Tv Ysc 0 0 0 0 0 0 0 0
40 05 40 05 40 5
Tv Ysc ¬30
Tv Ysc Tv Ysc 0
Yii
Tv Ysc ¬ À30
0 0 Tv Ysc 0 0 Tv Ysc 0 0
2 3 2 3 2 3
Ã
Tv Ysc 0 0 0 0 0 0 0 0
40 05 4 0 Tv Ysc 05 4 0 Tv Ysc ¬ À 30 5
Ã Ã Ã
Tv Ysc 0
Yiii
Ã Ã
Tv Ysc ¬30
Ã
0 0 Tv Ysc 0 0 Tv Ysc 0 0
2 3 2 3 2 3
2 0 0 0 0 0 0
Tv Ysc 0 0
40 05 4 0 Tv Ysc 05 4 0 Tv Ysc 05
2 2
2
Yiv Tv Ysc
2 2
2 0 0 Tv Ysc 0 0 Tv Ysc
0 0 Tv Ysc
The sequence domain representation of a transformer, in compact form, is:
!" # !
°iÞ °iiÞ
ÀYi
I012p V012p
Yi
¼ ; °3:104Þ
°iiiÞ °ivÞ
ÀY
I012s V012s
Yi
i
where the subscripts 0, 1, and 2 refer to zero, positive, and negative sequence quantities,
respectively. It has been emphasised in various points.
Careful examination of the sequence domain parameters indicates that three independent
transfer admittance matrix equations, leading to three independent circuits, are generated for
a threephase transformer. This is more easily realised if the transformer taps are taken to be
81
POWER TRANSFORMER MODELLING
real as opposed to complex, yielding symmetrical matrix equations and, hence, reciprocal
circuits. The star“star, delta“delta, and star“delta connections share the same positive and
negative sequence equivalent circuits, given in Figure 3.16.
Ip Is
TYsc
(T2T)Ysc
Vp Vs
(1T)Ysc
Figure 3.16 Positive and negative sequence equivalent circuit for the star“star, delta“delta, and star“
delta connections
In contrast, the zero sequence equivalent circuits for the three connections differ from one
another. The equivalent circuits are shown in Figures 3.17(a), 3.17(b), and 3.17(c) for the
star“star, delta“delta, and star“delta connections, respectively.
I0p
I0p I0s Ysc
TYsc
(1’T)Ysc V0p
V0p V0s
V0s
V0p V0s
(T’T2)Ysc
(c)
(a) (b)
Zero sequence equivalent circuits for: (a) star“star, (b) delta“delta, and (c) star“delta
Figure 3.17
connections
It should be noted that for the star“delta transformer connection the primary and
secondary terminals of the zero sequence equivalent circuit are not electrically connected.
However, the primary terminal contains an admittance Ysc connected between this terminal
and the reference. It is also interesting to note that the positive and negative transfer
admittances contain an asymmetrical phase shift of 30 between the primary and secondary
terminals giving rise to nonreciprocal equivalent circuits. The asymmetrical phase shift is
entirely attributable to the star“delta connection and it is present even when no taps are
availabe in the transformer. It is common practice in application studies, such as positive
sequence power ¬‚ow and sequence domainbased fault levels to ignore the phase shift
during the calculations and then to account for it during the analysis of results.
82 MODELLING OF CONVENTIONAL POWER PLANT
3.4 ROTATING MACHINERY MODELLING
In general, synchronous machines are grouped into two main types, according to their rotor
structure: round rotor and salient pole machines (Grainger and Stevenson, 1994). Steam
turbine driven generators (turbogenerators) work at high speed and have cylindrical rotors.
The rotor carries a DC excited ¬eld winding. Hydro units work at low speed and have salient
pole rotors. They normally have damper windings in addition to the ¬eld winding. Damper
windings consist of bars placed in slots on the pole faces and connected together at both
ends. Turbogenerators contain no damper windings but the solid steel of the rotor offers a
path for eddy currents, which have similar damping effects.
For simulation purposes, the currents circulating in the solid steel or in the damping
windings can be treated as currents circulating in two closed circuits. Accordingly, a three
phase synchronous machine may be assumed to have three stator windings and three rotor
windings. This is illustrated in Figure 3.18, where all six windings are magnetically coupled.
ia
ea
a
q (t)
d
ifd
ikq
ikd
ib
eb
q
ec
ic
Figure 3.18 Schematic representation of a threephase synchronous generator. Redrawn by
permission of the Institution of Electrical Engineers from R.G. Wasley and M.A. Shlash, ˜Steady
state Phasevariable Model of the Synchronous Machine for Use in 3phase Load¬‚ow Studies™,
Proceedings of the IEEE 121(10) 1155“1164 # 1974 IEEE
The relative position of the rotor with respect to the stator is given by the angle between
the rotor™s direct axis and the stator™s phase a axis, termed the d axis and a axis, respectively.
In the stator, the axis of phases a, b, and c are displaced from each other by 120 electrical
degrees. In the rotor, the d axis is magnetically centred in the north pole of the machine. A
83
ROTATING MACHINERY MODELLING
second axis, located 90 electrical degrees behind the d axis is termed the quadrature axis or
q axis.
Three main control systems directly affect the turbinegenerator set, namely the boiler™s
¬ring control, the governor control, and the excitation system control. The excitation system
consists of the exciter and absolute value recti¬er (AVR). The latter regulates the generator
terminal voltage by controlling the amount of current supplied to the ¬eld winding by the
exciter. For the purpose of steadystate analysis, it is assumed that the three control systems
act in an idealised manner, enabling the synchronous generator to produce constant power
output, to run at synchronous speed, and to regulate voltage magnitude at the generator™s
terminal with no delay and up to its reactive power design limits.
3.4.1 Machine Voltage Equation
The objective of this section is to derive a steadystate expression for the stator threephase
voltages and currents of the synchronous generator (Wasley and Shlash, 1974b). The rotor
emfs (electromagnetic forces) and saliency are accounted for in the resulting voltage
equation, which may form the basis for connecting the machine model to a given three
phase bus of an unbalanced power system representation.
With reference to Figure 3.18, using stator and rotor quantities expressed in frames of
reference attached to their respective physical circuits, namely stator and rotor circuits, the
instantaneous voltages of the machine may be expressed as:
v ¼ Ri þ pLi; °3:105Þ
where R and L are the machine resistance and inductance matrices, respectively, and p is the
time derivative operator.
Furthermore, expanding Equation (3.105) into stator and rotor subsets, we obtain:
! !! !! ! !
Gss Gsr is Lss Lsr pis
vs Rs 0 is
¼ þ !r þ ; °3:106Þ
Gtsr 0 Ltsr Lrr pir
vr 0 Rr ir ir
where G ¼ dL=d; !r , equal to d=dt, is the rotor speed; and ¼ !r t þ .
The submatrix coef¬cients L, G, and R are:
2 3
Laa0 þ La2 cos°2Þ ÀLab0 À La2 cos°2 þ 60Þ ÀLab0 À La2 cos°2 À 60Þ
Lss ¼ 4 ÀLab0 À La2 cos°2 þ 60Þ Laa0 þ La2 cos°2 þ 120Þ ÀLab0 À La2 cos°2 À 180Þ 5;
ÀLab0 À La2 cos°2 À 60Þ ÀLab0 À La2 cos°2 À 180Þ Laa0 þ La2 cos°2 À 120Þ
°3:107Þ
2 3
Laf cos°Þ Laf cos°Þ ÀLaf sin°Þ
4 Laf cos° À 120Þ ÀLaf sin° À 120Þ 5;
Laf cos° À 120Þ
Lsr ¼ °3:108Þ
Laf cos° þ 120Þ Laf cos° þ 120Þ ÀLaf sin° þ 120Þ
2 3
Lfd Lmkd 0
4 Lmkd 0 5;
Lrr ¼ °3:109Þ
Lkd
0 0 Lkq
84 MODELLING OF CONVENTIONAL POWER PLANT
2 3
À2La2 sin°2Þ 2La2 sin°2 þ 60Þ 2La2 sin°2 À 60Þ
4 2La2 sin°2 þ 60Þ À2La2 sin°2 þ 120Þ 2La2 sin°2 À 180Þ 5;
Gss ¼ °3:110Þ
2La2 sin°2 À 60Þ 2La2 sin°2 À 180Þ À2La2 sin°2 À 120Þ
2 3
ÀLaf sin°Þ ÀLaf sin°Þ ÀLaf cos°Þ
Gsr ¼ 4 ÀLaf sin° À 120Þ ÀLaf sin° À 120Þ ÀLaf cos° À 120Þ 5; °3:111Þ
ÀLaf sin° þ 120Þ ÀLaf sin° þ 120Þ ÀLaf cos° þ 120Þ
2 3
Ra 0 0
40 0 5;
Rr ¼ °3:112Þ
Rb
0 0 Rc
2 3
Rfd 0 0
Rr ¼ 4 0 0 5: °3:113Þ
Rkd
0 0 Rkq
Since the rotor circuits are represented by a ¬eld winding on the d axis and two short
circuited damper windings on the d axis and q axis, respectively, the rotor voltage vector
may be written as
23
v fd
vr ¼ 4 0 5 ; °3:114Þ
0
where vfd is the applied direct ¬eld voltage.
For the purpose of steadystate analysis, it will be assumed that the applied direct ¬eld
voltage equals the voltage drop across the ¬eld resistance owing to the DC component of the
¬eld current and that additional voltages from Rf if can be neglected. Using such a
simpli¬cation, the relevant part of Equation (3.106) is solved for pir :
pir ¼ ÀLÀ1 ½!Grs is þ Lrs pis : °3:115Þ
rr
Assuming the following set of unbalanced stator currents:
2 3
I1 sin°!t þ 1 Þ
is ¼ 4 I2 sin°!t þ 2 Þ 5; °3:116Þ
I3 sin°!t þ 3 Þ
and the fact that the rotor runs at synchronous speed (i.e. !r ¼ !), we have,
2 3
Im cos° m Þ
X3
4 Im cos° m Þ 5;
pir ¼ À!Laf LÀ1 °3:117Þ
rr
m¼1 ÀIm sin° m Þ
where
9
¼ 2!t þ þ 1 ; >
=
1
¼ 2!t þ þ 2 À 120; °3:118Þ
2
>
;
¼ 2!t þ þ 3 þ 120:
3
85
ROTATING MACHINERY MODELLING
The expression for pir can be further simpli¬ed by substituting the inverse relation of
Equation (3.109) into Equation (3.117):
2 3
X k1 Im cos° m Þ
3
4 Àk2 Im cos° m Þ 5;
pir ¼ °3:119Þ
k2 Im sin° m Þ
m¼1
where
!9
!Laf Lmkd >
;>
k1 ¼ À 1À =
Lfd Lkd
°3:120Þ
>
!Laf >
;
k2 ¼ :
Lkd
It should be mentioned that the following practical simpli¬cations have been made while
substituting the inverse relation of Equation (3.109) into Equation (3.117): Lfd is much
greater than Lmkd, and Lkq ¼ Lkd.
Equation (3.119) is now integrated:
2 323
k1 Im sin° m Þ
1 X4
ifd
3
5 þ 4 0 5;
ir ¼ Àk2 Im sin° m Þ °3:121Þ
2! m¼1
Àk2 Im cos° m Þ 0
where ifd is the DC component of the ¬eld current.
Inspection of Equation (3.118) and (3.121) reveals that the presence of negative sequence
currents at the machine terminals gives rise to rotor currents of double the supply frequency.
In contrast, positive sequence currents are associated with zero frequency rotor currents,
other than the direct ¬eld current. Also, owing to balanced machine design considerations,
there is no contribution from zero sequence currents:
Substituting Equations (3.119) and (3.121) into Equation (3.106) we obtain a reduced
expression for the stator voltage vector:
2 3 2 32 3 2 32 3
va I1 sin°!t þ 1 Þ À1 À1 I1 cos°!t þ 1 Þ
Ra 0 0 1
676 76 7 6 76 7
4 vb 5 ¼ 4 0 0 54 I2 sin°!t þ 2 Þ 5 þ !Laa0 4 À1 1 À1 54 I2 cos°!t þ 2 Þ 5
Rb
vc I3 sin°!t þ 3 Þ À1 À1 I3 cos°!t þ 3 Þ
0 0 Rc 1
2 32 3
cos°!t þ 1 Þ cos°!t þ 2 À 120Þ cos°!t þ 3 þ 120Þ I1
°k1 À 2k2 ÞLaf 6 76 7
þ 4 cos°!t þ 1 þ 120Þ cos°!t þ 2 Þ cos°!t þ 3 À 120Þ 54 I2 5
4
cos°!t þ 1 À 120Þ cos°!t þ 2 þ 120Þ cos°!t þ 3 Þ I3
2 3
cos°!t À 1 þ 2Þ À cos°!t À 2 þ 2 þ 60Þ À cos°!t À 3 þ 2 À 60Þ
!La2 6 7
À 4 À cos°!t À 1 þ 2 þ 60Þ cos°!t À 2 þ 2 þ 120Þ À cos°!t À 3 þ 2 À 180Þ 5
2
À cos°!t À 1 þ 2 À 60Þ À cos°!t À 2 þ 2 À 180Þ cos°!t À 3 þ 2 À 120Þ
2 3 2 3
sin°!tÞ
I1
67 6 7
‚ 4 I2 5 À !Laf ifd 4 sin°!t þ À 120Þ 5:
°3:122Þ
sin°!t þ þ 120Þ
I3
86 MODELLING OF CONVENTIONAL POWER PLANT
The last term in Equation (3.122) may be interpreted as an array of rotor emfs. Moreover,
taking the stator a phase as reference,
232 3
va V1 sin°!tÞ
4 v b 5 ¼ 4 V2 sin°!t þ 2 Þ 5: °3:123Þ
vc V3 sin°!t þ 3 Þ
It is seen that a root mean square (rms) form of Equation (3.122) may be established very
readily. Also, by negating the stator currents to correspond to generator operating
conditions, we have
ES ¼ ½RS þ j°X1 þ X2 ÞIS þ jX3 IÃ þ VS ; °3:124Þ
S
where
2 3
ÀLab0 ÀLab0
Laa0
4 ÀLab0 ÀLab0 5;
X1 ¼ ! °3:125Þ
Laa0
ÀLab0 ÀLab0 Laa0
2 3
h2
Laf °k1 À 2k2 Þ 4 1 h
h2 5;
X2 ¼ °3:126Þ
h 1
4
h2 h 1
2 3
h2
1 h
!La2 e
j2
4 h2 1 5:
X3 ¼ À °3:127Þ
h
2
h2
h 1
It is observed that the term Laf °k1 À 2k2 Þ=4 reduces to À!L2 =4Lfd if damper windings are
af
not present.
As a means of evaluating the reactance elements in Equation (3.124), it is noted that the
usually available dq0 reactances of the machine may be used in the following expressions:
Á9
1À
!Laa0 ¼ Xd þ Xq þ X0 ; > >
>
>
>
3
Á=
À
1
°3:128Þ
!Lab0 ¼ Xd þ Xq À 2X0 ;
>
>
6 >
>
1À Á >
;
!La2 ¼ Xd À Xq :
3
Equation Á (3.124) includes the effect of machine saliency through matrix X3, where
À
Xd À Xq expresses the degree of saliency. Notice that if saliency can be ignored (i.e.
Xd ¼ Xq ) matrix X3 plays no role in machine performance. Also, X3 is dependent on
external circuit conditions through the machine angle . Matrix X2 contributes negative
sequence impedance, impairing the balanced behaviour of the machine.
3.5 SYSTEM LOAD
In general, power system loads can be classi¬ed into rotating and static loads (Weedy,
1987). A third category corresponds to power electronicbased loads. Rotating loads consist
mainly of induction and synchronous motors, and their steadystate operation is affected by
87
SYSTEM LOAD
voltage and frequency variations in the supply. Power electronicbased loads are also
affected by voltage and frequency variations in the supply. There is general agreement that
such loads are more dif¬cult to operate because, in addition to being susceptible to supply
variations, they inject harmonic current distortion back into the supply point (Acha and
Madrigal, 2001).
Detailed representation of a synchronous motor load in a threephase power ¬‚ow study
requires use of Equation (3.124), with changed signs to re¬‚ect the motoring action. An
expression of comparable detail can be derived for the induction motor load. However,
owing to the large number and diversity of loads that exist in power networks, it is
preferable to group loads and to treat them as bulk load points. It is only very important
loads that are singled out for detailed representation. It is interesting to note that a group of
rotating loads operating at constant torque may be adequately represented as a static load
that exhibits the characteristic of a constant current sink (Weedy, 1987).
In steadystate applications, most system loads are adequately represented by a three
phase power sink, which may be connected either in a star or delta con¬guration, depending
on requirements (Chen and Dillon, 1974). Figure 3.19(a) shows the schematic rep
resentation of a starconnected load with the star point solidly grounded, whereas
Figure 3.19(b) shows a schematic represantation of a deltaconnected load.
Va Va
Vb Vb
Vc Vc
SLc SLb SLa SLc SLb SLa
(a) (b)
Figure 3.19 System load representation: (a) starconnected load with star point solidly grounded and
(b) deltaconnected load
In threephase power ¬‚ow studies it is normal to represent bulk power load points as
complex powers per phase, on a perunit basis:
9
SLa ¼ PLa þ jQLa ; >
=
SLb ¼ PLb þ jQLb ; °3:129Þ
>
;
SLc ¼ PLc þ jQLc :
88 MODELLING OF CONVENTIONAL POWER PLANT
Re¬nements can be applied to the above equations to make the power characteristic more
responsive to voltage performance:
!9
!
1 1 >
;>
0
>
SLa ¼ PLa þjQLa
Va > >
>
Va
! >>
! =
1 1
0
°3:130Þ
SLb ¼ PLb þjQLb ;
Vb > >
Vb
! >>
! >
>
>
:>
1 1
0
;
SLc ¼ PLc þjQLc
Vc Vc
In Equations (3.130), and take values in the range 0“2 and Va, Vb, and Vc are the perunit
threephase nodal voltage magnitudes at the load point. Notice that when ¼ ¼ 0 the
complex power expressions in Equations (3.130) coincide with those in Equations (3.129).
However, if ¼ ¼ 1, Equations (3.130) resemble complex current characteristics more
than complex power characteristics. Also, if ¼ ¼ 2, the complex powers in Equa
tions (3.130) would behave like complex admittances.
The admittancelike characteristic in Equations (3.130) may be expressed in matrix form
for both kinds of load connections, star and delta, respectively:
2
2 3
SLa Va 0 0
2
40 0 5; °3:131Þ
SLb Vb 2
0 0 SLc Vc
2 2 2 2
2 3
SLa Va þ SLb Vb ÀSLb Vb ÀSLa Va
2 2 2 2
16 7
ÀSLb Vb SLb Vb þ SLc Vc ÀSLc Vc 5: °3:132Þ
4
2 2 2 2
3
ÀSLa Va ÀSLc Vc SLc Vc þ SLa Va
Moreover, if it is assumed that the load powers and voltage magnitudes are taken to be
balanced, SLa ¼ SLb ¼ SLc ¼ SL, and Va ¼ Vb ¼ Vc ¼ V, then application of the following
symmetrical component operation, Y012 ¼ TÀ1 Yabc TS , leads to the load model representa
S
tion for zero, positive, and negative (0, 1, 2) sequences:
2 2 3
SL V 0 0
40 05
2; °3:133Þ
SL V 2
0 0 SL V
2 3
0 0 0
2
40 05
2: °3:134Þ
SL V
0 0 SL V
Notice that no zero sequence loads exist for the case of a threephase deltaconnected load,
only positive and negative sequences.
As an extension of the above result, the positive, negative, and zero sequence expression
of a starconnected load with its star point solidly grounded may be expressed as
! !
1 1
0 0 0
SL °1Þ ¼ SL °2Þ ¼ SL °0Þ ¼ PL þjQL ; °3:135Þ
V V
89
SUMMARY
whereas for the case of a deltaconnected load we have
! !
1 1
0 0 0
SL °1Þ ¼ SL °2Þ ¼ PL þjQL ; SL °0Þ ¼ 0: °3:136Þ
and
V V
It should be remarked that the exponents and are not con¬ned to integer values and that
a wide range of load characteristics can be achieved by judicious selection of and ,
depending on the group of loads present in the study.
Also, a threephase delta connected load can always be transformed into an equivalent
star circuit by using a delta“star transformation. However, notice that the transform
ation will generate an extra bus in the form of the star point, which yields no physical
meaning.
3.6 SUMMARY
The chapter has addressed the mathematical modelling of the most common elements found
in conventional electrical power systems, namely, transmission lines, transformers,
generators, loads, and shunt and series passive compensation. The tools and methods
covered in the book are limited to fundamental frequency steadystate phenomena, and the
modelling approach followed in this chapter re¬‚ects this fact. Notwithstanding this, the
overall modelling philosophy is quite general in the sense that all plant component models
are formulated in the frame of reference of the phases, which is closely associated with the
physical structure of the equipment and its actual steadystate electrical operation. It is
shown throughout the chapter that simpler models do exist to represent a given plant
component but that these models are based on the assumption of perfect geometric balance
conditions. These models are realised with the help of the symmetrical component
transform, leading to plant component representation in the frame of reference of the
sequences.
Multiphase transmission line parameters are calculated with great accuracy, incorporating
all key effects that affect fundamental frequency operation such as geometric imbalances,
ground return loops, and even longline effects. Practical transmission lines include several
conductors per phase and ground wires as well as more than one threephase circuit sharing
the same right of way, giving rise to a large number of electromagnetically coupled
conductors. The chapter has presented a methodology for handling all these effects in a
systematic and ef¬cient manner. A comprehensive computer program in Matlab1 has been
written to calculate multiconductor transmission line parameters.
Threephase power transformers have been modelled in the frame of reference of the
phases, with particular reference to complex offnominal tapping positions. This caters for
the possibility of the transformer acting as a tap changer or as a phase shifter. The most
popular transformer connections used in highvoltage transmission have been addressed
and, under the assumption of perfect geometric conditions, transformer models in the frame
of reference of the sequences have been derived. The thrust of these models is fundamental
frequency, steadystate operation, and there is little loss of accuracy in representing the
threephase transformer as a threephase bank of transformers. A detailed model of the
synchronous generator, based on its physical windings arrangement, has been presented.
90 MODELLING OF CONVENTIONAL POWER PLANT
The effects of saliency and the generator load angle are explicitly represented in the model.
This model also serves the purpose of representing a synchronous motor, by suitable
modi¬cation of signs to conform to motoring action. Static loads suitable for bulk load
representation have also received attention.
The models of conventional power plant components developed in this chapter interface
quite naturally with the models of FACTS components developed in Chapter 2. Together,
they provide a very sophisticated tool with which to represent power system networks
containing a vast array of power electronic controllers of various kinds. These are the
power systems that may be in operation tomorrow. Two different modelling ¬‚avours
emerge from this modelling exercise, the frame of reference of the phases and the frame of
reference of the sequences, each one having its own time and space. Chapters 4, 5, and 7,
dealing with positive sequence power ¬‚ow and optimal power ¬‚ow, use the positive
sequence models derived in this and Chapter 2. Chapter 6 covers the topic of threephase
power ¬‚ow and uses the comprehensive models developed in the frame of reference of
the phases.
REFERENCES
Acha, E., Madrigal, M., 2001, Power System Harmonics: Computer Modelling and Analysis, John
Wiley & Sons, Chichester.
´
Acha, E., AmbrizPerez, H., FuerteEsquivel, C.R., 2000, ˜Advanced Transformer Control Modelling in
an Optimal Power Flow using Newton™s Method™, IEEE Trans. Power Systems 15(1) 290“298.
Anderson, P.M., 1973, Analysis of Power Systems, Iowa State University Press, Ames IA.
Arrillaga, J., Acha, E., Densem, T.J., Bodger, P.S., 1986, ˜Ineffectiveness of Transmission Line
Transpositions at Harmonic Frequencies™, Proceedings of the IEE: Part C 133(2) 99“104.
Arrillaga, J., Smith, B.C., Watson, N.R., Wood, A.R., 1997, Power System Harmonic Analysis,
John Wiley & Sons, Chichester.
Bowman, W.I., McNamee, J.M., 1964, ˜Development of Equivalent PI and T Matrix Circuits for
Long Untransposed Transmission Lines™, IEEE Trans. Power Apparatus and Systems PAS84(6)
625“632.
Carson, J.R., 1926, ˜Wave Propagation in Overhead Wires with Ground Return™, Bell System Technical
Journal, 5 (October) 539“554.
Chen, M.S., Dillon, W.E. 1974, ˜Power System Modelling™, Proceeding of the IEE 62(7) 901“915.
Chen, B.K., Chen, M.S., Shoults, R.R., Liang, C.C., 1990, ˜Hydrid Threephase Load Flow™, IEE
Proceedings on Generation, Transmission and Distribution: Part C 137(3) 177“185.
Deri, A., Tevan, G., Semlyen, A., Castanheira, A., 1981, ˜The Complex Ground Return Plane, Simpli¬ed
Model for Homogeneous and Multilayer Earth Return™, IEEE Trans. Power Apparatus and Systems,
PAS100(8) 3686“3693.
Elgerd, O.I., 1982, Electric Energy System Theory: An Introduction, McGrawHill, New York.
Grainger, J.J., Stevenson, W.D., 1994, Power System Analysis, McGrawHill, New York.
Harker, B.J., Arrillaga, J., 1979, ˜3Phase a.c./d.c. Load Flows™, Proceedings of the IEE 126(12) 1275“
1281.
Hesse, M.H., 1966, ˜Circulating Currents in Paralleled Untransposed Multicircuit Lines™, IEEE Trans.
Power Apparatus and Systems, PAS85(7) 802“820.
Laughton, M.A., 1968, ˜Analysis of Unbalanced Polyphase Networks by the Method of Phase Co
ordinates™, Proceedings of the IEE 115(8) 1163“1172.
91
REFERENCES
Shlash, M.A., 1974, ThreePhase Analysis of Unbalanced Power System Networks, PhD thesis,
Department of Electrical Engineering and Electronics, University of Manchester Institute of Science
and Technology, Manchester, UK.
Wasley, R.G., Shlash, M.A., 1974a, ˜Newton“Raphson Algorithm for 3Phase Load Flow™, Proceedings
of the IEE 121(7) 630“638.
Wasley, R.G., Shlash, M.A., 1974b, ˜Steadystate Phasevariable Model of the Synchronous Machine