Ee 1404 Power System Lab Manual

DHANALAKSHMI COLLEGE OF ENGINEERING
DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING

EE 1404 POWER SYSTEM SIMULATION LABORATORY
LAB MANUAL/ OBSERVATION
2009 2010 ODD SEMESTER
NAME
:
REG
:
BRANCH
:
PREPARED BY V.BALAJI, M.Tech, (Ph.D)
Asst.Professor /EEE
PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE
Page 1
LIST OF EXPERIMENTS
1. COMPUTATION OF PARAMETERS AND MODELLINGOF TRANSMISSION LINES.
2.
FORMATION OF ADMITTANCE MATRICES.
3.
FORMATION OF IMPEDANCE MATRICES.
4.
SOLUTION OF POWER FLOW USING GAUSS-SEIDEL METHOD.
5.
SHORT CIRCUIT ANALYSIS.
6.
SOLUTION OF POWER FLOW USING NEWTON-RAPHSON METHOD.
7.
LOAD FREQUENCY DYNAMICS OF SINGLE AREA POWER SYSTEMS.
8.
LOAD FREQUENCY DYNAMICS OF TWO AREA POWER SYSTEMS.
9.
TRANSIENT AND SMALL SIGNAL STABILITY ANALYSIS SINGLE MACHINE INFINITE BUS SYSTEM
.
10.
ECONOMIC DISPATCH IN POWER SYSTEMS
Page 2
I CYCLE
1.
COMPUTATION OF PARAMETERS AND MODELLINGOF TRANSMISSION LINES.
2.
FORMATION OF ADMITTANCE MATRICES.
3.
FORMATION OF IMPEDANCE MATRICES.
4.
SOLUTION OF POWER FLOW USING GAUSS-SEIDEL METHOD.
5.
SHORT CIRCUIT ANALYSIS.
II CYCLE
6.
SOLUTION OF POWER FLOW USING NEWTON-RAPHSON METHOD.
7.
LOAD FREQUENCY DYNAMICS OF SINGLE AREA POWER SYSTEMS.
8.
LOAD FREQUENCY DYNAMICS OF TWO AREA POWER SYSTEMS.
9.
.
10.
ECONOMIC DISPATCH IN POWER SYSTEMS
Page 3
CONTENTS
S.NO
Date
Name of the Experiment
Page NO
Marks
Remarks
out of 10 /Signature
Page 4
S.NO
Date
Name of the Experiment
Page NO
Marks
Remarks
out of 10 /Signature
Page 5
COMPUTATION OF PARAMETERS AND MODELLING OF TRANSMISSION LINES

Expt.No : Date : AIM : (i) To determine the positive sequence line parameters L
and C per phase per kilometre of a three phase single and double circuit transmi
ssion lines for different conductor arrangements. (ii) To understand modeling an
d performance of medium lines. SOFTWARE REQUIRED: MATLAB 5.3 THEORY : Transmissi
on line has four parameters resistance, inductance, capacitance and conductance.
The inductance and capacitance are due to the effect of magnetic and electric f
ields around the conductor. The resistance of the conductor is best determined f
rom the manufactures data, the inductances and capacitances can be evaluated usi
ng the formula. Inductance: The general formula: L = 0.2 ln (Dm / Ds) Where, Dm
= geometric mean distance (GMD) Ds = geometric mean radius (GMR) I. Single phase
2 wire system GMD = D GMR = re-1/4 = r Where, r = radius of conductor II. Three
phase symmetrical spacing GMD = D GMR = re-1/4 = r Where, r = radius of conductor
III. Three phase Asymmetrical Transposed GMD = geometric mean of the three dist
ance of the symmetrically place conductors = 3DABDBCDCA
Page 6
GMR = re-1/4 = r Where, r = radius of conductors Composite conductor lines The in

ductance of composite conductor X., is given by Lx = 0.2 ln (GMD/GMR) where, GMD
=
mn
(Daa Dab).(Dna.Dnm)
GMR = n2 (Daa Dab.Dan ).(DnaDnb.Dnn) where, ra = ra e(-1/ 4) Bundle Conductors: The G

