1

A Library of Second-Order Models
for Synchronous Machines

arXiv:1803.09707v1 [] 21 Mar 2018

Olaoluwapo Ajala, Student Member, IEEE, Alejandro Domı́nguez-Garcı́a, Member, IEEE,
Peter Sauer, Life Fellow, IEEE, and Daniel Liberzon, Fellow, IEEE

Abstract—This paper presents a library of two-dimensional
state-space models for salient-pole and round-rotor synchronous
machines. The dynamic states of the models are the power angle,
and the angular speed, and the fidelity of the models is verified
through a simulation-based comparison of the response of the
proposed models with that of a high-order model, and that of
the so-called classical model.
Index Terms—Synchronous machines, Reduced-order modeling, Singular perturbation analysis.

I. I NTRODUCTION

R

EDUCED-ORDER models of synchronous machines
find applications in power systems analysis, control
systems design and education, with each application requiring
simplified models that capture dynamical phenomena relevant
to the intended use. This has led to the proliferation of
synchronous machine models in the literature [1, 2, 3, 4], with
varying degrees of complexity, computational cost, and statespace dimension. One of such models is the so-called classical
model advocated in [5] and [6], a two-dimensional state-space
model which only captures dynamics of the machines output
voltage phase and angular speed. The main contribution of
this paper is the development of synchronous machine secondorder models—the elemental model, the damped model and
the semi-damped model—that are significantly more accurate
than the classical model and are therefore useful for a broader
range of applications.
In this work, we show that the classical model has certain
limitations that restrict its applications, and we propose a
library of high-fidelity second-order models for synchronous
machines with no such limitations. Using singular perturbation
analysis as our main tool [2, 7, 8, 9, 10], the second-order
models presented in this paper are derived from a high-order
model by identifying the fastest dynamic states in the highorder model, developing approximate manifold equations for
them, which are algebraic equations, and replacing the differential equations for these states with the algebraic counterparts.
Our approach to developing the proposed second-order
synchronous machine models is based on developments in
[2, 10, 11], where zero-order and first-order approximations of
manifolds for fast dynamic states are used to develop reducedorder models. In [10, 11], the use of integral manifolds for
model order reduction is introduced with some applications

The authors are with the Department of Electrical and Computer
Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801
USA. E-mail:{ooajala2, aledan, psauer, liberzon}@ILLINOIS.EDU.

presented, and in [2], the technique is used to develop the socalled two-axis model, a ninth-order model, the so-called oneaxis model, an eighth-order model, and the classical model, a
second-order model.
The remainder of the paper is organized as follows. In
Section II, we present the high-order model and the classical
model of a synchronous machine. In Section III, we develop
a library of high-fidelity second-order models from the highorder model, using singular perturbation analysis. Finally, in
Section IV we validate the second-order models developed,
using numerical examples, and in Section V we comment on
some implications of the presented results.
II. P RELIMINARIES
We begin this section by presenting a high-order model of
the synchronous machine adopted in this work. In addition, the
time-scale properties of the model are discussed. Afterwards,
we introduce the so-called classical model and describe how
it can be developed from the high-order model.
A. High-Order Synchronous Machine Model
The high-order synchronous machine model we describe
in this section is based on the developments in [2, 3]. The
components included in the model are: (i) three damper
windings, (ii) a wound-rotor synchronous machine, (iii) an
IEEE type DC1A excitation system [12], and (iv) a Woodward
diesel governor (DEGOV1) [13] coupled to a diesel engine,
which acts as the prime mover. Next, we provide mathematical expressions that describe the dynamic behavior of these
components. [Note that the model is presented utilizing the
qd0 transformation, with all parameters and variables scaled,
and normalized using the per-unit system].
Assumption 1. The synchronous machine is connected to an
electrical network bus through a short transmission line.
1) Damper windings model: Let Φq2 (t) and Ed0 (t) denote
the flux linkages of two damper windings aligned with the
quadrature axis (q-axis) of the synchronous machine, let
Φd1 (t) and Eq0 (t) denote the flux linkages of a damper
winding and a field winding, respectively, aligned with the
direct axis (d-axis) of the synchronous machine, and let Iq and
Id denote the q-axis and d-axis components of the stator output
current, respectively. Then, the damper winding dynamics can
be described as follows:

τq00 Φ̇q2 = − Φq2 − Xq0 − Xk Iq − Ed0 ,
(1)
τd00 Φ̇d1 = − Φd1 − (Xd0 − Xk ) Id + Eq0 (t),

2

τq0 Ėd0 = − Ed0 + Xq − Xq0



+(Xq0 − Xk )Iq − Ed0

Iq −


Xq0 − Xq00
Φq2
(Xq0 − Xk )2

,

(2)

where Xk denotes the machine leakage reactance, Xq denotes the machine stator reactance, Xq0 and Xd0 denote
machine transient reactances, and Xq00 denotes the machine

1
Xkq2 + Xmq ,
sub-transient reactance and τq00 = ω0 R
q
2


Xkq1 +Xmq
Xmd Xkf
1
0 =
X
+
and
τ
are
τd00 = ω0 R
k
q
1
X
+X
ω0 Rq1
d1
md
kf
time constants, with Xkq2 , Xkd1 , Xkf , Xkq1 denoting leakage
reactances, Xmq and Xmd denoting mutual reactances, and
Rq2 , Rd1 , Rq1 denoting winding resistances.
(s)

2) Stator windings and network model: Let Φq (t) and
(s)
Φd (t) denote the q-axis and d-axis components of flux
(e)
linkages for the stator windings, respectively, let Φq (t) =
(e)
−X (e) Iq , and Φd (t) = −X (e) Id denote the q-axis and
d-axis components of flux linkages for the electrical line,
respectively, let ω (s) (t) denote the machine angular speed, in
electrical radians per second, and let δ (s) (t) denote the power
angle of the synchronous machine in electrical radians. At
the electrical network bus, let V (l) and δ (l) denote the voltage
magnitude, in per unit, and the voltage phase relative to a reference frame rotating at the nominal frequency,
in electrical
radi

