initial commit / add all books

20 files changed, 346 insertions, 0 deletions

diff --git a/1076/CH2/EX2.1/2_1.sce b/1076/CH2/EX2.1/2_1.sce
new file mode 100755
index 000000000..ec4a758c8
--- /dev/null
+++ b/1076/CH2/EX2.1/2_1.sce
@@ -0,0 +1,10 @@
+clear

+clc;

+dia=1.213;

+dist=1.25*100;

+f=50;

+rad=dia/2;

+effrad=.7788*rad;

+L=4d-7 * log (dist/effrad)*1d3;

+X=2*%pi*f*L;

+mprintf("L=%.2f *1e-4 H/km, X=%.1f ohm/km",L*1e4, X);

diff --git a/1076/CH2/EX2.10/2_10.sce b/1076/CH2/EX2.10/2_10.sce
new file mode 100755
index 000000000..cd9847588
--- /dev/null
+++ b/1076/CH2/EX2.10/2_10.sce
@@ -0,0 +1,23 @@
+clearglobal;

+clc;

+

+V=220e3;

+f=50;

+l=200;

+d1=6;

+d2=12;

+rad=1.81e-2;

+

+Deq=(d1*d1*d2)^(1/3);

+Cn=.02412/log10(Deq/rad);

+mprintf("Capacitance per phase per km = %.1f *1e-3 e-6F/phase/km \n",Cn *1000);

+C=l*Cn;

+C=round(C*100)/100

+mprintf("Capacitance per phase = %.2f e-6F/phase\n",C);

+Xc=1/(2 * %pi * f * C *1e-6);

+mprintf("Capacitive reactance per phase = %.0f ohms/phase\n",Xc);

+I=2 * %pi * f * C *1e-6 * V / sqrt(3);

+I=round(I*100)/100

+mprintf("Charging current = %.2f A/phase\n",I);

+MVA=sqrt(3)*V *I *1e-6;

+mprintf("Charging MVA = %.2f MVA\n",fix(MVA*100)/100);

diff --git a/1076/CH2/EX2.11/2_11.sce b/1076/CH2/EX2.11/2_11.sce
new file mode 100755
index 000000000..8c1e3b7c6
--- /dev/null
+++ b/1076/CH2/EX2.11/2_11.sce
@@ -0,0 +1,10 @@
+clear;

+clc;

+

+r= 1.6e-2;

+d= 45e-2;

+D=12;

+Dscb=sqrt(r*d);

+Deq=(D*D*(2*D))^(1/3);

+Cn= .02412/(log10 (Deq/Dscb));

+mprintf("Capacitance per phase per km= %.4f e-6 F/km\n",Cn);

diff --git a/1076/CH2/EX2.12/2_12.sce b/1076/CH2/EX2.12/2_12.sce
new file mode 100755
index 000000000..16de703ed
--- /dev/null
+++ b/1076/CH2/EX2.12/2_12.sce
@@ -0,0 +1,23 @@
+clear

+clc;

+

+D=350;

+r=1.09;

+f=50;

+V=132e3;

+d=100;

+

+//(b)

+Cn=.04824 / log10((sqrt(3)*D)/(2*r));

+Cn=round(Cn*10000)/1e4

+C=Cn/2;

+mprintf("Capacitance per conductor per km = %.5f e-6 F/conductor/km\n",C);

+

+//(c)

+w=2 * %pi * f;

+Vn=V/sqrt(3);

+Ic= w * Cn * Vn * 1e-6;

+Ic=round(Ic *1e3)/1e3

+I=Ic * d;

+MVA= sqrt(3)*V * I / 1e6;

+mprintf("Charging MVA = %.2f MVA\n",MVA);

diff --git a/1076/CH2/EX2.13/2_13.sce b/1076/CH2/EX2.13/2_13.sce
new file mode 100755
index 000000000..f469cf816
--- /dev/null
+++ b/1076/CH2/EX2.13/2_13.sce
@@ -0,0 +1,23 @@
+clear

+

+clc;

+

+r=.9e-2;

+d=6;

+D1=6;

+D2=7;

+

+Daa=sqrt ((d*d)+(D1*D1));

+Daa=round(Daa*1e3)/1e3

+Dbb=D2;

+Dcc=Daa;

+Deq=5.074;

+

+Dsc_a = round(sqrt(r*Daa)*1e4)/1e4;

+Dsc_b = round(sqrt(r*Dbb)*1e4)/1e4;

