initial commit / add all books

7 files changed, 261 insertions, 0 deletions

diff --git a/1067/CH23/EX23.03/23_03.sce b/1067/CH23/EX23.03/23_03.sce
new file mode 100755
index 000000000..843625968
--- /dev/null
+++ b/1067/CH23/EX23.03/23_03.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+vf=1;

+r=1250e3;

+V=600;

+x1=.5;

+x2=.5;

+x3=.02;

+ia2=vf/(x1+x2+x3);

+ia=3*ia2;

+ia1=ia2;

+ia0=ia1;

+iab=r/(sqrt(3)*V);

+iab=round(iab/10)*10;

+ia=round(ia*100)/100;

+If=ia*iab;//the difference in result is due to erroneous calculation in textbook.

+printf("fault current If=%fA",If);

+disp("the difference in result is due to erroneous calculation in textbook.")

diff --git a/1067/CH23/EX23.04/23_04.sce b/1067/CH23/EX23.04/23_04.sce
new file mode 100755
index 000000000..9c2228f53
--- /dev/null
+++ b/1067/CH23/EX23.04/23_04.sce
@@ -0,0 +1,33 @@
+clear;

+clc;

+v=1;

+r=1250e3;

+V=600;

+x1=.05*%i;

+x2=.05*%i;

+x0=.02*%i;

+a=1*%e^(%i*(120*%pi/180));

+b=a^2;

+ia1=v/(x1+inv(inv(x2)+inv(x0)));

+ibase=1200;

+va1=v-(ia1*x1);

+ia2=-va1/x2;

+ia0=-va1/x0;

+ia=(ia0+ia1+ia2);

+ib=(ia0+(b*ia1)+(a*ia2));

+ic=(ia0+(a*ia1)+(b*ia2));

+ia0r=real(ia);

+ia0i=imag(ia);

+iam=sqrt((ia0r^2)+(ia0i^2));

+ia1r=real(ib);

+ia1i=imag(ib);

+ibm=sqrt((ia1r^2)+(ia1i^2));//the difference in result is due to erroneous calculation in textbook.

+ia2r=real(ic);

+ia2i=imag(ic);

+icm=sqrt((ia2r^2)+(ia2i^2));

+iaa=0;

+iba=atand(ia1i/ia1r);

+ica=atand(ia2i/ia2r);

+im=ibm*ibase;

+mprintf("fault current for double line to ground fault=%fA",im)

+disp("the difference in result is due to erroneous calculation in textbook.")

diff --git a/1067/CH23/EX23.05/23_05.sce b/1067/CH23/EX23.05/23_05.sce
new file mode 100755
index 000000000..0ad05ca7f
--- /dev/null
+++ b/1067/CH23/EX23.05/23_05.sce
@@ -0,0 +1,33 @@
+clear;

+clc;

+v=1;

+r=1250e3;

+V=600;

+x1=.05*%i;

+x2=.05*%i;

+x0=.02*%i;

+ia1=v/(x1+x2);

+ia2=-ia1;

+ia=ia1+ia2;

+ia0=0;

+a=1*%e^(%i*(120*%pi/180));

+b=a^2;

+ia=(ia0+ia1+ia2);

+ib=(ia0+(b*ia1)+(a*ia2));

+ic=(ia0+(a*ia1)+(b*ia2));

+ia0r=real(ia);

+ia0i=imag(ia);

+iam=sqrt((ia0r^2)+(ia0i^2));

+ia1r=real(ib);

+ia1i=imag(ib);

+ibm=sqrt((ia1r^2)+(ia1i^2));

+ia2r=real(ic);

+ia2i=imag(ic);

+icm=sqrt((ia2r^2)+(ia2i^2));

+iaa=0;

+iba=atand(ia1i/ia1r);

+ica=atand(ia2i/ia2r);

+ibase=r/(sqrt(3)*V);

+ibm=ibm*ibase;

+ibm=round(ibm/100)*100;

+mprintf("fault current for double line to ground fault=%dA",ibm);

diff --git a/1067/CH23/EX23.06/23_06.sce b/1067/CH23/EX23.06/23_06.sce
new file mode 100755
index 000000000..dd24eb573
--- /dev/null
+++ b/1067/CH23/EX23.06/23_06.sce
@@ -0,0 +1,20 @@
+clear;

