diff options
Diffstat (limited to '3872/CH5/EX5.1/Ex5_1.sce')
| -rw-r--r-- | 3872/CH5/EX5.1/Ex5_1.sce | 60 |
1 files changed, 60 insertions, 0 deletions
diff --git a/3872/CH5/EX5.1/Ex5_1.sce b/3872/CH5/EX5.1/Ex5_1.sce new file mode 100644 index 000000000..78ddea092 --- /dev/null +++ b/3872/CH5/EX5.1/Ex5_1.sce @@ -0,0 +1,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);
|