From 79c59acc7af08ede23167b8455de4b716f77601f Mon Sep 17 00:00:00 2001 From: hardythe1 Date: Thu, 11 Jun 2015 17:31:11 +0530 Subject: add books --- .../Chapter18.ipynb | 326 +++++++++++++++++++++ 1 file changed, 326 insertions(+) create mode 100755 Electrical_Power_System_by_C.L._Wadhwa/Chapter18.ipynb (limited to 'Electrical_Power_System_by_C.L._Wadhwa/Chapter18.ipynb') diff --git a/Electrical_Power_System_by_C.L._Wadhwa/Chapter18.ipynb b/Electrical_Power_System_by_C.L._Wadhwa/Chapter18.ipynb new file mode 100755 index 00000000..379c3e7a --- /dev/null +++ b/Electrical_Power_System_by_C.L._Wadhwa/Chapter18.ipynb @@ -0,0 +1,326 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 18 : Load Flows" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 18.1, Page No 605" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "import numpy\n", + "#initialisation of variables\n", + "Y=[[complex(3,-12),complex(-2,8),complex(-1,4),0],\n", + " [complex(-2,8), complex(3.666,-14.664), complex(-0.666,2.6664),complex(-1,4)],\n", + " [complex(-1,4), complex(-0.666,2.6664), complex(3.666,-14.664), complex(-2,8)],\n", + " [0,complex(-1,4), complex(-2,8),complex(3,-12)]]\n", + "P2=-.5\n", + "P3=-.4\n", + "P4=-.3\n", + "Q4=-.1\n", + "Q3=-.3\n", + "Q2=-.2\n", + "V2=1.0\n", + "V3=1.0\n", + "V4=1.0\n", + "V10=1.06\n", + "V30=1.0\n", + "V40=1\n", + "\n", + "#Calculations\n", + "V21=(((complex(P2,-Q2))/V2)-Y[2][1]*V10-Y[1][2]*V30-Y[1][3]*V40)/(Y[1][1])\n", + "V21acc=1+1.6*(V21-1)\n", + "print(\"V21acc= {0:.5f}{1:.5f}i\".format(V21acc.real, V21acc.imag))\n", + "V31=((complex(P3,-Q3)/V3)-Y[2][0]*V10-Y[2][1]*V21acc-Y[2][3]*V40)/(Y[2][2])\n", + "V31acc=1+1.6*(V31-1)\n", + "print(\"V31acc= {0:.5f}{1:.5f}i\".format(V31acc.real, V31acc.imag))\n", + "V41=((complex(P4,-Q4)/V4)-Y[3][1]*V21acc-Y[3][2]*V31acc)/(Y[3][3])\n", + "V41acc=1+1.6*(V41-1)\n", + "\n", + "#Results\n", + "print(\"V41acc= {0:.5f}{1:.5f}i\".format(V41acc.real, V41acc.imag))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "V21acc= 0.40236-0.04634i\n", + "V31acc= 0.81149-0.04688i\n", + "V41acc= 0.45822-0.10923i\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 18.2, Page No 606" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#initialisation of variables\n", + "\n", + "Y=[[complex(3,-12),complex(-2,8),complex(-1,4),0],[complex(-2,8),complex(3.666,14.664),complex(-0.666,2.6664),complex(-1,4)],[complex(-1,4), complex(-.666,2.6664),complex(3.666,-14.664),complex(-2,8)],[0,complex(-1,4),complex(-2,8), complex(3,-12)]]\n", + "P2=.5\n", + "P3=-.4\n", + "P4=-.3\n", + "Q4=-.1\n", + "Q3=-.3\n", + "V3=1\n", + "V4=1\n", + "V1=1.06\n", + "V2=1.04\n", + "V30=1\n", + "V40=1\n", + "\n", + "#Calculations\n", + "Q2=-((V2*(Y[1][0]*V1+Y[1][1]*V2+Y[1][2]*V3+Y[1][3]*V4))).imag\n", + "V21=((complex(P2,-Q2)/V2)-Y[1][0]*V1-Y[1][2]*V30-Y[1][3]*V40)/(Y[1][1])\n", + "d=math.degrees(math.atan(0.0291473/1.0472868))\n", + "V21=1.04*complex(math.cos(math.radians(d)),math.sin(math.radians(d)))\n", + "\n", + "#Results\n", + "print(\"The value of V21 is= {0:.5f}+{1:.5f}i\".format(V21.real, V21.imag))\n", + "V31=((complex(P3,-Q3)/V3)-Y[2][0]*V1-Y[2][1]*V21-Y[2][3]*V40)/(Y[2][2])\n", + "print(\"The value of V31 is= {0:.5f}+{1:.5f}i\".format(V31.real, V31.imag))\n", + "V41=((complex(P4,-Q4)/V4)-Y[3][1]*V21-Y[3][2]*V31)/Y[3][3]\n", + "print(\"The value of V41 is= {0:.5f}+{1:.5f}i\".format(V41.real, V41.imag))\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The value of V21 is= 1.03960+0.02893i\n", + "The value of V31 is= 0.99805+-0.01564i\n", + "The value of V41 is= 0.99817+-0.02235i\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 18.3, Page No 607" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#initialisation of variables\n", + "\n", + "Y=[[complex(3,-12),complex(-2,8),complex(-1,4),0],[complex(-2,8),complex(3.666,14.664),complex(-0.666,2.6664),complex(-1,4)],[complex(-1,4), complex(-.666,2.6664),complex(3.666,-14.664),complex(-2,8)],[0,complex(-1,4),complex(-2,8), complex(3,-12)]]\n", + "P2=.5\n", + "P3=-.4\n", + "P4=-.3\n", + "Q4=-.1\n", + "Q3=-.3\n", + "V3=1\n", + "V4=1\n", + "V1=1.06\n", + "V2=1\n", + "V30=1\n", + "V40=1\n", + "Q2=.2\n", + "V3=1\n", + "\n", + "#Calculations\n", + "V21=((complex(P2,-Q2)/V2)-Y[1][0]*V1-Y[1][2]*V30-Y[1][3]*V40)/(Y[1][1])\n", + "V31=((complex(P3,-Q3)/V3)-Y[2][0]*V1-Y[2][1]*V21-Y[2][3]*V40)/(Y[2][2])\n", + "V41=((complex(P4,-Q4)/V4)-Y[3][1]*V21-Y[3][2]*V31)/Y[3][3]\n", + "\n", + "#Results\n", + "print(\"The value of V21 is= {0:.5f}+{1:.5f}i\".format(V21.real, V21.imag))\n", + "print(\"The value of V31 is= {0:.5f}+{1:.5f}i\".format(V31.real, V31.imag))\n", + "print(\"The value of V41 is= {0:.5f}+{1:.5f}i\".format(V41.real, V41.imag))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The value of V21 is= -0.91620+-0.52133i\n", + "The value of V31 is= 0.64242+-0.11561i\n", + "The value of V41 is= 0.10915+-0.27242i\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 18.4 Page No 615" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#initialisation of variables\n", + "V1=1.06\n", + "G11=6.25\n", + "G12=-1.25\n", + "G21=G12\n", + "G13=-5\n", + "G31=G13\n", + "G22=2.916\n", + "G23=-1.666\n", + "G32=G23\n", + "G33=6.666\n", + "B11=18.75\n", + "B12=-3.75\n", + "B21=B12\n", + "B13=-15\n", + "B31=B13\n", + "B22=8.75\n", + "B23=-5\n", + "B32=B23\n", + "B33=20\n", + "e1=1.06\n", + "e2=1\n", + "e3=1\n", + "f1=0\n", + "f2=0\n", + "f3=0\n", + "\n", + "#Calculations\n", + "P2=e2*(e1*G21+f1*B21) +f2*(f1*G21-e1*B21) +e2*(e2*G22+f2*B22)+f2*(f2*G22-e2*B22)+e2*(e3*G23+f3*B23)+f2*(f3*G23-e3*B23)\n", + "P3=-.3\n", + "Q2=-.225\n", + "Q3=-.9\n", + "dP2=.2-(-.225)\n", + "dP3=-.6-(-.3)\n", + "dQ2=0-(-.225)\n", + "dQ3=-.25-(-.9)\n", + "a1=2*e2*G22+e1*G21+f1*B21+e3*G23+f3*B23#a1=dP2/de2\n", + "a2=2*e3*G33+e1*G31+f1*B31+e3*G32+f2*B32#a2=dP3/de3\n", + "b1=2*f2*G22 +f1*G21-e1*B21+f3*G23-e3*B23#b1=dP2/df2\n", + "b2=20.9#dP3/df3\n", + "a3=e2*G23-f2*B23#dP2/de3\n", + "a4=-1.666\t\t#dP3/de2\n", + "b3=-5\t\t\t#dP2/df3\n", + "b4=-5\t\t\t#dP3/df2\n", + "c1=2*e2*B22-f1*G21+e1*B21-f3*G23+e3*B23#dQ2/de2\n", + "c2=19.1\t\t\t#dQ3/de3\n", + "c3=-2.991\t\t#dQ2/df2\n", + "c4=-6.966\t\t#dQ3/df3\n", + "\n", + "#Results\n", + "print(\"set of linear equations at the end of first iteration are\")\n", + "print(\"%.3fde2 %.3fde3+ %.3fdf2 %.3fdf3 = %.3f\" %(2.846,-1.666,8.975,-5,2.75))\n", + "print(\"%.3fde2 +%.3fde3 %.3fdf2 +%.3fdf3 = %.3f\" %(-1.666,6.366,-5,20.90,-.3))\n", + "print(\"%.3fde2 %.3fde3 %.3fdf2 +%.3fdf3 = %.3f\" %(8.525,-5,-2.991,1.666,.225))\n", + "print(\"%.3fde2 +%.3fde3+ %.3fdf2 %.3fdf3 = %.3f\" %(-5,19.1,1.666,-6.966,.65))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "set of linear equations at the end of first iteration are\n", + "2.846de2 -1.666de3+ 8.975df2 -5.000df3 = 2.750\n", + "-1.666de2 +6.366de3 -5.000df2 +20.900df3 = -0.300\n", + "8.525de2 -5.000de3 -2.991df2 +1.666df3 = 0.225\n", + "-5.000de2 +19.100de3+ 1.666df2 -6.966df3 = 0.650\n" + ] + } + ], + "prompt_number": 4 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 18.5, Page No 617" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "#initialisation of variables\n", + "Q2=-0.225\n", + "dP2=0.2-(-.075)\n", + "dP3=-0.6-(-0.3)\n", + "dQ3=-0.25-(-0.9)\n", + "\n", + "#Calculations\n", + "dV2=1.04**2 - 1**2 #dV2=|dV2|^2\n", + "\n", + "#Results\n", + "print(\"set of linear equations at the end of first iteration are\")\n", + "print(\"%.3fde2 %.3fde3+ %.3fdf2 %.3fdf3 = %.3f\" %(2.846,-1.666,8.975,-5,2.75))\n", + "print(\"%.3fde2 +%.3fde3 %.3fdf2 +%.3fdf3 = %.3f\" %(-1.666,6.366,-5,20.90,-.3))\n", + "print(\"%.3fde2 %.3fde3 %.3fdf2 +%.3fdf3 = %.3f\" %(8.525,-5,-2.991,1.666,.225))\n", + "print(\"%.3fde2 +%.3fde3+ %.3fdf2 +%.3fdf3 = %.5f\" %(2,0,0,0,dV2))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "set of linear equations at the end of first iteration are\n", + "2.846de2 -1.666de3+ 8.975df2 -5.000df3 = 2.750\n", + "-1.666de2 +6.366de3 -5.000df2 +20.900df3 = -0.300\n", + "8.525de2 -5.000de3 -2.991df2 +1.666df3 = 0.225\n", + "2.000de2 +0.000de3+ 0.000df2 +0.000df3 = 0.08160\n" + ] + } + ], + "prompt_number": 5 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file -- cgit