of bundle conductor is normally calculated GMR for two sub conductor c = (Ds * d
)1/2 GMR for three sub conductor Dsb =(Ds * d2)1/3 GMR for four sub conductor Ds
b = 1.09 (Ds * d3)1/4 where, Ds is the GMR of each subconductor and d is bundle
spacing Three phase Double circuit transposed: The inductance per phase in milli
henries per km is L = 0.2 ln (GMD / GMRL) mH/km where, GMRL is equivalent geome
tric mean radius and is given by GMRL = (DSADSBDSC)1/3 where, DSADSB and DSC are
GMR of each phase group and given by DSA = 4(Dsb Da1a2)2 = [Dsb Da1a2]1/2 DSB =
4(Dsb Db1b2)2 = [Dsb Db1b2]1/2 DSC = 4(Dsb Dc1c2)2 = [Dsb Dc1c2]1/2 where, Dsb =GM
R of bundle conductor if conductor a1, a2.. are bundle conductor. Dsb = ra1= rb1=
ra2 = rb2 = rc2 if a1, a2. are bundle conductor
Page 7
GMD is the equivalent GMD per phase & is given by GMD = [DAB * DBC * DCA]1/3 wher
e, DAB, DBC&DCA are GMD between each phase group A-B, B-C, C-A which are given b
y DAB = [Da1b1 * Da1b2 * Da2b1 * Da2b2]1/4 DBC = [Db1c1 * Db1c2 * Db2c1 * Db2c2]
1/4 DCA = [Dc1a1 * Dc2a1 * Dc2a1 * Dc2a2]1/4 Capacitance A general formula for e
valuating capacitance per phase in micro farad per km of a transmission line is
given by C = 0.0556/ ln (GMD/GMR) F/km Where, GMD is the Geometric mean distance wh
ich is same as that defined for inductance under various cases.
PROCEDURE: 1. Enter the command window of the MATLAB. 2. Create a new M file by
selecting File - New M File 3. Type and save the program in the editor window. 4
. Execute the program by either pressing Tools Run. 5. View the results.
EXERCISES: 1 A three-phase transposed line composed of one ACSR, 1,43,000 cmil,
47/7 Bobolink conductor per phase with flat horizontal spacing of 11m between ph
ases a and b and between phases b and c. The conductors have a diameter of 3.625
cm and a GMR of 1.439 cm. The line is to be replaced by a three-conductor bundl
e of ACSR 477,000-cmil, 26/7 Hawk conductors having the same cross sectional are
a of aluminum as the single-conductor line. The conductors have a diameter of 2.
1793 cm and a GMR of 0.8839 cm. The new line will also have a flat horizontal co
nfiguration, but it is to be operated at a higher voltage and therefore the phas
e spacing is increased to 14m as measured from the center of the bundles. The sp
acing between the conductors in the bundle is 45 cm. (a) Determine the inductanc
e and capacitance per phase per kilometer of the above two lines. (b) Verify the
results using the MATLAB program.
Page 8
PROGRAM :
[GMD, GMRL, GMRC] = gmd; L = 0.2*log(GMD/GMRL) C = 0.0556/log(GMD/GMRC)
MANUAL CALCULATIONS :
Page 9

Page 10
2. A three phase overhead line 200km long R = 0.16 ohm/km and Conductor diameter
of 2cm with spacing 4,5,6m transposed.Find A,B,C,D constants ,sending end volta
ge,current ,power factor and power when the line is delivering full load of 50MW
at 132kV ,0.8 pf lagging , transmission efficiency , receiving end voltage and
regulation. PROGRAM :
ab=input( value of ab ); bc=input( value of bc ); ca=input( value of ca ); pr=in
put( receving end power in mw ); vr=input( receving end voltage in kv ); pfr=inp
ut( receving end powerfactor ); l=input( length of the line in km ); r=input( re
sistance/ph/km ); f=input( frequency ); D=input( diameter in m ); rad=D/2; newra
d=(0.7788*rad); deq=(ab*bc*ca)^(1/3); L=2*10^(-7)*log(deq/newrad); C=(2*pi*8.854
*10^-12)/log(deq/rad); XL=2*pi*f*L*l*1000; rnew=r*l; Z=rnew+i*(XL); Y=i*(2*pi*f*
C*l*1000); A=1+((Y*Z)/2); D=A; B=Z; C=Y*(1+(Y*Z)/4); vrph=(vr*10^3)/1.732; irold
=(pr*10^6)/(1.732*vr*10^3*.8); k=sin(acos(pfr)); ir=irold*(pfr-(j*k)); vs=((A*vr
ph)+(B*ir)); is=((C*vrph)+(D*ir)); angle(vs); angle(is); f=angle(vs); u=angle(is
); PFS=cos(f-u); eff=((pr*10^6)/(3*abs(vs)*abs(is)*PFS))*100; reg=(((abs(vs)/abs
(A))-abs(vrph))/abs(vrph))*100; L C rnew A B C abs(vs) abs(is) angle(vs)*180/pi
angle(is)*180/pi PFS eff reg
Page 11
MANUAL CALCULATIONS:
Page 12
RESULT :
Page 13
FORMATION OF BUS ADMITTANCE MATRICES

Expt.No : Date : AIM: To determine the admittance matrices for the given power s
ystem network. SOFTWARE REQUIRED: MATLAB THEORY: Bus admittance is often used in
power system studies. In most of the power system studies it is required to for
m y- bus matrix of the system by considering certain power system parameters dep
ending upon the type of analysis. Y-bus may be formed by inspection method only
if there is no mutual coupling between the lines. Every transmission line should
be represented by - equivalent. Shunt impedances are added to diagonal element c
orresponding to the buses at which these are connected. The off diagonal element
s are unaffected. The equivalent circuit of Tap changing transformers is include
d while forming Y-bus matrix.
FORMATION OF Y-BUS MATRIX Generalised Y-bus = yii .. yid ydi ydd
where, Yii = Self admittance Ydi = Transfer admittance
PROCEDURE: 1. 2. 3. 4. 5. Enter the command window of the MATLAB. Create a new M
file by selecting File - New M File Type and save the program in the editor win
dow. Execute the program by either pressing Tools Run. View the results.
I. EXERCISE: (i) Determine the Y bus matrix and Z bus matrix for the power syste
m network shown in fig. (ii) Check the results obtained in using MATLAB.
Page 14
PROGRAM : z = [0 1 0 1.0 0 2 0 0.8 1 2 0 0.4 1 3 0 0.2 2 3 0 0.2 3 4 0 0.08]; Y

