{ "metadata": { "name": "" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Chapter 9: Single- and Two-Phase Motors" ] }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "Example 9.1, Page number: 459" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from __future__ import division\n", "from sympy import *\n", "import math\n", "\n", "#Variable declaration:\n", "Zmain=4.5+3.7j #main winding impedance(ohm)\n", "Zaux=9.5+3.5j #auxilliary winding impedance(ohm)\n", "f=60 #frequency(Hz)\n", "\n", "\n", "#Calculations:\n", "phy_main=math.degrees(math.atan(Zmain.imag/Zmain.real))\n", "phy=phy_main-90\n", "w=2*pi*60\n", "Xc=symbols('Xc')\n", "a=solve((3.5+Xc)/9.5-math.tan(math.radians(float(phy))), Xc)\n", "C=-1/(w*a[0])\n", "\n", "\n", "#Results:\n", "print \"The starting capacitance:\",round(float(C)*10**6,0), \"uF\"" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "The starting capacitance: 176.0 uF\n" ] } ], "prompt_number": 1 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "Example 9.2, Page number: 467" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from __future__ import division\n", "import cmath\n", "from math import *\n", "\n", "\n", "#Variable Declaration:\n", "R1_m=2.02 #resistance of main winding(ohm)\n", "X1_m=2.79 #resistance of main\n", "R2_m= 4.12 #Rotor resistance ref. to stator(ohm)\n", "X2_m=2.12 #Rotor reactance ref. to stator(ohm)\n", "Xm=66.8 #Magnetising reactance(ohm)\n", "s=0.05 #slip\n", "Pcu=24 #copper loss(W)\n", "Pw=13 #friction & windage loss(W)\n", "V=110 #line-to-line voltage(V)\n", "p=4 #no.of poles\n", "fc=60 #frequency(Hz)\n", "\n", "#Calculations:\n", "X22=X2_m+Xm\n", "Q2_m=X22/R2_m\n", "Rf=(Xm**2/X22)*(1/(s*Q2_m+1/(s*Q2_m)))\n", "Xf=(X2_m*Xm/X22)+Rf/(s*Q2_m)\n", "Zf=Rf+1j*Xf #forward field impedance(ohm)\n", "\n", "Rb=R2_m*(Xm/X22)**2/(2-s)\n", "Xb=(X2_m*Xm/X22)+Rb/((2-s)*Q2_m)\n", "Zb=Rb+1j*Xb #bachward field impedance\n", "Zt=0.5*(Zf+Zb)+R1_m+1j*X1_m\n", "I=V/abs(Zt) #Stator current(A)\n", "pf=cos(cmath.phase(Zt)) #power factor\n", "Pin=V*I*pf\n", "Pg_f=I**2*0.5*Rf #power absorbed by forward field(W)\n", "Pg_b=I**2*0.5*Rb #power absorbed by backward field(W)\n", "Pmech=(1-s)*(Pg_f-Pg_b)\n", "Pshaft=Pmech-(Pcu+Pw)\n", "ws=(2/p)*120*pi\n", "ns=(120/p)*fc\n", "n=(1-s)*ns #Rotor speed(rpm)\n", "wm=(1-s)*ws\n", "Tshaft=Pshaft/wm #shaft torque(Nm)\n", "eff=Pshaft/Pin\n", "\n", "#Results:\n", "print \"Stator current:\",round(I),\"A\", \"\\nPower factor:\",round(pf,3)\n", "print \"Power output:\",round(Pshaft),\"W\", \"\\nSpeed:\",n,\"rpm\"\n", "print \"Shaft torque:\",round(Tshaft,3),\"Nm\",\"Efficiency\",round(eff*100),\"%\"" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Stator current: 4.0 A \n", "Power factor: 0.621\n", "Power output: 147.0 W \n", "Speed: 1710.0 rpm\n", "Shaft torque: 0.823 Nm Efficiency 60.0 %\n" ] } ], "prompt_number": 2 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "Example 9.3, Page number: 474" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from __future__ import division\n", "import math\n", "import cmath\n", "\n", "\n", "#Variable declaration:\n", "f=60 #freq(Hz)\n", "omeag=2*pi*f\n", "s=0.05 #slip\n", "R1=0.534 #resistance of main winding(ohm)\n", "X1=2.45\n", "Xm=70.1\n", "R2=0.956\n", "X2=2.96\n", "Valpha=230\n", "Vbeta=210*cmath.exp(1j*80*pi/180)\n", "\n", "#Calculations:\n", "Vf = 0.5*(Valpha - 1j*Vbeta)\n", "Vb = 0.5*(Valpha + 1j*Vbeta)\n", "Zf=R1+1j*X1+1j*Xm*(R2/s+1j*X2)/(R2/s+1j*(X2+Xm))\n", "If=Vf/Zf\n", "Zb=R1+1j*X1+1j*Xm*(R2/(2-s)+1j*X2)/(R2/(2-s)+1j*(X2+Xm))\n", "Ib = Vb/Zb\n", "Ialpha=If+Ib\n", "Ibeta=1j*(If-Ib)\n", "Pgf=2*((Vf*(If.conjugate())).real-R1*abs(If)**2)\n", "Pgb=2*((Vb*(Ib.conjugate())).real-R1*abs(Ib)**2)\n", "Pmech=(1-s)*(Pgf-Pgb)\n", "\n", "\n", "#Results:\n", "print \"(a) Positive seq components:\", round(Vf.real,1)+1j*round(Vf.imag,1),\"V\"\n", "print\" Negative seq. components:\", round(Vb.real,1)+1j*round(Vb.imag,1),\"V\"\n", "\n", "print\"\\n(b) Positive stator currents:\",round(If.real,1)+1j*round(If.imag,1),\"A\"\n", "print\" Negative stator currnets:\",round(Ib.real,1)+1j*round(Ib.imag,1),\"A\"\n", "\n", "print\"\\n(c) Positive currents:\",round(Ialpha.real,1)+1j*round(Ialpha.imag,1),\"A\"\n", "print\" Negative currnets:\",round(Ibeta.real,1)+1j*round(Ibeta.imag,1),\"A\"\n", "\n", "print \"\\n(d) Power to forward field:\",round(Pgf,0),\"W\"\n", "print \" Power to backward field:\",round(Pgb,0),\"W\"\n", "print \" Pmech:\",round(Pmech,0),\"W\"" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "(a) Positive seq components: (218.4-18.2j) V\n", " Negative seq. components: (11.6+18.2j) V\n", "\n", "(b) Positive stator currents: (9.3-6.3j) A\n", " Negative stator currnets: (3.7-1.5j) A\n", "\n", "(c) Positive currents: (13-7.8j) A\n", " Negative currnets: (4.8+5.6j) A\n", "\n", "(d) Power to forward field: 4149.0 W\n", " Power to backward field: 15.0 W\n", " Pmech: 3928.0 W\n" ] } ], "prompt_number": 4 }, { "cell_type": "heading", "level": 3, "metadata": {}, "source": [ "Example 9.5, Page number: 483" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from __future__ import division\n", "%pylab inline\n", "import cmath\n", "from math import *\n", "\n", "#Variable declaration:\n", "Lmain=0.0806 #main winding inductance(H)\n", "Rmain = 0.58 #main winding resistance(ohm)\n", "Laux = 0.196 #auxilliary winding inductance(H)\n", "Raux = 3.37 #auxilliary winding resistance(ohm)\n", "Lr=4.7*10**-6 #rotor inductance(H)\n", "Rr=37.6*10**-6 #rotor resistance(ohm)\n", "Lmain_r=0.588*10**-3 #main inductance ref. to rotor(H)\n", "Laux_r = 0.909*10**-3 #aux inductance ref. to rotor(H)\n", "p=2 #poles\n", "Vo=230 #terminal voltage(V)\n", "w=120*pi #angular frequency(Hz)\n", "C=35*10**-6\n", "Prot=40 #Windage losses(W)\n", "Pcore=105 #Core loss(W)\n", "n=3500 #rpm\n", "\n", "\n", "#calculations and Results:\n", "Xc=-1/(w*C)\n", "speed=[0]*102\n", "for cal in range(1,3,1):\n", " if cal==1:\n", " mmax=2\n", " else:\n", " mmax=102\n", " for m in range(1,mmax,2):\n", " if cal==1:\n", " speed[m-1]=3500\n", " else:\n", " speed[m-1]=3599*(m-1)/100\n", " \n", " ns=(2/p)*3600\n", " s=(ns-speed[m-1])/ns\n", "\n", "#for part (a):\n", " Kplus=s*w/(2*(Rr+1j*s*w*Lr))\n", " Kminus=(2-s)*w/(2*(Rr+1j*(2-s)*w*Lr))\n", " A1=Lmain-1j*Lmain_r**2*(Kplus+Kminus)\n", " A2=Lmain_r*Laux_r*(Kplus-Kminus)\n", " A3=Laux-1j*Laux_r**2*(Kplus+Kminus)\n", " M=[[0]*2,[0]*2]\n", " M[0][0]=Rmain + 1j*w*A1\n", " M[0][1] = 1j*w*A2;\n", " M[1][0] = -1j*w*A2;\n", " M[1][1] = Raux + 1j*Xc+ 1j*w*A3\n", " V=[[Vo],[-Vo]]\n", " M1=inv(M)\n", " I=dot(M1,V)\n", " Imain=I[0][0]\n", " Iaux=I[1][0]\n", " Is=Imain-Iaux\n", " magImain=abs(Imain)\n", " angleImain=math.degrees(cmath.phase(Imain))\n", " magIaux=abs (Iaux)\n", " angleIaux=math.degrees(cmath.phase(Iaux))\n", " magIs=abs(Is)\n", " angleIs=math.degrees(cmath.phase(Is))\n", " Vcap=Iaux*Xc\n", " magVcap=abs(Vcap)\n", " \n", " #for part (b):\n", " Tmech=[0]*102\n", " Pshaft=[0]*102\n", " Tmechl = (Kplus-Kminus).conjugate()\n", " Tmechl=Tmechl*(Lmain_r**2*Imain*((Imain).conjugate())+Laux_r**2*Iaux*((Iaux).conjugate()))\n", " Tmech2 = 1j*Lmain_r*Laux_r*((Kplus+Kminus).conjugate())\n", " Tmech2 = Tmech2*((Imain).conjugate()*Iaux-Imain*((Iaux).conjugate()));\n", " Tmech[m-1] = (p/2)*(Tmechl+Tmech2).real\n", " Pshaft=((2/p)*(1-s)*w*Tmech[m-1])-Prot\n", " \n", " #for part (c):\n", " Pmech=[0]*102\n", " Pmain = (Vo*(Imain.conjugate())).real\n", " Paux = (-Vo*(Iaux.conjugate())).real\n", " Pin = Pmain+Paux+Pcore\n", " eta = Pshaft/Pin;\n", " if cal==1:\n", " print \"part (a):\"\n", " print \"\\nImain=\",round(magImain,1),\"A at an angle\",round(angleImain,1),\"degrees\"\n", " print \"\\nImain=\",round(magIaux,1),\"A at an angle\",round(angleIaux,1),\"degrees\"\n", " print \"\\nImain=\",round(magIs,1),\"A at an angle\",round(angleIs,1),\"degrees\"\n", " print \"\\nVcap=\",round(magVcap,0),\"V\"\n", " print \"\\npart (b):\"\n", " print \"\\nTmech=\",round(Tmech[0],2),\"Nm\"\n", " print \"\\nPshaft=\",round(Pshaft),\"W\"\n", " print \"\\npart (c):\"\n", " print \"\\nPmain=\",round(Pmain,0),\"W\"\n", " print \"\\nPaux=\",round(Paux,0),\"W\"\n", " print \"\\nPin=\",round(Pin,0),\"W\"\n", " print \"\\nEfficiency=\",round(eta*100,1),\"%\"\n", " else:\n", " \n", " plot(speed,Tmech,'g.')\n", " xlabel('speed (rpm)')\n", " ylabel('Tmech (Nm)')\n", " title('Electromagnetic torque vs speed')\n", " show()" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Populating the interactive namespace from numpy and matplotlib\n", "part (a):" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "\n", "\n", "Imain= 15.9 A at an angle -37.6 degrees\n", "\n", "Imain= 5.2 A at an angle -150.8 degrees\n", "\n", "Imain= 18.5 A at an angle -22.7 degrees\n", "\n", "Vcap= 394.0 V\n", "\n", "part (b):\n", "\n", "Tmech= 9.75 Nm\n", "\n", "Pshaft= 3532.0 W\n", "\n", "part (c):\n", "\n", "Pmain= 2893.0 W\n", "\n", "Paux= 1043.0 W\n", "\n", "Pin= 4041.0 W\n", "\n", "Efficiency= 87.4 %\n" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "WARNING: pylab import has clobbered these variables: ['Polygon', 'poly', 'sign', 'ldexp', 'hypot', 'flatten', 'conjugate', 'diff', 'tan', 'Circle', 'roots', 'plot', 'isnan', 'eye', 'trace', 'fabs', 'floor', 'diag', 'invert', 'nan', 'modf', 'sqrt', 'frexp', 'source', 'add', 'degrees', 'take', 'var', 'zeros', 'pi', 'log10', 'plotting', 'product', 'seterr', 'power', 'multinomial', 'copysign', 'transpose', 'expm1', 'ceil', 'test', 'beta', 'ones', 'isinf', 'sinh', 'vectorize', 'cosh', 'mod', 'cos', 'prod', 'e', 'f', 'tanh', 'det', 'radians', 'sin', 'binomial', 'solve', 'log', 'fmod', 'reshape', 'exp', 'trunc', 'log1p', 'gamma', 'interactive']\n", "`%pylab --no-import-all` prevents importing * from pylab and numpy\n" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEZCAYAAAB8culNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtYlGX+P/D3ABIqZ5QRRcWQRI6SWH09Nba6yphZKSmK\na2qx62YHzdr8WYLuplm7lXnIzNO1i6ZkbbumspmG5pYoqHnIpRVFwUTBQ6Kc4fP7g4tZxoezc3iG\neb+ui0vmnmee+cyNzIf7/tz3MxoRERAREdXhYO0AiIhIfZgciIhIgcmBiIgUmByIiEiByYGIiBSY\nHIiISIHJoY3YuHEjhgwZYu0wbJper8ff/vY3a4dBzcD/7+bH5GBDAgIC0KFDB7i5uRm+XnjhBZOd\nPycnBw4ODqiurjbZOdUqKSkJU6ZMMWrbuXOnoq05HBwccPbsWVOFRqQKTtYOgJpPo9Hgyy+/xCOP\nPGLW52lsX2RVVRUcHR3N+vy2qLV7SSsrK+HkxF9DUh+OHNqoY8eOYciQIXB3d0fPnj3x17/+1XDf\n7du38bvf/Q6+vr5wd3fHoEGDUFpaiqFDhwIAPD094e7ujoMHD2Ljxo0YNGgQ5syZA19fXyxatAjX\nr1/H+PHj4eHhAa1Wi/nz5xveHOse7+3tjd69e+O7777Dhg0bEBAQAC8vL3z88ceGWLZv346IiAi4\nu7tDq9XitddeM3odH330Ebp06QKtVos//elPCAgIwN69ewHU/PX/1FNPYerUqfDw8EDv3r3x/fff\nGx6bk5MDvV4PT09P+Pn5YenSpQCA1NRULFmyBFu3boWbmxuioqIAADqdDuvWrTM8/r333kOvXr3g\n5uaG4OBgHDlyRNHPtX0WGRkJNzc3fPrppwCA999/H/7+/nB3d8evf/1rnD9/3vAYBwcHrFq1Cn36\n9EFwcDAAYOHChfD29kaPHj2wfv16o9HInXHdOaXS2M+6rq1bt2LAgAFGbe+99x7Gjh0LAPjiiy8Q\nFBQEV1dXdO3aFW+//Xa95zl9+jQGDhwIV1dX+Pj4IDY21ui1LV++HEFBQXB3d8esWbOMRqLLly9H\nQEAA3N3d8fDDDyM7O7tZr+Py5csYPnw43Nzc8NBDDxk9jsxEyGYEBATI119/Xe99GzZskMGDB4uI\nyPXr18XX11eSk5NFROTUqVPi4+MjmZmZIiIydepUGTVqlBQWFoqIyOHDh6WsrExycnJEo9FIVVWV\n0XmdnJxk7dq1IiJSWloq48aNkwkTJkhJSYn8/PPPEhISIsuXLzc6vva5FyxYIN26dZOXXnpJKisr\nZc+ePdK+fXspKioSEZH9+/dLVlaWiIicPn1aunbtKp988omIiGRmZoq7u7tkZGRIVVWVzJ8/X9q1\nayd79uwREZHExERxcXEx9Mm8efPk/vvvFxGRyspKCQ4OliVLlkhVVZXk5ubKvffeK3//+99FRCQp\nKUmmTJli1Ic6nU7WrVtneB09e/aUkydPiojI+fPn5cKFC/X2vUajkezsbMPt7du3i6+vr5w+fVoq\nKytl7ty50r9/f6Pjx4wZI0VFRVJWViaff/65dO3aVbKzs6WsrEymTp1qdM66cbX0Z11XcXGxuLm5\nyX//+19DW3R0tGzdulVERLy9veXAgQMiIlJUVCQ//PBDva/3ySeflMWLF4uISEVFhaSnpxu9tlGj\nRklRUZHk5+dLSEiIfPDBByIismnTJgkKCpKzZ8+KiMiSJUukX79+zXodjz32mEyZMkXKy8vlv//9\nr3Tv3l2GDBlSb3xkGkwONqRnz57i6uoqnp6ehq/aN+26bxgbN25U/OIkJCTIvHnzpKSkRJydneU/\n//mP4vznzp2rNzkEBgYabhcXF0u7du3kzJkzhrb169fLQw89ZDg+KCjIcN/JkydFo9HIlStXDG2d\nO3eWjIyMel/jyy+/LDNnzhSRmjf7qVOnGu4rKyuT9u3bGyWHESNGGO4/deqUODk5iYhIWlqa9OjR\nw+jcixcvlri4OMNj4+Pjje6v+yY8ZMgQWb16db0x3unO5DBp0iR5/fXXDbdLSkrExcXFkAQ1Go3h\nTVhEJC4uThYsWGC4XZukm5McGvtZ1yc+Pl4WLVokIiI//fSTuLm5SUlJiYiI9OjRQ9asWSM3b95s\n9PX+5je/kd/+9rdy8eLFevti7969httr166VQYMG1fs6qqqqpEOHDpKVldXo6yguLhYnJydDUhGp\nSe61fUDmwWklG6LRaPCPf/wD169fN3zNmDFDcVxeXh7S09Ph5eVl+Nq8eTOuX7+Oa9euoaKiAvfe\ne2+zn9fPz8/w/dWrV1FZWYkePXoY2rp3747Lly8bbmu1WsP399xzDwCgc+fORm1lZWUAgG+//RaD\nBg2Ct7c3vLy8sHLlSty+fRsAUFBQgK5duxoe5+zsjE6dOhnFVve5OnTogKqqKlRXVyMvLw8///yz\nUR8sWbIEN27caNZrzs/Pb1Ef1XXlyhWj/nFxcUGnTp2M+qhunxYUFMDf399wu1u3bs1+rsZ+1vWZ\nNGkSPvnkEwDA5s2b8cQTT8DFxQUAkJKSgn/+85/o2bMnBg8ejG+//bbec7z11lsoLy/HgAED0Ldv\nX6xZs8bo/jtfS+3rzsvLw4svvmiI08fHx/D6m/o/W1VV1eo+otZhJawN8vPzw/Dhw7Fjxw7FfaWl\npXB2dsbZs2fRp08fo/s0Gk2T5/bx8YGjoyPOnz+P3r17AwByc3PRpUuXVsUaFxeH+fPn49lnn4WT\nkxNeeeUVw5uJr68vLl68aDi2rKwMhYWFzTpvly5dcN999+HUqVP13t/Ua+3atWurVyBptVqjGkNp\naSkKCwuNElldvr6+yMvLM9yu+z1QkxRrEyZQk6BrNfazrs/w4cNRUFCAH374AVu2bMH7779vuO/B\nBx/E9u3bUV1djeXLl+Opp57CpUuXFOfw8/PD+vXrAQDff/89hg0bBp1Oh/vuu88Qf1BQkOH72v8b\nfn5+WLJkCcaPH684Z1ZWVoOvo6SkBI6OjsjLy0OvXr3q7SMyPY4cbIw0Y1XM448/jmPHjmHbtm2G\nv6SPHj2KrKwsuLi4IC4uDnPmzMHVq1chIjh8+DDKy8vh6ekJjUaDc+fONXju9u3b47HHHsMbb7yB\n0tJSXLp0Ce+++y7i4uJa9XqKi4vRsWNHODk54ejRo9i0aZPhvieffBJ///vfceTIEVRVVeHNN99E\nZWVls8778MMPo7q6GitWrEB5eTlEBFlZWYaiso+PD3Jzcxvsz2nTpmHp0qX48ccfAdQUt3Nzc+s9\n1tvb26jPJkyYgLVr1+I///kPKisrsWDBAoSGhhrePO80fvx4rFu3DmfPnkVZWRkWLVpkdH9kZCQ+\n//xzlJSU4Pz580YF/cZ+1vVp164dYmNjMXfuXFy/fh0jRowAAFRUVCAlJQW3b9+Gg4MDXF1d4eBQ\n/9vDF198gfz8fACAu7s7HBwcjJLtn//8Z9y6dQuXL1/GsmXL8NRTTwEAEhISsHjxYpw5cwYAcOvW\nLXzxxRcAgCeeeKLB19G+fXvo9XosXLgQ5eXlyM7OxoYNG5r1xwy1HpODjRkzZozRPodx48YBqPlL\nuPaXxdvbG6mpqVi9ejW8vb3h4+OD2bNno7S0FACwcuVK+Pv7o0+fPvD09MScOXMgIvDw8MCcOXMQ\nHR0Nb29vpKenG5231po1a1BeXg6tVovIyEg8+uijeP755xVx1Grsl3jFihWYN28ePDw8sGDBAqO/\nKvv374+lS5ciJiYGXbt2hbOzM/z8/AxLaRt7LicnJ/zrX//Cnj17oNVq4enpid/85jeG6ZbY2FiU\nlJTAw8MD0dHRirimTp2KWbNmISYmBm5ubtDr9UZ/sdf1+uuvY8KECfDy8sK2bdswZswYvPrqq/jV\nr34FLy8vHD16FJ999lmD/fHEE09gxowZ6N+/P4KCghSbu1555RVUVVWhU6dOiI+PR1xcXLN/1vWZ\nNGkS9uzZg9jYWKMEsHbtWvj7+6Njx45YsWKFUaKu68CBA4iKikLHjh2h1+vx9ttvG0YKABATE4N+\n/fohKCgIQ4cOxaxZswAA8fHxSEhIQExMDNzd3dGnTx9DcvDy8mr0dXz00UfIzc2Fj48PJk+ejKlT\npzb4+sg0NNKcP0VbITc3F5MnT8b169dRXl6OGTNm4NVXX8W1a9cwYcIEXL58GX5+fti6dSs8PT3N\nEQK1MSUlJfDy8sIPP/ygmBJraxwcHHDmzJlW1z2sxVbjJiWzjRycnZ2xatUqnDhxApmZmVi7di1+\n+OEHJCYmYvTo0Th+/DhiYmKQmJhorhCoDUhNTUVpaSnKysowb948BAQEtPnEQKQGZksOWq0WYWFh\nAABXV1dERETg4sWLRpcoiI+Pb3YhjexTSkoKunTpAm9vb2RmZmLbtm3WDskibHU+3VbjJiWzTSvV\nlZOTg4cffhgnTpyAv78/bt68abjP3d3d6DYREVmf2QvSt27dwvjx47Fs2TK4u7ub++mIiMgEzLrP\noaKiAuPGjcPkyZPx+OOPA6jZDFVYWIhOnTqhoKAAvr6+isf17t2b104hImqhwMBAw1Lhu2W2kYOI\nYMaMGQgJCcHs2bMN7Xq9HsnJyQCA5ORk6PV6xWOzs7MhNZf2UPVXYmKi1WNgnIyTcTLG2i9T/lFt\ntpHDv//9byQnJyMiIsJw1cslS5Zg4cKFmDBhAtavX48uXbogJSXFXCEQEVErmS05DB48uMEPjdm9\ne7e5npaIiEyAO6Tvgk6ns3YIzcI4TYtxmpYtxGkLMZqaRZaytpRGo4EKwyIiUjVTvndy5EBERApM\nDkREpMDkQERECkwORESkwORAREQKTA5ERKTA5EBERApMDkREpMDkQERECma9ZDcRkSUlbE/AT1d/\nQod2HbB53GZ4uvDz6VuLyYGIbE5DSeCnqz9h3/l9hmNSYnnV59bitBIR2ZzaJLDrzC4kbE8wtHdo\n1wEAEN01GmvGrLFWeG0CkwMR2ZyGksDmcZsRGxKL3VN2c0rpLvGqrERkc26U3kDC9gSsGbOGSaAO\nU753MjkQkWqxwNwyvGQ3EdmFhmoLZH5MDkSkWiwwWw+nlYhItUxVW7CX6SnWHIiIWkC3UWfY/xAb\nEttm9z+w5kBE1AKcnmo5jhyIqM2zl6WvnFYiIiIFU7538tpKRGRV9lIstjWsORCRVXEvgzoxORCR\nVbFYrE6sORCRVdlLsdgSWJAmIiIF7nMgIiKzYnIgIiIFJgciIlLgPgcisgjuZ7AtHDkQkUVwP4Nt\nYXIgIovgfgbbwqWsRGQR3M9gftznQERECtznQEREZsXVSkRkt7iCqmEcORCR3eIKqoYxORCR3eIK\nqoaxIE1EdqutraDiaiUiUi3O41sPVysRkWpxHr9tYHIgIpPiPH7bwGklIjKptjaPb0tYcyAiIgWb\nqTlMnz4dWq0W4eHhhrakpCT4+/sjKioKUVFRSE1NNWcIRETUCmZNDtOmTVO8+Ws0GsyZMwdHjx7F\n0aNHMWrUKHOGQERErWDW5DBkyBB4eXkp2jllRESkblZZrbRy5Ur07dsX8fHxuHbtmjVCICKiRlj8\nwnvPPfccFixYAKCm/vDCCy8gOTlZcVxSUpLhe51OB51OZ6EIiag5uNnN+tLS0pCWlmaWc5t9tVJO\nTg7GjBmDEydOKO77+eefMWzYMGRlZRkHxdVKRKqn26jDvvP7AACxIbFIiU2xckRkM6uV6nPlyhXD\n95999hlCQ0MtHQIRmQA3u7VtZh05xMXFYd++fSgsLIRWq8XChQvxzTff4Pjx4ygvL0fPnj2xbt06\ndOvWzTgojhyIVI+b3dSHm+CIiEjBpqeViIhI/ZgciIhIgcmBiIgUmByIiEjB4pvgiMi2cLObfeLI\ngYgaxU92s09MDkTUKG52s0/c50BEjeJmN9vBTXBERKTATXBERGRWTA5ERKTA5EBERApMDkREpMBN\ncEQEgJvdyBhHDkQEgJvdyBiTAxEB4GY3MsZ9DkQEgJvd2gJugiMiIgVugiMiIrPiaiUiojtw5RZH\nDkRECly5xZEDkd3hX8VN48qtJgrSV65cwaeffor9+/cjJycHGo0GPXv2xNChQxEbGwtfX1/zBMWC\nNJHZ6DbqsO/8PgBAbEgsUmJTrByR+tjqyi2LrFaaMWMGsrOzERMTgwceeAB+fn4QEVy6dAmHDh1C\namoqevfujbVr15okEKOgmByIzEa/SY9dZ3Yhums0dk/ZbVNvftQ4iySH48ePIyIiotEHN+eYVgXF\n5EBkNrb6VzE1jfsciIhIwaL7HD7//HOEhYXBw8MDbm5ucHNzg7u7u0menIiI1KnJkUOPHj3w5Zdf\nIiwsDA4Olln5ypEDEVHLmfK9s8mlrL1790Z4eDg0Go1JnpCIiNSvyZFDeno6FixYAJ1OB2dn55oH\naTSYM2eO+YLiyIHornE/g/2x6Mjh9ddfh5ubG0pLS1FeXm6SJyUi86vd5QvUJAruZ6CWaDI5XL58\nGbt377ZELERkQtzlS3ejyQpzTEwMkwORDdo8bjNiQ2K50Y1apcmag6urK4qLi+Hs7Ix27drVPEij\nwc2bN80XFGsOREQtxk1wRESkYJGC9IULFxp9YI8ePUwSABERqU+DI4ewsLB69zYUFBSgoKAAVVVV\n5guKIweiZuOSVaplkZHDyZMnjW7n5OTgrbfewtdff4358+eb5MmJ6O5xySqZQ5OrlX766Sc8/fTT\nGDVqFPr374/Tp0/j+eeft0RsRNQMXLJK5tDgtNKJEyfw5ptv4tSpU3j11VcxadIkODo6WiYoTisR\nNRsvwU21LLJaydHREf7+/nj00UcVF9zTaDT44IMPTBJAvUExORARtZhFag7r1q0zPFldIsKL8BER\ntXHc50BE1EZY5MN+pk+fjsOHDzf4wPT0dEybNs0kQRBR0xK2J0C3UQf9Jj1ulN6wdjjUxjU4rTR7\n9my88847OHjwIPr06QM/Pz+ICPLz85GVlYWBAwdi7ty5loyVyK5xySpZUpPTSmVlZTh69CjOnz8P\njUaDnj17IjIyEi4uLuYLitNKRAr6TXrsOrML0V2jeTE9qhevrURkh7hklZpikZqDKUyfPh1arRbh\n4eGGtmvXrmHEiBGIiIjAyJEjceMG506JmsPTxRMpsSlMDGQRZk0O06ZNQ2pqqlFbYmIiRo8ejePH\njyMmJgaJiYnmDIGIiFrB7NNKOTk5GDNmDE6cOAEACAwMxKFDh+Dj44PCwkI89NBDOHPmjHFQnFYi\nO8YL6VFrWfQzpE+ePIk///nPyM3NRXV1tSGAvXv3tuoJCwoK4OPjAwDo1KkTrly50qrzELVVXJVE\natBkchg3bhxeeuklzJw503BtJUvskE5KSjJ8r9PpoNPpzP6cRGrAC+lRc6WlpSEtLc0s525yWumB\nBx7AoUOHWv0E9U0rpaeno1OnTigoKMD//d//cVqJqA6uSqLWsshqpWvXruHq1avQ6/VYvXo1Ll26\nhGvXrhm+Wkuv1yM5ORkAkJycDL1e3+pzEbVFXJVEatDgyCEgIKDR6aNz5841efK4uDjs27cPhYWF\n0Gq1WLRoEcaOHYsJEybg8uXL6NKlC1JSUuDpafxLwJEDEVHLcRMcURvAVUlkahbdBPfBBx/gl19+\nMdz+5ZdfsGLFCpM8OZE9q12VtOvMLiRsT7B2OERGmkwO69atg4eHh+G2h4cH1q5da9agiOwBVyWR\nmjWZHMrLy41uiwhKS0vNFhCRvdg8bjNiQ2J5ET1SpSZrDs899xyuXr2KZ599FiKCjz/+GD4+Pli1\napX5gmLNgYioxSxakK6srMTy5cuxZ88eAMCIESMwa9Ysw4Y4c2ByICJqOYuvVioqKsKFCxcQGhpq\nkidtCpMDtSVclUSWYtHVSp9++imioqIwevRoADXXWqr9noiaxlVJZIuaTA5JSUnIyMiAl5cXACAs\nLAy5ublmD4yoreCqJLJFTSYHJycnxQ7myspKswVE1NZwVRLZoiavyhoSEoJNmzahsrIS586dw6pV\nqzBgwABLxEbUJtReK4nIljQ5cvj444+RmZkJEcGYMWNQXV2NDz/80BKxERGRlfDaSkQmwlVJZG0W\n/SS47777DosXL1Z8Etzx48dNEgBRW8FPcKO2pMnkMHnyZCxbtgxhYWFwcGhyForIbnFVErUlTU4r\nDR06FPv377dUPAA4rUS2iZ/gRtZm0R3Su3fvRkpKCh555BE4OzsbAnjyySdNEkC9QTE5EBG1mEVr\nDhs3bkRWVhbKy8uNppXMmRyI1IyFZ7IHTSaHzMxMnD59utGPDCWyJyw8kz1ossI8aNAgZGVlWSIW\nIpvAwjPZgwZrDpWVlXByckJwcDCys7PRq1cv3HPPPTUPMvNSVtYcSM1YeCa1skhB+v7778eRI0eQ\nk5NT7wMDAgJMEkC9QTE5kAqwtkC2xiIF6donMGcSIFIz1hbInjWYHAoKCvDuu+/Wm4U0Gg3mzJlj\n1sCIrI21BbqTPY0mG0wOVVVVKCoqsmQsRKqyedxm1hbIiD2NJhusOURFReHo0aOWjgcAaw5EpE76\nTXrsOrML0V2jVfn5HBb9mFCiti5hewJ0G3XQb9LjRukNa4dDKmZPH9zU4Mjh6tWr8PHxsXQ8ADhy\nIMvSbdQZpgpiQ2Lb9FQBtW0WGTlYKzEQWRoLz0RK/LAfsnvc1EZthUWvymoNTA5kDva0DJHsEwvS\nRK1Quwxx15ldSNieYO1wiFSNyYHsBmsLRM3HaSWyG6wtUFvHmgNRI1hbIHvFmgNRI1hbILp7TA7U\n5rC2QHT3OK1EbQ5rC2SvWHMgAmsLRHdizYEIrC0QmROTA9ks1haIzIfTSmSzWFsgMsaaA9kV1haI\nmoc1B7IrrC0QWR6TA6keawtElsdpJVI91haImoc1B2qTWFsgujusOVCbxNoCkXo4WeuJAwIC4O7u\nDkdHR7Rr1w6HDh2yViikEqwtEKmH1aaVevXqhczMTHh7eyvu47RS29bQ9BFrC0R3p81MKzEB2KeG\npo88XTyREpvCxECkAlZLDhqNBiNGjEBERARWrFhhrTDICjh9RKR+Vqs5HDx4EL6+vigoKMCoUaMQ\nHByM4cOHG+5PSkoyfK/T6aDT6SwfJN2VhqaPNo/bzOkjIhNIS0tDWlqaWc6tiqWsS5YsAQDMmzcP\nAGsObYVuow77zu8DAMSGxCIlNsXKERG1bTZfcyguLkZxcTEA4Pbt20hNTUVoaKg1QiEz4vQRke2y\nysjh3LlzePzxx6HRaFBcXIyJEydi0aJF/wuKI4c2gauPiCyLO6RJVbizmUgdbH5aidoW7mwmanuY\nHOiusbZA1PZwWomajTubidSNNQeyCi5NJVI31hzIKjh9RGQ/OHIgBU4fEdkmTiuRWXH6iMg2mfK9\n02rXViLra2iEwOkjImLNwY41tD9h87jNiA2Jxe4puzl9RGSnOHKwYw2NEGo/V4GI7BdrDnaABWYi\n+8CCNLUIC8xE9oH7HKhFWGAmopbiyKEN4fQRkX3jtBLVi9NHRPaN+xzsHPcnEJG5seZgg7g/gYjM\njSMHG8T9CURkbqw5qBgLzETUEixI2wkWmImoJViQbmNYYCYitWFBWgVYYCYiteHIQQVYYCYitWHN\nwYJYYCYic2JB2kaxwExE5sSCtMqxwExEto4FaTNggZmIbB1HDmbAAjMR2TrWHO4CC8xEpCYsSKsE\nC8xEpCYsSFsYC8xEZG9YkG4GFpiJyN5w5FBHS0cILDATUVvFkUMdHCEQEdWwy5EDRwhERI2zy5ED\nRwhERI2zy5EDRwhERI1r0/scuEmNiOwJN8E1EzepEZE94Sa4OhoaHQDcpEZE1Fo2X5BuqLgMsMBM\nRNRaNjNyaM0lLFhgJiJqHZsZOXD5KRGR5ag2Oeg36XGj9IbhdlPLT5kYiIhMR7XJgSMEIiLrUe1S\n1ug10UwEREQtYMqlrFYZOaSmpiI8PBwhISFYunRpvccwMRARWY/Fk0NZWRlmzpyJ1NRUHD9+HNu2\nbcPRo0cVx9lCYkhLS7N2CM3COE2LcZqWLcRpCzGamsWTQ3p6OkJDQ9GtWzc4OTlhwoQJ2LFjh6XD\nMAlb+Q/DOE2LcZqWLcRpCzGamsWTQ15eHrp372647e/vj7y8PEuHQURkMgnbE6DbqFOssrRlFi9I\nf/LJJ9i/fz8+/PBDAMCWLVuQlpaG1atX/y8ojQZIqvleEmvC0yzUGO53hCMqEysV7U5wQkVixf/O\nU+e+dmiH8sRyRbsznFGWWKZovwf3oDSxVNHuAheUJJbUtOs0wLCa9vZoj+LEYsXxHdABtxNvK9o7\najri1oJbinZXB1cUvVGkaHdzdMPN128q2t2d3PHL/F8U7R7tPHDj/91QxOnl7IVr864pjve5xweF\nrxUq2ju7dMaVP1xRtGvba5H/ar6i3a+DH35+5WdFe7eO3ZA3N0/R3t21Oy68fEERZ4BbAM7NOac4\nPtA9EGdmn1G03+dxH7JeylK09/Xqix9f+FHRHuYThhOzTijaIztH4tjvjyna+2v7I+N3GYo4H/R7\nEAcTDiqOH9RtEA48c0DRruuuwzfTv1G0jwgYga+mfqVo1wfqsSN+h6J97H1j8UXcF4r22L6xSHkq\nRRHn5LDJSB6XrDh+euR0rHt8naJ9Zv+ZWPXoKkX77Adn491R7yra5w2ah8XDFyvaFz28CG/o3lC0\nvzP8HcwdNFcR58qYlfj9A783HFf3MRse24Cno55WtKeMT0FsaCwAwGGhAwQ17xk7J+1ETFAMAMD5\nj86oqK6ABhrsn7Yfg3sMBgB4vuWJW+W34KBxQEZCBiK0EQCA4BXByL+Vj3aO7ZDxbAY2vL8BSUlJ\nDW7GVct13Gz6wnvffvstli5dii+//BIA8M4776C8vBzz58//X1DeGuC6JaMiIrJ9gYGBOHPmjEnO\nZfHkUFpaiuDgYPz73/+Gr68vBg4ciI8++gj333+/JcMgIqJGWPzaSi4uLvjwww8xcuRIVFdXY8qU\nKUwMREQqo8pNcEREZF2qu3xGczbIWVJAQAAiIiIQFRWFBx54AABw7do1jBgxAhERERg5ciRu3Pjf\n6oQlS5aWhx3EAAAJXklEQVQgJCQE4eHh+Oqrr8wS0/Tp06HVahEeHm5oa01MmZmZiIqKQmhoKF58\n8UWLxJmUlAR/f39ERUUhKioKu3btsnqcubm5GDp0KMLDw9GnTx+8/fbbANTXpw3FqbY+LS0txYAB\nAxAVFYX77rsPs2fPBqCu/mwoRrX1Za2qqipERUVhzJgxACzUl6IipaWlEhAQIHl5eVJRUSHR0dFy\n5MgRq8YUEBAgV69eNWqbNWuWvPfeeyIi8t5778kLL7wgIiIZGRkSHR0tlZWVkpeXJwEBAVJWVmby\nmPbv3y9HjhyRsLCwVsVUXl4uIiLh4eGG/h07dqx8/vnnZo8zKSlJ/vKXvyiOtWac+fn5cuLECRER\nKSoqkqCgIDl27Jjq+rShONXYp8XFxSIiUlFRIQ8++KDs3btXdf1ZX4xq7EsRkb/85S8yadIkGTNm\njIhY5vddVSMHtW6Qkztm3nbu3IkpU6YAAOLj4w0x7tixAxMnToSjoyO6deuG0NBQHDp0yOTxDBky\nBF5eXq2OKT09HRcuXEB1dTWioqIUjzFnnICyP60dp1arRVhYGADA1dUVERERuHjxour6tKE4AfX1\nafv27QEA5eXlqKqqgq+vr+r6884YtVotAPX1ZV5eHnbu3IlnnnnGEJsl+lJVyUGNG+Q0Go1h+LZi\nxQoAQEFBAXx8fAAAnTp1wpUrNfsALl68CH9/f8NjLRl/S2O6ePGiUV9369bNYrGuXLkSffv2RXx8\nPK5du6aqOHNycnD48GEMHjxY1X1aG+eQIUMAqK9Pq6ur0a9fP2i1WgwbNgyhoaGq6887YwwJCQGg\nvr6cPXs23nnnHTg4/O/t2hJ9qarkoNFomj7Iwg4ePIgjR45gz5492LBhA77++mtrh2TTnnvuOWRn\nZ+PHH39EYGAgXnjhBWuHZHDr1i2MHz8ey5Ytg7u7u7XDadCtW7cQGxuLZcuWwc3NTZV96uDggGPH\njiEvLw/79+/HN998Y+2QFO6MMS0tTXV9+eWXX8LX1xdRUVEm29zWXKpKDv7+/sjNzTXczs3NNcp2\n1uDr6wsA6Ny5M8aPH4/Dhw+jc+fOKCys2U1cUFBgOObO+O8cCZlTS2Oqr73uXxzm0qlTJ2g0Gmg0\nGvz2t7/F4cOHVRFnRUUFxo0bh8mTJ+Pxxx8HoM4+rY1z0qRJhjjV2qcA4OHhgdGjRyM9PV2V/Vk3\nxoMHD6quL7/77jv885//RK9evRAXF4e9e/diypQpFulLVSWHAQMG4OTJk7h48SIqKiqQkpKCmJgY\nq8VTXFyM4uKaS2Lcvn0bqampCA0NhV6vR3JyzaUIkpOTodfrAQB6vR5bt25FZWUl8vLycPLkScMK\nJ3NraUzdu3eHg4OD4Yq4mzZtMjzGnGqHvwDw2WefITQ01OpxighmzJiBkJAQw6qV2pjU1KcNxam2\nPr169SqKimouAVNSUoLdu3cjPDxcVf3ZUIwFBQWGY9TQl4sXL0Zubi7OnTuHLVu24JFHHsHf/vY3\ny/Slqavqd2vnzp0SGhoqffv2lcWLF1s1lrNnz0pERIRERkZKUFCQvPHGGyIicvXqVRk+fLiEh4fL\niBEj5Pr164bHvPnmm9K3b18JDQ2V1NRUs8Q1ceJE8fPzk3bt2om/v7+sX7++VTFlZGRIv379JCQk\nRJ5//nmzx7lu3TqJj4+XiIgICQ4OlpEjR0peXp7V4/z2229Fo9FIZGSk9OvXT/r16ye7du1SXZ/W\nF+fOnTtV16fHjx+Xfv36SWRkpPTp00cWLlwoIq37vTFXnA3FqLa+rCstLc2wWskSfclNcEREpKCq\naSUiIlIHJgciIlJgciAiIgUmByIiUmByICIiBSYHIiJSYHIgaiadTofMzMx675swYQKys7PN8ry/\n+tWvDBu2iCyFyYGomWovq3CnM2fO4Pbt2wgMDFTcV11dfdfPO3HiRHz88cd3fR6ilmByIJt18+ZN\n6PV6REZGIjw8HCkpKQBqPqDpD3/4A6KjoxEZGYmsrCwAQH5+Ph599FFERkaiX79+2LdvH4CaC9nF\nxcUhMjISoaGh+PTTTwHUXD5l7NixCA0Nxfjx41FSUlLvxc+2bNmCxx57zHDb1dUVc+fORXR0NA4e\nPNhgPE8//TR+//vfY/DgwQgMDERaWhqmTZuG4OBgTJo0yXC+xx57DFu2bDFPJxI1xMQ7vIksZuvW\nrTJz5kzD7aKiIhGp+YCmpUuXiojIpk2b5Ne//rWIiDzxxBNy4MABERE5f/68BAYGiojI7NmzJTk5\nWURErl+/LoGBgXLz5k1ZvHixJCQkiIjIqVOnxMnJSTIzMxVxjBo1yqhdo9EYfZBKQ/FMnTpVJk+e\nLCIi//jHP8TNzU1Onz4t1dXV0r9/fzl8+LDhHL169ZJbt261uq+IWoojB7JZUVFR+Ne//oXXXnsN\n+/fvh6urq+G+p556CgAQGxuL77//HgDw9ddfY9asWYiKisLYsWNRVlaGmzdv4quvvsJbb72FqKgo\nDBs2DJWVlbhw4QIOHDiAuLg4AEBISAgiIiLqjeP8+fPw8/Mz3HZ0dDRcMbWxeDQaDUaPHg0ACAsL\nQ5cuXRAcHAyNRoPQ0FCjq2hqtVqj20Tm5mTtAIhaKygoCJmZmdixYwcSExMxbNgwLFiwoMHjNRoN\nDh8+DCcn5X/72ssi33m8NPPSY3WPc3FxafZnkzg7OwOo+WyBe+65x9Du4OBgVK8QEVV+3gm1XRw5\nkM3Kz89Hhw4dMHnyZLz88svIyMgw3Ldt2zbDvwMHDgQADB8+HKtXrzYcc/LkSQDAyJEjsWrVKkX7\n4MGDsXXrVgDA6dOncfz48Xrj6NmzJy5dutRorPXF0xKXL1+2yOdtENXiyIFs1vHjxzF37lw4OTnB\nycnJ8DGuAFBYWIjo6GhUVlYaCtWrV6/GM888g48++ggigoEDB2LNmjX44x//iJkzZyIkJAROTk7o\n3r07duzYgRdffBETJ05EaGgoQkJCEB0dXW8cgwcPRkZGBvr37w+g/k80rC+eO4+983G1t/Pz8+Hj\n44OOHTu2sqeIWo6X7KY2p1evXsjMzIS3t7dFnu/s2bN4/vnnG/zA9ruNZ82aNbh9+7bRB/wQmRun\nlajNsfTc/L333gs3N7cGN8HdbTxbt27Fs88+e1fnIGopjhyIiEiBIwciIlJgciAiIgUmByIiUmBy\nICIiBSYHIiJSYHIgIiKF/w/ondA2j7ZbKwAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 5 } ], "metadata": {} } ] }