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);