From 6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d Mon Sep 17 00:00:00 2001 From: kinitrupti Date: Fri, 12 May 2017 18:53:46 +0530 Subject: Removed duplicates --- .../Chapter_4.ipynb | 1168 ++++++++++++++++++++ 1 file changed, 1168 insertions(+) create mode 100755 Elements_of_Power_system_by_J.B._Gupta/Chapter_4.ipynb (limited to 'Elements_of_Power_system_by_J.B._Gupta/Chapter_4.ipynb') diff --git a/Elements_of_Power_system_by_J.B._Gupta/Chapter_4.ipynb b/Elements_of_Power_system_by_J.B._Gupta/Chapter_4.ipynb new file mode 100755 index 00000000..db67cdfa --- /dev/null +++ b/Elements_of_Power_system_by_J.B._Gupta/Chapter_4.ipynb @@ -0,0 +1,1168 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:818ef23f34385ee54f8d5672494c60336bc18fcc46f244fceb4d0fa90c810f66" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 4 - INDUCTANCE AND CAPACITANCE OF TRANSMISSION LINES" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E1 - Pg 85" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Loop inductance,Reactance of transmission line\n", + "#Given data :\n", + "import math\n", + "f=50.##Hz\n", + "d=1.*100.##cm\n", + "r=1.25/2.##cm\n", + "r_dash=r*0.7788##cm\n", + "L=0.4*math.log(d/r_dash)##mH\n", + "print '%s %.2f' %(\"Loop inductance per km(mH)\",L)#\n", + "XL=2.*math.pi*f*L*10.**-3.##ohm/Km\n", + "print '%s %.3f' %(\"Reactance of transmission line(ohm/km)\",XL)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Loop inductance per km(mH) 2.13\n", + "Reactance of transmission line(ohm/km) 0.669\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E2 - Pg 85" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Loop inductance\n", + "#Given data :\n", + "import math\n", + "f=50.##Hz\n", + "a=10.##cm**2\n", + "l=500./1000.##km\n", + "r=math.sqrt(a/math.pi)##cm\n", + "d=5.*100.##cm\n", + "r_dash=r*0.7788##cm\n", + "L=0.4*math.log(d/r_dash)*l##mH\n", + "print '%s %.3f' %(\"Loop inductance per km(mH)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Loop inductance per km(mH) 1.177\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E3 - Pg 85" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Loop inductance per km of copper conductor line,Loop inductance per km of steel conductor line\n", + "#Given data :\n", + "import math\n", + "r=1./2.##cm\n", + "d=1.*100.##cm\n", + "mu=50.##relative permeability\n", + "r_dash=r*0.7788##cm\n", + "L_cu=.1+0.4*math.log(d/r)##mH\n", + "print '%s %.2f' %(\"Loop inductance per km of copper conductor line(mH)\",L_cu)#\n", + "L_steel=(mu+4.*math.log(d/r))*10.**-7.*10.**3.##mH\n", + "print '%s %.2f' %(\"Loop inductance per km of steel conductor line(mH)\",L_steel*10**3)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Loop inductance per km of copper conductor line(mH) 2.22\n", + "Loop inductance per km of steel conductor line(mH) 7.12\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E4 - Pg 86" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Geometric mean radius\n", + "#Given data :\n", + "import math\n", + "r=3.##mm\n", + "d11=r##mm\n", + "d12=2.*r##mm\n", + "d34=2.*r##mm\n", + "d16=2.*r##mm\n", + "d17=2.*r##mm\n", + "d14=4.*r##mm\n", + "d13=math.sqrt(d14**2.-d34**2.)##mm\n", + "d15=d13##mm\n", + "Ds1=(0.7788*d11*d12*d13*d14*d15*d16*d17)**(1./7.)##mm\n", + "Ds2=Ds1##mm\n", + "Ds3=Ds1##mm\n", + "Ds4=Ds1##mm\n", + "Ds5=Ds1##mm\n", + "Ds6=Ds1##mm\n", + "Ds7=(2.*r*0.7788*d11*d12*d13*2.*r*2.*r)**(1./7.)##mm\n", + "Ds=(Ds1*Ds2*Ds3*Ds4*Ds5*Ds6*Ds7)**(1./7.)##mm\n", + "print '%s %.3f' %(\"Geometric mean radius(mm)\",Ds)#\n", + "#Answer in the book is wrong\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Geometric mean radius(mm) 6.367\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E5 - Pg 86" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Loop inductance of line\n", + "#Given data :\n", + "import math\n", + "r=1.2##cm\n", + "rdash=0.7788*r##cm\n", + "d12=0.12*100.##cm\n", + "d11dash=(0.2+1.2)*100.##cm\n", + "d22dash=(0.2+1.2)*100.##cm\n", + "d12dash=(0.2+1.2+0.2)*100.##cm\n", + "d21dash=(1.2)*100.##cm\n", + "Dm=(d11dash*d12dash*d21dash*d22dash)**(1./4.)##cm\n", + "d11=0.93456##cm\n", + "d22=0.93456##cm\n", + "d12=20.##cm\n", + "d21=20.##cm\n", + "Ds=(d11*d12*d21*d22)**(1./4.)##cm\n", + "L=0.4*math.log(Dm/Ds)##mH/km\n", + "print '%s %.3f' %(\"Loop inductance of line(mH/km)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Loop inductance of line(mH/km) 1.389\n" + ] + } + ], + "prompt_number": 5 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E6 - Pg 87" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Loop inductance of line\n", + "#Given data :\n", + "import math\n", + "r=2./2.##cm\n", + "rdash=0.7788*r##cm\n", + "d12=0.12*100##cm\n", + "d11dash=300.##cm\n", + "d12dash=math.sqrt(300.**2.+100.**2.)##cm\n", + "d21dash=d12dash##cm\n", + "d22dash=d11dash##cm\n", + "d11=rdash##cm\n", + "d22=rdash##cm\n", + "d12=100.##cm\n", + "d21=100.##cm\n", + "Dm=(d11dash*d12dash*d21dash*d22dash)**(1./4.)##cm\n", + "Ds=(d11*d12*d21*d22)**(1./4.)##cm\n", + "L=0.4*math.log(Dm/Ds)##mH/km\n", + "print '%s %.3f' %(\"Loop inductance of line(mH/km)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Loop inductance of line(mH/km) 1.421\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + " Example E7 - Pg 89" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Inductance per phase \n", + "#Given data :\n", + "import math\n", + "r=1.24/2##cm\n", + "rdash=0.7788*r##cm\n", + "d=2.*100.##cm\n", + "L=0.2*math.log(d/rdash)##mH\n", + "print '%s %.3f' %(\"Inductance per phase per km(mH)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Inductance per phase per km(mH) 1.205\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E8 - Pg 89" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Inductance per phase\n", + "#Given data :\n", + "import math\n", + "r=(20./2.)/10.##cm\n", + "d1=4.*100.##cm\n", + "d2=5.*100.##cm\n", + "d3=6.*100.##cm\n", + "rdash=0.7788*r##cm\n", + "L=0.2*math.log((d1*d2*d3)**(1./3.)/rdash)##mH\n", + "print '%s %.2f' %(\"Inductance per phase(mH)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Inductance per phase(mH) 1.29\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E9 - Pg 90" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Inductance per km of phase1,phase2,phase3\n", + "#Given data :\n", + "import math\n", + "r=4./2.##cm\n", + "rdash=0.7788*r##cm\n", + "d=300.##cm\n", + "d3=6.*100.##cm\n", + "LAr=0.2*(math.log(d/rdash)+1./2.*math.log(2.))\n", + "LAi=0.866*math.log(2.)##mH\n", + "print '%s' %(\"Inductance per km of phase1(mH)\")#\n", + "print '%.4f %s %.2f' %(LAr,'-j',LAi)\n", + "LB=0.2*math.log(d/rdash)##mH\n", + "print '%s %.3f' %(\"Inductance per km of phase2(mH)\",LB)#\n", + "LCr=0.2*(math.log(d/rdash)+1./2.*math.log(2.))\n", + "LCi=0.866*math.log(2.)##mH\n", + "print '%s' %(\"Inductance per km of phase3(mH)\")#\n", + "print '%.4f %s %.2f' %(LCr,'+j',LCi)\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Inductance per km of phase1(mH)\n", + "1.1214 -j 0.60\n", + "Inductance per km of phase2(mH) 1.052\n", + "Inductance per km of phase3(mH)\n", + "1.1214 +j 0.60\n" + ] + } + ], + "prompt_number": 9 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E10 - Pg 90" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Spacing between adjacent conductors\n", + "import math\n", + "#Given data :\n", + "r=1.2/2.*10.##mm\n", + "rdash=0.7788*r##mm\n", + "d=3.5*1000.##mm\n", + "L=2.*10.**-7.*math.log(d/rdash)##H/m\n", + "Lav=1./3.*(L+L+L)##H/m\n", + "d=rdash*math.exp(Lav/(2.*10.**-7.)-1./3.*math.log(2.))##mm\n", + "print '%s %.4e' %(\"Spacing between adjacent conductors(m)\",d/1000.)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Spacing between adjacent conductors(m) 2.7780e+00\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E11 - Pg 94" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Spacing between adjacent conductors\n", + "import math\n", + "#Given data :\n", + "r=20.##mm\n", + "rdash=0.7788*r##mm\n", + "d=7.*1000.##mm\n", + "L=10**-7*math.log(math.sqrt(3.)/2.*d/rdash)##H/m\n", + "print '%s %.4f' %(\"Spacing between adjacent conductors(mH)\",L*10.**3./10.**-3.)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Spacing between adjacent conductors(mH) 0.5964\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E12 - Pg 94" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Inductance per phase\n", + "import math \n", + "#Given data :\n", + "r=0.9##cm\n", + "rdash=0.7788*r*10.**-2.##m\n", + "daa_dash=math.sqrt(6.**2.+6.**2.)##m\n", + "dbb_dash=7.##m\n", + "dcc_dash=daa_dash##m\n", + "daa=rdash##m\n", + "d_adash_adash=rdash##m\n", + "d_adash_a=daa_dash##m\n", + "Dsa=(daa*daa_dash*d_adash_adash*d_adash_a)**(1./4.)##m\n", + "Dsb=(daa*7.)**(1./2.)##m\n", + "Dsc=(daa*daa_dash)**(1./2.)##m\n", + "Ds=(Dsa*Dsb*Dsc)**(1./3.)##m\n", + "dab=math.sqrt(3.**2.+0.5**2.)##m\n", + "dab_dash=math.sqrt(3.**2.+6.5**2.)##m\n", + "d_adash_b=math.sqrt(3.**2.+6.5**2.)##m\n", + "d_adash_bdash=math.sqrt(3.**2.+0.5**2.)##m\n", + "Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)**(1./4.)##m\n", + "Dbc=((dab*dab_dash)**2.)**(1./4.)##m\n", + "Dca=((6.*6.)**2.)**(1./4.)##m\n", + "Dm=(Dab*Dbc*Dca)**(1./3.)##m\n", + "L=0.2*math.log(Dm/Ds)##mH/km\n", + "print '%s %.4f' %(\"Inductance per phase(mH/km)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Inductance per phase(mH/km) 0.6135\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E13 - Pg 95" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate GMD,Deq or Dm,Inductance of 100 km line,Using Alternate method, Inductance of 100 km line\n", + "import math\n", + "#Given data :\n", + "r=5./2.##mm\n", + "rdash=2.176*r*10.**-3.##m\n", + "daa_dash=math.sqrt(6**2+16**2)##m\n", + "dbb_dash=6.##m\n", + "dcc_dash=daa_dash##m\n", + "dab=8.##m\n", + "dab_dash=math.sqrt(6.**2.+8.**2.)##m\n", + "dbc=8.##m\n", + "dbc_dash=math.sqrt(6.**2.+8.**2.)##m\n", + "dca=16.##m\n", + "dca_dash=6.##m\n", + "Dsa=math.sqrt(rdash*daa_dash)##m\n", + "Dsb=math.sqrt(rdash*dbb_dash)##m\n", + "Dsc=math.sqrt(rdash*dcc_dash)##m\n", + "Ds=(Dsa*Dsb*Dsc)**(1/3)##m\n", + "print '%s %.2f' %(\"GMD(m) : \",Ds)#\n", + "Dab=(dab*dab_dash)**(1./2.)##m\n", + "Dbc=(dbc*dbc_dash)**(1./2.)##m\n", + "Dca=(dca*dca_dash)**(1./2.)##m\n", + "Dm=(Dab*Dbc*Dca)**(1./3.)##m\n", + "print '%s %.2f' %(\"Deq or Dm(m) : \",Dm)#\n", + "L=0.2*math.log(Dm/Ds)##mH/km\n", + "L=L*10.**-3.*100.##H(for 100 km line)\n", + "print '%s %.4f' %(\"Inductance of 100 km line(H)\",L)#\n", + "#/Alternate method is given below\n", + "d1=dab##m\n", + "d2=dca_dash##m\n", + "L=0.2*math.log(2.**(1./6.))*math.sqrt(d1/rdash)*((d1**2.+d2**2.)/(4.*d1**2.+d2**2.))**(1./6.)##mH\n", + "L=L*10.**-3.*100.##H(for 100 km line)\n", + "print '%s %.4f' %(\"Using Alternate method, Inductance of 100 km line(H)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "GMD(m) : 1.00\n", + "Deq or Dm(m) : 9.22\n", + "Inductance of 100 km line(H) 0.0444\n", + "Using Alternate method, Inductance of 100 km line(H) 0.0741\n" + ] + } + ], + "prompt_number": 13 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E14 - Pg 97" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Inductance per phase\n", + "import math\n", + "#Given data :\n", + "r=5./2.##cm\n", + "rdash=0.7788*r*10.**-2.##m\n", + "d=6.5##m\n", + "s=0.4##m\n", + "Ds=math.sqrt(rdash*s)##m\n", + "dab=6.5##m\n", + "dab_dash=6.9##m\n", + "d_adash_b=6.1##m\n", + "d_adash_bdash=6.5##m\n", + "Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)**(1./4.)##m\n", + "Dbc=Dab##m\n", + "dca=13.##m\n", + "dca_dash=12.6##m\n", + "d_cdash_a=13.4##m\n", + "d_cdash_adash=13.##m\n", + "Dca=(dca*dca_dash*d_cdash_a*d_cdash_adash)**(1./4.)##m\n", + "Dm=(Dab*Dbc*Dca)**(1./3.)##m\n", + "L=0.2*math.log(Dm/Ds)##mH/km\n", + "print '%s %.3f' %(\"Inductance per phase(mH/km)\",L)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Inductance per phase(mH/km) 0.906\n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E15 - Pg 98" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Inductive reactance of bundled conductor line,Inductive reactance with single conductor\n", + "import math\n", + "#Given data :\n", + "f=50.##Hz\n", + "r=3.5/2.##cm\n", + "rdash=0.7788*r*10.**-2.##m\n", + "d=7.##m\n", + "s=40./100.##m\n", + "Ds=math.sqrt(rdash*s)##m\n", + "dab=7.##m\n", + "dab_dash=7.4##m\n", + "d_adash_b=6.6##m\n", + "d_adash_bdash=7.##m\n", + "Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)**(1./4.)##m\n", + "Dbc=Dab##m\n", + "dca=14.##m\n", + "dca_dash=13.6##m\n", + "d_cdash_a=14.4##m\n", + "d_cdash_adash=14.##m\n", + "Dca=(dca*dca_dash*d_cdash_a*d_cdash_adash)**(1./4.)##m\n", + "Dm=(Dab*Dbc*Dca)**(1./3.)##m\n", + "L=0.2*math.log(Dm/Ds)##mH/km\n", + "XL=2*math.pi*f*L*10.**-3.##ohm/km\n", + "print '%s %.2f' %(\"Inductive reactance of bundled conductor line(ohm/km)\",XL)#\n", + "#Equivalent single conductor\n", + "n=2#\n", + "r1=math.sqrt(n*math.pi*r**2./math.pi)##m\n", + "r1dash=0.7788*r1*10.**-2.##m\n", + "Dm1=(Dab*Dbc*Dca)**(1./3.)##m\n", + "L1=0.2*math.log(Dm1/r1dash)##mH/km\n", + "XL1=2*math.pi*f*L1*10.**-3.##ohm/km\n", + "print '%s %.3f' %(\"Inductive reactance with single conductor(ohm/km)\",XL1)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Inductive reactance of bundled conductor line(ohm/km) 0.30\n", + "Inductive reactance with single conductor(ohm/km) 0.385\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E16 - Pg 102" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance of line\n", + "import math\n", + "#Given data :\n", + "r=15./2.##mm\n", + "d=1.5*1000.##mm\n", + "l=30.##km\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "C=math.pi*epsilon_o/math.log(d/r)*l*1000.##F\n", + "print '%s %.4f' %(\"Capacitance of line(micro F)\",C*10**6)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance of line(micro F) 0.1575\n" + ] + } + ], + "prompt_number": 16 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E17 - Pg 105" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance of line\n", + "#Given data :\n", + "import math\n", + "r=2./2.##cm\n", + "d=2.5*100.##cm\n", + "l=100.##km\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "C=2*math.pi*epsilon_o/math.log(d/r)*l*1000.##F\n", + "print '%s %.4f' %(\"Capacitance of line(micro F)\",C*10.**6.)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance of line(micro F) 1.0075\n" + ] + } + ], + "prompt_number": 17 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E18 - Pg 105" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance of line\n", + "import math\n", + "#Given data :\n", + "r=2./2./100.##m\n", + "d1=3.5##m\n", + "d2=5.##m\n", + "d3=8.##m\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "CN=2.*math.pi*epsilon_o*1000/math.log((d1*d2*d3)**(1./3.)/r)##F\n", + "print '%s %.4f' %(\"Capacitance of line(micro F)\",CN*10.**6.)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance of line(micro F) 0.0089\n" + ] + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E19 - Pg 105" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance per phase per meter line,Charging current per phase\n", + "import math\n", + "#Given data :\n", + "f=50.##Hz\n", + "VL=220.##KV\n", + "r=20./2./1000.##m\n", + "d1=3.##m\n", + "d2=3.##m\n", + "d3=6.##m\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "CN=2.*math.pi*epsilon_o/math.log((d1*d2*d3)**(1./3.)/r)##F\n", + "print '%s %.4e' %(\"Capacitance per phase per meter line(F)\",CN)#\n", + "Vph=VL*1000./math.sqrt(3.)##V\n", + "Ic=2.*math.pi*f*CN*Vph##A\n", + "print '%s %.3f' %(\"Charging current per phase(mA) : \",Ic*1000.)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance per phase per meter line(F) 9.3737e-12\n", + "Charging current per phase(mA) : 0.374\n" + ] + } + ], + "prompt_number": 19 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E20 - Pg 106" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance per phase per meter line,Charging current per phase\n", + "import math\n", + "#Given data :\n", + "f=50.##Hz\n", + "VL=110.##kV\n", + "r=1.05/2.##cm\n", + "d1=3.5##m\n", + "d2=3.5##m\n", + "d3=7.##m\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "CN=2.*math.pi*epsilon_o/math.log((d1*d2*d3)**(1./3.)*100./r)##F\n", + "print '%s %.2e' %(\"Capacitance per phase per meter line(F)\",CN)#\n", + "Vph=VL*1000./math.sqrt(3.)##V\n", + "Ic=2.*math.pi*f*CN*Vph##A/m\n", + "print '%s %.3f' %(\"Charging current per phase(A/km) : \",Ic/10**-3)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance per phase per meter line(F) 8.26e-12\n", + "Charging current per phase(A/km) : 0.165\n" + ] + } + ], + "prompt_number": 20 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E21 - Pg 108" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitive reactance too neutral,Charging current\n", + "import math\n", + "#Given data :\n", + "r=2.5/2.*10.**-2.##m\n", + "VL=132.##KV\n", + "epsilon_o=8.85*10.**-12.##permitivity\n", + "f=50.##Hz\n", + "dRRdash=math.sqrt(7.**2.+(4.+4.)**2.)##m\n", + "dBBdash=dRRdash##m\n", + "dYYdash=9.##m\n", + "DSR=math.sqrt(r*dRRdash)##m\n", + "DSY=math.sqrt(r*dYYdash)##m\n", + "DSB=math.sqrt(r*dBBdash)##m\n", + "Ds=(DSR*DSB*DSY)**(1./3.)##m\n", + "dRY=math.sqrt(4.**2.+(4.5-3.5)**2.)##m\n", + "dRYdash=math.sqrt((9.-1.)**2.+4.**2.)##m\n", + "dRdashY=math.sqrt((9.-1.)**2.+4.**2.)##m\n", + "dRdashYdash=math.sqrt(4.**2.+(4.5-3.5)**2.)##m\n", + "DRY=(dRY*dRYdash*dRdashY*dRdashYdash)**(1./4.)##m\n", + "DYB=((dRY*dRYdash)**2.)**(1./4.)##m\n", + "DBR=((8.*7.)**2.)**(1./4.)##m\n", + "Dm=(DRY*DYB*DBR)**(1./3.)##m\n", + "C=2*math.pi*epsilon_o/math.log(Dm/Ds)##F/m\n", + "C=C/10.**-3.##F/km\n", + "X=1./(2.*math.pi*f*C)##ohm\n", + "print '%s %.3f' %(\"Capacitive reactance too neutral(kohm) : \",X/1000.)#\n", + "Vph=VL*1000./math.sqrt(3.)##Volt\n", + "Ic=2.*math.pi*f*C*Vph##A\n", + "print '%s %.4f' %(\"Charging current(A/km)\",Ic)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitive reactance too neutral(kohm) : 166.599\n", + "Charging current(A/km) 0.4574\n" + ] + } + ], + "prompt_number": 21 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E22 - Pg 109" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance of 100 km line\n", + "import math\n", + "#Given data :\n", + "d1=8.##m\n", + "d2=6.##m\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "r=3.*5./2.*10.**-3.##m\n", + "C=4*math.pi*epsilon_o/math.log(2.**(1./3.)*d1/r*((d1**2.+d2**2.)/(4.*d1**2.+d2**2.)**(1./3.)))##F/m\n", + "C100=C*100.*1000.*10.**6.##microF\n", + "print '%s %.3f' %(\"Capacitance of 100 km line(micro Farad) : \",C100)#\n", + "#answer in the textbook is wrong.\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance of 100 km line(micro Farad) : 1.122\n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E23 - Pg 110" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance,Charging currentimport math\n", + "#Given data :\n", + "import math\n", + "VL=132.##kV\n", + "f=50.##Hz\n", + "r=5./2.##cm\n", + "rdash=0.7788*r*10.**-2.##m\n", + "d=6.5##m\n", + "s=0.4##m\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "Ds=math.sqrt(rdash*s)##m\n", + "dab=6.5##m\n", + "dab_dash=6.9##m\n", + "d_adash_b=6.1##m\n", + "d_adash_bdash=6.5##m\n", + "Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)**(1./4.)##m\n", + "Dbc=Dab##m\n", + "dca=13.##m\n", + "dca_dash=12.6##m\n", + "d_cdash_a=13.4##m\n", + "d_cdash_adash=13.##m\n", + "Dca=(dca*dca_dash*d_cdash_a*d_cdash_adash)**(1./4.)##m\n", + "Dm=(Dab*Dbc*Dca)**(1./3.)##m\n", + "L=0.2*math.log(Dm/Ds)##mH/km\n", + "C=2*math.pi*epsilon_o/math.log(Dm/Ds)##F/m\n", + "C=C/10.**-3.##F/km\n", + "print '%s %.2e' %(\"Capacitance per km(F/km) : \",C)#\n", + "Vph=VL*1000./math.sqrt(3)##Volt\n", + "Ic=2*math.pi*f*C*Vph##A/km\n", + "print '%s %.1f' %(\"Charging current per km(A/km) : \",Ic)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance per km(F/km) : 1.23e-08\n", + "Charging current per km(A/km) : 0.3\n" + ] + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E24 - Pg 112" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate One conductor ACSR moose conductor line,Inductive reactance per Km per phase(ohm) :,Capacitivetive reactance per Km per phase(Mohm),Three conductor bundled line\n", + "#calculate Inductive reactance per km per phase(ohm)\n", + "#calculate Capacitivetive reactance per km per phase(Mohm)\n", + "import math\n", + "#Given data :\n", + "VL=132.##kV\n", + "f=50.##Hz\n", + "r=31.8/2##mm\n", + "rdash=0.7788*r##mm\n", + "d=10.*1000.##mm\n", + "epsilon_o=8.854*10**-12##permitivity\n", + "print '%s' %(\"One conductor ACSR moose conductor line : \")#\n", + "LAr=0.2*(math.log(d/rdash)+1./2.*math.log(2.))\n", + "LAi=0.866*math.log(2.)##mH/km\n", + "LB=0.2*math.log(d/rdash)##mH/km\n", + "LCr=0.2*(math.log(d/rdash)+1./2.*math.log(2.))\n", + "LCi=0.866*math.log(2.)##mH/km\n", + "Lav=(LAr+LAi+LB+LCr+LCi)/3.##mH/km\n", + "XL=2*math.pi*f*Lav*10.**-3.##ohm\n", + "print '%s %.3f' %(\"Inductive reactance per Km per phase(ohm) : \",XL)#\n", + "d1=10.##m\n", + "d2=10.##m\n", + "d3=20.##m\n", + "CN=2.*math.pi*epsilon_o/math.log((d1*d2*d3)**(1./3.)/(rdash*10.**-3.))/10.**3.##F/km\n", + "XC=1./(2.*math.pi*f*CN*10.**6.)##ohm\n", + "print '%s %.3f' %(\"Capacitivetive reactance per Km per phase(Mohm) : \",XC/10**6)#\n", + "print '%s' %(\"Three conductor bundled line : \")#\n", + "S=40./100.##m\n", + "Ds=(rdash*10.**-3.*S**2.)**(1./3.)##m\n", + "Deq=(d1*d2*d3)**(1./3.)##m\n", + "Ldash=0.2*math.log(Deq/Ds)##mH/km\n", + "XLdash=2*math.pi*f*Ldash*10.**-3.##ohm\n", + "print '%s %.3f' %(\"Inductive reactance per km per phase(ohm) : \",XLdash)#\n", + "Ds=(r*10.**-3.*S**2.)**(1./3.)##m\n", + "Cdash=2.*math.pi*epsilon_o*10.**3./math.log(Deq/Ds)##microF/km\n", + "XC=1./(2.*math.pi*f*Cdash)/10.**6.##Mohm\n", + "print '%s %.3f' %(\"Capacitivetive reactance per km per phase(Mohm) : \",XC)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "One conductor ACSR moose conductor line : \n", + "Inductive reactance per Km per phase(ohm) : 0.561\n", + "Capacitivetive reactance per Km per phase(Mohm) : 0.396\n", + "Three conductor bundled line : \n", + "Inductive reactance per km per phase(ohm) : 0.290\n", + "Capacitivetive reactance per km per phase(Mohm) : 0.259\n" + ] + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E25 - Pg 114" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Capacitance per km of line\n", + "import math\n", + "#Given data :\n", + "r=1.5/2.##cm\n", + "d=3.*100.##cm\n", + "h=6.*100.##cm\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "C=math.pi*epsilon_o/math.log(d/(1.+d**2./4./h**2.)**r)*10.**3.##F\n", + "print '%s %.3e' %(\"Capacitance per km of line(F) : \",C)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance per km of line(F) : 4.916e-09\n" + ] + } + ], + "prompt_number": 25 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example E26 - Pg 114" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "#calculate Part(i) Capacitance per phase per meter length,Part(ii) Capacitance per phase per meter length\n", + "import math\n", + "#Given data :\n", + "r=2./100.##m\n", + "d1=4.##m\n", + "d2=4.##m\n", + "d3=8.##m\n", + "epsilon_o=8.854*10.**-12.##permitivity\n", + "CN=2.*math.pi*epsilon_o/math.log((d1*d2*d3)**(1./3.)/r)##F\n", + "print '%s %.3e' %(\"Part(i) Capacitance per phase per meter length(F) : \",CN)#\n", + "h1=20.##m\n", + "h2=20.##m\n", + "h3=20.##m\n", + "h12=math.sqrt(20.**2.+4.**2.)##m\n", + "h23=math.sqrt(20.**2.+4.**2.)##m\n", + "h31=math.sqrt(20.**2.+8.**2.)##m\n", + "Deq=(d1*d2*d3)**(1./3.)##m\n", + "CN=2.*math.pi*epsilon_o/(math.log(Deq/r)-math.log((h12*h23*h31/h1/h2/h3)**(1./3.)) )##F\n", + "print '%s %.3e' %(\"Part(ii) Capacitance per phase per meter length(F) : \",CN)#\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Part(i) Capacitance per phase per meter length(F) : 1.006e-11\n", + "Part(ii) Capacitance per phase per meter length(F) : 1.013e-11\n" + ] + } + ], + "prompt_number": 26 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file -- cgit