International Journal of Electrical and Electronics Engineering 1:1 2008

Simulation of Synchronous Machine in Stability

Study for Power System
Tin Win Mon, and Myo Myint Aung

systems the state-space model has been used more frequently

Abstract—This paper presents a linear mathematical model of in connection with system described by linear differential
the synchronous generator with excitation system for power system equations. By the MATLAB modeling results of synchronous
stability. For large systems the state space model has been used more machine-infinite system stability is determined by root-locus
frequently in connection with system described by linear differential
plot.
equations. The flux linkage of each circuit in the machine depends
upon the exciter output voltage. The excitation system models This overall system consists of synchronous machine,
described use a pu system. The complete system is mainly included exciter, amplifier portion and the components of regulator.
synchronous machine, exciter and transmission line. The system is The mathematical modeling of this system is described in this
whether stable or unstable is determined by eigen-values of the paper and then this system is transferred into transfer function
system coefficient matrix A. These eigen-values are identified with Simulation result (time response) obtains from this transfer
the parameters of the machine and are not depend on the exciter
function.
parameters. Complete state equations are obtained first order
differential equation from the system model by identifying
appropriate state variables. The conversion of state space equation to II. MODELING OF SYNCHRONOUS MACHINE-INFINITE BUS
transfer function is applied in the MATLAB/SIMULINK program. SYSTEM
The first-order approximations are made for the system equations. To be modeling a Synchronous machine-infinite bus
The current and flux linkage are used to as state variables the overall
system, simple circuit of its electrical diagram as shown in
system. Stability characteristics of the system are determined by
examining the eigen-values. It can be prove that root locus plot. The Fig. 1 is considered. To be Modeling and Simulation of
current state-space model is specifically concerned with the 35 Synchronous Machine-infinite bus system, the following steps
MVA, 11 kV synchronous generator modeling aspects of the boost- to be made step by step;
buck excitation system. MATLAB/SIMULINK software is used to Step 1: Represent Synchronous Machine Circuit diagram
compute the state-space model of the system. Step 2: Represent System Mathematical equations
Step 3: Calculate the Transfer Function of overall system
Keywords—Park's transformation, Per-unit conversion, state- Step 4: Convert to model block
space model, transfer function.
Step 5: Create the M file to simulate the model
I. INTRODUCTION

C ONTROL system design and analysis technologies are

very useful to be applied in real time development. A
stable power system is one in which the synchronous
machines, when perturbed will either return to their original
state if there is no net change of power or will acquire a new
state without losing synchronism. If the system equations are
linear (have been linearized ), the techniques of linear system
analysis are used to study dynamic behavior. The most
common method is to simulate each component by its transfer
function. The system performance may then be analyzed by
such methods as root-locus plots, frequency domain analysis Fig. 1 A close loop system that representing the Synchronous
(Nyquist criterion), and Routh-Hurwitz criterion. These Machine and its excitation System
methods have been frequently used in studies pertaining to
small systems or a small number of machines. For large A. Description of Synchronous Machine
Large scale power is generated by three phase synchronous
generator, known as alternator, driven either by steam
Tin Win Mon is with the Mandalay Technological University Myanmar,
(corresponding author to provide phone: 095-02-88702; fax: 095-02-88702; e- turbines, hydro turbines, or gas turbines. The armature
mail: monmon.23@ gmail.com and nyinyi.29@ gmail.com). windings are placed on the stationary part called stator. The
Myo Myint Aung, was with Yangon Technological University, Myanmar. armature windings are designed for generation of balance
He is now with the Head of Department of Electrical Power Engineering,
Mandalay Technological University, Myanmar. three phase voltages and are arranged to develop the same

49
 International Journal of Electrical and Electronics Engineering 1:1 2008

number of magnetic poles as the field winding that is on the Ld = Equivalent direct-axis Reactance
rotor. The rotor is also equipped with one or more short- LF = Filed winding Self –inductance
circuited windings known as damper winding. The rotor is LD = Self-inductance damper winding
driven by prime mover at constant speed and its field circuit is Lq = Equivalent quadrature axis reactance
excited by direct current flow. The excitation may be provided LQ= Self inductance of quadrature
through slip rings and brushes by means of AC generators reactance
mounted on the same shaft as the rotor of the synchronous kMD = Stator to damper winding
machine. Modern excitation system usually use AC generator resistance
with rotating rectifiers, are known as brushless excitation. kMQ = Stator to quadrature winding resistance
r = Stator winding current
III. MATHEMATICAL MODELING OF A MACHINE WITH rF = Field winding resistance
INFINITE-BUS SYSTEM rD = resistance of d axis damper winding
The mathematical description of the synchronous machine rQ = resistance of q axis damper winding
is obtained if a certain transformation of variables is iq = armature current in the q direction
performed. Park's transformation is simply to transform all iF = Field current
stator quantities from phase a, b and c into equivalent dq axis iD = d axis damper winding current
new variables. The Park's transformation P is defined as (1); id = q axis damper winding current

