00 WindingsAxes
00 WindingsAxes
00 WindingsAxes
1.0 Introduction
In these notes, we will describe the different windings on
a synchronous machine. We will confine our analysis to
two-pole machines of the salient pole rotor construction.
Results will be generalizable because
• A machine with p>2 poles will have the same
phenomena, except p times/cycle.
• Round rotor machines can be well approximated using a
salient pole model and proper designation of the
machine parameters.
We will also define an important coordinate frame that we
will use heavily in the future.
2.0 Defined axes
The magnetic circuit and all rotor winding circuits (which
we will describe shortly) are symmetrical with respect to
the polar and inter-polar (between-poles) axes. This
proves convenient, so we give these axes special names:
1
• Polar axis: Direct, or d-axis
• Interpolar axis: Quadrature, or q-axis.
The q-axis is 90° from the d-axis, but which way?
• Ahead?
• Or behind?
Correct modeling can be achieved either way, and some
books do it one way, and some another. We will remain
consistent with your text and choose the q-axis to lag the
d-axis by 90°, which is “consistent with the
recommendation and rationale of [15]1” (p. 92, VMAF).
1
IEEE Committee Report, Harrington, D. (chair), “Recommended phasor diagram for synchronous machines, IEEE
Trans. Power App. Syst. PAS-88: pp. 1593-1610, 1969.
2
Fig. 1 is from your text, and shows the q-axis lagging the d-
axis, consistent with our assumption.
Fig. 2 is from Kundur, and shows the q-axis leading the d-
axis, which we will NOT do.
Fig. 1 Fig. 2
Example 5.1 in Chapter 5 is worked with the q-axis lagging
the d-axis; Example 5.2 reworks the same problem with
the q-axis leading the d-axis. Example 5.2 states that ‘”The
objective of this example is to illustrate, though the
calculations differ slightly, that the results are identical;
thus, the choice of which coordinate reference frame is
assumed is arbitrary, as long as calculations are consistent
with the approach.”
3
3.0 Physical windings
There are typically 5 physical windings on a synchronous
machine:
• 3 stator windings (a-phase, b-phase, and c-phase)
• 1 main field winding
• Amortissuer windings on the pole-faces
The stator windings and the field winding are familiar to
you based on the previous notes. The amortissuer winding
might not be, so we briefly describe it.
Amortissuer is a French word for “shock absorber.” These
windings provide positive sequence damping under
transient conditions and so are good for stability2. Under
steady state, they are dead, i.e., they have no influence.
Amortissuer windings are not usually used on smooth-
rotor machines, because the solid steel rotor cores of such
machines provide paths for eddy currents and thus
produce the same effects as amortissuer windings3.
Amortissuer windings are often used in salient-pole
machines, but even when not, eddy currents in pole
2
E. Kimbark, “Power system stability, Vol. III: Synchronous Machines,” IEEE Press, 1995 (orig pub 1956), p. 219.
3
Ibid, p. 215.
4
faces contribute the same effect, although greatly
diminished.
Amortissuers have a number of other good effects, as
articulated by Kimbark in his Volume III book on
synchronous machines4 (he lists 10 such effects). For our
purposes, they affect machine subtransient data5.
Amortissuer windings are embedded in the pole-face (or
“shoe” of the pole) and consist of copper or brass rods
connected to end rings. They are similar in construction to
the squirrel cage of an induction motor.
Figures 3 (from Sarma) and 4 (from Kundur) illustrate
amortissuer windings. Note that they may be continuous
(Fig. 3a and Fig. 4) or noncontinuous (Fig. 3b).
Fig. 3
4
Ibid, p. 216.
5
Ibid, p. 217-219.
5
Fig. 4
4.0 Modeled windings and currents
Although there are typically 5 physical windings on a
machine, we will model a total of 7, with associated
currents as designated below.
• 3 stator windings: ia, ib, ic
• Field windings: There are 2: one physical; one fictitious
o Main field winding: carrying current iF and
producing flux along the d-axis.
o G-winding: carrying current iG and producing flux
along the q-axis. This is the fictitious one, but it
serves to improve the model accuracy of the round-
rotor machine (by modeling the q-axis flux
produced by the eddy-current effects in the rotor
6
during the transient period), and it can simply be
omitted when modeling the salient pole machine (in
salient pole machines, there is little q-axis flux
produced by the eddy current effect in the rotor).
The G-winding is like the F-winding of the main field,
except it has no source voltage in its circuit. Kimbark
suggests modeling it in his Vol. III, pg. 73.
• Amortissuer winding: This one represents a physical
winding for salient-pole machines with dampers, and a
fictitious winding otherwise. Because these produce flux
along both the d-axis and the q-axis, we model two
windings:
o d-axis: amortissuer winding carrying current iD
o q-axis: amortissuer winding carrying current iQ
It is important to understand the difference between the F
and G windings and the D and Q windings, respectively,
driven by the fact that D and Q windings have higher
resistance than F and G windings. Therefore:
• Both the F and D produce flux along the D-axis, but D is
“faster” (lower time constant or L/R ratio) than F.
• Both the G and Q produce flux along the Q-axis, but Q is
“faster” (lower time constant or L/R ratio) than G.
7
5.0 Flux linkages and currents
So we have seven windings (circuits) in our synchronous
machine. The flux linkage seen by any winding i will be a
function of
• Currents in all of the windings and
• Magnetic coupling between winding i and winding j, as
characterized by Lij, where j=1,…,7.
That is
7
i = Lij ij (1)
j =1
For example, the flux linking the main field winding is:
F = LFa ia + LFb ib + LFc ic + LFF iF + LFD iD + LFQ iQ + LFG iG (2)
8
a Laa Lab Lac LaF LaG LaD LaQ i
a
L Lbb Lbc LbF LbG LbD LbQ ib
stator b ba
c Lca
Lbc Lcc LcF LcG LcD LcQ ic
F = LFa LFb LFc LFF LFG LFD LFQ iF
G L
Ga
LGb LGc LGF LGG LGD LGQ iG (3)
rotor
D LDa i
LDb LDc LDF LDG LDD LDQ D
i
Q LQa LQb LQc LQF LQG LQD LQQ Q
9
iN i2
Lii = = (4a)
ii Ri
where Ri is the reluctance of the path seen by λi, given by
l
Ri = (4b)
A
where l is the mean length of the path, μ is the
permeability of the path’s material, and A is the cross-
sectional area of the path.
φF
N S
Fig. 5
10
Thus, since Lii=(Ni)2/R, rotor winding self-inductances are
constants, and we define the following nomenclature,
consistent with eq. (4.13) in VMAF.
• d-axis field winding LFF = LF (5)
• q-axis field winding LGG = LG (6)
• d-axis amortissuer winding: LDD = LD (7)
• q-axis amortissuer winding: LQQ = LQ (8)
Note your text’s convention of using only a single subscript
for constant terms.
6.1b Rotor-rotor terms: mutual inductances
Recall (eq. (15) in “Preliminary Fundamentals”) that:
i Ni N j
Lij = =
ij Rij (9)
11
Therefore, we have the following:
Fig. 6
12
6.2a Stator-stator terms: self inductances
We can derive these rigorously (see Kundur pp. 61-65) but
the insight gained in this effort may not be great. Rather,
we may be better served by gaining a conceptual
understanding of four ideas, as follows:
1. Sinusoidal dependence of permeance on θ: Due to
saliency of the poles (and to field winding slots in a smooth
rotor machine), the path reluctance seen by the stator
windings depends on θ, as illustrated in Fig. 7.
Rotation Rotation
φa φa
a' a'
a a
N S
13
a flux in this path is at a minimum, and permeance is at a
maximum.
From Fig. 7b, we observe that when θ=90°, the path of
phase-a flux contains more air that at any other angle
0→180°, and therefore the reluctance seen by the phase-
a flux in this path is at a maximum, and permeance is at a
minimum. This suggests a sinusoidal variation of
permeance with θ.
2. Constant permeance component: There will be a
constant permeance component due to the amount of
permeance seen by the phase-a flux at any angle. This will
include the iron in the middle part of the rotor (indicated
by a box in Figs. 7a and 7b), the stator iron, and the air gap.
Denote the corresponding component as Ps.
3. Double angle dependence: Because the effects
described in 1 and 2 above depend on permeance (or
reluctance), and not on rotor polarity, the maximum
permeance occurs twice each cycle, and not once.
Taking (1), (2), and (3) together, we may write that
P = Ps + Pm cos 2 (12)
14
4. Inductance: Because L=N2/R=N2P, the self inductance
of the a-phase winding can be written as
Laa = Ls + Lm cos 2 (13)
Likewise, we will obtain:
Lbb = Ls + Lm cos 2( − 120) (14)
Lcc = Ls + Lm cos 2( − 240) (15)
Equations (13), (14), (15) are denoted (4.12) in your text.
Note that because θ is a function of t, then (13), (14), (15)
imposes that stator self-inductances are functions of t!
(Recall our discussion in “Preliminary Fundamentals”.) This
means in Faraday’s law, e=d(Li)/dt, L is not constant.
6.2b Stator-stator terms: mutual inductances
We will identify 3 important concepts for understanding
mutual terms of stator-stator inductances.
1. Sign:
First, we need to remind ourselves of a preliminary fact:
➔ For any circuits i and j, Lij is positive if positive currents
in the two circuits produce fluxes in the same direction.
15
With this fact, we can state important concept 1:
➔As a result of defined stator current directions, the
stator-stator mutual inductance is always negative.
To see this, we can observe that the flux produced by
positive currents of a and b phases are in opposite
directions, as indicated in Fig. 8.
6 The thinking here is that if θ=0° results in max flux (min reluctance path) from a-current seen by a-winding, and if θ=-60° results
in max flux (min reluctance path) from b-current seen by b-winding, then halfway between the two will result in max flux (min
reluctance path) from a-(b-) current seen by b-(a-) winding.
17
Rotation Rotation
φa φa
b' b'
N
a' a'
a● a
●
S
φb S φb
●b ●b
19
• the rotor moves within the path of mutual flux and
causes the iron in the path to vary, and for this reason,
the path permeance varies.
Now, in this case, for stator-rotor terms (all mutuals),
• the rotor winding locations vary, the stator winding
locations are fixed, and so…
• the path of mutual flux varies, and so…
• the iron in the path of mutual flux varies, and for this
reason, the path permeance varies.
To illustrate, consider the permeance between the a-phase
winding and the main field winding (F).
• When the main field winding and the stator winding are
aligned, as in Fig. 10a, the permeance is maximum, and
therefore inductance is maximum.
20
φF
N
φa φa
a' φF a'
a a S
N
S
φa
a'
a
N
φF
Fig. 11
21
This discussion results in a conclusion that the mutual
inductance between a-phase winding and the main field
winding should have the form:
LaF = M F cos (20a)
The d-axis damper (amortissuer) winding is positioned
concentric with the main field winding, both producing flux
along the d-axis. Therefore, the reasoning about the
mutual inductance between the a-phase winding and the
d-axis damper winding will be similar to the reasoning
about the mutual inductance between the a-phase
winding and the main field (F) winding, leading to
LaD = M D cos (21a)
Now consider the mutuals between the a-phase winding
and the windings on the q-axis, i.e., the G-winding and the
Q damper (amortissuer) winding.
The only difference in reasoning about these mutuals and
the mutuals between the a-phase winding and the
windings on the d-axis (the F-winding and the D damper
winding) is that the windings on the q-axis are 90° behind
the windings on the d-axis. Therefore, whereas the a-
22
phase/d-axis mutuals were cosine functions, these
mutuals will be sine functions, i.e.,
LaQ = M Q sin (22a)
23
LaF = M F cos (20a)
LbF = M F cos( − 120) (20b)
LcF = M F cos( − 240) (20c)
LaD = M D cos (21a)
LbD = M D cos( − 120) (21b)
LcD = M D cos( − 240) (21c)
LaQ = M Q sin (22a)
LbQ = M Q sin( − 120) (22b)
LbQ = M Q sin( − 240) (22c)
24
7.0 Summary
Summarizing all of our needed expressions:
Rotor-rotor self terms: 5, 6, 7, 8
Rotor-rotor mutuals: 10a, 10b, 11a, 11b, 11c, 11d
Stator-stator self terms: 13, 14, 15
Stator-stator mutuals: 17, 18, 19
Rotor-stator mutuals: 20a, 20b, 20c, 21a, 21b, 21c, 22a,
22b, 22c, 23a, 23b, 23c
Counting the above expressions, we see that we have 28.
But let’s look back at our original flux linkage relation (3):
26