ASSIGNMENT -02

NAME B. SAI KISHORE

VTU NO 10572

Year/Sem III YEAR / VI SEM

Course Code 1151EE110

Course Title POWER SYSTEM ANALYSIS

SCHOOL OF ELECTRICAL AND COMMUNICATION

DEPARTMENT OF ELECTRICAL AND ELECTRONICS

ENGINEERING

DECEMBER 2019 –MAY 2020

ACADEMIC YEAR 2019-2020
EVEN SEM
Problem 1: A synchronous generator & synchronous-motor each rated 25MVA,
11kV having 15% sub-transient reactance are connected through transformer & a line
as shown. The transformers are rated at 25MVA, 11/66kV & 66/11 kV with leakage
reactance of 10% each. The line has a reactance of 10% on a base of 25MVA and
66kV. The motor is drawing 15MW at 0.8pf leading & terminal voltage at 10.6kV
when a symmetrical 3 phase fault occurs at motor terminals.

- Find the sub transient current in generator, motor & at the fault location using
hand calculations.
Problem 2: Draw Reactance Diagram for the system shown below,

The ratings of various components are,

TL = 50+j200 Ω/ph

Choose the system base MVA as100 MVA and base voltage for G1 to be 2400V.
With the help of any programming technique like MATLAB /C++ verify the per
unit reactance values of the power system.

Gen-1 10kVA, 2400V, x”=0.2 pu

Gen-2 20kVA, 2400V, x”=0.2 pu
Motor 25kVA, 4 kV 2400/9600V x=0.15 pu
T1 40kVA, x=0.1 pu

T2 80kVA, 10kV/5kV, x=0.1 pu

PYTHON CODE :

n=int(input('Enter No of buses '))

g=int(input('Enter No of Generators '))

t=int(input('Enter No of transformers '))

tr=int(input('Enter No of transmission lines '))

m=int(input('Enter No of motors/loads '))

M=int(input('Enter the new Base MVA '))

K=float(input('Enter the new Base KV '))

for j in range(1,t+1):

print('Transformer:',j)

Y=float(input('Transformer is connected with which bus '))

Tm=float(input('Enter the Base MVA of transformer '))

Tk=float(input('Enter the Base KV of bus ' ))

Tx=float(input('Enter the given reactance of transformer ' ))

LT=float(input('Enter LT rating '))

HT=float(input('Enter HT rating '))

