EXPERIMENT-1
Aim: TO calculate the parameters of transmission line
Matlab Code for the Calculation of the parameters of transmission line
%analysis of transmission line
clc;
display('----------------------------------------------------------')
display('output for medium transmission line parameters-13g21a0246')
display('----------------------------------------------------------')
l=input('enter length of medium 3 phase transmission line in km:');
f=input('enter operating frequency in Hz:');
d=input('enter diameter of conductor in meters:');
r=d/2;
d12=input('enter distance between lines 1&2:');
d23=input('enter distance between lines 2&3:');
d31=input('enter distance between lines 3&1:');
R=input('enter resistance per km of a line:');
RT=R*l;
pr=input('enter recieving power in MW:');
vr=input('enter recieving end line voltage in kv:');
pf=input('enter power factor at output:');
w=2*pi*f;
L=2*(10^-7)*l*(10^3)*log((nthroot((d12*d23*d31),3))/(0.7788*r))
C=(2*pi*l*(10^3)*8.854*(10^-12))/log((nthroot((d12*d23*d31),3))/r)
Z=complex(RT,(w*L))
Y=complex(0,(w*C))
A=(1+((Y*Z)/2))
B=Z
C=Y*(1+((Y*Z)/4))
D=A
irold=(pr*(10^3))/((nthroot(3,2)*vr*pf))
ir=irold*(pf-1i*sin(acos(pf)))
vp=(vr*10^3)/(nthroot(3,2))
vs=(A*vp)+(B*ir)
is=(C*vp)+(D*ir)
anglevs=angle(vs)*(180/pi)
angleis=angle(is)*(180/pi)
inputpf=cos(angle(vs)-angle(is))
ps=3*abs(vs)*abs(is)*inputpf
efficiency=((pr*10^6)/ps)*100
regulation=((abs(vs)-abs(vp))/abs(vp))*100
Alternate Program :
clc
clear all
close all
L=input('Enter the length of the conductor:');
phase=input('Enter Phase:');
Pr = input('Receiving end Power: ');
Vr = input('Receiving end Volatge: ');
PF = input('Power Factor: ');
r = input('Resistance per Km: ');
x = input('Reactance Per Km: ');
Susceptance=input('The given value of susceptance:');
R = r*L;
X = x*L;
Z = R + j*X;
S=Susceptance*L;
Ic=j*S*Vr;
Method=input('Enter the desirable Method: \n 0-End Condenser Method \n 1-
Nominal-T \n 2-Nominal pi:\n')
if (Method==0);
if(phase==1)
I = Pr/(Vr*PF);
I = I*(PF - j*sin(acos(PF)));
Is=I+Ic;
Vs = Vr + Is*Z;
reg = (abs(Vs)-Vr)*100/Vr
Ps = Pr+(abs(Is)^2*R);
Efficiency = Pr*100/Ps
elseif(phase==3)
Vrph = Vr/sqrt(3);
I = Pr/(3*Vrph*PF);
I = I*(PF - j*sin(acos(PF)));
Is=I+Ic;
Vsph = Vrph + Is*Z;
Vs = sqrt(3)*Vsph;
reg = (abs(Vs)-Vr)*100/Vr
Ps = Pr + 3*(abs(Is)^2)*R;
Efficiency = Pr*100/Ps
else
disp('invalid phase');
end
elseif (Method==1);
if (phase==3)
Vrph = Vr/sqrt(3);
I = Pr/(3*Vrph*PF);
I = I*(PF - j*sin(acos(PF)));
V1=Vrph+(I*(Z/2));
Ic=j*S*V1;
Is=I+Ic;
Vsph=V1+Is*(Z/2);
Vs=sqrt(3)*Vsph;;
reg = (abs(Vs)-Vr)*100/Vr
Losses=3*(abs(Is)^2)*(R/2)+3*(abs(I)^2)*(R/2);
Efficiency= ((Pr)/((Pr)+Losses))*100
elseif (phase==1)
I = Pr/(Vr*PF);
I = I*(PF - j*sin(acos(PF)));
V1=Vr+(I*(Z/2));
Ic=j*S*V1;
Is=I+Ic;
Vs=V1+Is*(Z/2);
reg = (abs(Vs)-Vr)*100/Vr
Losses=(abs(Is)^2)*(R/2)+(abs(I)^2)*(R/2);
Efficiency= ((Pr)/((Pr)+Losses))*100
else
disp('Invalid Phase');
end
elseif(Method==2)
if (phase==3)
Vrph = Vr/sqrt(3);
I = Pr/(3*Vrph*PF);
I = I*(PF - j*sin(acos(PF)));
Ic1=j*(S/2)*Vrph;
Il=I+Ic1;
Vs=Vrph+(Il*(Z));
Ic2=j*(S/2)*Vs;
Is=Il+Ic2;
reg = (abs(Vs)-Vr)*100/Vr
Losses=3*(abs(Is)^2)*(R/2)+3*(abs(I)^2)*(R/2);
Efficiency= ((Pr)/((Pr)+Losses))*100
elseif(phase==1)
I=Pr/(Vr*PF);
I = I*(PF - j*sin(acos(PF)));
Ic1=j*(S/2)*Vr;
Il=I+Ic1;
Vs=Vr+Il*Z;
Ic2=j*(S/2)*Vs;
Is=Ic2+Il;
reg = (abs(Vs)-Vr)*100/Vr
Losses=(abs(Is)^2)*(R/2)+(abs(I)^2)*(R/2);
Efficiency= ((Pr)/((Pr)+Losses))*100
end
else
disp('Invalid Method')
end
