10 files changed, 164 insertions, 0 deletions
diff --git a/3793/CH3/EX3.1/exp_3_1.sce b/3793/CH3/EX3.1/exp_3_1.sce
new file mode 100644
index 000000000..34ae376bd
--- /dev/null
+++ b/3793/CH3/EX3.1/exp_3_1.sce
@@ -0,0 +1,9 @@
+clear;
+clc;
+//there are 9 distances and self GMD of conducto will be the ninth root of the product of the nine distance
+d12=2;//consider 2 as 2r r=radius
+d23=2;//consider 2 as 2r r=radius
+d13=2; //consider 2 as 2r r=radius
+r1=.7788; // onsider .7788 as .7788r
+SGMD=((((.7788)^3)*(d12^2)*(d23^2)*(d13^2))^(1/9))
+mprintf("Self-GMD = %.4fr",SGMD);
diff --git a/3793/CH3/EX3.10/exp_3_10.sce b/3793/CH3/EX3.10/exp_3_10.sce
new file mode 100644
index 000000000..227a4e800
--- /dev/null
+++ b/3793/CH3/EX3.10/exp_3_10.sce
@@ -0,0 +1,19 @@
+clear;
+clc;
+//geometric constants
+a=sqrt((5^2)+1);
+b=sqrt(((4+1)^2)+(5^2));
+h=10;
+r=.02;
+s1=4;
+s2=6;
+c=10.773;
+ca=4*%pi*8.854*(10^(-12+9))*150;
+cb=(((a^2)*(b^2)*10*4)/((r^3)*(c^2)*s2))^(1/3);
+Cn=ca/log(cb);
+mprintf(" capacitance to neutral is %.3f microF\n",Cn);
+I=(2*%pi*50*Cn*10^(-3)*220)/(3^(1/2));
+mprintf(" Line Charging Current is %.3f A\n",I);
+Ip=I/2;
+mprintf("Charging Current per conductor is %.3f A",Ip);
+
diff --git a/3793/CH3/EX3.11/exp_3_11.sce b/3793/CH3/EX3.11/exp_3_11.sce
new file mode 100644
index 000000000..03405cdcd
--- /dev/null
+++ b/3793/CH3/EX3.11/exp_3_11.sce
@@ -0,0 +1,7 @@
+clear;
+clc;
+D=(6*6*12)^(1/3);
+ca=2*%pi*8.854*10^(-9);
+cb=log(D/sqrt(.013*.25));
+Cn=ca/cb;
+mprintf("Capacitance to neutral is %.14f F/Km",Cn);
diff --git a/3793/CH3/EX3.2/exp_3_2.sce b/3793/CH3/EX3.2/exp_3_2.sce
new file mode 100644
index 000000000..f4454aba7
--- /dev/null
+++ b/3793/CH3/EX3.2/exp_3_2.sce
@@ -0,0 +1,22 @@
+clear;
+clc;
+Aln=6;
+Dacsr=6;
+f=50;
+DAl=2;
+d1=120//distance
+Dst=(Dacsr-(2*DAl));
+d12=2;
+d16=2
+d=2;
+d13=sqrt(3)*d;
+d15=sqrt(3)*d;
+d14=2*d;
+Ds=((((.7788*d)*(d^2)*((sqrt(3)*d)^2)*(2*d))^6)^(1/36));// in book answer is misprinted of ds
+L=(((2*(10^(-7)))*(log1p(d1/Ds)))*(10^6));
+mprintf("L=%.5f mH/km\n",L);
+Li=2*L;
+mprintf(" Loop inductance = %.5f mH/km\n",Li);
+Xl=(2*(3.14)*(f*Li)*(10^(-3)))
+mprintf(" Inductive Reactance= %.5f ohm/km",Xl);
+
diff --git a/3793/CH3/EX3.3/exp_3_3.sce b/3793/CH3/EX3.3/exp_3_3.sce
new file mode 100644
index 000000000..6e0c9b12a
--- /dev/null
+++ b/3793/CH3/EX3.3/exp_3_3.sce
@@ -0,0 +1,22 @@
+clear;
+clc;
+D=1.200;
+r=.75*(10^2);
+rd=.7788*r;
+Ir=complex(25,-30);
+Iy=complex(35,-50);
+Ib=complex(-60,80);
+fln=2*(10^(-7))*(Ir*log(1/(3*D))+Iy*log(1/(2*D))+Ib*log(1/D));
+mprintf("flux linkage of the neutral = %.9f + i%.9f Wb-T/m\n",real(fln),imag(fln));
+Dn=%i*2*%pi*50*fln*10000;
+LMAT=zeros(3);
+LMAT(1,1)=log((2*D)/rd);
+LMAT(1,2)=log(2);
+LMAT(2,2)=log(D/rd);
+LMAT(3,2)=log(2);
+LMAT(3,3)=LMAT(1,1);
+I=[Ir;Iy;Ib];
+Vryb= %i*200*%pi*(10^(-7))*LMAT*I;
+mprintf(" Voltage drop per unit length = ");
+disp(Vryb);
+
diff --git a/3793/CH3/EX3.4/exp_3_4.sce b/3793/CH3/EX3.4/exp_3_4.sce
new file mode 100644
index 000000000..4bb4f4649
--- /dev/null
+++ b/3793/CH3/EX3.4/exp_3_4.sce
@@ -0,0 +1,14 @@
+clear;
+clc;
+
+//given parameter
+r=1; //radius
+d=3
+h=300;
+p=500;
+q=(100*sqrt((6^2)+(5^2)));
+l=(2*(10^-7)*log(((2^(1/6))*((d/(.7788*r))^(1/2))*((p/q)^(2/3)))));
+L=l*h*1000*1000;
+Xl= 2*%pi*50*L*(10^-3);
+mprintf(" Inductance = %.3f mH\n",L);
+mprintf(" Rectance = %.3f ohm",Xl);
diff --git a/3793/CH3/EX3.5/exp_3_5.sce b/3793/CH3/EX3.5/exp_3_5.sce
new file mode 100644
index 000000000..b49c3ad49
--- /dev/null
+++ b/3793/CH3/EX3.5/exp_3_5.sce
@@ -0,0 +1,20 @@
+clear;
+clc;
+
+f=50;
+d=.05;     //Diameter Of each Conductor in m
+r=d/2;        //radius Of each conductor
+D=.5;     //Space Between Two Conductor in m
+ln=200;      //Distance Of the Line in km
+Dsb=sqrt(.7788*r*D);          
+Dyb=10;
+Dry=10;
+Dbr=20;
+D_eq=((Dyb*Dry*Dbr)^(1/3));         
+L=(2*(10^-7)*log(D_eq/Dsb)*(10^(6)));          
+mprintf(" Inductaance = %.3f mH/Km \n",L);
+l=L*ln*(10^-3)           
+ mprintf("the Inductaance of the line= %.3f H\n",l);
+Xl=2*%pi*l*f;            
+mprintf("the reactance of the line= %.3f ohm",Xl);
+
diff --git a/3793/CH3/EX3.6/exp_3_6.sce b/3793/CH3/EX3.6/exp_3_6.sce
new file mode 100644
index 000000000..a0dee9de9
--- /dev/null
+++ b/3793/CH3/EX3.6/exp_3_6.sce
@@ -0,0 +1,20 @@
+clear;
+clc;
+pe=(%pi*(8.85*10^(-12)));
+v=33000;
+d=1.5//distance
+r=.005 //radius in m
+h=5;
+l=5;
+Cab= ((.01206)/log10((d/r)))*(l);
+Cn=2*Cab;
+I=(2*%pi*50*Cab*v*10^(-6));
+ab=sqrt(1+((d^2)/(4*(h^2))))
+f=log(d/.050559);
+Cbb=((pe/f)*((10^9)*5)/2);
+
+mprintf(" capacitance = %.3f microF\n",Cab);
+mprintf(" capacitance to neutral = %.3f microF\n",Cn);
+mprintf(" charing current = j%.3f A\n",I);
+mprintf(" capacitance with earth cnsideration = %.3f microF\n",Cbb);// capacitance when effect of earth is considered
+
diff --git a/3793/CH3/EX3.8/exp_3_8.sce b/3793/CH3/EX3.8/exp_3_8.sce
new file mode 100644
index 000000000..1dc63a9d4
--- /dev/null
+++ b/3793/CH3/EX3.8/exp_3_8.sce
@@ -0,0 +1,12 @@
+clear;
+clc;
+D=6;
+r=.015;
+e=8.85*10^(-12);
+a=complex(0.5,-0.866);
+ca=log(r/(2*D))-(log(D/(2*r))*a);
+cb=(log((D^2)/(r^2))*log(r/(2*D)))-(log((2*D)/r)*log(D/(2*r)));
+Cyr= (2*%pi*e*(10^9)*ca)/cb;
+mprintf(" capacitance in microF/Km is %.3f %.3fi\n",real(Cyr),imag(Cyr));
+I=%i*2*%pi*50*Cyr*110000*10^-6;
+mprintf(" Current is %.3f+%.3fi A",real(I),imag(I));
diff --git a/3793/CH3/EX3.9/exp_3_9.sce b/3793/CH3/EX3.9/exp_3_9.sce
new file mode 100644
index 000000000..35dd1bca8
--- /dev/null
+++ b/3793/CH3/EX3.9/exp_3_9.sce
@@ -0,0 +1,19 @@
+clear;
+clc;
+Hr1y=15.93;
+Hyb1=12.65;
+Hbr1=15.93;
+Hrr1=18.93;
+Hyy1=12;
+Hbb1=12;
+Deq=4;
+Cb=(2*%pi*(8.854*10^(-12)));
+
+cb1=log(Deq/.02);
+cb2=log(((Hr1y*Hyb1*Hbr1)^(1/3))/((Hrr1*Hyy1*Hbb1)^(1/3)));
+Cn=(Cb/(cb1-cb2));
+
+Cln=Cn*100*10^(9);
+I=((2*%pi*50*Cln*(10^(-3))*220)/sqrt(3));
+mprintf("line capacitance to neutral = %.3f microF\n",Cln);
+mprintf("line current = %.3f A",I);