From fba055ce5aa0955e22bac2413c33493b10ae6532 Mon Sep 17 00:00:00 2001 From: hardythe1 Date: Tue, 5 May 2015 14:21:39 +0530 Subject: add books --- Power_Electronics/Power_electronics_ch_6_1.ipynb | 852 +++++++++++++++++++++++ 1 file changed, 852 insertions(+) create mode 100755 Power_Electronics/Power_electronics_ch_6_1.ipynb (limited to 'Power_Electronics/Power_electronics_ch_6_1.ipynb') diff --git a/Power_Electronics/Power_electronics_ch_6_1.ipynb b/Power_Electronics/Power_electronics_ch_6_1.ipynb new file mode 100755 index 00000000..9456765e --- /dev/null +++ b/Power_Electronics/Power_electronics_ch_6_1.ipynb @@ -0,0 +1,852 @@ +{ + "metadata": { + "name": "" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 6: AC Regulators" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.1, Page No.286" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# input/output power, input power factor, average and rms vales of thyristor\n", + "\n", + "\n", + "import math\n", + "from scipy.integrate import quad\n", + "#variable declaration\n", + "nc = 36.0 # conduction cycles\n", + "no = 64.0 # off cycles\n", + "V = 150.0 # input voltage\n", + "Rl = 8.0 # load resistance\n", + "\n", + "#calculations\n", + "#(a)\n", + "alfa = nc/(nc+no)\n", + "Vl = V*math.sqrt(alfa)\n", + "#(b)\n", + "P = (Vl**2)/Rl\n", + "#(d)\n", + "Il = Vl/Rl\n", + "va_i = V*Il\n", + "pf_i = P/va_i\n", + "#(e)\n", + "sqrt_2 = math.floor(math.sqrt(2)*1000)/1000\n", + "Ip = sqrt_2*V/Rl\n", + "\n", + "def f(x):\n", + " return (alfa/(2*math.pi))*Ip*(math.sin(x))\n", + "wt_lower=0\n", + "wt_upper =math.pi\n", + "val = quad(f,wt_lower,wt_upper)\n", + "\n", + "val2 =Ip*math.sqrt(alfa)/2.0\n", + "\n", + "#result\n", + "print(\"(a) rms output voltage, Vl = %.0f V\\n(b) Power output = %.1f W\"%(Vl,P))\n", + "print(\"(c) Since losses are neglected, input power = Output power = %.1f W\"%(P))\n", + "print(\"(d) input power factor = %.1f lagging\\n(e)\\tPeak thyristor current = %.4f A \"%(pf_i,Ip))\n", + "print(\"\\tAverage thyristor current = %.3f A\\n\\tRMS thyristor current = %.2f A\"%(val[0],val2))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) rms output voltage, Vl = 90 V\n", + "(b) Power output = 1012.5 W\n", + "(c) Since losses are neglected, input power = Output power = 1012.5 W\n", + "(d) input power factor = 0.6 lagging\n", + "(e)\tPeak thyristor current = 26.5125 A \n", + "\tAverage thyristor current = 3.038 A\n", + "\tRMS thyristor current = 7.95 A\n" + ] + } + ], + "prompt_number": 26 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.2, Page No.288" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# single phase half wave regulator parameter\n", + "\n", + "import math\n", + "#vaariable declaration\n", + "V = 150.0 # input voltage\n", + "Rl = 8.0 # load resistance\n", + "theta = 60.0 # thyristor firing angle in degrees\n", + "\n", + "# Calculations\n", + "Vm =math.sqrt(2)*V\n", + "Vm = math.floor(Vm*10)/10\n", + "#(a)\n", + "Vo = Vm*(math.cos(theta*math.pi/180)-1)/(2*math.pi)\n", + "#(b)\n", + "Vl = Vm*math.sqrt(((2*math.pi)-(theta*math.pi/180)+(math.sin(2*theta*math.pi/180)/2))/(4*math.pi))\n", + "Vl = math.ceil(Vl*100)/100\n", + "#(c)\n", + "Po = (Vl**2)/Rl\n", + "#(d)\n", + "I = Vl/Rl\n", + "va = V*I\n", + "pf = Po/va\n", + "#(e)\n", + "Iavg = Vo/Rl\n", + "\n", + "#Result\n", + "print(\"(a) Average output voltage = %.1f V\\n This is negative because only a part of positive half cycle appears at the output\"%Vo)\n", + "print(\" whereas the whole negative cycle appears at the output.\")\n", + "print(\"\\n(b) Vl = %.2f V\\n\\n(c) Power output = %.2f W\\n\\n(d) Input pf = %.2f lagging\"%(Vl,Po,pf))\n", + "print(\"\\n(e) Average input current = Average output current = %.2f A \"%(Iavg))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) Average output voltage = -16.9 V\n", + " This is negative because only a part of positive half cycle appears at the output\n", + " whereas the whole negative cycle appears at the output.\n", + "\n", + "(b) Vl = 142.46 V\n", + "\n", + "(c) Power output = 2536.86 W\n", + "\n", + "(d) Input pf = 0.95 lagging\n", + "\n", + "(e) Average input current = Average output current = -2.11 A \n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.3, Page No.293" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Single phase full wave regulator parameters\n", + "\n", + "import math\n", + "from scipy.integrate import quad\n", + "#variable declaration\n", + "V = 150.0 # input voltage\n", + "Rl = 8.0 # load resistance\n", + "theta = 60.0 # thyristor firing angle in degrees\n", + "\n", + "#calculations\n", + "Vm =math.sqrt(2)*V\n", + "Vm = math.floor(Vm*10)/10\n", + "#(a)\n", + "Vav = Vm*(math.cos(theta*math.pi/180)+1)/(math.pi)\n", + "#(b)\n", + "Vl = Vm*math.sqrt(((3.14)-(3.14/3)+(math.sin(2*theta*3.14/180)/2))/(2*3.141))\n", + "Vl = math.ceil(Vl*100)/100\n", + "#(c)\n", + "Io = Vl/Rl\n", + "Po = (Io**2)*Rl\n", + "#(d)\n", + "va = V*Io\n", + "pf = Po/va\n", + "#(e)\n", + "def f(x):\n", + " return math.sin(x)\n", + "wt_lower = math.pi/3 \n", + "wt_upper = math.pi\n", + "val = quad(f,wt_lower,wt_upper)\n", + "Iavg = val[0]*Vm/(2*math.pi*Rl)\n", + "\n", + "def g(x):\n", + " return math.sin(x)**2\n", + "wt_lower1 = math.pi/3\n", + "wt_upper1 = math.pi\n", + "val1 = quad(g,wt_lower1,wt_upper1)\n", + "Irms = (Vm/(Rl))*math.sqrt(val1[0]/(math.pi*2))\n", + "\n", + "#Result\n", + "print(\"(a) Average output voltage = %.2f V\\n \"%Vav)\n", + "print(\"\\n(b) Vl = %.2f V\\n\\n(c) Power output = %.2f W\\n\\n(d) Input pf = %.3f lagging\"%(Vl,Po,pf))\n", + "print(\"\\n(e) Average thyristor current= %.2f A\\n RMS thyristor current = %.2f A \"%(Iavg,Irms))\n", + "\n", + "#(b) For Vl calculation, value of pi is used different at different place. if math.pi is used then Vl = 134.53V\n", + "#(c) answer in the book is misprinted" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) Average output voltage = 101.27 V\n", + " \n", + "\n", + "(b) Vl = 134.52 V\n", + "\n", + "(c) Power output = 2261.95 W\n", + "\n", + "(d) Input pf = 0.897 lagging\n", + "\n", + "(e) Average thyristor current= 6.33 A\n", + " RMS thyristor current = 11.89 A \n" + ] + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.4, Page No. 294" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Single phase full wave ac voltage regulator\n", + "\n", + "import math\n", + "from scipy.integrate import quad\n", + "#variable declaration\n", + "V = 120.0 # input voltage\n", + "Rl = 10.0 # load resistance\n", + "alfa = 90.0 # thyristor firing angle in degrees\n", + "\n", + "#Calculations\n", + "Vm =math.sqrt(2)*V\n", + "#Vm = math.floor(Vm*10)/10\n", + "#(a)\n", + "Vo = Vm*math.sqrt(((math.pi)-(alfa*math.pi/180)+(math.sin(2*alfa*math.pi/180)/2.0))/(2*math.pi))\n", + "#Vo = 60*math.sqrt(2) \n", + "#(b)\n", + "Il = Vo/Rl\n", + "Po = (Il**2)*Rl\n", + "VA = Il*V\n", + "pf = Po/VA\n", + "#(c)\n", + "def f(x):\n", + " return math.sin(x)\n", + "wt_lower = math.pi/2\n", + "wt_upper = math.pi\n", + "val = quad(f,wt_lower,wt_upper)\n", + "Iav = Vm*val[0]/(2*math.pi*Rl)\n", + " \n", + "#(d)\n", + "Irms = Il/math.sqrt(2)\n", + "#(e)\n", + "Irmsl = Vo/Rl\n", + "\n", + "#Result\n", + "print(\"(a) RMS output voltage = %f V\\n(b) Input p.f. = %.3f lagging\\n(c) Average thyristor current = %.1f A\"%(Vo,pf,Iav))\n", + "print(\"(d) RMS thyristor current = %.1f A\\n(e) RMS load current =%f A\"%(Irms,Irmsl))\n", + "# Answer for \"RMS output voltage\" and \"RMS load current\" is wrong in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) RMS output voltage = 84.852814 V\n", + "(b) Input p.f. = 0.707 lagging\n", + "(c) Average thyristor current = 2.7 A\n", + "(d) RMS thyristor current = 6.0 A\n", + "(e) RMS load current =8.485281 A\n" + ] + } + ], + "prompt_number": 48 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.5, Page No. 294" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# rms output voltage and average power\n", + "\n", + "import math\n", + "R = 400.0 # load resistance\n", + "V = 110.0 # input voltage\n", + "alfa = 60.0 # firing angle\n", + "\n", + "# calculations\n", + "sqrt_2 =math.floor(math.sqrt(2)*1000)/1000\n", + "Vm = sqrt_2*V\n", + "Vo = Vm*math.sqrt(((math.pi)-(alfa*math.pi/180)+(math.sin(2*alfa*math.pi/180)/2.0))/(2*math.pi))\n", + "P = (Vo**2)/R\n", + "\n", + "#Result\n", + "print(\"RMS output voltage = %.2f V\\nAverage power = %.2f W\"%(Vm,P))\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "RMS output voltage = 155.54 V\n", + "Average power = 24.33 W\n" + ] + } + ], + "prompt_number": 58 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.6, Page No.294 " + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Firing angle\n", + "\n", + "import math\n", + "from numpy import poly1d\n", + "# variable declaration\n", + "a1 = 0.80 # % of maximum power value1\n", + "a2 = 0.30 # % of maximum power value2\n", + "\n", + "#Calculation\n", + "#After solving equation using taylor seris of X, we got following coefficient. \n", + "P1 = poly1d([128.0/math.factorial(7),0,-32.0/math.factorial(5),0,8.0/math.factorial(3),0,0,(-2*math.pi*(1-a1))], variable = 'x')\n", + "P2 = poly1d([128.0/math.factorial(7),0,-32.0/math.factorial(5),0,8.0/math.factorial(3),0,0,(-math.pi*1.525)], variable = 'x')\n", + "# hardcoded value used to match the answer to the book\n", + "x1 = P1.r[(P1.order+1)/2]*180/math.pi\n", + "x2 = P2.r[(P1.order+1)/2]*180/math.pi\n", + "alfa1 = x1.real\n", + "alfa2 = x2.real\n", + "print(\"(a) alfa1 = %.1f\u00b0\"%(x1.real))\n", + "print(\"(b) alfa2 = %.1f\u00b0\"%(x2.real))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) alfa1 = 60.5\u00b0\n", + "(b) alfa2 = 108.6\u00b0\n" + ] + } + ], + "prompt_number": 98 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.7, Page No.295" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# single phase AC regulator\n", + "\n", + "import math\n", + "#variable declaration\n", + "V = 230.0 # Single phase supply voltage\n", + "f = 50.0 # supply frequency\n", + "R = 15.0 # load resistance\n", + "alfa = 30.0 # firing angle\n", + "\n", + "#calculations\n", + "sqrt_2 =math.floor(math.sqrt(2)*1000)/1000\n", + "Vm = V*sqrt_2\n", + "alfa = alfa*math.pi/180 #degree to radians\n", + "#(a)\n", + "#-Theoretical\n", + "Iavg = Vm*(1+ math.cos(alfa))/(2*math.pi*R)\n", + "Irms = (Vm/R)*math.sqrt(((math.pi - alfa)/(4*math.pi))+(math.sin(2*alfa)/(8*math.pi)))\n", + "#(b)\n", + "a1 = (Vm)*(math.cos(2*alfa)-1)/(2*math.pi)\n", + "b1 = (Vm)*((math.pi-alfa)+((math.sin(2*alfa))/2.0))/(math.pi)\n", + "Va = math.sqrt(a1**2 + b1**2)\n", + "Vrms = Va/sqrt_2\n", + "Vrms = math.floor(Vrms*1000)/1000\n", + "#(e)\n", + "P = (Vrms**2)/R\n", + "#(f)\n", + "Vl = Vm*math.sqrt((math.pi -alfa+((math.sin(2*alfa))/2.0))/(2*math.pi))\n", + "Vl =math.floor(Vl*1000)/1000\n", + "#(g)\n", + "Pt = (Vl**2)/R\n", + "\n", + "#Result\n", + "print(\"(a) Average thyristor current = %.2f A\\n RMS thyristor current = %.3f A\"%(Iavg,Irms))\n", + "print(\"\\n(b) Amplitude of fundamental component of load voltage = %.1f V\\n RMS value = %.3f V\"%(math.floor(Va*10)/10,Vrms))\n", + "print(\"\\n(c) Since load is purely resistive, the load current rises from 0 to peak value instantaneously at alfa.\")\n", + "print(\" Hence max.(di/dt) is infinite\")\n", + "print(\"\\n(d) Maximum forward or reverse voltage across thyristor = %.2f V \"%(Vm))#math.pi =1.414\n", + "print(\"\\n\\n(e) Power delevered to load by fundamental component of load voltage = %.2f W\"%P)\n", + "print(\"\\n\\n(f) Load voltage = %.3f V\"%Vl)\n", + "print(\"\\n\\n(g) Total power output = %.2f W\"%Pt)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) Average thyristor current = 6.44 A\n", + " RMS thyristor current = 10.683 A\n", + "\n", + "(b) Amplitude of fundamental component of load voltage = 316.9 V\n", + " RMS value = 224.116 V\n", + "\n", + "(c) Since load is purely resistive, the load current rises from 0 to peak value instantaneously at alfa.\n", + " Hence max.(di/dt) is infinite\n", + "\n", + "(d) Maximum forward or reverse voltage across thyristor = 325.22 V \n", + "\n", + "\n", + "(e) Power delevered to load by fundamental component of load voltage = 3348.53 W\n", + "\n", + "\n", + "(f) Load voltage = 226.625 V\n", + "\n", + "\n", + "(g) Total power output = 3423.93 W\n" + ] + } + ], + "prompt_number": 26 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.8, Page No. 300" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Single phase full wave regulator\n", + "\n", + "import math\n", + "#Variable declaration\n", + "V = 150 # input voltage\n", + "R = 4 # resistance of the circuit\n", + "L = 22*10**-3 # Inductance\n", + "f = 50 # frequency \n", + "alfa = 60 # firing angle\n", + "\n", + "#Calculation\n", + "#(a)\n", + "w = 2*math.pi*f\n", + "theta = (180/math.pi)*math.atan(w*L/R)\n", + "beta = 180+alfa\n", + "#(b)\n", + "alfa = alfa*math.pi/180\n", + "Vm = V*math.sqrt(2)\n", + "Vm = math.floor(Vm*10)/10\n", + "Vl = Vm*math.sqrt((math.pi-((math.sin(2*beta*math.pi/180))/2)+((math.sin(2*alfa))/2.0))/(2*math.pi))\n", + "\n", + "#Result\n", + "print(\"(a)\\ntheta = %.0f\u00b0\\nAs alfa = theta, (beta-alfa) = Conduction angle= pi\\ntherefore, Beta = %d\u00b0\"%(theta,beta))\n", + "print(\"\\n(b)\\n Vl = %.0f V\"%Vl)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a)\n", + "theta = 60\u00b0\n", + "As alfa = theta, (beta-alfa) = Conduction angle= pi\n", + "therefore, Beta = 240\u00b0\n", + "\n", + "(b)\n", + " Vl = 150 V\n" + ] + } + ], + "prompt_number": 41 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.9, Page No. 300" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Single phase full wave regulator parameters\n", + "\n", + "import math\n", + "#Variable declaration\n", + "V = 230 # input voltage\n", + "f = 50 # frequency\n", + "R = 5 # resistance of the circuit\n", + "L = 20*10**-3 # inductance\n", + "\n", + "#Calculations\n", + "w = 2*math.pi*f\n", + "theta = (180/math.pi)*math.atan(R/(2*math.pi*f*L))\n", + "theta = math.ceil(theta*100)/100\n", + "Il = V/math.sqrt((R**2)+((w**2)*(L**2)))\n", + "Il = math.floor(Il*100)/100\n", + "P =R*Il**2\n", + "ipf = (P)/(V*Il)\n", + "\n", + "#Result\n", + "print(\"theta = %.2f\u00b0\"%theta)\n", + "print(\"\\n(a) The minimum value of firing angle(alfa) is theta. Therefore, range of firing angle is %.2f\u00b0 to 180\u00b0.\"%theta)\n", + "print(\"\\n(b) Conduction period of each thyristor is 180\u00b0.\")\n", + "print(\"\\n(c)\\n Load current = %.2f A\\n Power output = %.2f W\\n Input power factor = %.3f lagging\"%(Il,P,ipf)) " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "theta = 38.52\u00b0\n", + "\n", + "(a) The minimum value of firing angle(alfa) is theta. Therefore, range of firing angle is 38.52\u00b0 to 180\u00b0.\n", + "\n", + "(b) Conduction period of each thyristor is 180\u00b0.\n", + "\n", + "(c)\n", + " Load current = 28.64 A\n", + " Power output = 4101.25 W\n", + " Input power factor = 0.623 lagging\n" + ] + } + ], + "prompt_number": 52 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.10, Page No. 300" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Singal phase full wave regulator circuit parameters\n", + "\n", + "import math\n", + "from scipy.integrate import quad\n", + "#Variable declaration\n", + "R = 9 # resistance of the circuit\n", + "L = 0.03 # inductance\n", + "V = 240 # input voltage\n", + "f = 50 # frequency\n", + "alfa1 = 0 # firing angle case a\n", + "alfa2 = 60 # firing angle case b\n", + "alfa3 = 90 # firing angle case c\n", + "\n", + "#calculations\n", + "RbyL = R/L\n", + "alfa1 = 0*math.pi/180\n", + "alfa2 = 60*math.pi/180\n", + "alfa3 = 90*math.pi/180 \n", + "alfa4 = 13.7*math.pi/180 \n", + "sqrt_2 =math.floor(math.sqrt(2)*1000)/1000\n", + "w = 2*math.pi*f\n", + "Z_mag = math.sqrt((R**2)+(w**2)*(L**2))\n", + "Z_angle = math.floor((math.atan(w*L/R))*1000)/1000\n", + "x=math.floor(math.cos(Z_angle)*10000)/10000\n", + "\n", + "#(a)\n", + "Il = (V*math.sqrt(2)/Z_mag)#*math.sin((w/f)-(Z_angle*math.pi/180))\n", + "Il = math.floor(Il*100)/100\n", + "P = V*Il*x/sqrt_2\n", + "#b\n", + "k1 = math.ceil((Il*math.sin(alfa2-Z_angle))*100)/100\n", + "def h(t):\n", + " return ((V*sqrt_2*math.sin(w*t*alfa2))*((Il*math.sin((w*t)+(alfa2-Z_angle)))-(k1*math.e**(-RbyL*t))))\n", + "t2 = 0.00919\n", + "t1 = 0\n", + "val = quad(h,t1,t2)\n", + "P2 = val[0]/0.01\n", + "#c\n", + "k2 = math.floor((Il*math.sin(alfa3-Z_angle))*100)/100\n", + "def g(t):\n", + " return ((V*sqrt_2*math.sin(w*t*alfa3))*((Il*math.sin((w*t)+(alfa3-Z_angle)))-(k2*math.e**(-RbyL*t))))\n", + "t2 = 0.00733\n", + "t1 = 0\n", + "val2 = quad(g,t1,t2)\n", + "P3 = val2[0]/0.01\n", + "\n", + "angle1 = math.ceil((Z_angle*180/math.pi)*10)/10\n", + "angle2 = (alfa2-Z_angle)*180/math.pi\n", + "angle3 = (alfa3-Z_angle)*180/math.pi\n", + "\n", + "#Result\n", + "print(\"(a)\\nCurrent waveform--> i = %.2fsin(2pi*%dt-%.1f\u00b0) A\\nPower delivered = %.1f W\"%(Il,f,angle1,P))\n", + "print(\"(b)\\nCurrent waveform--> i = %.2fsin(2pi*%dt+%.1f\u00b0)-%.2fe^(-%dt) A\\nPower delivered = %.1f W\"%(Il,f,angle2,k1,RbyL,P2))\n", + "print(\"(c)\\nCurrent waveform--> i = %.2fsin(2pi*%dt+%.1f\u00b0)-%.2fe^(-%dt) A\\nPower delivered = %.1f W\"%(Il,f,angle3,k2,RbyL,P3))\n", + "#Power answer for (b) and (c) do not match with the book" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a)\n", + "Current waveform--> i = 26.04sin(2pi*50t-46.3\u00b0) A\n", + "Power delivered = 3053.6 W\n", + "(b)\n", + "Current waveform--> i = 26.04sin(2pi*50t+13.7\u00b0)-6.17e^(-300t) A\n", + "Power delivered = 3885.4 W\n", + "(c)\n", + "Current waveform--> i = 26.04sin(2pi*50t+43.7\u00b0)-17.99e^(-300t) A\n", + "Power delivered = 2138.7 W\n" + ] + } + ], + "prompt_number": 200 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.11, Page No. 304" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# Current and voltage ratings\n", + "\n", + "import math\n", + "# variable declaration\n", + "V = 415 # 3-phase input voltage\n", + "f = 50 # frequency\n", + "P = 20*10**3 # output power\n", + "sf = 1.5 # safety factor\n", + "\n", + "#Calculations\n", + "#(a)\n", + "Il = P/(math.sqrt(3)*V)\n", + "Irms = Il*sf\n", + "Irms = math.floor(Irms*100)/100\n", + "Vrms = V*sf\n", + "#(b)\n", + "Irms_thyristor = Il*sf/math.sqrt(2)\n", + "\n", + "#Result\n", + "print(\"(a)\\n Line Current = %.2f A\\n RMS rating of each triac = %.2f A\\n rms voltage rating = %.1f V\"%(Il,Irms,Vrms))\n", + "print(\"\\n(b)\\n Line Current = %.2f A\\n RMS rating of each thyristor = %.2f A\\n rms voltage rating = %.1f V\"%(Il,Irms_thyristor,Vrms))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a)\n", + " Line Current = 27.82 A\n", + " RMS rating of each triac = 41.73 A\n", + " rms voltage rating = 622.5 V\n", + "\n", + "(b)\n", + " Line Current = 27.82 A\n", + " RMS rating of each thyristor = 29.51 A\n", + " rms voltage rating = 622.5 V\n" + ] + } + ], + "prompt_number": 70 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.12, Page No.305" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# output parameters of 3-phase AC regulator\n", + "\n", + "import math\n", + "# variable declaration\n", + "V = 415 # 3-phase input voltage\n", + "R = 15 # load resistance per phase\n", + "alfa = 30 # firing angle\n", + "\n", + "#Calculations\n", + "#(a)\n", + "Vrms = V/math.sqrt(3)\n", + "Vl = math.sqrt(3)*(math.sqrt(2)*Vrms)*math.sqrt(((math.pi/6)-((alfa*math.pi)/(180*4))+(math.sin(2*alfa*math.pi/180)/8))/(math.pi))\n", + "Vl = math.floor(Vl*10)/10\n", + "#(b)\n", + "P = (3*Vl**2)/R\n", + "#(c)\n", + "Il = Vl/R\n", + "#(d)\n", + "VA = 3*Vrms*Il\n", + "ipf = P/VA\n", + "\n", + "#Result\n", + "print(\"(a) rms value of input phase voltage = %.1f V\\n Vl = %.1f V\"%(Vrms,Vl))\n", + "print(\"\\n(b) Output power = %.1f W or %.4f kW\\n\\n(c) Line Current = %.2fA\\n\\n(d) input p.f. = %.3f lagging\"%(P,P/1000,Il,ipf))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) rms value of input phase voltage = 239.6 V\n", + " Vl = 234.3 V\n", + "\n", + "(b) Output power = 10979.3 W or 10.9793 kW\n", + "\n", + "(c) Line Current = 15.62A\n", + "\n", + "(d) input p.f. = 0.978 lagging\n" + ] + } + ], + "prompt_number": 88 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "example 6.13, Page No. 305" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "# output parameters of 3-phase AC regulator\n", + "\n", + "import math\n", + "# variable declaration\n", + "V = 415 # 3-phase input voltage\n", + "R = 15 # load resistance per phase\n", + "alfa = 60 # firing angle\n", + "\n", + "#Calculations\n", + "#(a)\n", + "Vrms = V/math.sqrt(3)\n", + "#Vrms = math.floor(Vrms*10)/10\n", + "Vl = math.sqrt(3)*(1.414*Vrms)*math.sqrt(((3.141/6)-((alfa*3.141)/(180*4))+(math.sin(2*alfa*3.141/180)/8))/(math.pi))\n", + "#pi value = 3.141 to match the answer in the book\n", + "Vl = math.floor(Vl*100)/100\n", + "#(b)\n", + "P = (3*Vl**2)/R\n", + "#(c)\n", + "Il = Vl/R\n", + "#(d)\n", + "VA = 3*Vrms*Il\n", + "ipf = P/VA\n", + "\n", + "#Result\n", + "print(\"(a) rms value of input phase voltage = %f V\\n Vl = %f V\"%(Vrms,Vl))\n", + "print(\"\\n(b) Output power = %.2f W or %.5f kW\\n\\n(c) Line Current = %.3fA\\n\\n(d) input p.f. = %.2f lagging\"%(P,P/1000,Il,ipf))" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "(a) rms value of input phase voltage = 239.600362 V\n", + " Vl = 201.390000 V\n", + "\n", + "(b) Output power = 8111.59 W or 8.11159 kW\n", + "\n", + "(c) Line Current = 13.426A\n", + "\n", + "(d) input p.f. = 0.84 lagging\n" + ] + } + ], + "prompt_number": 109 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file -- cgit