diff options
Diffstat (limited to '3872/CH8')
| -rw-r--r-- | 3872/CH8/EX8.1/EX8_1.jpg | bin | 0 -> 23162 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.1/EX8_1.sce | 17 | ||||
| -rw-r--r-- | 3872/CH8/EX8.2/EX8_2.jpg | bin | 0 -> 21631 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.2/EX8_2.sce | 17 | ||||
| -rw-r--r-- | 3872/CH8/EX8.3/EX8_3.jpg | bin | 0 -> 34548 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.3/EX8_3.sce | 20 | ||||
| -rw-r--r-- | 3872/CH8/EX8.4/EX8_4.jpg | bin | 0 -> 28492 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.4/EX8_4.sce | 19 | ||||
| -rw-r--r-- | 3872/CH8/EX8.5/EX8_5.jpg | bin | 0 -> 21453 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.5/EX8_5.sce | 18 | ||||
| -rw-r--r-- | 3872/CH8/EX8.6/EX8_6.jpg | bin | 0 -> 27458 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.6/EX8_6.sce | 34 | ||||
| -rw-r--r-- | 3872/CH8/EX8.7/EX8_7.jpg | bin | 0 -> 17889 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.7/EX8_7.sce | 40 | ||||
| -rw-r--r-- | 3872/CH8/EX8.8/EX8_8.jpg | bin | 0 -> 18129 bytes | |||
| -rw-r--r-- | 3872/CH8/EX8.8/EX8_8.sce | 24 | ||||
| -rw-r--r-- | 3872/CH8/EX8.9/EX8_9.sce | 42 | ||||
| -rw-r--r-- | 3872/CH8/EX8.9/Ex8_9.JPG | bin | 0 -> 13483 bytes |
18 files changed, 231 insertions, 0 deletions
diff --git a/3872/CH8/EX8.1/EX8_1.jpg b/3872/CH8/EX8.1/EX8_1.jpg Binary files differnew file mode 100644 index 000000000..b30dc65ca --- /dev/null +++ b/3872/CH8/EX8.1/EX8_1.jpg diff --git a/3872/CH8/EX8.1/EX8_1.sce b/3872/CH8/EX8.1/EX8_1.sce new file mode 100644 index 000000000..07519b4cc --- /dev/null +++ b/3872/CH8/EX8.1/EX8_1.sce @@ -0,0 +1,17 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.1
+//Scilab Version - 6.0.0 ; OS - Windows
+clc;
+clear;
+Vp = [277; 277*(cos(-120*%pi/180)+%i*sin(-120*%pi/180)); 277*(cos(120*%pi/180)+%i*sin(120*%pi/180))]; //given column vector of phase voltage in volts
+function [Vp1]=phaseshift(x1,x2) //Function for shifting the phase
+ [r theta]=polar(x1);
+ Vp1=r*(cos(theta+x2*%pi/180)+%i*sin(theta+x2*%pi/180));
+endfunction
+V0 = 1*(Vp(1,1)+Vp(2,1)+Vp(3,1))/3; //zero sequence voltage in V
+V1 = 1*(Vp(1,1)+phaseshift(Vp(2,1),120)+phaseshift(Vp(3,1),240))/3; //positive sequence voltage in V
+V2 = 1*(Vp(1,1)+phaseshift(Vp(2,1),240)+phaseshift(Vp(3,1),120))/3; //negative sequence voltage in V
+printf('\nThe zero sequence voltage V0 = %f V',V0);
+printf('\nThe positive sequence voltage V1 = %f V',V1);
+printf('\nThe negative sequence voltage V2 = %f V',V2);
diff --git a/3872/CH8/EX8.2/EX8_2.jpg b/3872/CH8/EX8.2/EX8_2.jpg Binary files differnew file mode 100644 index 000000000..21c8bac3f --- /dev/null +++ b/3872/CH8/EX8.2/EX8_2.jpg diff --git a/3872/CH8/EX8.2/EX8_2.sce b/3872/CH8/EX8.2/EX8_2.sce new file mode 100644 index 000000000..8b3d83108 --- /dev/null +++ b/3872/CH8/EX8.2/EX8_2.sce @@ -0,0 +1,17 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.2
+//Scilab Version - 6.0.0 ; OS - Windows
+clc;
+clear;
+Ip = [10; 10*(cos(120*%pi/180)+%i*sin(120*%pi/180)); 10*(cos(-120*%pi/180)+%i*sin(-120*%pi/180))]; //given column vector of phase current in A
+function [Ip1]=phaseshift(x1,x2) //Function for shifting the phase
+ [r theta]=polar(x1);
+ Ip1=r*(cos(theta+x2*%pi/180)+%i*sin(theta+x2*%pi/180));
+endfunction
+I0 = 1*(Ip(1,1)+Ip(2,1)+Ip(3,1))/3; //zero sequence current in A
+I1 = 1*(Ip(1,1)+phaseshift(Ip(2,1),120)+phaseshift(Ip(3,1),240))/3; //positive sequence current in A
+I2 = (Ip(1,1)+phaseshift(Ip(2,1),240)+phaseshift(Ip(3,1),120))/3; //negative sequence current in A
+printf('\nThe zero sequence current V0 = %f A',I0);
+printf('\nThe positive sequence current V1 = %f A',I1);
+printf('\nThe negative sequence current V2 = %f A',I2);
diff --git a/3872/CH8/EX8.3/EX8_3.jpg b/3872/CH8/EX8.3/EX8_3.jpg Binary files differnew file mode 100644 index 000000000..80f456e4e --- /dev/null +++ b/3872/CH8/EX8.3/EX8_3.jpg diff --git a/3872/CH8/EX8.3/EX8_3.sce b/3872/CH8/EX8.3/EX8_3.sce new file mode 100644 index 000000000..89a919222 --- /dev/null +++ b/3872/CH8/EX8.3/EX8_3.sce @@ -0,0 +1,20 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.3
+//Scilab Version - 6.0.0 ; OS - Windows
+clc;
+clear;
+Ip = [10; 0; 10*(cos(120*%pi/180)+%i*sin(120*%pi/180))];; //given column vector of phase current in A
+function [Ip1]=phaseshift(x1,x2) //Function for shifting the phase
+ [r theta]=polar(x1);
+ Ip1=r*(cos(theta+x2*%pi/180)+%i*sin(theta+x2*%pi/180));
+endfunction
+
+I0 = (Ip(1,1)+Ip(2,1)+Ip(3,1))/3; //zero sequence current in A
+I1 = 1*(Ip(1,1)+(Ip(2,1)+phaseshift(Ip(3,1),240)))/3; //positive sequence current in A
+I2 = (Ip(1,1)+Ip(2,1)+phaseshift(Ip(3,1),120))/3; //negative sequence current in A
+In = (Ip(1,1)+Ip(2,1)+Ip(3,1)); //neutral current in A
+printf('\nThe magnitude of zero sequence current I0 in Ampere is %0.3f and its angle is %0.3f degree',abs(I0), atand(imag(I0), real(I0)));
+printf('\nThe magnitude of positive sequence current in Ampere is %0.3f and its angle is %0.3f degree ',abs(I1), atand(imag(I1), real(I1)));
+printf('\nThe magnitude of negative sequence current in Ampere is %0.3f and its angle is %0.3f degree',abs(I2), atand(imag(I2), real(I2)));
+printf('\nThe magnitude of neutral current in Ampere is %0.3f and its angle is %0.3f degree',abs(In), atand(imag(In), real(In)));
diff --git a/3872/CH8/EX8.4/EX8_4.jpg b/3872/CH8/EX8.4/EX8_4.jpg Binary files differnew file mode 100644 index 000000000..678d8b830 --- /dev/null +++ b/3872/CH8/EX8.4/EX8_4.jpg diff --git a/3872/CH8/EX8.4/EX8_4.sce b/3872/CH8/EX8.4/EX8_4.sce new file mode 100644 index 000000000..e906a446c --- /dev/null +++ b/3872/CH8/EX8.4/EX8_4.sce @@ -0,0 +1,19 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.4
+//Scilab Version - 6.0.0 ; OS - Windows
+clc;
+clear;
+Zy = (3+(%i*4)); //Y load impedance per phase
+Xn = 2; //inductive reactance in ohm per phase
+Xc = -%i*30; //capacitor bank reactance in ohm per phase
+Zn = %i*2 //neutral impedance in ohm per phase
+Zdel = Xc/3;
+
+Z0 = Zy+(3*Zn); //zero load sequence impedane in ohm
+Z1 = 1/(1/Zy+1/Zdel); //positive load sequence impedane in ohm
+Z2 =Z1; //negativa load sequence impedane in ohm
+printf('\nThe zero load sequence impedance Z0 is %0.4f + %0.4fi ohm',real(Z0), imag(Z0));
+printf('\nThe amplitude of positive load sequence impedance Z1 is %.4f ohm and its angle is %.4f degree ',abs(Z1), atand(imag(Z1), real(Z1)));
+printf('\nThe amplitude of negative load sequence impedance Z2 is %.4f ohm and its angle is %.4f degree ',abs(Z2), atand(imag(Z2), real(Z2)));
+
diff --git a/3872/CH8/EX8.5/EX8_5.jpg b/3872/CH8/EX8.5/EX8_5.jpg Binary files differnew file mode 100644 index 000000000..1585ee12c --- /dev/null +++ b/3872/CH8/EX8.5/EX8_5.jpg diff --git a/3872/CH8/EX8.5/EX8_5.sce b/3872/CH8/EX8.5/EX8_5.sce new file mode 100644 index 000000000..c2d59db7d --- /dev/null +++ b/3872/CH8/EX8.5/EX8_5.sce @@ -0,0 +1,18 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.5
+//Scilab Version - 6.0.0 ; OS - Windows
+
+clc;
+clear;
+
+Zn = %i*10; //generator neutral impedance in ohm
+Zgo = %i*1; //generator zero sequence impedance in ohm
+Zg1 = %i*15; //generator positive sequence impedance in ohm
+Zg2 = %i*3; //generator negative sequence impedance in ohm
+Zl1 = 0.087+(%i*0.99); //line impedace in ohm
+Zdel = 22.98+%i*(19.281); //impedance of the delta load in ohm
+V1=(415.69-(%i*240))/sqrt(3); //RMS line to neutral phase voltage of AC supply in Volts
+I1 = V1/(Zl1+((1/3)*Zdel)); //sequence component of line current in A
+
+printf('\nThe sequence component of the line current Ia is %.4f amperes and its angle is %.4f degree ',abs(I1), atand(imag(I1), real(I1)));
diff --git a/3872/CH8/EX8.6/EX8_6.jpg b/3872/CH8/EX8.6/EX8_6.jpg Binary files differnew file mode 100644 index 000000000..619fb7388 --- /dev/null +++ b/3872/CH8/EX8.6/EX8_6.jpg diff --git a/3872/CH8/EX8.6/EX8_6.sce b/3872/CH8/EX8.6/EX8_6.sce new file mode 100644 index 000000000..1a90e57b4 --- /dev/null +++ b/3872/CH8/EX8.6/EX8_6.sce @@ -0,0 +1,34 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.6
+//Scilab Version - 6.0.0 ; OS - Windows
+
+clc;
+clear;
+
+Vp = [277; 260*(cos(-120*%pi/180)+%i*sin(-120*%pi/180)); 295*(cos(115*%pi/180)+%i*sin(115*%pi/180))]; //given column vector of phase voltage in volts
+Zl1 = 0.087+%i*(0.99); //impedace of line 1 in ohm
+Zdel = 22.98+%i*(19.281); //impedance of the delta load in ohm
+Zl2 = 0.087+%i*(0.99); //impedance of line 2 in ohm
+function [Vp1]=phaseshift(x1,x2) //Function for shifting the phase
+ [r theta]=polar(x1);
+ Vp1=r*(cos(theta+x2*%pi/180)+%i*sin(theta+x2*%pi/180));
+endfunction
+
+V0 = (Vp(1,1)+Vp(2,1)+Vp(3,1))/3; //zero sequence voltage in V
+V1 = (Vp(1,1)+phaseshift(Vp(2,1),120)+phaseshift(Vp(3,1),240))/3; //positive sequence voltage in V
+V2 = (Vp(1,1)+phaseshift(Vp(2,1),240)+phaseshift(Vp(3,1),120))/3; //negative sequence voltage in V
+I0 = 0; //zero sequence current in A
+I1 = V1/(Zl1+(Zdel/3)); //positive sequence current in A
+I2 = V2/(Zl2+(Zdel/3)); //negative sequence current in A
+Ia = I0+I1+I2; //zero source current in A
+Ib = I0+phaseshift(I1,240)+phaseshift(I2,120); //positive source current in A
+Ic = I0+phaseshift(I1,120)+phaseshift(I2,240); //negative source current in A
+printf('The zero source current Ia is %.4f amperes and its angle is %.4f degree ',abs(Ia), atand(imag(Ia), real(Ia)));
+printf('\nThe positive source current Ib is %.4f amperes and its angle is %.4f degree ',abs(Ib), atand(imag(Ib), real(Ib))+360);
+printf('\nThe negative source current Ic is %.4f amperes and its angle is %.4f degree ',abs(Ic), atand(imag(Ic), real(Ic)));
+
+
+
+
+
diff --git a/3872/CH8/EX8.7/EX8_7.jpg b/3872/CH8/EX8.7/EX8_7.jpg Binary files differnew file mode 100644 index 000000000..82a8c2e7c --- /dev/null +++ b/3872/CH8/EX8.7/EX8_7.jpg diff --git a/3872/CH8/EX8.7/EX8_7.sce b/3872/CH8/EX8.7/EX8_7.sce new file mode 100644 index 000000000..81b1b21de --- /dev/null +++ b/3872/CH8/EX8.7/EX8_7.sce @@ -0,0 +1,40 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.7
+//Scilab Version - 6.0.0 ; OS - Windows
+clc;
+clear;
+Q = 75; //rated power in kVA
+Vprtr = 480; //primary side voltage of transformer in volts
+Vsectr = 208; //secondary side voltage of transformer in volts
+Xeq = 0.10; //leakage reactance in per unit
+Sbase = Q/3; //base quantity of rated power in single phase in kVA
+VbaseHLN = Vprtr/(sqrt(3)); //base quantity of primary side voltage of transformer in volts
+VbaseXLN = Vsectr/(sqrt(3)); //base quantity of secondary side voltage of transformer in volts
+ZbaseX = 0.5770; //base quantity of impedance in ohm
+Vp = [277; 260*(cos(-120*%pi/180)+%i*sin(-120*%pi/180)); 295*(cos(115*%pi/180)+%i*sin(115*%pi/180))]; //given column vector of phase voltage in volts
+function [Vp1]=phaseshift(x1,x2)
+ [r theta]=polar(x1);
+ Vp1=r*(cos(theta+x2*%pi/180)+%i*sin(theta+x2*%pi/180));
+endfunction
+
+V0 = (Vp(1,1)+Vp(2,1)+Vp(3,1))/3; //zero sequence voltage in V
+V1 = (Vp(1,1)+phaseshift(Vp(2,1),120)+phaseshift(Vp(3,1),240))/3; //positive sequence voltage in V
+V2 = (Vp(1,1)+phaseshift(Vp(2,1),240)+phaseshift(Vp(3,1),120))/3; //negative sequene voltage in v
+V0 = V0/VbaseHLN; //zero sequence voltage in per unit
+V1 = V1/VbaseHLN; //positive sequence voltage in per unit
+V2 = V2/VbaseHLN; //negative sequene voltage in per unit
+Zline0 = 0.087+%i*(0.99); //line impedance in ohm
+Zload1 = 22.98+%i*(19.281); //load impedance in ohm
+Zline0 = Zline0/ZbaseX; //line impedance in per unit
+Zload1 = Zload1/(3*ZbaseX); //line impedance in per unit
+I0 = 0; //zero sequence component of source current in per unit
+I1 = V1/((%i*Xeq)+Zline0+Zload1); //positive sequence component of source current in per unit
+I2 = V2/((%i*Xeq)+Zline0+Zload1); //negative sequence component of source current in per unit
+Ia = I0+I1+I2; //phase 'a' source current in per unit
+IbaseH=(Q*10^3)/(Vprtr*sqrt(3)); //base current in A
+Ia = Ia*IbaseH; //phase 'a' source current in A
+printf('The magnitude of phase a source current Ia is %.4f Ampere and its angle is %.4f degree',abs(Ia),atand(imag(Ia),real(Ia)));
+
+
+
diff --git a/3872/CH8/EX8.8/EX8_8.jpg b/3872/CH8/EX8.8/EX8_8.jpg Binary files differnew file mode 100644 index 000000000..7a682e3f6 --- /dev/null +++ b/3872/CH8/EX8.8/EX8_8.jpg diff --git a/3872/CH8/EX8.8/EX8_8.sce b/3872/CH8/EX8.8/EX8_8.sce new file mode 100644 index 000000000..e11b07266 --- /dev/null +++ b/3872/CH8/EX8.8/EX8_8.sce @@ -0,0 +1,24 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.8
+//Scilab Version - 6.0.0 ; OS - Windows
+
+clc;
+clear;
+
+Q = 900; //rated power in MVA
+Vg = 13.8; //generator voltage in kV
+Vt = 345; //transmission line voltage in kV
+Vd = 34.5; //distribution line voltage in kV
+V1 = 13.8; //voltage at the winding X in kV
+V2 = 199.2; //voltage at the winding H in kV
+V3 = 19.92; //voltage at the winding M in kV
+Zn = %i*0.10; //neutral impedance in ohm
+VbaseX = 13.8; //rated line to line voltage of terminal X in kV
+VbaseM = sqrt(3)*V3; //rated line to line voltage of terminal M in kV
+ZbaseM = (Vd^2)/Q; //base impedance of medium line voltage in ohm
+Zn = Zn/ZbaseM; //neutral impedance in per unit
+
+printf('\n The base impedance of medium voltage terminal ZbaseM is %f ohm',ZbaseM);
+printf('\n The per unit neutral impedance is Zn is i%0.4f per unit', imag(Zn));
+
diff --git a/3872/CH8/EX8.9/EX8_9.sce b/3872/CH8/EX8.9/EX8_9.sce new file mode 100644 index 000000000..7f21a866a --- /dev/null +++ b/3872/CH8/EX8.9/EX8_9.sce @@ -0,0 +1,42 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.8
+//Scilab Version - 6.0.0 ; OS - Windows
+
+clc;
+clear;
+
+Vp = [277; 260*(cos(-120*%pi/180)+%i*sin(-120*%pi/180)); 295*(cos(115*%pi/180)+%i*sin(115*%pi/180))]; //given column vector of phase voltage in volts
+Zl1 = 0.087+%i*(0.99); //impedace of line 1 in ohm
+Zdel = 22.98+%i*(19.281); //impedance of the delta load in ohm
+Zl2 = 0.087+%i*(0.99); //impedance of line 2 in ohm
+function [Vp1]=phaseshift(x1,x2)
+ [r theta]=polar(x1);
+ Vp1=r*(cos(theta+x2*%pi/180)+%i*sin(theta+x2*%pi/180));
+endfunction
+
+V0 = (Vp(1,1)+Vp(2,1)+Vp(3,1))/3; //zero sequence voltage in V
+V1 = (Vp(1,1)+phaseshift(Vp(2,1),120)+phaseshift(Vp(3,1),240))/3; //positive sequence voltage in V
+V2 = (Vp(1,1)+phaseshift(Vp(2,1),240)+phaseshift(Vp(3,1),120))/3; //negative sequence voltage in V
+I0 = 0; //zero sequence current in A
+I1 = V1/(Zl1+(Zdel/3)); //positive sequence current in A
+I2 = V2/(Zl2+(Zdel/3)); //negative sequence current in A
+Ia = I0+I1+I2; //zero source current in A
+Ib = I0+phaseshift(I1,240)+phaseshift(I2,120); //positive source current in A
+Ic = I0+phaseshift(I1,120)+phaseshift(I2,240); //negative source current in A
+Sp = (Vp(1,1)*(conj(Ia)))+(Vp(2,1)*(conj(Ib)))+(Vp(3,1)*(conj(Ic))); //total complex power delivered to load in VA
+Ss = (V0*conj(I0))+(V1*conj(I1))+(V2*conj(I2)); //total complex power delivered to the sequence networks in VA
+SS = 3*Ss;
+printf('\n 3Ss = %0.2f , Sp = %0.2f ',abs(SS), abs(Sp));
+if (ceil(real(SS))==ceil(real(Sp))) then
+ printf('\n Sp is equal to 3Ss');
+else
+ printf('\n Sp is not equal to 3Ss');
+end
+
+
+
+
+
+
+
diff --git a/3872/CH8/EX8.9/Ex8_9.JPG b/3872/CH8/EX8.9/Ex8_9.JPG Binary files differnew file mode 100644 index 000000000..b3b58f960 --- /dev/null +++ b/3872/CH8/EX8.9/Ex8_9.JPG |
