Scribd
Upload
2 views

MATLAB cod1

The MATLAB code calculates the electrical design criteria for a transmission line connecting a hydro power plant to a substation in the Integrated Nepal Power System. It evaluates parameters such as ampacity, efficiency, voltage regulation, and corona inception voltage for different conductor configurations. The calculations are based on given variables including voltage, power, resistance, and geometric factors of the conductors.

Uploaded by

aavanishk
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
2 views

MATLAB cod1

The MATLAB code calculates the electrical design criteria for a transmission line connecting a hydro power plant to a substation in the Integrated Nepal Power System. It evaluates parameters such as ampacity, efficiency, voltage regulation, and corona inception voltage for different conductor configurations. The calculations are based on given variables including voltage, power, resistance, and geometric factors of the conductors.

Uploaded by

aavanishk
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

MATLAB code:

%Check the electrical design criterion of any transmission line of your preference
connecting a hydro power to a substation of Integrated Nepal Power System. Repeat
this calculation for “two more better conductors than one currently used (See Table
1)” and compare the electrical performance for all conductor cases. Use single
circuit configuration for voltage level up to 132kV and double circuit configuration
for voltage level from 220kV and above,

%the variables are


v = 133 X 10^3; % receiving voltage
pf = 0.95;
p = 30 X 10^6; % power transmitted
alpha = 0.004; % alpha
t_max = 40;
t_0 = 20; %
r_20 = 0.45650;
n = 1; % number of line circuits
line_length = 28;

i = p / (sqrt(3) * v * pf);

% Resistance at 40°C
r_40 = r_20 * (1 + alpha * (t_max - t_0));

% Total power loss and efficiency


total_loss = n * 3 * i^2 * r_40 * line_length;
efficiency = (1 - total_loss / p) * 100;

% Voltage regulation parameters


s = 7; % number of strands
d = 10.98;
s7 = 0.768; % GMR value
GMRL = s7 * d;
GMRC = d / 2;

%diagonal
da = 8.97; % diagonal distance between conductors
db = 8.728; % distance between two conductors
dc = 2.055; % height of the two conductors

GMD = (da * db * dc)^(1 / 3);


L = (2e-7) * log((GMD * 1000) / GMRL) * 1000; % inductance
C = 0.02411 / log10((GMD * 1000) / GMRC); % capacitance

% Impedance and admittance


w = 314.16; % angular frequency
Z = (r_40 + 1i * w * L) * line_length; % impedance
Y = 1i * (2 * w * C * 1e-6) * line_length; % admittance
A = 1 + (Y * Z) / 2;
B = Z;
C1 = Y * (1 + (Y * Z) / 4);
D = A;
Vr = v / sqrt(3); % receiving voltage
Ir = complex(62.3312, -20.4814); % receiving current
Vs = A * Vr + B * Ir; % sending voltage
Is = C1 * Vr + D * Ir; % sending current

Vnl = abs(Vs) / abs(A);


VR = ((Vnl - Vr) / Vr) * 100;

% Corona effect
m = 0.9;
ro = 0.9;
Vci = sqrt(3) * 21.21 * (GMRC / 10) * m * ro * log(GMD / (GMRC / 1000)); % corona
voltage

% Display results
fprintf('Ampacity is %.2f A\n', i);
fprintf('Efficiency is %.2f %%\n', efficiency);
fprintf('Voltage Regulation is %.2f %%\n', VR);
fprintf('Corona Inception Voltage is %.2f kV\n', Vci);

Result:

You might also like