diff options
Diffstat (limited to '3872/CH8/EX8.7/EX8_7.sce')
| -rw-r--r-- | 3872/CH8/EX8.7/EX8_7.sce | 40 |
1 files changed, 40 insertions, 0 deletions
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)));
+
+
+
|