initial commit / add all books

8 files changed, 123 insertions, 0 deletions

diff --git a/3872/CH14/EX14.1/EX14_1.jpg b/3872/CH14/EX14.1/EX14_1.jpg
new file mode 100644
index 000000000..5a2d83caa
--- /dev/null
+++ b/3872/CH14/EX14.1/EX14_1.jpg
diff --git a/3872/CH14/EX14.1/EX14_1.sce b/3872/CH14/EX14.1/EX14_1.sce
new file mode 100644
index 000000000..d1a189498
--- /dev/null
+++ b/3872/CH14/EX14.1/EX14_1.sce
@@ -0,0 +1,28 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 14 ; Example 14.1

+//Scilab Version - 6.0.0 ; OS - Windows

+

+

+clc;

+clear;

+Vdelpri=230;                                //The rated RMS line voltage of the primary winding in kV

+Vwyesec=34.5;                               //The rated RMS line voltage of the secondary winding in kV

+MVAoa=75;                                   //The MVA OA rating of the transformer

+MVAfa=100;                                  //The MVA FA rating of the transformer

+MVAfoa=133;                                 //The MVA FOA rating of the transformer

+Zpu=0.07;                                   //The percentage impedance of the transformer in terms of MVA OA ratings

+MVAbase=100;                                //The MVA base in MVA

+kVbase=34.5;                                //The KV base in kV

+Ioa=(MVAoa/3)/(Vwyesec/sqrt(3));            //The OA transformer current in the low voltage side in kA

+Ifa=(MVAfa/3)/(Vwyesec/sqrt(3));            //The FA transformer current in the low voltage side in kA

+Ifoa=(MVAfoa/3)/(Vwyesec/sqrt(3));          //The FOA transformer current in the low voltage side in kA

+Zbasepu=Zpu*MVAbase/MVAoa;                  //The per unit impedance of the system in ohm pu

+Isc3ph=(1/Zpu)*Ioa;                         //The transformer current during three phase bolted faultin kA

+printf('\nThe OA transformer current in the low voltage side is %f kA',Ioa);

+printf('\nThe FA transformer current in the low voltage side is %f kA\n',Ifa);

+printf('The FOA transformer current in the low voltage side is %f kA\n',Ifoa);

+printf('The per unit impedance of the system is %f pu\n',Zbasepu);

+printf('The transformer current during three phase bolted fault is %f kA',Isc3ph);

+

+//The answer in the book is not correct. Eg.75/(sqrt(3)*34.5)---Actual result is 1.255, but it is given as 7.372 in the book.

diff --git a/3872/CH14/EX14.2/EX14_2.jpg b/3872/CH14/EX14.2/EX14_2.jpg
new file mode 100644
index 000000000..af7df6527
--- /dev/null
+++ b/3872/CH14/EX14.2/EX14_2.jpg
diff --git a/3872/CH14/EX14.2/EX14_2.sce b/3872/CH14/EX14.2/EX14_2.sce
new file mode 100644
index 000000000..df7fb91dd
--- /dev/null
+++ b/3872/CH14/EX14.2/EX14_2.sce
@@ -0,0 +1,20 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 14 ; Example 14.2

+//Scilab Version - 6.0.0 ; OS - Windows

+

+

+clc;

+clear;

+MVAtr1=40;                                                 //MVA FOA rating of transformer 1

+MVAtr2=40;                                                 //MVA FOA rating of transformer 2

+normal=1.28;                                               //Factor for normal summer operation

+emergency2hr=1.70;                                         //Factor for two hour emergency operation

+emergency30day=1.55;                                       //Factor for thirty days emergency operation

+unequalloadingfactor=0.95;                                 //Factor to account for unequal transformer loading

+MVAstation=normal*(MVAtr1+MVAtr2)*unequalloadingfactor;    //MVA rating of thr station

+MVAstationemergency2hr=emergency2hr*MVAtr1;                //MVA rating of a single transformer for two hour emergency

+MVAstationemergency30day=emergency30day*MVAtr1;            //MVA rating of a single transformer for thirty days emergency

+printf('\nThe summer normal rating of the station is %f MVA',MVAstation);

+printf('\nThe emergency rating of the single transformer for two hours is %f MVA',MVAstationemergency2hr);

+printf('\nThe emergency rating of the single transformer for thirty days is %f MVA',MVAstationemergency30day)

