initial commit / add all books

2 files changed, 103 insertions, 0 deletions

diff --git a/83/CH9/EX9.6/example_9_6.sce b/83/CH9/EX9.6/example_9_6.sce
new file mode 100755
index 000000000..2ab3fd629
--- /dev/null
+++ b/83/CH9/EX9.6/example_9_6.sce
@@ -0,0 +1,73 @@
+//Chapter 9
+//Example 9.6
+//page 352
+//To calculate short circuit solution using algorithm for short circuit studies
+clear;clc;
+
+Y11=1/(0.15*%i)+1/(0.15*%i)+1/(0.1*%i)+1/(0.2*%i);
+Y12=-1/(0.2*%i);
+Y21=Y12;
+Y13=-1/(0.15*%i);
+Y31=Y13;
+Y14=-1/(0.1*%i);
+Y41=Y14;
+Y22=1/(0.15*%i)+1/(0.15*%i)+1/(0.1*%i)+1/(0.2*%i);
+Y23=-1/(0.1*%i);
+Y32=Y23;
+Y24=-1/(0.15*%i);
+Y42=Y24;
+Y33=1/(0.15*%i)+1/(0.1*%i);
+Y34=0;
+Y43=Y34;
+Y44=1/(0.15*%i)+1/(0.1*%i);
+
+//Ybus matrix can be written as
+
+Ybus=[Y11 Y12 Y13 Y14;Y21 Y22 Y23 Y24;Y31 Y32 Y33 Y34;Y41 Y42 Y43 Y44];
+
+Zbus=inv(Ybus);
+
+//preault voltages
+V10=1;V20=1;V30=1;V40=1;
+
+//post fault voltages
+V1f=V10-(Zbus(1,4)/Zbus(4,4))*V40;
+V2f=V20-(Zbus(2,4)/Zbus(4,4))*V40;
+V3f=V30-(Zbus(3,4)/Zbus(4,4))*V40;
+V4f=V40-(Zbus(4,4)/Zbus(4,4))*V40;
+
+//to calculate fault current through Zf=0
+If=V40/(Zbus(4,4)+0);
+
+//short circuit current in lines 1-3,1-2,1-4,2-4 and 2-3
+
+I13f=(V1f-V3f)/(0.15*%i);
+I12f=(V1f-V2f)/(0.2*%i);
+I14f=(V1f-V4f)/(0.1*%i);
+I24f=(V2f-V4f)/(0.15*%i);
+I23f=(V2f-V3f)/(0.1*%i);
+
+//If at all fault occurs on bus1 or bus2
+If12=1/Zbus(1,1);
+
+//displaying the results
+printf('\n Ybus=');
+disp(Ybus);
+
+printf('\n Zbus=');
+disp(Zbus);
+
+printf('\nV1f= %0.4f pu',V1f);
+printf('\nV2f= %0.4f pu',V2f);
+printf('\nV3f= %0.4f pu',V3f);
+printf('\nV4f= %0.1f pu\n',V4f);
+
+printf('\nFault current=-j%0.5f pu\n',abs(If));
+
+printf('\nI13f=j%0.3f pu',abs(I13f));
+printf('\nI12f=j%0.3f pu',abs(I12f));
+printf('\nI14f=-j%0.3f pu',abs(I14f));
+printf('\nI24f=-j%0.3f pu',abs(I24f));
+printf('\nI23f=-j%0.3f pu\n',abs(I23f));
+
+printf('\n Fault current for a fault on bus 1 (or bus 2)\n If=-j%0.6f pu\n\n',abs(If12));
\ No newline at end of file
diff --git a/83/CH9/EX9.6/result_example_9_6.txt b/83/CH9/EX9.6/result_example_9_6.txt
new file mode 100755
index 000000000..81e1467ae
--- /dev/null
+++ b/83/CH9/EX9.6/result_example_9_6.txt
@@ -0,0 +1,30 @@
+ 
+
+ Ybus= 
+  - 28.333333i    5.i           6.6666667i    10.i        
+    5.i         - 28.333333i    10.i          6.6666667i  
+    6.6666667i    10.i        - 16.666667i    0           
+    10.i          6.6666667i    0           - 16.666667i  
+
+ Zbus= 
+    0.0903061i    0.0596939i    0.0719388i    0.0780612i  
+    0.0596939i    0.0903061i    0.0780612i    0.0719388i  
+    0.0719388i    0.0780612i    0.1356122i    0.0743878i  
+    0.0780612i    0.0719388i    0.0743878i    0.1356122i  
+
+V1f= 0.4244 pu
+V2f= 0.4695 pu
+V3f= 0.4515 pu
+V4f= 0.0 pu
+
+Fault current=-j7.37397 pu
+
+I13f=j0.181 pu
+I12f=j0.226 pu
+I14f=-j4.244 pu
+I24f=-j3.130 pu
+I23f=-j0.181 pu
+
+ Fault current for a fault on bus 1 (or bus 2)
+ If=-j11.073446 pu
+