+Dsc_c = round(sqrt(r*Dcc)*1e4)/1e4;

+

+Dsc=(Dsc_a*Dsc_b*Dsc_c)^(1/3);

+Cn=.02412/log10(Deq/Dsc);

+

+mprintf("Capacitance per conductor per km = %f e-6 F/conductor/km\n",Cn);

diff --git a/1076/CH2/EX2.14/2_14.sce b/1076/CH2/EX2.14/2_14.sce
new file mode 100755
index 000000000..414ccca91
--- /dev/null
+++ b/1076/CH2/EX2.14/2_14.sce
@@ -0,0 +1,13 @@
+clear;

+clc;

+

+h=5.5;

+Dia=1.213e-2;

+d=1.25;

+l=10;

+

+r=Dia/2;

+Cn=(2*.01206)/(log10 (d/(Dia* 0.5 * sqrt(1 + (d*d*.25/h^2)))));

+C=Cn * l;

+

+mprintf("Capacitance per conductor= %.5f e-6 F/conductor\n",C);

diff --git a/1076/CH2/EX2.15/2_15.sce b/1076/CH2/EX2.15/2_15.sce
new file mode 100755
index 000000000..c9caf8dfb
--- /dev/null
+++ b/1076/CH2/EX2.15/2_15.sce
@@ -0,0 +1,18 @@
+clear

+clc;

+

+dia=4.22e-3;

+n=6;

+T1=20;

+T2=50;

+rho1 = 2.826e-8;

+a=.004;

+dL=1.5e-2;

+L=1000;

+

+A=%pi * n *dia *dia / 4;

+R1= rho1 * L * (1+dL) / A;

+R2= R1 * (1+ (a * (T2-T1)));

+

+mprintf("Temperature at %d = %.4f ohm/km\n", T1, R1);

+mprintf("Temperature at %d = %.4f ohm/km", T2, R2);

diff --git a/1076/CH2/EX2.16/2_16.sce b/1076/CH2/EX2.16/2_16.sce
new file mode 100755
index 000000000..b68216c3b
--- /dev/null
+++ b/1076/CH2/EX2.16/2_16.sce
@@ -0,0 +1,25 @@
+clear;

+clc;

+

+A= 1.5e-4;

+Deff=39.8e-3;

+D= 8;

+rho1=1.73e-6 / 100;

+l=1e3;

+f=50;

+V=132e3;

+

+//(a)

+R= rho1 * l / A;

+r=.5 * Deff;

+L= .4605 * log10 (D/(.7788 *r));

+mprintf("L = %.2f mH/km\n",L);

+

+C= .02412/(log10 (D/r));

+mprintf("C = %.5f e-6 F/km\n",C);

+

+//(b)

+Ic = 2 * %pi * f * C *1e-6 * V / sqrt(3);

+mprintf("Charging current = %.4f A/km/phase\n",Ic);

+MVA=sqrt(3)*V *Ic *1e-6;

+mprintf("Charging MVA = %.4f MVA/km\n",MVA);

diff --git a/1076/CH2/EX2.17/2_17.sce b/1076/CH2/EX2.17/2_17.sce
new file mode 100755
index 000000000..80df46ec5
--- /dev/null
+++ b/1076/CH2/EX2.17/2_17.sce
@@ -0,0 +1,11 @@
+clear;

+clc;

+

+r=.9e-2;

+d=3.5;

+

+

+//(b)

+

+L= .4605 * log10(1.375 * sqrt(d/r));

+mprintf("(b) L= %.2f mH/km", L);

diff --git a/1076/CH2/EX2.18/2_18.sce b/1076/CH2/EX2.18/2_18.sce
new file mode 100755
index 000000000..6f5f3ca8d
--- /dev/null
+++ b/1076/CH2/EX2.18/2_18.sce
@@ -0,0 +1,19 @@
+clear;

+clc;

+clear;

+

+

+dia=26.88;

+d=450

+d1=15.25*1e3;

+

+R=dia/2;

+R1= .7788 * R;

+R1=round(R1*100)/100

+GMR= 1.09 * (R1 * d*d*d)^(.25);

+GMR=round(GMR*10)/10

+GMD=(d1* d1 * (2*d1))^(1/3);

+GMD=round(GMD/10)*10

+L=.4605 * log10 (GMD/GMR);

+

+mprintf(" L= %.4f mH/km", L);

