Power Systems Simulation Lab Assignment-4: H R Sai Kiran Pandit (116ee0328)
Power Systems Simulation Lab Assignment-4: H R Sai Kiran Pandit (116ee0328)
ASSIGNMENT-4
H R SAI KIRAN PANDIT(116EE0328)
Question-1: A 3Φ generator is connected to a Y-∆ transformer as shown in figure. The
positive negative and zero sequence impedances of alternator are 0.1j, 0.1j and 0.05j, for
transformer are 0.05j each and for the transmission line 0.4j, 0.4j and 0.8j respectively.
Determine the fault current and voltage at fault point when a single line to ground fault
occur at mid-point of line and alternator neutral is – (a) grounded, (b) isolated.
Code:
clc
clear
Xgen=[0.05i 0.1i 0.1i];
Xtransf=[0.05i 0.05i 0.05i];
Xtrline=[0.8i 0.4i 0.4i];
XTrFault=Xtrline/2;
E=1;
Z0=0.225i;
Z1=Xgen(2)+XTrFault(2);
Z2=Xgen(3)+XTrFault(3);
Neutral Grounded
fprintf('For grounded system\n\n')
If=(E*3)/(Z0+Z1+Z2)
If1=If/3;
ifmag=abs(If);
sequence voltages
Vf1=1-If1*Z1;
Vf2=-If1*Z2;
Vf0=-If1*Z0;
a=-0.5+(sqrt(3)/2)*1i;
Va=Vf0+Vf1+Vf2
Vb=Vf0+a^2*Vf1+a*Vf2
Vc=Vf0+a*Vf1+a^2*Vf2
Va = -2.2204e-16
Vb = -0.4091 - 0.8660i
Vc = -0.4091 + 0.8660i
For isolated system
fprintf('For isolated system\n\n')
eg=1; %generator emf
z1=0.3j;
z2=z1;
z0=0.45j; %considering zero sequence circuit of transformer
I1=eg/(z1+z1+z0);
If=3*I1;%fault current in first phase
%other phase currents are zero because no fault there
vseq=[-z0*I1;eg-z1*I1;-z2*I1];
alpha=exp(1j*2*pi/3);
conv_mat=[1,1,1;1,alpha^2,alpha;1,alpha,alpha^2];
vph=conv_mat*vseq;
fprintf('Fault current = %f + j(%f) A\n',real(If),imag(If));
fprintf('Va = %f + j(%f) V\nVb = %f + j(%f) V\nVc = %f + j(%f)
V\n',real(vph(1)),imag(vph(1)),real(vph(2)),imag(vph(2)),real(vph(3)),imag(vph(3)));
clc
clear
Z1=0.25i;
Z2=0.35i;
Z0=0.1i;
E=1;
Zeq=(Z2*Z0)/(Z2+Z0);
I1=E/(Z1+Zeq);
I2=-I1*(Z0/(Z0+Z2));
I0=-I1*(Z2/(Z0+Z2));
a=-0.5+(sqrt(3)/2)*1i;
Ia=I1+I2+I0;
Ib=I0+a^2*I1+a*I2;
Ic=I0+a*I1+a^2*I2;
Vbase=(13.8/sqrt(3))*10^3;
Sbase=10*10^6;
Ibase=Sbase/Vbase;
Ib1=Ib*Ibase
Ic1=Ib*Ibase
[theta_b,rhob]=cart2pol(-3.2292,3.5593);
[theta_c,rhoc]=cart2pol(3.2292,3.5593);
theta_b=theta_b*180/pi
theta_c=theta_c*180/pi
Ibf=rhob*Ibase
Icf=rhoc*Ibase
theta_b = 132.2161
theta_c = 47.7839
Ibf = 6.0319e+03
Icf = 6.0319e+03