diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3872/CH9/EX9.3/Ex9_3.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '3872/CH9/EX9.3/Ex9_3.sce')
| -rw-r--r-- | 3872/CH9/EX9.3/Ex9_3.sce | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/3872/CH9/EX9.3/Ex9_3.sce b/3872/CH9/EX9.3/Ex9_3.sce new file mode 100644 index 000000000..b5b8775a5 --- /dev/null +++ b/3872/CH9/EX9.3/Ex9_3.sce @@ -0,0 +1,37 @@ +//Book - Power system: Analysisi & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J.Overbye
+//Chapter-9 ;Example 9.3
+//Scilab Version - 6.0.0; OS - Windows
+
+clc;
+clear;
+
+Xn=0.05 //motor neutral is grounded through reactance in per unit
+Sb=100 //Base value of system in MVA
+Vb=13.8 //Base voltage of system in kV
+Vf=1.05 //Prefault voltage in per unit
+Z0=%i*0.250 //Zero sequence impedance in per unit
+Z1=%i*0.13893 //Positive sequence impedance in per unit
+Z2=%i*0.14562 //Negative sequence impedance in per unit
+Zf=0 //Fault through impedance in per unit
+
+If0=Vf/(Z0+Z1+Z2+(3*Zf)) //sequence line to ground fault current in per unit
+If1=If0;If2=If0; // Since If0=If1=If2
+If=[If0;If1;If2]
+Isf=3*If0 //subtransient fault current in per unit
+Ib2=Sb/(Vb*sqrt(3)) //base current at bus 2 in kA
+Ib22=Isf*Ib2
+Vsf=[0;Vf;0]-([Z0 0 0;0 Z1 0;0 0 Z2]*If) //sequence componenets of the voltages at the fault in per unit
+a=exp(%i*(120)*(%pi/180)); //operator a
+Vlg2=[1 1 1;1 (a^2) a;1 a (a^2)]*Vsf //line to ground voltages at faulted bus 2 in per unit
+for i=1:3 //This loop is included to avoid discrepancies in angle values when the voltage value is near to zero or zero
+ if abs(Vlg2(i))<1e-6 //For example, atand(0,0) gives 0 degree and atand(0,-0) gives 180 degree
+ Vlg2(i)=0;
+ end
+end
+
+printf('The magnitude of subtransient at bus 2 in is %.4f kA and its angle is %.4f degrees\n',abs(Ib22),atand(imag(Ib22),real(Ib22)));
+disp(abs(Vlg2),'The magnitude of line to ground voltages at faulted bus 2 in per unit is:');
+disp(atand(imag(Vlg2),real(Vlg2)),'The angle of line to ground voltages at faulted bus 2 is :');
+
+
|