(l)
(l)
ans, respectively. Let Vq := V (l) cos δ (s) − δ (l) , Vd :=


(e)
(s)
V (l) sin δ (s) − δ (l) , Φq (t) := Φq (t) + Φq (t), Φd (t) :=
(s)

(e)

Φd (t) + Φd (t). Then, the stator winding and network
dynamics are described by:
δ̇ (s) = ω (s) (t) − ω0 ,


ω (s) (t)
1
Φ̇q = −
Φd + Vq(l) + Rs + R(e) Iq ,
ω0
ω0


1
ω (s) (t)
(l)
Φ̇d =
Φq + Vd + Rs + R(e) Id ,
ω0
ω0
Xq00 − Xk
Xq0 − Xq00
(e)
Φq2 −
Ed0 ,
Φq = − Xq00 Iq +
Xq0 − Xk
Xq0 − Xk
Xd0 − Xd00
Xd00 − Xk
(e)
Φd = − Xd00 Id +
Φd1 +
Eq0 (t),
0
Xd − Xk
Xd0 − Xk
1 (e)
ω (s) (t) (e)
Φ̇q = R(e) Iq −
Φd − Vq(s) + Vq(l) ,
ω0
ω0
1 (e)
ω (s) (t) (e)
(s)
(l)
Φ̇d = R(e) Id +
Φq − V d + V d ,
ω0
ω0
(3)
(e)
(e)
where Xq00 := Xq00 + X (e) , and Xd00 := Xd00 + X (e) ,
X (e) denotes the per-phase line reactance, Xd00 denotes a
machine sub-transient reactance, R(e) denotes the per-phase
line resistance, Rs denotes the per-phase stator resistance, and
ω0 denotes the nominal frequency in electrical radians per
second.
3) Excitation system model: Let Ef (t) denote the output
voltage of the machines excitation system, let Uf (t) denote
the exciter control input, let Ūf (t) denote the rate feedback
variable of the voltage regulator (see [2], pp. 71–72 for

details), and let V (s) denoterthe magnitude of the machine

2 
2
(s)
(s)
Vq
+ Vd
.
output voltage, with V (s) :=
Assumption 2. The effects of magnetic saturation on the
machines excitation system are negligible.
Then, the dynamics of the machines excitation system can be
described as follows:

Xd0 − Xd00
(Φd1
τd0 Ėq0 = − Eq0 − (Xd − Xd0 ) Id −
(Xd0 − Xk )2


+(Xd0 − Xk )Id − Eq0 + Ef ,
τf Ėf = − Kf Ef + Uf ,
τu U̇f = − Uf + Ku Ūf −



Ku K̄u
Ef + Ku Vr(s) − V (s) ,
τ̄u

K̄u
Ef ,
τ̄u Ū˙ f = − Ūf +
τ̄u
(4)
(s)
Vr

denotes the reference voltage magnitude, τd0 =
where
Xf
R̄f
Lf
Lt +Lm
t2 Lm
, K̄u = N
ω0 Rf , τf = Kg , Kf = Kg , τ̄u =
Rt
Nt2 Rt , Xd
denotes the machine stator reactance, τu denotes the amplifier
time constant, Ku denotes the amplifier gain, Xf denotes
the field winding reactance, Rf denotes the field winding
resistance, Lf denotes the unsaturated field inductance, Kg
denotes the slope of the unsaturated portion of the exciter
saturation curve, R̄f denotes the exciter circuit resistance, Lt
and Lm denote series and magnetizing inductances of the
stabilizing transformer, which is used to stabilize the excitation
system through voltage feedback [2], respectively, Rt denotes
Nt2
the series resistance of a stabilizing transformer, and N
t1
denotes the turns ratio of the stabilizing transformer.
4) Prime mover and speed governor model: Let Tm (t)
denote the mechanical torque output of the machine, and let
Pu (t) denote the valve position of the diesel engine, which acts
as the prime mover. For the speed governor system, let Pa2
denote the output of its actuator, with Ṗa1 = Pa2 , and let Pb2
denote the output of its electric control box, with Ṗb1 = Pb2 .
Then, the speed control system of the synchronous machine
can be expressed as follows:
M ω̇ (s) = Tm − Φd (t)Iq + Φq (t)Id − D̃0 ω (s) ,
τm Ṫm = − Tm + Pu ,
Ṗu = Pa1 + τ4 Pa2 ,

1
Pa1 − κ (Pb1 + τ3 Pb2 ) − Pa2 ,
τa2 Ṗa2 = −
τ5 +
 τ6

1  (s)
1
1
τ2 Ṗb2 = −
(Pc − Pu ) +
Pb1 −
ω − ω0
τ1
ω0
D̄0 ω0
− Pb 2 ,
(5)
where τ2 , τ3 , τ4 , τ5 and τ6 denote time constants of the
τ6
control system, τa2 = ττ55+τ
, κ denotes a controller gain
6
for the actuator, Pc denotes the power change setting of the
machine, M denotes the inertia of the machine, D̃0 denotes
the friction and windage damping coefficient of the machine,
τm denotes the time constant of the engine, and D̄0 = RD1ω0 ,
with RD denoting the droop coefficient.

3

Note that for salient pole machines, Xq = Xq0 , so that
Ed0 (t) = 0, and for round-rotor machines, Xq = Xd .
B. High-Order Model Time-Scale Properties
The following observations are based on standard parameter
values obtained from synchronous machine models in [1, 2,
3, 13], and an eigenvalue analysis of these models.
(e)

