Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Parameters Calculation For Transmission Line

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 3
At a glance
Powered by AI
The document discusses the calculation of parameters like resistance, reactance and susceptance of a transmission line using MATLAB code. It also presents an alternative method to model the transmission line and calculate its parameters.

The aim of the experiment is to calculate the parameters of a transmission line like resistance, reactance and susceptance.

The parameters calculated for the transmission line using MATLAB code are: resistance (R), inductance (L), capacitance (C), impedance (Z), admittance (Y) and other derived parameters.

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

You might also like