= ybus(z)
Page 15

Page 16
RESULT:
Page 17
FORMATION OF BUS IMPEDANCE MATRICES

Expt.No: Date : AIM : To determine the bus impedance matrices for the given powe
r system network. SOFTWARE REQUIRED: MATLAB 5.3 THEORY: FORMATION OF Z-BUS MATRI
X In bus impedance matrix the elements on the main diagonal are called driving p
oint impedance and the off-diagonal elements are called the transfer impedance o
f the buses or nodes. The bus impedance matrix are very useful in fault analysis
. The bus impedance matrix can be determined by two methods. In one method we ca
n form the bus admittance matrix and than taking its inverse to get the bus impe
dance matrix. In another method the bus impedance matrix can be directly formed
from the reactance diagram and this method requires the knowledge of the modific
ations of existing bus impedance matrix due to addition of new bus or addition o
f a new line (or impedance) between existing buses.
selecting File - New M File 3. Type and save the program in the editor Window 4.
Execute the program by either pressing Tools Run. 5. View the results. EXERCISE
: (i) Determine the and Z bus matrix for the power system network shown in fig.
(ii) Check the results obtained in using MATLAB.
Page 18
PROGRAM : z = [0 1 0 1.0 0 2 0 0.8 1 2 0 0.4 1 3 0 0.2 2 3 0 0.2 3 4 0 0.08]; Y

= ybus(z) Ibus = [-j*1.1; -j*1.25; 0; 0]; Zbus = inv(Y) Vbus = Zbus*Ibus MANUAL
CALCULATIONS:
Page 19

Page 20
RESULT:
Page 21
SOLUTION OF POWER FLOW USING GAUSS-SEIDEL METHOD Expt.No : Date :

AIM : To understand, in particular, the mathematical formulation of power flow m
odel in complex form and a simple method of solving power flow problems of small
sized system using GaussSeidel iterative algorithm
SOFTWARE REQUIRED: MATLAB 5.3 THEORY: The GAUSS SEIDEL method is an iterative al
gorithm for solving a set of non-linear load flow equations. The non-linear load
flow equation is given by 1 Vp
k+1
Pp j Qp (Vp )
k *
p-1
n k+1
= Ypp
- Ypq Vq
q=1
- Vq k
q=p+1
The reactive power of bus-p is given by
p-1 n q=p
QPk+1 =(-1) x Im
(Vpk)*
Ypq Vqk+1 + Ypq Vqk
q=1
selecting File - New M File 3. Type and save the program in the editor Window 4.
Execute the program by either pressing Tools Run. 5. View the results.
Page 22
EXERCISE : The figure shows the single line diagram of a simple 3 buses power sy
stem with generator at bus 1.The magnitude at bus 1 is adjusted to 1.05pu. The s
cheduled loads at buses 2 and 3 are marked on the diagram. Line impedance are ma
rked in pu. The base value is 100kVA. The line charging susceptances are neglect
ed. Determine the phasor values of the voltage at the load buses 2 and 3. Find t
he slack bus real and reactive power. Verify the result using MATLAB.
Program %Gauss Sedial clc; data=[1 1 2 10-j*20 2 1 3 10-j*30 3 2 3 16-j*32] elem
ents=max(data(:,1)); bus=max(max(data(:,2)),max(data(:,3))); y=zeros(bus,bus); f
or p=1:bus, for q=1:elements, if(data(q,2)==p|data(q,3)==p) y(p,p)=y(p,p)+data(q
,4); end end end
Page 23
for p=1:bus, for q=1:bus, if (p~=q) for r=1:elements if((data(r,2)==p&data(r,3)=

=q)|(data(r,2)==q&data(r,3)==p)) y(p,q)=-(data(r,4)); end end end end end a1=inp
ut( enter p2 in MW: ); b1=input( enter q2 in MVAR: ); a2=input( enter p3 in MW:
); b2=input( enter q3 in MVAR ); pu=input( enter the base value in MVA ); p2=(a1
/pu); q2=(b1/pu); p3=(a2/pu); q3=(b2/pu); dx1=1+j*0; dx2=1+j*0; v1=1.05; v2=1+j*
0; v3=1+j*0; iter=0; disp( iter v2 v3 ); while(abs(dx1)&abs(dx2)>=0.00001)&iter<
7; iter=iter+1; g1=(((p2-j*q2)/conj(v2))+(-y(1,2)*v1)+(-y(2,3)*v3))/y(2,2); g2=(
((p3-j*q3)/conj(v3))+(-y(1,3)*v1)+(-y(2,3)*g1))/y(3,3); dx1=g1-v2; dx2=g2-v3; v2
=v2+dx1; v3=v3+dx2; fprintf ( %g ,iter),disp([v2,v3]); end MANUAL CALCULATION
Page 24