diff --git a/1076/CH2/EX2.19/2_19.sce b/1076/CH2/EX2.19/2_19.sce
new file mode 100755
index 000000000..3e03557ad
--- /dev/null
+++ b/1076/CH2/EX2.19/2_19.sce
@@ -0,0 +1,25 @@
+clear

+clc

+

+dia=2.5e-2

+d=3

+

+r=dia/2

+r1=.7788*r

+

+c=exp(%i *2*%pi/3)

+b=exp(%i *-2*%pi/3)

+k=2

+

+Dab=d

+Dac=2*d

+Dbc=d

+Dca=2*d

+

+La=round(k*(log(1/r1) + log(1/Dab)*b + log(1/Dac)*c)*1e3)/1e4

+Lb=round(k*(log(1/Dab) + log(1/r1)*b + log(1/Dbc)*c)/b*1e3)/(1e4)

+Lc=round(k*(log(1/Dac) + log(1/Dbc)*b + log(1/r1)*c)/c*1e3)/(1e4)

+

+disp(La, "La= (mH/km)")

+disp(Lb, "Lb= (mH/km)")

+disp(Lc, "Lc= (mH/km)")

diff --git a/1076/CH2/EX2.20/2_20.sce b/1076/CH2/EX2.20/2_20.sce
new file mode 100755
index 000000000..72ad4368c
--- /dev/null
+++ b/1076/CH2/EX2.20/2_20.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+

+dia=26.88e-3;

+d=45e-2;

+d1=15.25;

+

+r=dia/2;

+

+GMR= 1.09 * (r * d*d*d)^(.25);

+GMD=(d1* d1 * (2*d1))^(1/3);

+

+C= 0.02412 / log10 (GMD/GMR);

+mprintf("C= %.4f e-6 F/km", C);

diff --git a/1076/CH2/EX2.21/2_21.sce b/1076/CH2/EX2.21/2_21.sce
new file mode 100755
index 000000000..c510227f5
--- /dev/null
+++ b/1076/CH2/EX2.21/2_21.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+f=50;

+dia=25e-3;

+sp=.3;

+D=6;

+

+rad=dia/2;

+r=.7788*rad;

+GMR= ((r*r*sp*sp)^.25);

+Dab= (D* (D+sp)* D *(D-sp))^.25;

+Dbc=Dab;

+Dac= ((2*D)* ((2*D) +sp)* (2*D) *((2*D)-sp))^.25;

+Deq=(Dab * Dbc *Dac)^(1/3);

+L=.4605 * log10 (Deq/GMR);

+mprintf("\nL= %.3f mH/phase/km", L);

+XL=2 * %pi * f * L *1e-3;

+mprintf("\nXL= %.2f ohm/phase/km", XL);

diff --git a/1076/CH2/EX2.3/2_3.sce b/1076/CH2/EX2.3/2_3.sce
new file mode 100755
index 000000000..7066c2ca4
--- /dev/null
+++ b/1076/CH2/EX2.3/2_3.sce
@@ -0,0 +1,9 @@
+clear

+clc;

+Dia=1.63;

+Dist=3*100;

+Reff=.7788*.5*Dia;

+Reff=round(Reff*1e3)*1e-3

+L = .4605 * log10(Dist/Reff);

+

+mprintf("L= %.4f mH/km",L)

diff --git a/1076/CH2/EX2.4/2_4.sce b/1076/CH2/EX2.4/2_4.sce
new file mode 100755
index 000000000..e12af0c3a
--- /dev/null
+++ b/1076/CH2/EX2.4/2_4.sce
@@ -0,0 +1,10 @@
+clear

+clc;

+rad=1.81;

+Dist1=6*100;

+Dist2=Dist1;

+Dist3=12*100;

+Reff=.7788*rad;

+Dist=(Dist1*Dist2*Dist3)^(1/3);

+L = .4605 * log10(Dist/Reff);

+mprintf("L=%.4f mH/km",L);

diff --git a/1076/CH2/EX2.5/2_5.sce b/1076/CH2/EX2.5/2_5.sce
new file mode 100755
index 000000000..843a0c8ca
--- /dev/null
+++ b/1076/CH2/EX2.5/2_5.sce
@@ -0,0 +1,7 @@
+clear

+clc;

+d=25;

+rad=1.573;

+Ds=.7788*rad;

+Dsb=((Ds*d*d*sqrt(2)*d)^4)^(1/16);

+mprintf("GMR= %.3f cm",fix(Dsb*1000)/1000);

