0% found this document useful (0 votes)
17 views3 pages

Code

This MATLAB program calculates line parameters and performance for a 3-phase transmission line by inputting length, conductor properties, and bundled/unbundled configuration. It then models the line as short, medium, or long and outputs voltage, current, power, efficiency, and other metrics.

Uploaded by

tarunsharma.amu
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
17 views3 pages

Code

This MATLAB program calculates line parameters and performance for a 3-phase transmission line by inputting length, conductor properties, and bundled/unbundled configuration. It then models the line as short, medium, or long and outputs voltage, current, power, efficiency, and other metrics.

Uploaded by

tarunsharma.amu
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

% MATLAB Program to calculate line parameters and performance of the line using 3

phase transmission line system


clc
fprintf('----Enter the values----\n\n');
Length =input('Enter the length of the Transmission Line (in km): ');
f = 50; % in Hz
R1 = input('Enter the radius of Conductor (in m): ');
GMR = 0.7788*R1;
Dab = input('distance between conductors 1 and 2 (in m) ');
Dbc = input('distance between conductors 2 and 3 (in m) ');
Dca = input('distance between conductors 3 and 1 (in m) ');
GMD = (Dab*Dbc*Dca)^(1/3);
r = input('Enter the resistance of Transmission line (in Ohm/km): ');
bundling = input("Is the line bundled? '0' for NO '1' for YES : ");
fprintf("Bundling option: %d\n", bundling);

if bundling == 0
fprintf("GMR: %.6f\n", GMR);
GMR = 0.7788*Radius;
GMD = (Dab*Dbc*Dca)^(1/3);
l = (2e-4)*log(GMD/GMR); % in H/km
c = (2*pi*8.85e-9)/log(GMD/GMR); % in F/km and line to neutral
fprintf("l: %.6f H/km\n", l);
fprintf("c: %.6f F/km (line to neutral)\n", c);

elseif bundling == 1
Conductor = input("Enter the number of conductors (2 or 3) : ");
fprintf("Number of conductors: %d\n", Conductor);

conductor_spacing = input("Enter conductor spacing (in m) : ");


fprintf("Conductor spacing: %.2f m\n", conductor_spacing);