print('Base MVA of transformer ',Tm,'Base KV of bus ',Tk,'Reactance of transformer ',Tx,'LT

rating ',LT,'HT rating ',HT)

d=int(input('If you want to find transformer Base KV of LT side enter 1 if HT side enter
2 '))

s=int(input('If the LT side of the transformer is connected with bus press 1 if not press 0
'))

if s==1 & d==1:

q=float(input('Enter the Base KV of bus '))

Kt=(LT/q)*(LT/q)

XT=Tx*Kt*(M/Tm)

print('Per Unit Reactance of transformer ',XT,' per unit')

BH=q*(HT/LT)

print('Base KV of HT side =',BH)

elif s==0 & d==1:

Q=float(input('Enter the Base KV of LT side '))

BL=Q*(LT/HT)
 print('Base KV of HT side =',BL)

Kt=(LT/BL)*(LT/BL)

XT=Tx*Kt*(M/Tm)

print('Per Unit Reactance of transformer ',XT,' per unit')

elif d==0 & s==1:

Q=float(input('Enter the Base KV of LT side '))

BL=Q*(LT/HT)

print('Base KV of HT side =',BL)

Kt=(HT/BL)*(HT/BL)

XT=Tx*Kt*(M/Tm)

print('Per Unit Reactance of transformer ',XT,' per unit')

elif d==0 & s==0:

q=float(input('Enter the Base KV of bus '))

Kt=(HT/q)*(HT/q)

XT=Tx*Kt*(M/Tm)

print('Per Unit Reactance of transformer ',XT,' per unit')

BH=q*(LT/HT)

print('Base KV of LT side =',BH)

for i in range(1,g+1):

print('Generator:',i)

X=float(input('Generator is connected with which bus ' ))

Gm=float(input('Enter the Base MVA of generator '))

Gk=float(input('Enter the Base KV of Generator ' ))

Gx=float(input('Enter the given reactance of generator ' ))

print('Base MVA of generator ',Gm,'Base KV of Generator ',Gk,'reactance of generator ',Gx)

w=int(input('If any transformer is present on this bus press 1 if not press 0 '))

if w==1:

LT=float(input('Enter LT rating '))

HT=float(input('Enter HT rating '))

s=int(input('If the LT side of the transformer is connected with bus press 1 if not press 0 '))
 if s==1:

q=float(input('Enter the Base KV of LT side '))

Kg=(Gk/q)*(Gk/q)

XG=Gx*Kg*(M/Gm)

print('Per Unit Reactance of Generator ',XG,' per unit')

BH=q*(HT/LT)

print('Base KV of HT side = ',BH)

elif s==0:

Q=float(input('Enter the Base KV of HT side '))

Kg=(Gk/Q)*(Gk/Q)

XG=Gx*Kg*(M/Gm)

print('Per Unit Reactance of Generator ',XG,' per unit')

BL=Q*(LT/HT)

print('Base KV of LT side = ',BL)

if w==0:

Kg=(Gk/K)*(Gk/K)

XG=Gx*Kg*(M/Gm)

print('Per Unit Reactance of Generator ',XG,' per unit')

for k in range(1,tr+1):

print('Transmission line:',k)

Z=int(input('Transmission line is connected with which bus '))

trlxr=float(input('Enter the given Actual resistance of transmission line ' ))

trlxl=float(input('Enter the given Actual reactance of transmission line ' ))

print('Actual Reactance of transmission line ',trlxr,' + j ',trlxl)

w=int(input('If any transformer is present on this bus press 1 If not press 0 '))

if w==1:

LT=float(input('Enter LT rating '))

HT=float(input('Enter HT rating '))

s=int(input('If the LT side of the transformer is connected with bus press 1 If not press 0 '))

if s==1:

q=float(input('Enter the Base KV of HT side '))

 Ktr=(q*q)/M

XTrl1=trlxr/Ktr

XTrl2=trlxl/Ktr

print('Per Unit Reactance of transmission line ',XTrl1,' + j ',XTrl2,' per unit')

BH=q*(LT/HT)

print('Base KV of HT side = ',BH)

elif s==0:

Q=float(input('Enter the Base KV of LT side '))

BL=Q*(HT/LT)

print('Base KV of lT side = ',BL)

Ktr=(BL*BL)/M

XTrl1=trlxr/Ktr

XTrl2=trlxl/Ktr

print('Per Unit Reactance of transmission line ',XTrl1,' + j ',XTrl2,' per unit')

if w==0:

r=('Enter the Base KV of transmission line ')

Ktr=(r*r)/m

XTrl1=trlxr/Ktr

XTrl2=trlxl/Ktr

print('Per Unit Reactance of transmission line ',Xtrl,' + j ',XTrl2,' per unit')

for l in range(1,m+1):

print('Loads/Motors:',l)

G=int(input('Motor is connected with which bus '))

Mm=float(input('Enter the Base MVA of Motor '))

Mk=float(input('Enter the Base KV of Motor ' ))

Mx=float(input('Enter the given reactance of Motor ' ))

print('Base MVA of Motor ',Mm,'Base KV of Motor ',Mk,'reactance of Motor ',Mx)

w=int(input('If any transformer is present on this bus press 1 If not press 0 '))

if w==1:

LT=float(input('Enter LT rating '))

HT=float(input('Enter HT rating '))

 s=int(input('If the LT side of the transformer is connected with bus press 1 If not press 0
'))

if s==1:

q=float(input('Enter the Base KV of LT side '))

Kt=(Mk/q)*(Mk/q)

XM=Mx*Kt*(M/Mm)

print('Per Unit Reactance of load/motor ',XM,'per unit')

BH=q*(HT/LT)

print('Base KV of HT side = ',BH)

elif s==0:

Q=float(input('Enter the Base KV of HT side '))

Kt=(Mk/Q)*(Mk/Q)

XM=Mx*Kt*(M/Mm)

print('Per Unit Reactance of load/motor ',XM,' per unit')

BL=Q*(LT/HT)

print('Base KV of LT side = ',BL)

if w==0:

R=('Enter the Base KV Load/Motor ')

Kt=(R/K)*(R/K)

XM=Mx*Kt*(M/Mm)

print('Per Unit Reactance of load/motor ',XM,' per unit')

print('Please Draw the Reactance Diagram using Above Vlues')

Problem 3: With the help of MATLAB /C++ program determine the bus
impedance matrix of the power system shown below in which impedances are shown
in p.u. using bus building algorithm.

Solution :
Matlab code :

% % % ------Formulation Zbus matrix-------

clear all;

clc

disp('------Formulation Z-bus Matrix-------')

n=input('enter total number of busses including Ref busses = ');

zbus=zeros(n,n);

t=1;

while t==1;

zbus

s=menu('specify case no','new bus to reference bus','existing bus to new bus','between

existing busses','existing bus to reference bus','print','quit');

switch s

case{1}

zb=input('enter impedance value = ');

zbus=zb;

case{2}

k=input('enter starting bus number = ');

n=input('enter new bus number = ');

zb=input('Enter impedance value = ');

for i=1:n

if i==n

zbus(n,n)=zbus(k,k)+zb;

else

zbus(i,n)=zbus(i,k);

zbus(n,i)=zbus(k,i);

end

end

case{3}

a=input('enter first bus number = ');

b=input('enter second bus number = ');

 zb=input('enter impedance value = ');

m1=zb+zbus(a,a)+zbus(b,b)-(2*zbus(a,b));

ztemp=(1/m1)*((zbus(:,a))-(zbus(:,b)))*((zbus(a,:))-(zbus(b,:)));

zbus=zbus-ztemp;

case{4}

k=input('enter the old bus no = ');

zb=input('enter the impedance value = ');

m2=zbus(k,k)+zb;

ztemp=(1/m2)*zbus(:,k)*zbus(k,:);

zbus=zbus-ztemp;

case{5}

zbus

case{6}

'end program'

t=0;

end

end