diff --git a/3872/CH14/EX14.3/EX14_3.jpg b/3872/CH14/EX14.3/EX14_3.jpg
new file mode 100644
index 000000000..08973c988
--- /dev/null
+++ b/3872/CH14/EX14.3/EX14_3.jpg
diff --git a/3872/CH14/EX14.3/EX14_3.sce b/3872/CH14/EX14.3/EX14_3.sce
new file mode 100644
index 000000000..44b090683
--- /dev/null
+++ b/3872/CH14/EX14.3/EX14_3.sce
@@ -0,0 +1,58 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 14 ; Example 14.3

+//Scilab Version - 6.0.0 ; OS - Windows

+

+

+clc;

+clear;

+kV=13.8;                                          //The sending end line voltage in kVolts

+Vsln=1.05*kV/sqrt(3);                             //The sending end voltage with 5% above rated in kVolts

+Rload=20;                                         //The Wye connected load resistance in Ohm

+Xload=40*%i;                                      //The Wye connected load inductive reactance in Ohm

+Xc=-40*%i;                                        //The Wye connected capacitive reactance in Ohm

+Rline=3;                                          //The line resistance in Ohm

+Xline=6*%i;                                       //The line inductive reactance in Ohm

+Ztot1=Rline+Xline+(Rload*Xload/(Rload+Xload));    //The total impedance seen by source without capacitance in Ohm

+Iline1=Vsln/Ztot1;                                //The line current without shunt capacitor in kA

+Vdrop1=(Rline+Xline)*Iline1;                      //The voltage drop across the line without shunt capacitor in KVolts

+Vload1=Vsln-Vdrop1;                               //The voltage drop across the load without shunt capacitor in KVolts

+Pload1=3*abs(Vload1)^2/Rload;                     //The real power delivered to the load without shunt capacitor in MW

+Qload1=3*abs(Vload1)^2/abs(Xload);                //The reactive power delivered to the load without shunt capacitor in MVAR

+pf1=cos((atan(Qload1/Pload1)));                   //The power factor of the load without shunt capacitor

+Pline1=3*abs(Iline1)^2*Rline;                     //The real power loss in the line without shunt capacitor in MW

+Qline1=3*abs(Iline1)^2*abs(Xline);                //The reactive power loss in the line without shunt capacitor in MVAR

+Psource1=Pload1+Pline1;                           //The real power delivered by the source without shunt capacitor in MW

+Qsource1=Qload1+Qline1;                           //The reactive power delivered by the source without shunt capacitor in MVAR

+Ssource1=sqrt(Psource1^2+Qsource1^2);             //The apparent power delivered by the source without shunt capacitor in MVA

+Ztot2=Rline+Xline+(1/(1/Rload+1/Xload+1/Xc));     //The total impedance seen by source with capacitance in Ohm

+Iline2=Vsln/Ztot2;                                //The line current with shunt capacitor in kA

+Vdrop2=(Rline+Xline)*Iline2;                      //The voltage drop across the line with shunt capacitor in KVolts

+Vload2=Vsln-Vdrop2;                               //The voltage drop across the load with shunt capacitor in KVolts

+Pload2=3*abs(Vload2)^2/Rload;                     //The real power delivered to the load with shunt capacitor in MW

+Qload2=3*abs(Vload2)^2/abs(Xload);                //The reactive power delivered to the load with shunt capacitor in MVAR

+pf2=cos((atan(Qload2/Pload2)));                   //The power factor of the load with shunt capacitor

+Pline2=3*abs(Iline2)^2*Rline;                     //The real power loss in the line with shunt capacitor in MW

+Qline2=3*abs(Iline2)^2*abs(Xline);                //The reactive power loss in the line with shunt capacitor in MVAR

+Qc=3*abs(Vload2)^2/abs(Xc);                       //The reactive power delivered by the shunt capacitor inb MVAR

+Psource2=Pload2+Pline2;                           //The real power delivered by the source with shunt capacitor in MW

+Qsource2=Qload2+Qline2-Qc;                        //The reactive power delivered by the source with shunt capacitor in MVAR

+Ssource2=sqrt(Psource2^2+Qsource2^2);             //The apparent power delivered by the source with shunt capacitor in MVA

+printf('a. Without Capacitor');

+printf('\nThe magnitude of line current is %f kA and %f degree',abs(Iline1),atand(imag(Iline1)/real(Iline1)));