Page 25
RESULT
Page 26
SOLUTION OF POWER FLOW USING NEWTON-RAPHSON METHOD Expt.No : Date : AIM : To det
ermine the power flow analysis using Newton Raphson method SOFTWARE REQUIRED : M
ATLAB THEORY : The Newton Raphson method of load flow analysis is an iterative m
ethod which approximates the set of non-linear simultaneous equations to a set o
f linear simultaneous equations using Taylors series expansion and the terms are
limited to first order approximation. The load flow equations for Newton Raphson
method are non-linear equations in terms of real and imaginary part of bus volt
ages. PP =
n q=1 n
ep(eq Gpq + fq Bpq) + fp (fq Gpq - eq Bpq)
QP =
q=1
f (eq Gpq + fq Bpq) - ep (fq Gpq - eq Bpq)
Vp 2 = ep2 + fp2 where, ep = Real part of Vp fp = Imaginary part of Vp Gpq, Bpq = C

onductance and Susceptance of admittance Ypq respectively.
EXERCISE
Page 27
1. Consider the 3 bus system each of the 3 line bus a series impedance of 0.02 +
j0.08 p.u and a total shunt admittance of j0.02 pu. The specified quantities at
the buses are given below : Bus Real load Reactive demand, QD 1 0 0.6 Load Real
power Reactive Power Voltage Specified V1=1.04 Unspecified V3 = 1.04
demand, PD 1 2 3 2 0 1.5
generation,PG Generation, QG 0.5 0 1 QG3 = ?
2. Verify the result using MATLAB
PROGRAM :
%NEWTON RAPHSON METHOD clc; gbus = [1 2.0 1.0 0.0 0.0 2 0.0 0.0 0.5 1.0 3 1.5 0.
6 0.0 0.0]; ybus = [5.882-j*23.528 -2.941+j*11.764 -2.941+j*11.764 -2.941+j*11.7
64 5.882-j*23.528 -2.941+j*11.764 -2.941+j*11.764 -2.941+j*11.764 5.882-j*23.528
]; t= 0.001 v1=1.04+j*0; v2=1+j*0; v3=1.04+j*0; del3=angle(v3); del1=angle(v1);
del2=angle(v2); %abs(ybus(2,1)) %abs(v2) for i=1:10 p2=(abs(v2)*abs(v1)*abs(ybus
(2,1))*cos((angle(ybus(2,1)))+del1del2))+abs(v2)* abs(v2)*abs(ybus(2,2))*cos((an
gle(ybus(2,2))))+(abs(v2)*abs(v3)* abs(ybus(2,3))*cos((angle(ybus(2,3))+del3-del
2)); q2=-(abs(v2)*abs(v1)*abs(ybus(2,1))*sin((angle(ybus(2,1)))+del1-del2))abs(v
2)*abs(v2)*abs(ybus(2,2))*sin((angle(ybus(2,2))))-(abs(v2)*abs(v3)* abs(ybus(2,3
))*sin((angle(ybus(2,3)))+del3-del2)); p3=(abs(v3)*abs(v1)*abs(ybus(3,1))*cos((a
ngle(ybus(3,1)))+del1del3))+abs(v3)*abs(v3)*abs(ybus(3,3))*cos((angle(ybus(3,3))
))+(abs(v2)*abs(v3)* abs(ybus(3,2))*cos((angle(ybus(3,2)))+del2-del3)); delp20=g
bus(2,4)-gbus(2,2)-p2; delp30=gbus(3,4)-gbus(3,2)-p3; delq20=gbus(2,5)-gbus(2,3)
-q2; J(1,1)=(abs(v2)*abs(v1)*abs(ybus(2,1))*sin((angle(ybus(2,1)))+del1del2))+(a
bs(v2)*abs(v3)*abs(ybus(2,3))*sin((angle(ybus(2,3)))+del3del2)); J(1,2)=-(abs(v2
)*abs(v3)*abs(ybus(2,3))*sin((angle(ybus(2,3)))+del3del2)); J(1,3)=(abs(v1)*abs(
ybus(2,1))*cos((angle(ybus(2,1)))+del1PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE
Page 28
del2))+2*(abs(v2)*abs(ybus(2,2))*cos((angle(ybus(2,2))))+(abs(v3)*abs(ybus(2,3)
)* cos((angle(ybus(2,3)))+del3-del2)); J(2,1)=-(abs(v3)*abs(v2)*abs(ybus(3,2))*s
in((angle(ybus(3,2)))+del2del3)); J(2,2)=(abs(v3)*abs(v1)*abs(ybus(3,1))*sin((an
gle(ybus(3,1)))+del1del3))+(abs(v3)*abs(v2)*abs(ybus(3,2))*sin((angle(ybus(3,2))
)+del2del3)); J(2,3)=(abs(v3)*abs(ybus(3,2))*cos((angle(ybus(3,2)))+del2-del3));
J(3,1)=(abs(v2)*abs(v1)*abs(ybus(2,1))*cos((angle(ybus(2,1)))+del1del2))-(abs(v
2)*abs(v3)*abs(ybus(2,3))*cos((angle(ybus(2,3)))+del2del3)); J(3,2)=(abs(v2)*abs
(v3)*abs(ybus(2,3))*cos((angle(ybus(2,3)))+del2del3)); J(3,3)=-(abs(v2)*abs(ybus
(2,1))*sin((angle(ybus(2,1)))+del1del2))-2*(abs(v2)*abs(ybus(2,2))*sin((angle(yb
us(2,2))))(abs(v3)*abs(ybus(2,3))* sin((angle(ybus(2,3)))+del3-del2)); end J inv
(J); A=[del2;del3;abs(v2)]; delA0=[delp20;delp30;delq20]; delA1=inv(J)*delA0; de
lA1; b0=abs(v2); A1=[del2;del3;b0]+delA1; A1-delA0; if((A1-delA0)<=t) break; del
2=A1(1,1); del3=A1(2,1); abs(v2)=A1(3,1); end A1
MANUAL CALCULATIONS :
Page 29

Page 30

Page 31
RESULT:
Page 32
SHORT CIRCUIT ANALYSIS Expt.No : Date : AIM : To become familiar with modelling
and analysis of power systems under faulted condition and to compute the fault l
evel, post-fault voltages and currents for different types of faults, both symme
tric and unsymmetric. PROGRAM REQUIRED: MATLAB 5.3 THEORY : Symmetrical Fault :
II. Three phase fault : From the thevenins equivalent circuit Vth Fault current If
= Z th Where Vth = Thevenins Voltage Z th = Thevenins Impedance Unsymmetrical Fau
lt : Single line to ground fault : Fault current If = Ia = 3Ia1 Ia1 = Ea Z1+Z2+Z
0
Line to line fault: Fault current If = Ia1(a2 a ) Ea Ia1 = Z1+Z2 Double Line to
ground fault : Fault current If = 2 Ia0 +( Ia1+ Ia2) (a2 + a ) Ea Ia1 = Z1 + Z0Z
2 Page 33
Z0 + Z2 Ia2 = (- Ia1) * Z0 Z0 + Z2 Ia0 = ( Ia1 Ia2) Fault MVA = 3 * If * Vpu wher

e, Ia1, Ia2 and Ia0 are positive, negative and zero phase sequence currents Z1 ,
Z2 and Zo are positive, negative and zero phase sequence impedances PROCEDURE: 1
. Enter the command window of the MATLAB. 2. Create a new M file by selecting Fi
le - New M File 3. Type and save the program. 4. Execute the program by either p
ressing Tools Run. View the results. EXERCISE :
Page 34
The one line diagram of a simple power system is shown in figure. The neutral of
each generator is grounded through a current limiting reactor of 0.25/3 per uni
t on a 100MVA base. The system data expressed in per unit on a common 100 MVA ba
se is tabulated below. The generators are running on no load at their rated volt
age and rated frequency with their emfs in phase. Determine the fault current fo
r the following faults. (a) A balanced three phase fault at bus 3 through a faul
t impedance Zf = j0.1 per unit. (b) A single line to ground fault at bus3 throug
h a fault impedance Zf = j0.1 per unit. (c) A line to line fault at bus3 through
a fault impedance Zf = j0.1 per unit. (d) A double line to ground fault at bus3
through a fault impedance Zf = j0.1 per unit. Item G1 G2 T1 T2 L12 L13 L23 Base
MVA 100 100 100 100 100 100 100 Voltage Rating kV 20 20 20/220 20/220 220 220 2
20 X1 0.15 0.15 0.10 0.10 0.125 0.15 0.25 X2 0.15 0.15 0.10 0.10 0.125 0.15 0.25
X0 0.05 0.05 0.10 0.10 0.30 0.35 0.7125
Verify the result using MATLAB program. PROGRAM : zdata1 = [0 0 1 1 2 1 0 0.25 2
0 0.25 2 0 0.125 3 0 0.15 3 0 0.25];
zdata0 = [0 1 0 0.40 0 2 0 0.10 1 2 0 0.30 1 3 0 0.35 2 3 0 0.7125]; zdata2 = zd
ata1; Zbus1 = zbuild(zdata1) Zbus0 = zbuild(zdata0) Zbus2 = Zbus1; symfault(zdat
a1,Zbus1) lgfault(zdata0, Zbus0, zdata1, Zbus1, zdata2, Zbus2) llfault(zdata1, Z
bus1, zdata2, Zbus2) dlgfault(zdata0, Zbus0, zdata1, Zbus1, zdata2, Zbus2)
Page 35
Page 36

Page 37

Page 38
RESULT
Page 39
LOAD FREQUENCY DYNAMICS OF SINGLE AREA POWER SYSTEMS Expt . No : Date AIM : To b
ecome familiar with modelling and analysis of the frequency and tie-line flow dy
namics of a power system without and with load frequency controllers (LFC) and t
o design better controllers for getting better responses. THEORY : Active power
control is one of the important control actions to be perform to be normal opera
tion of the system to match the system generation with the continuously changing
system load in order to maintain the constancy of system frequency to a fine to
lerance level. This is one of the foremost requirements in proving quality power
supply. A change in system load cases a change in the speed of all rotating mas
ses ( Turbine generator rotor systems) of the system leading to change in system
frequency. The speed change form synchronous speed initiates the governor contr
ol (primary control) action result in all the participating generator turbine un
its taking up the change in load, stabilizing system frequency. Restoration of f
requency to nominal value requires secondary control action which adjust the loa
d - reference set points of selected ( regulating) generator turbine units. The
primary objectives of automatic generation control (AGC) are to regulate system
frequency to the set nominal value and also to regulate the net interchange of e
ach areas to the scheduled value by adjusting the outputs of the regulating unit
s. This function is referred to as load frequency control(LFC). :
PROCEDURE : 1. Enter the command window of the MATLAB. 2. Create a new Model by
selecting File - New Model 3. Pick up the blocks from the simulink library brows
er and form a block diagram. 4. After forming the block diagram , save the block
diagram. 5. Double click the scope and view the result.
Page 40
EXERCISE: 1.An isolated power station has the following parameters Turbine time
constant T = 0.5sec Governor time constant g = 0.2sec Generator inertia constant H
= 5sec Governor speed regulation = R per unit The load varies by 0.8 percent fo
r a 1 percent change in frequency , i.e,D = 0.8 (a) Use the Routh Hurwitz array
to find the range of R for control system stability. (b) Use MATLAB to obtain th
e root locus plot. (c) The governor speed regulation is set to R = 0.05 per unit
.The turbine rated output is 250MW at nominal frequency of 60Hz. A sudden load c
hange of 50MW( PL = 0.2 per unit) occurs. (i) (ii) Find the steady state frequen
cy deviation in Hz. Use MATLAB to obtain the time domain performance specificati
ons and the frequency deviation step response.
Page 41
Page 42

Page 43
RESULT:
Page 44
LOAD FREQUENCY DYNAMICS OF TWO AREA POWER SYSTEMS Expt . No : Date AIM : To beco
me familiar with modelling and analysis of the frequency and tie-line flow dynam
ics of a two area power system without and with load frequency controllers (LFC)
and to design better controllers for getting better responses. THEORY: Active p
ower control is one of the important control actions to be perform to be normal
operation of the system to match the system generation with the continuously cha
nging system load in order to maintain the constancy of system frequency to a fi
ne tolerance level. This is one of the foremost requirements in proving quality
power supply. A change in system load cases a change in the speed of all rotatin
g masses ( Turbine generator rotor systems) of the system leading to change in s
ystem frequency. The speed change form synchronous speed initiates the governor
control (primary control) action result in all the participating generator turbi
ne units taking up the change in load, stabilizing system frequency. Restoration
of frequency to nominal value requires secondary control action which adjust th
e load - reference set points of selected ( regulating) generator turbine units.
The primary objectives of automatic generation control (AGC) are to regulate sy
stem frequency to the set nominal value and also to regulate the net interchange
of each areas to the scheduled value by adjusting the outputs of the regulating
units. This function is referred to as load frequency control(LFC). :
PROCEDURE: 1. Enter the command window of the MATLAB. 2. Create a new Model by s
electing File - New Model 3. Pick up the blocks from the simulink library browse
r and form a block diagram. 4. After forming the block diagram , save the block
diagram. 5. Double click the scope and view the result.
Page 45
EXERCISE : 1. A two area system connected by a tie line has the following parame
ters on a 1000MVA common base Area Speed Regulation Frequency sens.load coeff. In
ertia Constant Base Power Governor Time Constant Turbine Time Constant 1 R1=0.05
D1=0.6 H1=5 1000MVA g1 = 0.2sec T1 =0.5sec 2 R2=0.0625 D2=0.9 H2=4 1000MVA g1 = 0.
3sec T1 =0.6sec
The units are operating in parallel at the nominal frequency of 60Hz. The synchr
onizing power coefficient is computed from the initial operating condition and i
s given to be Ps = 2 p.u. A load change of 187.5 MW occurs in area1. (a) Derterm
ine the new steady state frequency and the change in the tie-line flow. (b) Cons
truct the SIMULINK block diagram and obtain the frequency deviation response for
the condition in part(a). SIMULINK BLOCK DIAGRAM :
Page 46
MANUAL CALCULATION:
Page 47

Page 48

Page 49
RESULT:
Page 50
Expt.No : Date :
AIM : To become familiar with various aspects of the transient and small signal
stability analysis of Single-Machine-Infinite Bus (SMIB) system
PROGRAM REQUIRED : MATLAB 5.3 THEORY : Stability : Stability problem is concerne
d with the behaviour of power system when it is subjected to disturbance and is
classified into small signal stability problem if the disturbances are small and
transient stability problem when the disturbances are large. Transient stabilit
y: When a power system is under steady state, the load plus transmission loss eq
uals to the generation in the system. The generating units run a synchronous spe
ed and system frequency, voltage, current and power flows are steady. When a lar
ge disturbance such as three phase fault, loss of load, loss of generation etc.,
occurs the power balance is upset and the generating units rotors experience ei
ther acceleration or deceleration. The system may come back to a steady state co
ndition maintaining synchronism or it may break into subsystems or one or more m
achines may pull out of synchronism. In the former case the system is said to be
stable and in the later case it is said to be unstable. Small signal stability:
When a power system is under steady state, normal operating condition, the syst
em may be subjected to small disturbances such as variation in load and generati
on, change in field voltage, change in mechanical toque etc., The nature of syst
em response to small disturbance depends on the operating conditions, the transm
ission system strength, types of controllers etc. Instability that may result fr
om small disturbance may be of two forms, (i) Steady increase in rotor angle due
to lack of synchronising torque. (ii) Rotor oscillations of increasing magnitud
e due to lack of sufficient damping torque. FORMULA : Reactive power Qe = sin(co
s-1(p.f)) S* Stator Current It = = Et* Voltage behind transient condition E1 = E
t + j Xd1It Et* Pe - jQe
Page 51
Voltage of infinite bus EB = Et - j( X3 + Xtr )It X1 X2 where, X3 = X1 + X2 Angu

lar separation between E1 and EB o = E1 EB Prefault Operation: X = j Xd1+ jXtr
+ X1 X2 X1 + X2 E1 x EB Power Pe = X o = sin Pe * X
-1
sino
E1 * EB
During Fault Condition: Pe = PEii = 0 Find out X from the equivalent circuit dur
ing fault condition Post fault Condition: Find out X from the equivalent circuit
during post fault condition E1 x EB Power Pe = X max = - o Pm Pe = sinmax Critical
Clearing Angle: Pm(max - o ) + P3maxcosmax - P2maxcoso Coscr = P3max - P2max sin
Page 52
Critical Clearing Time: 2H (cr - o) tcr = fo Pm Sec

PROCEDURE : 1. Enter the command window of the MATLAB. 2. Create a new M file by
selecting File - New M File 3. Type and save the program. 4. Execute the progra
m by either pressing Tools Run 5. View the results. EXERCISE : 1. A 60Hz synchro
nous generator having inertia constant H = 5 MJ/MVA and a direct axis transient
reactance Xd1 = 0.3 per unit is connected to an infinite bus through a purely re
active circuit as shown in figure. Reactances are marked on the diagram on a com
mon system base. The generator is delivering real power Pe = 0.8 per unit and Q
= 0.074 per unit to the infinite bus at a voltage of V = 1 per unit.
a) A temporary three-phase fault occurs at the sending end of the line at point
F.When the fault is cleared, both lines are intact. Determine the critical clear
ing angle and the critical fault clearing time. . b) Verify the result using MAT
LAB program. PROGRAM : Pm = 0.8; E = 1.17; V = 1.0; X1 = 0.65; X2 = inf; X3 = 0.
65; eacfault(Pm, E, V, X1, X2, X3) For b) Pm = 0.8; E = 1.17; V = 1.0; X1 = 0.65
; X2 = 1.8; X3 = 0.8; eacfault(Pm, E, V, X1, X2, X3)
Page 53
MANUAL CALCULATION:
Page 54
RESULT:
Page 55
ECONOMIC DISPATCH IN POWER SYSTEMS Expt.No : Date

