EE6501 Dec
EE6501 Dec
Visit OF ENGINEERING
: www.Easyengineering.net
DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING
EE6501 - POWER SYSTEM ANALYSIS
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
UNIT - II
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Scanned by CamScanner
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
EE6501 - Power System Analysis
Part - A Unit - III
1
Part A
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
4.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
5.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
7.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
3. In what type of fault the +ve sequence component of current is equal in magnitude but opposite
in phase to negative sequence components of current?
Line to line fault.
4. In which fault the negative and zero sequence currents are absent?
In three phase fault the negative and zero sequence currents are absent.
8. Name the fault in which positive, -ve and zero sequence component currents are equal.(May
2012)
In single line to ground fault the +ve, -ve and zero sequence component currents are equal.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
17. What are the symmetrical components of a three phase system? (May 2011)(Nov 2012) (Nov
2014) (Nov 2015)(May 2016)
1) Positive sequence 2) negative sequence 3) Zero sequence
19. List the various symmetrical and unsymmetrical faults in a power system.(May 2012)
Symmetricalfault:3 phase short circuit fault.
Unsymmetrical fault: i) single line to ground fault ii) line to line fault iii) double line to ground fault
iv)open conductor fault
21. What are the observations made from the analysis of various faults?(Nov 2013)
i) To check the MVA ratings of the existing circuit breakers, when new generation are added into a
system; ii) To select the rating for fuses, circuit breaker and switch gear in addition to setting up of
protective relays; iii) To determine the magnitudes of currents flowing throughout the power system at
various time intervals after a fault occurs.
22. Write the boundary conditions for single line to ground fault.(Nov 2013)
The boundary conditions are Va = 0; Ib=Ic=0
23.What are the features of zero sequence current?(May 2014)
As zero sequence currents in three phases are equal and of same phase, three systems operate like single
phase as regards zero sequence currents.Zero sequence currents flow only if return path is available
through which circuit is completed.
24.Write the symmetrical component current of phase ‘a’ in terms of 3ϕ currents.( May
2016).
25. State the reason why, the negative sequence impedance of a transmission line is taken as equal
to positive sequence impedance of the line. (May 2015).
A transmission line is a passive and bilateral device. By passive, we mean there are no voltage or current
sources present in the equivalent model of a transmission line. Bilateral means the line behaves the same
way regardless of the direction of the current. Note that although a single transmission line is bilateral.
Because of a transmission line’s passive and bilateral properties, the phase sequence of the applied
voltage makes no difference, as a-b-c (positive-sequence) voltages produce the same voltage drops as a-c-
b (negative-sequence) voltages. This means that the positive- and negative-sequence impedances of a
transmission line are identical, provided that the line is transposed. Transposition means physically
exchanging the position of each phase conductor along the length of the line such that conductor #1
occupies: position #1 for 1/3 of the line length, position #2 for 1/3 of the line length, and position #3 for
1/3 of the line length. Conductors #2 and #3 are rotated similarly.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
PART - B
1. A 50Hz, 13.2 KV, 15MVA alternator has X1 = X2 = 20% and X0 = 8% and the neutral is
grounded through a reactor of 0.5 ohm. Determine the initial symmetrical rms current in the
ground reactor when a double line to ground fault occurs at the generator terminals at a time
when the generator voltage was 12 kV.
Solution:
Base MVA=15MVA,Base KV=13.2KV
Ea(p.u)=0.909p.u
Z1=j0.2 p.u,Z2=j0.2p.u,Zg0=j0.08p.u
Xn=0.5Ω,Xn(p.u)=j0.043 p.u
Z0=Zg0+3Zn=j0.08+3xj0.043=j0.21p.u
Ia0=j1.47p.u Ib=-3.938+j2.197
Ic=3.938+j2.197 In=Ib+Ic=4.394∟90o
Base current=656 A
There fore In=2882.46A
2. A 3-phase, 10 MVA, 11KV, generator with solidity earthed neutral point supplies a feeder. The
relevant impedances of the generator and feeder in ohm are as below:
Generator Feeder
(a) +ve sequence j1.2 j1.0
(b )-ve sequence j0.8 j1.0
(c) zero sequence j0.4 j3.0
If the line to line fault occurs at the far end of the feeder, calculate the fault current.
Solution:
Ea=6350V
Z1=j2.2Ω
Z2=j1.8Ω
Z3=j3.4Ω
I1=I2=I0=-j858.10A
Fault Current,Ia=3I0=-j2574.32A.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
3. A salient pole generator is rated 20 MVA, 13.8 kV and has X1=0.25p.u X2=0.35p.u and
X0=0.1p.u. The neutral of the generator is solidly grounded. Compute fault current in the
generator and line to line to ground fault at its terminals. Neglect initial load on the generator.
The reactance of a generator are X’’=X2=0.15p.u and X0=0.05p.u. The generator ratings are 10
MVA, 6KV. The generator is star connected with neutral point grounded through a reactor of
0.5ohm reactance. Compute fault current in amps when a single line to ground fault occurs at
the generator terminals.
Solution:
Ia1= -j1.667p.u
Ia0=0 Ia=Ia0+Ia1+Ia2=0
Therefore
Ia=0
Ib= 2416∟180o
Ic= 2416∟0o
4. A11 kV, 30MVA alternator has Z1=Z2=-j0.2 pu and Z0=-j0.05 pu. A line to ground fault occurs
on the generator terminals. Determine the fault current and line to line voltages during faulted
conditions. Assume that the generator neutral is solidly grounded and the generator is operating
at no load and at the rated voltage during the occurrence of the fault. (May 2012)
Solution:
Base MVA=30MVA,Base Voltage=11KV
Base Current=1574.6A
Zf=0
Ia1= Ia2= Ia0=-j2.222pu
Fault Current=3Ia1=-j6.666p.u=10496.3∟-900A
Line to line voltages are
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Vab=5.6∟79.1oKV
Vbc=11∟270oKV
Vca=5.6∟100.9oKV
Actual value of line current=7873∟90oA
5. A single line diagram of a power network is shown in the figure. (May 2013)
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
6. A 25 MVA , 13.2 KV alternator with solidly grounded neutral has a sub transient reactance of
0.25 p.u. The negative and zero sequence reactance are 0.35 and 0.01 p.u. respectively. If a
double line-to-ground fault occurs at the terminals of the alternator, determine the fault
current and line to line voltage at the fault. (May 2014)
Solution:
Base MVA=25MVA,Base KV=13.2KV
Ea(p.u)=0.909p.u
Z1=j0.25 p.u,Z2=j0.35p.u,Zg0=j0.1p.u
Xn=0.1p.u,
Ia1= -j2.29p.u
Ia2=j1.225p.u
Ia0=j1.072p.u
Ib=-3.04+j1.6045
Ic=3.04+j1.6045
In=Ib+Ic=j3.208 p.u
In=4209A
7. Obtain the expression for fault current for a line to line fault taken place through an impedance
Zb in a power system.(Nov 2013)(May 2014)
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
8. A 30 MVA, 11Kv, 3ϕ synchronous generator has a direct sub transient reactance of 0.25 pu.
The negative and zero sequence reactance are 0.35 and 0.1pu respectively. The neutral of the
generator is solidly grounded. Determine the sub transient current in the generator and the line
to line voltages for sub transient conditions when a single line to ground fault occurs at the
generator terminals with the generator operating unloaded at rated voltages. (Nov 2015)
Solution:
Z1 = j0.25 p.u Z2 = j0.35 p.u Z0 = j0.1 p.u Prefault voltage V0=1∟00
Symmetrical components of Fault Current Ia1 = I a2 = Ia0 = -j1.4286 p.u
1
Fault Current in p.u = 3 Ia = -j4.2857 p.u
Base Current = 1574.59 Amp
Subtransient Current or Phase Currrent:
Ia= -j4.2858 x 1574.59 =-j6748.37 Amp Ib=0 Ic=0
Symmetrical Component of Bus Voltages for Phase a:
Va1 = -0.1429 p.u
Va2 = 0.6429 p.u
Va0 = -0.5 p.u
Subtransient phase voltages:
Va =0
Vb = -0.2144 – j0.9898
Vc = -0.2144 + j0.9898
Line to Line Volage :
Vab=0.2144 +j0.9898p.u
Vbc = -j1.9796p.u
Vac = -0.2144 –j 0.9898 p.u
9. Derive the expression to calculate the fault current for single-line-to-ground fault.
Let a 1LG fault has occurred at node k of a network. The faulted segment is then as shown in Fig. 8.2 where it
is assumed that phase-a has touched the ground through an impedance Zf. Since the system is unloaded before the
occurrence of the fault we have
I fb I fc 0 (8.1)
Vka Z f I fa (8.2)
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
1 1 1 I fa
1
I fa 012 1 a a 2 0 (8.3)
3
1 a 2 a 0
I fa
I fa 0 I fa1 I fa 2 (8.4)
3
This implies that the three sequence currents are in series for the 1LG fault. Let us denote the zero, positive and
negative sequence Thevenin impedance at the faulted point as Zkk0, Zkk1 and Zkk2 respectively. Also since the Thevenin
voltage at the faulted phase is Vf we get three sequence circuits that are similar to the ones shown in Fig. 7.7. We can
then write
Vka 0 Z kk 0 I fa 0
Vka1 V f Z kk 1I fa1 (8.5)
Vka 2 Z kk 2 I fa 2
Again since
Vka Z f I fa Z f I fa 0 I fa1 I fa 2 3Z f I fa 0
Vf
I fa 0 (8.7)
Z kk 0 Z kk1 Z kk 2 3Z f
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
10. A three-phase Y-connected synchronous generator is running unloaded with rated voltage when a 1LG fault
occurs at its terminals. The generator is rated 20 kV, 220 MVA, with subsynchronous reactance of 0.2 per unit.
Assume that the subtransient mutual reactance between the windings is 0.025 per unit. The neutral of the
generator is grounded through a 0.05 per unit reactance. The equivalent circuit of the generator is shown in Fig.
8.4. We have to find out the negative and zero sequence reactances.
Since no current flows in phases b and c, once the fault occurs, we have from Fig. 8.4
1
I fa j 4.0
j 0.2 0.05
Vn X n I fa 0.2
Va 0
Vb Ebn Vn j 0.025I fa 0.6 j 0.866 1.0536 124.72
Vc Ecn Vn j 0.025I fa 0.6 j 0.866 1.0536124.72
Therefore
0 0.4
Va 012 C 1.0536 124.72 0.7
1.0536124.72 0.3
From (7.38) we can write Z1 = j(Ls + Ms) = j0.225. Then from Fig. 7.7 we have
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
I fa 0 I fa1 I fa 2
Comparing the above two values with (7.37) and (7.39) we find that Z0 indeed is equal to j(Ls 2Ms) and Z2 is equal
to j(Ls + Ms). Note that we can also calculate the fault current from (8.7) as
1
I fa 0 j1.3333
j 0.225 0.225 0.15 3 0.05
11. Derive the expression to calculate the fault current for line-to-line fault.
The faulted segment for an L-L fault is shown in Fig. 8.5 where it is assumed that the fault has occurred at
node k of the network. In this the phases b and c got shorted through the impedance Zf. Since the system is unloaded
before the occurrence of the fault we have
I fa 0 (8.8)
I fb I fc (8.9)
0 0
1
I fa 012 C I fb a a 2 I fb (8.10)
I fb
3 2
a a I fb
We can then summarize from (8.10)
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
I fa 0 0
(8.11)
I fa1 I fa 2
Therefore no zero sequence current is injected into the network at bus k and hence the zero sequence remains a dead
network for an L-L fault. The positive and negative sequence currents are negative of each other.
Now from Fig. 8.5 we get the following expression for the voltage at the faulted point
Again
a 2
a V V
ka1 ka 2
Moreover since Ifa0 = Ifb0 = 0 and Ifa1 = Ifb2, we can write
I fb I fb1 I fb 2 a 2 I fa1 aI fb 2 a 2 a I fa1 (8.14)
Equations (8.12) and (8.15) indicate that the positive and negative sequence networks are in parallel. The sequence
network is then as shown in Fig. 8.6. From this network we get
Vf
I fa1 I fa 2 (8.16)
Z kk1 Z kk 2 Z f
12. Let us consider the same generator as given in Example 8.1. Assume that the generator is unloaded when a
bolted (Zf = 0) short circuit occurs between phases b and c. Then we get from (8.9) Ifb = Ifc. Also since the
generator is unloaded, we have Ifa = 0. Therefore from (7.34) we get
Also since Vbn = Vcn, we can combine the above two equations to get
1 120 1120
I fb I fc 3.849
j 0.45
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Then
0 0
I fa 012 C 3.849 j 2.2222
3.849 j 2.2222
1
I fa1 I fa 2 j 2.2222
j 0.225 j 0.225
Also since the neutral current In is zero, we can write Va = 1.0 and
Vb Vc Vbn 0.5
1.0 0
Va 012 C 0.5 0.5
0.5 0.5
13. Derive the expression to calculate the fault current for double-line-to-ground fault.
The faulted segment for a 2LG fault is shown in Fig. 8.7 where it is assumed that the fault has occurred at
node k of the network. In this the phases b and c got shorted through the impedance Zf to the ground. Since the system
is unloaded before the occurrence of the fault we have the same condition as (8.8) for the phase-a current. Therefore
1
I fa 0
3
I fa I fb I fc I fb I fc
1
3 (8.17)
3I fa 0 I fb I fc
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Therefore
I fa 0 I fa1 I fa 2 0 (8.23)
The Thevenin equivalent circuit for 2LG fault is shown in Fig. 8.8. From this figure we get
Vf Vf
I fa1
Z kk1 Z kk 2 Z kk 0 3Z f Z kk 2 Z kk 0 3Z f
l (8.24)
Z kk1
Z kk 2 Z kk 0 3Z f
The zero and negative sequence currents can be obtained using the current divider principle as
Z kk 2
I fa 0 I fa1 (8.25)
Z Z 3Z
kk 2 kk 0 f
Z kk 0 3Z f
I fa 2 I fa1 (8.26)
Z Z 3Z
kk 2 kk 0 f
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
14. Let us consider the same generator as given in Examples 8.1 and 8.2. Let us assume that the generator is
operating without any load when a bolted 2LG fault occurs in phases b and c. The equivalent circuit for this
fault is shown in Fig. 8.9. From this figure we can write
Fig. 8.9 Equivalent circuit of the generator in Fig. 8.4 for a 2LG fault in phases b and c.
Combining the above three equations we can write the following vector-matrix form
I b 3.849 j1.8182
I c 3.849 j1.8182
Hence
0 j1.2121
i fa 012 C 3.849 j1.8182 j 2.8283
3.849 j1.8182 j1.6162
We can also obtain the above values using (8.24)-(8.26). Note from Example 8.1 that
Then
1.0
I fa1 j 2.8283
j 0.225 j 0.3
j 0.225
j 0.525
j 0.3
I fa 2 I fa1 j1.6162
j 0.525
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
j 0.225
I fa 0 I fa1 j1.2121
j 0.525
and Vb = Vc = 0. Therefore
1.0909 0.3636
Va 012 C 0 0.3636
0 0.3636
15. Consider the network shown in Fig. 8.10. The system parameters are given below:
We shall find the fault current for when a (a) 1LG, (b) LL and (c) 2LG fault occurs at bus-2.
Let us choose a base in the circuit of the generator. Then the per unit impedances of the generator are:
X G1 X G 2 0.2, X G 0 0.075
X T 1 X T 2 0.1
The MVA base of the motor is 40, while the base MVA of the total circuit is 50. Therefore the per unit impedances of
the motor are
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
50 50 50
X M 1 X M 2 0.2 0.25, X M 0 0.1 0.125, X n 0.05 0.0625
40 40 40
For the transmission line
1102
Zbase 242
50
Therefore
24.2 60.5
X L1 X L 2 0.1, X L0 0.25
242 242
Let us neglect the phase shift associated with the Y/ transformers. Then the positive, negative and zero sequence
networks are as shown in Figs. 8.11-8.13.
Fig. 8.11 Positive sequence network of the power system of Fig. 8.10.
Fig. 8.12 Negative sequence network of the power system of Fig. 8.10.
Fig. 8.13 Zero sequence network of the power system of Fig. 8.10.
From Figs. 8.11 and 8.12 we get the following Ybus matrix for both positive and negative sequences
15 10 0 0
10 20 10 0
Ybus1 Ybus2 j
0 10 20 10
0 0 10 14
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Again from Fig. 8.13 we get the following Ybus matrix for the zero sequence
13.3333 0 0 0
0 14 4 0
Ybus0 j
0 4 14 0
0 0 0 3.2
0.075 0 0 0
0 0.0778 0.0222 0
Z bus0 j
0 0.0222 0.0778 0
0 0 0 0.3125
Z1 Z 2 j 0.18, Z0 j 0.0778
(a) Single-Line-to-Ground Fault: Let a bolted 1LG fault occurs at bus-2 when the system is unloaded with
bus voltages being 1.0 per unit. Then from (8.7) we get
1
I fa 0 I fa1 I fa 2 j 2.2841 per unit
j 2 0.18 0.0778
Also Ifb = Ifc = 0. From (8.5) we get the sequence components of the voltages as
V2 a 0 j 0.0778I fa 0 0.1777
V2 a1 1 j 0.18I fa1 0.5889
V2 a 2 j 0.18I fa 2 0.4111
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Va V2 a 0 0
V C 1 V 0.9061 107.11
b 2 a1
Vc V2 a 2 0.9061107.11
(b) Line-to-Line Fault: For a bolted LL fault, we can write from (8.16)
1
I fa1 I fa 2 j 2.7778 per unit
j 2 0.18
I fa 0 0
1
I fb C I fa1 4.8113
I fc I fa 2 4.8113
Finally the sequence components of bus-2 voltages are
V2 a 0 0
V2 a1 1 j 0.18I fa1 0.5
V2 a 2 j 0.18I fa 2 0.5
(c)Double-Line-to-Ground Fault: Let us assumes that a bolted 2LG fault occurs at bus-2. Then
1
I fa1 j 4.2676 per unit
j 0.18 Z eq
The zero and negative sequence currents are then computed from (8.25) and (8.26) as
j 0.18
I fa 0 I fa1 j 2.9797 per unit
j 0.18 0.0778
j 0.0778
I fa 2 I fa1 j1.2879 per unit
j 0.18 0.0778
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
I fa I fa 0 0
1
I fb C I fa1 6.657137.11
I fc I fa 2 6.65742.89
Furthermore the sequence components of bus-2 voltages are
V2 a 0 j 0.0778I fa 0 0.2318
V2 a1 1 j 0.18I fa1 0.2318
V2 a 2 j 0.18I fa 2 0.2318
16. Let us now assume that a 2LG fault has occurred in bus-4 instead of the one in bus-2.
Therefore
X1 X 2 j 0.1667, X 0 j 0.3125
Also we have
Hence
1
I fa1 j3.631 per unit
j 0.1667 Z eq
Also
j 0.1667
I fa 0 I fa1 j1.2631 per unit
j 0.1667 0.3125
j 0.3125
I fa 2 I fa1 j 2.3678 per unit
j 0.1667 0.3125
I fa I fa 0 0
1
I fb C I fa1 5.5298159.96
I fc I fa 2 5.529820.04
We shall now compute the currents contributed by the generator and the motor to the fault. Let us denoted the
current flowing to the fault from the generator side by Ig, while that flowing from the motor by Im. Then from Fig. 8.11
using the current divider principle, the positive sequence currents contributed by the two buses are
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
j 0.25
I ga1 I fa1 j1.2103 per unit
j 0.75
j 0.5
I ma1 I fa1 j 2.4206 per unit
j 0.75
Similarly from Fig. 8.12, the negative sequence currents are given as
j 0.25
I ga2 I fa 2 j 0.7893 per unit
j 0.75
j 0.5
I ma 2 I fa 2 j1.5786 per unit
j 0.75
Finally notice from Fig. 8.13 that the zero sequence current flowing from the generator to the fault is 0. Then
we have
I ga0 0
I ma 0 j1.2631 per unit
Therefore the fault currents flowing from the generator side are
I ma I ma 0 0.421090
I C 1 I 3.8512154.07
mb ma1
I mc I ma 2 3.851225.93
In the above two examples we have neglected the phase shifts of the Y/ transformers. However according to
the American standard, the positive sequence components of the high tension side lead those of the low tension side by
30, while the negative sequence behavior is reverse of the positive sequence behavior. Usually the high tension side of
a Y/ transformer is Y-connected. Therefore as we have seen in Fig. 7.16, the positive sequence component of Y side
leads the positive sequence component of the side by 30 while the negative sequence component of Y side lags that
of the side by 30. We shall now use this principle to compute the fault current for an unsymmetrical fault.
16. Let us consider the same system as given in Example 8.5. Since the phase shift does not alter the zero
sequence, the circuit of Fig. 8.13 remains unchanged. The positive and the negative sequence circuits must
however include the respective phase shifts. These circuits are redrawn as shown in Figs. 8.14 and 8.15.
Note from Figs. 8.14 and 8.15 that we have dropped the 3 vis-à-vis that of Fig. 7.16. This is because the per
unit impedances remain unchanged when referred to the either high tension or low tension side of an ideal transformer.
Therefore the per unit impedances will also not be altered.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Fig. 8.14 Positive sequence network of the power system of Fig. 8.10 including transformer phase shift.
Fig. 8.15 Negative sequence network of the power system of Fig. 8.10 including transformer phase shift.
Since the zero sequence remains unaltered, these currents will not change from those computed in Example
8.6. Thus
Now the positive sequence fault current from the generator Iga1, being on the Y-side of the Y/ transformer will lead
Ima1 by 30. Therefore
Finally the negative sequence current Iga2 will lag Ima2 by 30. Hence we have
Therefore
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
17. Let us consider the same power system as given in Example 1.2, the sequence diagrams of which are given in
Figs. 7.18 to 7.20. With respect to Fig. 7.17, let us define the system parameters as:
Let us choose the circuit of Generator 3 as the base, the base MVA for the circuit is 300. The base voltages are
then same as those shown in Fig. 1.23. Per unit reactances are then computed as shown below.
300
Generator G1: X 0.2 0.3 , X0 = 0.15
200
2
18
Generator G2: X 0.2 0.1312 , X0 = 0.0656
22.22
2
220
Transformer T1: X 0.1 0.121
200
2
25
Transformer T2: X 0.1 0.1266
22.22
2
22
Transformer T3: X 0.1 0.121
20
75 100
Line B-C: X1 X 2 0.5625 , X 0 0.75
133.33 133.33
75 100
Line C-D: X1 X 2 0.5625 , X 0 0.75
133.33 133.33
50 75
Line C-F: X1 X 2 0.375 , X 0 0.5625
133.33 133.33
Neglecting the phase shifts of Y/ connected transformers and assuming that the system is unloaded, we shall find the
fault current for a 1LG fault at bus-1 (point C of Fig. 7.17).
From Figs. 7.18 and 7.19, we can obtain the positive and negative sequence Thevenin impedance at point C as
(verify)
Similarly from Fig. 7.20, the Thevenin equivalent of the zero sequence impedance is
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
1
I fa 0 j1.0188 per unit
j 2 0.2723 0.4369
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
PART – A
1. Define Dynamic stability of a power system.
Dynamic stability is the stability given to an inherently unstable system by automatic control devices and
this dynamic stability is concerned with small disturbances lasting for times of the order of 10 to 30
seconds.
6. If two machines with inertia’s H1, H2 are swinging together, what will be the inertia of the
equivalent machine?
H1 and H2 is the Inertia constant of M1 and M2; G1 and G2 is the capacity of M1 and M2.
Hs is the equivalent inertia of M1 and M2 ; Gs is the equivalent capacity of M1 and M2.
7. On what basis do you conclude that the given synchronous machine has lost stability?
Following a sudden disturbance on a power system rotor speeds, rotor angular differences and power
transfer undergo fast changes whose magnitude is dependent on the severity of the disturbance. If these
disturbances leads to growing oscillations in the power system even after some period of time say more
than 30 seconds then system said are in asynchronous state and it has lost synchronism.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
11. A 50Hz, 4 pole turbo alternator rated at 20 MVA, 13.2 KV has as inertia constant H = 4 KW –
sec/ KVA. Find the K.E stored in the rotor at synchronous speed.
F = 50Hz. P = 4, G = 20 MVA, H = 4 KW – Sec/ KVA. Stored K.E = 4 x 20 = 80MJ.
12. Mention the methods used for the solution of swing equation.
Methods used for solution of swing equation are: Point by point method, Modified Euler‟s method
andRunge-kutta method.
17. Define critical clearing angle and time? (May 2011)(May 2012)(Nov 2012) (Nov 2014) (May15)
Critical clearing angle „c ‟corresponds to critical clearing timetc,in which the fault in the line is cleared
by the circuit breaker above which the system goes out of synchronism.
20. Distinguish between steady state and transient state stability.(Nov 2012)
Steady state stability is basically concerned with the ability of the system to restore back to its stable state
upon a small disturbance whereas the transient stability is concerned with large disturbances.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
23. Differentiate between voltage stability and rotor angle stability.(Nov 2013, Nov 2016)
Voltage stability is the ability of a power system to maintain steady acceptable voltage at all buses in the
system under normal operating conditions and after being subjected to a disturbance.
Rotor angle stability is the ability of interconnected synchronous machines of a power system to remain in
synchronism.
PART – B
1. Explain the equal area criterion for the stability studies in the power system.
Equal Area Criterion
1.0 Development of equal area criterion
As in previous notes, all powers are in per-unit.
I want to show you the equal area criterion a little differently than the book does it.
Let‟s start from Eq. (2.43) in the book.
2 H d 2
Pm Pe Pa
Re dt 2 (1)
Note in (1) that the book calls ωRe as ωR; however, we need to use 377 (for a 60 Hz system).
We can also write (1) as
2 H d
Pm Pe Pa
Re dt (2)
Now multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:
H d d
2 P P
Re dt m e
dt (3)
Note:
d (t ) 2
2 (t )
d (t )
dt dt (4)
Substitution of (4) into the left-hand-side of (3) yields:
H d 2 d
P P
Re dt m e
dt (5)
Multiply by dt to obtain:
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
d 2 Pm Pe d
H
Re (6)
Now consider a change in the state such that the angle goes from δ1 to δ2 while the speed goes from
ω1 to ω2. Integrate (6) to obtain:
22 2
P Pe d
H
2
d
Re 2
1 1
m
(7)
2 S BRe 1
(8a)
2
R2 1 2 1 2
J2 J1 Pm Pe d
S BRe 2 2 1
(8b)
Returning to (8), let ω1 be the speed at the initial moment of the fault (t=0+, δ=δ1), and ω2 be the
speed at the maximum angle (δ=δr), as shown in Fig. 1 below.
Note that the fact that we identify a maximum angle δ=δr indicates an implicit assumption that the
performance is stable. Therefore the following development assumes stable performance.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Pm1
Pe
Pm3
Pm
Pe1
Pm2
Pe3
Pe2
δ
δ 1 δ3 δc 90° δr δm 180°
Fig. 1
Since speed is zero at t=0, it remains zero at t=0+. Also, since δr is the maximum angle, the speed is
zero at this point as well. Therefore, the angle and speed for the two points of interest to us are
(note the dual meaning of δ1: it is lower variable of integration; it is initial angle):
δ1=δ1 δ=δr
ω1=0 ω2=0
P m Pe d 0
(9b)
1
Recalling that Pa=Pm-Pe, we see that (9b) says that for stable performance, the integration of the
accelerating power from initial angle to maximum angle must be zero. Recalling again (8b), which
indicated the left-hand-side was proportional to the change in the kinetic energy between the two
states, we can say that (9b) indicates that the accelerating energy must exactly counterbalance the
decelerating energy.
Inspection of Fig. 1 indicates that the integration of (9b) includes a discontinuity at the moment
when the fault is cleared, at angle δ=δc. Therefore we need to break up the integration of (9b) as
follows:
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
c r
P m Pe 2 d P e3 Pm d
(12)
1 c
Observing that these two terms each represent areas on the power-angle curve, we see that we have
developed the so-called equal-area criterion for stability. This criterion says that stable
performance requires that the accelerating area be equal to the decelerating area, i.e.,
A1 A2 (13)
where
c
A1 P m Pe 2 d (13a)
1
r
A2 P e3 Pm d
(13b)
c
Figure 2 illustrates.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Pm1
Pe
Pm3
A2
Pm
Pe1
A1
Pm2
Pe3
Pe2
δ
δ 1 δ3 δc 90° δr δm 180°
Fig. 2
Figure 2 indicates a way to identify the maximum swing angle, δr. Given a particular clearing angle
δc, which in turn fixes A1, the machine angle will continue to increase until it reaches an angle δr
such that A2=A1.
Let us consider a three-phase synchronous alternator that is driven by a prime mover. The equation of motion
of the machine rotor is given by
d 2
J Tm Te Ta (9.6)
dt 2
where
Neglecting the losses, the difference between the mechanical and electrical torque gives the net
accelerating torque Ta. In the steady state, the electrical torque is equal to the mechanical torque,
and hence the accelerating power will be zero. During this period the rotor will move at
synchronous speed s in rad/s.
The angular position is measured with a stationary reference frame. To represent it with
respect to the synchronously rotating frame, we define
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
st (9.7)
where is the angular position in rad with respect to the synchronously rotating reference frame.
Taking the time derivative of the above equation we get
d d
s (9.8)
dt dt
d
r
dt
d
r s (9.9)
dt
We can therefore conclude that the rotor angular speed is equal to the synchronous speed only
when ddt is equal to zero. We can therefore term ddt as the error in speed. Taking derivative of
(9.8), we can then rewrite (9.6) as
d 2
J 2 Tm Te Ta (9.10)
dt
d 2
J r 2 Pm Pe Pa (9.11)
dt
where Pm, Pe and Pa respectively are the mechanical, electrical and accelerating power in MW.
In steady state, the machine angular speed is equal to the synchronous speed and hence we can
replace r in the above equation by s. Note that in (9.13) Pm, Pe and Pa are given in MW.
Therefore dividing them by the generator MVA rating Srated we can get these quantities in per unit.
Hence dividing both sides of (9.13) by Srated we get
2 H d 2
Pm Pe Pa per unit (9.14)
s dt 2
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Equation (7.14) describes the behaviour of the rotor dynamics and hence is known as the swing
equation. The angle is the angle of the internal emf of the generator and it dictates the amount of
power that can be transferred. This angle is therefore called the load angle.
3. A 50 Hz, 4-pole turbogenerator is rated 500 MVA, 22 kV and has an inertia constant (H) of 7.5. Assume that
the generator is synchronized with a large power system and has a zero accelerating power while delivering a
power of 450 MW. Suddenly its input power is changed to 475 MW. We have to find the speed of the generator
in rpm at the end of a period of 10 cycles. The rotational losses are assumed to be zero.
We then have
d 2 s
2
Pm Pe 100 25 523.6 electrical deg/s 2
dt 2H 15
523.6
9.1385 electrical rad/s 2
180
Noting that the generator has four poles, we can rewrite the above equation as
d 2 9.1385
2
4.5693 mechanical rad/s 2
dt 2
4.5693
60 43.6332 rpm/s
2
The machines accelerates for 10 cycles, i.e., 2010 = 200 ms = 0.2 s, starting with a synchronous speed of 1500 rpm.
Therefore at the end of 10 cycles
The power-angle relationship has been discussed in Section 2.4.3. In this section we shall consider this
relation for a lumped parameter lossless transmission line. Consider the single-machine-infinite-bus (SMIB) system
shown in Fig. 9.1. In this the reactance X includes the reactance of the transmission line and the synchronous reactance
or the transient reactance of the generator. The sending end voltage is then the internal emf of the generator. Let the
sending and receiving end voltages be given by
We then have
The sending end real power and reactive power are then given by
This is simplified to
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
PS jQS
V1V2 sin j V12 V1V2 cos (9.3)
X
Since the line is loss less, the real power dispatched from the sending end is equal to the real power received
at the receiving end. We can therefore write
V1V2
Pe PS PR sin Pmax sin (9.4)
X
where Pmax = V1V2/X is the maximum power that can be transmitted over the transmission line. The power-angle curve
is shown in Fig. 9.2. From this figure we can see that for a given power P0. There are two possible values of the angle
0 and max. The angles are given by
P0
0 sin 1
Pmax (9.5)
max 180 0
5. A generator is connected to a constant voltage bus through an external reactance of 0.3 per unit. The
synchronous reactance of the generator is 0.2 per unit and the voltage magnitude of the constant voltage bus is
1.0 per unit with its angle being 0. The generator delivers 0.9 per unit power to the constant voltage bus when
the angle of its terminal voltage is 15. We have to determine the magnitude and angle of its internal emf.
Let us denote the angle of the terminal voltage by t, while its magnitude by Vt. Since the generator delivers
0.9 per unit power to the constant voltage bus of magnitude 1.0 per unit through a 0.3 per unit reactance, we can write
from (9.4)
Vt 1
sin t t sin 15
V
0.9
0.3 0.3
Solving the above equation we get Vt = 1.0432 per unit. Therefore the current flowing through the line is
1.043215 1
Is 0.9 j 0.0255 per unit
j 0.3
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Then the magnitude of the internal emf is 1.1082 per unit angle its angle is 23.96. It can be readily verified that the
power delivered is
1.1082 1
sin 23.96 0.9 per unit
0.3 0.2
6. Consider the system of Example 9.1. Let us assume that the system is operating with Pm = Pe = 0.9 per unit
when a circuit breaker opens inadvertently isolating the generator from the infinite bus. During this period the
real power transferred becomes zero. From Example 9.1 we have calculated 0 = 23.96 = 0.4182 rad and the
maximum power transferred as
1.1082 1
Pmax 2.2164 per unit
0.5
We have to find the critical clearing angle.
From (9.15) the accelerating area is computed as by note that Pe = 0 during this time. This is then given by
cr
A1 0.9 d 0.9
0.4182
cr 0.9 0.4182 0.9 cr 0.3764
To calculate the decelerating area we note that m = 0.4182 = 2.7234 rad. This area is computed by noting that Pe =
2.2164 sin() during this time. Therefore
2.7234
A2 2.2164 sin 0.9d
cr
0.0493
cr cos 1 1.5486 rad 88.73
2.2164
7. Consider the system in which a generator is connected to an infinite bus through a double circuit
transmission line as shown in Fig. 9.5. The per unit system reactances that are converted in a common base, are
also shown in this figure. Let us assume that the infinite bus voltage is 10. The generator is delivering 1.0 per
unit real power at a lagging power factor of 0.9839 to the infinite bus. While the generator is operating in steady
state, a three-phase bolted short circuit occurs in the transmission line connecting buses 2 and 4 very near to
bus 4. The fault is cleared by opening the circuit breakers at the two ends of this line. We have to find the
critical clearing time for various values of H.
Let the current flowing through the line be denoted by I. Then the power delivered to the infinite bus is
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Pe 1 Re 1.0 I 0.9839 I
The total impedance during the time when both the lines are operational, the impedance between the generator and the
infinite bus is j(0.3+0.1+0.1) = j0.5 per unit. Then the generator internal voltage is
Therefore the machine internal voltage is E = 1.2 per unit its angle is 24.625 or 0.4298 rad.
The pre-fault equivalent circuit is shown in Fig. 9.6 (a). From this figure we can write the power transfer
equation as
1.2
Pre fault : Pe sin 2.4 sin
0.5
Once the fault is cleared by opening of the breakers connected near buses 2 and 4, only of one the two lines will be
operational. Therefore during the post fault period, the impedance between the generator and the infinite bus is
j(0.3+0.2+0.6) = j0.5 per unit as shown in Fig. 9.6 (b). Then the post-fault power transfer equation is given by
1.2
Post fault : Pe sin 2.0 sin
0.6
It is to be noted that since one of the two lines is functional during the fault, the power transfer during the
fault will not be zero. The equivalent circuit during the fault is shown in Fig. 9.7. Since the fault has occurred very near
to bus-4, we can assume that this bus has been short circuited. We shall find the Thevenin equivalent of the portion of
the circuit to the right of the points A and B.
The circuit between buses 2 and 3 has been converted into an equivalent Y using -Y transformation. This is
shown in Fig. 9.8 (a). From this figure we find the Thevenin impedance as
0.05 0.025
X th 0.45 0.4667 per unit
0.075
Also the Thevenin voltage is then given by
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
j 0.025
Vth 10 0.33330 per unit
j 0.075
The reduced circuit is shown in Fig. 9.8 (b). From this circuit we can write the following power transfer equation
during the fault
1.2 0.3333
During fault : Pe sin 0.857 sin
0.4667
Fig. 9.8 Equivalent circuit during the fault: (a) -Y transformed and (b) Thevenin equivalent.
Three power-angle curves are shown in Fig. 9.9. From this figure we find that
1
max sin 1 2.618 rad
2
Fig. 9.9 Power-angles curves for the three modes of operation of the system of Example 9.4.
cr
A1 1 0.857 sin d
0.4298
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
2.618
A2 2 sin 1d
cr
0.323
cr cos 1 1.8573 rad 106.41
1.1429
As mentioned earlier, the critical clearing angle depends on the system network
configuration. The critical clearing time, however, is dependent on H and will vary as this
parameter varies. Usually numerical methods are employed for finding out the clearing time. We
shall however demonstrate the effect of H through a MATLAB program. The program uses the
built-in ordinary differential equation solver though which the swing equations are solved. The
results obtained are listed in Table 9.1. It can be seen that as the value of H is increases, the
clearing time is also increases, even though the clearing angle remains the same. This is obvious as
the value of H increases, the response of the system becomes more sluggish due to larger inertia.
Hence, the rotor takes more time to accelerate.
Table 9.1 Effect of H on critical clearing time
Consider plotting our function f(x(t)) as a function of t. What we want to do, based on (11), is to
obtain the area under the curve of f(x(t)) from t=kT-T to t=kT, as illustrated in Fig. 1.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
f(x(t))
f(x(kT)
)
f(x(kT-T))
kT- kT
T Fig. 1
Here is a proposed approach: Assume that f is constant throughout the interval at f(kT-T). This
assumption is clearly not good if kT is large, but it might be reasonable if kT is made small
enough.
This approach approximates the integral as the area shaded by the vertical lines in Fig. 2. One
observes that it misses the area shaded by the horizontal lines in Fig. 2.
f(x(t))
f(x(kT)
) error
f(x(kT-
T))
kT- kT
T Fig. 2
Analytically, (10) becomes
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
x(t)
error
x(kT)
x(kT-T)
kT- kT
T Fig. 3
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
One major problem with this method is that if T is too large, the error will propagate from one
time step to another.
This means that error incurred at one time step kT will add to the error incurred at later time steps.
Solution methods where this can happen are said to be numerically unstable.
However, occurrence of this phenomenon does depend on our choice of T. I will provide you with
a way to consider this issue.
Two algorithms that improve on the forward (Euler) method are predictor-corrector and Runge-
Kutta. We will look at both briefly.
3.1 Predictor-corrector method
This method is called the modified Euler in your text. The idea here is that we will take a step to
compute x(kT) (a predictor) and then we will use that calculation to recalculate that same step (the
corrector).
Step 1: Predict x(kT) using Euler to get x(kT):
f c
1
2
f ( x(kT T )) f ( x p (kT )) (33)
x c (kT ) x(kT T ) Tf c
x(kT T )
T
2
f ( x(kT T )) f ( x p (kT )) (34)
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Average these
two slopes
4
x p(kT) 3 2
1
x(kT-T)
x(t)
kT- kT
T Fig. 5
Understanding the method is facilitated by observing the sequence of points in Fig. 5. The
derivative f(x(kT-T)) is obtained at point 1. The predicted point xp(kT) is obtained at point 2. The
derivative of the predicted point f(xp(kT)) is obtained at point 3. The final point, point 4, is
obtained from averaging the two obtained derivatives f(x(kT-T)) and f(xp(kT)).
One can observe intuitively from Fig. 5 that the error will be reduced. This method can be shown
to be equivalent to considering up to the second derivative term in the Taylor series, therefore the
error is of O(T3).
This is a significant improvement over the Euler method, but it is still a numerically unstable
algorithm. In other words, for the predictor-corrector method, for a given minimum eigenvalue, T
can be larger than it can be in the Euler method, but it is still true that the algorithm may be
unstable if T is too large.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
x K Tf ( x(kT T ))
Start-point derivative
1 1 only
K
First interior derivative
x2 K 2 Tf x(kT T ) 1
2
K2
Second interior
x3 K 3 Tf x(kT T )
derivative
2
x4 K 4 Tf x(kT T ) K 3
Approx. end-point
derivative
2. Each Ki is not a derivative but rather an increment in the integration variable, i.e.,
4. The derivatives are computed at three different locations within the interval:
Figure 6 below illustrates the sequence of calculations, which can be understood by following the
single arrows from point 1 to point 2 to point 3, and then the double arrows from point 3 to point 4
to point 5, and then the triple arrows from point 5 to point 6 to point 7 to point 8.
Visit : www.Easyengineering.net
Visit : www.Easyengineering.net
Fig. 6
Once each of the increments K1-K4 are computed, then the final increment is obtained by taking a
weighted average of the four increments, where the middle increments are weighted heaviest,
according to (36).
x K1 2 K 2 2 K 3 K 4
1
(36)
6
The middle increments (K2 and K3) are weighted heaviest as they are computed based on slopes
that will be more representative of the slope in the interval than the beginning (K1) and final (K4)
increments.
The R-K method can be shown to be equivalent to considering up to the fourth derivative term in
the Taylor series, therefore the error is O(T5). Although this is a significant improvement over the
Euler or the P/C method, R-K is also a numerically unstable algorithm therefore the stability
domain.
Visit : www.Easyengineering.net