+printf('\nThe magnitude of voltage drop in the line is %f kV and %f degree',abs(Vdrop1),atand(imag(Vdrop1)/real(Vdrop1)));

+printf('\nThe magnitude of voltage drop in the load is %f kV and %f degree',abs(Vload1),atand(imag(Vload1)/real(Vload1)));

+printf('\nThe real and reactive power delivered to the three phase load is %f MW and %f MVAR',Pload1,Qload1);

+printf('\nThe load power factor is %f lagging',pf1);

+printf('\nThe real and reactive power losses in the line is %f MW and %f MVAR',Pline1,Qline1);

+printf('\nThe real power, reactive power and Apparent power delivered by the source is %f MW , %f MVAR and %f MVA',Psource1,Qsource1,Ssource1);

+printf('\n\n\nb. With Capacitor');

+printf('\nThe magnitude of line current is %f kA and %f degree',abs(Iline2),atand(imag(Iline2)/real(Iline2)));

+printf('\nThe magnitude of voltage drop in the line is %f kV and %f degree',abs(Vdrop2),atand(imag(Vdrop2)/real(Vdrop2)));

+printf('\nThe magnitude of voltage drop in the load is %f kV and %f degree',abs(Vload2),atand(imag(Vload2)/real(Vload2)));

+printf('\nThe real and reactive power delivered to the three phase load is %f MW and %f MVAR',Pload2,Qload2);

+printf('\nThe load power factor is %f lagging',pf2);

+printf('\nThe real and reactive power losses in the line is %f MW and %f MVAR',Pline2,Qline2);

+printf('\nThe reactive power delivered by the shunt capacitor bank is %f MVAR',Qc);

+printf('\nThe real power, reactive power and Apparent power delivered by the source is %f MW , %f MVAR and %f MVA',Psource2,Qsource2,Ssource2);

+                                                  //The third part of this question cannot be executed in SCILAB because of its theoritical nature

diff --git a/3872/CH14/EX14.4/EX14_4.jpg b/3872/CH14/EX14.4/EX14_4.jpg
new file mode 100644
index 000000000..caad930cb
--- /dev/null
+++ b/3872/CH14/EX14.4/EX14_4.jpg
diff --git a/3872/CH14/EX14.4/EX14_4.sce b/3872/CH14/EX14.4/EX14_4.sce
new file mode 100644
index 000000000..9a302c3eb
--- /dev/null
+++ b/3872/CH14/EX14.4/EX14_4.sce
@@ -0,0 +1,17 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 14 ; Example 14.1

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+time_interruptions=[8.17 200;71.3 600; 30.3 25; 267.2 90; 120 700; 10 1500; 40 100];                                                                                                                                               //Number of custumers interrupted for a time duration in minutes

+total_customers=2000;                //Total number of customers

+SAIFI=(time_interruptions(1,2)+time_interruptions(2,2)+time_interruptions(3,2)+time_interruptions(4,2)+time_interruptions(5,2)+time_interruptions(6,2)+time_interruptions(7,2))/total_customers;                                   //System average interruption frequency index

+SAIDI=(time_interruptions(1,2)*time_interruptions(1,1)+time_interruptions(2,2)*time_interruptions(2,1)+time_interruptions(3,2)*time_interruptions(3,1)+time_interruptions(4,2)*time_interruptions(4,1)+time_interruptions(5,2)*time_interruptions(5,1)+time_interruptions(6,2)*time_interruptions(6,1)+time_interruptions(7,2)*time_interruptions(7,1))/total_customers;                                                         //System average interruption duration index

+CAIDI=SAIDI/SAIFI;                   //Customer average interruption duration index

+ASAI=(365*24*total_customers-(time_interruptions(1,2)*time_interruptions(1,1)+time_interruptions(2,2)*time_interruptions(2,1)+time_interruptions(3,2)*time_interruptions(3,1)+time_interruptions(4,2)*time_interruptions(4,1)+time_interruptions(5,2)*time_interruptions(5,1)+time_interruptions(6,2)*time_interruptions(6,1)+time_interruptions(7,2)*time_interruptions(7,1))/60)*100/(365*24*total_customers);                //Average service availability index

+printf('\nSAIFI = %f interruptions/year',SAIFI);

+printf('\nSAIDI = %f minutes/year',SAIDI);

+printf('\nCAIDI = %f minutes/year',CAIDI);

+printf('\nASAI = %f percentage',ASAI);