 Open Access
 Total Downloads : 3051
 Authors : K.Himaja, Dr T S Surendra, Dr.S. Tara Kalyani
 Paper ID : IJERTV1IS7309
 Volume & Issue : Volume 01, Issue 07 (September 2012)
 Published (First Online): 25092012
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Steady State Stability Analysis Of A Single Machine Power System Using Matlab
Steady State Stability Analysis Of A Single Machine Power System Using Matlab
K.Himaja1, Dr T S Surendra2, Dr.S. Tara Kalyani3
1,2 Dept.of.EEE, PDBVRIT, JNT University, Hyderabad, Andhra Pradesh, India.
3 Dept.of.EEE, JNTH University, Hyderabad, Andhra Pradesh, India.
ABSTRACT— This paper provides the stability of an interconnected power system is its ability to return to normal or stable operation after having been subjected to some form of disturbance. We have taken a single machine connected to a infinite busbar example to demonstrate the features of Steady state stability and the swing equation and its solution using MATLAB.
Keywords—: MATLAB, Modeling of a Power system, steady state Stability.

INTRODUCTION
The stability of power systems has been and continues to be of major concern in system operation. Modern electrical power systems have grown to a large complexity due to increasing interconnections, installation of large generating units and extrahigh voltage tielines etc.
Successful operation of a power system depends largely on the engineer's ability to provide reliable and uninterrupted service to the loads. The reliability of the power supply implies much more than merely being available. Ideally, the loads must be fed at constant voltage and frequency at all times. The first requirement of reliable service is to keep the synchronous generators running in parallel and with adequate capacity to meet the load demand.
Synchronous machines do not easily fall out of step under normal conditions. If a machine tends to speed up or slow down, synchronizing forces tend to keep it in step.
Conditions do arise, however, such as a fault on the network, failure in a piece of equipment, sudden application of a major load such as a steel mill, or loss of a line or generating Unit., in which operation is such that the synchronizing forces for one or more machines may not be adequate and small impacts in the system may cause these machines to lose synchronism.
There are two forms of instability in power system, the stalling of asynchronous loads and the loss of synchronism between synchronous machines. The synchronous stability is again divided into two regimes (i) Steady state stability and (ii) Transient state stability. The steady state stability is the stability of the power.
System under conditions of gradual or relatively slow change in load. The transient state stability refers to the maximum flow of power possible through a point without losing the stability with sudden and large changes in the network conditions such as brought about by faults, by sudden large increment of loads [1].
Transient stability is the ability of the power system to maintain synchronism when subjected to a severe transient disturbance, such as a fault on transmission facilities, sudden loss of generation, or loss of a large load. The system response to such disturbances involves large excursions of generator rotor angles, power flows, bus voltages, and other system variables. It is important that, while steadystate stability is a function only of operating conditions, transient stability is a function of both the operating conditions and the disturbance(s). [2]
This complicates the analysis of transient stability considerably. Repeated analysis is required for different disturbances that are to be considered. In the transient stability studies, frequently considered disturbances are the short circuits of different types. Out of these, normally the threephase short circuit at the generator bus is the most severe type, as it causes maximum acceleration of the connected machine. [3]
Power systems over the global are becomingly increasingly complex and the requirements for providing stable, secure, controlled, economic, quality power, especially so in the deregulated environment are becoming vitally important. An overwhelming need is being felt for increasing the transmission capacity on transmission lines and controlling power flow in specific corridors, while assuring system reliability during any fault scenarios. Power transfer capacity along a transmission corridor is limited by several factors such as thermal limit, steady state stability limit, transient stability limit and system damping. In certain situations, a power system may have inadequate damping or even negative damping.
All mutual inductances between stator and rotor circuits are periodic functions of rotor angular position. In addition, because of the rotor saliency, the self inductances of the stator phases and the mutual inductances between any two stator phases are also periodic functions of rotor angular position. Therefore, the resulting circuit equations are awkward to handle. However, if certain reasonable

The stator slots cause no appreciable variation of any of the rotor inductances with change in rotor position.

Saturation may be neglected. The effect of saturation can be accounted for separately.
If the magnitude of the disturbance is sufficiently small so that the system response in the initial stage is essentially linear, the stability may be classified as small disturbance stability. Small disturbance stability is assured if the eigen values of the appropriate dynamic model, linearized about the equilibrium point, have negative real parts. If there is an eigenvalue with positive real part the system is unstable. Complex eighenvalues occur in conjugate pairs. They signify the oscillations. With negative real parts oscillations damp out.
2. STEADY STATE STABILITY
2.1 .Swing equation
Under normal operating conditions, the relative position of the rotor axis and the resultant magnetic field axis is fixed. The angle between the two is known as the power angle or torque angle. During any disturbance, rotor will decelerate or accelerate with respect to the synchronously rotating air gap mmf, a relative motion begins. The equation describing the relative motion is known as the swing equation.
By the law of rotation
d 2
assumptions are made, a relatively simple
J m T T T
(1)
transformation of variables will eliminate the
dt 2
a m e
troublesome functions of angle from the equations.
The assumptions are:
Where J is the combined movement of inertia of the prime mover and generator
1. The windings are sinusoidally distributed along the air gap as far as all mutual effects with the rotor are concerned.
H
f0
d 2
dt 2 ) Pm( pu) Pe( pu) Pm Pmax sin
(2)
P dp P
cos
We have considered the dynamic behavior of a one
s d
max 0
(3)
machine system connected to an infinite bus bar as shown in Fig1. These the system and the given data widely used in
Swing equation in terms of
H d 2
the literature. [4]
The generator is delivering real power 0.6 p.u, unity power factor lagging to the infinite bus at a voltage V=1 per unit.
f0
dt 2
EV
Pm cos 0 0
(4)
The system data are given in Appendix I. The equations written in Matlab commands.
Pmax X
(5)
Synchronizing Coefficient is
Ps Pmax cos0
Undammed natural frequency is
(6)
n
f0 p H s
(7)
Fig 1 One Line diagram for Steady State stability
Damping ratio is
Before applying the disturbance the rotor angle and
D
2
f0
HPs
frequency of the synchronous machine is 9.80 and 60 Hz in Fig 2.
0 1 
2 
3 
4 
5 6 t, sec 
7 
8 
9 
10 