+clc;

+r=1250e3;

+v=600;

+z1=.15*%i;

+z2=.3*%i;

+z3=.05*%i;

+z4=.55*%i;

+x1=inv(inv(z2)+inv(z1));

+x2=x1;

+x0=inv(inv(z3)+inv(z4));

+e=1;

+ia1=e/(x1+x2+x0);

+ia2=ia1;

+ia0=ia2;

+ia=3*ia1;//the difference in result is due to erroneous calculation in textbook.

+base=r/(sqrt(3)*v);

+ita=ia*base;

+mprintf("the fault current=%fA",-imag(ita));

+disp("the difference in result is due to erroneous calculation in textbook.");

diff --git a/1067/CH23/EX23.07/23_07.sce b/1067/CH23/EX23.07/23_07.sce
new file mode 100755
index 000000000..d0323eba9
--- /dev/null
+++ b/1067/CH23/EX23.07/23_07.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+e=1;

+r=1500e3;

+v=11e3;

+x1=.1;

+ia=3*e/(x1*3);

+ibase=r/(sqrt(3)*v);

+i=ia*ibase;

+mprintf("the single line to ground fault = %dA",i);

+ia1=e/(2*x1);

+ib=sqrt(3)*ia1;

+ib=ibase*ib;

+mprintf("\nline to line fault current=%dA",ib);

diff --git a/1067/CH23/EX23.08/23_08.sce b/1067/CH23/EX23.08/23_08.sce
new file mode 100755
index 000000000..a32e29757
--- /dev/null
+++ b/1067/CH23/EX23.08/23_08.sce
@@ -0,0 +1,50 @@
+clear;

+clc;

+X1=6.6*%i;

+X2=6.3*%i;

+X0=12.6*%i;

+r=37.5e6;

+v=33e3;

+e=1;

+zb=v^2/r;

+x1=X1/zb;

+x2=X2/zb;

+x0=X0/zb;

+x1g=.18*%i;

+x2g=.12*%i;

+x0g=.1*%i;

+x1=x1+x1g;

+x2=x2+x2g;

+x0=x0+x0g;

+ia=3*e/(x1+x2+x0);

+ia1=ia/3;

+a=1*%e^(%i*(120*%pi/180));

+b=a^2;

+ibase=r/(sqrt(3)*v);

+ian=ia*ibase;

+printf("fault current=%djAmp",imag(ian));

+va=e-(ia1*x1g);

+vb=-ia1*x2g;

+vc=-ia1*x0g;

+va0=(va+vb+vc);

+va1=(va+(b*vb)+(a*vc));

+va2=(va+(a*vb)+(b*vc));

+v=v/sqrt(3);

+va0=va0*v;

+va1=va1*v;

+va2=va2*v;

+va0r=real(va0);

+va0i=imag(va0);

+va0m=sqrt((va0r^2)+(va0i^2));

+va0a=atand(va0i/va0r);

+va1r=real(va1);

+va1i=imag(va1);

+va1m=sqrt((va1r^2)+(va1i^2));

+va1a=atand(va1i/va1r)-120;

+va2r=real(va2);

+va2i=imag(va2);

+va2m=sqrt((va2r^2)+(va2i^2));

+va2a=atand(va2i/va2r)+120;

+mprintf("\nthe voltage levels are \n va=%f+j%f V \tor\t %d/_%d kV",va0r/1e3,va0i/1e3,va0m/1e3,va0a);

+mprintf("\n vb=%f+j(%f) kV \tor\t %d/_%d kV",va1r/1e3,va1i/1e3,va1m/1e3,va1a);

+mprintf("\n vc=%f+j(%f) kV \tor\t %d/_%d kV",va2r/1e3,va2i/1e3,va2m/1e3,va2a);

diff --git a/1067/CH23/EX23.09/23_09.sce b/1067/CH23/EX23.09/23_09.sce
new file mode 100755
index 000000000..d8505bb11
--- /dev/null
+++ b/1067/CH23/EX23.09/23_09.sce
@@ -0,0 +1,93 @@
+clear;

+clc;

+e=100/75;

+r=100e6;

+v=66e3;

