1 files changed, 42 insertions, 0 deletions
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
+
+
+
+
+
+
+