Scribd
Upload
292 views4 pages

PSS Lab Manual

This document contains summaries of multiple experiments involving power system analysis calculations. Experiment 1 involves calculating line inductance (L) and capacitance (C) values given geometric parameters. Experiment 2 extends this to bundled conductors. Experiment 3 calculates generalized ABCD constants for a transmission line given system data, then uses this to determine voltage regulation and efficiency. Experiment 4 shows the bus admittance (YBUS) and inverse bus impedance (ZBUS) matrix calculations for a power system. Experiment 5 demonstrates the bus building algorithm to calculate the ZBUS matrix through sequential modifications. Experiment 6 applies the Gauss-Seidel method to iteratively solve a set of bus voltages given admittances and

Uploaded by

Gokul Chandrasekaran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
292 views4 pages

PSS Lab Manual

This document contains summaries of multiple experiments involving power system analysis calculations. Experiment 1 involves calculating line inductance (L) and capacitance (C) values given geometric parameters. Experiment 2 extends this to bundled conductors. Experiment 3 calculates generalized ABCD constants for a transmission line given system data, then uses this to determine voltage regulation and efficiency. Experiment 4 shows the bus admittance (YBUS) and inverse bus impedance (ZBUS) matrix calculations for a power system. Experiment 5 demonstrates the bus building algorithm to calculate the ZBUS matrix through sequential modifications. Experiment 6 applies the Gauss-Seidel method to iteratively solve a set of bus voltages given admittances and

Uploaded by

Gokul Chandrasekaran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

-------------------------------------------------------EXPERIMENT 2 Problem-1 -------------------------------------------------

