Fault Analysis (Short Circuit Analysis)

FAULT ANALYSIS (SHORT CIRCUIT ANALYSIS)

PROGRAM:

clear all

clc

z data[0 1 0  0.15 0 2 0 0.075 1 2 0 0.1 2 3 0 0.1 3 1 0 0.1];

nl=zdata(:,1);nr=zdata(:,2);x=zdata(:,4);

nbr=length(zdata(:,1));nbus=max(max(n1),max(nr));

z=j*x;

y=ones(nbr,1)./z;

ybus=zeros(nbus,nbus);

for k=1:nbr;

ifnl(k)>0&&nr(k)>0

ybus(nl(k),nr(k))=ybus(nl(k),nr(k))-y(k);

ybus(nr(k),nl(k))=ybus(nl(k),nr(k));

end

end

for n=1:nbus

for k=1:nbr

ifnl(k)==n/nr(k)==n

ybus(n,n)=ybus(n,n)+y(k);

else

end

end

end

ybus;

zbus=inv(ybus);

zf=0;

v=[1 1 1];

vfa=[0 0 0];

selected bus=input('Select bus 1,2,3=');

I=selectedbus

If=v(I)/(zf+zbus(I,I));

for k=1:3

ifk~I

vf(k)=v(k)-(zbus(I,k)*v(k)/(zf+zbus(I,I)));

vfa(k)=vfa(k)+vf(k);

end

vf(I)=zf*v(I)/(zf+zbus(I,I));

vfa(I)=vfa((I)+vf(k);

    I12f=abs(vfa(1)-vfa(2))/x(3);

    I23f=abs(vfa(1)-vfa(3))/x(4);

    I23f=abs(vfa(2)-vfa(3))/x(5);

end

vfa

I12f

I13f

I23f

If=abs(If)