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
48
49
50
51
52
53
54
55
56
57
58
59
60
|
// Book - Power System: Analysis & Design 5th Edition
// Authors - J. Duncan Glover, Mulukutla S. Sharma, Thomas J. Overbye
// Chapter - 5 : Example 5.1
// Scilab Version 6.0.0 : OS - Windows
clc;
clear;
f = 60; // Frequency in Hz
N = 2; // Number of Conductors
V = 345; // Voltage in kV
L = 200; // Line length
S = 795000; // Size of the conductor
z = 0.032+(%i*0.35); // Impedance in Ohm/km
y = (%i)*(4.2*10^-6); // Admitance in S/km
Pr = 700; // Full load Power in MW
pf = 0.99; // Power factor
v = 95/100; // rated voltage
Z = z*L; // Total series impedance
Y = y*L; // Total shunt Admitance
A = 1 + ((Y*Z)/2); // Line Paramater A in per unit
D = A; // Line Paramater D in per unit
B = Z; // Line Paramater B in Ohm
C = Y*(1+(Y*Z)/4); // Line Paramater C in Siemens
VrLL = V*v; // Receiving end Line to Line Voltage in kVLL
VrLN = VrLL/sqrt(3); // Receiving end Line to Neutral Voltage in kVLN
theta = acos(pf);
Ir = (((Pr)*exp(%i*theta))/(sqrt(3)*(v*V)*(pf))); // Receiving end current in kA
VsLN = ((A*VrLN)+(B*Ir)); // Sending end Line to Neutral Voltage in kVLN
VsLL = abs(VsLN)*sqrt(3); // sending end Line to Line Voltage in kVLL
Is = ((C*VrLN)+(D*Ir)); // sending end current in kA
[r theta1] = polar(VsLN);
[r theta2] = polar(Is);
Ps = sqrt(3)*abs(VsLL)*abs(Is)*cos(theta1-theta2); // Power delivered to the sending end in MW
VrNL = abs(VsLL)/abs(A); // No load receiving end voltage in kVLL
PercentVR = ((abs(VrNL)-abs(VrLL))/abs(VrLL))*100; // Full load voltage in percent
J = N*0.9; // Approximate Current carrying capacity of 2 ACSR conductors in kA taken From table A.4
P = Ps-Pr; // Full load line losses in MW
PercentEFF = (Pr/Ps)*100; // Full load transmission efficiency in percent
printf('The magnitude of Transmission line parameter A in per unit is %0.4f and its angle is %0.3f degree', abs(A), atand(imag(A), real(A)));
printf('\nThe magnitude of Transmission line parameter B in Ohm is %0.2f and its angle is %0.2f degree', abs(B), atand(imag(B), real(B)));
printf('\nThe magnitude of Transmission line parameter C in Siemens is %0.2e and its angle is %0.2f degree', abs(C), atand(imag(C), real(C)));
printf('\nThe magnitude of Transmission line parameter D in per unit is %0.4f and its angle is %0.3f degree', abs(D), atand(imag(D), real(D)));
printf('\n\nSending end Line to Neutral Voltage in kVLN is : %0.1f and its angle is : %0.2f degree', abs(VsLN), atand(imag(VsLN), real(VsLN)));
printf('\nSending end Line to Line Voltage is (VsLL) = %0.1f kV', abs(VsLL));
printf('\nThe magnitude of sending end current in kA is (Is) : %0.3f and its angle is : %0.2f degree', abs(Is), atand(imag(Is), real(Is)));
printf('\nPower delivered to the sending end is (Ps) = %0.1f MW', Ps);
printf('\n\nNo load receiving end voltage is (VrNL) = %0.1f kVLL', VrNL);
printf('\nFull load voltage is (Percent VR) = %0.1f Percent', PercentVR);
printf('\n\nApproximate Current carrying capacity of 2 ACSR conductors is (J) = %0.1f kA', J);
printf('\nFull load line losses is (P) = %0.1f MW', P);
printf('\nFull load transmission efficiency is (Percent EFF) = %0.1f Percent', PercentEFF);
|
