From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 1076/CH2/EX2.1/2_1.sce | 10 ++++++++++ 1076/CH2/EX2.10/2_10.sce | 23 +++++++++++++++++++++++ 1076/CH2/EX2.11/2_11.sce | 10 ++++++++++ 1076/CH2/EX2.12/2_12.sce | 23 +++++++++++++++++++++++ 1076/CH2/EX2.13/2_13.sce | 23 +++++++++++++++++++++++ 1076/CH2/EX2.14/2_14.sce | 13 +++++++++++++ 1076/CH2/EX2.15/2_15.sce | 18 ++++++++++++++++++ 1076/CH2/EX2.16/2_16.sce | 25 +++++++++++++++++++++++++ 1076/CH2/EX2.17/2_17.sce | 11 +++++++++++ 1076/CH2/EX2.18/2_18.sce | 19 +++++++++++++++++++ 1076/CH2/EX2.19/2_19.sce | 25 +++++++++++++++++++++++++ 1076/CH2/EX2.20/2_20.sce | 14 ++++++++++++++ 1076/CH2/EX2.21/2_21.sce | 18 ++++++++++++++++++ 1076/CH2/EX2.3/2_3.sce | 9 +++++++++ 1076/CH2/EX2.4/2_4.sce | 10 ++++++++++ 1076/CH2/EX2.5/2_5.sce | 7 +++++++ 1076/CH2/EX2.6/2_6.sce | 23 +++++++++++++++++++++++ 1076/CH2/EX2.7/2_7.sce | 11 +++++++++++ 1076/CH2/EX2.8/2_8.sce | 41 +++++++++++++++++++++++++++++++++++++++++ 1076/CH2/EX2.9/2_9.sce | 13 +++++++++++++ 20 files changed, 346 insertions(+) create mode 100755 1076/CH2/EX2.1/2_1.sce create mode 100755 1076/CH2/EX2.10/2_10.sce create mode 100755 1076/CH2/EX2.11/2_11.sce create mode 100755 1076/CH2/EX2.12/2_12.sce create mode 100755 1076/CH2/EX2.13/2_13.sce create mode 100755 1076/CH2/EX2.14/2_14.sce create mode 100755 1076/CH2/EX2.15/2_15.sce create mode 100755 1076/CH2/EX2.16/2_16.sce create mode 100755 1076/CH2/EX2.17/2_17.sce create mode 100755 1076/CH2/EX2.18/2_18.sce create mode 100755 1076/CH2/EX2.19/2_19.sce create mode 100755 1076/CH2/EX2.20/2_20.sce create mode 100755 1076/CH2/EX2.21/2_21.sce create mode 100755 1076/CH2/EX2.3/2_3.sce create mode 100755 1076/CH2/EX2.4/2_4.sce create mode 100755 1076/CH2/EX2.5/2_5.sce create mode 100755 1076/CH2/EX2.6/2_6.sce create mode 100755 1076/CH2/EX2.7/2_7.sce create mode 100755 1076/CH2/EX2.8/2_8.sce create mode 100755 1076/CH2/EX2.9/2_9.sce (limited to '1076/CH2') 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); -- cgit