if Conductor == 2
GMR = sqrt(conductor_spacing * Radius * 0.7788);
GMD = (conductor_spacing);
fprintf("GMR: %.6f\n", GMR);
l = (2e-4)*log(GMD/GMR); % in H/km
c = (2*pi*8.85e-9)/log(GMD/GMR); % in F/km and line to neutral
fprintf("l: %.6f H/km\n", l);
fprintf("c: %.6f F/km (line to neutral)\n", c);
elseif Conductor == 3
GMR = sqrt(conductor_spacing^2 * Radius * 0.7788);
GMD = (conductor_spacing*conductor_spacing*conductor_spacing)^(1/3);
fprintf("GMR: %.6f\n", GMR);
l = (2e-4)*log(GMD/GMR); % in H/km
c = (2*pi*8.85e-9)/log(GMD/GMR); % in F/km and line to neutral
fprintf("l: %.6f H/km\n", l);
fprintf("c: %.6f F/km (line to neutral)\n", c);
else
disp("Invalid number of conductors. Please enter 2 or 3.");
return;
end
else
disp("Invalid bundling option. Please enter '0' or '1'.");
return;
end
l = (2e-4)*log(GMD/GMR); % in H/km
c = (2*pi*8.85e-9)/log(GMD/R1); % in F/km and line to neutral
L = l*Length;
R = r*Length;
C = c*Length;
% per phase series impedance and Shunt Capacitance
Xl = 2*pi*f*L; % in Ohm
Xc = 1/(2*pi*f*C); % in Ohm
Z = R+1i*Xl;
Y = 1i/Xc;
% Receiving end parameters
Pr1 = input('Receiving end Power (in MW): ');
Pr = Pr1*(1e+6);
Vr1 = input('Receiving end Voltage (in kV): ');
Vr = (Vr1*(1e+3))*(1+1i*0);
PF = input('Power Factor: ');
sint = (1-PF^2)^(1/2);
Ir1 = Pr/(3*Vr*PF);
Ir = Ir1*(PF-1i*sint);
%String Efficiency
m = input('Enter the value of m (i.e. mutual Capacitance/capacitance to ground)
:');
v1 = Vr1/((1+m)*(3+m));
v2 = v1*(1+m);
v3 = v2+(v1+v2)*m;
se = (Vr1*100)/(3*v3);
% Transmission Line model
if (Length<=80)
fprintf('\n------Short Transmission Line Model-----\n');
Is = Ir;
Vs = Vr+Ir*Z; % per phase
Vs1 = abs(Vs);
Vsl = sqrt(3)*Vs;
Is1 = abs(Is);
Vreg = ((Vs1-Vr)/Vr)*100;
Loss = 3*(R)*(Is1^2);
Ps = Pr+Loss;
PFs = Ps/(3*Vs1*Is1);
Efficiency = (Pr/Ps)*100;
elseif (Length>=80 && Length<=200)
fprintf('\n------Medium Transmission Line Model-----\n');
Model=input('Enter the desirable Model:\n For Nominal-T press 1 \n For Nominal-Pi
press 2');
if (Model==1)
V1 = Vr+Ir*(Z/2);
Ic = Y*V1;
Is = Ir+Ic;
Vs = V1+Is*(Z/2);
Vsl = sqrt(3)*Vs;
Vs1 = abs(Vs);
Is1 = abs(Is);
Vro = Vs*(-1i*Xc)/(Z/2-1i*Xc);
Vreg = ((abs(Vro)-Vr)/Vr)*100;
Loss = 3*(R/2)*(Is1^2+Ir1^2);
Ps = Pr+Loss;
PFs = Ps/(3*Vs1*Is1);
Efficiency = (Pr/Ps)*100;
elseif (Model==2)
Icr = (Y/2)*Vr;
Il = Ir+Icr;
Vs = Vr+Il*Z;
Ics = (Y/2)*Vs;
Is = Il+Ics;
Vsl = sqrt(3)*Vs;
Vs1 = abs(Vs);
Is1 = abs(Is);
Vro = Vs*(-1i*(Xc/2))/(Z/-1i*(Xc/2));
Vreg = ((abs(Vro)-Vr)/Vr)*100;
Loss = 3*R*(abs(Il)^2);
Ps = Pr+Loss;
PFs = Ps/(3*Vs1*Is1);
Efficiency = (Pr/Ps)*100;
end
else
fprintf('\n------Long Transmission Line Model-----\n');
Gamma = sqrt(Y*Z);
Zc = sqrt(Y/Z);
A = cosh(Gamma);
B = sinh(Gamma)/Zc;
C = Zc*sinh(Gamma);
D = A;
Vs = A*Vr+B*Ir;
Is = C*Vr+D*Ir;
Vsl = sqrt(3)*Vs;
Vs1 = abs(Vs);
Is1 = abs(Is);
Vreg = (((Vs1/abs(A))-Vr)/Vr)*100;
avs = angle(Vs)*(180/pi);
ais = angle(Is)*(180/pi);
PFs = cos(avs-ais);
Ps = 3*abs(Vs)*abs(Is)*PFs;
Efficiency = (Pr/Ps)*100;
end
fprintf('\n-------Transmission line Parameters-------\n\n');
fprintf('\nResistance of transmission line is %d Ohm\n',R);
fprintf('Inductance of transmission line is %d H\n',L);
fprintf('Capacitance of transmission line is %d F\n',C);
fprintf('Per phase series Impedance is %d Ohm \n',Z);
fprintf('Shunt capacitance is %d Ohm \n\n',Xc);
fprintf('\n-------Performance of transmission line-------\n');
fprintf('Sending end current: %d Ampere\n',Is);
fprintf('Sending end voltage: %d Volt\n',Vsl);
fprintf('Sending end power factor: %d \n',PFs);
fprintf('Sending end power: %d Watt\n',Ps);
fprintf('Efficiency: %d \n',Efficiency);
fprintf('Voltage Regulation: %d \n',Vreg);
fprintf('String Efficiency for 3 insulator disc is: %d \n',se);

You might also like