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
|
// 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)
|