diff --git a/1076/CH2/EX2.6/2_6.sce b/1076/CH2/EX2.6/2_6.sce
new file mode 100755
index 000000000..c55ac29e4
--- /dev/null
+++ b/1076/CH2/EX2.6/2_6.sce
@@ -0,0 +1,23 @@
+clear

+clc;

+d=45;

+D=12e2;

+rad=1.6;

+

+//(a)

+reff=.7788*rad;

+GMR=(reff*d*reff*d)^(1/4);

+Dab=(D*(D+d)*D*(D-d))^(1/4);

+Dbc=Dab;

+Dca=((D*2)*((2*D)+d)*(D*2)*((2*D)-d))^(1/4);

+Deq=(Dab*Dbc*Dca)^(1/3);

+L=.4605 *log10 (Deq/GMR);

+mprintf("\n(a)\tL=%.2f mH/km",L);

+

+//(b)

+R=sqrt(2)*rad;

+reff=R*.7788;

+Deq=(D*D*(2*D))^(1/3);

+L=.4605 *log10 (Deq/reff);

+

+mprintf("\n(b)\tL=%.3f mH/km",L);

diff --git a/1076/CH2/EX2.7/2_7.sce b/1076/CH2/EX2.7/2_7.sce
new file mode 100755
index 000000000..eda50a1b8
--- /dev/null
+++ b/1076/CH2/EX2.7/2_7.sce
@@ -0,0 +1,11 @@
+clear

+clc;

+D=3e2;

+h=5e2;

+rad=1.1;

+

+m=(h^2  +  D^2)^(1/2);

+n=(h^2  +  (D*2)^2)^(1/2);

+reff=.7788* rad;

+L=2e-7 * log ((2)^(1/6)* (D/reff)^(1/2) * (m/n)^(1/3)) *1e3;

+mprintf("L=%.2f * 1e-4 H/phase/km",L*1e4);

diff --git a/1076/CH2/EX2.8/2_8.sce b/1076/CH2/EX2.8/2_8.sce
new file mode 100755
index 000000000..ef3f51e7d
--- /dev/null
+++ b/1076/CH2/EX2.8/2_8.sce
@@ -0,0 +1,41 @@
+clear;

+clc;

+D=3e2;

+d1=6e2;

+d2=7e2;

+rad=.9;

+

+reff=.7788* rad;

+Daa=(d1^2  +  d1^2)^(1/2);

+Dcc=Daa;

+Dbb=d2;

+GMRa=sqrt(reff*Daa);

+GMRb=sqrt(reff*Dbb);

+GMRc=sqrt(reff*Dcc);

+Ds=(GMRa*GMRb*GMRc)^(1/3);

+Ds=round(Ds*10)/10

+

+Dab=(D^2  +  ((d2-d1)/2)^2)^(1/2);

+Dcb=Dab;

+Dc1b1=Dab;

+Da1b1=Dab;

+

+Dab1=(D^2  +  (((d2-d1)/2)+d1)^2)^(1/2);

+Da1b=Dab1;

+Dc1b=Dab1;

+Dcb1=Dab1;

+

+Dac=2*D;

+Da1c1=Dac;

+Da1c=(d1);

+Dac1=Da1c;

+

+GMRab=(Dab*Da1b1*Da1b*Dab1)^(1/4);

+GMRbc=(Dcb*Dc1b1*Dc1b*Dcb1)^(1/4);

+GMRac=(Dac*Da1c1*Da1c*Dac1)^(1/4);

+

+Deq=(GMRab*GMRbc*GMRac)^(1/3);

+Deq=round(Deq*10)/10

+

+L=2e-7 * log (Deq/Ds) * 1e3;

+mprintf("L=%.3f *1e-4 H/phase/km",L*1e4);

diff --git a/1076/CH2/EX2.9/2_9.sce b/1076/CH2/EX2.9/2_9.sce
new file mode 100755
index 000000000..3f230030a
--- /dev/null
+++ b/1076/CH2/EX2.9/2_9.sce
@@ -0,0 +1,13 @@
+clear;

+clc;

+

+l=10;

+D=1.25;

+dia=1.213e-2;

+

+Cab=.01206/log10(D/(dia/2))

+C=l*Cab;

+Cn=2*C;

+

+mprintf("Capacitance between 2 conductors= %.4f e-6F/km\n", C);

+mprintf("Capacitance between conductor and neutral= %.4f e-6F/km\n", Cn);