summaryrefslogtreecommitdiff
path: root/3793/CH10
diff options
context:
space:
mode:
Diffstat (limited to '3793/CH10')
-rw-r--r--3793/CH10/EX10.1/exp_10_1.sce8
-rw-r--r--3793/CH10/EX10.10/exp_10_10.sce62
-rw-r--r--3793/CH10/EX10.2/exp_10_2.sce8
-rw-r--r--3793/CH10/EX10.3/exp_10_3.sce10
-rw-r--r--3793/CH10/EX10.4/exp_10_4.sce63
-rw-r--r--3793/CH10/EX10.5/exp_10_5.sce33
-rw-r--r--3793/CH10/EX10.6/exp_10_6.sce60
-rw-r--r--3793/CH10/EX10.7/exp_10_7.sce57
-rw-r--r--3793/CH10/EX10.8/exp_10_8.sce31
-rw-r--r--3793/CH10/EX10.9/exp_10_9.sce41
10 files changed, 373 insertions, 0 deletions
diff --git a/3793/CH10/EX10.1/exp_10_1.sce b/3793/CH10/EX10.1/exp_10_1.sce
new file mode 100644
index 000000000..a2cc2b3f7
--- /dev/null
+++ b/3793/CH10/EX10.1/exp_10_1.sce
@@ -0,0 +1,8 @@
+clear;
+clc;
+v=[complex(100,0);complex(-75,129.90);complex(-105,-181.865)];
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+vs=inv(A)*v;
+mprintf("symmetrical components of phase voltages are \n");
+disp(vs);
diff --git a/3793/CH10/EX10.10/exp_10_10.sce b/3793/CH10/EX10.10/exp_10_10.sce
new file mode 100644
index 000000000..07622ccf9
--- /dev/null
+++ b/3793/CH10/EX10.10/exp_10_10.sce
@@ -0,0 +1,62 @@
+clear;
+clc;
+Sb=30;
+vb=11;
+sg=20;
+p=10;
+R=6.6;
+//generator
+X1=complex(0,.1);
+X2=complex(0,.1);
+X0=complex(0,.15);
+x1=X1*(Sb/sg);
+x2=X2*(Sb/sg);
+x0=X0*(Sb/sg);
+//transformer12
+xt1=complex(0,.12);
+xt2=complex(0,.12);
+xt0=complex(0,.12);
+//transmission line
+vtr=22;
+Ztr=vtr^2/Sb;
+Z=complex(1,5);
+Zpu=Z/Ztr;
+//transformer34
+Xt1=complex(0,.05);
+Xt2=complex(0,.05);
+Xt0=complex(0,.05);
+xtt1=Xt1*(Sb/sg);
+xtt2=Xt2*(Sb/sg);
+xtt0=Xt0*(Sb/sg);
+Vf3=1;
+Rpu=((Vf3^2)*Sb)/p;
+Rf=(R*Sb)/vtr^2;
+Il=p/Sb;
+Vf4=Vf3+(Il*xtt1);
+Zfp=((x1+xt1+Zpu)*(Rpu+xtt1))/(x1+xt1++Zpu+Rpu+xtt1);
+Zfn=Zfp;
+zf0=Zpu+xt1;
+Ia1=Vf3/complex(1.611,1.5); //from figure
+Ia0=Ia1;
+Ia2=Ia1;
+Ia=3*Ia1;
+Ip23=(Rpu+xtt1)*Ia1/((x1+xt1+Zpu)+Rpu+xtt1);
+Ip34=Ia1-Ip23;
+In23=Ip23;
+In34=Ip34;
+Iz23=Ia1;
+Iz34=0;
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+I23=[Iz23;Ip23;In23];
+II=A*I23;
+mprintf("Current for node 2-3 is Ia=%f%f A Ib=%f+%f A Ic=%f+%f A\n",real(II(1,1)),imag(II(1,1)),real(II(2,1)),imag(II(2,1)),real(II(3,1)),imag(II(3,1)));
+I34=[Iz34;Ip34;In34];
+III=A*I34;
+mprintf("Current for node 3-4 is Ia=%f+%f A Ib=%f%f A Ic=%f%f A\n",real(III(1,1)),imag(III(1,1)),real(III(2,1)),imag(III(2,1)),real(III(3,1)),imag(III(3,1)));
+Vp2=Vf3+(Zpu*Ip23);
+Vn2=Zpu*Ip23;
+Vz2=Zpu*Ia1;
+Vz=[Vz2;Vp2;Vn2];
+V=A*Vz;
+mprintf("Phase voltages at node 2 are Va=%f+%f Kv Vb=%f%f Kv Vc=%f+%f Kv",real(V(1,1)),imag(V(1,1)),real(V(2,1)),imag(V(2,1)),real(V(3,1)),imag(V(3,1)));
diff --git a/3793/CH10/EX10.2/exp_10_2.sce b/3793/CH10/EX10.2/exp_10_2.sce
new file mode 100644
index 000000000..8610c3784
--- /dev/null
+++ b/3793/CH10/EX10.2/exp_10_2.sce
@@ -0,0 +1,8 @@
+clear;
+clc;
+I=[complex(0,0);complex(5.006,-2.89);complex(4.087,4.087)];
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+Ip=A*I;
+mprintf("line currents are \n");
+disp(Ip);
diff --git a/3793/CH10/EX10.3/exp_10_3.sce b/3793/CH10/EX10.3/exp_10_3.sce
new file mode 100644
index 000000000..463404ddd
--- /dev/null
+++ b/3793/CH10/EX10.3/exp_10_3.sce
@@ -0,0 +1,10 @@
+clear;
+clc;
+Zabc=[4 1 1;1 4 1;1 1 4];
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+Z012=inv(A)*Zabc*A;
+mprintf("Sequence impedence matrix\n");
+disp(Z012);
+
+
diff --git a/3793/CH10/EX10.4/exp_10_4.sce b/3793/CH10/EX10.4/exp_10_4.sce
new file mode 100644
index 000000000..75707c17e
--- /dev/null
+++ b/3793/CH10/EX10.4/exp_10_4.sce
@@ -0,0 +1,63 @@
+clear;
+clc;
+Sb=100;
+//generator 1
+x0=complex(0,.08);
+s=100;
+v=11;
+z=3;
+Zb=v^2/s;
+X=z/Zb;
+x=z*X;
+mprintf("Magnitutde of grounding reactor for G1 %.3f\n",x);
+//generator 2
+s1=50;
+x1=complex(0,.05);
+sm=25;
+xm=complex(0,.05);
+z1=x1*(Sb/s1);
+zm=xm*(Sb/sm);
+mprintf("pu reactance for G2 is %.3f\n",imag(z1));
+mprintf("pu reactance for M is %.3f\n",imag(zm));
+//transformers
+s12=100;
+x12=.1;
+z12=x12*(Sb/s12);
+mprintf("pu reactance for T12 is %.3f\n",imag(z12));
+s34=50;
+x34=complex(0,.075);
+z34=x34*(Sb/s34);
+mprintf("pu reactance for T34 is %.3f\n",imag(z34));
+s45=50;
+x45=complex(0,.08);
+z45=x45*(Sb/s45);
+mprintf("pu reactance for T45 is %.3f\n",imag(z45));
+s67=50;
+x67=complex(0,.076);
+z67=x67*(Sb/s67);
+mprintf("pu reactance for T67 is %.3f\n",imag(z67));
+s78=75;
+x78=complex(0,.06);
+z78=x78*(Sb/s78);
+mprintf("pu reactance for T78 is %.3f\n",imag(z78));
+vb=132;
+zb=vb^2/s;
+X11=z/zb;
+x11=z*X11;
+mprintf("pu reactance for zero sequence is %.3f\n",x11);
+Zt23=vb^2/s;
+xt23=complex(0,300);
+zt23=xt23/Zt23;
+mprintf("pu reactance for Tr23 is %.3f\n",imag(zt23));
+Zt28=vb^2/s;
+xt28=complex(0,250);
+zt28=xt28/Zt28;
+mprintf("pu reactance for Tr28 is %.3f\n",imag(zt28));
+vt=66;
+Zt56=vt^2/s;
+xt56=complex(0,200);
+zt56=xt56/Zt56;
+mprintf("pu reactance for Tr56 is %.3f\n",imag(zt56));
+
+
+
diff --git a/3793/CH10/EX10.5/exp_10_5.sce b/3793/CH10/EX10.5/exp_10_5.sce
new file mode 100644
index 000000000..62c8e4af8
--- /dev/null
+++ b/3793/CH10/EX10.5/exp_10_5.sce
@@ -0,0 +1,33 @@
+clear;
+clc;
+s=30;
+v=11;
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+xs=complex(0,.5);
+xn=complex(0,.3);
+xz=complex(0,.08);
+Z=5;
+X=(Z*s)/v^2;
+Xn=complex(0,3*X);
+Ea=complex(1,0);
+Ia1=Ea/(Xn+xn+xs+xz);
+Ia=3*Ia1;
+Va0=-Ia1*(xz+Xn);
+Va1=Ea-(Ia1*xs);
+Va2=-Ia1*xn;
+V1=[Va0;Va1;Va2];
+V=A*V1;
+Vab=V(1,1)-V(2,1);
+Vbc=V(2,1)-V(3,1);
+Vca=V(3,1)-V(1,1);
+Ibase=(s*1000)/(sqrt(3)*v);
+If=Ia1*Ibase;
+Vab1=Vab*(v/sqrt(3));
+Vbc1=Vbc*(v/sqrt(3));
+Vca1=Vca*(v/sqrt(3));
+mprintf("Sub transient fault current is %.3f A \n",imag(If));
+mprintf("Actual Line Voltages are Vab=%f+%f Kv Vbc=%f%f Kv Vca=%f+%f Kv",real(Vab1),imag(Vab1),real(Vbc1),imag(Vbc1),real(Vca1),imag(Vca1));
+
+
+
diff --git a/3793/CH10/EX10.6/exp_10_6.sce b/3793/CH10/EX10.6/exp_10_6.sce
new file mode 100644
index 000000000..e0b2e8bf7
--- /dev/null
+++ b/3793/CH10/EX10.6/exp_10_6.sce
@@ -0,0 +1,60 @@
+clear;
+clc;
+Sb=30;
+vb=11;
+sg=20;
+p=10;
+R=6.6;
+//generator
+X1=complex(0,.1);
+X2=complex(0,.1);
+X0=complex(0,.15);
+x1=X1*(Sb/sg);
+x2=X2*(Sb/sg);
+x0=X0*(Sb/sg);
+//transformer12
+xt1=complex(0,.12);
+xt2=complex(0,.12);
+xt0=complex(0,.12);
+//transmission line
+vtr=22;
+Ztr=vtr^2/Sb;
+Z=complex(1,5);
+Zpu=Z/Ztr;
+//transformer34
+Xt1=complex(0,.05);
+Xt2=complex(0,.05);
+Xt0=complex(0,.05);
+xtt1=Xt1*(Sb/sg);
+xtt2=Xt2*(Sb/sg);
+xtt0=Xt0*(Sb/sg);
+Vf3=1;
+Rpu=((Vf3^2)*Sb)/p;
+Rf=(R*Sb)/vtr^2;
+Il=p/Sb;
+Vf4=Vf3+(Il*xtt1);
+Zfp=((x1+xt1+Zpu)*(Rpu+xtt1))/(x1+xt1++Zpu+Rpu+xtt1);
+Zfn=Zfp;
+zf0=Zpu+xt1;
+Ia1=Vf3/complex(1.611,1.5); //from figure
+Ia0=Ia1;
+Ia=3*Ia1;
+Va0=-zf0*Ia0;
+Va1=Vf3-(Zfn*Ia1);
+Va2=-(Zfn*Ia1);
+Ibase=(Sb*1000)/(sqrt(3)*vtr);
+If=Ibase*Ia;
+mprintf("fault current is %f%f A \n",real(If),imag(If));
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+V1=[Va0;Va1;Va2];
+V=A*V1;
+Vab=V(1,1)-V(2,1);
+Vbc=V(2,1)-V(3,1);
+Vca=V(3,1)-V(1,1);
+mprintf("Actual Line Voltages are Vab=%f+%f Kv Vbc=%f%f Kv Vca=%f+%f Kv",real(Vab),imag(Vab),real(Vbc),imag(Vbc),real(Vca),imag(Vca));
+
+
+
+
+
diff --git a/3793/CH10/EX10.7/exp_10_7.sce b/3793/CH10/EX10.7/exp_10_7.sce
new file mode 100644
index 000000000..39e02844a
--- /dev/null
+++ b/3793/CH10/EX10.7/exp_10_7.sce
@@ -0,0 +1,57 @@
+clear;
+clc;
+Sb=30;
+vb=11;
+sg=20;
+p=10;
+R=6.6;
+Ea=1;
+//generator
+X1=complex(0,.1);
+X2=complex(0,.1);
+X0=complex(0,.15);
+x1=X1*(Sb/sg);
+x2=X2*(Sb/sg);
+x0=X0*(Sb/sg);
+//transformer12
+xt1=complex(0,.12);
+xt2=complex(0,.12);
+xt0=complex(0,.12);
+//transmission line
+vtr=22;
+Ztr=vtr^2/Sb;
+Z=complex(1,5);
+Zpu=Z/Ztr;
+//transformer34
+Xt1=complex(0,.05);
+Xt2=complex(0,.05);
+Xt0=complex(0,.05);
+xtt1=Xt1*(Sb/sg);
+xtt2=Xt2*(Sb/sg);
+xtt0=Xt0*(Sb/sg);
+Vf3=1;
+Rpu=((Vf3^2)*Sb)/p;
+Rf=(R*Sb)/vtr^2;
+Il=p/Sb;
+Vf4=Vf3+(Il*xtt1);
+Zfp=((x1+xt1+Zpu)*(Rpu+xtt1))/(x1+xt1++Zpu+Rpu+xtt1);
+Zfn=Zfp;
+zf0=Zpu+xt1;
+Ia1=Vf3/(Zfp+Zfn+Rf);
+Ia2=-Ia1;
+Va0=0;
+Va1=Ea-(Zfn*Ia1);
+Va2=-Zfn*Ia2;
+Ibase=(Sb*1000)/(sqrt(3)*vtr);
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+V1=[Va0;Va1;Va2];
+V=A*V1;
+Vab=V(1,1)-V(2,1);
+Vbc=V(2,1)-V(3,1);
+Vca=V(3,1)-V(1,1);
+mprintf("Actual Line Voltages are Vab=%f+%f Vbc=%f%f Vca=%f+%f \n",real(Vab),imag(Vab),real(Vbc),imag(Vbc),real(Vca),imag(Vca));
+
+Ib1=complex(0,-sqrt(3))*Ia1;
+Ib=Ib1*Ibase;
+mprintf("Phase current is %f%f A",real(Ib),imag(Ib));
diff --git a/3793/CH10/EX10.8/exp_10_8.sce b/3793/CH10/EX10.8/exp_10_8.sce
new file mode 100644
index 000000000..973e0229f
--- /dev/null
+++ b/3793/CH10/EX10.8/exp_10_8.sce
@@ -0,0 +1,31 @@
+clear;
+clc;
+s=30;
+v=22;
+E=1;
+Zf=6.6;
+Zg=13.2;
+Z0=complex(.062,.43);
+Z1=complex(.161,.535);
+Z2=Z1;
+Zff=.409;
+Zgp=(Zg*s)/v^2;
+Ia1=E/(Z1+Zff+((Z2+Zff)*(Z0+Zff+3*Zgp)/(Z2+Z0+2*Zff+3*Zgp)));
+Ia2=-((Z0+Zff+3*Zgp)/(Z2+Z0+2*Zff+3*Zgp))*(Ia1);
+Ia0=-((Z2+Zff)/(Z2+Z0+2*Zff+3*Zgp))*(Ia1);
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+I=[Ia0;Ia1;Ia2];
+Ia=A*I;
+Va0=-Z0*Ia0;
+Va1=E-Z1*Ia1;
+Va2=-Z1*Ia2;
+V1=[Va0;Va1;Va2];
+V=A*V1;
+Vab=V(1,1)-V(2,1);
+Vbc=V(2,1)-V(3,1);
+Vca=V(3,1)-V(1,1);
+mprintf("Actual Phase Voltages are Va0=%f+%f Va1=%f%f Va2=%f+%f \n",real(Va0),imag(Va0),real(Va1),imag(Va1),real(Va2),imag(Va2));
+mprintf("Actual Phase currents are \n");
+
+disp(Ia);
diff --git a/3793/CH10/EX10.9/exp_10_9.sce b/3793/CH10/EX10.9/exp_10_9.sce
new file mode 100644
index 000000000..2c39d58a6
--- /dev/null
+++ b/3793/CH10/EX10.9/exp_10_9.sce
@@ -0,0 +1,41 @@
+clear;
+clc;
+s=30;
+v=22;
+E=1;
+Zf=6.6;
+Zg=13.2;
+Z0=complex(.062,.43);
+Z1=complex(.161,.535);
+Z2=Z1;
+Zff=.409;
+Zgp=(Zg*s)/v^2;
+Ia=E/(Z1+Zff);
+absIa=abs(Ia);
+mprintf("Absolute Ia=%.3f\n",absIa);
+angleIa=atand(imag(Ia)/real(Ia));
+mprintf("Abngle of Ia=%.3f\n",angleIa);
+a=complex((-.5),.866);
+aa=a^2;
+anglea=atand(imag(aa)/real(aa));
+angleaa=atand(imag(a)/real(a));
+angleIb=angleIa+(anglea);
+disp(anglea)
+
+mprintf("Angle of Ib=%.3f\n",angleIb);
+angleIc=angleIa+(angleaa);
+
+mprintf("Angle of Ic=%.3f\n",angleIc);
+Va=E-(Zff*Ia);
+absVa=abs(Va);
+mprintf("Absolute Va=%.3f\n",absVa);
+angleVa=atand(imag(Va)/real(Va));
+mprintf("Angle of Va=%.3f\n",angleVa);
+angleVb=angleVa+(anglea);
+mprintf("Angle of Vb=%.3f\n",angleVb);
+angleVc=angleVa+(angleaa);
+mprintf("Angle of Vc=%.3f",angleVc);
+
+
+
+