Observation 1. The dynamics of Φq2 , Φd1 , Ed0 , Φq , Φd , Φq ,
(e)
Φd , Eq0 , Ef , Uf , Ūf , Tm , Pu , Pa2 , Pb2 , Pa1 and Pb1 , are
much faster than those of ω (s) and δ (s) .
Observation 2. For  = 0.1 denoting a constant, the parameters Rs , τq00 , τq0 , ω10 , τf , τu , τ̄u , τm , τa2 , τ2 , τ1 , (τ5 + τ6 ),
τ5 τ6
1
1
(τ5 +τ6 ) , κRD are O () .
Based on these observations, the fifteenth order machine model
described by (1)–(5) can be expressed compactly as:

ẋ(t) = f x(t), z(t),  ,
x(0) = x0 ,

(6)
ż(t) = g x(t), z(t),  ,
z(0) = z 0 ,
h
i>

where x(t) = δ (s) ω (s) , and z(t) = Φq Φd Ed0
(e)

(e)

Φq2 Φd1 Φq
Φd
E 0 Ef Uf Ūf Tm Pu
> q
Pa2 Pb2 Pa1 Pb1 . In the remainder of this paper, we
refer to the elements of z(t) as the fast states, and the elements
of x(t) as the slow states. Other observations, which will prove
useful in Sections III-B and III-C are:
(e)

(e)

Observation 3. The dynamics of Φq , Φd , Φq and Φd
much faster than those of Φq2 , Φd1 , Ed0 and Eq0 .

III. L IBRARY OF S ECOND -O RDER M ODELS
In this section, a library of two-dimensional state-space
models for synchronous machines are developed from the
high-order model presented in Section II-A. By employing the
time-scale properties described in Section II-B, and singular
perturbation analysis, the fifteenth order machine model is
reduced to the elemental model, the damped model, and the
semi-damped model. The models are based on the following
simplifying assumption:
Assumption 3. The angular speed of the machine, ω (s) (t), is
sufficiently close to the nominal speed of the machine so that
ω (s) (t)
= 1 + O ().
ω0
A. The Elemental Model
The elemental model is formulated by replacing the differential equations for the fast states with algebraic counterparts
called zero-order approximate manifolds. The manifolds are
developed by setting τd00 , τd0 , and all O () parameters in (6)
to zero, to give:

Φq2 ,0 (t) = − Xq − Xk Iq ,

Ed0 ,0 (t) = Xq − Xq0 Iq ,


Φq,0 (t) = − Rs(e) Id − V (l) sin δ (s) (t) − δ (l) ,


Φd,0 (t) = Rs(e) Iq + V (l) cos δ (s) (t) − δ (l) ,
(e)

Φq,0 (t) = − X (e) Iq ,

are

(e)

Φd,0 (t) = − X (e) Id ,


(s)
Ku Vr − V (s)
,
Ef,0 (t) =
Kf
Eq0 ,0 (t) = − (Xd − Xd0 ) Id + Ef,0 (t),

Observation 4. The dynamics of Φq2 and Φd1 are much faster
than those of Ed0 and Eq0 .

Φd1 ,0 (t) = − (Xd − Xk ) Id + Ef,0 (t),

C. Classical Model
The classical model of a synchronous machine is a secondorder model whose formulation is based on the following
assumptions [14]: (i) the machine can be modeled as a constant
magnitude voltage source with a series reactance, (ii) the
mechanical rotor angle of the machine can be represented by
the angle of the voltage source, (iii) damping can be neglected,
and (iv) the machines mechanical power input is constant.
Thus, the classical model can be obtained from the high-order
(s)
model by setting τq00 = 0, τd00 = 0, ω10 = 0, ω ω0(t) = 1,
Rs = 0, Re = 0, Xq0 = Xd0 , τq0 = ∞, τd0 = ∞, τm = ∞ to
give:
δ̇ (s) = ω (s) − ω0 ,
M ω̇ (s) = Tm (0) −

(8)

E0



(l)
(s)
(l)
V
sin
δ
−
δ
− D̃0 ω (s) ,
(e)

Uf,0 (t) = Kf Ef,0 (t),
K̄u
Ūf,0 (t) =
Ef,0 (t),
τ̄u



Pu,0 (t) = Pc − D̄0 ω (s) (t) − ω0 ,
Tm,0 (t) = Pu,0 (t),
also, Pa1 ,0 = 0, Pa2 ,0 = 0, Pb1 ,0 = 0, and Pb2 ,0 = 0,
where ‘0’ subscripts denote a zero-order approximation, and
(e)
(e)
(e)
Rs := Rs + R(e) , Xq := Xq + X (e) , Xd := Xd + X (e) .
(s)
Also, the output 
voltage is described
by: Vq = R(e) Iq +

(s)
X (e) Id +V (l) cos δ (s) (t) − δ (l) , Vd = R(e) Id −X (e) Iq +


V (l) sin δ (s) (t) − δ (l) , and the output current is described


Xd0

(7)
q

2
2
(e)
where E0 =
Eq0 (0) + Eq0 (0) and Xd0 := Xd0 +
X (e) denote constants.
1 Consider a positive constant , where  < 1, and a function f (), defined

on some subset of the real numbers. We write f () = O i if and only if
i
there exists a positive real number k, such that: f () ≤ k , as  → 0.

Rs(e) 

by Iq =





(s)
Ku Vr −V (s)


(e) 2

Rs



Kf
(e)

(e)

+Xq Xd

−




Rs(e) V (l) cos δ (s) (t)−δ (l)


(e) 2
(e)
(e)
Rs
+Xq Xd



(s)

+

Ku Vr −V (s)




Xq(e) 
Kf
V (l) sin δ (s) (t)−δ (l)




