blob: 960dd8d183e5734d25c5964a23dd81ab71a0cee7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
// Variable Declaration
R = 0.16 //Resistance(ohm)
L = 1.26*10**(-3) //Inductance(H)
C = 8.77*10**(-9) //Capacitance(F)
l = 200.0 //Length of line(km)
P = 50.0 //Power(MVA)
pf = 0.8 //Lagging power factor
V_r = 132000.0 //Receiving end voltage(V)
f = 50.0 //Frequency(Hz)
// Calculation Section
w = 2 * %pi * f
z = complex(R, w*L) //Series impedance per phase per km(ohm)
y = complex(0, w*C) //Shunt admittance per phase per km(mho)
g = (y*z)**(0.5) //propagation constant(/km)
gl = g * l
Z_c = (z/y)**(0.5) //Surge impedance(ohm)
cosh_gl = cosh(gl)
sinh_gl = sinh(gl)
A = cosh_gl
B = Z_c * sinh_gl
C = (sinh_gl/Z_c)
D = cosh_gl
fi = acos(pf) //Power factor angle(radians)
V_R = V_r/(3**0.5) //Receiving end voltage(V)
I_R = (P*10**6/((3**0.5)*V_r))*(pf - complex(0,sin(fi)))//Receiving end current(A)
V_S = (A*V_R + B*I_R) //Sending end voltage(V/phase)
V_S_L = V_S * (3**0.5)*10**-3 //Sending end line voltage(kV)
I_S = C*V_R + D*I_R //Sending end current(A)
pf_S = cos((phasemag(I_S)*%pi/180) - (phasemag(V_S)*%pi/180)) //Sending end power factor
P_S = abs(V_S*I_S)*pf_S*10**-6 //Sending end power/phase(MW)
P_R = (P/3)*pf //Receiving end power/phase(MW)
P_L = 3*(P_S - P_R) //Total line loss(MW)
// Result Section
printf('Sending end voltage , V_S = %.2f∠%.2f° kV/phase' ,abs(V_S*10**-3),phasemag(V_S))
printf('Sending end line voltage = %.2f kV' ,abs(V_S_L))
printf('Sending end current , I_S = %.2f∠%.2f° A' ,abs(I_S),phasemag(I_S))
printf('Sending end power factor = %.2f lagging' ,pf_S)
printf('Total transmission line loss = %.3f MW' ,P_L)
printf('NOTE : Answers are slightly different because of rounding error.')
|