This equation is in the state space form and can be written

in compact form as (5) and (6).
v =-(R+ωN)i - Li˙ (5)

i˙ = -L-1 (R+ωN)I - L-1 v (6)

Where,
(1)
Ld = Ls + Ms +3/2 Lm
Angle θ is given by,
Lq = Ls + Ms -3/2 Lm
θ = ωR t + δ + π/2 (2) R = resistance matrix
where, N = matrix of speed voltage inductance coefficients
ωR = rated angular frequency of L= matrix of constant inductances
synchronous machine The excitation system controls the terminal voltage of
δ =Synchronous torque angle synchronous machine. This system is also plays an important
role in the machine's mechanical oscillations, since its effects
Equation (3) is the generator voltage equation in the rotor the electrical power. Two transfer functions are of vital
frame of reference is described in Per-unit. The machine importance. The first of these is the generator transfer
equation in the rotor frame of reference becomes (4). function. The generator equations are nonlinear and the
transfer function is a linearized approximation of the
V= -ri - λ˙ + Vn (3) behaviour of the generator terminal voltage near a quiescent
operating point or equilibrium state. The load equations are
also nonlinear and reflect changes in the electrical output
quantities due to changes in terminal voltage. Fig. 2 (a) and
(b) show synchronous generator stator, field, and dq
equivalent circuit.

(4) (a) Stator and Field Circuit

Where,

50
 International Journal of Electrical and Electronics Engineering 1:1 2008

vqΔ = − KSin(δ0 − α)δΔ + ReiqΔ − ω0 Lei − id0 LeωΔ

dΔ
(9)
⋅
+ Lei
qΔ

Since all variables are now small displacement. These

equations are written in the following matrix form (10).
v =- Kx - Mx. (10)

These matrix M and K are included the transmission line

(b) dq axis equivalent circuit constants, synchronous machine parameters and the infinite
Fig. 2 Schematic diagram of Synchronous Machine bus voltage. Therefore, state space model of the system is as
(11) and (12);
IV. SYNCHRONOUS MACHINE CONNECTED TO AN INFINITE x˙ = -M-1 Kx – M-1v (11)
BUS TROUGH A TRANSMISSION LINE
The differential equations for one machine connected to an which is the same form as
infinite bus through a transmission line with impedance Ze =
Re + jωRLe. In equation (7), id and iq are dc currents obtained x˙ = Ax + Bu (12)
from the modified Park's transformation.
Stability characteristics may be determined by examining
vd = − 3Vα Sin(δ − α) + Reid + ωLeiq (11). x is an n vector denoting the states of the system and A
vq= 3Vα Cos(δ − α ) + Rei − ωLe i
is a coefficient matrix. The system inputs are represented by
(7)
q d the r vector u, and these input are related mathematically to
If machine connected to an infinite bus with local load at differential equations by an n × r matrix B. This description
machine terminal, the currents id and iq should be replaced by has the advantage that A may be time varying and u may be
the currents itd and itq. (8) is used to add the transmission line used to represent several inputs if necessary.
in the synchronous machine.
V. STATE-SPACE DESCRIPTION OF OVERALL SYSTEM
vd= − 3 VαSin(δ − α) + Rei + ωLeitq
td The excitation system controls the generated EMF of the
vq = 3 Vα Cos(δ − α ) + Rei − ωLeitd (8) generator and therefore controls not only the output voltage
tq but also the power factor and current magnitude as well. The
If the conditions at the machine terminals are given, the following equations (13) through (17) of the excitation control
voltage of the infinite bus is then determined by subtracting system transformed into S plane are given below. The
the appropriate voltage drops to the machine terminal voltage. corresponding model block diagram is shown in Fig. 3.
Otherwise, if the terminal conditions at the infinite bus are
given as the boundary conditions, the position of the d and q ⋅ ⎛k ⎞ ⎛ 1 ⎞
v1 = ⎜ R ⎟Vt −⎜ ⎟V (13)
axis currents, voltages and the machine terminal voltage can ⎜τ ⎟ ⎜τ ⎟ 1
then be determined. ⎝ R ⎠ ⎝ R ⎠