, Id =
−
(e) 2
(e)
(e)
(e) 2
(e)
(e)
Rs
+Xq Xd
Rs
+Xq Xd






Xq(e) V (l) cos δ (s) (t)−δ (l)
Rs(e) V (l) sin δ (s) (t)−δ (l)

2


−
. Substi(e)
(e)
(e)
(e) 2
(e)
(e)
Rs
+Xq Xd
Rs
+Xq Xd
(e)

Xd

4

tuting the zero-order approximate manifolds formulated in this
section into (1)–(5), the elemental model is given by:
δ̇ (s) = ω (s) − ω0 ,

2
2
M ω̇ (s) = Pr(s) − D0 ω (s) − Rs(e) (I) + Cr V (l)




− Ck Cr Vr(s) − V (s) V (l) cos δ (s) − δ (l) (9)




− Ck C̃x Vr(s) − V (s) V (l) sin δ (s) − δ (l)


Cx  (l) 2
−
sin 2 δ (s) − δ (l) ,
V
2
where Cr , Ck , C̃x and Cx are constants, with Cr =
X (e)
R(e)
u
 s
 q


, Ck = K
, Cx =
(e)
(e)
(e)
(e)
(e) 2
(e) 2
Kf , C̃x =
+Xq Xd
+Xq Xd
Rs
Rs
(s)
Xd −Xq
(s)
,
and
P
=
P
+
D̄
ω
,
D
=
D̄
=
r
c
0 0
0
0 + D̃0 , V
(e)
(e)
Xq Xd
r
q

2 
2
2
(s)
(s)
2
Iq + (Id ) . The dyVq
+ Vd
, and I =

namic circuit of the elemental model is depicted in Fig. 1. For
(e)
the special case where Rs and R(e) are O (), we set Rs = 0,
1
from where it follows that Cr = 0, and C̃x = (e) .
Xd

I

q

 jI d  e j 

(s)

(t )

jX d

R(e)

jX ( e )

Rs
+


 j X

Ck Vr

(s)

d

V

(s)

e

j ( s ) ( t )

 X q  I q e j

+
(s)

V

(s)
q

 jVd( s )

+

e

j

(s)

(t )

V ( l ) e j

(l )

(t )

(t )
-

-

Fig. 1: Dynamic circuit of synchronous machine elemental
model.

B. The Damped Model
The damped model is formulated by replacing (1) and
(2) with first-order approximate manifolds, and replacing the
differential equations for other fast states with zero-order
approximate manifolds. By using a first-order approximation
for the damper windings manifolds, the effects of damper
windings on the machine response are captured by the resulting reduced model. The following simplifying assumption
is employed:
Assumption 4. The per-phase line resistance, R(e) is O ().
Starting with the states observed to have the fastest dynamics,
(e)
(e)
Φq (t), Φd (t), Φq (t) and Φd (t), we formulate the following
zero-order approximations by setting Rs = 0, R(e) = 0, ω10 =
0 and

ω (s) (t)
ω0

= 1:



Φq,0 (t) = − V (l) sin δ (s) (t) − δ (l) ,


Φd,0 (t) = V (l) cos δ (s) (t) − δ (l) ,


(e)
(s)
Φq,0 (t) = Vd − V (l) sin δ (s) (t) − δ (l) ,


(e)
Φd,0 (t) = − Vq(s) + V (l) cos δ (s) (t) − δ (l) ,

(10)

(Xq0 −Xq00 ) 
Φq2 (t)
X
( q0 −Xk ) Xq(e)
00

V (l) sin δ (s) (t)−δ (l)

from where it follows that: Iq =

−

(Xq00 −X
k)

 Ed0 (t)
,
and
+
(e)
Xq00
(Xq0 −Xk ) Xq(e)
00
(Xd0 −Xd00 ) 
(Xd00 −X
k)