AIM : To understand the fundamentals of economic dispatch and solve the problem
using classical method with and without line losses.
:
PROGRAM REQUIRED : MATLAB 5.3 THEORY : Mathematical Model for Economic Dispatch
of Thermal Units Without Transmission Loss: Statement of Economic Dispatch Probl
em In a power system, with negligible transmission loss and with N number of spi
nning thermal generating units the total system load PD at a particular interval
can be met by different sets of generation schedules {PG1(k) , PG2(k) , PGN(K) }; k =
1,2,..NS
Out of these NS set of generation schedules, the system operator has to choose t
he set of schedules, which minimize the system operating cost, which is essentia
lly the sum of the production cost of all the generating units. This economic di
spatch problem is mathematically stated as an optimization problem. Given : The
number of available generating units N, their production cost functions, their o
perating limits and the system load PD, To determine : The set of generation sch
edules, PGi ; i = 1,2N Which minimize the total production cost, Min ; FT = Fi (PGi
)
i=1 N
(1)
(2)
and satisfies the power balance constraint = PGi PD = 0
i=1 N
(3)
and the operating limits
Page 56
PGi,min PGi PGi, ,max The units production cost function is usually approximated
by quadratic function Fi (PGi) = ai PG2i + bi PGi + ci ; i = 1,2,.N where ai , bi
and ci are constants Necessary conditions for the existence of solution to ED p
roblem
(4)
(5)
The ED problem given by the equations (1) to (4). By omitting the inequality con
straints (4) tentatively, the reduce ED problem (1),(2) and (3) may be restated
as an unconstrained optimization problem by augmenting the objective function (1
) with the constraint multiplied by LaGrange multiplier, to obtained the LaGrang
e function, L as Min : L (PG1 ..PGN , ) = Fi(PGi) - [ PGi PD]
i=1 i=1 N N
(6)
The necessary conditions for the existence of solution to (6) are given by L / PGi
= 0 = dFi (PGi) / dPGi - ; i = 1, 2,..N L / = 0 = PGi PD
i=1 N
(7) (8)
The solution to ED problem can be obtained by solving simultaneously the necessa
ry conditions (7) and (8) which state that the economic generation schedules not
only satisfy the system power balance equation (8) but also demand that the inc
remental cost rates of all the units be equal be equal to which can be interpret
ed as incremental cost of received power.
When the inequality constraints(4) are included in the ED problem the necessary
condition (7) gets modified as dFi (PGi) / dPGi = PGi,min PGi PGi, ,max ____(9)
(10)
for
for PGi = PGi, ,max ` PGi = PGi, ,mi for Economic Schedule PGi = ( -bi)/ 2ai ; i=1
,2.N
Page 57
Incremental fuel cost = PD + ( bi/2ai ) / (1/2ai)

i=1 i=1 N N
(11)
PROCEDURE : 1. Enter the command window of the MATLAB. 2. Create a new M file by
selecting File - New M File 3. Type and save the program. 4. Execute the progra
m by either pressing Tools Run. 5. View the results.
EXERCISE :
1.The fuel cost functions for three thermal plants in $/h are given by C1 = 500
+ 5.3 P1 + 0.004 P12 ; C2 = 400 + 5.5 P2 + 0.006 P22 ; C3 = 200 +5.8 P3 + 0.009
P3 ;
2
P1 in MW P2 in MW P3 in MW
The total load , PD is 800MW.Neglecting line losses and generator limits, find t
he optimal dispatch and the total cost in $/h by analytical method. Verify the r
esult using MATLAB program. PROGRAM :
alpha = [500; 400; 200]; beta = [5.3; 5.5; 5.8]; gamma = [0.004; 0.006; 0.009];
PD = 800; DelP = 10; lamda = input( Enter estimated value of Lamda = ); fprintf
( ) disp([ Lamda P1 P2 P3 DP ... grad Delamda ]) iter = 0; while abs(DelP) >
= 0.001
Page 58
iter = iter + 1; P = (lamda - beta)./(2*gamma); DelP = PD - sum(P); J = sum(ones

(length(gamma),1)./(2*gamma)); Delamda = DelP/J; disp([lamda,P(1),P(2),P(3),DelP
,J,Delamda]) lamda = lamda + Delamda; end totalcost = sum(alpha + beta.*P + gamm
a.*P.^2)
MANUAL CALCULATION:
Page 59
2. The fuel cost functions for three thermal plants in $/h are given by C1 = 500
+ 5.3 P1 + 0.004 P12 ; C2 = 400 + 5.5 P2 + 0.006 P22 ; C3 = 200 +5.8 P3 + 0.009
P32 ; The total load , PD is 975MW. Generation limits: 200 P1 450 MW 150 P2 350
MW 100 P3 225 MW Find the optimal dispatch and the total cost in $/h by analyti
cal method. Verify the result using MATLAB program. PROGRAM : cost = [500 5.3 0.
004 400 5.5 0.006 200 5.8 0.009]; mwlimits = [200 450 150 350 100 225]; Pdt = 97
5; dispatch gencost P1 in MW P2 in MW P3 in MW
MANUAL CALCULATION:
Page 60

Page 61

Page 62
RESULT :
Page 63

Ee 1404 Power System Lab Manual

Uploaded by

Copyright:

Available Formats

Ee 1404 Power System Lab Manual

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Ee 1404 Power System Lab Manual

Uploaded by

Copyright:

Available Formats

What are the different experiments covered in the lab manual?

What are the different experiments covered in the lab manual?

What are the objectives of the experiment on computation of parameters and modelling of transmission lines?

What are the objectives of the experiment on computation of parameters and modelling of transmission lines?

DHANALAKSHMI COLLEGE OF ENGINEERING

DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING

COMPUTATION OF PARAMETERS AND MODELLING OF TRANSMISSION LINES

GMR = re-1/4 = r Where, r = radius of conductors Composite conductor lines The in

GMR = n2 (Daa Dab.Dan ).(DnaDnb.Dnn) where, ra = ra e(-1/ 4) Bundle Conductors: The G

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

FORMATION OF BUS ADMITTANCE MATRICES

PROGRAM : z = [0 1 0 1.0 0 2 0 0.8 1 2 0 0.4 1 3 0 0.2 2 3 0 0.2 3 4 0 0.08]; Y

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

FORMATION OF BUS IMPEDANCE MATRICES

PROGRAM : z = [0 1 0 1.0 0 2 0 0.8 1 2 0 0.4 1 3 0 0.2 2 3 0 0.2 3 4 0 0.08]; Y

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

SOLUTION OF POWER FLOW USING GAUSS-SEIDEL METHOD Expt.No : Date :

for p=1:bus, for q=1:bus, if (p~=q) for r=1:elements if((data(r,2)==p&data(r,3)=

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

Vp 2 = ep2 + fp2 where, ep = Real part of Vp fp = Imaginary part of Vp Gpq, Bpq = C

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

Z0 + Z2 Ia2 = (- Ia1) * Z0 Z0 + Z2 Ia0 = ( Ia1 Ia2) Fault MVA = 3 * If * Vpu wher

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

Voltage of infinite bus EB = Et - j( X3 + Xtr )It X1 X2 where, X3 = X1 + X2 Angu

Critical Clearing Time: 2H (cr - o) tcr = fo Pm Sec

ECONOMIC DISPATCH IN POWER SYSTEMS Expt.No : Date

Incremental fuel cost = PD + ( bi/2ai ) / (1/2ai)

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

iter = iter + 1; P = (lamda - beta)./(2*gamma); DelP = PD - sum(P); J = sum(ones

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

PREPARED BY V.BALAJI, M.Tech, (Ph.D) Asst.Professor /EEE

You might also like