0 1 
2 
3 
4 
5 6 t, sec 
7 
8 
9 
10 
Damped frequency of oscillation is 11
Delta, degrees
10
d
Rotor is suddenly is perturbed by a small angle
0
1 2
ent sin( t )
9
(10) 8
Rotor frequency is
0 e
n 0
1 2
nt
sin d t
(11)
61
Frequency, Hz
60.5
60
The motion of rotor relative to the synchronously revolving field is
59.5
59
0
0
1 2
e nt sin( t )
d
(12)
Fig 2. The Rotor angle and frequency for machine
2.2. Illustrative system example
before disturbance.
By applying a small disturbance of per unit damping coefficient D=0.13 and a small disturbance of =120
.The breakers open and quickly close to the Line diagram. The Synchronous machine will be followed by a slowly damped oscillation and comes to the steady state operation in7.714 seconds in Fig 3. The rotor locks back into synchronous speed, the generator will maintain its stability at 7.714 seconds. It is correlated with theoretical results.
20
Delta, degrees
15
10
5
0
0 1 2 3 4 5 6 7 8 9 10
t, sec
By applying the input power is increased by a small amount P =0.3 per unit and the breakers open and quickly close to the Line diagram. The Synchronous machine will be followed by a slowly damped oscillation and comes to the steady state operation in7.714 seconds and a new steady state operating point is attained at =150 in Fig 4. The rotor locks back into synchronous speed, the generator will maintain its stability at 7.714 seconds. The stability is entirely independent of the input, and for a positive damping coefficient the system is always stable as long as the synchronizing power coefficient is positive.
20
Delta, degrees
15
10
5
0 1 2 3 4 5 6 7 8 9 10
t, sec
60.1
Frequency, Hz
60
59.9
59.8
0 1 2 3 4 5 6 7 8 9 10
t, sec
60.1
Frequency, Hz
60.05
60
59.95
0 1 2 3 4 5 6 7 8 9 10
t, sec
Fig 3. The Rotor angle and frequency for synchronous machine after disturbance
The motion of rotor relative to the synchronously revolving field in electrical radians becomes
Fig 4. The Rotor angle and frequency for Synchronous machine after small power input
f P
1 t
The qualitative conclusions regarding system stability
2
0 0 1 e
Hn 1 2
n sin(d t )
(13)
drawn from a twomachine or an equivalent onemachine infinite bus system can be easily extended to a multi machine system. Any parameter within Matlab
Rotor angular frequency in radian per second is
programming can be easily modified through simple MATLAB commands to suit the changes in the original
0
Hn
f0 P
1 2
e nt sin( t)
d
(14)
power system network due to fault or a corrective action. By applying a small disturbance the breakers open and quickly close to the Line diagram. The Synchronous machine will be followed by a slowly damped oscillation
and comes to the steady state operation the rotor locks back into synchronous speed, the generator will maintain its stability it is correlated with theoretical results.
References:
[1]C.L.Wadhwa, Electrical power systems (New Age International (P) Limited, 1994). [2]P. Kundur, Power System Stability and Control, EPRI Power System Engineering Series (Mc GrawHill, New York 1994). [3].I.J. Nagrath and D. P. Kothari, Power System Engineering (Tata McGrawHill, New Delhi, 1994). [4].Hadi Saadat, Power System Analysis (McGrawHill, New York, 1999). [5]. LouisA Dessaint et al., Power system simulation tool based on Simulink, IEEE Trans. Industrial Electronics, 46(6) (1999), 12521254. [6]. B.Stott, Power System Dynamic Response Calculations,Proc. IEEE, Vol.67, pp.210241, February 1979
[7].Khatir, M.; Zidi, S.A.; Fellah, M.K.; Hadjeri, S.; Dahou, O., “ HVDC Transmission Line Models For Simulink Environment , Proc.IEEE, IECON 2006, pp.436441, November 2006.Appendix I (generator data)
Generator 
RATING 
Power Factor 
1.0 
H(S) 
23.64 
X d 
0.0608 
Xt 
0.2 
X12 
0.0608 
12 Degrees 

f0 
60 