0
Id =
(e) Φd1 (t) +
(e) Eq (t) −
Xd0 −Xk ) Xd00
Xd0 −Xk ) Xd00
(
(


V (l) cos δ (s) (t)−δ (l)
(e)

Xd00

. Next, for the subsequent fastest states,

Φq2 (t) and Φd1 (t), which are damper winding states, we
derive a first-order approximation of its manifold. Manifolds
for Φq2 (t) and Φd1 (t), can be expressed as power series in
τq00 and τd00 , respectively, to give:
Φq2 (t) = Φq2 ,0 (t) + τq00 Φq2 ,1 (t) + (τq00 )2 Φq2 ,2 (t) + · · · ,
Φd1 (t) = Φd1 ,0 (t) + τd00 Φd1 ,1 (t) + (τd00 )2 Φd1 ,2 (t) + · · · ,
(11)
from where it follows that first-order approximations are given
by:
Φq2 (t) ≈ Φq2 ,0 (t) + τq00 Φq2 ,1 (t),
Φd1 (t) ≈ Φd1 ,0 (t) + τd00 Φd1 ,1 (t).

(12)

Expressions for Φq2 ,0 (t), Φq2 ,1 (t), Φd1 ,0 (t) and Φd1 ,1 (t) are
derived using the following steps:
• Substitute (11) into (1) to give:

d
Φq2 ,0 (t) + τq00 Φq2 ,1 (t) + · · · = − Φq2 ,0 (t)
τq00
dt


+τq00 Φq2 ,1 (t) + · · · − Xq0 − Xk Iq − Ed0 (t),

d
τd00
Φd1 ,0 (t) + τd00 Φd1 ,1 (t) + · · · = − Φd1 ,0 (t)
dt

+τd00 Φd1 ,1 (t) + · · · − (Xd0 − Xk )Id + Eq0 (t).
(13)
• Using the zero-order approximations in (10), substitute
0
expressions for Iq and Id into (13) and equate the τq00
0
and (τd00 ) terms to give:
(e)


X
Xq0 − Xk (l)
V sin δ (s) (t) − δ (l) ,
Φq2 ,0 (t) = − k(e) Ed0 (t) −
(e)
Xq0
Xq0
(e)


X
Xd0 − Xk (l)
(s)
(l)
Φd1 ,0 (t) = k(e) Eq0 (t) +
V
cos
δ
(t)
−
δ
,
(e)
Xd0
Xd0
(e)

where Xk := Xk + X (e) .
1
1
• Also equate the τq00 and (τd00 ) terms to give:



(e)
(e)

Xq00
Xk

Φq2 ,1 (t) = −
Xq(e) Ed0 (t) − Xq − Xq0 V (l)

3
(e)
τq0 Xq0



(e)


Xq00
Xq0 − Xk
(l)
· sin δ (s) (t) − δ (l)
+
V̇d ,

2
(e)
Xq0



(e)
(e)

Xd00
Xk
(e)
Φd1 ,1 (t) =
Xd Eq0 (t) − (Xd − Xd0 ) V (l)

3
(e)
τd0 Xd0


(e)


Xd00 (Xd0 − Xk )
· cos δ (s) (t) − δ (l)
−
V̇q(l)

2
(e)
Xd0

5





(e)
(e)
Xd00
Xk
−

2 Ef (t),
(e)
τd0 Xd0



(l)
where V̇d
= V (l) cos δ (s) (t) − δ (l) δ̇ (s) (t) − δ̇ (l) +




(l)
V̇ (l) sin δ (s) (t) − δ (l) , V̇q = V̇ (l) cos δ (s) (t) − δ (l) −



(e)
V (l) sin δ (s) (t) − δ (l) δ̇ (s) (t) − δ̇ (l) , and Xq0 := Xq0 +

following zero-order approximations by setting τd0 , and all
O () parameters except τq00 and τq0 , to zero:


(s)
Ku Vr − V (s)
Ef,0 (t) =
,
Kf
(e)

X (e) . Next, for the damper winding state observed to have the
slower dynamics, Ed0 (t), we derive a first-order approximation
of its manifold. A manifold for Ed0 (t) can be expressed as a
power series in τq0 to give:
Ed0 (t) = Ed0 ,0 (t) + τq0 Ed0 ,1 (t) + (τq0 )2 Ed0 ,2 (t) + · · · , (14)
from where it follows that a first-order approximation is given
by:
Ed0 (t) ≈ Ed0 ,0 (t) + τq0 Ed0 ,1 (t).

(15)

Expressions for Ed0 ,0 (t) and Ed0 ,1 (t) can be derived using the
following steps:
• Substitute (12) and (14) into (2) to give:

d
Ed0 ,0 (t) + τq0 Ed0 ,1 (t) + · · · =
τq0
dt

− Ed0 ,0 (t) + τq0 Ed0 ,1 (t) + · · ·




Xq0 Xq0 − Xq00
+ Xq − Xq0 Iq −
2 τq00 Φq2 ,1 (t) .
Xq00 Xq0 − Xk
(16)
• Using the zero-order approximations in (10), substitute the
expressions for Iq and Φq2 ,1 (t), into (16), and equate the
0
τq0 terms to give:
Ed0 ,0 (t) =

Xq − Xq0
(e)

Xq


 N
q (l)
V̇ ,
V (l) sin δ (s) (t) − δ (l) −
Dq d

(e) (e)
τq0 τq00 Xq0 Xk



where Nq =
Xq − Xq0 Xq0 − Xq00

2

2
(e)
(e)
· Xq0 − Xk , Dq = τq0 Xq
Xq0
Xq0 − Xk −

2


(e)
(e)
τq00 Xq
Xk
Xq − Xq0 Xq0 − Xq00 .
1
• Also equate the τq0 terms to give:

Nq0 (l)
Ed0 ,1 (t) = −
V̇d + O τq0 ,
D̃q

3

2
(e)
where Nq0 = τq0 Xq0
Xq − Xq0 Xq0 − Xk and
(e)

D̃q = Xq Dq
Finally, for other states observed to have fast dynamics, Eq0 ,
Ef , Uf , Ūf , Tm , Pu , Pa2 , Pb2 , Pa1 , Pb1 , we derive the

Xd0

Nd (l)
V̇
Dd q


Xd − Xd0 (l)
(s)
(l)
,
V
cos
δ
(t)
−
δ
+
(e)
Xd
Uf,0 (t) = Kf Ef,0 (t),
K̄u
Ūf,0 (t) =
Ef,0 (t),
τ̄u



Eq0 ,0 (t) =

(e)
Xd

Ef,0 (t) −

(17)

Pu,0 (t) = Pc − D̄0 ω (s) (t) − ω0 ,
Tm,0 (t) = Pu,0 (t),

also, Pa1 ,0 = Pa2 ,0 = Pb1 ,0 = Pb2 ,0 = 0, where Nd =
(e) (e)
τd0 τd00 Xd0 Xk (Xd − Xd0 ) (Xd0 − Xd00 ) (Xd0 − Xk ) , and

2

2
(e)
(e)
(e)
(e)
2
Dd = τd0 Xd Xd0
(Xd0 − Xk ) − τd00 Xd Xk
· (Xd − Xd0 ) (Xd0 − Xd00 ) .
Substituting the first-order and zero-order approximate manifolds in (12), (15), (10) and (17) into (1)–(5), and setting
O (τq0 )2 terms to zero, the damped model for a non-salient
pole machine is given by:
δ̇ (s) = ω (s) − ω0 ,


Cx  (l) 2
V
sin 2 δ (s) − δ (l)
M ω̇ (s) = Pr(s) − D0 ω (s) −
2



Ck  (s)
(s)
− (e) Vr − V
V (l) sin δ (s) − δ (l)
Xd

2



− Cq V (l) cos2 δ (s) − δ (l) δ̇ (s) − δ̇ (l)

2



− Cd V (l) sin2 δ (s) − δ (l) δ̇ (s) − δ̇ (l)



Cq − Cd (l) (l)
V̇ V sin 2 δ (s) − δ (l) ,
−
2
(18)
Ck ,

where
Ck =

Ku
Kf ,

Cx ,
Cq ,
and
Cd are constants,

2
(Xd −Xq )
Cx = (e) (e) , Cq = Cq00 + Cq0 + C̃q00 C̃q ,

with Cq00 =
(e)

τq0 Xq0

Xq Xd
τq00 (Xq0 −Xq00 )


,
(e) 2
Xq 0

Cd

Cd00

=

(Xq −Xq0 )
D̃q

C̃q00

Xd0

(e)

, Cq 0 =

(e)
τq00 Xq(e) Xk (Xq0 −Xq00 )
,
=
(e)

 Xq 0
=
Cd00 +
Cd0 + C̃d00 C̃d00 C̃d , with
τd00 (Xd0 −Xd00 )
(Xd −Xd0 )


, C̃d =
, Cd0 =
(e) 2
D̃


Xq0 − Xk ,

and

C̃q =

τd0 Xd0 (Xd0 − Xk ) ,

d

(e)

C̃d00

=

(e)

τd00 Xd Xk

(Xd0 −Xd00 )

(e)

Xd0

.The

dynamic circuit of the damped
model is depicted in Fig.

(l)
2, with V̇d = V (l) cos δ (s) (t) − δ (l) δ̇ (s) (t) − δ̇ (l) +




(l)
V̇ (l) sin δ (s) (t) − δ (l) , V̇q = V̇ (l) cos δ (s) (t) − δ (l) −



V (l) sin δ (s) (t) − δ (l) δ̇ (s) (t) − δ̇ (l) . Note that for salient
pole machines, C̃q = 0, whereas for round-rotor machines,
Cx = 0.

6

I

C V
 j X
k

(s)

r

d

V

(s)

 X

(e)
d

q

(l )
q

Cd V

 X q  I q  X q( e )CqVd( l )

 jI d  e j 

e
e

(s)

(t )

jX ( e )

jX d

R(e)

+

j ( s ) ( t )
+
-



(s)
q

V

+

 jVd( s )

 e j

(s)

( l ) j ( l ) ( t )

(t )

V e

j ( s ) ( t )

-

-

Fig. 2: Dynamic circuit of synchronous machine damped
model.

C. The Semi-Damped Model
This model is only applicable to round-rotor machines. Employing Assumption 4, the semi-damped model is developed
by replacing (2) with a first-order approximate manifold, and
replacing the differential equations for other fast states with
zero-order approximate manifolds. Starting with fastest states
(e)
(e)
Φq (t), Φd (t), Φq (t), Φd (t), Φq2 (t) and Φd1 (t), we develop
the following zero-order approximations by setting Rs = 0,
(s)
R(e) = 0, ω10 = 0, ω ω0(t) = 1, τq00 = 0 and τd00 = 0:


Φq,0 (t) = − V (l) sin δ (s) (t) − δ (l) ,


Φd,0 (t) = V (l) cos δ (s) (t) − δ (l) ,


(e)
(s)
Φq,0 (t) = Vd − V (l) sin δ (s) (t) − δ (l) ,
(19)


(e)
Φd,0 (t) = − Vq(s) + V (l) cos δ (s) (t) − δ (l) ,

Φq2 ,0 (t) = − Xq0 − Xk Iq − Ed0 (t),
Φd1 ,0 (t) = − (Xd0 − Xk ) Id + Eq0 (t),
from where it follows that: Iq


V (l) sin δ (s) (t)−δ (l)
(e)

Xq0

, Id =

1
(e) Eq 0 (t)
Xd0

1
(e) Ed0 (t) +
Xq0


V (l) cos δ (s) (t)−δ (l)

=
−

−

(e)

Xd0

.

Next, we derive a first-order approximate manifold for Ed0 (t)
having the form Ed0 (t) ≈ Ed0 ,0 (t) + τq0 Ed0 ,1 (t). Substituting
the expression for Φq2 ,0 (t) in (19), and the power series
expansion in (14) into (2), it follows that:

d
τq 0
Ed0 ,0 (t) + τq0 Ed0 ,1 (t) + · · · =
dt
 Xq(e)
− Ed0 ,0 (t) + τq0 Ed0 ,1 (t) + · · ·
(20)
(e)
Xq0


Xq − Xq0 (l)
(s)
(l)
+
V
sin
δ
(t)
−
δ
.
(e)
Xq0
0
Equating the τq0 terms in (20), we have that:


Xq − Xq0 (l)
(s)
(l)
Ed0 ,0 (t) =
V
sin
δ
(t)
−
δ
,
(21)
(e)
Xq
1
and equating the τq0 terms, we have that:

(e)
Xq0 Xq − Xq0 (l)
Ed0 ,1 (t) = −
V̇d ,
(22)

2
(e)
Xq



(l)
where V̇d
= V (l) cos δ (s) (t) − δ (l) δ̇ (s) (t) − δ̇ (l) +


V̇ (l) sin δ (s) (t) − δ (l) . Finally, for other states observed to

have fast dynamics, i.e., Eq0 , Ef , Uf , Ūf , Tm , Pu , Pa2 , Pb2 ,
Pa1 and Pb1 , the following zero-order approximations, derived
by setting τd00 , τd0 , and all O () parameters except τq0 , to zero
are used:


(s)
Ku Vr − V (s)
,
Ef,0 (t) =
Kf
(e)


X 0
Xd − Xd0 (l)
(s)
(l)
,
Eq0 ,0 (t) = d(e) Ef,0 (t) +
V
cos
δ
(t)
−
δ
(e)
Xd
Xd
Uf,0 (t) = Kf Ef,0 (t),
K̄u
Ūf,0 (t) =
Ef,0 (t),
τ̄u


Pu,0 (t) = Pc − D̄0 ω (s) (t) − ω0 ,

Tm,0 (t) = Pu,0 (t),
(23)
with, Pa1 ,0 = Pa2 ,0 = Pb1 ,0 = Pb2 ,0 = 0. Substituting the
zero-order approximate manifolds in (23), and the first-order
approximate manifold Ed0 ,0 (t) + τq0 Ed0 ,1 (t), into (1)–(5), the
semi-damped model is given by:
δ̇ (s) = ω (s) − ω0 ,


C̃q0 (l) (l)
V̇ V sin 2 δ (s) − δ (l)
M ω̇ (s) = Pr(s) − D0 ω (s) −
2




2
2
(l)
− C̃q0 V
cos δ (s) − δ (l) δ̇ (s) − δ̇ (l)



Ck 
− (e) Vr(s) − V (s) V (l) sin δ (s) − δ (l) ,
Xd
(24)
Ck

where
Ku
Kf ,

C̃q0

and
C̃q0

and

are

constants,

τq0 (Xq −Xq0 )


.
(e) 2
Xq

=

Ck

with

=

The dynamic circuit

of the semi-damped model
is depicted
in Fig. 3,


(l)
with V̇d
= V (l) cos δ (s) (t) − δ (l) δ̇ (s) (t) − δ̇ (l) +


V̇ (l) sin δ (s) (t) − δ (l) .

I

q

 jI d  e j 

(s)

(t )

jX ( e )

jX d
+


 j X

Ck Vr

(s)

(e)
q

V

(s)

(l )
q' d

C V


e
e

j ( s ) ( t )
+
-

V

(s)
q

 jV

(s)
d

R(e)
+

e

j ( s ) ( t )

V ( l ) e j

(l )

(t )

j ( s ) ( t )

-

-

Fig. 3: Dynamic circuit of synchronous machine semi-damped
model.
IV. N UMERICAL VALIDATION
In this section, simulation results comparing the high-order
model, the classical model, the elemental model, the semidamped model, and the damped model, of a round-rotor
synchronous machine, are presented. We consider a two-bus
power system with a synchronous machine connected to a
constant power load through a short electrical transmission
line. See Fig. 4 for a one-line diagram, and Table II for the
system parameters.

7

(s)

V e

(s)
j  ( s ) ( t )t  ( t ) 



R

(e)

(l )

V e

jX

(l )
j  ( s ) ( t )t  



response, the damped model and the semi-damped model have
a higher accuracy for machine angular frequency response.
TABLE II: System parameters for a salient pole synchronous
machine

(e)

where:
 V ( s ) 
 ( s ) (t )  arctan  (ds )    ( s ) (t )
 Vq 



P ( l )  jQ ( l )
Damper windings

Fig. 4: One line diagram of a power system with a synchronous
machine connected to a constant power load through a short
transmission line.
Stator windings

A. Case 1
This case is used to highlight the high-fidelity of the secondorder models in comparison to the classical model, and we
consider the system response to an increase in real power
demand by the load. A round-rotor synchronous machine is
considered, and a stable equilibrium point for the high-order
model is chosen as the common initial condition for all the
models. The real power demand by the load is increased from
0.05 [pu] to 0.25 [pu] at time t = 30, and the reference voltage
(s)
magnitude Vr is changed at time t = 30 to keep the bus
voltage magnitude at V (l) = 1 [pu]. Numerical results are
depicted in Fig. 5.
Case 1 numerical results show that the elemental model,
the semi-damped model and the damped model have an overall
better accuracy than the classical model, and that after a certain
time range, the error of the classical model response increases
exponentially.

IEEE DC1A exciter

DEGOV1 speed governor

Transmission line

parameter
τq00
τd00
τq 0
Xq00
Xq0
Xq
Xk
ω0
Rs
Xd00
Xd0
τd 0
τf
τu
τ̄u
Xd
Kf
Ku
K̄u

value
0.9453 [s]
0.042 [s]
3.6123 [s]
0.2388 [pu]
0.7299 [pu]
1.7997 [pu]
0.19 [pu]
376.99 [rad/s]
0.003 [pu]
0.24 [pu]
0.32 [pu]
5.0141 [s]
1 × 10−8 [s]
0.002 [s]
1 × 10−12 [s]
1.7997 [pu]
1 [pu]
200
0 [s]

τ1
τ2
τ3
τ4
τ5
τ6
τm
κ
(s)
Pr
M
D̃0
D̄0

1 × 10−4 [s]
0 [s]
0.5001 [s]
25 × 10−3 [s]
9 × 10−4 [s]
5.74 × 10−3 [s]
24 × 10−3 [s]
10
0 [pu]
0.1188 [s2 ]
2.5825 × 10−7 [s/rad]
0.0531 [s/rad]

R(e)
X (e)

0.004 [pu]
0.0595 [pu]

B. Case 2
This case is used to compare the fidelity of the elemental
model, the semi-damped model and the damped model. The
machine whose parameters are described in Table II is employed. The real power demand by the load is increased from
0.05 [pu] to 0.25 [pu] at time t = 30, from 0.25 [pu] to 0.35
[pu] at time t = 1530, from 0.35 [pu] to 0.3 [pu] at time
t = 3030, and from 0.3 [pu] to 0.15 [pu] at time t = 4530.
(s)
For each load change, the reference voltage magnitude Vr
is changed to keep the bus voltage magnitude at V (l) = 1
[pu]. The root mean square errors of the models, relative to
the high-order model, are outlined in Table I, and numerical
results are presented in Figs. 6.
TABLE I: Root Mean Square Error (RMSE)
damped model
semi-damped model
elemental model

ω (s)
2.9093[rpm]
2.9093[rpm]
80.753[rpm]

V (s)
0.0042 [pu]
0.0042 [pu]
0.0042 [pu]

δ (s)
1065.6 [deg]
1065.6 [deg]
1065.6 [deg]

The RMSE results from Case 2 show that although the
elemental model, the semi-damped model and the damped
model match in accuracy for machine voltage magnitude

V. C ONCLUDING R EMARKS
In this paper, we introduced a library of second-order
synchronous machine models, comprising of the elemental
model, the damped model, and the semi-damped model. We
also described how these models, and the so-called classical
model, can be obtained from a high-order machine model.
While the classical model is obtained by identifying small
and large parameters in the high-order model, and setting
them to zero and infinity, respectively, the library of secondorder models are obtained by identifying fast and slow states
in the high-order model, and replacing differential equations
for the fast states with algebraic counterparts, referred to as
approximate manifolds (zero-order or first-order). The library
of second-order models were validated by comparing their
responses to those of a high-order model, and the classical
model, for given test cases.
We note that the second-order models presented in this
paper are not exhaustive. By developing higher-order approximations of manifolds for the fast states of the high-order
model, other second-order synchronous machine models can
be formulated. However, in this work we have focused on

8

ω (s) (rpm)

V (s) (per-unit)

δ (s) (degrees)

0

3500

high-order model
damped model
semi-damped model
elemental model
classical model

1.0045

3000

high-order model
damped model
semi-damped model
elemental model
classical model

2500
2000

1

1.004

-500

1.0035

0.995

1500 3600.05

1

3600

1000

0.99
3535
3530
3525

500
0

-1000

high-order model
damped model
semi-damped model
elemental model
classical model

0.985

0.12
0.1
0.08
0.06
0.04
0.02

-1500

1 2 3 4 5

-500

-2000
0

1

2

3

4

5

0

1

2

3

4

5

0

1

2

Time (Minutes)

Time (Minutes)

3

4

5

Time (Minutes)

Fig. 5: Case 1 numerical results: machine angular frequency, voltage magnitude and phase.
V (s) (per-unit)
1.005

0
high-order model
damped model
semi-damped model
elemental model

1.0045

δ (s) (degrees)

×10 4

×10 4

-1

-3.7219336543
-3.72193365435

-2
1.004
-3
high-order model
damped model
semi-damped model
elemental model

1.0036416

1.0035

1.0036414

0

10

20

30

-4

40

50

60

Time (Minutes)

70

80

90

100

0

10

20

30

40

50

60

70

80

90

100

Time (Minutes)

Fig. 6: Case 2 numerical results: machine angular frequency, voltage magnitude and phase.

first-order approximations of manifolds for damper winding
flux linkages.
[11]
R EFERENCES
[1] P. Kundur, N. J. Balu, and M. G. Lauby, Power system
stability and control. McGraw-Hill, 1994.
[2] P. Sauer and A. Pai, Power System Dynamics and Stability. Stipes Publishing L.L.C., 2006.
[3] P. Krause, O. Wasynczuk, S. Sudhoff, and S. Pekarek,
Analysis of Electric Machinery and Drive Systems, ser.
IEEE Press Series on Power Engineering. Wiley, 2013.
[4] L. Wang, J. Jatskevich, and H. W. Dommel, “Reexamination of synchronous machine modeling techniques for electromagnetic transient simulations,” IEEE
Transactions on Power Systems, vol. 22, no. 3, pp. 1221–
1230, Aug Aug. 2007.
[5] S. Crary, Power System Stability: Transient stability, ser.
General Electric series. John Wiley, 1947.
[6] E. Kimbark, Power Systems Stability. Vol. 3. Synchronous
Machines. Wiley, 1956.
[7] P. Kokotović, H. K. Khalil, and J. O’Reilly, Singular
Perturbation Methods in Control: Analysis and Design,
ser. Classics in Applied Mathematics.
Society for
Industrial and Applied Mathematics, 1986.
[8] H. K. Khalil, Nonlinear Systems. Pearson Education,
Limited, 2013.
[9] J. H. Chow, Time-Scale Modeling of Dynamic Networks
with Applications to Power Systems, B. A.V. and T. M.,
Eds., 1982.
[10] P. W. Sauer, S. Ahmed-Zaid, and P. V. Kokotovic, “An
integral manifold approach to reduced order dynamic
modeling of synchronous machines,” IEEE Transactions

[12]

[13]

[14]

on Power Systems, vol. 3, no. 1, pp. 17–23, Feb Feb.
1988.
P. V. Kokotovic and P. W. Sauer, “Integral manifold as a
tool for reduced-order modeling of nonlinear systems: A
synchronous machine case study,” IEEE Transactions on
Circuits and Systems, vol. 36, no. 3, pp. 403–410, March
Mar. 1989.
“IEEE Recommended Practice for Excitation System
Models for Power System Stability Studies,” IEEE Std
421.5-2016 (Revision of IEEE Std 421.5-2005), pp. 1–
207, Aug Aug. 2016.
PowerWorld corporation. (2017) Woodward diesel
governor model. [Online]. Available: https://www.
powerworld.com/WebHelp/Content/TransientModels
HTML/Governor%20DEGOV1.htm?tocpath=Transient%
20Stability%20Add-On%20(TS)%7CTransient%
20Models%7CGenerator%7CGovernor%7C
14
A. Fouad and V. Vittal, Power System Transient Stability
Analysis Using the Transient Energy Function Method.
Prentice Hall, 1992.