A. Linearization of Load Equation for a Synchronous ⎛k ⎞ ⋅ ⎛ 1 ⎞

Machine with Infinite Bus System v3⋅ = ⎜⎜ F ⎟⎟ E FD − ⎜⎜ ⎟⎟V3 (14)
When the Machine-Infinite bus system is subjected to a ⎝τ F ⎠ ⎝τ F ⎠
small load change, it tends to acquire a new operating state.
The first order approximations for these have been illustrated ⎛K ⎞ ⎛ 1 ⎞
in previous (7) and (8). Let the state variables xi and xj have v R⋅ = ⎜⎜ A ⎟⎟Ve − ⎜⎜ ⎟⎟VR (15)
the initial values xio and xjo. The changes in these variables are ⎝ τA ⎠ ⎝τ A ⎠
xiΔ and xjΔ. Therefore, (7) is used to linearize, the linearized
system equations are as (9). ⎛ 1 ⎞ ⎛ S + KE ⎞
⋅
E FD = ⎜⎜ ⎟⎟VR − ⎜⎜ E ⎟⎟ E FD (16)
⋅ ⎝τ E ⎠ ⎝ τE ⎠
v dΔ = − KCos(δ 0 − α)δ Δ + R e i dΔ + ω0 L e i qΔ + i q0 L e ω Δ

+ L e i dΔ VE = VREF + VS – V1 – V3 (17)

51
 International Journal of Electrical and Electronics Engineering 1:1 2008

3. The complete linear state space equation (4) is

transferred to output stator and rotor current by using
MATLAB programming.
4. Root-locus method is used to prove the system
performance is either satisfaction or not.
6. The process is continued until all the system roots are
located in the left hand s-plane.
Flow diagram of MATLAB programming step is shown in
Fig. 4.

VI. MODELING AND SIMULATION

Fig. 3 Closed Loop System that representing Synchronous Machine MATLAB/SIMULINK program was used in testing of the
and Excitation System
stability condition of the synchronous machine infinite bus
system. The M-file and SIMULINK model can be combined
The amplifier portion of excitation control system has an
by the commands [num, den]=ss2tf(A, B,C, D, 1);c=step(num,
amplification factor KA and time constant τA, supplying an
den, t); plot(t, c); title; x label and y label. The parameters for
output voltage VR which is limited by VRmax and VRmin. In the
35 MVA, 11kV generator, exciter and transmission line are
linearized model, the generator terminal voltage to its field
listed in Table I and II. This input parameters are used to
voltage can be represented by a gain kG and a time constant
simulation program. The simulation results are shown in Fig.5
τG. Typical values of KA are in the range of 10 to 400. The
to Fig. 9. The out put field current of synchronous machine is
amplifier time constant τA is in the range of 0.02 to 0.1
shown in Fig. 5. The d axis stator current and damper winding
second, and often is neglected. Gain KG may vary between 0.7
current of the machine is shown in Fig. 6 and 7. The response
to 1, and τG is between 1.0 and 2.0 second from full load to
is highly oscillatory, with a very large overshoot and a long
no-load. τR is very small, its range is between 0.01 to 0.06
settling time. It cannot have a small steady state error and a
second.
satisfactory transient response at the same time. The system is
A. Flow Diagram for Simulation for overall System whether stable or unstable is determined not only eigen-values
of overall system matrix A but also root locus plot. In Fig. 8,
all the poles of the system transfer function have negative real
part. Therefore, the system is stable. Fig. 9 is terminal voltage
step response (with (or) without feed back stabilizer).
Transmission line impedance is 0.02 + j0.4. The real parts
of all the eigen-values are negative, which means that the
system is stable under the conditions assumed in the
development of this model.
TABLE I
INPUT PARAMETERS OF MACHINE BASE VALUE
Rated MVA (MVA) 35
Rated voltage (kV) 11
Excitation voltage (V) 400
stator current (kA) 1.837
field current (A) 385
power factor 0.85
TABLE II
INPUT PARAMETERS OF EXCITER AND REGULATOR
proportionality constant KR 1.0
Regulator input fitter time τR 0.01
(sec)
Fig. 4 Flow Diagram of Simulation Step regulator gain KA 40
regulator amplifier time constant τA 0.05
The procedure for modeling of Synchronous machine with (sec)
feedback excitation system is as follows: exciter constant related to self excited -0.05
1. Mathematical models of synchronous generator with field KE
excitation systems are essential for the assessment of desired exciter time constant τE 0.5
performance requirements. (sec)
2. The per-unit value of the system transfer function matrix generator stabilizing circuit gain KF 0.04
is calculated from Equation (1 to 17). regulator stabilizing time constant τF 0.715

