diff options
Diffstat (limited to '3872/CH5/EX5.2/Ex5_2.sce')
| -rw-r--r-- | 3872/CH5/EX5.2/Ex5_2.sce | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/3872/CH5/EX5.2/Ex5_2.sce b/3872/CH5/EX5.2/Ex5_2.sce new file mode 100644 index 000000000..c1ad005a7 --- /dev/null +++ b/3872/CH5/EX5.2/Ex5_2.sce @@ -0,0 +1,32 @@ +// Book - Power System: Analysis & Design 5th Edition
+// Authors - J. Duncan Glover, Mulukutla S. Sharma, Thomas J. Overbye
+// Chapter - 5 : Example 5.2
+// Scilab Version 6.0.0 : OS - Windows
+
+clc;
+clear;
+
+V = 765; // Line voltage in kV
+f = 60; // frequency in Hz
+L = 300; // line length in km
+z = 0.0165+(%i*0.3306); // Positive sequence impedance in Ohm/km
+y = %i*4.674e-6; // Positive sequence admitance in S/km
+Zc = sqrt(z/y); // Characteristic impedance in Ohm
+GammaL = sqrt(z*y)*L; // Propagation constant in per unit
+eGammaL = exp(0.00930)*exp(%i*0.3730);
+eNegGammaL = exp(-0.00930)*exp(-%i*0.3730);
+coshGammaL = (eGammaL+eNegGammaL)/2; // Hyperbolic function
+sinhGammaL = (eGammaL-eNegGammaL)/2; // Hyperbolic function
+A = cosh(GammaL); // line parameter in per unit
+D = A; // line parameter in per unit
+B = Zc*sinh(GammaL); // line parameter in Ohm
+C = (1/Zc)*sinh(GammaL); // Line parameter in S
+Bnominalpi = z*L; // The B parameter for the nominal pi circuit in Ohm
+Bdiff=100-(abs(B)/abs(Bnominalpi))*100; //The difference in B parameter in percentage
+
+printf('\n\The 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 line parameter B in Ohm is %0.1f and its angle is : %0.1f degree', abs(B), atand(imag(B), real(B)));
+printf('\nThe line parameter C in Siemens is %0.2e and its angle is : %0.2f degree', abs(C), atand(imag(C), real(C)));
+printf('\n\The line parameter D in per unit is %0.4f and its angle is : %0.3f degree', abs(A), atand(imag(A), real(A)));
+printf('\nThe B parameter for the nominal pi circuit in Ohm is (Bnominalpi) : %0.1f and its angle is : %0.2f degree', abs(Bnominalpi), atand(imag(Bnominalpi), real(Bnominalpi)));
+printf('\nThe difference in B parameter for the nominal pi circuit is %d percentage',Bdiff)
|