diff options
Diffstat (limited to 'Introduction_to_Electric_Drives_by_J._S._Katre/chapter8.ipynb')
| -rwxr-xr-x | Introduction_to_Electric_Drives_by_J._S._Katre/chapter8.ipynb | 806 |
1 files changed, 0 insertions, 806 deletions
diff --git a/Introduction_to_Electric_Drives_by_J._S._Katre/chapter8.ipynb b/Introduction_to_Electric_Drives_by_J._S._Katre/chapter8.ipynb deleted file mode 100755 index d3af51e4..00000000 --- a/Introduction_to_Electric_Drives_by_J._S._Katre/chapter8.ipynb +++ /dev/null @@ -1,806 +0,0 @@ -{ - "metadata": { - "name": "", - "signature": "sha256:9acc72237a22bafac745af0975ba25d35e2bd8c4f39083e9c7def597741576b5" - }, - "nbformat": 3, - "nbformat_minor": 0, - "worksheets": [ - { - "cells": [ - { - "cell_type": "heading", - "level": 1, - "metadata": {}, - "source": [ - "Chapter8, Control of DC drives" - ] - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.12.1: page 8-26" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from math import pi\n", - "#back emf ,Required armature voltage and Rated armatuer current\n", - "#given data :\n", - "TL=45 # in N-M\n", - "N=1200 #in rpm\n", - "Rf=147 #in ohm\n", - "Ra=25 # in ohm\n", - "Kv=0.7032 \n", - "w=(2*pi*N)/60 \n", - "Vf=220 #in volts\n", - "Kt=Kv \n", - "If=Vf/Rf \n", - "T=TL \n", - "Ia=T/(Kt*If) \n", - "Eg=Kv*w*If \n", - "print \"part (a)\"\n", - "print \"Back emf,Eg = %0.2f Volts\" %Eg\n", - "print \"part (b)\"\n", - "Ea=(Ia*(Ra/100))+Eg \n", - "print \"Required armature voltage, Ea = %0.2f V\"%Ea\n", - "print \"part (c)\"\n", - "rac=11191.4/Vf #\n", - "print \"rated armature current = %0.2f A\" %rac" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "Back emf,Eg = 132.25 Volts\n", - "part (b)\n", - "Required armature voltage, Ea = 142.94 V\n", - "part (c)\n", - "rated armature current = 50.87 A\n" - ] - } - ], - "prompt_number": 22 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.12.2: page 8-27" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from math import sqrt, cos, pi, acos\n", - "#the field current,Evaluation of alfa,Evaluation of power factor\n", - "#given data :\n", - "TL=50 # in N-M\n", - "N=1000 #in rpm\n", - "Rf=150 #in ohm\n", - "Ra=.25 # in ohm\n", - "Kv=0.7032 \n", - "alfa=0 \n", - "Vm=230 # in volts\n", - "Ef=((Vm*sqrt(2))/pi)*(1+cos(pi/180*alfa)) \n", - "If=Ef/Rf \n", - "print \"part (a)\"\n", - "print \"Field current, If = %0.2f A\" %If\n", - "print \"part (b)\"\n", - "w=(2*pi*N)/60 \n", - "Ia=TL/(Kv*If) \n", - "Eg=Kv*w*If \n", - "Ea=Eg+(Ra*Ia) \n", - "alfa_a=acos(((Ea*pi)/(Vm*sqrt(2)))-1)*180/pi \n", - "print \"angle = %0.2f degree\" %alfa_a\n", - "print \"part (c)\"\n", - "Ismax=Ia*((180-alfa_a)/180)**(1/2) #in amperes\n", - "PF=((Ea*Ia)/(Vm*Ismax)) #lagging\n", - "print \"power factor = %0.3f lagging\" %PF" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "Field current, If = 1.38 A\n", - "part (b)\n", - "angle = 83.90 degree\n", - "part (c)\n", - "power factor = 0.682 lagging\n" - ] - } - ], - "prompt_number": 24 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.12.3: page 8-29" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#torque \n", - "#given data :\n", - "Ia=50 # in A\n", - "Rf=150 #in ohm\n", - "Ra=.25 # in ohm\n", - "Kv=1.4 # in V/A-rad/sec\n", - "alfa_f=0 \n", - "alfa_a=45 # in degree\n", - "Vm=230*sqrt(2) # in volts\n", - "Vs=230 # in volts\n", - "Ef=((2*Vm)/pi)*(cos(pi/180*alfa_f)) \n", - "If=Ef/Rf \n", - "T=Kv*Ia*If \n", - "print \"part (a)\"\n", - "print \"Torque developed by the motor, T = %0.2f N/m\" %T\n", - "Ea=((2*Vm)/pi)*(cos(pi/180*alfa_a)) \n", - "Eg=Ea-(Ia*Ra) \n", - "w=Eg/(Kv*If) \n", - "N=(w/(2*pi))*60 \n", - "print \"part (b)\"\n", - "print \"Speed, N = %0.2f rpm \"%N\n", - "print \"part (c)\"\n", - "Ea=Eg+(Ra*Ia) \n", - "alfa_a=180/pi*acos(((Ea*pi)/(Vm*sqrt(2)))-1) \n", - "Ismax=Ia*((180-alfa_a)/180)**(1/2) #in amperes\n", - "PF=((Ea*Ia)/(Vm*Ismax)) #lagging\n", - "print \"power factor = %0.4f lagging\" %PF\n", - "#supply power factor is calculated wrong in the textbook" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "Torque developed by the motor, T = 96.63 N/m\n", - "part (b)\n", - "Speed, N = 661.71 rpm \n", - "part (c)\n", - "power factor = 0.6366 lagging\n" - ] - } - ], - "prompt_number": 26 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.12.4: page 8-32" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#Motor torque\n", - "#given data :\n", - "Vs_rms=230 # in volts\n", - "N=1200 # in rpm\n", - "Ia=40 # in A\n", - "Ra=0.25 #in ohm\n", - "Ka_fi1=0.182 # in V/rpm\n", - "Ka_fi=(0.182*60)/(2*pi) \n", - "alfa_a=30 \n", - "T=Ka_fi*Ia \n", - "print \"part (a)\"\n", - "print \"Motor torque, T = %0.1f N-m\" %T\n", - "print \"part (b)\"\n", - "Ea=((2*sqrt(2)*Vs_rms)/pi)*(cos(alfa_a*pi/180)) \n", - "N=(Ea-(Ra*Ia))/Ka_fi1 \n", - "print \"Speed of the motor, N = %0.1f rpm \"%N\n", - "print \"part (c)\"\n", - "Is_rms=Ia \n", - "PF=(Ea*Ia)/(Vs_rms*Is_rms) \n", - "print \"Power factor, PF = %0.4f lagging\" %PF" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "Motor torque, T = 69.5 N-m\n", - "part (b)\n", - "Speed of the motor, N = 930.4 rpm \n", - "part (c)\n", - "Power factor, PF = 0.7797 lagging\n" - ] - } - ], - "prompt_number": 30 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.12.6 : page " - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from numpy import arange, nditer\n", - "#Torque speed charaterstics\n", - "#given data :\n", - "v=230 #in volts\n", - "vm=sqrt(2)*v #in clts\n", - "Ka=1 \n", - "QR=1 #\n", - "ra=0.05 #\n", - "alpha=30 #in degree\n", - "y=(60/(2*pi)) #\n", - "z=((vm/pi)*(1+cos(pi/180*alpha))) #\n", - "x=(ra/(0.5)**2)\n", - "i = arange(0,9)\n", - "def func(i):\n", - " it = nditer([i, None])\n", - " for a, b in it:\n", - " b[...] = (z-a*x)*y\n", - " return it.operands[1]\n", - "\n", - "\n", - "wm = func(i)\n", - "print \"varoius values of speed in RPM is\"\n", - "for x in nditer([wm]):\n", - " print x,'\\t',\n", - "T=arange(0,9)\n", - "\n", - "###############PLOT#############\n", - "\n", - "\n", - "%matplotlib inline\n", - "import matplotlib.pyplot as plt\n", - "plt.plot(T, wm)\n", - "plt.xlabel(\"Torque ,N-m\")\n", - "plt.ylabel(\"Speed (rpm) for alpha=30 degree\")\n", - "plt.show()\n", - "\n", - "\n" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "varoius values of speed in RPM is\n", - "1844 \t1843 \t1841 \t1839 \t1837 \t1835 \t1833 \t1831 \t1829 \t" - ] - }, - { - "metadata": {}, - "output_type": "display_data", - "png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEPCAYAAABcA4N7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuYXXV97/H3ByIYbsmDWFG5BJUCRcAQiUW5DAgWbLUc\nK1K8Huw5tEUPNmgRUCAeoVgQby1UW5MIVvGOSgmXIzKIlxokgYQAAirl0ohggRBBSsjn/LF+A9th\nZs+embVn7z3783qe/WSt395rr2/EzHe+63eTbSIiIkazUacDiIiI7pZEERERTSVRREREU0kUERHR\nVBJFREQ0lUQRERFNtS1RSFos6T5Jqxra5ktaJmmFpOsk7TPsmh0krZP03hG+79uN3xUREVOjnRXF\nEuCwYW1nA6fangucVs4bfQy4dPgXSXoD8AiQSR8REVOsbYnC9rXAg8Oa1wCzyvFs4N6hNyQdAfwc\nuLnxAklbAAuAMwC1K96IiBjZjCm+30nA9yV9lCpJvRKeSgYnAocAfzvsmg8DHwUencI4IyKimOrO\n7EXA8bZ3oKoSFpX2hcDHbT9KQ9Ug6WXAi2x/i1QTEREdoXau9SRpDnCJ7T3K+VrbW5VjAQ/ZniXp\ne8D25bLZwAaqPowngVOB/6aqfn4P+IHtg0e4V/ovIiImwHbTX8SnuqK4Q9KB5fhg4DYA2wfY3sn2\nTsAngDNtn2f707ZfWNr3A24bKUkMsc2jj5qLLzZvfrOZNcsMDJjzzjNr1hi786/TTz+94zEkzsSY\nOBPn0KsV7RweexHwQ2AXSXdLOgY4Fjhb0g1UndPHjucraWHU08yZcMQR8IUvwJo18J73wA9+ALvu\nCgcdBOefD7/85UT+RhER/altndm2jx7lrVeMcd2HRmm/E9hzPDEMJY0jjoDHHoMrroCvfhVOOQXm\nzoUjj4Q3vAG23XY83xoR0V/6ZmZ2t1QaAwMD7b9JDRJnfXohRkicdeuVOFvR1s7sqSTJE/m7NFYa\nl16aSiMi+oskPEZndt8nikZJGhHRb5IoJiFJIyL6QRJFTZI0ImK6SqJogySNiJhOkijaLEkjInpd\nEsUUStKIiF6URNEhSRoR0SuSKLrA8KRx9NFw9tmw5ZadjiwiorVE0TczszulcUb4nXfC44/DnnvC\nd7/b6cgiIlqTiqIDli6Fv/xLeN3rqupiiy06HVFE9KtUFF3qta+FVauqx1J77JHqIiK6WyqKDkt1\nERGdlIqiB6S6iIhul4qii6S6iIiploqix6S6iIhulIqiS6W6iIipkIqih6W6iIhu0dZEIWmxpPsk\nrWpomy9pmaQVkq6TtM+wa3aQtE7Se8v5ZpIulXSLpJskndXOmLvJ7NmwZAmcdx684x1w3HGwbl2n\no4qIftPuimIJcNiwtrOBU23PBU4r540+BlzacG7gbNu7AXOBV0ka/p3TWqqLiOiktiYK29cCDw5r\nXgPMKsezgXuH3pB0BPBz4OaG73jM9jXl+AlgOfDCNobdlVJdRESndKKP4iTgXEl3AecApwBI2gI4\nEVg42oWSZgOvA65qf5jdKdVFREy1MROFpI0kvU3SaeV8B0nzJ3HPRcDxtncAFpRzqBLEx20/Cjyj\nB17SDOAi4JO275zE/XteqouImEpjDo+V9GlgA3Cw7V0lbQ1cafvlLd1AmgNcYnuPcr7W9lblWMBD\ntmdJ+h6wfblsdrnnqbbPL59dDKy1/Tej3Menn376U+cDAwMMDAy0EmJPe+ghWLAABgdh0SI4+OBO\nRxQR3WxwcJDBwcGnzj/0oQ9Nfj8KSStszx36s7TdaHuvVoIaIVEsBxbYvkbSq4GP2B4+8ul04BHb\nHyvnZwC7AkeONllius2jGK/Mu4iIiahrHsV/S9q44UufS/XbfisBXAT8ENhF0t2SjgGOBc6WdANw\nRjlv9h3bUfVj7AYsL8Nq39nK/ftJ+i4iol1aqSjeCrwJmAdcALwR+KDtr7Q/vNb1e0XRKNVFRLSq\ntq1QJe0GvLqcXmX7lhriq1USxe9K30VEtKLORLE/8BLbS8qjpy1s/6KmOGuRRDGyVBcR0UwtfRSS\nFlLNbzi5NG0C/Ouko4spkb6LiJisVjqz/wfwp8BvAGzfC2zZzqCiXpl3ERGT0UqieNz2U6OcJG3e\nxniijVJdRMREtJIovirpM8BsScdSLZ/x2faGFe2S6iIixqtpoigzp78MfL28fp9qtvSnpiC2aKNU\nFxHRqqajnkqiWGX7pVMX0sRk1NPEZWRURP+a9Kin8pP3+kkuAhhdLtVFRDTTyszsnwIvAf6DMvKJ\nKofs2ebYxiUVRT0aq4u//3vYMuPbIqa1WibclUX9nqHblvpOoqhPZnVH9I+6EsXWIzQ/Unab6xpJ\nFPVL30XE9FfX6rHLgQeA28vrAeA/JC2XNG/yYUa3St9FREBrieL/AYfbfo7t5wCHAf8GvAv4p3YG\nF52XeRcR0Uqi2Nf2FUMntq8sbT+iWvcp+kCqi4j+1UqiWCPp/ZJ2lDRH0onAfWUzo5Y2MIrpIdVF\nRH9qJVG8mWov628CFwM7AEcDG1NtaBR9JtVFRH9paT8KqBYDtP2bsT/ZGRn11BkZGRXR2+raj+KV\nkm4Gbi3ne0k6v6YYo8eluoiY/lp59PQJqpFODwDYvhE4cKyLJC2WdJ+kVQ1t8yUtk7RC0nWS9hl2\nzQ6S1kl6b0PbPEmrJN0u6ZOt/sVi6qTvImJ6ayVRYPuuYU3rW7hsCVWCaXQ21eqzc4HTynmjjwGX\nDmv7J+AvbO8M7Cxp+HdGl0h1ETE9tZIo7pL0KgBJm0h6H3DLWBfZvhZ4cFjzGmBWOZ4N3Dv0hqQj\ngJ8DNze0PR/Y0vay0nQhcEQLMUeHpLqImH5aSRR/TTW57oVUP9jnlvOJOAk4V9JdwDnAKQCStqDa\nl3vhsM+/ELin4fze0hZdLtVFxPQxY6wP2L6faohsHRYBx9u+WNKR5fxQqgTxcduPlj0wJmThwoVP\nHQ8MDDAwMDCpYGNyhqqLpUur6iIjoyI6b3BwkMHBwXFdM+rwWEn/0HBqQOXPqsE+fswvr1aevcT2\nHuV8re2tyrGAh2zPkvQ9qrkaUD2S2gCcCnwDuNr2buWao4EDbf/VCPfK8NgulhVpI7rTZIfHXl9e\nmwJ7A7cBd1A9epro0h13SBoaMXVw+U5sH2B7J9s7UY2yOtP2+bZ/CayV9IqSWN5GNfEvekz6LiJ6\nVyvLjP8Y2G9oWXFJzwK+b/sVY1x3EdUw2m2A+6hGOa0CzqNKPo8Bx9leMey606mWMf9YOZ8HfA6Y\nCSwdrZJJRdE7Ul1EdI+69qP4KfBK278u51sDP7K9S22R1iCJovdkVndE59W1H8VHgOWSPifpAqr9\nKc6qI8DobxkZFdEbWlrrqcxneAVVZ/Yy22vaHdh4paLobakuIjqjrooC22tsf9P2t7oxSUTvS3UR\n0b1aXj2226WimD5SXURMndoqioiplOoioru0MurpecB2VP0T99q+byoCG69UFNNTqouI9ppURSFp\nrqR/B64B/p5qpddrJP27pL3rDTViZKkuIjqv2RIeNwLH2v7xsPY/BD5je68piK9lqSimv1QXEfWb\nbB/FZsOTBIDtfwc2n2xwEeOV6iKiM5pVFJ8CXgJcANxNtSjg9sDbgZ/bfvdUBdmKVBT9JdVFRD0m\nvYSHpNcCfwq8oDTdC3zb9tLaoqxJEkX/yZpREZNXy1pPvSKJon+luoiYuMmOepol6YOS/krSxpJO\nl/Rvkj4saWb94UZMTPouItqrWR/FxVR7WG8G7AGsBL4EvB7Y2vY7pyrIVqSiCEh1ETFek3r0JGml\n7T3LhkFrgBfY3lDOb7S9Z/0hT1wSRQxJ30VE6yY7PHYDQPnpe5ntxvOIrpXd9CLq1XQrVElbAtg+\nZqhR0kuAte0OLGKy0ncRUY8JjXpSFz7n6cKQoouk7yJiZLWtHitp1/LnbpDHT9F7Ul1ETFyry4x/\ncdifY5K0WNJ9klY1tM2XtEzSCknXSdqnoX1Fea2UdFTDNcdIWiXpRkmXSXpOqzFENErfRcTEtJoo\nmpYlo1gCHDas7WzgVNtzgdPKOcAqYF5pfw1wXpm7sQnwUeDAsgjhSqCrlg6J3pPqImJ82rZxke1r\ngQeHNa8BZpXj2VRLgmD7saFRVcBM4GHbTwLry3dsUYblbjV0TcRkpLqIaN1U73B3EnCupLuAc4BT\nht4oj59WA6uBEwBK8ngPcBNVgtgNWDzFMcc0luoiYmwzpvh+i4DjbV8s6chyfiiA7WXA7qXj/HJJ\nV5drPgXsZfsXkv4BOBk4c6QvX7hw4VPHAwMDDAwMtOvvEdPIUHWxdGlVXWRkVExng4ODDA4Ojuua\nlobHSlphe66kG2y/rOUvl+YAl9jeo5yvtb1VORbwkO1ZI1x3FXAiVSI70/Yhpf0A4P22/3iEazIY\nKyYts7qj39Q2PBY4oPy5/+RC4g5JB5bjg4HboEookmaU4x2BnYHbqdaa2lXSNuWaQ4GbJxlDxKjS\ndxHxTC0lCtuPNP7ZCkkXAT8EdpF0t6RjgGOBsyXdAJxRzgH2A26QtAL4KtUWrGtt30/Vj3F12Zp1\nT+DvWo0hYqLSdxHxtDEfPUnal6qf4A+ATYCNgXVDj5C6RR49RbtkVndMZ3U9evpH4M1Uj4meDfwF\ncP7kw4voDakuot+1+ujpdmBj20/aHmkiXcS0lr6L6GetJIrfSNoUuFHS2ZJOYGIztSN6XqqL6Eet\n9FHMAe6j6p9YQDU7+nzbd7Q7uPFIH0VMtfRdxHQwqR3uek0SRXRC5l1Er6slUUjaDzgdmMPTM7lt\n+0V1BFmXJIropFQX0avqGvW0CPgY1VyHfcpr/uTDi5g+0ncR01krFcWPbb9iiuKZsFQU0S1SXUQv\nmVRFIWmepHlUs6LPkbSvpL2HXrVHGzFNpLqI6WbUikLSIDD0phqOAbB9UFsjG6dUFNGNUl1Et8uo\np4gukJFR0c3qGvW0DdWop/2oqoprgf9r+9d1BVqHJIrodqkuohvVNerpS8CvgDcAbwTuB748+fAi\n+kv6LqJXtVJR3GT7pcPaVg1tRtQtUlFEL0l1Ed2iroriSklHS9qovI4CrqwnxIj+lOoiekkrFcU6\nYDNgQ2naCPhNOXa37EuRiiJ6VaqL6KRaKgrbW9jeyPaM8trI9pbl1RVJIqKXpbqIbtdsHsU8hs2d\naGR7ebuCmohUFDEdpLqIqTap4bHDJtw9QybcRbRH5l3EVOrohDtJi4E/Bn41NEJK0nyqrVWfBawH\njrN9XWn/TLl0Y+BM218u12xSrjmQqp/kA7a/McL9kihiWkl1EVOhtkQhaQ9gN6o9swGwfeEY1+wP\nrAMubEgUg8BZtq+QdDhwou2DJM0EHre9QdK2wE3A82w/KelDJc7Tync8Z6TJfkkUMR2luoh2ayVR\nzGj2ZvmShVS/ze8OXAocDnwfaJoobF9bdsdrtAaYVY5nA/eWzz7W8JmZwMO2nyznxwC7NHxvV80I\nj2inob26ly6t9upOdRGd0Mo8ijcChwBrbB8D7EX1Q34iTgLOlXQXcA5wytAbkuZLWg2sBk4obUP3\nOUPS9ZK+Iun3JnjviJ6VkVHRSWNWFMBj5RHQekmzqJbz2H6C91sEHG/7YklHlvNDAWwvA3aXtCtw\nuaSrS3zbAT+w/V5JC4CPAm8f6csXLlz41PHAwAADAwMTDDOi+6S6iDoMDg4yODg4rmtamXB3PvAB\n4CjgvVST7VaU6mKsa+cAlzT0UawdmnshScBDtmeNcN1VwInAcuAR21uU9u2By4YvKVLeSx9F9I30\nXURd6ppwd5ztB21/GngN8I5WksQo7pB0YDk+GLitBDpH0oxyvCOwM3B7+cl/iaShobivpno0FdHX\nhqqL886rqovjjoN16zodVUxX7RweexFVJ/g2wH3AacAq4DxgU+AxquGxKyS9lar/4onyOs325eV7\ndgA+T9Uv8ivgGNv3jHC/VBTRl1JdxGRk46KIPpJ5FzERda0eGxE9ICOjol2aVhSl32C17V1G/VCX\nSEUR8bRUF9GqSVcUttcDt5YO5ojoEakuok6tDI+9FpgLLON396F4fZtjG5dUFBEjS3URzdTSmS1p\noBwOfVBUieKaSUdYoySKiNFlZFSMps5FAbcF9qFKFsts/6qeEOuTRBExtlQXMVwto54kvQn4MXAk\n8CZgWVl+IyJ6TPouYiJaefS0EjhkqIqQ9FzgKtt7TkF8LUtFETE+qS4C6ptHIeD+hvNfl7aI6GGp\nLqJVrVQU51AtLf5FqgRxFLDS9ontD691qSgiJi7VRf+aVEUh6dkAtv+WapvSPYE9gM90W5KIiMlJ\ndRHNjFpRSFpue29Jn7f9timOa9xSUUTUI9VFf5lsH8Wmkt4CvErSGxpefybpDfWGGhHdItVFDNes\notgfeAvVsNhvD39/EntStEUqioj6pbqY/uqamf2/bH+21sjaIIkioj0yq3t6y34UEVGbVBfTU/aj\niIjapO+if6WiiIhxS3UxfaSiiIi2SHXRX9qWKCQtlnSfpFUNbfMlLZO0QtJ1kvZpaF9RXislHTXC\n93278bsiorNmz4YlS+C88+Ad74DjjoN16zodVbRDOyuKJcBhw9rOBk61PRc4rZwDrALmlfbXAOdJ\n2njoojJv4xGe3hMjIrpEqovpr6VEIWlzSbtK2kXS5q1cY/ta4MFhzWuAWeV4NnBv+exjtjeU9pnA\nw7afLPfeAlgAnEEWI4zoSqkuprdmaz1tKekEScuofuNfAlwA3CTpJ5IWlB/i43EScK6ku4BzgFMa\n7jdf0mpgNXBCwzUfBj4KPDrOe0XEFEt1MT3NaPLeN4EvAa+zfV/jG2XHu9cD3wJePY77LQKOt31x\n2fxoEXAogO1lwO6SdgUulzQI7AS8yPYCSXPG+vKFCxc+dTwwMMDAwMA4QouIOgxVF0uXVtVFRkZ1\nl8HBQQYHB8d1TVuHx5Yf7pfY3qOcr7W9VTkW8JDtWSNcdxXwfuDlwKnAf1Mltd8DfmD7GXNDMzw2\novtkVnf3a2V4bLOKovGL9gLmNHzetr8xgZjukHSg7WuAg4HbyvfPAe6xvV7SjsDOwG22fwJ8unxm\nR+DfRkoSEdGdUl1MD2MmCklLqPahWA1saHiraaKQdBFwILCNpLupRjkdSzWiaVPgsXIOsB9wkqQn\ngCeAY22vHf6VZNRTRE8a6rtYsKDqu0h10VtaWRTwZmD3bn+uk0dPEb0hs7q7S10zs68D/qCekCKi\n3w1VF7/9bUZG9YpWKooBqv0ofgk8Xppte8/2hjY+qSgiek+qi86raz+Kn1FNeLuJhj4K23fWEGNt\nkigietNDD8EJJ8DVV6fvohPqShQ/sr1vrZG1QRJFRG9LddEZdfVRrJD0RUlHl/2ys2d2RNQufRfd\nq5WK4nPl8Hc+mD2zI6JdUl1MnWyFGhE9K30XU6OWR0+SXizpEkkPSLpf0rckvai+MCMinmn2bFi8\nOCvSdoNW+ii+CHwFeD7wAuCrwEXtDCoiYkj6LjqvlT6KlcPnTEi60fZebY1snPLoKWL6u+wyOPbY\n9F3Uqa5RT5dJOlnSnPJ6f2nbWtLW9YQaETG2ww9PddEJrVQUdzL6Yny23RX9FakoIvpLRkbVY9IV\nhaSNgLfY3mmUV1ckiYjoP+m7mDqtVBQ32H7ZFMUzYakoIvpXqouJq6uP4juS3lh2pIuI6DqpLtqr\nlYpiHbAZ8CTw29LsoS1Nu0UqioiAVBfjVUtFYXsL2xvZfpbtLcurq5JERMSQVBf1G7WikPRi2z9r\nenELn5kqqSgiYrhUF2Ob1FpPkr4MbE61adFPgDVU+1Y/H3g58HrgEdt/XmfQE5VEEREjyZpRzU16\nUUBJLwH+HHgVsGNp/g/g+8BFtn8+RgCLgT8GfmV7j9I2H/hH4FnAeuA429eV9s+USzcGzrT9ZUmb\nUS0b8iKqfpJLbJ88wr2SKCJiVKkuRtbx1WMl7Q+sAy5sSBSDwFm2r5B0OHCi7YMkzQQet71B0rZU\nO+o9D9gEmG/7GknPAq4C/s725cPulUQREU2lunimuobHTpjta4EHhzWvAWaV49nAveWzj9ke2mp1\nJvCw7SdL+zXlM08Ay4EXtjPuiJiesiLtxLQ1UYziJOBcSXcB5wCnDL0hab6k1cBq4IThF0qaDbyO\nqqqIiJiQjIwanxkduOci4HjbF0s6spwfCmB7GbC7pF2ByyUN2n4YQNIMquXNP2n7zpG+eOHChU8d\nDwwMMDAw0Ma/RkT0sqHqYunSqrrol76LwcFBBgcHx3VNs1FP86gWAxQjLApoe3lLN5DmUHVAD/VR\nrB2ah1Fmez9ke9YI111F1X9xfTlfDKy1/Tej3Cd9FBExIf3cd9FKH0WziuJcqgQxE5gHrCzte1IN\nl913gnHdIenA0u9wMHBbCXYOcI/t9ZJ2BHYGbi/vnQFsBfzFBO8ZETGqfq0uWjVqH4XtAdsHAf8J\n7G17nu15wNzSNiZJFwE/BHaRdLekY4BjgbMl3QCcUc4B9gNukLSCajjssbbXStqOqh9jN2C5pBWS\n3jmhv21ERBPpuxhZK2s93Wz7D8Zq67Q8eoqIOvXLvIu6hseulPRZSQOSDpL0L8CN9YQYEdGdUl08\nrZWKYibw18D+pel7wD/Z/u3oV029VBQR0S7TubqobWZ2WUZjB9u31hVc3ZIoIqKdpuvIqFoePUl6\nPbACuLycz5X07XpCjIjoDf08q7uVPoqFwCsoS3HYXkG1QF9ERN/px76LVhLFE7YfGta2YcRPRkT0\ngX6rLlpJFKslvQWYIWlnSf9ANTciIqKv9Ut10cqop82BDwCvKU1XAB/OqKeIiKf16sioWvejkLS5\n7d/UElkbJFFERKf14sioukY9vVLSzcCt5XwvSefXFGNExLQx1Hdx/vnTq++ilT6KTwCHAQ8A2L4R\nOLCdQUVE9LLDD59efRctbVxk+65hTevbEEtExLQxnUZGtZIo7pL0KgBJm0h6H3BLe8OKiJgepsPI\nqFZGPT0X+CRwCNUmRldS7VD36/aH17p0ZkdEt+vGkVG1jnrqdkkUEdELum1kVF2jnl4s6RJJD0i6\nX9K3JGUJj4iICejFvotW+ii+CHwFeD7wAqrd5y5qZ1AREdNdL/VdtNJHsdL2nsPabrS9V1sjG6c8\neoqIXtXJvou6dri7TNLJkuaU1/tL29aStm5y88WS7pO0qqFtvqRlZd/r6yTt09C+orxWSjqq4Zp5\nklZJul3SJ1uINyKip3R7ddFKRXEnMNqHbHvE/gpJ+wPrgAtt71HaBoGzbF8h6XDgRNsHlV30Hre9\nQdK2wE3A82w/KWkZ8G7byyQtBT5l+/IR7peKIiJ63lRXF7VUFLbn2N5plNeondq2r6XsYdFgDTCr\nHM8G7i2ffcz20NLlM4GHS5J4PrCl7WXlvQuBI8aKOSKiV3VjdTFqoiiPg57fcP4OSd+W9Klmj5zG\ncBJwrqS7gHOAU4bdbzWwGjihNL8QuKfh+ntLW0TEtNVtI6OaVRSfAR4HkHQA8BHgAmAt8M8TvN8i\nqsl6OwALyjkAtpfZ3h3YG/ikpFmjfEdERF/olupiRpP3NrL9X+X4KOAztr8OfF3SjRO833zbh5Tj\nrwGfHf4B27dK+hnwEqpqYruGt7ejPK4aycKFC586HhgYYGBgYIJhRkR0h6HqYunSqrqYbN/F4OAg\ng4OD47pm1M5sSTcBc20/IemnwLG2rynvrS6//Tf/cmkOcElDZ/ZyYIHtayS9GviI7X3K5+6xvV7S\njsC1wEttr5X0Y+B4YBlwKenMjog+1Y5Z3a10ZjerKC4CrpH0APAo1Q9vJO0MDN9De6SbX0S1HPk2\nku4GTgOOBc6TtCnwWDkH2A84SdITwBNUSWltee844HNUndxLR0oSERH9oO7qolVNh8dK2hfYFrhy\naHc7Sb8PbGF7eXtDG59UFBHRT+qqLrIoYETENDfZeRd1zcyOiIguNRUjo1JRRERMExOpLlJRRET0\nkXZVF6koIiKmoVari1QUERF9qs7qIhVFRMQ016y6SEURERGTri5SUURE9JHh1cWWW6aiiIiIBsOr\ni1akooiI6FNXXgl/9EdZwiMiIppIZ3ZERExaEkVERDSVRBEREU0lUURERFNJFBER0VQSRURENJVE\nERERTbUtUUhaLOk+Sasa2uZLWiZphaTrJO1T2g+V9BNJK8ufBzVcc4ykVZJulHSZpOe0K+aIiHim\ndlYUS4DDhrWdDZxqey5wWjkHuB/4E9t7Au8APg8gaRPgo8CBtvcCVgLvbmPMbTc4ONjpEFqSOOvT\nCzFC4qxbr8TZirYlCtvXAg8Oa14DzCrHs4F7y2dvsP3L0n4zMFPSs4D15Tu2kCRgq6FrelWv/J8n\ncdanF2KExFm3XomzFTOm+H4nAd+X9FGqJLXvCJ/5M+B6208ASHoPcBOwDrgNeNcUxRoREUx9Z/Yi\n4HjbOwALgMWNb0raHfgI8JflfCvgU8Betl8ArAJOntKIIyL6ne22vYA5wKqG87UNxwIebjjfDvgp\nsG9D2yuA7zScHwBcOsq9nFdeeeWV1/hfY/0sn+pHT3dIOtD2NcDBVI+SkDQbuBR4v+0fNXz+58Cu\nkrax/QBwKFUfxjOMtfphRERMTNuWGZd0EXAgsA1wH9Uop1XAecCmwGPAcbZXSPogVf/F7Q1fcajt\nByS9HfhbYANwJ/A/bQ/vJI+IiDaZNvtRREREe/T8zGxJh0m6VdLtkt7f6XhGM9IExG4jaXtJV0ta\nLekmScd3OqaRSHq2pB9LukHSzZLO6nRMzUjauEwyvaTTsYxG0p1lwusKScs6Hc9oJM2W9DVJt5T/\n9n/Y6ZiGk7RL+d9x6PVwN/5bknRy+be+StIXJW066md7uaKQtDFVB/ghVPMrrgOOtn1LRwMbgaT9\nqYb4Xmi7xZ1qp5akbYFtbd8gaQvgeuCILv3fczPbj0qaAXwfeJ/t73c6rpFIOgGYB2xp+/Wdjmck\nkn4BzLP9X52OpRlJFwDX2F5c/ttvbvvhTsc1GkkbUf1smm/77k7HM0TSHOC7wG62H5f0ZWCp7QtG\n+nyvVxRjHUmgAAAFHklEQVTzgTts31nmXXwJ+NMOxzSiUSYgdhXbv7R9QzleB9wCvKCzUY3M9qPl\ncBNgY6Arf8BJ2g54LfBZqpF+3ayr45M0C9jf9mIA2+u7OUkUhwA/66YkUawFngA2Kwl3M5pMZu71\nRPFCoPE/wD2lLSap/MYxF/hxZyMZmaSNJN1ANVDiatsjjobrAh/n6cEY3czAd8paa/+708GMYifg\nfklLJC2X9C+SNut0UGP4c+CLnQ5iuFI5ngvcBfwn8JDt74z2+V5PFL373KyLlcdOXwPeUyqLrmN7\ng+2XUc2/OUDSQIdDegZJfwL8yvYKuvy3deBVZQ22w4F3lUel3WYGsDdwvu29gd9QjZbsSmWtutcB\nX+10LMNJejHwN1Rz3V5AtUzSW0b7fK8ninuB7RvOt6eqKmKCyhpbXwf+1fY3Ox3PWMqjh0uBl3c6\nlhG8Enh9ef5/EXCwpAs7HNOIbK8pf94PXEz1WLfb3APcY/u6cv41qsTRrQ6nWo7o/k4HMoKXAz+0\n/Wvb64FvUP3/dUS9nih+AuwsaU7J3kcB3+5wTD2rLLy4CLjZ9ic6Hc9oJG1TJmkiaSbVRMwVnY3q\nmWyfYnt72ztRPYL4ru23dzqu4SRtJmnLcrw58BqqOU9dpSwcerek3y9NhwCrOxjSWI6m+gWhG90K\n/KGkmeXf/SGMMpkZpn5RwFrZXi/p3cAVVB2ai7pxhA78zgTE50i6GzjN9pIOhzXcq4C3AislDf3g\nPdn25R2MaSTPBy4oI0o2Aj5v+6oOx9SKbn1U+jzg4urnBTOAL9i+srMhjer/AF8ovxj+DDimw/GM\nqCTcQ4Cu7O+xfWOpbn9C1X+2HPjn0T7f08NjIyKi/Xr90VNERLRZEkVERDSVRBEREU0lUURERFNJ\nFBER0VQSRURENJVEEX1F0nMaln9eI+mecry8LI7W6fg+V2LapJxvU2Z2R3RMx/9hREwl27+mWuwQ\nSacDj9j+2FjXSdrI9lQt7LceeCfw6Sm6X0RTqSii30nSq0tVsVLSoobf5u+U9BFJ1wNHlk2ybpF0\nvaRPDW1EJGmhpPc2fOFNknYox28tmyytkPTpMpu8GQOfBBaM9dlSfZwv6UeSfiZpQNIFZUOfbpv1\nHz0siSL63bOBJcCRtvekqrL/urxn4AHb84BvUS1x8Cfl/Hk8vSTH8OUNDCBpN+BNwCvLyqwbgFFX\n6GxwF9VmTG8f4buH32e27X2BBVTrnJ0N7A7sIWmvFu4VMaYkiuh3GwM/t31HOb8AOKDh/S+XP3cF\nfmH7Z+X8X2m+dLiAV1PtbPeTsnbWwVR7KozFwFlU+1iM9W90aHvVm4Bf2l7tal2e1VRLSEdMWvoo\nIn73B7743d/if9PCNev53R/oz244vsD2KeMNyPYdZWOmo566oXQm1W55LvsxAPx3+XMD8HjDV2wg\n/76jJqkoot89CcwpG7kAvA24ZoTP3Vo+96JyfjRPJ5Q7KfsiSNqbqmowcBXwRknPLe9tPdR30UjS\nWZKOaGwqf54JvG+o0fYHbM9tSBIRUyKJIvrdY1RLVX9V0kqq6mBotNFTlYXt3wLHApeWzu37ePoH\n+teBrSXdBLwL+Gm55hbgg8CVkm4ErgS2HSGGlwJrGs5drr8ZuJ6x+ylGOh7pPGJCssx4xARIOhB4\nn+3X1fBdl9s+rIawItoiFUXExNXyW1aSRHS7VBQREdFUKoqIiGgqiSIiIppKooiIiKaSKCIioqkk\nioiIaCqJIiIimvr/ge16EIi/jlkAAAAASUVORK5CYII=\n", - "text": [ - "<matplotlib.figure.Figure at 0x7f3de8ed3d90>" - ] - } - ], - "prompt_number": 39 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.18.1: page 8-52" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from math import sqrt, degrees, pi, acos, cos\n", - "from __future__ import division\n", - "#No load speed ,firing angle ,Power Factor and speed regulation\n", - "#given data :\n", - "Ra=0.075 #in ohm\n", - "alfa1=0 # in degree\n", - "alfa2=30 # in degree\n", - "VL_rms=480 # in volts\n", - "Ka_fi=0.3 # in V/rms\n", - "Vs_rms=round(VL_rms/sqrt(3)) \n", - "Vm=sqrt(2)*Vs_rms \n", - "Ea=round((3*sqrt(3)*Vm*cos(pi/180*alfa1))/pi) \n", - "Ea1=((3*sqrt(3)*Vm*cos(alfa2*pi/180))/pi) \n", - "Ia=(10/100)*160 # in A\n", - "N_0=(Ea-Ia*Ra)/Ka_fi \n", - "N_30=(Ea1-Ia*Ra)/Ka_fi \n", - "print \"part (a)\"\n", - "print \"No load speed at alfa=0 degree = %0.2f rpm \"%N_0\n", - "print \"No load speed at alfa=30 degree = %0.2f rpm \"%N_30\n", - "print \"part (b)\"\n", - "Ia=160 # in A\n", - "N=1800 # in rpm\n", - "Eg=540 # in volts\n", - "Ea=(Eg+(Ia*Ra)) \n", - "alfa=degrees(acos((Ea*pi)/(3*sqrt(3)*Vm))) \n", - "print \"the firng angel, alfa = %0.1f degree \"%alfa\n", - "print \"part (c)\"\n", - "Is_rms=sqrt(2/3)*Ia \n", - "Sva=3*Vs_rms*Is_rms \n", - "PF=(Ea*Ia)/(Sva) \n", - "print \"Power Factor, PF = %0.4f lagging\"%PF\n", - "print \"part (d)\"\n", - "Ra=0.075 #in ohm\n", - "Ia=160 # in A\n", - "Ia1=16 # in A\n", - "Eg=540 # in volts\n", - "Ka_fi=0.3 # in V/rms\n", - "N=1800 # in rpm\n", - "Ea=(Eg+(Ia*Ra)) \n", - "Eg1=Ea-(Ia1*Ra) \n", - "N_0=Eg1/Ka_fi \n", - "SR=((N_0-N)/N)*100 \n", - "print \"Speed Regulation, SR = %0.2f %%\" %SR" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "No load speed at alfa=0 degree = 2156.00 rpm \n", - "No load speed at alfa=30 degree = 1866.41 rpm \n", - "part (b)\n", - "the firng angel, alfa = 31.6 degree \n", - "part (c)\n", - "Power Factor, PF = 0.8135 lagging\n", - "part (d)\n", - "Speed Regulation, SR = 2.00 %\n" - ] - } - ], - "prompt_number": 41 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.18.2: page 8-54" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#Delay Angel of Armature,No load speed and speed regulation\n", - "#given data :\n", - "VL_rms=208 # in volts\n", - "Kv=1.2 # in V/A-rad/sec\n", - "Vs_rms=round(VL_rms/sqrt(3)) \n", - "Vm=sqrt(2)*Vs_rms \n", - "Rf=240 # in ohm\n", - "Ra=0.25 # in ohm\n", - "alfa_f=0 # in degree\n", - "V=280 # in volts\n", - "Twenty_HP=20*746 #in watt\n", - "Ia=Twenty_HP/V\n", - "Ef=round((3*sqrt(3)*Vm*cos(pi/180*alfa_f))/pi) \n", - "N=1800 \n", - "w=(N*2*pi)/60 \n", - "If=Ef/Rf \n", - "Eg=Kv*w*If \n", - "Ea=round(Eg+(Ia*Ra)) \n", - "alfa_a=degrees(acos((Ea*pi)/(3*sqrt(3)*Vm))) \n", - "print \"part (a)\"\n", - "print \"Delay Angel Of Armature, alfa_a = %0.2f degree\"%alfa_a\n", - "print \"part (b)\"\n", - "Ia1=(Ia*10)/100\n", - "Eg_noL=Ea-(Ia1*Ra) \n", - "w_0=(Eg_noL/(1.2*1.17)) # rad/sec\n", - "N_0=(w_0*60)/(2*pi) \n", - "print \"NO load speed at alfa|_a, = %0.2f \"%N_0\n", - "# no load speed is calculated wrong in textbook\n", - "print \"part (c)\"\n", - "SR=((N_0-N)/N)*100 \n", - "print \"Speed Regulation, SR = %0.2f %% \"%SR\n", - "# speed regulation is calculated wrong in the textbook" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "Delay Angel Of Armature, alfa_a = 7.94 degree\n", - "part (b)\n", - "NO load speed at alfa|_a, = 1881.75 \n", - "part (c)\n", - "Speed Regulation, SR = 4.54 % \n" - ] - } - ], - "prompt_number": 42 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.18.3: page 8-56" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#alphas,speed and delay angle\n", - "#given data :\n", - "v1=208 #\n", - "vsrms=v1/sqrt(3) #\n", - "n=1000 #rpm\n", - "w=n*(pi/30) #in rad/s\n", - "ang=0 #\n", - "ef=((3*sqrt(3)*sqrt(2)*vsrms*cos(pi/180*ang))/pi) #in volts\n", - "rf=140 #in ohms\n", - "If=ef/rf #in amperes\n", - "t=120 #N-m\n", - "kv=1.2 #\n", - "ia=(t)/(kv*If) #in amperes\n", - "eg=kv*If*w #in volts\n", - "ra=0.25 #in ohms\n", - "ea=eg+(ia*ra) #\n", - "alpha=degrees(acos((ea*pi)/(3*sqrt(3)*sqrt(2)*vsrms)))\n", - "print \"part (a)\"\n", - "print \"alpha = %0.2f degree\"%round(alpha)\n", - "print \"part (b)\"\n", - "rf=140 #in ohms\n", - "If=ea/rf #in amperes\n", - "t=120 #N-m\n", - "kv=1.2 #\n", - "ia=(t)/(kv*If) #in amperes\n", - "ra=0.25 #in ohms\n", - "eg=ea-(ia*ra) #\n", - "w=(eg/(kv*If)) #in rad/s\n", - "N=w*(30/pi) #rpm\n", - "print \"speed = %0.2f rpm\" %N\n", - "#speed is calculated wrong in the textbook\n", - "print \"part (c)\"\n", - "n1=1000 #rpm\n", - "w=n1*(pi/30) #in rad/s\n", - "v1=208 #\n", - "vsrms=v1/sqrt(3) #\n", - "w1=(1800*(pi/30)) #\n", - "n=1800 #rpm\n", - "ang=0 #\n", - "T=120 #n-m\n", - "alphas=0 #\n", - "ang=0 #\n", - "ea=((3*sqrt(3)*sqrt(2)*vsrms*cos(pi/180*ang))/pi) #in volts\n", - "rf=140 #in ohms\n", - "If=ea/rf #in amperes\n", - "t=120 #N-m\n", - "kv=1.2 #\n", - "ia=(t)/(kv*If) #in amperes\n", - "ra=0.25 #in ohms\n", - "eg=ea-(ia*ra) #\n", - "if1=eg/(kv*w1) #in amperese\n", - "ef1=if1*rf #in volts\n", - "alphaf=degrees(acos((ef1*pi)/(3*sqrt(3)*120*sqrt(2)))) \n", - "print \"delay angle = %0.2f degree\"%alpha\n", - "# Ans in the textbook are not accurate." - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "alpha = 20.00 degree\n", - "part (b)\n", - "speed = 1058.39 rpm\n", - "part (c)\n", - "delay angle = 19.62 degree\n" - ] - } - ], - "prompt_number": 44 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.19.1: page 8-58" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#Firing angle to keep the motor current and Power fed back \n", - "#given data :\n", - "Vs_rms=260 # in volts\n", - "Ia=40 # in A\n", - "Eg=192 #in volts\n", - "kv=0.182 # in V/rpm\n", - "Ra=0.3 # in ohm\n", - "Ea=Eg+(Ia*Ra) \n", - "alfa_a=degrees(acos((Ea*pi)/(2*Vs_rms*sqrt(2))) )\n", - "print \"part (a)\"\n", - "print \"Firing angle to keep motor current, alfa_a = %0.2f degree\" %alfa_a\n", - "Ea1=-Eg+(Ia*Ra) \n", - "alfa_b=degrees(acos((Ea1*pi)/(2*Vs_rms*sqrt(2))) )\n", - "print \"Firing angle, alfa_a = %0.2f degree\"%alfa_b\n", - "print \"part (b)\"\n", - "Ia=40 # in A\n", - "Eg=192 #in volts\n", - "Ra=0.3 # in ohm\n", - "Ea=-Eg+(Ia*Ra) \n", - "P=abs(Ea)*Ia \n", - "print \"Power fed back, P = %0.2f Watt\" %P" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "Firing angle to keep motor current, alfa_a = 29.37 degree\n", - "Firing angle, alfa_a = 140.26 degree\n", - "part (b)\n", - "Power fed back, P = 7200.00 Watt\n" - ] - } - ], - "prompt_number": 45 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.19.2 : page 8-58" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#Average armature voltage ,back emf ,speed of the motor , motor torque and supply power factor\n", - "#given data :\n", - "Vm=230 # in volts\n", - "Ia=40 # in A\n", - "Ra=0.5 # in ohm\n", - "Ka_fi=0.2 # in V/rpm\n", - "alfa=30 \n", - "Ea=(Vm*sqrt(2)*(1+cos(pi/180*alfa)))/pi \n", - "print \"part (a)\"\n", - "print \"Average armature current, Ea = %0.2f V\"%Ea\n", - "print \"part (b)\"\n", - "Eb=Ea-(Ia*Ra) \n", - "print \"Back emf, Eb = %0.2f V\"%Eb\n", - "print \"part (c)\"\n", - "N=Eb/Ka_fi \n", - "print \"Speed of the motor, N = %0.2f rpm\" %round(N)\n", - "print \"part (d)\"\n", - "Ka_fi1=(Ka_fi*60)/(2*pi) \n", - "T=Ka_fi1*Ia \n", - "print \"Torque, T = %0.1f N/m\" %T\n", - "print \"part (e)\"\n", - "alfa=pi/6 \n", - "PF=(2*sqrt(2)*cos(alfa/2)**2)/(sqrt(pi*(pi-alfa))) \n", - "print \"power factor = %0.2f lagging\" %PF" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "Average armature current, Ea = 193.20 V\n", - "part (b)\n", - "Back emf, Eb = 173.20 V\n", - "part (c)\n", - "Speed of the motor, N = 866.00 rpm\n", - "part (d)\n", - "Torque, T = 76.4 N/m\n", - "part (e)\n", - "power factor = 0.92 lagging\n" - ] - } - ], - "prompt_number": 47 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.19.3: page 8-59" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#torque developed,speed and input power factor\n", - "#given data :\n", - "v=208 #in volts\n", - "f=50 #in Hz\n", - "ra=0.5 #in ohms\n", - "rf=345 #in ohms\n", - "kv=0.71 #in V/A-rad/sec\n", - "alpha=45 #in degree\n", - "ia=55 #in amperes\n", - "If=((2*sqrt(2)*v*cos(0))/(pi*rf)) #in amperes\n", - "t=kv*If*ia #in N/m\n", - "print \"part (a)\"\n", - "print \"torque = %0.2f N/m\"%t\n", - "print \"part (b)\"\n", - "eb=((2*sqrt(2)*v*cos(pi/180*alpha))/pi)-(ia*ra) #in volts\n", - "w=eb/(kv*If) #in rad/sec\n", - "N=w/(2*pi) #rps\n", - "N*=60 # rpm\n", - "print \"speed = %0.2f rpm\" %N\n", - "#speed is calculated wrong in the textbook\n", - "print \"part (c)\"\n", - "ea=132.4 #in volts\n", - "ef=187.3 #in volts\n", - "pi=(ea*ia)+(ef*If) #in watts\n", - "Isrms=sqrt((ia)**2+(If)**2) #in amperes\n", - "va1=Isrms*v #in VA\n", - "Pf=pi/va1 #\n", - "print \"Power factor = %0.4f lagging\" %Pf" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "torque = 0.01 N/m\n", - "part (b)\n", - "speed = -681.49 rpm\n", - "part (c)\n", - "Power factor = 0.6365 lagging\n" - ] - } - ], - "prompt_number": 52 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.19.4: page 8-60" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "from math import sqrt, pi\n", - "#develepoed back emf,required armature voltage and firing angle and rated armature current\n", - "#given data :\n", - "hp=20 #\n", - "v=230 #volts\n", - "n=1000 #rpm\n", - "lt=50 #load torque in N-m\n", - "s=1000 #speed in rpm\n", - "ra=0.2 #in ohms\n", - "rf=150 #in ohms\n", - "la=10 #in mH\n", - "kv=0.7 #\n", - "vf=(2*sqrt(2)*v)/(pi) #\n", - "If=vf/rf #in amperes\n", - "ia=(lt/(kv*If)) #in amperes\n", - "eg=((kv*2*pi*n*If))/(60) #in volts\n", - "print \"part (a)\"\n", - "print \"back emf = %0.1f V\"%eg\n", - "print \"part (b)\"\n", - "ea=eg+(ia*ra) #in volts\n", - "\n", - "alpha=degrees(acos((ea*pi)/(2*sqrt(2)*v))) #\n", - "print \"armature voltage = %0.2f V\" %ea\n", - "print \"firing angle = %0.1f degree\" %alpha\n", - "print \"part (c)\"\n", - "ea1=220 #in volts\n", - "ha20=746*20 #\n", - "iar=(ha20/ea1) #in amperes\n", - "print \"rated armature current = %0.1f A\" %iar" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "part (a)\n", - "back emf = 101.2 V\n", - "part (b)\n", - "armature voltage = 111.54 V\n", - "firing angle = 57.4 degree\n", - "part (c)\n", - "rated armature current = 67.8 A\n" - ] - } - ], - "prompt_number": 68 - }, - { - "cell_type": "heading", - "level": 2, - "metadata": {}, - "source": [ - "Example 8.21.1 : page 8-65" - ] - }, - { - "cell_type": "code", - "collapsed": false, - "input": [ - "#Average armature current\n", - "#given data :\n", - "V=200 # in volts\n", - "D=50/100 # duty cycle\n", - "VL_dc=V*D \n", - "Eb=75 # in volts\n", - "Ra=1 # in ohm\n", - "Ia=(VL_dc-Eb)/Ra \n", - "print \"Average armature current, Ia = %0.2f A\" %Ia" - ], - "language": "python", - "metadata": {}, - "outputs": [ - { - "output_type": "stream", - "stream": "stdout", - "text": [ - "Average armature current, Ia = 25.00 A\n" - ] - } - ], - "prompt_number": 69 - } - ], - "metadata": {} - } - ] -}
\ No newline at end of file |