+xg1=.175*%i*e;

+xg2=.135*%i*e;

+X1=.1*%i*e;

+zn=3*58;

+ibase=r/(sqrt(3)*v);

+vbase=v/sqrt(3);

+zb=vbase/ibase;

+zg0=zn/zb;

+f=70e3;

+e=f/v;

+x1=.367*%i;

+x2=.313*%i;

+z0=zg0+(.133*%i);

+a=1*%e^(%i*(120*%pi/180));

+b=a^2;

+ia1=e/x1;

+mprintf("%f",real(vbase));

+ia=ia1;

+ib=b*ia;

+ic=a*ia;

+ia=ibase*ia;

+ib=ibase*ib;

+ic=ibase*ic;

+ia0r=real(ia);

+ia0i=imag(ia);

+iam=sqrt((ia0r^2)+(ia0i^2));

+ia1r=real(ib);

+ia1i=imag(ib);

+ibm=sqrt((ia1r^2)+(ia1i^2));

+ia2r=real(ic);

+ia2i=imag(ic);

+icm=sqrt((ia2r^2)+(ia2i^2));

+iaa=-90;

+iba=180+atand(ia1i/ia1r);

+ica=atand(ia2i/ia2r);

+mprintf("the symmetric components for three phase fault are \n ia0=%f+j%f A \tor\t %f/_%d A",ia0r,ia0i,iam,iaa);

+mprintf("\n ib=%f+j%f A \tor\t %f/_%d A",ia1r,ia1i,ibm,iba);

+mprintf("\n ic=%f+j(%f) A \tor\t %f/_%d A",ia2r,ia2i,icm,ica);

+ia1=e/(x1+x2);

+ia2=-ia1;

+ia0=0;

+ia=(ia0+ia1+ia2);

+ib=(ia0+(b*ia1)+(a*ia2));

+ic=(ia0+(a*ia1)+(b*ia2));

+i=r/(sqrt(3)*v);

+ia=ia*i;

+ib=ib*i;

+ic=ic*i;

+ia0r=real(ia);

+ia0i=imag(ia);

+iam=sqrt((ia0r^2)+(ia0i^2));

+ia1r=real(ib);

+ia1i=imag(ib);

+ibm=sqrt((ia1r^2)+(ia1i^2));

+ia2r=real(ic);

+ia2i=imag(ic);

+icm=sqrt((ia2r^2)+(ia2i^2));

+iaa=0;

+iba=180+atand(ia1i/ia1r);

+ica=atand(ia2i/ia2r);

+icm=round(icm/10)*10;

+ibm=round(ibm/10)*10;

+mprintf("\nthe symmetric components for line to line fault are \n ia0=%f+j%f A \tor\t %f/_%f A",ia0r,ia0i,iam,iaa);

+mprintf("\n ib=%f+j%f A \tor\t %f/_%f A",ia1r,ia1i,ibm,iba);

+mprintf("\n ic=%f+j(%f) A \tor\t %f/_%f A",ia2r,ia2i,icm,ica);

+ia1=e/(x1+x2+z0);

+ia2=ia1;

+ia0=ia2;

+ia=(ia0+ia1+ia2);

+ib=(ia0+(b*ia1)+(a*ia2));

+ic=(ia0+(a*ia1)+(b*ia2));

+i=r/(sqrt(3)*v);

+ia=ia*874;

+ia0r=real(ia);

+ia0i=imag(ia);

+iam=sqrt((ia0r^2)+(ia0i^2));

+ia1r=real(ib);

+ia1i=imag(ib);

+ibm=sqrt((ia1r^2)+(ia1i^2));

+ia2r=real(ic);

+ia2i=imag(ic);

+icm=sqrt((ia2r^2)+(ia2i^2));

+iaa=atand(ia0i/ia0r);

+iba=0;

+ica=0;

+mprintf("\nthe symmetric components for single line to ground fault are \n ia0=%f+j%f A \tor\t %f/_%f A",ia0r,ia0i,iam,iaa);

+mprintf("\n ib=%f+j%f A \tor\t %f/_%f A",ia1r,ia1i,ibm,iba);

+mprintf("\n ic=%f+j(%f) A \tor\t %f/_%f A",ia2r,ia2i,icm,ica);