Dab= input('Enter Dab value);


Dbc= input('Enter Dbc value);
Dca= input('Enter Dac value);
GMD=(Dab*Dbc*Dca)^(1/3);
GMRL= input('Enter GMRL(GMR) Value in ft (inch/12): ');
disp('Value of L in mH/KM:);
L = 0.2*log(GMD/GMRL)
dia= input('Enter conductor diameter in in);
r=dia/2; GMRC=r/12;
disp('Value of C in F/KM:);
C = 0.0556/log(GMD/GMRC)
----------------------------------------------------------EXPERIMENT 2 Problem-2 -----------------------------------------------
Dab= input('Enter Dab value:);
Dbc= input('Enter Dbc value:);
Dca= input('Enter Dac value:);
GMD=(Dab*Dbc*Dca)^(1/3);
dia=input('Enter conductor diameter in in);
r=dia/2;
d=input('Enter Bundle Spacing in inch);
Ds=input('Enter GMR(Ds) in inch);
GMRL=(sqrt(d*Ds))/12;
GMRC=(sqrt(d*r))/12;
disp('Value of L in mH/KM:);
L = 0.2*log(GMD/GMRL)
disp('Value of C in F/KM:);
C = 0.0556/log(GMD/GMRC)
--------------------------------------------------------------EXPERIMENT 2 Problem-3 -------------------------------------------
S = input('Enter row vector [S11, S22, S33] = ');
H = input('Enter row vector [H12, H23] = ');
d = input('Bundle spacing in inch = ');
dia = input('Conductor diameter in inch = '); r=dia/2;
Ds = input('Geometric Mean Radius in inch = ');
S11 = S(1); S22 = S(2); S33 = S(3); H12 = H(1); H23 = H(2);
a1 = -S11/2 + j*H12;
b1 = -S22/2 + j*0;
c1 = -S33/2 - j*H23;
a2 = S11/2 + j*H12;
b2 = S22/2 + j*0;
c2 = S33/2 - j*H23;
Da1b1 = abs(a1 - b1); Da1b2 = abs(a1 - b2);
Da1c1 = abs(a1 - c1); Da1c2 = abs(a1 - c2);
Db1c1 = abs(b1 - c1); Db1c2 = abs(b1 - c2);
Da2b1 = abs(a2 - b1); Da2b2 = abs(a2 - b2);
Da2c1 = abs(a2 - c1); Da2c2 = abs(a2 - c2);
Db2c1 = abs(b2 - c1); Db2c2 = abs(b2 - c2);
Da1a2 = abs(a1 - a2);
Db1b2 = abs(b1 - b2);
Dc1c2 = abs(c1 - c2);
DAB=(Da1b1*Da1b2* Da2b1*Da2b2)^0.25;
DBC=(Db1c1*Db1c2*Db2c1*Db2c2)^.25;
DCA=(Da1c1*Da1c2*Da2c1*Da2c2)^.25;
GMD=(DAB*DBC*DCA)^(1/3)
Ds = 2.54*Ds/100; r = 2.54*r/100; d = 2.54*d/100;
Dsb = (d*Ds)^(1/2); rb = (d*r)^(1/2);
DSA=sqrt(Dsb*Da1a2); rA = sqrt(rb*Da1a2);
DSB=sqrt(Dsb*Db1b2); rB = sqrt(rb*Db1b2);
DSC=sqrt(Dsb*Dc1c2); rC = sqrt(rb*Dc1c2);
GMRL=(DSA*DSB*DSC)^(1/3)
GMRC = (rA*rB*rC)^(1/3)
L=0.2*log(GMD/GMRL) % mH/km
C = 0.0556/log(GMD/GMRC) % micro F/km
-----------------------EXP 3 PGM -1 short transmission line --------------
VRLL=220;
VR=VRLL/sqrt(3);
Z=[0.15+j*2*pi*60*1.3263e-3]*40;
disp=('(a)')
SR=304.8+j*228.6;
IR=conj(SR)/(3*conj(VR));
IS=IR;
VS=VR+Z*IR;
VSLL=sqrt(3)*abs(VS)
SS=3*VS*conj(IS)
REG=((VSLL-VRLL)/(VRLL))*10
EFF=(real(SR)/real(SS))*100
disp=('(b)')
SR=304.8-j*228.6;
IR=conj(SR)/(3*conj(VR));
IS=IR;
VS=VR+Z*IR;
VSLL=sqrt(3)*abs(VS)
SS=3*VS*conj(IS)
REG=((VSLL-VRLL)/(VRLL))*100
EFF=(real(SR)/real(SS))*100

---------------------------------------------------------------------------
-----------------------EXP 3 PGM -2 MEDIUM transmission line --------------
clc
clear all
z=[];
z=input('Enter the length of line(km),Receiving voltage(v),Power(W),Powerfactor,L,R&C(per phase per km)');
l=z(1);vr=z(2);p=z(3);pf=z(4);L=z(5);R=z(6);c=z(7);
z=32+i*81;
y=i*0.00056;
disp('Medium Transmission line by nominal T method:')
A=1+(z*y)/2;C=y;B=z*(1+(z*y)/4);D=1+(z*y)/2;
disp('Generalized ABCD Constants are:')
disp(A);
disp(B);
disp(C);
disp(D);
ir=p/(1.732*vr*0.8)
vs=(A*vrph)+(ir*B);
fprintf('Sending end voltage Vs is %f\n',vs);
is=(vrph*C)+(D*ir);
fprintf('Sending end current Is is %f\n',is);
reg=((3*conj(vs)/conj(A))-v)*100/v;
fprintf('Voltage Regulation is %f\n',reg);
pows=3* conj(vs)* conj(is)*0.808;
eff=p*100/pows;
fprintf('Efficiency is %f\n',eff);
--------------------------------------EXP 4,5-------------YBUS,ZBUS(INVERSE YBUS)----------------------
clc
cler all
disp('');
b = input('Enter no. of buses:');
s = input('Enter no. of impedences:');
for i=1:s
sb(i)=input('Enter starting bus no.');
rb(i)=input('Enter receiving bus no.');
imp(i)=input('Enter impendance of bus:');
lc(i)=input('Enter line charging admittance:');
ybus=diag(0,b-1);
end
for i=1:s
k1=sb(i);
k2=rb(i);
adm(i)=1/imp(i);
ybus(k1,k1)=ybus(k1,k1)+adm(i)+lc(i);
ybus(k2,k2)=ybus(k2,k2)+adm(i)+lc(i);
ybus(k1,k2)=-adm(i);
ybus(k2,k1)=ybus(k1,k2);
end
disp('bus admittance matrix is:)
ybus
disp('bus impendance matrix is:)
z=inv(ybus)
----------------------------------------------------------------------------------------------------------------------------------------
----------------------------------------------------EXP 6 -------GAUSS SEIDAL Method--------------------------------------
y12=10-j*20;
y13=10-j*30;
y23=16-j*32;
V1=1.05+j*0;
V2=1+j*0;
V3=1+j*0;
S2=-2.566-j*1.102;
S3=-1.386-j*.452;
iter =0;
for I=1:10;
iter=iter+1;
V2 = (conj(S2)/conj(V2)+y12*V1+y23*V3)/(y12+y23);
V3 = (conj(S3)/conj(V3)+y13*V1+y23*V2)/(y13+y23);
disp([iter, V2, V3])
end
V2= .98-j*.06; V3= 1-j*.05;
I12=y12*(V1-V2); I21=-I12;
I13=y13*(V1-V3); I31=-I13;
I23=y23*(V2-V3); I32=-I23;
S12=V1*conj(I12); S21=V2*conj(I21);
S13=V1*conj(I13); S31=V3*conj(I31);
S23=V2*conj(I23); S32=V3*conj(I32);
I1221=[I12,I21]
I1331=[I13,I31]
I2332=[I23,I32]
S1221=[S12, S21 (S12+S13) S12+S21]
S1331=[S13, S31 (S31+S32) S13+S31]
S2332=[S23, S32 (S23+S21) S23+S32]
---------------------EXP 5 -------ZBUS using bus building algorithm-----------
clc
clear all
close all
g=[1.2 0.2 0.15 1.5 0.3];
z1=[g(1)];
disp(z1);
disp('TYPE I MODIFICATION')
z2=[g(1) g(1)
g(1) g(1)+g(2)];
disp(z2);
disp('TYPE II MODIFICATION')
z3=[g(1) g(1) g(1)
g(1) g(1)+g(2) g(1)+g(2)
g(1) g(1)+g(2) g(1)+g(2)+g(3)];
disp(z3);
disp('TYPE III MODIFICATION')
z4=[g(1) g(1) g(1) g(1)
g(1) g(1)+g(2) g(1)+g(2) g(1)+g(2)
g(1) g(1)+g(2) g(1)+g(2)+g(3) g(1)+g(2)+g(3)
g(1) g(1)+g(2) g(1)+g(2)+g(3) g(1)+g(2)+g(3)+g(4)]
z4=[z4];
disp(z4);
disp('Actual Zbus matrix is:')
n=4;
for i=1:1:n
for j=1:1:n
z4(i,j)=z4(i,j)-(z4(i,n)*z4(n,j)/z4(n,n));
end
end
z4(:,4)=[];
z4(4,:)=[];
disp(z4);
disp('TYPE IV MODIFICATION MATRIX')
l=4;
p=2;q=4;
for i=1:l-1
z4(l,i)=z4(p-1,i)-z4(q-1,i);
z4(i,l)=z4(l,i);
end
z4(l,l)=g(5)+z4(1,1)+z4(3,3)-2*z4(1,3);
disp(z4);
n=4;
for i=1:1:n
for j=1:1:n
z4(i,j)=z4(i,j)-(z4(i,n)*z4(n,j)/z4(n,n));
end
end
z4(:,4)=[];
z4(4,:)=[];
disp('THE REQUIRED ZBUS MATRIX IS:');
ZBUS=z4*1i;
disp(ZBUS);

You might also like