52
 International Journal of Electrical and Electronics Engineering 1:1 2008

Fig. 8 Root Locus Plot for overall System

Fig. 5 Output Field Current of Machine

Fig. 6 d axis stator current of Synchronous Machine

(a) without stabilizer feedback

Fig. 7 Output d axis Damper Current of Machine

(b) with stabilizer feedback

Fig. 9 Terminal Voltage Step Response of Synchronous Machine-

Infinite bus System

53
 International Journal of Electrical and Electronics Engineering 1:1 2008

VII. CONCLUSION Tin Win Mon made her first publication of International Paper at this paper
“Simulation of Synchronous Machine in Stability Study for Power System”.
The aim of this paper is to introduce Technicians to the
Modeling of Synchronous machine with excitation system on
stability studies and to use computer simulation as a tool for
conducting transient and control studies. Next to having an
actual system to experiment on, simulation is often chosen by
engineers to study transient and control performance or to test
conceptual designs. MATLAB/SIMULINK is used because of
the short learning curve that most students require to start
using it, its wide distribution, and its general-purpose nature.
This will demonstrate the advantages of using MATLAB
for analyzing steady state power system stability and its
capabilities for simulating transients in power systems,
including control behavior. This paper mentioned only a part
of the author’s research and approaching of her studies in
MATLAB software. The real application of this research
paper is stability of synchronous machine-infinite bus system
as mentioned in above.

ACKNOWLEDGMENT
Firstly the author would like to thank his parents: U Pann
Nyunt and Daw Nyein Nu for their best wishes to join the
PhD research. Special thanks are due to her supervisor/head of
Electrical Power Engineering Department, MTU, Union of
Myanmar: Dr Myo Myint Aung. The author would like to
express her thank to her teachers: Mr. Than Zaw Htwe and
Mr. Hla Myo Aung. The author greatly expresses her thanks
to all persons whom will concern to support in preparing this
paper.

REFERENCES
[1] Paul M. Aerson, A. A. Fouad, “Power System Control and Stability,”
IEEE Power Systems Engineering series , 1997.
[2] M. K. EI. Sherbiny and. D. M. Mehta, “Dynamic system stability,” IEEE
Trans; PAS- 92:1538- 46, 1973.
[3] Hadi Saadat, “Power System analysis” McGraw-Hill companies, Inc,
1999.
[4] Hadi Saadat, “Computational Aids in Control Systems using
MATLAB,” McGraw-Hill Companies, Inc, 1993.
[5] G. T. Heydt, “Identification and Tracking of Parameters for a Large
Synchronous Generator,” Final project Report, Department of Electrical
Power Engineering, Arizona State University, 2002.
[6] P. Kundur, “Power System Stability and Control,” McGraw-Hill
companies, Inc, 1994.
[7] Courtesy by Kinda Hydro Power Plant.
[8] Math Works, 2001, Introduction to MATLAB, the Math Works, Inc.
[9] Math Works, 2001, SIMULINK, the Math Works, Inc.
[10] Math Works, 2001, What is SIMULINK, the Math Works, Inc.

Tin Win Mon was born in 1980, November 29. A.GTI Certificate in 2000,
November. Graduated in 3rd November, 2003 with B.E (Electrical Power) and
finished Master degree on March, 2006 with M.E (Electrical Power). At
present, she is a Ph.D candidate of Electrical Power Engineering Department,
MTU, Mandalay, Myanmar. I’m trying to accomplish my Ph.D. thesis.
She served as a Demonstrator at Mandalay GTC from May, 2002 to
January, 2004 when she was attending the Special Engineering Course in
MTU. From February, 2004 to now, she promoted as an Assistant Lecturer of
Hpa-an Technological University, Department of Technical and Vocational
Education, Myanmar. Now, she is making her PhD research at Mandalay
Technological University (MTU), Myanmar with the title of “Simulation of
Synchronous Machine in Stability Study for Power System’.

54