diff options
Diffstat (limited to 'Electrical_Machines_II')
| -rwxr-xr-x | Electrical_Machines_II/chapter_1.ipynb | 1297 | ||||
| -rwxr-xr-x | Electrical_Machines_II/chapter_1_1.ipynb | 1297 | ||||
| -rwxr-xr-x | Electrical_Machines_II/chapter_2.ipynb | 2810 | ||||
| -rwxr-xr-x | Electrical_Machines_II/chapter_2_1.ipynb | 2810 | ||||
| -rwxr-xr-x | Electrical_Machines_II/chapter_3.ipynb | 68 | ||||
| -rwxr-xr-x | Electrical_Machines_II/chapter_3_1.ipynb | 68 | ||||
| -rwxr-xr-x | Electrical_Machines_II/chapter_4.ipynb | 183 | ||||
| -rwxr-xr-x | Electrical_Machines_II/chapter_4_1.ipynb | 183 | ||||
| -rwxr-xr-x | Electrical_Machines_II/screenshots/capacitance_3.png | bin | 0 -> 53826 bytes | |||
| -rwxr-xr-x | Electrical_Machines_II/screenshots/capacitance_3_1.png | bin | 0 -> 53826 bytes | |||
| -rwxr-xr-x | Electrical_Machines_II/screenshots/flux_per_pole.png | bin | 0 -> 68471 bytes | |||
| -rwxr-xr-x | Electrical_Machines_II/screenshots/flux_per_pole_1.png | bin | 0 -> 68471 bytes | |||
| -rwxr-xr-x | Electrical_Machines_II/screenshots/rotor_current_2.png | bin | 0 -> 43696 bytes | |||
| -rwxr-xr-x | Electrical_Machines_II/screenshots/rotor_current_2_1.png | bin | 0 -> 43696 bytes |
14 files changed, 8716 insertions, 0 deletions
diff --git a/Electrical_Machines_II/chapter_1.ipynb b/Electrical_Machines_II/chapter_1.ipynb new file mode 100755 index 00000000..0ea2a442 --- /dev/null +++ b/Electrical_Machines_II/chapter_1.ipynb @@ -0,0 +1,1297 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:e4feee2966703522c4e009f16556422b345ba2813bb92ad5cc7b8ab8058895bb" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 1: Synchronous Machines\n" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.1, Page 24" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "slots=24;#no. of slotes\n", + "NoOfPhase=3;#no of phase\n", + "MotorSpeed=1450;#in rpm\n", + "N=1500;#Synchonous Speed in rpm\n", + "f=50;#in Hz\n", + "\n", + "#Calculations&Results\n", + "print \"As the winding is in single layer, each slot contains one coil slide only.\";\n", + "CoilSlidePerSlot=1;#coil slide per slot\n", + "CoilSlidePerCoil=2;#coil slide per Coil\n", + "TotalCoils=slots*CoilSlidePerSlot/CoilSlidePerCoil;#no. of coils\n", + "print \"Total no. of Coils : \",TotalCoils;\n", + "P=120*f/N;#no. of poles\n", + "print \"No. of poles : \",P;\n", + "PolesPitch=slots/P;#unitless\n", + "print \"Poles Pitch = \",PolesPitch;\n", + "print \"In case of single layer winding, the pole ppitch is generally taken in odd numbers only\"\n", + "print \"let pole pitch = 5(for short pitch winding)\";\n", + "PolesPitch=5;#for short pitch winding\n", + "print \"Coil Span or coil through = 1-6\";\n", + "CoilsPerPolePerPhase=TotalCoils/(P*NoOfPhase);#No. of Coils/Pole/Phase\n", + "print \"No. of Coils/Pole/Phase = \",CoilsPerPolePerPhase;\n", + "pair_of_poles=2;#no. of pair of poles\n", + "TotalElectricalDegree=360*pair_of_poles;#in degree \n", + "ElectricalDegreesPerSlot=TotalElectricalDegree/slots;#in degree electrical\n", + "print \"Electrical Degrees/Slot = \",ElectricalDegreesPerSlot;\n", + "Slots_required=120/ElectricalDegreesPerSlot;#No. of slotes required for proper phase displacement\n", + "print \"No. of slotes required for proper phase displacement = \",Slots_required;\n", + "print \"Winding Table is as follows :\";\n", + "print \"Coil No. Connection Lead from Coil Span Phase and Group No.\";\n", + "print \" 1 A1 1-6 A1\";\n", + "print \" 2 3-8 C4\";\n", + "print \" 3 B1 5-10 B1\";\n", + "print \" 4 7-12 A2\";\n", + "print \" 5 C1 9-14 C1\";\n", + "print \" 6 11-16 B2\";\n", + "print \" 7 13-18 A3\";\n", + "print \" 8 15-20 C2\";\n", + "print \" 9 17-22 B3\";\n", + "print \" 10 19-24 A4\";\n", + "print \" 11 21-2 C3\";\n", + "print \" 12 23-4 B4\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "As the winding is in single layer, each slot contains one coil slide only.\n", + "Total no. of Coils : 12\n", + "No. of poles : 4\n", + "Poles Pitch = 6\n", + "In case of single layer winding, the pole ppitch is generally taken in odd numbers only\n", + "let pole pitch = 5(for short pitch winding)\n", + "Coil Span or coil through = 1-6\n", + "No. of Coils/Pole/Phase = 1\n", + "Electrical Degrees/Slot = 30\n", + "No. of slotes required for proper phase displacement = 4\n", + "Winding Table is as follows :\n", + "Coil No. Connection Lead from Coil Span Phase and Group No.\n", + " 1 A1 1-6 A1\n", + " 2 3-8 C4\n", + " 3 B1 5-10 B1\n", + " 4 7-12 A2\n", + " 5 C1 9-14 C1\n", + " 6 11-16 B2\n", + " 7 13-18 A3\n", + " 8 15-20 C2\n", + " 9 17-22 B3\n", + " 10 19-24 A4\n", + " 11 21-2 C3\n", + " 12 23-4 B4\n" + ] + } + ], + "prompt_number": 52 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.2, Page 28" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "slots=24;#no. of slotes\n", + "P=4;#no. of poles\n", + "CoilPitch=5;#unitless\n", + "NoOfPhase=3;#no of phase\n", + "\n", + "#Calculations&Results\n", + "print \"As the winding is in double layer, each slot contains 2 coils per slide.\";\n", + "CoilSlidePerSlot=2;#coil slide per slot\n", + "CoilSlidePerCoil=2;#coil slide per Coil\n", + "TotalCoils=slots*CoilSlidePerSlot/CoilSlidePerCoil;#no. of coils\n", + "print \"Total no. of Coils : \",TotalCoils;\n", + "PolesPitch=slots/P;#unitless\n", + "print \"Pole Pitch = \",PolesPitch;\n", + "CoilsPerPolePerPhase=TotalCoils/(P*NoOfPhase);#No. of Coils/Pole/Phase\n", + "print \"No. of Coils/Pole/Phase = \",CoilsPerPolePerPhase;\n", + "pair_of_poles=2;#no. of pair of poles\n", + "TotalElectricalDegree=360*pair_of_poles;#in degree \n", + "ElectricalDegreesPerSlot=TotalElectricalDegree/slots;#in degree electrical\n", + "print \"Electrical Degrees/Slot = \",ElectricalDegreesPerSlot;\n", + "Slots_required=120/ElectricalDegreesPerSlot;#No. of slotes required for proper phase displacement\n", + "print \"No. of slots required for proper phase displacement = \",Slots_required;\n", + "print \"ie. Phase A1 is brought out from slot no. = 1\";\n", + "print \"Phase B1 at slot no. = 1+4 = 5\";\n", + "print \"Phase C1 at slot no. = 5+4 = 9\";\n", + "print \"Col Connection - end to start to start\";\n", + "print \"Winding Table is as follows :\";\n", + "print \"Coil No. Connection Lead from Coil Span Phase and Group No.\";\n", + "print \" 1 A1 1-6 A1\";\n", + "print \" 2 2-7 \";\n", + "print \" 3 3-8 C4\";\n", + "print \" 4 4-9 \";\n", + "print \" 5 B1 5-10 B1\";\n", + "print \" 6 6-11 \";\n", + "print \" 7 7-12 A2\";\n", + "print \" 8 8-13 \";\n", + "print \" 9 C1 9-14 C1\";\n", + "print \" 10 10-15 \";\n", + "print \" 11 11-16 B2\";\n", + "print \" 12 12-17 \";\n", + "print \" 13 13-18 A3\";\n", + "print \" 14 14-19 \";\n", + "print \" 15 15-20 C2\";\n", + "print \" 16 16-21 \";\n", + "print \" 17 17-22 B3\";\n", + "print \" 18 18-23 \";\n", + "print \" 19 19-24 A4\";\n", + "print \" 20 20-1 \";\n", + "print \" 21 21-2 C3\";\n", + "print \" 22 22-3 \";\n", + "print \" 23 23-4 \";\n", + "print \" 24 24-5 B4\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "As the winding is in double layer, each slot contains 2 coils per slide.\n", + "Total no. of Coils : 24\n", + "Pole Pitch = 6\n", + "No. of Coils/Pole/Phase = 2\n", + "Electrical Degrees/Slot = 30\n", + "No. of slots required for proper phase displacement = 4\n", + "ie. Phase A1 is brought out from slot no. = 1\n", + "Phase B1 at slot no. = 1+4 = 5\n", + "Phase C1 at slot no. = 5+4 = 9\n", + "Col Connection - end to start to start\n", + "Winding Table is as follows :\n", + "Coil No. Connection Lead from Coil Span Phase and Group No.\n", + " 1 A1 1-6 A1\n", + " 2 2-7 \n", + " 3 3-8 C4\n", + " 4 4-9 \n", + " 5 B1 5-10 B1\n", + " 6 6-11 \n", + " 7 7-12 A2\n", + " 8 8-13 \n", + " 9 C1 9-14 C1\n", + " 10 10-15 \n", + " 11 11-16 B2\n", + " 12 12-17 \n", + " 13 13-18 A3\n", + " 14 14-19 \n", + " 15 15-20 C2\n", + " 16 16-21 \n", + " 17 17-22 B3\n", + " 18 18-23 \n", + " 19 19-24 A4\n", + " 20 20-1 \n", + " 21 21-2 C3\n", + " 22 22-3 \n", + " 23 23-4 \n", + " 24 24-5 B4\n" + ] + } + ], + "prompt_number": 53 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.3, Page 35" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import *\n", + "\n", + "#given data\n", + "StatorSlots=36;#No. of stator slots\n", + "Poles=4;#No. of poles\n", + "#coilSpan=1:8;#unitless\n", + "\n", + "#Calculations&Results\n", + "SlotsPerPole=StatorSlots/Poles;# no. of slots per pole\n", + "print \"Slots/Pole = \",SlotsPerPole;\n", + "print \"Here the coil span falls short by, \",((2./9)*180),\" Degree\";\n", + "alfa=40*pi/180;#short pitch angle in degree\n", + "Kp=cos(alfa/2);#Coil span Factor\n", + "print \"Pitch Factor or coil span factor : \",round(Kp,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slots/Pole = 9\n", + "Here the coil span falls short by, 40.0 Degree\n", + "Pitch Factor or coil span factor : 0.9\n" + ] + } + ], + "prompt_number": 62 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.4, Page 38" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import *\n", + "\n", + "#given data\n", + "NoOfPhase=3;#no of phase\n", + "P=16;#No. of pole alternator\n", + "Slots=144;#No. of slots\n", + "Conductors=10;#per slot\n", + "fi=0.03;#in Weber\n", + "N=375;#machine speed in rpm\n", + "\n", + "#Calculations\n", + "f=P*N/120;#in Hz\n", + "SlotsPerPole=Slots/P;#unitless\n", + "m=Slots/(P*NoOfPhase);#unitless\n", + "Beta=180/(SlotsPerPole);#in Degree\n", + "Kd=sin(m*Beta/2*pi/180)/(m*sin(Beta/2*pi/180));#unitless\n", + "TotalConductors=Conductors*Slots;#no. of conductors\n", + "TotalConductorsPerPhase=Conductors*Slots/NoOfPhase;#no. of conductors/phase\n", + "TurnsPerPhase=TotalConductorsPerPhase/2;#No. of turns per phase\n", + "EMFPerPhase=4.44*Kd*fi*f*TurnsPerPhase;#in Volt\n", + "LineVoltage=sqrt(3)*EMFPerPhase;#in Volt\n", + "\n", + "#Results\n", + "print \"Frequency in Hz : \",f;\n", + "print \"Phase Electromotive force in Volt : \",round(EMFPerPhase,1);\n", + "print \"Line Electromotive force in Volt : \",round(LineVoltage,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Frequency in Hz : 50\n", + "Phase Electromotive force in Volt : 1534.1\n", + "Line Electromotive force in Volt : 2657.202\n" + ] + } + ], + "prompt_number": 76 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.5, Page 39" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "NoOfPhase=3;#no of phase\n", + "SlotsPerPhase=3;#o. of slots\n", + "\n", + "#Calculations&Results\n", + "m=SlotsPerPhase;#no. of slots\n", + "SlotsPerPolePerPhase=SlotsPerPhase*NoOfPhase;#unitless\n", + "Beta=180./SlotsPerPolePerPhase;#in degree\n", + "print \"The phase difference between the induced emf in two coils lying in adjacent slots = \",(Beta),\" degree\";\n", + "Kb=sin((m*Beta*pi)/(2*180))/(m*sin((Beta*pi)/(2*180)));#unitless\n", + "print \"The breadth factor : \",round(Kb,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The phase difference between the induced emf in two coils lying in adjacent slots = 20.0 degree\n", + "The breadth factor : 0.96\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.6, Page 39" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "StatorSlots=24#No. of stator slots\n", + "Poles=4;#No. of poles\n", + "SlotsPerPole=StatorSlots/Poles;# no. of slots per pole\n", + "#coilSpan=1:6;#unitless\n", + "\n", + "#Calculations&Results\n", + "print \"If the sides of the coil are placed in slots 1 and 7, then it is full pitched, \\nIf the coil slides are placed\" \\\n", + " + \" in 1 and 6 then it is short pitched and the distance equal to 5/6th of pole-pitch.\";\n", + "print \"Since it falls short by 1/6th of the pole-pitch, hence it is short by :\"\n", + "theta=180/SlotsPerPole;#in Degree\n", + "print (theta),\" Degree.\"\n", + "Kp=cos(theta/2*pi/180);#unitless\n", + "print \"Pitch-factor Kp : \",round(Kp,4);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "If the sides of the coil are placed in slots 1 and 7, then it is full pitched, \n", + "If the coil slides are placed in 1 and 6 then it is short pitched and the distance equal to 5/6th of pole-pitch.\n", + "Since it falls short by 1/6th of the pole-pitch, hence it is short by :\n", + "30 Degree.\n", + "Pitch-factor Kp : 0.9659\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.7, Page 40" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "NoOfPhase=3;#no of phase\n", + "Eph=3300/sqrt(3);#in Volts\n", + "f=50;#in Hz\n", + "Poles=12;#No. of poles\n", + "StatorSlots=144#No. of stator slots\n", + "\n", + "#Calculations\n", + "SlotsPerPhase=StatorSlots/NoOfPhase;#no. of slots/phase\n", + "Conductors=5;#per slot\n", + "ConductorsPerphase=SlotsPerPhase*Conductors;#Conductors/Phase\n", + "S=ConductorsPerphase;#Conductors/phase\n", + "SlotsPerPolePerPhase=SlotsPerPhase/Poles;#no. of slots/phase\n", + "Kf=1.11;#Form Factor\n", + "Kb=0.96;#Breadth Factor\n", + "Kp=1;#For concentric winding\n", + "fi=Eph/(2*Kf*Kb*Kp*S*f);#in weber\n", + "\n", + "#Result\n", + "print \"The Flux per pole in weber : \",round(fi,4);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The Flux per pole in weber : 0.0745\n" + ] + } + ], + "prompt_number": 93 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.8, Page 41" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "Phase=3;#no. of phase\n", + "f=50;#in Hz\n", + "P=16;#No. of pole alternator\n", + "Slots=144.;#No. of slots\n", + "conductors=10;#conductors per slot\n", + "\n", + "#Calculations&Results\n", + "fi=2.48*10**-2;#in weber\n", + "n=Slots/P;#No. of slots/pole\n", + "Zr=Slots*conductors/Phase;#No. of conductors/Phase\n", + "T=Zr/2;#N. of turns/phase\n", + "Beta=180/n;#Angular displacement between slots in degree\n", + "m=n/Phase;#No. of slots/pole/Phase\n", + "Kd=sin(m*Beta/2*pi/180)/(m*sin(Beta/2*pi/180));#Distribution factor :unitless \n", + "print \"The coil span falls short of 2 slots i.e. \",(2*180/9),\" degree.\";\n", + "alfa=40;#short pitch angle in degree\n", + "Kp=(cos(alfa/2*pi/180));#Unitless\n", + "#Formula : f=P*N/120;#in Hz\n", + "N=120*f/P;#in rpm\n", + "Ep=4.44*Kd*Kp*fi*f*T;#in Volts\n", + "LineVoltage=sqrt(3)*Ep;#in Volts\n", + "print \"The speed is \",(N),\" in rpm\";\n", + "print \"The line emf is \",round(LineVoltage,2),\" Volts\";\n", + "#Note : Answer in the book is not accurate for last part due to rounding off errors" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The coil span falls short of 2 slots i.e. 40 degree.\n", + "The speed is 375 in rpm\n", + "The line emf is 2064.15 Volts\n" + ] + } + ], + "prompt_number": 109 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.9, Page 58" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "RatedPower=100;#in KVA\n", + "RatedPower=100*1000;#in VA\n", + "VL=1040;#in Volt\n", + "Phase=3;#Machine phase\n", + "If=40;#in Ampere\n", + "Isc=200;#in Ampere\n", + "EL=1040;#in Volt\n", + "\n", + "#Calculations&Results\n", + "Eph=EL/sqrt(3);#in Volt\n", + "Zs=Eph/Isc;#in Ohm\n", + "Rs=0.2;#in Ohm\n", + "Xs=sqrt(Zs**2-Rs**2);#in Ohm\n", + "IL=19.25;#in Ampere\n", + "V=3000/sqrt(3);#in Volt\n", + "#At 0.8 power factor lagging\n", + "IRa=IL*0.2;#in Volt\n", + "IXs=IL*Xs;#in Volt\n", + "Vsin_fi=V*0.6;#in Volt\n", + "Vcos_fi=V*0.8;#in Volt\n", + "Eo=sqrt((Vcos_fi+IRa)**2+(Vsin_fi+IXs)**2);#in Volts\n", + "Regulation=((Eo-V)/V)*100;#in %\n", + "print \"Full load percentage regulation at a power factor of 0.8 lagging : \",round(Regulation,1);\n", + "\n", + "#At 0.8 power factor leading\n", + "Eo=sqrt((Vcos_fi+IRa)**2+(Vsin_fi-IXs)**2);#in Volts\n", + "Regulation=((Eo-V)/V)*100;#in %\n", + "print \"Full load percentage regulation at a power factor of 0.8 leading : \",round(Regulation,2);\n", + "print \"Negative regulation due to leading power factor.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Full load percentage regulation at a power factor of 0.8 lagging : 2.2\n", + "Full load percentage regulation at a power factor of 0.8 leading : -1.78\n", + "Negative regulation due to leading power factor.\n" + ] + } + ], + "prompt_number": 110 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.10, Page 60" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "RatedPower=50;#in KVA\n", + "RatedPower=50*1000;#in VA\n", + "VL=173;#in Volts\n", + "Ra=0.1;#in Ohm\n", + "\n", + "#Calculations&Results\n", + "VP=VL/sqrt(3);#in Volts\n", + "print \"Some exciting curent on short circuit produces a current of 100 A.\";\n", + "OC_PhaseVoltage=100;#in Volt\n", + "SC_Current=100;#in Ampere\n", + "Zs=OC_PhaseVoltage/SC_Current;#n ohm\n", + "Xs=sqrt(Zs**2-Ra**2);#in Ohm\n", + "print \"Impedence of the alternator in Ohm : \",round(Xs,2);\n", + "V=400;#in Volts\n", + "I_FL=RatedPower/(sqrt(3)*V);#in Ampere\n", + "V=400/sqrt(3);#in Volts\n", + "Eo=sqrt((V+I_FL*Ra)**2+(I_FL*Xs)**2);#in Volts\n", + "Regulation=(Eo-V)*100/V;#in %\n", + "print \"Regulation at U.P.F. in % :\",round(Regulation,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Some exciting curent on short circuit produces a current of 100 A.\n", + "Impedence of the alternator in Ohm : 0.99\n", + "Regulation at U.P.F. in % : 7.71\n" + ] + } + ], + "prompt_number": 112 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.11, Page 61" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=500;#in KVA\n", + "OutputPower=500*1000;#in VA\n", + "VL=3300;#in Volts\n", + "Ra=0.3;#in Ohm\n", + "Xs=4;#in Ohm\n", + "PF=0.8;#Lagging Power factor\n", + "\n", + "#Calculations&Results\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "IL=OutputPower/(sqrt(3)*VL);#in Ampere\n", + "print \"For a star connected alternator, line current is equal to phase current. Therefore Ia=IL\";\n", + "Ia=IL;#in Ampere\n", + "#PF=cosd(fi)=0.8 and sind(fi)=0.6\n", + "cos_fi=0.8;#Power factor\n", + "sin_fi=0.6;#Unitless\n", + "VPerPhase=VL/sqrt(3);#in Volts\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);#in Volts/Phase\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "print \"Voltage Regulation at Full Load in % :\",round(Regulation,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For a star connected alternator, line current is equal to phase current. Therefore Ia=IL\n", + "Voltage Regulation at Full Load in % : 12.98\n" + ] + } + ], + "prompt_number": 113 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.12, Page 62" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "V=2000;#in Volt\n", + "Ia=100;#in Ampere\n", + "OC_Voltage=500;#in Volt\n", + "SC_Current=100;#in Ampere\n", + "\n", + "#Calculations&Results\n", + "Zs=OC_Voltage/SC_Current;#in ohm\n", + "Ra=0.8;#in Ohm\n", + "Xs=sqrt(Zs**2-Ra**2);#in Ohm\n", + "#formula : Induced EMF, E=sqrt((V*cos_fi+Ia*Ra)^2+(V*sin_fi+Ia*Xs)^2)\n", + "#Part (a) : at unity pf\n", + "cos_fi=1;#Unitless\n", + "sin_fi=0;#Unitless\n", + "E=sqrt((V*cos_fi+Ia*Ra)**2+(V*sin_fi+Ia*Xs)**2)\n", + "Regulation=(E-V)*100/V;#in %\n", + "print \"Regulation at U.P.F. in % :\",round(Regulation,2);\n", + "\n", + "#Part (b) : at 0.71 pf lagging\n", + "cos_fi=0.71;#Unitless\n", + "sin_fi=0.704;#Unitless\n", + "E=sqrt((V*cos_fi+Ia*Ra)**2+(V*sin_fi+Ia*Xs)**2)\n", + "Regulation=(E-V)*100/V;#in %\n", + "print \"Regulation at 0.71 pf lagging in % :\",round(Regulation,2);\n", + "\n", + "#Part (c) : at 0.8 pf leading\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((V*cos_fi+Ia*Ra)**2+(V*sin_fi-Ia*Xs)**2)\n", + "Regulation=(E-V)*100/V;#in %\n", + "print \"Regulation at 0.8 pf leading in % :\",round(Regulation,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Regulation at U.P.F. in % : 6.89\n", + "Regulation at 0.71 pf lagging in % : 21.1\n", + "Regulation at 0.8 pf leading in % : -8.88\n" + ] + } + ], + "prompt_number": 114 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.13, Page 63" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "Ia=60;#in Ampere\n", + "\n", + "#Calculations&Results\n", + "print \"The value of synchronous impedence at this excitation :\" ;\n", + "OC_Voltage=900;#in Volt\n", + "SC_Current=150;#in Ampere\n", + "Zs=OC_Voltage/SC_Current;#in ohm\n", + "print \"Zs equals to \",(Zs),\" Ohm\";\n", + "print \"Internal Voltage drop when tthe load current is 60A=Ia*Zs=\",(Ia*Zs),\" Volts\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The value of synchronous impedence at this excitation :\n", + "Zs equals to 6 Ohm\n", + "Internal Voltage drop when tthe load current is 60A=Ia*Zs= 360 Volts\n" + ] + } + ], + "prompt_number": 118 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.14, Page 63" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "V=6600;#in Volts\n", + "OutputPower=2000;#in KVA\n", + "\n", + "#Calculations\n", + "OutputPower=2000*1000;#in VA\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "IL=OutputPower/(sqrt(3)*V);#in Ampere\n", + "Ia=IL;#in Ampere\n", + "Ra=0.4;#in Ohm\n", + "Xs=4.5;#in Ohm\n", + "#PF=cosd(fi)=0.8 and sind(fi)=0.6\n", + "cos_fi=0.8;#Power factor\n", + "sin_fi=0.6;#Unitless\n", + "VPerPhase=V/sqrt(3);#in Volts\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);#in Volts/Phase\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "\n", + "#Result\n", + "print \"Percentage Change in Terminal Voltage :\",round(Regulation,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Percentage Change in Terminal Voltage : 14.9\n" + ] + } + ], + "prompt_number": 119 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.15, Page 63" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=1200;#in KVA\n", + "OutputPower=1200*1000;#in VA\n", + "V=3300;#in Volt\n", + "Ra=0.25;#in Ohm\n", + "\n", + "#Calculations&Results\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "IL=OutputPower/(sqrt(3)*V);#in Ampere\n", + "Ia=IL;#in Ampere\n", + "VPerPhase=V/sqrt(3);#in Volts\n", + "OC_Voltage=1100;#in Volt\n", + "SC_Current=200;#in Ampere\n", + "Zs=OC_Voltage/(sqrt(3)*SC_Current);#in ohmRa\n", + "Xs=sqrt(Zs**2-Ra**2);#in Ohm\n", + "#formula : Induced EMF, E=sqrt((V*cos_fi+Ia*Ra)^2+(V*sin_fi+Ia*Xs)^2)\n", + "\n", + "#Part (a) : For lagging pf load\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "print \"Regulation at U.P.F. in % :\",round(Regulation,2);\n", + "\n", + "#Part (b) : For leading pf load\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi-Ia*Xs)**2)\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "print \"Regulation at 0.71 pf lagging in % :\",round(Regulation,2);\n", + "#Note: For leading power factor load, the regulation s negative. " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Regulation at U.P.F. in % : 25.9\n", + "Regulation at 0.71 pf lagging in % : -13.52\n" + ] + } + ], + "prompt_number": 120 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.16, Page 64" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=1500;#in KVA\n", + "OutputPower=1500*1000;#in VA\n", + "V=6600;#in Volt\n", + "Ra=0.4;#in Ohm\n", + "Xs=6;#in Ohm per phase\n", + "pf=0.8;#lagging power factor\n", + "\n", + "#Calculations&Results\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "Ia=OutputPower/(sqrt(3)*V);#in Ampere\n", + "VPerPhase=V/sqrt(3);#in Volts\n", + "#formula : Induced EMF, E=sqrt((V*cos_fi+Ia*Ra)^2+(V*sin_fi+Ia*Xs)^2)\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);#in volt\n", + "print \"Induced emf in volt : \",round(E,2);\n", + "print \"As excitation remains constant, E at 4364 volt remains constant.\";\n", + "E=4364;#in Volt\n", + "V=4743;#in Volts\n", + "TerminalVoltage=sqrt(3)*V;#in Volts\n", + "print \"Terminal voltage line to line in Volts : \",round(TerminalVoltage);\n", + "#Note ans of 1st part is wrong in the books" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Induced emf in volt : 4366.07\n", + "As excitation remains constant, E at 4364 volt remains constant.\n", + "Terminal voltage line to line in Volts : 8215.0\n" + ] + } + ], + "prompt_number": 123 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.17, Page 81" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=2500;#in KVA\n", + "OutputPower=2500*1000;#in VA\n", + "V=6600;#in Volt\n", + " \n", + "#Calculations\n", + "#For first load i.e. Lighting load : \n", + "KW1=2500;#in KWatts\n", + "KVAR1=0;#Kwatts\n", + "#For second load i.e. Motor load : \n", + "KW2=5000;#in KWatts\n", + "cos_fi=0.707;#unitless\n", + "sin_fi=0.707;#unitless\n", + "KVAR2=KW2*sin_fi/cos_fi;#Kwatts\n", + "#For total load\n", + "TotalKW=KW1+KW2;#in KWatts\n", + "TotalKVAR=KVAR1+KVAR2;#Kwatts \n", + "#For first Machine\n", + "KWm=4000;#in KWatts\n", + "cos_fi=0.8;#unitless\n", + "sin_fi=0.6;#unitless\n", + "KVARm=KWm*sin_fi/cos_fi;#Kwatts\n", + "#so, second machine will be supplying\n", + "KW=TotalKW-KWm;#in Kwatts\n", + "print \"KW output of second machine : \",KW;\n", + "KVAR=TotalKVAR-KVARm;#in KWatts\n", + "tan_fi=KVAR/KW;#unitless\n", + "fi=atan(tan_fi);#ib degree\n", + "#Power factor of other machine\n", + "pf=cos(fi);#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of other machine : \",round(pf,4),\" lagging\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "KW output of second machine : 3500\n", + "Power factor of other machine : 0.8682 lagging\n" + ] + } + ], + "prompt_number": 126 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.18, Page 82" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "\n", + "#Load1 : \n", + "KW1=500;#in KWatts\n", + "KVAR1=0;#Kwatts\n", + "\n", + "#Load2 : \n", + "KW2=1000;#in KWatts\n", + "pf=0.9;#lagging\n", + "cos_fi=0.9;#unitless\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR2=KW2*sin_fi/cos_fi;#Kwatts\n", + "\n", + "#Load3 : \n", + "KW3=800;#in KWatts\n", + "pf=0.8;#lagging\n", + "cos_fi=0.8;#unitless\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR3=KW3*sin_fi/cos_fi;#Kwatts\n", + "\n", + "#Load4 : \n", + "KW4=500;#in KWatts\n", + "pf=0.9;#lagging\n", + "cos_fi=0.9;#unitless\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR4=-KW4*sin_fi/cos_fi;#Kwatts\n", + "#TOtalKW and TotalKVAR\n", + "TotalKW=KW1+KW2+KW3+KW4;#in KWatts\n", + "TotalKVAR=KVAR1+KVAR2+KVAR3+KVAR4;#in KWAtts\n", + "#For the first Machine :\n", + "KW=1500;#n Kwatts\n", + "cos_fi=0.95;#unitless\n", + "sin_fi=0.3123;#unitless\n", + "KVAR=KW*sin_fi/cos_fi;#Kwatts\n", + "\n", + "KW1=TotalKW-KW;#in KWatts\n", + "KVAR1=TotalKVAR-KVAR;#in Volts\n", + "print \"KW supplied by other machine : \",(TotalKW-KW);\n", + "print \"KVAR supplied by other machine : \",round(TotalKVAR-KVAR); #answer differs due to rounding off the digits\n", + "tan_fi=KVAR1/KW1;#unitless\n", + "#fi=atand(tan_fi);#in degree\n", + "cos_fi=cos(atan(tan_fi));#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of the other machine : \",round(cos_fi,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "KW supplied by other machine : 1300\n", + "KVAR supplied by other machine : 349.0\n", + "Power factor of the other machine : 0.966\n" + ] + } + ], + "prompt_number": 129 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.19, Page 82" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "\n", + "#Lighting Load : \n", + "MW1=20;#load in Mwatts\n", + "KW1=MW1*1000;#in KWatts\n", + "KVAR1=0;#Kwatts\n", + "\n", + "#Motor Load : \n", + "MW2=40;#load in Mwatts\n", + "KW2=MW2*1000;#in KWatts\n", + "pf=0.8;#unitless\n", + "cos_fi=0.8\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR2=KW2*sin_fi/cos_fi;#Kwatts\n", + "\n", + "#For Total Load : \n", + "TotalKW=KW1+KW2;#load in Mwatts\n", + "TotalKVAR=KVAR1+KVAR2;#in KWatts\n", + "#For first machine : \n", + "MWm=32;#load in Mwatts\n", + "KWm=MWm*1000;#in KWatts\n", + "cos_fi=0.866;#unitless\n", + "fi=acos(cos_fi);\n", + "tan_fi=tan(fi);#unitless\n", + "KVARm=KWm*tan_fi;#in KWatts\n", + "#so, load supplied by the second machine\n", + "KW2=TotalKW-KWm;#in Kwatts\n", + "print \"Load of other machine,KW : \",(KW2);\n", + "KVAR2=TotalKVAR-KVARm;#in Kwatts\n", + "tan_fi=KVAR2/KW2;#unitless\n", + "fi=atan(tan_fi);#in degree\n", + "cos_fi=cos(atan(tan_fi));#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of the other machine : \",round(cos_fi,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Load of other machine,KW : 28000\n", + "Power factor of the other machine : 0.925\n" + ] + } + ], + "prompt_number": 131 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.20, Page 84" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "cos_fi=0.8;#unitless\n", + "fi=acos(cos_fi);\n", + "tan_fi=tan(fi);#unitless\n", + "\n", + "#For Alternator A : \n", + "cos_fi_A=0.9;#unitless\n", + "fi_A=acos(cos_fi_A);\n", + "tan_fi_A=tan(fi_A);#unitless\n", + "#Formula : Active load, KW=V*I*cos_fi\n", + "#Formula : Reactive load, KVAR=V*I*sin_fi\n", + "ActiveLoad=8000;#in KW\n", + "ReactiveLoad=ActiveLoad*tan_fi;#in KVAR\n", + "\n", + "#For A:\n", + "ActiveLoadA=5000;#in KW\n", + "ReactiveLoadA=ActiveLoadA*tan_fi_A;#in KVAR\n", + "\n", + "#For B :\n", + "ActiveLoadB=ActiveLoad-ActiveLoadA;#in KW\n", + "ReactiveLoadB=ReactiveLoad-ReactiveLoadA;#in KVAR\n", + "tan_fi_B=ReactiveLoadB/ActiveLoadB;#unitless\n", + "fi_B=atan(tan_fi_B);#in degree\n", + "cos_fi=cos(atan(tan_fi_B));#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of the other machine : \",round(cos_fi,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Power factor of the other machine : 0.642\n" + ] + } + ], + "prompt_number": 133 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.21, Page 97" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "V=6600;#in Volts\n", + "KW=6000.;#in KWatts\n", + "pf=0.8;#unitless\n", + "cos_fi=pf;#unitless\n", + "Eff=90.;#in %\n", + "\n", + "#Calculations&Results\n", + "#Part (a) : \n", + "KVA=KW/cos_fi;#in KVAR\n", + "print \"KVA rating of the alternator : \",(KVA),\" KVA\";\n", + "#Part (b) : \n", + "TotalRating=KVA;#in KVA\n", + "VA=TotalRating*1000;#in VA\n", + "I=VA/(sqrt(3)*KW);#in Ampere\n", + "print \"Current Rating in Ampere : \",round(I,2);\n", + "#Part (c) :\n", + "Input=KW/(Eff/100);#in KW\n", + "print \"Power Input(in KW) :\",round(Input,2);\n", + "Input=Input*1000/735.5;#in hp\n", + "print \"Power Input(in hp) :\",round(Input,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "KVA rating of the alternator : 7500.0 KVA\n", + "Current Rating in Ampere : 721.69\n", + "Power Input(in KW) : 6666.67\n", + "Power Input(in hp) : 9064.13\n" + ] + } + ], + "prompt_number": 139 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.22, Page 97" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "Ecoil=8000;#in Volts\n", + "Icoil=418;#in Ampere\n", + "pf=80.;#in % lgging\n", + "\n", + "#Calculations&Results\n", + "pf=pf/100;#in fraction\n", + "cos_fi=pf;#unitless\n", + "#Part (i) : \n", + "EL=sqrt(3)*Ecoil;#in volt\n", + "print \"Line volts(in V): \",round(EL);\n", + "#Part (ii) : \n", + "IL=Icoil;#in Ampere\n", + "print \"Line Current in Ampere : \",IL;\n", + "#Part (iii) :\n", + "Rating=sqrt(3)*EL*IL/1000;#in KVA\n", + "print \"Rating (in KVA) :\",Rating;\n", + "#Part (iv) :\n", + "FullLoadPower=sqrt(3)*EL*IL*cos_fi/1000;#in KW\n", + "print \"Full Load Power in KW :\",FullLoadPower;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Line volts(in V): 13856.0\n", + "Line Current in Ampere : 418\n", + "Rating (in KVA) : 10032.0\n", + "Full Load Power in KW : 8025.6\n" + ] + } + ], + "prompt_number": 140 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/chapter_1_1.ipynb b/Electrical_Machines_II/chapter_1_1.ipynb new file mode 100755 index 00000000..0ea2a442 --- /dev/null +++ b/Electrical_Machines_II/chapter_1_1.ipynb @@ -0,0 +1,1297 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:e4feee2966703522c4e009f16556422b345ba2813bb92ad5cc7b8ab8058895bb" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 1: Synchronous Machines\n" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.1, Page 24" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "slots=24;#no. of slotes\n", + "NoOfPhase=3;#no of phase\n", + "MotorSpeed=1450;#in rpm\n", + "N=1500;#Synchonous Speed in rpm\n", + "f=50;#in Hz\n", + "\n", + "#Calculations&Results\n", + "print \"As the winding is in single layer, each slot contains one coil slide only.\";\n", + "CoilSlidePerSlot=1;#coil slide per slot\n", + "CoilSlidePerCoil=2;#coil slide per Coil\n", + "TotalCoils=slots*CoilSlidePerSlot/CoilSlidePerCoil;#no. of coils\n", + "print \"Total no. of Coils : \",TotalCoils;\n", + "P=120*f/N;#no. of poles\n", + "print \"No. of poles : \",P;\n", + "PolesPitch=slots/P;#unitless\n", + "print \"Poles Pitch = \",PolesPitch;\n", + "print \"In case of single layer winding, the pole ppitch is generally taken in odd numbers only\"\n", + "print \"let pole pitch = 5(for short pitch winding)\";\n", + "PolesPitch=5;#for short pitch winding\n", + "print \"Coil Span or coil through = 1-6\";\n", + "CoilsPerPolePerPhase=TotalCoils/(P*NoOfPhase);#No. of Coils/Pole/Phase\n", + "print \"No. of Coils/Pole/Phase = \",CoilsPerPolePerPhase;\n", + "pair_of_poles=2;#no. of pair of poles\n", + "TotalElectricalDegree=360*pair_of_poles;#in degree \n", + "ElectricalDegreesPerSlot=TotalElectricalDegree/slots;#in degree electrical\n", + "print \"Electrical Degrees/Slot = \",ElectricalDegreesPerSlot;\n", + "Slots_required=120/ElectricalDegreesPerSlot;#No. of slotes required for proper phase displacement\n", + "print \"No. of slotes required for proper phase displacement = \",Slots_required;\n", + "print \"Winding Table is as follows :\";\n", + "print \"Coil No. Connection Lead from Coil Span Phase and Group No.\";\n", + "print \" 1 A1 1-6 A1\";\n", + "print \" 2 3-8 C4\";\n", + "print \" 3 B1 5-10 B1\";\n", + "print \" 4 7-12 A2\";\n", + "print \" 5 C1 9-14 C1\";\n", + "print \" 6 11-16 B2\";\n", + "print \" 7 13-18 A3\";\n", + "print \" 8 15-20 C2\";\n", + "print \" 9 17-22 B3\";\n", + "print \" 10 19-24 A4\";\n", + "print \" 11 21-2 C3\";\n", + "print \" 12 23-4 B4\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "As the winding is in single layer, each slot contains one coil slide only.\n", + "Total no. of Coils : 12\n", + "No. of poles : 4\n", + "Poles Pitch = 6\n", + "In case of single layer winding, the pole ppitch is generally taken in odd numbers only\n", + "let pole pitch = 5(for short pitch winding)\n", + "Coil Span or coil through = 1-6\n", + "No. of Coils/Pole/Phase = 1\n", + "Electrical Degrees/Slot = 30\n", + "No. of slotes required for proper phase displacement = 4\n", + "Winding Table is as follows :\n", + "Coil No. Connection Lead from Coil Span Phase and Group No.\n", + " 1 A1 1-6 A1\n", + " 2 3-8 C4\n", + " 3 B1 5-10 B1\n", + " 4 7-12 A2\n", + " 5 C1 9-14 C1\n", + " 6 11-16 B2\n", + " 7 13-18 A3\n", + " 8 15-20 C2\n", + " 9 17-22 B3\n", + " 10 19-24 A4\n", + " 11 21-2 C3\n", + " 12 23-4 B4\n" + ] + } + ], + "prompt_number": 52 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.2, Page 28" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "slots=24;#no. of slotes\n", + "P=4;#no. of poles\n", + "CoilPitch=5;#unitless\n", + "NoOfPhase=3;#no of phase\n", + "\n", + "#Calculations&Results\n", + "print \"As the winding is in double layer, each slot contains 2 coils per slide.\";\n", + "CoilSlidePerSlot=2;#coil slide per slot\n", + "CoilSlidePerCoil=2;#coil slide per Coil\n", + "TotalCoils=slots*CoilSlidePerSlot/CoilSlidePerCoil;#no. of coils\n", + "print \"Total no. of Coils : \",TotalCoils;\n", + "PolesPitch=slots/P;#unitless\n", + "print \"Pole Pitch = \",PolesPitch;\n", + "CoilsPerPolePerPhase=TotalCoils/(P*NoOfPhase);#No. of Coils/Pole/Phase\n", + "print \"No. of Coils/Pole/Phase = \",CoilsPerPolePerPhase;\n", + "pair_of_poles=2;#no. of pair of poles\n", + "TotalElectricalDegree=360*pair_of_poles;#in degree \n", + "ElectricalDegreesPerSlot=TotalElectricalDegree/slots;#in degree electrical\n", + "print \"Electrical Degrees/Slot = \",ElectricalDegreesPerSlot;\n", + "Slots_required=120/ElectricalDegreesPerSlot;#No. of slotes required for proper phase displacement\n", + "print \"No. of slots required for proper phase displacement = \",Slots_required;\n", + "print \"ie. Phase A1 is brought out from slot no. = 1\";\n", + "print \"Phase B1 at slot no. = 1+4 = 5\";\n", + "print \"Phase C1 at slot no. = 5+4 = 9\";\n", + "print \"Col Connection - end to start to start\";\n", + "print \"Winding Table is as follows :\";\n", + "print \"Coil No. Connection Lead from Coil Span Phase and Group No.\";\n", + "print \" 1 A1 1-6 A1\";\n", + "print \" 2 2-7 \";\n", + "print \" 3 3-8 C4\";\n", + "print \" 4 4-9 \";\n", + "print \" 5 B1 5-10 B1\";\n", + "print \" 6 6-11 \";\n", + "print \" 7 7-12 A2\";\n", + "print \" 8 8-13 \";\n", + "print \" 9 C1 9-14 C1\";\n", + "print \" 10 10-15 \";\n", + "print \" 11 11-16 B2\";\n", + "print \" 12 12-17 \";\n", + "print \" 13 13-18 A3\";\n", + "print \" 14 14-19 \";\n", + "print \" 15 15-20 C2\";\n", + "print \" 16 16-21 \";\n", + "print \" 17 17-22 B3\";\n", + "print \" 18 18-23 \";\n", + "print \" 19 19-24 A4\";\n", + "print \" 20 20-1 \";\n", + "print \" 21 21-2 C3\";\n", + "print \" 22 22-3 \";\n", + "print \" 23 23-4 \";\n", + "print \" 24 24-5 B4\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "As the winding is in double layer, each slot contains 2 coils per slide.\n", + "Total no. of Coils : 24\n", + "Pole Pitch = 6\n", + "No. of Coils/Pole/Phase = 2\n", + "Electrical Degrees/Slot = 30\n", + "No. of slots required for proper phase displacement = 4\n", + "ie. Phase A1 is brought out from slot no. = 1\n", + "Phase B1 at slot no. = 1+4 = 5\n", + "Phase C1 at slot no. = 5+4 = 9\n", + "Col Connection - end to start to start\n", + "Winding Table is as follows :\n", + "Coil No. Connection Lead from Coil Span Phase and Group No.\n", + " 1 A1 1-6 A1\n", + " 2 2-7 \n", + " 3 3-8 C4\n", + " 4 4-9 \n", + " 5 B1 5-10 B1\n", + " 6 6-11 \n", + " 7 7-12 A2\n", + " 8 8-13 \n", + " 9 C1 9-14 C1\n", + " 10 10-15 \n", + " 11 11-16 B2\n", + " 12 12-17 \n", + " 13 13-18 A3\n", + " 14 14-19 \n", + " 15 15-20 C2\n", + " 16 16-21 \n", + " 17 17-22 B3\n", + " 18 18-23 \n", + " 19 19-24 A4\n", + " 20 20-1 \n", + " 21 21-2 C3\n", + " 22 22-3 \n", + " 23 23-4 \n", + " 24 24-5 B4\n" + ] + } + ], + "prompt_number": 53 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.3, Page 35" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import *\n", + "\n", + "#given data\n", + "StatorSlots=36;#No. of stator slots\n", + "Poles=4;#No. of poles\n", + "#coilSpan=1:8;#unitless\n", + "\n", + "#Calculations&Results\n", + "SlotsPerPole=StatorSlots/Poles;# no. of slots per pole\n", + "print \"Slots/Pole = \",SlotsPerPole;\n", + "print \"Here the coil span falls short by, \",((2./9)*180),\" Degree\";\n", + "alfa=40*pi/180;#short pitch angle in degree\n", + "Kp=cos(alfa/2);#Coil span Factor\n", + "print \"Pitch Factor or coil span factor : \",round(Kp,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slots/Pole = 9\n", + "Here the coil span falls short by, 40.0 Degree\n", + "Pitch Factor or coil span factor : 0.9\n" + ] + } + ], + "prompt_number": 62 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.4, Page 38" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import *\n", + "\n", + "#given data\n", + "NoOfPhase=3;#no of phase\n", + "P=16;#No. of pole alternator\n", + "Slots=144;#No. of slots\n", + "Conductors=10;#per slot\n", + "fi=0.03;#in Weber\n", + "N=375;#machine speed in rpm\n", + "\n", + "#Calculations\n", + "f=P*N/120;#in Hz\n", + "SlotsPerPole=Slots/P;#unitless\n", + "m=Slots/(P*NoOfPhase);#unitless\n", + "Beta=180/(SlotsPerPole);#in Degree\n", + "Kd=sin(m*Beta/2*pi/180)/(m*sin(Beta/2*pi/180));#unitless\n", + "TotalConductors=Conductors*Slots;#no. of conductors\n", + "TotalConductorsPerPhase=Conductors*Slots/NoOfPhase;#no. of conductors/phase\n", + "TurnsPerPhase=TotalConductorsPerPhase/2;#No. of turns per phase\n", + "EMFPerPhase=4.44*Kd*fi*f*TurnsPerPhase;#in Volt\n", + "LineVoltage=sqrt(3)*EMFPerPhase;#in Volt\n", + "\n", + "#Results\n", + "print \"Frequency in Hz : \",f;\n", + "print \"Phase Electromotive force in Volt : \",round(EMFPerPhase,1);\n", + "print \"Line Electromotive force in Volt : \",round(LineVoltage,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Frequency in Hz : 50\n", + "Phase Electromotive force in Volt : 1534.1\n", + "Line Electromotive force in Volt : 2657.202\n" + ] + } + ], + "prompt_number": 76 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.5, Page 39" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "NoOfPhase=3;#no of phase\n", + "SlotsPerPhase=3;#o. of slots\n", + "\n", + "#Calculations&Results\n", + "m=SlotsPerPhase;#no. of slots\n", + "SlotsPerPolePerPhase=SlotsPerPhase*NoOfPhase;#unitless\n", + "Beta=180./SlotsPerPolePerPhase;#in degree\n", + "print \"The phase difference between the induced emf in two coils lying in adjacent slots = \",(Beta),\" degree\";\n", + "Kb=sin((m*Beta*pi)/(2*180))/(m*sin((Beta*pi)/(2*180)));#unitless\n", + "print \"The breadth factor : \",round(Kb,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The phase difference between the induced emf in two coils lying in adjacent slots = 20.0 degree\n", + "The breadth factor : 0.96\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.6, Page 39" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "StatorSlots=24#No. of stator slots\n", + "Poles=4;#No. of poles\n", + "SlotsPerPole=StatorSlots/Poles;# no. of slots per pole\n", + "#coilSpan=1:6;#unitless\n", + "\n", + "#Calculations&Results\n", + "print \"If the sides of the coil are placed in slots 1 and 7, then it is full pitched, \\nIf the coil slides are placed\" \\\n", + " + \" in 1 and 6 then it is short pitched and the distance equal to 5/6th of pole-pitch.\";\n", + "print \"Since it falls short by 1/6th of the pole-pitch, hence it is short by :\"\n", + "theta=180/SlotsPerPole;#in Degree\n", + "print (theta),\" Degree.\"\n", + "Kp=cos(theta/2*pi/180);#unitless\n", + "print \"Pitch-factor Kp : \",round(Kp,4);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "If the sides of the coil are placed in slots 1 and 7, then it is full pitched, \n", + "If the coil slides are placed in 1 and 6 then it is short pitched and the distance equal to 5/6th of pole-pitch.\n", + "Since it falls short by 1/6th of the pole-pitch, hence it is short by :\n", + "30 Degree.\n", + "Pitch-factor Kp : 0.9659\n" + ] + } + ], + "prompt_number": 2 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.7, Page 40" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "NoOfPhase=3;#no of phase\n", + "Eph=3300/sqrt(3);#in Volts\n", + "f=50;#in Hz\n", + "Poles=12;#No. of poles\n", + "StatorSlots=144#No. of stator slots\n", + "\n", + "#Calculations\n", + "SlotsPerPhase=StatorSlots/NoOfPhase;#no. of slots/phase\n", + "Conductors=5;#per slot\n", + "ConductorsPerphase=SlotsPerPhase*Conductors;#Conductors/Phase\n", + "S=ConductorsPerphase;#Conductors/phase\n", + "SlotsPerPolePerPhase=SlotsPerPhase/Poles;#no. of slots/phase\n", + "Kf=1.11;#Form Factor\n", + "Kb=0.96;#Breadth Factor\n", + "Kp=1;#For concentric winding\n", + "fi=Eph/(2*Kf*Kb*Kp*S*f);#in weber\n", + "\n", + "#Result\n", + "print \"The Flux per pole in weber : \",round(fi,4);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The Flux per pole in weber : 0.0745\n" + ] + } + ], + "prompt_number": 93 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.8, Page 41" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "Phase=3;#no. of phase\n", + "f=50;#in Hz\n", + "P=16;#No. of pole alternator\n", + "Slots=144.;#No. of slots\n", + "conductors=10;#conductors per slot\n", + "\n", + "#Calculations&Results\n", + "fi=2.48*10**-2;#in weber\n", + "n=Slots/P;#No. of slots/pole\n", + "Zr=Slots*conductors/Phase;#No. of conductors/Phase\n", + "T=Zr/2;#N. of turns/phase\n", + "Beta=180/n;#Angular displacement between slots in degree\n", + "m=n/Phase;#No. of slots/pole/Phase\n", + "Kd=sin(m*Beta/2*pi/180)/(m*sin(Beta/2*pi/180));#Distribution factor :unitless \n", + "print \"The coil span falls short of 2 slots i.e. \",(2*180/9),\" degree.\";\n", + "alfa=40;#short pitch angle in degree\n", + "Kp=(cos(alfa/2*pi/180));#Unitless\n", + "#Formula : f=P*N/120;#in Hz\n", + "N=120*f/P;#in rpm\n", + "Ep=4.44*Kd*Kp*fi*f*T;#in Volts\n", + "LineVoltage=sqrt(3)*Ep;#in Volts\n", + "print \"The speed is \",(N),\" in rpm\";\n", + "print \"The line emf is \",round(LineVoltage,2),\" Volts\";\n", + "#Note : Answer in the book is not accurate for last part due to rounding off errors" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The coil span falls short of 2 slots i.e. 40 degree.\n", + "The speed is 375 in rpm\n", + "The line emf is 2064.15 Volts\n" + ] + } + ], + "prompt_number": 109 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.9, Page 58" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "RatedPower=100;#in KVA\n", + "RatedPower=100*1000;#in VA\n", + "VL=1040;#in Volt\n", + "Phase=3;#Machine phase\n", + "If=40;#in Ampere\n", + "Isc=200;#in Ampere\n", + "EL=1040;#in Volt\n", + "\n", + "#Calculations&Results\n", + "Eph=EL/sqrt(3);#in Volt\n", + "Zs=Eph/Isc;#in Ohm\n", + "Rs=0.2;#in Ohm\n", + "Xs=sqrt(Zs**2-Rs**2);#in Ohm\n", + "IL=19.25;#in Ampere\n", + "V=3000/sqrt(3);#in Volt\n", + "#At 0.8 power factor lagging\n", + "IRa=IL*0.2;#in Volt\n", + "IXs=IL*Xs;#in Volt\n", + "Vsin_fi=V*0.6;#in Volt\n", + "Vcos_fi=V*0.8;#in Volt\n", + "Eo=sqrt((Vcos_fi+IRa)**2+(Vsin_fi+IXs)**2);#in Volts\n", + "Regulation=((Eo-V)/V)*100;#in %\n", + "print \"Full load percentage regulation at a power factor of 0.8 lagging : \",round(Regulation,1);\n", + "\n", + "#At 0.8 power factor leading\n", + "Eo=sqrt((Vcos_fi+IRa)**2+(Vsin_fi-IXs)**2);#in Volts\n", + "Regulation=((Eo-V)/V)*100;#in %\n", + "print \"Full load percentage regulation at a power factor of 0.8 leading : \",round(Regulation,2);\n", + "print \"Negative regulation due to leading power factor.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Full load percentage regulation at a power factor of 0.8 lagging : 2.2\n", + "Full load percentage regulation at a power factor of 0.8 leading : -1.78\n", + "Negative regulation due to leading power factor.\n" + ] + } + ], + "prompt_number": 110 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.10, Page 60" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "RatedPower=50;#in KVA\n", + "RatedPower=50*1000;#in VA\n", + "VL=173;#in Volts\n", + "Ra=0.1;#in Ohm\n", + "\n", + "#Calculations&Results\n", + "VP=VL/sqrt(3);#in Volts\n", + "print \"Some exciting curent on short circuit produces a current of 100 A.\";\n", + "OC_PhaseVoltage=100;#in Volt\n", + "SC_Current=100;#in Ampere\n", + "Zs=OC_PhaseVoltage/SC_Current;#n ohm\n", + "Xs=sqrt(Zs**2-Ra**2);#in Ohm\n", + "print \"Impedence of the alternator in Ohm : \",round(Xs,2);\n", + "V=400;#in Volts\n", + "I_FL=RatedPower/(sqrt(3)*V);#in Ampere\n", + "V=400/sqrt(3);#in Volts\n", + "Eo=sqrt((V+I_FL*Ra)**2+(I_FL*Xs)**2);#in Volts\n", + "Regulation=(Eo-V)*100/V;#in %\n", + "print \"Regulation at U.P.F. in % :\",round(Regulation,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Some exciting curent on short circuit produces a current of 100 A.\n", + "Impedence of the alternator in Ohm : 0.99\n", + "Regulation at U.P.F. in % : 7.71\n" + ] + } + ], + "prompt_number": 112 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.11, Page 61" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=500;#in KVA\n", + "OutputPower=500*1000;#in VA\n", + "VL=3300;#in Volts\n", + "Ra=0.3;#in Ohm\n", + "Xs=4;#in Ohm\n", + "PF=0.8;#Lagging Power factor\n", + "\n", + "#Calculations&Results\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "IL=OutputPower/(sqrt(3)*VL);#in Ampere\n", + "print \"For a star connected alternator, line current is equal to phase current. Therefore Ia=IL\";\n", + "Ia=IL;#in Ampere\n", + "#PF=cosd(fi)=0.8 and sind(fi)=0.6\n", + "cos_fi=0.8;#Power factor\n", + "sin_fi=0.6;#Unitless\n", + "VPerPhase=VL/sqrt(3);#in Volts\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);#in Volts/Phase\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "print \"Voltage Regulation at Full Load in % :\",round(Regulation,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For a star connected alternator, line current is equal to phase current. Therefore Ia=IL\n", + "Voltage Regulation at Full Load in % : 12.98\n" + ] + } + ], + "prompt_number": 113 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.12, Page 62" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "V=2000;#in Volt\n", + "Ia=100;#in Ampere\n", + "OC_Voltage=500;#in Volt\n", + "SC_Current=100;#in Ampere\n", + "\n", + "#Calculations&Results\n", + "Zs=OC_Voltage/SC_Current;#in ohm\n", + "Ra=0.8;#in Ohm\n", + "Xs=sqrt(Zs**2-Ra**2);#in Ohm\n", + "#formula : Induced EMF, E=sqrt((V*cos_fi+Ia*Ra)^2+(V*sin_fi+Ia*Xs)^2)\n", + "#Part (a) : at unity pf\n", + "cos_fi=1;#Unitless\n", + "sin_fi=0;#Unitless\n", + "E=sqrt((V*cos_fi+Ia*Ra)**2+(V*sin_fi+Ia*Xs)**2)\n", + "Regulation=(E-V)*100/V;#in %\n", + "print \"Regulation at U.P.F. in % :\",round(Regulation,2);\n", + "\n", + "#Part (b) : at 0.71 pf lagging\n", + "cos_fi=0.71;#Unitless\n", + "sin_fi=0.704;#Unitless\n", + "E=sqrt((V*cos_fi+Ia*Ra)**2+(V*sin_fi+Ia*Xs)**2)\n", + "Regulation=(E-V)*100/V;#in %\n", + "print \"Regulation at 0.71 pf lagging in % :\",round(Regulation,2);\n", + "\n", + "#Part (c) : at 0.8 pf leading\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((V*cos_fi+Ia*Ra)**2+(V*sin_fi-Ia*Xs)**2)\n", + "Regulation=(E-V)*100/V;#in %\n", + "print \"Regulation at 0.8 pf leading in % :\",round(Regulation,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Regulation at U.P.F. in % : 6.89\n", + "Regulation at 0.71 pf lagging in % : 21.1\n", + "Regulation at 0.8 pf leading in % : -8.88\n" + ] + } + ], + "prompt_number": 114 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.13, Page 63" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#given data\n", + "Ia=60;#in Ampere\n", + "\n", + "#Calculations&Results\n", + "print \"The value of synchronous impedence at this excitation :\" ;\n", + "OC_Voltage=900;#in Volt\n", + "SC_Current=150;#in Ampere\n", + "Zs=OC_Voltage/SC_Current;#in ohm\n", + "print \"Zs equals to \",(Zs),\" Ohm\";\n", + "print \"Internal Voltage drop when tthe load current is 60A=Ia*Zs=\",(Ia*Zs),\" Volts\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The value of synchronous impedence at this excitation :\n", + "Zs equals to 6 Ohm\n", + "Internal Voltage drop when tthe load current is 60A=Ia*Zs= 360 Volts\n" + ] + } + ], + "prompt_number": 118 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.14, Page 63" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "V=6600;#in Volts\n", + "OutputPower=2000;#in KVA\n", + "\n", + "#Calculations\n", + "OutputPower=2000*1000;#in VA\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "IL=OutputPower/(sqrt(3)*V);#in Ampere\n", + "Ia=IL;#in Ampere\n", + "Ra=0.4;#in Ohm\n", + "Xs=4.5;#in Ohm\n", + "#PF=cosd(fi)=0.8 and sind(fi)=0.6\n", + "cos_fi=0.8;#Power factor\n", + "sin_fi=0.6;#Unitless\n", + "VPerPhase=V/sqrt(3);#in Volts\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);#in Volts/Phase\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "\n", + "#Result\n", + "print \"Percentage Change in Terminal Voltage :\",round(Regulation,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Percentage Change in Terminal Voltage : 14.9\n" + ] + } + ], + "prompt_number": 119 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.15, Page 63" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=1200;#in KVA\n", + "OutputPower=1200*1000;#in VA\n", + "V=3300;#in Volt\n", + "Ra=0.25;#in Ohm\n", + "\n", + "#Calculations&Results\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "IL=OutputPower/(sqrt(3)*V);#in Ampere\n", + "Ia=IL;#in Ampere\n", + "VPerPhase=V/sqrt(3);#in Volts\n", + "OC_Voltage=1100;#in Volt\n", + "SC_Current=200;#in Ampere\n", + "Zs=OC_Voltage/(sqrt(3)*SC_Current);#in ohmRa\n", + "Xs=sqrt(Zs**2-Ra**2);#in Ohm\n", + "#formula : Induced EMF, E=sqrt((V*cos_fi+Ia*Ra)^2+(V*sin_fi+Ia*Xs)^2)\n", + "\n", + "#Part (a) : For lagging pf load\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "print \"Regulation at U.P.F. in % :\",round(Regulation,2);\n", + "\n", + "#Part (b) : For leading pf load\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi-Ia*Xs)**2)\n", + "Regulation=(E-VPerPhase)*100/VPerPhase;#in %\n", + "print \"Regulation at 0.71 pf lagging in % :\",round(Regulation,2);\n", + "#Note: For leading power factor load, the regulation s negative. " + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Regulation at U.P.F. in % : 25.9\n", + "Regulation at 0.71 pf lagging in % : -13.52\n" + ] + } + ], + "prompt_number": 120 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.16, Page 64" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=1500;#in KVA\n", + "OutputPower=1500*1000;#in VA\n", + "V=6600;#in Volt\n", + "Ra=0.4;#in Ohm\n", + "Xs=6;#in Ohm per phase\n", + "pf=0.8;#lagging power factor\n", + "\n", + "#Calculations&Results\n", + "#Formula : outputPower=sqrt(3)*VL*IL\n", + "Ia=OutputPower/(sqrt(3)*V);#in Ampere\n", + "VPerPhase=V/sqrt(3);#in Volts\n", + "#formula : Induced EMF, E=sqrt((V*cos_fi+Ia*Ra)^2+(V*sin_fi+Ia*Xs)^2)\n", + "cos_fi=0.8;#Unitless\n", + "sin_fi=0.6;#Unitless\n", + "E=sqrt((VPerPhase*cos_fi+Ia*Ra)**2+(VPerPhase*sin_fi+Ia*Xs)**2);#in volt\n", + "print \"Induced emf in volt : \",round(E,2);\n", + "print \"As excitation remains constant, E at 4364 volt remains constant.\";\n", + "E=4364;#in Volt\n", + "V=4743;#in Volts\n", + "TerminalVoltage=sqrt(3)*V;#in Volts\n", + "print \"Terminal voltage line to line in Volts : \",round(TerminalVoltage);\n", + "#Note ans of 1st part is wrong in the books" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Induced emf in volt : 4366.07\n", + "As excitation remains constant, E at 4364 volt remains constant.\n", + "Terminal voltage line to line in Volts : 8215.0\n" + ] + } + ], + "prompt_number": 123 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.17, Page 81" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "OutputPower=2500;#in KVA\n", + "OutputPower=2500*1000;#in VA\n", + "V=6600;#in Volt\n", + " \n", + "#Calculations\n", + "#For first load i.e. Lighting load : \n", + "KW1=2500;#in KWatts\n", + "KVAR1=0;#Kwatts\n", + "#For second load i.e. Motor load : \n", + "KW2=5000;#in KWatts\n", + "cos_fi=0.707;#unitless\n", + "sin_fi=0.707;#unitless\n", + "KVAR2=KW2*sin_fi/cos_fi;#Kwatts\n", + "#For total load\n", + "TotalKW=KW1+KW2;#in KWatts\n", + "TotalKVAR=KVAR1+KVAR2;#Kwatts \n", + "#For first Machine\n", + "KWm=4000;#in KWatts\n", + "cos_fi=0.8;#unitless\n", + "sin_fi=0.6;#unitless\n", + "KVARm=KWm*sin_fi/cos_fi;#Kwatts\n", + "#so, second machine will be supplying\n", + "KW=TotalKW-KWm;#in Kwatts\n", + "print \"KW output of second machine : \",KW;\n", + "KVAR=TotalKVAR-KVARm;#in KWatts\n", + "tan_fi=KVAR/KW;#unitless\n", + "fi=atan(tan_fi);#ib degree\n", + "#Power factor of other machine\n", + "pf=cos(fi);#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of other machine : \",round(pf,4),\" lagging\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "KW output of second machine : 3500\n", + "Power factor of other machine : 0.8682 lagging\n" + ] + } + ], + "prompt_number": 126 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.18, Page 82" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "\n", + "#Load1 : \n", + "KW1=500;#in KWatts\n", + "KVAR1=0;#Kwatts\n", + "\n", + "#Load2 : \n", + "KW2=1000;#in KWatts\n", + "pf=0.9;#lagging\n", + "cos_fi=0.9;#unitless\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR2=KW2*sin_fi/cos_fi;#Kwatts\n", + "\n", + "#Load3 : \n", + "KW3=800;#in KWatts\n", + "pf=0.8;#lagging\n", + "cos_fi=0.8;#unitless\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR3=KW3*sin_fi/cos_fi;#Kwatts\n", + "\n", + "#Load4 : \n", + "KW4=500;#in KWatts\n", + "pf=0.9;#lagging\n", + "cos_fi=0.9;#unitless\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR4=-KW4*sin_fi/cos_fi;#Kwatts\n", + "#TOtalKW and TotalKVAR\n", + "TotalKW=KW1+KW2+KW3+KW4;#in KWatts\n", + "TotalKVAR=KVAR1+KVAR2+KVAR3+KVAR4;#in KWAtts\n", + "#For the first Machine :\n", + "KW=1500;#n Kwatts\n", + "cos_fi=0.95;#unitless\n", + "sin_fi=0.3123;#unitless\n", + "KVAR=KW*sin_fi/cos_fi;#Kwatts\n", + "\n", + "KW1=TotalKW-KW;#in KWatts\n", + "KVAR1=TotalKVAR-KVAR;#in Volts\n", + "print \"KW supplied by other machine : \",(TotalKW-KW);\n", + "print \"KVAR supplied by other machine : \",round(TotalKVAR-KVAR); #answer differs due to rounding off the digits\n", + "tan_fi=KVAR1/KW1;#unitless\n", + "#fi=atand(tan_fi);#in degree\n", + "cos_fi=cos(atan(tan_fi));#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of the other machine : \",round(cos_fi,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "KW supplied by other machine : 1300\n", + "KVAR supplied by other machine : 349.0\n", + "Power factor of the other machine : 0.966\n" + ] + } + ], + "prompt_number": 129 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.19, Page 82" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "\n", + "#Lighting Load : \n", + "MW1=20;#load in Mwatts\n", + "KW1=MW1*1000;#in KWatts\n", + "KVAR1=0;#Kwatts\n", + "\n", + "#Motor Load : \n", + "MW2=40;#load in Mwatts\n", + "KW2=MW2*1000;#in KWatts\n", + "pf=0.8;#unitless\n", + "cos_fi=0.8\n", + "fi=acos(pf);\n", + "sin_fi=sin(fi);#unitless\n", + "KVAR2=KW2*sin_fi/cos_fi;#Kwatts\n", + "\n", + "#For Total Load : \n", + "TotalKW=KW1+KW2;#load in Mwatts\n", + "TotalKVAR=KVAR1+KVAR2;#in KWatts\n", + "#For first machine : \n", + "MWm=32;#load in Mwatts\n", + "KWm=MWm*1000;#in KWatts\n", + "cos_fi=0.866;#unitless\n", + "fi=acos(cos_fi);\n", + "tan_fi=tan(fi);#unitless\n", + "KVARm=KWm*tan_fi;#in KWatts\n", + "#so, load supplied by the second machine\n", + "KW2=TotalKW-KWm;#in Kwatts\n", + "print \"Load of other machine,KW : \",(KW2);\n", + "KVAR2=TotalKVAR-KVARm;#in Kwatts\n", + "tan_fi=KVAR2/KW2;#unitless\n", + "fi=atan(tan_fi);#in degree\n", + "cos_fi=cos(atan(tan_fi));#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of the other machine : \",round(cos_fi,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Load of other machine,KW : 28000\n", + "Power factor of the other machine : 0.925\n" + ] + } + ], + "prompt_number": 131 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.20, Page 84" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "cos_fi=0.8;#unitless\n", + "fi=acos(cos_fi);\n", + "tan_fi=tan(fi);#unitless\n", + "\n", + "#For Alternator A : \n", + "cos_fi_A=0.9;#unitless\n", + "fi_A=acos(cos_fi_A);\n", + "tan_fi_A=tan(fi_A);#unitless\n", + "#Formula : Active load, KW=V*I*cos_fi\n", + "#Formula : Reactive load, KVAR=V*I*sin_fi\n", + "ActiveLoad=8000;#in KW\n", + "ReactiveLoad=ActiveLoad*tan_fi;#in KVAR\n", + "\n", + "#For A:\n", + "ActiveLoadA=5000;#in KW\n", + "ReactiveLoadA=ActiveLoadA*tan_fi_A;#in KVAR\n", + "\n", + "#For B :\n", + "ActiveLoadB=ActiveLoad-ActiveLoadA;#in KW\n", + "ReactiveLoadB=ReactiveLoad-ReactiveLoadA;#in KVAR\n", + "tan_fi_B=ReactiveLoadB/ActiveLoadB;#unitless\n", + "fi_B=atan(tan_fi_B);#in degree\n", + "cos_fi=cos(atan(tan_fi_B));#unitless\n", + "\n", + "#Result\n", + "print \"Power factor of the other machine : \",round(cos_fi,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Power factor of the other machine : 0.642\n" + ] + } + ], + "prompt_number": 133 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.21, Page 97" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "V=6600;#in Volts\n", + "KW=6000.;#in KWatts\n", + "pf=0.8;#unitless\n", + "cos_fi=pf;#unitless\n", + "Eff=90.;#in %\n", + "\n", + "#Calculations&Results\n", + "#Part (a) : \n", + "KVA=KW/cos_fi;#in KVAR\n", + "print \"KVA rating of the alternator : \",(KVA),\" KVA\";\n", + "#Part (b) : \n", + "TotalRating=KVA;#in KVA\n", + "VA=TotalRating*1000;#in VA\n", + "I=VA/(sqrt(3)*KW);#in Ampere\n", + "print \"Current Rating in Ampere : \",round(I,2);\n", + "#Part (c) :\n", + "Input=KW/(Eff/100);#in KW\n", + "print \"Power Input(in KW) :\",round(Input,2);\n", + "Input=Input*1000/735.5;#in hp\n", + "print \"Power Input(in hp) :\",round(Input,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "KVA rating of the alternator : 7500.0 KVA\n", + "Current Rating in Ampere : 721.69\n", + "Power Input(in KW) : 6666.67\n", + "Power Input(in hp) : 9064.13\n" + ] + } + ], + "prompt_number": 139 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.22, Page 97" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#given data\n", + "Ecoil=8000;#in Volts\n", + "Icoil=418;#in Ampere\n", + "pf=80.;#in % lgging\n", + "\n", + "#Calculations&Results\n", + "pf=pf/100;#in fraction\n", + "cos_fi=pf;#unitless\n", + "#Part (i) : \n", + "EL=sqrt(3)*Ecoil;#in volt\n", + "print \"Line volts(in V): \",round(EL);\n", + "#Part (ii) : \n", + "IL=Icoil;#in Ampere\n", + "print \"Line Current in Ampere : \",IL;\n", + "#Part (iii) :\n", + "Rating=sqrt(3)*EL*IL/1000;#in KVA\n", + "print \"Rating (in KVA) :\",Rating;\n", + "#Part (iv) :\n", + "FullLoadPower=sqrt(3)*EL*IL*cos_fi/1000;#in KW\n", + "print \"Full Load Power in KW :\",FullLoadPower;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Line volts(in V): 13856.0\n", + "Line Current in Ampere : 418\n", + "Rating (in KVA) : 10032.0\n", + "Full Load Power in KW : 8025.6\n" + ] + } + ], + "prompt_number": 140 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/chapter_2.ipynb b/Electrical_Machines_II/chapter_2.ipynb new file mode 100755 index 00000000..c18790d9 --- /dev/null +++ b/Electrical_Machines_II/chapter_2.ipynb @@ -0,0 +1,2810 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:e3e9de8b11024daf372011e6906cf967e4f06ac156a70592b24845e8bc2b4696" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 2: Induction Motors" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.1, Page 130" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=2;#no. of poles\n", + "f=60;#in Hz\n", + "N=3460.;#in rpm\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "Slip=Ns-N;#in rpm\n", + "PercentageSlip=((Ns-N)/Ns)*100;#in %\n", + "\n", + "#Results\n", + "print \"Synchronous speed in rpm : \",Ns;\n", + "print \"Percentage Slip : \",round(PercentageSlip,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Synchronous speed in rpm : 3600\n", + "Percentage Slip : 3.89\n" + ] + } + ], + "prompt_number": 124 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.2, Page 131" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=8.;#no. of poles\n", + "f=50.;#in Hz\n", + "N=700.;#in rpm\n", + "\n", + "#Calculations&Results\n", + "#Part (a) :\n", + "Ns=120*f/P;#in rpm\n", + "print \"Synchronous speed in rpm : \",Ns;\n", + "#Part (b) :\n", + "S=(Ns-N)*100/Ns;#in %\n", + "print \"Slip(in %) : \",round(S,2);\n", + "#Part (c) :\n", + "#At the time of stsrt S=1;\n", + "fdash=f;#in Hz\n", + "print \"Rotor frequeny at the time of starting(in Hz) : \",fdash;\n", + "#Part (d) :\n", + "fdash=(S/100)*f;#in Hz\n", + "print \"Rotor frequeny at the given speed(in Hz) : \",round(fdash,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Synchronous speed in rpm : 750.0\n", + "Slip(in %) : 6.67\n", + "Rotor frequeny at the time of starting(in Hz) : 50.0\n", + "Rotor frequeny at the given speed(in Hz) : 3.3\n" + ] + } + ], + "prompt_number": 125 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.3, Page 131" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=10;#no. of poles\n", + "f=50;#in Hz\n", + "N=600;#in rpm\n", + "Nr=1440;#in rpm\n", + "\n", + "#Calculations&Results\n", + "f=P*N/120;#in Hz\n", + "#When P=2\n", + "P=2;#no. of poles\n", + "Ns=120*f/P;#in rpm\n", + "#When P=4\n", + "P=4;#no. of poles\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-Nr)*100/Ns;#Slip in %\n", + "print \"Percentage Slip(in %) : \",(S);\n", + "print \"No. of poles : \",(S);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Percentage Slip(in %) : 4\n", + "No. of poles : 4\n" + ] + } + ], + "prompt_number": 126 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.4, Page 132" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=8;#no. of poles\n", + "f=50;#in Hz\n", + "fr=1.5;#in Hz\n", + "\n", + "#Calculations\n", + "#Formula : fr=S*f\n", + "S=fr/f;#slip(unitless)\n", + "Ns=120*f/P;#in rpm\n", + "#Formula : S=(Ns-Nr)/Ns\n", + "Nr=Ns-S*Ns;#in rpm\n", + "\n", + "#Results\n", + "print \"Motor running speed in rpm : \",round(Nr);\n", + "print \"Slip(in %):\",(S*100);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Motor running speed in rpm : 728.0\n", + "Slip(in %): 3.0\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.5, Page 133" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=2;#no. of poles\n", + "f=50;#in Hz\n", + "S=4.;#in %\n", + "\n", + "#Calculations&Results\n", + "#Part (i) : \n", + "Ns=f/P;#in rps\n", + "Ns=Ns*60;#in rpm\n", + "print \"Synchronous speed in rpm : \",Ns;\n", + "#Part (ii) :\n", + "print \"Slip = \",(S),\"% or \",(S/100);\n", + "#Part (iii) :\n", + "N=Ns*(1-S/100);#in rpm\n", + "print \"Actual speed in rpm :\",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Synchronous speed in rpm : 1500\n", + "Slip = 4.0 % or 0.04\n", + "Actual speed in rpm : 1440.0\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.6, Page 133" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=4;#no. of poles\n", + "f1=50;#in Hz\n", + "fdash=1.5;#in Hz\n", + "\n", + "#Calculations\n", + "S=fdash/f1;#unitless\n", + "print \"Slip : \",S;\n", + "Ns=120*f1/P;#in rpm\n", + "N=Ns*(1-S);#in rpm\n", + "\n", + "#Result\n", + "print \"Speed of motor in rpm :\",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.03\n", + "Speed of motor in rpm : 1455.0\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.7, Page 138" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "N=1440;#in rpm\n", + "\n", + "#Calculations&Results\n", + "f=50;#in Hz\n", + "print \"For this speed of 1440 rpm the synchronous speed may be either 1500 rpm or 3000 rpm.\";\n", + "Ns=1500;#in rpm\n", + "P=4;#no. of poles\n", + "print \"No. of poles : \",P;\n", + "Slip=(Ns-N)*100/Ns;#in %\n", + "print Slip,\"Slip(in %) : \",Slip;\n", + "speed1=Ns-N;#Speed of rotor flux with respect to rotor\n", + "print \"Speed of rotor flux with respect to rotor in rpm : \",speed1;\n", + "speed2=(speed1)+N;#Speed of rotor flux with respect to stator\n", + "print \"Speed of rotor flux with respect to stator in rpm :\",speed2;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For this speed of 1440 rpm the synchronous speed may be either 1500 rpm or 3000 rpm.\n", + "No. of poles : 4\n", + "4 Slip(in %) : 4\n", + "Speed of rotor flux with respect to rotor in rpm : 60\n", + "Speed of rotor flux with respect to stator in rpm : 1500\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.8, Page 140" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "E2=100;#in volt\n", + "R2=0.05;#in ohm\n", + "X2=0.1;#im ohm\n", + "\n", + "#Calculations&Results\n", + "E2perphase=E2/sqrt(3);#in volt\n", + "#part (a) : \n", + "S=0.04;#slip\n", + "I2=S*E2perphase/sqrt(R2**2+(S*X2)**2);#in Ampere\n", + "print \"At 4% slip, Rotor current in Ampere : \",round(I2,2);\n", + "fi2=degrees(acos(R2/sqrt(R2**2+(S*X2)**2)));#in degree\n", + "print \"At 4% slip, Phase angle between rotor voltage and rotor current in degree :\",round(fi2,2);\n", + "#part (b) : \n", + "S=1;#slip\n", + "I2=S*E2perphase/sqrt(R2**2+(S*X2)**2);#in Ampere\n", + "print \"At 100% slip, Rotor current in Ampere : \",round(I2,2);\n", + "fi2=degrees(acos(R2/sqrt(R2**2+(S*X2)**2)));#in degree\n", + "print \"At 100% slip, Phase angle between rotor voltage and rotor current in degree :\",round(fi2,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "At 4% slip, Rotor current in Ampere : 46.04\n", + "At 4% slip, Phase angle between rotor voltage and rotor current in degree : 4.57\n", + "At 100% slip, Rotor current in Ampere : 516.4\n", + "At 100% slip, Phase angle between rotor voltage and rotor current in degree : 63.43\n" + ] + } + ], + "prompt_number": 137 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.9, Page 141" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "E2=100.;#in volt\n", + "R2=0.4;#in ohm\n", + "X2=2.25;#im ohm\n", + "\n", + "#Calculations\n", + "E2perphase=E2/sqrt(3);#in volt\n", + "S=4.;#in %\n", + "E=(S/100)*E2perphase;#rotor induced emf at a slip=4% in volt\n", + "Z2=sqrt(R2**2+((S/100)*X2)**2);\n", + "I2=E/Z2;#in Ampere\n", + "\n", + "#Result\n", + "print \"Rotor current in Ampere : \",round(I2,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor current in Ampere : 5.63\n" + ] + } + ], + "prompt_number": 138 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.10, Page 141" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "f=50;#in Hz\n", + "R2=0.2;#in ohm\n", + "X2=2;#im ohm\n", + "E2=60;#in volt\n", + "\n", + "#Calculations&Results\n", + "E2perphase=E2/sqrt(3);#in volt\n", + "#case (i) : S=1 \n", + "S=1;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"At standstill : Rotor current in Ampere : \",round(I2,2);\n", + "#case (ii) : S=0.1 \n", + "S=0.1;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"for 10% slip : Rotor current in Ampere : \",round(I2,2);\n", + "#case (iii) : S=0.05 \n", + "S=0.05;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"for 5% slip : Rotor current in Ampere : \",round(I2,2);\n", + "#case (iv) : S=0.01 \n", + "S=0.01;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"for 1% slip : Rotor current in Ampere : \",round(I2,2);\n", + "#Note : Answer in the book is wrong for S=0.05 and S=0.01" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "At standstill : Rotor current in Ampere : 17.23\n", + "for 10% slip : Rotor current in Ampere : 12.25\n", + "for 5% slip : Rotor current in Ampere : 7.75\n", + "for 1% slip : Rotor current in Ampere : 1.72\n" + ] + } + ], + "prompt_number": 139 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.11, Page 142" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "R2=0.5;#in ohm\n", + "X=4.5;#im ohm\n", + "E=50;#line voltage in volt\n", + "\n", + "#Calculations&Results\n", + "Eperphase=E/sqrt(3);#in volt\n", + "#part (a) : Sliprings are short circuited\n", + "Z2=sqrt(R2**2+X**2);#in ohm\n", + "I2=Eperphase/Z2;#in Ampere\n", + "print \"Sliprings are short circuited, Rotor current in Ampere : \",round(I2,2);\n", + "cosfi2=R2/Z2;#unitless\n", + "print \"Power factor : \",round(cosfi2,2);\n", + "#part (b) : Sliprings are connected \n", + "R2=0.4;#in ohm\n", + "Z2=sqrt(R2**2+X**2);#in ohm\n", + "I2=Eperphase/Z2;#in Ampere\n", + "print \"Sliprings are short circuited, Rotor current in Ampere : \",round(I2,2);\n", + "cosfi2=R2/Z2;#unitless\n", + "print \"Power factor : \",round(cosfi2,2);\n", + "#Note : answer of part (a) I2 is not curate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Sliprings are short circuited, Rotor current in Ampere : 6.38\n", + "Power factor : 0.11\n", + "Sliprings are short circuited, Rotor current in Ampere : 6.39\n", + "Power factor : 0.09\n" + ] + } + ], + "prompt_number": 140 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.12, Page 146" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "R2=0.02;#in ohm\n", + "X2=0.1;#im ohm\n", + "\n", + "#Calculations&Results\n", + "#let external resistance per phase = r then R2=R2+r ohm\n", + "S=1;#slip at starting\n", + "print \"Since at start speed is zero and slip is, therefore, unity or R2=X2\";\n", + "r=X2-R2;#in ohm\n", + "print \"External resistance per phase added to the rotor circuit in ohms : \",r;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Since at start speed is zero and slip is, therefore, unity or R2=X2\n", + "External resistance per phase added to the rotor circuit in ohms : 0.08\n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.13, Page 147" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "P=4;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.03;#in ohm\n", + "X2=0.12;#im ohm\n", + "\n", + "#Calculations&Results\n", + "Smax=R2/X2;#unitless\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-Smax);#in rpm\n", + "print \"Speed corresponding to maximum torque(in rpm) :\",N;\n", + "#alfa=(R2+r)/X2\n", + "#Ratio of starting torque to max torque=75/100=2*alfa/(alfa^2+1)\n", + "#It gives : 3*alfa^2-8*alfa+3=0\n", + "#coefficients :\n", + "a=3;b=-8;c=3;\n", + "alfa1=(-b+sqrt(b**2-4*a*c))/(2*a);alfa2=(-b-sqrt(b**2-4*a*c))/(2*a);\n", + "print \"rejecting higher values, alfa = \",round(alfa2,4);\n", + "fdash=1.5;#in Hz\n", + "r=alfa2*X2-R2;#in ohm;\n", + "print \"External rotor resistance per phase(in ohm) : \",round(r,4);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed corresponding to maximum torque(in rpm) : 1125.0\n", + "rejecting higher values, alfa = 0.4514\n", + "External rotor resistance per phase(in ohm) : 0.0242\n" + ] + } + ], + "prompt_number": 141 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.14, Page 147" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=4;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.024;#in ohm\n", + "X2=0.6;#in ohm \n", + "\n", + "#Calculations&Results\n", + "Sm=R2/X2;#slip corresponding to max torque\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-Sm);#in rpm\n", + "print \"Speed corresponding to maximum torque in rpm : \",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed corresponding to maximum torque in rpm : 1440.0\n" + ] + } + ], + "prompt_number": 26 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.15, Page 148" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=6;#no. of poles\n", + "f=50;#in Hz\n", + "Sf=4.;#in %\n", + "R2=0.01;#in ohm\n", + "X2=0.05;#im \n", + "\n", + "#Calculations&Results\n", + "Ratio=((R2/X2)**2+(Sf/100)**2)/(2*R2*(Sf/100)/X2);#ratio of max torque to full load torque\n", + "print \"Maximum torque, Tmax=\",(Ratio),\"Tf\";\n", + "Sm=R2/X2;#slip corresponding to max torque\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-Sm);#in rpm\n", + "print \"Speed corresponding to maximum torque in rpm : \",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum torque, Tmax= 2.6 Tf\n", + "Speed corresponding to maximum torque in rpm : 800.0\n" + ] + } + ], + "prompt_number": 28 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.16, Page 149" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=12;#no. of poles\n", + "f=50.;#in Hz\n", + "R2=0.03;#in ohm\n", + "X2=0.5;#im \n", + "\n", + "#Calculations&Results\n", + "S=R2/X2;#unitless\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-S);#in rpm\n", + "print \"Speed at maximum torque in rpm : \",N;\n", + "S=(Ns-495)/Ns;#slip at 495 rpm speed\n", + "Ratio=(2*R2*S/X2)/((R2/X2)**2+S**2);#ratio of max torque to full load torque\n", + "print \"Ratio of full load torque to max torque(in %) :\",round((Ratio*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed at maximum torque in rpm : 470.0\n", + "Ratio of full load torque to max torque(in %) : 32.4\n" + ] + } + ], + "prompt_number": 142 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.17, Page 149" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "P=12;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.5;#in ohm\n", + "N=475.;#in rpm\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#unitless\n", + "X2=R2/S;#in ohm\n", + "#\"At maximum torque, Tmax is proportional to 1/2*X2 or 1/20\";\n", + "# \"Ts is proportional to (R2+R)/((R2+R)^2+X2^2)\";\n", + "# \"Also, Ts is proportional to 1(20*3)\";\n", + "# \"Equating the two eqn we have : (0.5*R)/(R2+R)^2+X2^2\";\n", + "#R^2-59R+70.25=0 : \n", + "#coefficients :\n", + "a=1;b=-59;c=70.25;\n", + "R=(-b-sqrt(b**2-4*a*c))/(2*a);\n", + "print \"Resistance R to be inserted in series (in ohm) :\",round(R,1);\n", + "#Note : answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Resistance R to be inserted in series (in ohm) : 1.2\n" + ] + } + ], + "prompt_number": 143 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.18, Page 150" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=4;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.04;#in ohm\n", + "N=1200.;#in rpm(speed at max torque)\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm(synchronous speed)\n", + "S=(Ns-N)/Ns;#unitless\n", + "X2=R2/S;#in ohm\n", + "#Starting torque is the torque devloped when S=1: Tst=K*R2/(R2^2+X2^2)\n", + "#Let say, m=R2/(R2^2+X2^2) then Tst=K*m\n", + "m=R2/(R2**2+X2**2);#assumed\n", + "\n", + "#Results\n", + "print \"Starting torque, Tst=\",round(m,2),\"k\";\n", + "print \"Maximum torque, Tm=K/\",(2*X2);\n", + "print \"Thus, Tst in terms of Tm can be expressed as : \";\n", + "print \"Tst/Tm=\",(0.96*0.4),\" or Tst=\",(0.96*0.4),\"Tm\";\n", + "print \"Therefore, staring torque is \",(0.96*0.4*100),\"% of maximum torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Starting torque, Tst= 0.96 k\n", + "Maximum torque, Tm=K/ 0.4\n", + "Thus, Tst in terms of Tm can be expressed as : \n", + "Tst/Tm= 0.384 or Tst= 0.384 Tm\n", + "Therefore, staring torque is 38.4 % of maximum torque.\n" + ] + } + ], + "prompt_number": 144 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.19, Page 155" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=6;#no. of poles\n", + "f=50;#in Hz\n", + "fr=3.;#in Hz\n", + "R2=0.5;#in ohm\n", + "N=475;#in rpm\n", + "RotorInput=111.9;#in KW\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=fr*100/f;#unitless\n", + "print \"% Slip : \",S;\n", + "N=Ns*(1-S/100);#in rpm\n", + "print \"Speed of motor in rpm : \",N;\n", + "RotorCopperLoss=RotorInput*S/100;#in KW\n", + "RotorCopperLoss=RotorCopperLoss/3;#in KW/Phase\n", + "print \"Rotor Copper Loss per phase(in KW) : \",RotorCopperLoss;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "% Slip : 6.0\n", + "Speed of motor in rpm : 940.0\n", + "Rotor Copper Loss per phase(in KW) : 2.238\n" + ] + } + ], + "prompt_number": 146 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.20, Page 156" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3.;#no. of phase\n", + "Pin_stator=60.;#Power input of stator in KW\n", + "TotalStatorLosses=1.;#in KW\n", + "\n", + "#Calculations&Results\n", + "Pin_rotor=Pin_stator-TotalStatorLosses;#Power input of rotor in KW\n", + "S=3.;#slip in %\n", + "RotorCopperLosses=(S/100)*Pin_rotor;#in KW\n", + "RotorCopperLosses=RotorCopperLosses/phase;#in KW per phase\n", + "print \"Rotor Copper Losses per phase(in watts) : \",(RotorCopperLosses*10**3);\n", + "TotalMechPowerDev=Pin_rotor*(1-S/100);#in KW\n", + "print \"Total mechanial power developed(in KW) : \",TotalMechPowerDev;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor Copper Losses per phase(in watts) : 590.0\n", + "Total mechanial power developed(in KW) : 57.23\n" + ] + } + ], + "prompt_number": 40 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.21, Page 156" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=6;#no. of poles\n", + "S=4.;#slip in %\n", + "Output=20;#in KW\n", + "\n", + "#Calculations&Results\n", + "FrictionalLoss=250;#in watts\n", + "f=50;#in Hz\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns-Ns*(S/100);#in rpm\n", + "print \"Speed of motor(in rpm) : \",N;\n", + "MechPowerDeveloped=Output*10**3+FrictionalLoss;#in Watts\n", + "S=S/100;#unitless\n", + "RotorCopperLoss=(S/(1-S))*MechPowerDeveloped;#in watts\n", + "print \"Rotor Copper Loss(in watts) : \",RotorCopperLoss;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed of motor(in rpm) : 960.0\n", + "Rotor Copper Loss(in watts) : 843.75\n" + ] + } + ], + "prompt_number": 147 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.22, Page 157" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "Pin_rotor=80;#Power input of rotor in KW\n", + "f=50;#in Hz\n", + "fdash=100;#alternations per minute\n", + "\n", + "#Calculations&Results\n", + "fdash=100./60;#in Hz\n", + "S=fdash/f;#unitless\n", + "print \"Slip : \",round(S,4);\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "N=Ns-Ns*S;#in rpm\n", + "print \"Rotor Speed(in rpm) : \",round(N);\n", + "MechPowerDev=Pin_rotor*(1-S)*10**3/735.5;#in H.P.\n", + "print \"Total mechanial power developed(in KW) : \",round(MechPowerDev,2);\n", + "RotorCopperLoss=S*Pin_rotor*10**3;#in watts\n", + "print \"Rotor Copper Loss(in watts) : \",RotorCopperLoss;\n", + "CopperLoss=RotorCopperLoss/phase;#in watts/phase\n", + "print \"Copper Loss per phase(in watts) : \",round(CopperLoss);\n", + "I2=60;#in Ampere\n", + "R2=CopperLoss/I2**2;#in ohm\n", + "print \"Rotor resistance per phase(in ohm) : \",round(R2,3);\n", + "#Note : Some answers are not accurate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.0333\n", + "Rotor Speed(in rpm) : 967.0\n", + "Total mechanial power developed(in KW) : 105.14\n", + "Rotor Copper Loss(in watts) : 2666.66666667\n", + "Copper Loss per phase(in watts) : 889.0\n", + "Rotor resistance per phase(in ohm) : 0.247\n" + ] + } + ], + "prompt_number": 148 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.23, Page 158" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "UsefulOutput=55;#in H.P.\n", + "MechLosses=2;#in H.P.\n", + "Stator_Losses=13.5;#in KW\n", + "\n", + "#Calculations&Results\n", + "MechPowerDev=UsefulOutput+MechLosses;#in H.P.\n", + "MechPowerDev=MechPowerDev*735.5/1000;#in KW\n", + "S=50.;#in %\n", + "Pin_Rotor=MechPowerDev/(1-S/100);#in KW\n", + "RotorCopperLoss=(S/100)*Pin_Rotor;#in KW\n", + "print \"Rotor Copper Loss(in KW) : \",round(RotorCopperLoss,2);\n", + "Pin_Motor=Pin_Rotor+Stator_Losses;#in KW\n", + "Efficiency=UsefulOutput*0.7355/Pin_Motor;#/unitless\n", + "print \"Effiiency(in %) :\",round((Efficiency*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor Copper Loss(in KW) : 41.92\n", + "Effiiency(in %) : 41.6\n" + ] + } + ], + "prompt_number": 149 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.24, Page 159" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "Pin_rotor=80;#Power input of rotor in KW\n", + "f=50;#in Hz\n", + "fdash=100;#alternations per minute\n", + "\n", + "#Calculations&Results\n", + "fdash=100./60;#in Hz\n", + "S=fdash/f;#unitless\n", + "print \"Slip : \",round(S,3);\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "N=Ns-Ns*S;#in rpm\n", + "print \"Rotor Speed(in rpm) : \",round(N);\n", + "RotorCopperLoss=S*Pin_rotor*10**3;#in watts\n", + "#print \"Rotor Copper Loss(in watts) : \",RotorCopperLoss;\n", + "CopperLoss=RotorCopperLoss/phase;#in watts/phase\n", + "print \"Copper Loss per phase(in watts) : \",round(CopperLoss);\n", + "I2=60;#in Ampere\n", + "R2=CopperLoss/I2**2;#in ohm\n", + "print \"Rotorresistance per phase(in ohm) : \",round(R2,3);\n", + "MechPowerDev=Pin_rotor*(1-S)*10**3/735.5;#in H.P.\n", + "print \"Total mechanial power developed(in H.P.) : \",round(MechPowerDev,1);\n", + "#Note : Some answers are not accurate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.033\n", + "Rotor Speed(in rpm) : 967.0\n", + "Copper Loss per phase(in watts) : 889.0\n", + "Rotorresistance per phase(in ohm) : 0.247\n", + "Total mechanial power developed(in H.P.) : 105.1\n" + ] + } + ], + "prompt_number": 150 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.25, Page 160" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Efficiency=0.9;#unitless\n", + "Output=50;#in H.P.\n", + "\n", + "#Calculations\n", + "#formula : Efficiency=Output/(Output+Losses) \n", + "Losses=((1-Efficiency)/Efficiency)*Output*735.5;#in watts\n", + "Losses=round(Losses);#round\n", + "#Let, Ststor Cu loss = Rotor Cu los = Iron loss=K\n", + "#Mechanical Loss = Iron Loss/3 = K/3\n", + "#TotalLosses=k+K+K+K/3\n", + "K=Losses*3/10;#in watts\n", + "Pin_rotor=Output*735.5+Losses;#in watts\n", + "Slip=K/Pin_rotor;#unitless\n", + "\n", + "#Results\n", + "print \"Slip :\",round(Slip,3);\n", + "print \"or \",round((Slip*100)),\"%\"; #answer differes due to rouding" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.03\n", + "or 3.0 %\n" + ] + } + ], + "prompt_number": 154 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.26, Page 160" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Pin_rotor=20;#Power input of rotor in KW\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "f=50.;#in Hz\n", + "N=960;#in rpm(Actual speed of motor)\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "S=(Ns-N)/Ns;#unitless\n", + "RotorCuLoss=S*Pin_rotor*10**3;#in watts\n", + "RotorCuLoss=RotorCuLoss/phase;#in watts/phase\n", + "R2=1./3;#in ohm(Rotor resistance per phase)\n", + "I2=sqrt(RotorCuLoss/R2);#in Ampere\n", + "\n", + "#Result\n", + "print \"Motor current per phase(in Ampere) :\",round(I2,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Motor current per phase(in Ampere) : 28.3\n" + ] + } + ], + "prompt_number": 155 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.27, Page 161" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "V=500;#in volt\n", + "Pout_rotor=20;#Power output of rotor in H.P.\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "f=50;#in Hz\n", + "N=995.;#in rpm(Actual speed of motor)\n", + "cosfi=0.87;#powerfactor(unitless)\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "S=(Ns-N)/Ns;#unitless\n", + "print \"Slip : \",S;\n", + "RotorCuLoss=(S/(1-S))*Pout_rotor*735.5;#in watts\n", + "print \"Rotor Cu Loss(in watts) : \",round(RotorCuLoss,2);\n", + "Pin_rotor=RotorCuLoss/S;#in watts\n", + "print \"Power input to rotor(in KW) :\",round((Pin_rotor/10**3),2);\n", + "LineCurrent=Pin_rotor/(sqrt(3)*V*cosfi);#in Ampere\n", + "print \"Line Current(in A) :\",round(LineCurrent,2);\n", + "RotorFreq=S*f;#in Hz\n", + "print \"Rotor Frequency(in Hz) :\",RotorFreq;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.005\n", + "Rotor Cu Loss(in watts) : 73.92\n", + "Power input to rotor(in KW) : 14.78\n", + "Line Current(in A) : 19.62\n", + "Rotor Frequency(in Hz) : 0.25\n" + ] + } + ], + "prompt_number": 156 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.28, Page 162" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3;#no. of phase\n", + "Efficiency=90;#in %\n", + "Speed=480;#in rpm\n", + "VL=400;#in volt\n", + "IL=75;#in Ampere\n", + "cosfi=0.77;#powerfactor(unitless)\n", + "d=0.75;#diameter of pulley in meter\n", + "\n", + "#Calculations\n", + "Pin_motor=sqrt(3)*VL*IL*cosfi;#Power input of motor in watts\n", + "OutputPower=Pin_motor*Efficiency/100;#in watts\n", + "Omega=Speed*2*pi/60;#angular speed in radians/sec\n", + "Torque=OutputPower/Omega;#in N-meter\n", + "Torque=Torque/9.81;#in Kg-meter\n", + "PullOnBelt=Torque/(d/2);#in Kg\n", + "\n", + "#Result\n", + "print \"Pull On Belt(in Kg.) : \",round(PullOnBelt,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pull On Belt(in Kg.) : 194.74\n" + ] + } + ], + "prompt_number": 157 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.29, Page 162" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "#At 3% slip\n", + "OutputPower=24;#in KW(At 3% slip)\n", + "S=3.;#in %\n", + "\n", + "#Calculations\n", + "Efficiency=(1-S/100);#unitless\n", + "InputPower=OutputPower/Efficiency;#in KW\n", + "#At 5% slip\n", + "S=5;#in %\n", + "Efficiency=(1-S/100);#unitless\n", + "OutputPower=InputPower*Efficiency;#in KW\n", + "\n", + "#Result\n", + "print \"Mechanical power output at a slip of 5%(in KW) : \",round(OutputPower,1);\n", + "#answer differs due to rounding" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Mechanical power output at a slip of 5%(in KW) : 24.7\n" + ] + } + ], + "prompt_number": 158 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.30, Page 162" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "#At 50% speed of synchronous speed\n", + "S=50.;#in %\n", + "UsefulOutput=55;#in H.P.\n", + "MechLoss=2;#in H.P.\n", + "StatorLoss=13.5;#in KW\n", + "\n", + "#Calculations&Results\n", + "Pout_rotor=UsefulOutput+MechLoss;#in H.P.\n", + "Pout_rotor=Pout_rotor*0.7355;#in KW\n", + "Efficiency=(1-S/100);#unitless\n", + "RotorInputPower=Pout_rotor/Efficiency;#in KW\n", + "RotorCuLoss=RotorInputPower-Pout_rotor;#in KW\n", + "print \"Cu Loss in the rotor circuit(in KW) :\",round(RotorCuLoss,2);\n", + "TotalLosses=StatorLoss+RotorCuLoss+MechLoss*0.7355;#in KW\n", + "MotorEfficiency=UsefulOutput*0.7355/(UsefulOutput*0.7355+TotalLosses);\n", + "print \"Motor Efficiency(in %) : \",round((MotorEfficiency*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Cu Loss in the rotor circuit(in KW) : 41.92\n", + "Motor Efficiency(in %) : 41.6\n" + ] + } + ], + "prompt_number": 159 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.31, Page 163" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "N=1440.;#in rpm(Actual speed of motor)\n", + "Power_dev=10;#Power developed in H.P.\n", + "VL=400;#in volt\n", + "cosfi=0.8;#powerfactor(unitless)\n", + "f=50;#in Hz\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#fractional slip(unitless)\n", + "print \"Slip :\",S;\n", + "Omega=N*2*pi/60;#angular speed in radians/sec\n", + "Torque=Power_dev*735.5/Omega;#in N-meter\n", + "Torque=Torque/9.81;#in Kg-meter\n", + "print \"Torque(in Kg-meter) :\",round(Torque,2);\n", + "RotorCuLoss=(S/(1-S))*Power_dev*735.5;#in watts\n", + "print \"Rotor Cu Loss per phase(in watts) : \",round((RotorCuLoss/3),1);\n", + "Pin_rotor=RotorCuLoss/S;#in watts\n", + "print \"Power input to rotor(in KW) :\",round((Pin_rotor/10**3),2);\n", + "LineCurrent=Pin_rotor/(sqrt(3)*VL*cosfi);#in Ampere\n", + "print \"Line Current(in A) :\",round(LineCurrent,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.04\n", + "Torque(in Kg-meter) : 4.97\n", + "Rotor Cu Loss per phase(in watts) : 102.2\n", + "Power input to rotor(in KW) : 7.66\n", + "Line Current(in A) : 13.82\n" + ] + } + ], + "prompt_number": 160 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.32, Page 164" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=440;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin_rotor=80;#in KW\n", + "fr=100;#revolutions/min\n", + "\n", + "#Calculations&Results\n", + "fr=100./60;#in Hz\n", + "S=fr/f;#slip(unitless)\n", + "print \"Slip : \",round(S,3);\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-S);#in rpm(Speed of Rotor)\n", + "print \"Rotor speed(in rpm) : \",round(N);\n", + "RotorCuLoss=S*Pin_rotor*10**3;#in watts\n", + "print \"Rotor Cu Loss per phase(in watts) : \",round(RotorCuLoss/phase);\n", + "I2=65;#in Ampere\n", + "R2=(RotorCuLoss/phase)/I2**2;#in ohm\n", + "print \"Rotor resistance per phase(in ohm) : \",round(R2,2);\n", + "Pout_rotor=Pin_rotor-RotorCuLoss/10**3;#in KW\n", + "print \"Output power of rotor(in H.P.) : \",round((Pout_rotor/0.735),1);\n", + "#Note : answers of few part are not accurate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.033\n", + "Rotor speed(in rpm) : 967.0\n", + "Rotor Cu Loss per phase(in watts) : 889.0\n", + "Rotor resistance per phase(in ohm) : 0.21\n", + "Output power of rotor(in H.P.) : 105.2\n" + ] + } + ], + "prompt_number": 161 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.33, Page 165" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=440;#in volt\n", + "f=50.;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin_rotor=20;#in KW\n", + "N=960.;#in rpm(Speed of Motor)\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#slip(unitless)\n", + "RotorCuLoss=S*Pin_rotor*10**3;#in watts\n", + "RotorCuLoss=RotorCuLoss/phase;#in watts/phase\n", + "R2=1./3;#Rotor Resistance in ohm per phase\n", + "I2=sqrt(RotorCuLoss/R2)\n", + "\n", + "#Result\n", + "print \"Motor current per phase(in Ampere) :\",round(I2,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Motor current per phase(in Ampere) : 28.28\n" + ] + } + ], + "prompt_number": 162 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.34, Page 166" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "VL=500;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pout_rotor=20;#in H.P.\n", + "cosfi=0.87;#power factor\n", + "N=995.;#in rpm(Speed of Motor)\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#slip(unitless)\n", + "print \"Slip : \",S;\n", + "RotorCuLoss=(S/(1-S))*Pout_rotor*735.5;#in watts\n", + "print \"Rotor Cu Loss(in watts) :\",round(RotorCuLoss,2);\n", + "Pin_rotor=RotorCuLoss/S;#in watts\n", + "print \"Power input to rotor(in KW) :\",round((Pin_rotor/10**3),3);\n", + "IL=Pin_rotor/(sqrt(3)*VL*cosfi);#in Ampere\n", + "print \"Line current (in Ampere) :\",round(IL,2);\n", + "Rotorfreq=S*f;#in Hz\n", + "print \"Rotor frequency(in Hz) :\",round(Rotorfreq,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.005\n", + "Rotor Cu Loss(in watts) : 73.92\n", + "Power input to rotor(in KW) : 14.784\n", + "Line current (in Ampere) : 19.62\n", + "Rotor frequency(in Hz) : 0.25\n" + ] + } + ], + "prompt_number": 163 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.35, Page 166" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "Efficiency=85.;#in %\n", + "OutputPower=17.;#in KW\n", + "OutputPower=17*10**3;#in watts\n", + "StatorLosses=900.;#in watt\n", + "MechLosses=1100;#in watt\n", + "\n", + "#Calculations&Results\n", + "Pin=OutputPower/(Efficiency/100);#in watts\n", + "#Stator input or rotor input :\n", + "Pin_rotor=Pin-StatorLosses;#in watts\n", + "P_MechDev=OutputPower+MechLosses;#in watts\n", + "RotorCuLosses=Pin_rotor-P_MechDev;#in watts\n", + "S=RotorCuLosses/Pin_rotor;#slip(unitless)\n", + "print \"Slip :\",round(S,3);\n", + "Ns=120*f/P\n", + "N=Ns*(1-S);#in rpm(Speed of Motor)\n", + "Omega=2*pi*N/60;#angular speed in rad/sec\n", + "Torque=P_MechDev/Omega;#in Nm\n", + "print \"Torque developed(in Nm) :\",round(Torque,1);\n", + "Tm=OutputPower/Omega;#in Nm\n", + "Tm=Tm/9.81;#in Kg-meter\n", + "print \"Torque at the shaft(in Kg-meter) :\",round(Tm,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.052\n", + "Torque developed(in Nm) : 121.6\n", + "Torque at the shaft(in Kg-meter) : 11.64\n" + ] + } + ], + "prompt_number": 164 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.36, Page 168" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=500;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin_stator=50;#in KW\n", + "Statorlosses=2;#in KW\n", + "MechLosses=1;#in KW\n", + "\n", + "#Calculations&Results\n", + "Pin_rotor=Pin_stator-Statorlosses;#in KW\n", + "S=3.;#in %\n", + "RotorI2RLoss=(S/100)*Pin_rotor;#in KW\n", + "print \"Rotor I2R Loss(in KW) :\",RotorI2RLoss;\n", + "print \"Rotor core loss at 3% slip is very less and can be neglected.\";\n", + "Power_dev=Pin_rotor-RotorI2RLoss;#in KW\n", + "print \"Power developed by the rotor(in H.P.) : \",round((Power_dev/0.7355),1);\n", + "Outputpower=Power_dev-MechLosses;#in KW\n", + "Efficiency=Outputpower/Pin_stator;#unitless\n", + "print \"Efficiency of the motor(in %) :\",(Efficiency*100);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor I2R Loss(in KW) : 1.44\n", + "Rotor core loss at 3% slip is very less and can be neglected.\n", + "Power developed by the rotor(in H.P.) : 63.3\n", + "Efficiency of the motor(in %) : 91.12\n" + ] + } + ], + "prompt_number": 165 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.37, Page 168" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "MotorOutput=20;#in H.P.\n", + "MotorOutput=20*735.5;#in watts\n", + "S=4.;#full load slip in %\n", + "MechLosses=500;#in watts\n", + "\n", + "#Calculations&Results\n", + "Pdev_rotor=MotorOutput+MechLosses;#in watts\n", + "S=S/100;#fractional slip\n", + "RotorI2RLoss=(S/(1-S))*Pdev_rotor;#in watts\n", + "print \"Rotor I2R Loss(in watts) :\",RotorI2RLoss;\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "print \"Rotor speed(in rpm) :\",Nr;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor I2R Loss(in watts) : 633.75\n", + "Rotor speed(in rpm) : 1440.0\n" + ] + } + ], + "prompt_number": 80 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.38, Page 169" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "R2=0.02;#in ohm\n", + "X2=0.1;#in ohm\n", + "\n", + "#Calculations\n", + "S=R2/X2;#slip(unitless)\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "#At starting S=1\n", + "S=1;#slip\n", + "#Formula : T=K*S*R2/(R2^2+X2^2)\n", + "#Starting torque, Tst=K*R2/(R2^2+S^2*X2^2)\n", + "#Maximum torque, Tm=K/(2*X2)\n", + "#Tst=(2/3)*Tm : gives a equation\n", + "#100*R2^2-30*R2+1=0\n", + "a=100;b=-30;c=1;\n", + "R21=(-b+sqrt(b**2-4*a*c))/(2*a);#in ohm\n", + "R22=(-b-sqrt(b**2-4*a*c))/(2*a);#in ohm\n", + "#This R2 is the value of rotor circuit resistance.\n", + "RotWinResistance=0.02;#in ohm per phase\n", + "Extra_R1=R21-RotWinResistance;#in ohm\n", + "Extra_R2=R22-RotWinResistance;#in ohm\n", + "\n", + "#Result\n", + "print \"Extra Resistance(in ohm) : \",round(Extra_R1,3),\" ohm or \",round((Extra_R2),3),\" ohm.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Extra Resistance(in ohm) : 0.242 ohm or 0.018 ohm.\n" + ] + } + ], + "prompt_number": 166 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.39, Page 170" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "MotorShaftOutput=20;#in H.P.\n", + "\n", + "#Calculations&Results\n", + "MotorShaftOutput=20*735.5;#in watts\n", + "MechLosses=MotorShaftOutput*2/100;#in watts(2% of the output)\n", + "Pdev_rotor=MotorShaftOutput+MechLosses;#in watts\n", + "S=3.;#slip in %\n", + "S=S/100;#fractional slip\n", + "RotorI2RLoss=(S/(1-S))*Pdev_rotor;#in watts\n", + "print \"Rotor I2R Loss(in watts) :\",round(RotorI2RLoss);\n", + "print \"Rotor iron loss at 3% slip is very small and can be neglected.\";\n", + "RotorInput=Pdev_rotor+RotorI2RLoss;#in watts\n", + "print \"Rotor Input(in Watts) :\",round(RotorInput,1);\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "Nr=Nr/60;#in rps\n", + "T=MotorShaftOutput/(2*pi*Nr);#in Nm\n", + "print \"Output torque(in Nm) : \",round(T,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor I2R Loss(in watts) : 464.0\n", + "Rotor iron loss at 3% slip is very small and can be neglected.\n", + "Rotor Input(in Watts) : 15468.2\n", + "Output torque(in Nm) : 96.54\n" + ] + } + ], + "prompt_number": 167 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.40, Page 171" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=500;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Nr=975.;#in rpm\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-Nr)/Ns;#slip\n", + "print \"Slip :\",S;\n", + "Pin_stator=40;#in KW\n", + "StatorLosses=1;#in KW\n", + "Pin_rotor=Pin_stator-StatorLosses;#in KW\n", + "RotorCuLosses=S*Pin_rotor;#in KW\n", + "print \"Rotor Cu Losses(in KW) :\",RotorCuLosses;\n", + "RotorOutput=Pin_rotor-RotorCuLosses;#in KW\n", + "OutputHP=RotorOutput/0.735;#in H.P.\n", + "print \"Output Horse Power : \",round(OutputHP,2);\n", + "Efficiency=RotorOutput/Pin_stator;#unitless\n", + "print \"Efficiency(in %) :\",round(Efficiency*100);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.025\n", + "Rotor Cu Losses(in KW) : 0.975\n", + "Output Horse Power : 51.73\n", + "Efficiency(in %) : 95.0\n" + ] + } + ], + "prompt_number": 168 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.41, Page 172" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=440;#in volt\n", + "f=50.;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3.;#no. of phase\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "fr=120;#alternations per minute\n", + "fr=fr/60;#in Hz\n", + "S=fr/f;#slip\n", + "print \"Slip : \",S;\n", + "Nr=Ns-S*Ns;#in rpm\n", + "print \"Rotor speed(in rpm) :\",Nr;\n", + "Rotor_input=80;#in KW\n", + "RotorCuLoss=S*Rotor_input;#in KW\n", + "print \"Rotor Cu Loss per phase(in watts) :\",round((RotorCuLoss*10**3)/phase);\n", + "P_Mechdev=Rotor_input*10**3-RotorCuLoss*10**3;#in watts\n", + "P_Mechdev=P_Mechdev/735.5;#in H.P.\n", + "print \"Mechanical power devloped(in H.P.) :\",round(P_Mechdev,1);\n", + "Ir=60;#in Ampere\n", + "R2=(RotorCuLoss*10**3/phase)/Ir**2;#in ohm\n", + "print \"Rotor resistance per phase(in ohm) :\",round(R2,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.04\n", + "Rotor speed(in rpm) : 960.0\n", + "Rotor Cu Loss per phase(in watts) : 1067.0\n", + "Mechanical power devloped(in H.P.) : 104.4\n", + "Rotor resistance per phase(in ohm) : 0.296\n" + ] + } + ], + "prompt_number": 169 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.43, Page 175 " + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50.;#in Hz\n", + "P=6;#no. of poles\n", + "R2=0.2;#rotor resistance per phase in ohm\n", + "Nr=960;#in rpm\n", + "Nr_dash=600;#in rpm\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-Nr)/Ns;#slip(unitless)\n", + "Sdash=(Ns-Nr_dash)/Ns;#slip(unitless)\n", + "#Let the new value of resistance is R2dash=(R+0.2)\n", + "R=R2*(Sdash/S)-R2;#Resistance to be added in ohm\n", + "\n", + "#Result\n", + "print \"Resistance to be added(in ohm) : \",R;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Resistance to be added(in ohm) : 1.8\n" + ] + } + ], + "prompt_number": 91 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.44, Page 176" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "RotorIntResistance=0.1;#in ohm per phase\n", + "X2=1;#in ohm/phase\n", + "\n", + "#Calculation\n", + "Rext=X2-RotorIntResistance;#in ohm\n", + "\n", + "#Result\n", + "print \"External resistance to be included(in ohm/phase) : \",Rext;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "External resistance to be included(in ohm/phase) : 0.9\n" + ] + } + ], + "prompt_number": 92 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.45, Page 176" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50.;#in Hz\n", + "P=6.;#no. of poles\n", + "phase=3.;#no. of phase\n", + "R2=0.2;#rotor resistance per phase in ohm\n", + "N1=960.;#Full load speed in rpm\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "S1=(Ns-N1)/Ns;#Full load slip(unitless)\n", + "N2=N1*(1-10./100);#New speed in rpm(reduced 10%)\n", + "S2=(Ns-N2)/Ns;#New slip(unitless)\n", + "#Formula : S=RotorCuLoss/Pin_rotor=3*I2^2*R2/Pin_rotor\n", + "#Let the additional resistance is R\n", + "R=R2*S2/S1-R2;#Resistance to be added in ohm\n", + "\n", + "#Result\n", + "print \"Additional Rotor Resistance(in ohm) : \",R;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Additional Rotor Resistance(in ohm) : 0.48\n" + ] + } + ], + "prompt_number": 97 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.46, Page 180" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "R2inner=0.4;#in ohm\n", + "X2inner=2;#in ohm \n", + "R2outer=2;#in ohm\n", + "X2outer=0.4;#in ohm \n", + "#At standstill :\n", + "S=1;#unitless\n", + "\n", + "#Calculations&Results\n", + "Ratio=(R2outer*S/(R2outer**2+S**2*X2outer**2))/(R2inner*S/(R2inner**2+S**2*X2inner**2));#unitless\n", + "print \"Ratio of torque produced by two cages at standstill : \",Ratio;\n", + "#At Full load :\n", + "S=5;#in %\n", + "S=5./100;#fractional\n", + "Ratio=(R2outer*S/(R2outer**2+S**2*X2outer**2))/(R2inner*S/(R2inner**2+S**2*X2inner**2));#unitless\n", + "print \"Ratio of torque produced by two cages at full load(slip=5%) : \",round(Ratio,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ratio of torque produced by two cages at standstill : 5.0\n", + "Ratio of torque produced by two cages at full load(slip=5%) : 0.21\n" + ] + } + ], + "prompt_number": 171 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.47, Page 187" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "Output=10;#in H.P.\n", + "Output=Output*735.5;#in watts\n", + "cosfi=0.8;#unitless\n", + "ETA=0.83;#unitless\n", + "ISCbyIFL=3.5;#ratio of SC current to full load current\n", + "VL=500;#in volt\n", + "\n", + "#Calculations\n", + "Input=Output/ETA;#in watts\n", + "IFL=Input/(sqrt(3)*VL*cosfi);#in Ampere\n", + "ISC=IFL*ISCbyIFL;#in Ampere\n", + "Is=ISC/3;#in Ampere\n", + "\n", + "#Result\n", + "print \"Strting current(in Ampere) :\",round(Is,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Strting current(in Ampere) : 14.92\n" + ] + } + ], + "prompt_number": 172 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.48, Page 187" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=5;#ratio of SC current to full load current\n", + "S=0.04;#Full load slip\n", + "\n", + "#Calculations&Results\n", + "#(i) for star delta starter : \n", + "print \"for star delta starter : \";\n", + "Ratio=(1./3)*(ISCbyIFL)**2*S;#Ratio of starting torque to full load torque\n", + "print \"Starting torque is \",round((Ratio*100),2),\"% of full load torque.\";\n", + "#(i) for auto transformer starter :\n", + "print \"for auto transformer starter : \";\n", + "K=50./100;#tappings\n", + "Ratio=K**2*(ISCbyIFL)**2*S;#Ratio of starting torque to full load torque\n", + "print \"Starting torque is \",(Ratio*100),\"% of full load torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "for star delta starter : \n", + "Starting torque is 33.33 % of full load torque.\n", + "for auto transformer starter : \n", + "Starting torque is 25.0 % of full load torque.\n" + ] + } + ], + "prompt_number": 173 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.49, Page 188" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=4;#ratio of SC current to full load current\n", + "TsBYTf=1./4;#ratio of atarting torque to full load torque\n", + "Slip=3;#in %\n", + "Slip=3./100;#in fraction\n", + "\n", + "#Calculations\n", + "#Formula : TsBYTf=Percent_Tapping^2*ISCbyIFL^2*Slip\n", + "tapping=sqrt(TsBYTf/(Slip*ISCbyIFL**2));#in %\n", + "\n", + "#Result\n", + "print \"Tapping(in %) :\",round((tapping*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Tapping(in %) : 72.2\n" + ] + } + ], + "prompt_number": 174 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.50, Page 189" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=5.;#ratio of SC current to full load current\n", + "Slip=5;#in %\n", + "Slip=5./100;#in fraction\n", + "\n", + "#Calculations\n", + "#Formula : 3*IFL=K^2*ISC\n", + "K=sqrt(1./(ISCbyIFL/3));#unitless\n", + "TsBYTfl=K**2*(ISCbyIFL)**2*Slip;#ratio of starting torque to full load torque\n", + "\n", + "#Result\n", + "print \"Starting torque is \",(TsBYTfl*100),\"% of full load torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Starting torque is 75.0 % of full load torque.\n" + ] + } + ], + "prompt_number": 175 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.51, Page 189" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "print \"Star delta starter :\";\n", + "ISCbyIFL=6;#ratio of SC current to full load current\n", + "Slip=4;#in %\n", + "Slip=4./100;#in fraction\n", + "\n", + "#Calculations&Results\n", + "TsBYTfl=(1./3)*(ISCbyIFL)**2*Slip;#ratio of starting torque to full load torque\n", + "print \"Starting torque is \",(TsBYTfl*100),\"% of full load value.\";\n", + "print \"For an auto transformer :\";\n", + "K=70.70;#in %\n", + "K=70.70/100;#in fraction\n", + "TsBYTfl=K**2*(ISCbyIFL)**2*Slip;#ratio of starting torque to full load torque\n", + "print \"Starting torque is \",round(TsBYTfl*100),\"% of full load torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Star delta starter :\n", + "Starting torque is 48.0 % of full load value.\n", + "For an auto transformer :\n", + "Starting torque is 72.0 % of full load torque.\n" + ] + } + ], + "prompt_number": 176 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.52, Page 190" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "V=400.;#in volt\n", + "S=3.5;#in %\n", + "S=3.5/100;#in fraction\n", + "VL=92;#in volts\n", + "\n", + "#Calculations&Results\n", + "#ISC=(V/VL)*IFL;#in Ampere\n", + "ISCbyIFL=V/VL;#ratio of SC current to full load current\n", + "#2*IFL=K^2*ISC imples that 2*IFL=K^2*(V/VL)*IFL\n", + "K=sqrt(2./(V/VL));#in fraction\n", + "print \"Necessary tapping(in %) :\",round((K*100),1);\n", + "TsBYTfl=K**2*(ISCbyIFL)**2*S;#ratio of starting torque to full load torque\n", + "print \"Starting torque is \",round((TsBYTfl*100),1),\"% of full load value.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Necessary tapping(in %) : 67.8\n", + "Starting torque is 30.4 % of full load value.\n" + ] + } + ], + "prompt_number": 179 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.53, Page 190" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=4;#ratio of SC current to full load current\n", + "ISbyIFL=2;#ratio of Supply current to full load current\n", + "S=2.5;#in %\n", + "S=2.5/100;#in fraction\n", + "\n", + "#Calculations\n", + "#Formula : ISbyIFL=(1/K^2)*ISCbyIFL\n", + "K=sqrt(ISCbyIFL/ISbyIFL);#in fraction\n", + "print \"Auto Transormation ratio :\",round(K,2);\n", + "TsBYTfl=(1./K**2)*(ISCbyIFL)**2*S;#ratio of starting torque to full load torque\n", + "\n", + "#Results\n", + "print \"Ratio of starting torque to full load torque :\",TsBYTfl;\n", + "print \"or Starting torque is \",(TsBYTfl*100),\"% of full load value.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Auto Transormation ratio : 1.41\n", + "Ratio of starting torque to full load torque : 0.2\n", + "or Starting torque is 20.0 % of full load value.\n" + ] + } + ], + "prompt_number": 180 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.54, Page 191" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=4;#ratio of SC current to full load current\n", + "S=3;#in %\n", + "S=3./100;#in fraction\n", + "TsBYTfl=1./4;#ratio of starting torque to full load torque\n", + "\n", + "#Calculations\n", + "#Formula : TsBYTfl=ISCbyIFL^2*K^2*S\n", + "K=sqrt(TsBYTfl/(ISCbyIFL**2*S));#in fraction\n", + "\n", + "#Result\n", + "print \"Auto Transormation ratio(in %) :\",round((K*100),2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Auto Transormation ratio(in %) : 72.17\n" + ] + } + ], + "prompt_number": 181 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.55, Page 191" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "Output=3;#in H.P.\n", + "Output=3*735.5;#in watts\n", + "Efficiency=0.83;#unitless\n", + "cosfi=0.8;#power factor\n", + "\n", + "#Calculations\n", + "Vl=500;#in volt\n", + "Input=Output/Efficiency;#in watts\n", + "#Formula : Input=sqrt(3)*Vl*Il*cosfi\n", + "Il=Input/(sqrt(3)*Vl*cosfi);#in Ampere\n", + "ISCbyIFL=3.5;#ratio of SC current to full load current\n", + "ISC=ISCbyIFL*Il;#in Ampere\n", + "LineCurrent=ISC/3;#in Ampere(for star delta starter)\n", + "\n", + "#Result\n", + "print \"Line Current(in Ampere) :\",round(LineCurrent,2);\n", + "#Note : Ans in the book is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Line Current(in Ampere) : 4.48\n" + ] + } + ], + "prompt_number": 182 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.56, Page 192" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "Output=15;#in H.P.\n", + "Output=15*735.5;#in watts\n", + "Vl=400;#in volt\n", + "ISCat200=40;#in Ampere(at 200 volt)\n", + "Efficiency=0.88;#unitless\n", + "cosfi=0.85;#power factor\n", + "\n", + "#Calculations\n", + "ISCat400=ISCat200*(400./200);#in Ampere(at 400 volt)\n", + "Input=Output/Efficiency;#in watts\n", + "Ifl=Input/(sqrt(3)*Vl*cosfi);#in Ampere\n", + "#starting line current with star delta starter\n", + "Is=ISCat400/3;#in Ampere\n", + "Ratio=Is/Ifl;#ratio of starting current to full load current\n", + "\n", + "#Result\n", + "print \"Ratio of line current at starting to full load current :\",round(Ratio,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ratio of line current at starting to full load current : 1.25\n" + ] + } + ], + "prompt_number": 183 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.57, Page 192" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "#With star delta starter : \n", + "TstBYTfl=0.35;#ratio of starting torque to full load torque\n", + "IstBYIfl=1.75;#ratio of starting current to full load current\n", + "\n", + "#Calculations&Results\n", + "ISCBYIs=sqrt(3);#ratio of SC current to starting current\n", + "ISCBYIfl=sqrt(3)*IstBYIfl;#ratio of SC current to full load current\n", + "#Formula : TstBYTfl=(ISCBYIfl)^2*S\n", + "S=TstBYTfl/(ISCBYIfl)**2;#in fraction\n", + "print \"Full load Slip : \",round(S,3);\n", + "#With auto transformer with winding in delta : \n", + "Ip=sqrt(3)*1.750*0.8;#full voltage phase current in Ampere\n", + "IlBYIf=4.2;#ratio of Line current to full load current\n", + "Ratio=IlBYIf**2*S;#ratio of starting current to full load current\n", + "print \"Ratio of line current at starting to full load current :\",round(Ratio,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Full load Slip : 0.038\n", + "Ratio of line current at starting to full load current : 0.67\n" + ] + } + ], + "prompt_number": 188 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.58, Page 196" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "R2=0.25;#in ohm per phase\n", + "S=2;#in %\n", + "S=2./100;#in fraction\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "#When speed reduced to 10%\n", + "NewSpeed=Nr*90/100;#in rpm\n", + "Sdash=(Ns-NewSpeed)/Ns;#in fraction\n", + "R2dash=(Sdash/S)*R2;#in ohm per phase\n", + "R=R2dash-R2;#in ohm\n", + "\n", + "#Result\n", + "print \"Value of resistance to be added(in ohm) :\",R;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Value of resistance to be added(in ohm) : 1.225\n" + ] + } + ], + "prompt_number": 114 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.59, Page 197" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "R2=0.25;#in ohm per phase\n", + "N=1440.;#in rpm at full load\n", + "NewSpeed=1200;#in rpm\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#in %\n", + "NewS=(Ns-NewSpeed)/Ns;#in fraction\n", + "S1=S;S2=NewS;#slip new and old\n", + "#Torque remaining same : S1/R2=S2/(R2+R)\n", + "R=S2*R2/S1-R2;#in ohm\n", + "print \"External resistance per phase(in ohm) :\",R; #incorrect answer in textbook\n", + "Nr=Ns-Ns*S;#in rpm\n", + "#If S1 is taken as 0.03 and S2 as 0.127\n", + "S1=0.03;#slip in fraction\n", + "S2=0.127;#slip in fraction\n", + "R=S2*R2/S1-R2;#in ohm\n", + "print \"External resistance per phase(in ohm) :\",round(R,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "External resistance per phase(in ohm) : -0.25\n", + "External resistance per phase(in ohm) : 0.808\n" + ] + } + ], + "prompt_number": 194 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.60, Page 198" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "S=4;#in %\n", + "S=4./100;#in fraction\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "N1=Ns-Ns*S;#in rpm\n", + "#When speed reduced to 10%\n", + "N2=N1*85/100;#in rpm(NewSpeed)\n", + "print \"New speed(in rpm) :\",N2;\n", + "#New speed is reduced by 15 %\n", + "Sdash=(Ns-N2)/Ns;#in fraction\n", + "print \"New Slip : \",Sdash;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "New speed(in rpm) : 1224.0\n", + "New Slip : 0.184\n" + ] + } + ], + "prompt_number": 118 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.61, Page 202" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin=60;#in KW\n", + "StatorLoss=1.2;#in KW\n", + "WindingLoss=1.8;#in KW\n", + "S=4;#in %\n", + "S=4./100;#in fraction\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "print \"Rotor speed(in rpm) :\",Nr;\n", + "Pin_rotor=Pin-StatorLoss;#in KW\n", + "RotorCuLoss=S*Pin_rotor;#in KW\n", + "print \"Rotor Cu Loss(in KW) : \",RotorCuLoss;\n", + "Pout_rotor=Pin_rotor-WindingLoss-RotorCuLoss;#in KW\n", + "Efficiency=(Pout_rotor/Pin)*100;#in %\n", + "print \"Efficiency(in %) :\",round(Efficiency);\n", + "HP=Pout_rotor*10**3/735.5;#in H.P.\n", + "print \"Horse Power(inH.P.) :\",round(HP,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor speed(in rpm) : 1440.0\n", + "Rotor Cu Loss(in KW) : 2.352\n", + "Efficiency(in %) : 91.0\n", + "Horse Power(inH.P.) : 74.3\n" + ] + } + ], + "prompt_number": 195 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/chapter_2_1.ipynb b/Electrical_Machines_II/chapter_2_1.ipynb new file mode 100755 index 00000000..c18790d9 --- /dev/null +++ b/Electrical_Machines_II/chapter_2_1.ipynb @@ -0,0 +1,2810 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:e3e9de8b11024daf372011e6906cf967e4f06ac156a70592b24845e8bc2b4696" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 2: Induction Motors" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.1, Page 130" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=2;#no. of poles\n", + "f=60;#in Hz\n", + "N=3460.;#in rpm\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "Slip=Ns-N;#in rpm\n", + "PercentageSlip=((Ns-N)/Ns)*100;#in %\n", + "\n", + "#Results\n", + "print \"Synchronous speed in rpm : \",Ns;\n", + "print \"Percentage Slip : \",round(PercentageSlip,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Synchronous speed in rpm : 3600\n", + "Percentage Slip : 3.89\n" + ] + } + ], + "prompt_number": 124 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.2, Page 131" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=8.;#no. of poles\n", + "f=50.;#in Hz\n", + "N=700.;#in rpm\n", + "\n", + "#Calculations&Results\n", + "#Part (a) :\n", + "Ns=120*f/P;#in rpm\n", + "print \"Synchronous speed in rpm : \",Ns;\n", + "#Part (b) :\n", + "S=(Ns-N)*100/Ns;#in %\n", + "print \"Slip(in %) : \",round(S,2);\n", + "#Part (c) :\n", + "#At the time of stsrt S=1;\n", + "fdash=f;#in Hz\n", + "print \"Rotor frequeny at the time of starting(in Hz) : \",fdash;\n", + "#Part (d) :\n", + "fdash=(S/100)*f;#in Hz\n", + "print \"Rotor frequeny at the given speed(in Hz) : \",round(fdash,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Synchronous speed in rpm : 750.0\n", + "Slip(in %) : 6.67\n", + "Rotor frequeny at the time of starting(in Hz) : 50.0\n", + "Rotor frequeny at the given speed(in Hz) : 3.3\n" + ] + } + ], + "prompt_number": 125 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.3, Page 131" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=10;#no. of poles\n", + "f=50;#in Hz\n", + "N=600;#in rpm\n", + "Nr=1440;#in rpm\n", + "\n", + "#Calculations&Results\n", + "f=P*N/120;#in Hz\n", + "#When P=2\n", + "P=2;#no. of poles\n", + "Ns=120*f/P;#in rpm\n", + "#When P=4\n", + "P=4;#no. of poles\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-Nr)*100/Ns;#Slip in %\n", + "print \"Percentage Slip(in %) : \",(S);\n", + "print \"No. of poles : \",(S);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Percentage Slip(in %) : 4\n", + "No. of poles : 4\n" + ] + } + ], + "prompt_number": 126 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.4, Page 132" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=8;#no. of poles\n", + "f=50;#in Hz\n", + "fr=1.5;#in Hz\n", + "\n", + "#Calculations\n", + "#Formula : fr=S*f\n", + "S=fr/f;#slip(unitless)\n", + "Ns=120*f/P;#in rpm\n", + "#Formula : S=(Ns-Nr)/Ns\n", + "Nr=Ns-S*Ns;#in rpm\n", + "\n", + "#Results\n", + "print \"Motor running speed in rpm : \",round(Nr);\n", + "print \"Slip(in %):\",(S*100);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Motor running speed in rpm : 728.0\n", + "Slip(in %): 3.0\n" + ] + } + ], + "prompt_number": 8 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.5, Page 133" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=2;#no. of poles\n", + "f=50;#in Hz\n", + "S=4.;#in %\n", + "\n", + "#Calculations&Results\n", + "#Part (i) : \n", + "Ns=f/P;#in rps\n", + "Ns=Ns*60;#in rpm\n", + "print \"Synchronous speed in rpm : \",Ns;\n", + "#Part (ii) :\n", + "print \"Slip = \",(S),\"% or \",(S/100);\n", + "#Part (iii) :\n", + "N=Ns*(1-S/100);#in rpm\n", + "print \"Actual speed in rpm :\",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Synchronous speed in rpm : 1500\n", + "Slip = 4.0 % or 0.04\n", + "Actual speed in rpm : 1440.0\n" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.6, Page 133" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "P=4;#no. of poles\n", + "f1=50;#in Hz\n", + "fdash=1.5;#in Hz\n", + "\n", + "#Calculations\n", + "S=fdash/f1;#unitless\n", + "print \"Slip : \",S;\n", + "Ns=120*f1/P;#in rpm\n", + "N=Ns*(1-S);#in rpm\n", + "\n", + "#Result\n", + "print \"Speed of motor in rpm :\",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.03\n", + "Speed of motor in rpm : 1455.0\n" + ] + } + ], + "prompt_number": 11 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.7, Page 138" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "N=1440;#in rpm\n", + "\n", + "#Calculations&Results\n", + "f=50;#in Hz\n", + "print \"For this speed of 1440 rpm the synchronous speed may be either 1500 rpm or 3000 rpm.\";\n", + "Ns=1500;#in rpm\n", + "P=4;#no. of poles\n", + "print \"No. of poles : \",P;\n", + "Slip=(Ns-N)*100/Ns;#in %\n", + "print Slip,\"Slip(in %) : \",Slip;\n", + "speed1=Ns-N;#Speed of rotor flux with respect to rotor\n", + "print \"Speed of rotor flux with respect to rotor in rpm : \",speed1;\n", + "speed2=(speed1)+N;#Speed of rotor flux with respect to stator\n", + "print \"Speed of rotor flux with respect to stator in rpm :\",speed2;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "For this speed of 1440 rpm the synchronous speed may be either 1500 rpm or 3000 rpm.\n", + "No. of poles : 4\n", + "4 Slip(in %) : 4\n", + "Speed of rotor flux with respect to rotor in rpm : 60\n", + "Speed of rotor flux with respect to stator in rpm : 1500\n" + ] + } + ], + "prompt_number": 12 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.8, Page 140" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "E2=100;#in volt\n", + "R2=0.05;#in ohm\n", + "X2=0.1;#im ohm\n", + "\n", + "#Calculations&Results\n", + "E2perphase=E2/sqrt(3);#in volt\n", + "#part (a) : \n", + "S=0.04;#slip\n", + "I2=S*E2perphase/sqrt(R2**2+(S*X2)**2);#in Ampere\n", + "print \"At 4% slip, Rotor current in Ampere : \",round(I2,2);\n", + "fi2=degrees(acos(R2/sqrt(R2**2+(S*X2)**2)));#in degree\n", + "print \"At 4% slip, Phase angle between rotor voltage and rotor current in degree :\",round(fi2,2);\n", + "#part (b) : \n", + "S=1;#slip\n", + "I2=S*E2perphase/sqrt(R2**2+(S*X2)**2);#in Ampere\n", + "print \"At 100% slip, Rotor current in Ampere : \",round(I2,2);\n", + "fi2=degrees(acos(R2/sqrt(R2**2+(S*X2)**2)));#in degree\n", + "print \"At 100% slip, Phase angle between rotor voltage and rotor current in degree :\",round(fi2,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "At 4% slip, Rotor current in Ampere : 46.04\n", + "At 4% slip, Phase angle between rotor voltage and rotor current in degree : 4.57\n", + "At 100% slip, Rotor current in Ampere : 516.4\n", + "At 100% slip, Phase angle between rotor voltage and rotor current in degree : 63.43\n" + ] + } + ], + "prompt_number": 137 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.9, Page 141" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "E2=100.;#in volt\n", + "R2=0.4;#in ohm\n", + "X2=2.25;#im ohm\n", + "\n", + "#Calculations\n", + "E2perphase=E2/sqrt(3);#in volt\n", + "S=4.;#in %\n", + "E=(S/100)*E2perphase;#rotor induced emf at a slip=4% in volt\n", + "Z2=sqrt(R2**2+((S/100)*X2)**2);\n", + "I2=E/Z2;#in Ampere\n", + "\n", + "#Result\n", + "print \"Rotor current in Ampere : \",round(I2,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor current in Ampere : 5.63\n" + ] + } + ], + "prompt_number": 138 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.10, Page 141" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "f=50;#in Hz\n", + "R2=0.2;#in ohm\n", + "X2=2;#im ohm\n", + "E2=60;#in volt\n", + "\n", + "#Calculations&Results\n", + "E2perphase=E2/sqrt(3);#in volt\n", + "#case (i) : S=1 \n", + "S=1;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"At standstill : Rotor current in Ampere : \",round(I2,2);\n", + "#case (ii) : S=0.1 \n", + "S=0.1;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"for 10% slip : Rotor current in Ampere : \",round(I2,2);\n", + "#case (iii) : S=0.05 \n", + "S=0.05;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"for 5% slip : Rotor current in Ampere : \",round(I2,2);\n", + "#case (iv) : S=0.01 \n", + "S=0.01;#unitless\n", + "Z2=sqrt(R2**2+(S*X2)**2);\n", + "I2=E2perphase*S/Z2;#in Ampere\n", + "print \"for 1% slip : Rotor current in Ampere : \",round(I2,2);\n", + "#Note : Answer in the book is wrong for S=0.05 and S=0.01" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "At standstill : Rotor current in Ampere : 17.23\n", + "for 10% slip : Rotor current in Ampere : 12.25\n", + "for 5% slip : Rotor current in Ampere : 7.75\n", + "for 1% slip : Rotor current in Ampere : 1.72\n" + ] + } + ], + "prompt_number": 139 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.11, Page 142" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "R2=0.5;#in ohm\n", + "X=4.5;#im ohm\n", + "E=50;#line voltage in volt\n", + "\n", + "#Calculations&Results\n", + "Eperphase=E/sqrt(3);#in volt\n", + "#part (a) : Sliprings are short circuited\n", + "Z2=sqrt(R2**2+X**2);#in ohm\n", + "I2=Eperphase/Z2;#in Ampere\n", + "print \"Sliprings are short circuited, Rotor current in Ampere : \",round(I2,2);\n", + "cosfi2=R2/Z2;#unitless\n", + "print \"Power factor : \",round(cosfi2,2);\n", + "#part (b) : Sliprings are connected \n", + "R2=0.4;#in ohm\n", + "Z2=sqrt(R2**2+X**2);#in ohm\n", + "I2=Eperphase/Z2;#in Ampere\n", + "print \"Sliprings are short circuited, Rotor current in Ampere : \",round(I2,2);\n", + "cosfi2=R2/Z2;#unitless\n", + "print \"Power factor : \",round(cosfi2,2);\n", + "#Note : answer of part (a) I2 is not curate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Sliprings are short circuited, Rotor current in Ampere : 6.38\n", + "Power factor : 0.11\n", + "Sliprings are short circuited, Rotor current in Ampere : 6.39\n", + "Power factor : 0.09\n" + ] + } + ], + "prompt_number": 140 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.12, Page 146" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data : \n", + "R2=0.02;#in ohm\n", + "X2=0.1;#im ohm\n", + "\n", + "#Calculations&Results\n", + "#let external resistance per phase = r then R2=R2+r ohm\n", + "S=1;#slip at starting\n", + "print \"Since at start speed is zero and slip is, therefore, unity or R2=X2\";\n", + "r=X2-R2;#in ohm\n", + "print \"External resistance per phase added to the rotor circuit in ohms : \",r;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Since at start speed is zero and slip is, therefore, unity or R2=X2\n", + "External resistance per phase added to the rotor circuit in ohms : 0.08\n" + ] + } + ], + "prompt_number": 22 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.13, Page 147" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data : \n", + "P=4;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.03;#in ohm\n", + "X2=0.12;#im ohm\n", + "\n", + "#Calculations&Results\n", + "Smax=R2/X2;#unitless\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-Smax);#in rpm\n", + "print \"Speed corresponding to maximum torque(in rpm) :\",N;\n", + "#alfa=(R2+r)/X2\n", + "#Ratio of starting torque to max torque=75/100=2*alfa/(alfa^2+1)\n", + "#It gives : 3*alfa^2-8*alfa+3=0\n", + "#coefficients :\n", + "a=3;b=-8;c=3;\n", + "alfa1=(-b+sqrt(b**2-4*a*c))/(2*a);alfa2=(-b-sqrt(b**2-4*a*c))/(2*a);\n", + "print \"rejecting higher values, alfa = \",round(alfa2,4);\n", + "fdash=1.5;#in Hz\n", + "r=alfa2*X2-R2;#in ohm;\n", + "print \"External rotor resistance per phase(in ohm) : \",round(r,4);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed corresponding to maximum torque(in rpm) : 1125.0\n", + "rejecting higher values, alfa = 0.4514\n", + "External rotor resistance per phase(in ohm) : 0.0242\n" + ] + } + ], + "prompt_number": 141 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.14, Page 147" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=4;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.024;#in ohm\n", + "X2=0.6;#in ohm \n", + "\n", + "#Calculations&Results\n", + "Sm=R2/X2;#slip corresponding to max torque\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-Sm);#in rpm\n", + "print \"Speed corresponding to maximum torque in rpm : \",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed corresponding to maximum torque in rpm : 1440.0\n" + ] + } + ], + "prompt_number": 26 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.15, Page 148" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=6;#no. of poles\n", + "f=50;#in Hz\n", + "Sf=4.;#in %\n", + "R2=0.01;#in ohm\n", + "X2=0.05;#im \n", + "\n", + "#Calculations&Results\n", + "Ratio=((R2/X2)**2+(Sf/100)**2)/(2*R2*(Sf/100)/X2);#ratio of max torque to full load torque\n", + "print \"Maximum torque, Tmax=\",(Ratio),\"Tf\";\n", + "Sm=R2/X2;#slip corresponding to max torque\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-Sm);#in rpm\n", + "print \"Speed corresponding to maximum torque in rpm : \",N;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Maximum torque, Tmax= 2.6 Tf\n", + "Speed corresponding to maximum torque in rpm : 800.0\n" + ] + } + ], + "prompt_number": 28 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.16, Page 149" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=12;#no. of poles\n", + "f=50.;#in Hz\n", + "R2=0.03;#in ohm\n", + "X2=0.5;#im \n", + "\n", + "#Calculations&Results\n", + "S=R2/X2;#unitless\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-S);#in rpm\n", + "print \"Speed at maximum torque in rpm : \",N;\n", + "S=(Ns-495)/Ns;#slip at 495 rpm speed\n", + "Ratio=(2*R2*S/X2)/((R2/X2)**2+S**2);#ratio of max torque to full load torque\n", + "print \"Ratio of full load torque to max torque(in %) :\",round((Ratio*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed at maximum torque in rpm : 470.0\n", + "Ratio of full load torque to max torque(in %) : 32.4\n" + ] + } + ], + "prompt_number": 142 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.17, Page 149" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "P=12;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.5;#in ohm\n", + "N=475.;#in rpm\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#unitless\n", + "X2=R2/S;#in ohm\n", + "#\"At maximum torque, Tmax is proportional to 1/2*X2 or 1/20\";\n", + "# \"Ts is proportional to (R2+R)/((R2+R)^2+X2^2)\";\n", + "# \"Also, Ts is proportional to 1(20*3)\";\n", + "# \"Equating the two eqn we have : (0.5*R)/(R2+R)^2+X2^2\";\n", + "#R^2-59R+70.25=0 : \n", + "#coefficients :\n", + "a=1;b=-59;c=70.25;\n", + "R=(-b-sqrt(b**2-4*a*c))/(2*a);\n", + "print \"Resistance R to be inserted in series (in ohm) :\",round(R,1);\n", + "#Note : answer in the book is wrong." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Resistance R to be inserted in series (in ohm) : 1.2\n" + ] + } + ], + "prompt_number": 143 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.18, Page 150" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=4;#no. of poles\n", + "f=50;#in Hz\n", + "R2=0.04;#in ohm\n", + "N=1200.;#in rpm(speed at max torque)\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm(synchronous speed)\n", + "S=(Ns-N)/Ns;#unitless\n", + "X2=R2/S;#in ohm\n", + "#Starting torque is the torque devloped when S=1: Tst=K*R2/(R2^2+X2^2)\n", + "#Let say, m=R2/(R2^2+X2^2) then Tst=K*m\n", + "m=R2/(R2**2+X2**2);#assumed\n", + "\n", + "#Results\n", + "print \"Starting torque, Tst=\",round(m,2),\"k\";\n", + "print \"Maximum torque, Tm=K/\",(2*X2);\n", + "print \"Thus, Tst in terms of Tm can be expressed as : \";\n", + "print \"Tst/Tm=\",(0.96*0.4),\" or Tst=\",(0.96*0.4),\"Tm\";\n", + "print \"Therefore, staring torque is \",(0.96*0.4*100),\"% of maximum torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Starting torque, Tst= 0.96 k\n", + "Maximum torque, Tm=K/ 0.4\n", + "Thus, Tst in terms of Tm can be expressed as : \n", + "Tst/Tm= 0.384 or Tst= 0.384 Tm\n", + "Therefore, staring torque is 38.4 % of maximum torque.\n" + ] + } + ], + "prompt_number": 144 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.19, Page 155" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=6;#no. of poles\n", + "f=50;#in Hz\n", + "fr=3.;#in Hz\n", + "R2=0.5;#in ohm\n", + "N=475;#in rpm\n", + "RotorInput=111.9;#in KW\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=fr*100/f;#unitless\n", + "print \"% Slip : \",S;\n", + "N=Ns*(1-S/100);#in rpm\n", + "print \"Speed of motor in rpm : \",N;\n", + "RotorCopperLoss=RotorInput*S/100;#in KW\n", + "RotorCopperLoss=RotorCopperLoss/3;#in KW/Phase\n", + "print \"Rotor Copper Loss per phase(in KW) : \",RotorCopperLoss;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "% Slip : 6.0\n", + "Speed of motor in rpm : 940.0\n", + "Rotor Copper Loss per phase(in KW) : 2.238\n" + ] + } + ], + "prompt_number": 146 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.20, Page 156" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3.;#no. of phase\n", + "Pin_stator=60.;#Power input of stator in KW\n", + "TotalStatorLosses=1.;#in KW\n", + "\n", + "#Calculations&Results\n", + "Pin_rotor=Pin_stator-TotalStatorLosses;#Power input of rotor in KW\n", + "S=3.;#slip in %\n", + "RotorCopperLosses=(S/100)*Pin_rotor;#in KW\n", + "RotorCopperLosses=RotorCopperLosses/phase;#in KW per phase\n", + "print \"Rotor Copper Losses per phase(in watts) : \",(RotorCopperLosses*10**3);\n", + "TotalMechPowerDev=Pin_rotor*(1-S/100);#in KW\n", + "print \"Total mechanial power developed(in KW) : \",TotalMechPowerDev;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor Copper Losses per phase(in watts) : 590.0\n", + "Total mechanial power developed(in KW) : 57.23\n" + ] + } + ], + "prompt_number": 40 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.21, Page 156" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=6;#no. of poles\n", + "S=4.;#slip in %\n", + "Output=20;#in KW\n", + "\n", + "#Calculations&Results\n", + "FrictionalLoss=250;#in watts\n", + "f=50;#in Hz\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns-Ns*(S/100);#in rpm\n", + "print \"Speed of motor(in rpm) : \",N;\n", + "MechPowerDeveloped=Output*10**3+FrictionalLoss;#in Watts\n", + "S=S/100;#unitless\n", + "RotorCopperLoss=(S/(1-S))*MechPowerDeveloped;#in watts\n", + "print \"Rotor Copper Loss(in watts) : \",RotorCopperLoss;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Speed of motor(in rpm) : 960.0\n", + "Rotor Copper Loss(in watts) : 843.75\n" + ] + } + ], + "prompt_number": 147 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.22, Page 157" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "Pin_rotor=80;#Power input of rotor in KW\n", + "f=50;#in Hz\n", + "fdash=100;#alternations per minute\n", + "\n", + "#Calculations&Results\n", + "fdash=100./60;#in Hz\n", + "S=fdash/f;#unitless\n", + "print \"Slip : \",round(S,4);\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "N=Ns-Ns*S;#in rpm\n", + "print \"Rotor Speed(in rpm) : \",round(N);\n", + "MechPowerDev=Pin_rotor*(1-S)*10**3/735.5;#in H.P.\n", + "print \"Total mechanial power developed(in KW) : \",round(MechPowerDev,2);\n", + "RotorCopperLoss=S*Pin_rotor*10**3;#in watts\n", + "print \"Rotor Copper Loss(in watts) : \",RotorCopperLoss;\n", + "CopperLoss=RotorCopperLoss/phase;#in watts/phase\n", + "print \"Copper Loss per phase(in watts) : \",round(CopperLoss);\n", + "I2=60;#in Ampere\n", + "R2=CopperLoss/I2**2;#in ohm\n", + "print \"Rotor resistance per phase(in ohm) : \",round(R2,3);\n", + "#Note : Some answers are not accurate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.0333\n", + "Rotor Speed(in rpm) : 967.0\n", + "Total mechanial power developed(in KW) : 105.14\n", + "Rotor Copper Loss(in watts) : 2666.66666667\n", + "Copper Loss per phase(in watts) : 889.0\n", + "Rotor resistance per phase(in ohm) : 0.247\n" + ] + } + ], + "prompt_number": 148 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.23, Page 158" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "UsefulOutput=55;#in H.P.\n", + "MechLosses=2;#in H.P.\n", + "Stator_Losses=13.5;#in KW\n", + "\n", + "#Calculations&Results\n", + "MechPowerDev=UsefulOutput+MechLosses;#in H.P.\n", + "MechPowerDev=MechPowerDev*735.5/1000;#in KW\n", + "S=50.;#in %\n", + "Pin_Rotor=MechPowerDev/(1-S/100);#in KW\n", + "RotorCopperLoss=(S/100)*Pin_Rotor;#in KW\n", + "print \"Rotor Copper Loss(in KW) : \",round(RotorCopperLoss,2);\n", + "Pin_Motor=Pin_Rotor+Stator_Losses;#in KW\n", + "Efficiency=UsefulOutput*0.7355/Pin_Motor;#/unitless\n", + "print \"Effiiency(in %) :\",round((Efficiency*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor Copper Loss(in KW) : 41.92\n", + "Effiiency(in %) : 41.6\n" + ] + } + ], + "prompt_number": 149 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.24, Page 159" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "Pin_rotor=80;#Power input of rotor in KW\n", + "f=50;#in Hz\n", + "fdash=100;#alternations per minute\n", + "\n", + "#Calculations&Results\n", + "fdash=100./60;#in Hz\n", + "S=fdash/f;#unitless\n", + "print \"Slip : \",round(S,3);\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "N=Ns-Ns*S;#in rpm\n", + "print \"Rotor Speed(in rpm) : \",round(N);\n", + "RotorCopperLoss=S*Pin_rotor*10**3;#in watts\n", + "#print \"Rotor Copper Loss(in watts) : \",RotorCopperLoss;\n", + "CopperLoss=RotorCopperLoss/phase;#in watts/phase\n", + "print \"Copper Loss per phase(in watts) : \",round(CopperLoss);\n", + "I2=60;#in Ampere\n", + "R2=CopperLoss/I2**2;#in ohm\n", + "print \"Rotorresistance per phase(in ohm) : \",round(R2,3);\n", + "MechPowerDev=Pin_rotor*(1-S)*10**3/735.5;#in H.P.\n", + "print \"Total mechanial power developed(in H.P.) : \",round(MechPowerDev,1);\n", + "#Note : Some answers are not accurate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.033\n", + "Rotor Speed(in rpm) : 967.0\n", + "Copper Loss per phase(in watts) : 889.0\n", + "Rotorresistance per phase(in ohm) : 0.247\n", + "Total mechanial power developed(in H.P.) : 105.1\n" + ] + } + ], + "prompt_number": 150 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.25, Page 160" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Efficiency=0.9;#unitless\n", + "Output=50;#in H.P.\n", + "\n", + "#Calculations\n", + "#formula : Efficiency=Output/(Output+Losses) \n", + "Losses=((1-Efficiency)/Efficiency)*Output*735.5;#in watts\n", + "Losses=round(Losses);#round\n", + "#Let, Ststor Cu loss = Rotor Cu los = Iron loss=K\n", + "#Mechanical Loss = Iron Loss/3 = K/3\n", + "#TotalLosses=k+K+K+K/3\n", + "K=Losses*3/10;#in watts\n", + "Pin_rotor=Output*735.5+Losses;#in watts\n", + "Slip=K/Pin_rotor;#unitless\n", + "\n", + "#Results\n", + "print \"Slip :\",round(Slip,3);\n", + "print \"or \",round((Slip*100)),\"%\"; #answer differes due to rouding" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.03\n", + "or 3.0 %\n" + ] + } + ], + "prompt_number": 154 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.26, Page 160" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Pin_rotor=20;#Power input of rotor in KW\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "f=50.;#in Hz\n", + "N=960;#in rpm(Actual speed of motor)\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "S=(Ns-N)/Ns;#unitless\n", + "RotorCuLoss=S*Pin_rotor*10**3;#in watts\n", + "RotorCuLoss=RotorCuLoss/phase;#in watts/phase\n", + "R2=1./3;#in ohm(Rotor resistance per phase)\n", + "I2=sqrt(RotorCuLoss/R2);#in Ampere\n", + "\n", + "#Result\n", + "print \"Motor current per phase(in Ampere) :\",round(I2,1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Motor current per phase(in Ampere) : 28.3\n" + ] + } + ], + "prompt_number": 155 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.27, Page 161" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "V=500;#in volt\n", + "Pout_rotor=20;#Power output of rotor in H.P.\n", + "phase=3;#no. of phase\n", + "P=6;#no. of poles\n", + "f=50;#in Hz\n", + "N=995.;#in rpm(Actual speed of motor)\n", + "cosfi=0.87;#powerfactor(unitless)\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#synchronous speed in rpm\n", + "S=(Ns-N)/Ns;#unitless\n", + "print \"Slip : \",S;\n", + "RotorCuLoss=(S/(1-S))*Pout_rotor*735.5;#in watts\n", + "print \"Rotor Cu Loss(in watts) : \",round(RotorCuLoss,2);\n", + "Pin_rotor=RotorCuLoss/S;#in watts\n", + "print \"Power input to rotor(in KW) :\",round((Pin_rotor/10**3),2);\n", + "LineCurrent=Pin_rotor/(sqrt(3)*V*cosfi);#in Ampere\n", + "print \"Line Current(in A) :\",round(LineCurrent,2);\n", + "RotorFreq=S*f;#in Hz\n", + "print \"Rotor Frequency(in Hz) :\",RotorFreq;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.005\n", + "Rotor Cu Loss(in watts) : 73.92\n", + "Power input to rotor(in KW) : 14.78\n", + "Line Current(in A) : 19.62\n", + "Rotor Frequency(in Hz) : 0.25\n" + ] + } + ], + "prompt_number": 156 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.28, Page 162" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "phase=3;#no. of phase\n", + "Efficiency=90;#in %\n", + "Speed=480;#in rpm\n", + "VL=400;#in volt\n", + "IL=75;#in Ampere\n", + "cosfi=0.77;#powerfactor(unitless)\n", + "d=0.75;#diameter of pulley in meter\n", + "\n", + "#Calculations\n", + "Pin_motor=sqrt(3)*VL*IL*cosfi;#Power input of motor in watts\n", + "OutputPower=Pin_motor*Efficiency/100;#in watts\n", + "Omega=Speed*2*pi/60;#angular speed in radians/sec\n", + "Torque=OutputPower/Omega;#in N-meter\n", + "Torque=Torque/9.81;#in Kg-meter\n", + "PullOnBelt=Torque/(d/2);#in Kg\n", + "\n", + "#Result\n", + "print \"Pull On Belt(in Kg.) : \",round(PullOnBelt,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pull On Belt(in Kg.) : 194.74\n" + ] + } + ], + "prompt_number": 157 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.29, Page 162" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "#At 3% slip\n", + "OutputPower=24;#in KW(At 3% slip)\n", + "S=3.;#in %\n", + "\n", + "#Calculations\n", + "Efficiency=(1-S/100);#unitless\n", + "InputPower=OutputPower/Efficiency;#in KW\n", + "#At 5% slip\n", + "S=5;#in %\n", + "Efficiency=(1-S/100);#unitless\n", + "OutputPower=InputPower*Efficiency;#in KW\n", + "\n", + "#Result\n", + "print \"Mechanical power output at a slip of 5%(in KW) : \",round(OutputPower,1);\n", + "#answer differs due to rounding" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Mechanical power output at a slip of 5%(in KW) : 24.7\n" + ] + } + ], + "prompt_number": 158 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.30, Page 162" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "#At 50% speed of synchronous speed\n", + "S=50.;#in %\n", + "UsefulOutput=55;#in H.P.\n", + "MechLoss=2;#in H.P.\n", + "StatorLoss=13.5;#in KW\n", + "\n", + "#Calculations&Results\n", + "Pout_rotor=UsefulOutput+MechLoss;#in H.P.\n", + "Pout_rotor=Pout_rotor*0.7355;#in KW\n", + "Efficiency=(1-S/100);#unitless\n", + "RotorInputPower=Pout_rotor/Efficiency;#in KW\n", + "RotorCuLoss=RotorInputPower-Pout_rotor;#in KW\n", + "print \"Cu Loss in the rotor circuit(in KW) :\",round(RotorCuLoss,2);\n", + "TotalLosses=StatorLoss+RotorCuLoss+MechLoss*0.7355;#in KW\n", + "MotorEfficiency=UsefulOutput*0.7355/(UsefulOutput*0.7355+TotalLosses);\n", + "print \"Motor Efficiency(in %) : \",round((MotorEfficiency*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Cu Loss in the rotor circuit(in KW) : 41.92\n", + "Motor Efficiency(in %) : 41.6\n" + ] + } + ], + "prompt_number": 159 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.31, Page 163" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "N=1440.;#in rpm(Actual speed of motor)\n", + "Power_dev=10;#Power developed in H.P.\n", + "VL=400;#in volt\n", + "cosfi=0.8;#powerfactor(unitless)\n", + "f=50;#in Hz\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#fractional slip(unitless)\n", + "print \"Slip :\",S;\n", + "Omega=N*2*pi/60;#angular speed in radians/sec\n", + "Torque=Power_dev*735.5/Omega;#in N-meter\n", + "Torque=Torque/9.81;#in Kg-meter\n", + "print \"Torque(in Kg-meter) :\",round(Torque,2);\n", + "RotorCuLoss=(S/(1-S))*Power_dev*735.5;#in watts\n", + "print \"Rotor Cu Loss per phase(in watts) : \",round((RotorCuLoss/3),1);\n", + "Pin_rotor=RotorCuLoss/S;#in watts\n", + "print \"Power input to rotor(in KW) :\",round((Pin_rotor/10**3),2);\n", + "LineCurrent=Pin_rotor/(sqrt(3)*VL*cosfi);#in Ampere\n", + "print \"Line Current(in A) :\",round(LineCurrent,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.04\n", + "Torque(in Kg-meter) : 4.97\n", + "Rotor Cu Loss per phase(in watts) : 102.2\n", + "Power input to rotor(in KW) : 7.66\n", + "Line Current(in A) : 13.82\n" + ] + } + ], + "prompt_number": 160 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.32, Page 164" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=440;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin_rotor=80;#in KW\n", + "fr=100;#revolutions/min\n", + "\n", + "#Calculations&Results\n", + "fr=100./60;#in Hz\n", + "S=fr/f;#slip(unitless)\n", + "print \"Slip : \",round(S,3);\n", + "Ns=120*f/P;#in rpm\n", + "N=Ns*(1-S);#in rpm(Speed of Rotor)\n", + "print \"Rotor speed(in rpm) : \",round(N);\n", + "RotorCuLoss=S*Pin_rotor*10**3;#in watts\n", + "print \"Rotor Cu Loss per phase(in watts) : \",round(RotorCuLoss/phase);\n", + "I2=65;#in Ampere\n", + "R2=(RotorCuLoss/phase)/I2**2;#in ohm\n", + "print \"Rotor resistance per phase(in ohm) : \",round(R2,2);\n", + "Pout_rotor=Pin_rotor-RotorCuLoss/10**3;#in KW\n", + "print \"Output power of rotor(in H.P.) : \",round((Pout_rotor/0.735),1);\n", + "#Note : answers of few part are not accurate in the book." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.033\n", + "Rotor speed(in rpm) : 967.0\n", + "Rotor Cu Loss per phase(in watts) : 889.0\n", + "Rotor resistance per phase(in ohm) : 0.21\n", + "Output power of rotor(in H.P.) : 105.2\n" + ] + } + ], + "prompt_number": 161 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.33, Page 165" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=440;#in volt\n", + "f=50.;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin_rotor=20;#in KW\n", + "N=960.;#in rpm(Speed of Motor)\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#slip(unitless)\n", + "RotorCuLoss=S*Pin_rotor*10**3;#in watts\n", + "RotorCuLoss=RotorCuLoss/phase;#in watts/phase\n", + "R2=1./3;#Rotor Resistance in ohm per phase\n", + "I2=sqrt(RotorCuLoss/R2)\n", + "\n", + "#Result\n", + "print \"Motor current per phase(in Ampere) :\",round(I2,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Motor current per phase(in Ampere) : 28.28\n" + ] + } + ], + "prompt_number": 162 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.34, Page 166" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "VL=500;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pout_rotor=20;#in H.P.\n", + "cosfi=0.87;#power factor\n", + "N=995.;#in rpm(Speed of Motor)\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#slip(unitless)\n", + "print \"Slip : \",S;\n", + "RotorCuLoss=(S/(1-S))*Pout_rotor*735.5;#in watts\n", + "print \"Rotor Cu Loss(in watts) :\",round(RotorCuLoss,2);\n", + "Pin_rotor=RotorCuLoss/S;#in watts\n", + "print \"Power input to rotor(in KW) :\",round((Pin_rotor/10**3),3);\n", + "IL=Pin_rotor/(sqrt(3)*VL*cosfi);#in Ampere\n", + "print \"Line current (in Ampere) :\",round(IL,2);\n", + "Rotorfreq=S*f;#in Hz\n", + "print \"Rotor frequency(in Hz) :\",round(Rotorfreq,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.005\n", + "Rotor Cu Loss(in watts) : 73.92\n", + "Power input to rotor(in KW) : 14.784\n", + "Line current (in Ampere) : 19.62\n", + "Rotor frequency(in Hz) : 0.25\n" + ] + } + ], + "prompt_number": 163 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.35, Page 166" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "Efficiency=85.;#in %\n", + "OutputPower=17.;#in KW\n", + "OutputPower=17*10**3;#in watts\n", + "StatorLosses=900.;#in watt\n", + "MechLosses=1100;#in watt\n", + "\n", + "#Calculations&Results\n", + "Pin=OutputPower/(Efficiency/100);#in watts\n", + "#Stator input or rotor input :\n", + "Pin_rotor=Pin-StatorLosses;#in watts\n", + "P_MechDev=OutputPower+MechLosses;#in watts\n", + "RotorCuLosses=Pin_rotor-P_MechDev;#in watts\n", + "S=RotorCuLosses/Pin_rotor;#slip(unitless)\n", + "print \"Slip :\",round(S,3);\n", + "Ns=120*f/P\n", + "N=Ns*(1-S);#in rpm(Speed of Motor)\n", + "Omega=2*pi*N/60;#angular speed in rad/sec\n", + "Torque=P_MechDev/Omega;#in Nm\n", + "print \"Torque developed(in Nm) :\",round(Torque,1);\n", + "Tm=OutputPower/Omega;#in Nm\n", + "Tm=Tm/9.81;#in Kg-meter\n", + "print \"Torque at the shaft(in Kg-meter) :\",round(Tm,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.052\n", + "Torque developed(in Nm) : 121.6\n", + "Torque at the shaft(in Kg-meter) : 11.64\n" + ] + } + ], + "prompt_number": 164 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.36, Page 168" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=500;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin_stator=50;#in KW\n", + "Statorlosses=2;#in KW\n", + "MechLosses=1;#in KW\n", + "\n", + "#Calculations&Results\n", + "Pin_rotor=Pin_stator-Statorlosses;#in KW\n", + "S=3.;#in %\n", + "RotorI2RLoss=(S/100)*Pin_rotor;#in KW\n", + "print \"Rotor I2R Loss(in KW) :\",RotorI2RLoss;\n", + "print \"Rotor core loss at 3% slip is very less and can be neglected.\";\n", + "Power_dev=Pin_rotor-RotorI2RLoss;#in KW\n", + "print \"Power developed by the rotor(in H.P.) : \",round((Power_dev/0.7355),1);\n", + "Outputpower=Power_dev-MechLosses;#in KW\n", + "Efficiency=Outputpower/Pin_stator;#unitless\n", + "print \"Efficiency of the motor(in %) :\",(Efficiency*100);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor I2R Loss(in KW) : 1.44\n", + "Rotor core loss at 3% slip is very less and can be neglected.\n", + "Power developed by the rotor(in H.P.) : 63.3\n", + "Efficiency of the motor(in %) : 91.12\n" + ] + } + ], + "prompt_number": 165 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.37, Page 168" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "MotorOutput=20;#in H.P.\n", + "MotorOutput=20*735.5;#in watts\n", + "S=4.;#full load slip in %\n", + "MechLosses=500;#in watts\n", + "\n", + "#Calculations&Results\n", + "Pdev_rotor=MotorOutput+MechLosses;#in watts\n", + "S=S/100;#fractional slip\n", + "RotorI2RLoss=(S/(1-S))*Pdev_rotor;#in watts\n", + "print \"Rotor I2R Loss(in watts) :\",RotorI2RLoss;\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "print \"Rotor speed(in rpm) :\",Nr;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor I2R Loss(in watts) : 633.75\n", + "Rotor speed(in rpm) : 1440.0\n" + ] + } + ], + "prompt_number": 80 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.38, Page 169" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "R2=0.02;#in ohm\n", + "X2=0.1;#in ohm\n", + "\n", + "#Calculations\n", + "S=R2/X2;#slip(unitless)\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "#At starting S=1\n", + "S=1;#slip\n", + "#Formula : T=K*S*R2/(R2^2+X2^2)\n", + "#Starting torque, Tst=K*R2/(R2^2+S^2*X2^2)\n", + "#Maximum torque, Tm=K/(2*X2)\n", + "#Tst=(2/3)*Tm : gives a equation\n", + "#100*R2^2-30*R2+1=0\n", + "a=100;b=-30;c=1;\n", + "R21=(-b+sqrt(b**2-4*a*c))/(2*a);#in ohm\n", + "R22=(-b-sqrt(b**2-4*a*c))/(2*a);#in ohm\n", + "#This R2 is the value of rotor circuit resistance.\n", + "RotWinResistance=0.02;#in ohm per phase\n", + "Extra_R1=R21-RotWinResistance;#in ohm\n", + "Extra_R2=R22-RotWinResistance;#in ohm\n", + "\n", + "#Result\n", + "print \"Extra Resistance(in ohm) : \",round(Extra_R1,3),\" ohm or \",round((Extra_R2),3),\" ohm.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Extra Resistance(in ohm) : 0.242 ohm or 0.018 ohm.\n" + ] + } + ], + "prompt_number": 166 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.39, Page 170" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "MotorShaftOutput=20;#in H.P.\n", + "\n", + "#Calculations&Results\n", + "MotorShaftOutput=20*735.5;#in watts\n", + "MechLosses=MotorShaftOutput*2/100;#in watts(2% of the output)\n", + "Pdev_rotor=MotorShaftOutput+MechLosses;#in watts\n", + "S=3.;#slip in %\n", + "S=S/100;#fractional slip\n", + "RotorI2RLoss=(S/(1-S))*Pdev_rotor;#in watts\n", + "print \"Rotor I2R Loss(in watts) :\",round(RotorI2RLoss);\n", + "print \"Rotor iron loss at 3% slip is very small and can be neglected.\";\n", + "RotorInput=Pdev_rotor+RotorI2RLoss;#in watts\n", + "print \"Rotor Input(in Watts) :\",round(RotorInput,1);\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "Nr=Nr/60;#in rps\n", + "T=MotorShaftOutput/(2*pi*Nr);#in Nm\n", + "print \"Output torque(in Nm) : \",round(T,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor I2R Loss(in watts) : 464.0\n", + "Rotor iron loss at 3% slip is very small and can be neglected.\n", + "Rotor Input(in Watts) : 15468.2\n", + "Output torque(in Nm) : 96.54\n" + ] + } + ], + "prompt_number": 167 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.40, Page 171" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=500;#in volt\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "Nr=975.;#in rpm\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-Nr)/Ns;#slip\n", + "print \"Slip :\",S;\n", + "Pin_stator=40;#in KW\n", + "StatorLosses=1;#in KW\n", + "Pin_rotor=Pin_stator-StatorLosses;#in KW\n", + "RotorCuLosses=S*Pin_rotor;#in KW\n", + "print \"Rotor Cu Losses(in KW) :\",RotorCuLosses;\n", + "RotorOutput=Pin_rotor-RotorCuLosses;#in KW\n", + "OutputHP=RotorOutput/0.735;#in H.P.\n", + "print \"Output Horse Power : \",round(OutputHP,2);\n", + "Efficiency=RotorOutput/Pin_stator;#unitless\n", + "print \"Efficiency(in %) :\",round(Efficiency*100);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.025\n", + "Rotor Cu Losses(in KW) : 0.975\n", + "Output Horse Power : 51.73\n", + "Efficiency(in %) : 95.0\n" + ] + } + ], + "prompt_number": 168 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.41, Page 172" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "VL=440;#in volt\n", + "f=50.;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3.;#no. of phase\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "fr=120;#alternations per minute\n", + "fr=fr/60;#in Hz\n", + "S=fr/f;#slip\n", + "print \"Slip : \",S;\n", + "Nr=Ns-S*Ns;#in rpm\n", + "print \"Rotor speed(in rpm) :\",Nr;\n", + "Rotor_input=80;#in KW\n", + "RotorCuLoss=S*Rotor_input;#in KW\n", + "print \"Rotor Cu Loss per phase(in watts) :\",round((RotorCuLoss*10**3)/phase);\n", + "P_Mechdev=Rotor_input*10**3-RotorCuLoss*10**3;#in watts\n", + "P_Mechdev=P_Mechdev/735.5;#in H.P.\n", + "print \"Mechanical power devloped(in H.P.) :\",round(P_Mechdev,1);\n", + "Ir=60;#in Ampere\n", + "R2=(RotorCuLoss*10**3/phase)/Ir**2;#in ohm\n", + "print \"Rotor resistance per phase(in ohm) :\",round(R2,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Slip : 0.04\n", + "Rotor speed(in rpm) : 960.0\n", + "Rotor Cu Loss per phase(in watts) : 1067.0\n", + "Mechanical power devloped(in H.P.) : 104.4\n", + "Rotor resistance per phase(in ohm) : 0.296\n" + ] + } + ], + "prompt_number": 169 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.43, Page 175 " + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50.;#in Hz\n", + "P=6;#no. of poles\n", + "R2=0.2;#rotor resistance per phase in ohm\n", + "Nr=960;#in rpm\n", + "Nr_dash=600;#in rpm\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-Nr)/Ns;#slip(unitless)\n", + "Sdash=(Ns-Nr_dash)/Ns;#slip(unitless)\n", + "#Let the new value of resistance is R2dash=(R+0.2)\n", + "R=R2*(Sdash/S)-R2;#Resistance to be added in ohm\n", + "\n", + "#Result\n", + "print \"Resistance to be added(in ohm) : \",R;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Resistance to be added(in ohm) : 1.8\n" + ] + } + ], + "prompt_number": 91 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.44, Page 176" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=6;#no. of poles\n", + "phase=3;#no. of phase\n", + "RotorIntResistance=0.1;#in ohm per phase\n", + "X2=1;#in ohm/phase\n", + "\n", + "#Calculation\n", + "Rext=X2-RotorIntResistance;#in ohm\n", + "\n", + "#Result\n", + "print \"External resistance to be included(in ohm/phase) : \",Rext;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "External resistance to be included(in ohm/phase) : 0.9\n" + ] + } + ], + "prompt_number": 92 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.45, Page 176" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50.;#in Hz\n", + "P=6.;#no. of poles\n", + "phase=3.;#no. of phase\n", + "R2=0.2;#rotor resistance per phase in ohm\n", + "N1=960.;#Full load speed in rpm\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "S1=(Ns-N1)/Ns;#Full load slip(unitless)\n", + "N2=N1*(1-10./100);#New speed in rpm(reduced 10%)\n", + "S2=(Ns-N2)/Ns;#New slip(unitless)\n", + "#Formula : S=RotorCuLoss/Pin_rotor=3*I2^2*R2/Pin_rotor\n", + "#Let the additional resistance is R\n", + "R=R2*S2/S1-R2;#Resistance to be added in ohm\n", + "\n", + "#Result\n", + "print \"Additional Rotor Resistance(in ohm) : \",R;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Additional Rotor Resistance(in ohm) : 0.48\n" + ] + } + ], + "prompt_number": 97 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.46, Page 180" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "R2inner=0.4;#in ohm\n", + "X2inner=2;#in ohm \n", + "R2outer=2;#in ohm\n", + "X2outer=0.4;#in ohm \n", + "#At standstill :\n", + "S=1;#unitless\n", + "\n", + "#Calculations&Results\n", + "Ratio=(R2outer*S/(R2outer**2+S**2*X2outer**2))/(R2inner*S/(R2inner**2+S**2*X2inner**2));#unitless\n", + "print \"Ratio of torque produced by two cages at standstill : \",Ratio;\n", + "#At Full load :\n", + "S=5;#in %\n", + "S=5./100;#fractional\n", + "Ratio=(R2outer*S/(R2outer**2+S**2*X2outer**2))/(R2inner*S/(R2inner**2+S**2*X2inner**2));#unitless\n", + "print \"Ratio of torque produced by two cages at full load(slip=5%) : \",round(Ratio,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ratio of torque produced by two cages at standstill : 5.0\n", + "Ratio of torque produced by two cages at full load(slip=5%) : 0.21\n" + ] + } + ], + "prompt_number": 171 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.47, Page 187" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "Output=10;#in H.P.\n", + "Output=Output*735.5;#in watts\n", + "cosfi=0.8;#unitless\n", + "ETA=0.83;#unitless\n", + "ISCbyIFL=3.5;#ratio of SC current to full load current\n", + "VL=500;#in volt\n", + "\n", + "#Calculations\n", + "Input=Output/ETA;#in watts\n", + "IFL=Input/(sqrt(3)*VL*cosfi);#in Ampere\n", + "ISC=IFL*ISCbyIFL;#in Ampere\n", + "Is=ISC/3;#in Ampere\n", + "\n", + "#Result\n", + "print \"Strting current(in Ampere) :\",round(Is,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Strting current(in Ampere) : 14.92\n" + ] + } + ], + "prompt_number": 172 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.48, Page 187" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=5;#ratio of SC current to full load current\n", + "S=0.04;#Full load slip\n", + "\n", + "#Calculations&Results\n", + "#(i) for star delta starter : \n", + "print \"for star delta starter : \";\n", + "Ratio=(1./3)*(ISCbyIFL)**2*S;#Ratio of starting torque to full load torque\n", + "print \"Starting torque is \",round((Ratio*100),2),\"% of full load torque.\";\n", + "#(i) for auto transformer starter :\n", + "print \"for auto transformer starter : \";\n", + "K=50./100;#tappings\n", + "Ratio=K**2*(ISCbyIFL)**2*S;#Ratio of starting torque to full load torque\n", + "print \"Starting torque is \",(Ratio*100),\"% of full load torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "for star delta starter : \n", + "Starting torque is 33.33 % of full load torque.\n", + "for auto transformer starter : \n", + "Starting torque is 25.0 % of full load torque.\n" + ] + } + ], + "prompt_number": 173 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.49, Page 188" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=4;#ratio of SC current to full load current\n", + "TsBYTf=1./4;#ratio of atarting torque to full load torque\n", + "Slip=3;#in %\n", + "Slip=3./100;#in fraction\n", + "\n", + "#Calculations\n", + "#Formula : TsBYTf=Percent_Tapping^2*ISCbyIFL^2*Slip\n", + "tapping=sqrt(TsBYTf/(Slip*ISCbyIFL**2));#in %\n", + "\n", + "#Result\n", + "print \"Tapping(in %) :\",round((tapping*100),1);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Tapping(in %) : 72.2\n" + ] + } + ], + "prompt_number": 174 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.50, Page 189" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=5.;#ratio of SC current to full load current\n", + "Slip=5;#in %\n", + "Slip=5./100;#in fraction\n", + "\n", + "#Calculations\n", + "#Formula : 3*IFL=K^2*ISC\n", + "K=sqrt(1./(ISCbyIFL/3));#unitless\n", + "TsBYTfl=K**2*(ISCbyIFL)**2*Slip;#ratio of starting torque to full load torque\n", + "\n", + "#Result\n", + "print \"Starting torque is \",(TsBYTfl*100),\"% of full load torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Starting torque is 75.0 % of full load torque.\n" + ] + } + ], + "prompt_number": 175 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.51, Page 189" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "print \"Star delta starter :\";\n", + "ISCbyIFL=6;#ratio of SC current to full load current\n", + "Slip=4;#in %\n", + "Slip=4./100;#in fraction\n", + "\n", + "#Calculations&Results\n", + "TsBYTfl=(1./3)*(ISCbyIFL)**2*Slip;#ratio of starting torque to full load torque\n", + "print \"Starting torque is \",(TsBYTfl*100),\"% of full load value.\";\n", + "print \"For an auto transformer :\";\n", + "K=70.70;#in %\n", + "K=70.70/100;#in fraction\n", + "TsBYTfl=K**2*(ISCbyIFL)**2*Slip;#ratio of starting torque to full load torque\n", + "print \"Starting torque is \",round(TsBYTfl*100),\"% of full load torque.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Star delta starter :\n", + "Starting torque is 48.0 % of full load value.\n", + "For an auto transformer :\n", + "Starting torque is 72.0 % of full load torque.\n" + ] + } + ], + "prompt_number": 176 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.52, Page 190" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "V=400.;#in volt\n", + "S=3.5;#in %\n", + "S=3.5/100;#in fraction\n", + "VL=92;#in volts\n", + "\n", + "#Calculations&Results\n", + "#ISC=(V/VL)*IFL;#in Ampere\n", + "ISCbyIFL=V/VL;#ratio of SC current to full load current\n", + "#2*IFL=K^2*ISC imples that 2*IFL=K^2*(V/VL)*IFL\n", + "K=sqrt(2./(V/VL));#in fraction\n", + "print \"Necessary tapping(in %) :\",round((K*100),1);\n", + "TsBYTfl=K**2*(ISCbyIFL)**2*S;#ratio of starting torque to full load torque\n", + "print \"Starting torque is \",round((TsBYTfl*100),1),\"% of full load value.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Necessary tapping(in %) : 67.8\n", + "Starting torque is 30.4 % of full load value.\n" + ] + } + ], + "prompt_number": 179 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.53, Page 190" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=4;#ratio of SC current to full load current\n", + "ISbyIFL=2;#ratio of Supply current to full load current\n", + "S=2.5;#in %\n", + "S=2.5/100;#in fraction\n", + "\n", + "#Calculations\n", + "#Formula : ISbyIFL=(1/K^2)*ISCbyIFL\n", + "K=sqrt(ISCbyIFL/ISbyIFL);#in fraction\n", + "print \"Auto Transormation ratio :\",round(K,2);\n", + "TsBYTfl=(1./K**2)*(ISCbyIFL)**2*S;#ratio of starting torque to full load torque\n", + "\n", + "#Results\n", + "print \"Ratio of starting torque to full load torque :\",TsBYTfl;\n", + "print \"or Starting torque is \",(TsBYTfl*100),\"% of full load value.\";" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Auto Transormation ratio : 1.41\n", + "Ratio of starting torque to full load torque : 0.2\n", + "or Starting torque is 20.0 % of full load value.\n" + ] + } + ], + "prompt_number": 180 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.54, Page 191" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "ISCbyIFL=4;#ratio of SC current to full load current\n", + "S=3;#in %\n", + "S=3./100;#in fraction\n", + "TsBYTfl=1./4;#ratio of starting torque to full load torque\n", + "\n", + "#Calculations\n", + "#Formula : TsBYTfl=ISCbyIFL^2*K^2*S\n", + "K=sqrt(TsBYTfl/(ISCbyIFL**2*S));#in fraction\n", + "\n", + "#Result\n", + "print \"Auto Transormation ratio(in %) :\",round((K*100),2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Auto Transormation ratio(in %) : 72.17\n" + ] + } + ], + "prompt_number": 181 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.55, Page 191" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "Output=3;#in H.P.\n", + "Output=3*735.5;#in watts\n", + "Efficiency=0.83;#unitless\n", + "cosfi=0.8;#power factor\n", + "\n", + "#Calculations\n", + "Vl=500;#in volt\n", + "Input=Output/Efficiency;#in watts\n", + "#Formula : Input=sqrt(3)*Vl*Il*cosfi\n", + "Il=Input/(sqrt(3)*Vl*cosfi);#in Ampere\n", + "ISCbyIFL=3.5;#ratio of SC current to full load current\n", + "ISC=ISCbyIFL*Il;#in Ampere\n", + "LineCurrent=ISC/3;#in Ampere(for star delta starter)\n", + "\n", + "#Result\n", + "print \"Line Current(in Ampere) :\",round(LineCurrent,2);\n", + "#Note : Ans in the book is not accurate." + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Line Current(in Ampere) : 4.48\n" + ] + } + ], + "prompt_number": 182 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.56, Page 192" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "Output=15;#in H.P.\n", + "Output=15*735.5;#in watts\n", + "Vl=400;#in volt\n", + "ISCat200=40;#in Ampere(at 200 volt)\n", + "Efficiency=0.88;#unitless\n", + "cosfi=0.85;#power factor\n", + "\n", + "#Calculations\n", + "ISCat400=ISCat200*(400./200);#in Ampere(at 400 volt)\n", + "Input=Output/Efficiency;#in watts\n", + "Ifl=Input/(sqrt(3)*Vl*cosfi);#in Ampere\n", + "#starting line current with star delta starter\n", + "Is=ISCat400/3;#in Ampere\n", + "Ratio=Is/Ifl;#ratio of starting current to full load current\n", + "\n", + "#Result\n", + "print \"Ratio of line current at starting to full load current :\",round(Ratio,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Ratio of line current at starting to full load current : 1.25\n" + ] + } + ], + "prompt_number": 183 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.57, Page 192" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import * \n", + "\n", + "#Given data :\n", + "#With star delta starter : \n", + "TstBYTfl=0.35;#ratio of starting torque to full load torque\n", + "IstBYIfl=1.75;#ratio of starting current to full load current\n", + "\n", + "#Calculations&Results\n", + "ISCBYIs=sqrt(3);#ratio of SC current to starting current\n", + "ISCBYIfl=sqrt(3)*IstBYIfl;#ratio of SC current to full load current\n", + "#Formula : TstBYTfl=(ISCBYIfl)^2*S\n", + "S=TstBYTfl/(ISCBYIfl)**2;#in fraction\n", + "print \"Full load Slip : \",round(S,3);\n", + "#With auto transformer with winding in delta : \n", + "Ip=sqrt(3)*1.750*0.8;#full voltage phase current in Ampere\n", + "IlBYIf=4.2;#ratio of Line current to full load current\n", + "Ratio=IlBYIf**2*S;#ratio of starting current to full load current\n", + "print \"Ratio of line current at starting to full load current :\",round(Ratio,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Full load Slip : 0.038\n", + "Ratio of line current at starting to full load current : 0.67\n" + ] + } + ], + "prompt_number": 188 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.58, Page 196" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "R2=0.25;#in ohm per phase\n", + "S=2;#in %\n", + "S=2./100;#in fraction\n", + "\n", + "#Calculations\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "#When speed reduced to 10%\n", + "NewSpeed=Nr*90/100;#in rpm\n", + "Sdash=(Ns-NewSpeed)/Ns;#in fraction\n", + "R2dash=(Sdash/S)*R2;#in ohm per phase\n", + "R=R2dash-R2;#in ohm\n", + "\n", + "#Result\n", + "print \"Value of resistance to be added(in ohm) :\",R;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Value of resistance to be added(in ohm) : 1.225\n" + ] + } + ], + "prompt_number": 114 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.59, Page 197" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "R2=0.25;#in ohm per phase\n", + "N=1440.;#in rpm at full load\n", + "NewSpeed=1200;#in rpm\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "S=(Ns-N)/Ns;#in %\n", + "NewS=(Ns-NewSpeed)/Ns;#in fraction\n", + "S1=S;S2=NewS;#slip new and old\n", + "#Torque remaining same : S1/R2=S2/(R2+R)\n", + "R=S2*R2/S1-R2;#in ohm\n", + "print \"External resistance per phase(in ohm) :\",R; #incorrect answer in textbook\n", + "Nr=Ns-Ns*S;#in rpm\n", + "#If S1 is taken as 0.03 and S2 as 0.127\n", + "S1=0.03;#slip in fraction\n", + "S2=0.127;#slip in fraction\n", + "R=S2*R2/S1-R2;#in ohm\n", + "print \"External resistance per phase(in ohm) :\",round(R,3);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "External resistance per phase(in ohm) : -0.25\n", + "External resistance per phase(in ohm) : 0.808\n" + ] + } + ], + "prompt_number": 194 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.60, Page 198" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "S=4;#in %\n", + "S=4./100;#in fraction\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "N1=Ns-Ns*S;#in rpm\n", + "#When speed reduced to 10%\n", + "N2=N1*85/100;#in rpm(NewSpeed)\n", + "print \"New speed(in rpm) :\",N2;\n", + "#New speed is reduced by 15 %\n", + "Sdash=(Ns-N2)/Ns;#in fraction\n", + "print \"New Slip : \",Sdash;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "New speed(in rpm) : 1224.0\n", + "New Slip : 0.184\n" + ] + } + ], + "prompt_number": 118 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 2.61, Page 202" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "f=50;#in Hz\n", + "P=4;#no. of poles\n", + "phase=3;#no. of phase\n", + "Pin=60;#in KW\n", + "StatorLoss=1.2;#in KW\n", + "WindingLoss=1.8;#in KW\n", + "S=4;#in %\n", + "S=4./100;#in fraction\n", + "\n", + "#Calculations&Results\n", + "Ns=120*f/P;#in rpm\n", + "Nr=Ns-Ns*S;#in rpm\n", + "print \"Rotor speed(in rpm) :\",Nr;\n", + "Pin_rotor=Pin-StatorLoss;#in KW\n", + "RotorCuLoss=S*Pin_rotor;#in KW\n", + "print \"Rotor Cu Loss(in KW) : \",RotorCuLoss;\n", + "Pout_rotor=Pin_rotor-WindingLoss-RotorCuLoss;#in KW\n", + "Efficiency=(Pout_rotor/Pin)*100;#in %\n", + "print \"Efficiency(in %) :\",round(Efficiency);\n", + "HP=Pout_rotor*10**3/735.5;#in H.P.\n", + "print \"Horse Power(inH.P.) :\",round(HP,2);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Rotor speed(in rpm) : 1440.0\n", + "Rotor Cu Loss(in KW) : 2.352\n", + "Efficiency(in %) : 91.0\n", + "Horse Power(inH.P.) : 74.3\n" + ] + } + ], + "prompt_number": 195 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/chapter_3.ipynb b/Electrical_Machines_II/chapter_3.ipynb new file mode 100755 index 00000000..2c125106 --- /dev/null +++ b/Electrical_Machines_II/chapter_3.ipynb @@ -0,0 +1,68 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:d44762913f69b2b88f8f79d544b4c89ac298c3f68155a24b93b531ef54634cd3" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 3: Fractional K.W. \n", + "(FKW) Motors" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.1, Page 242" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import *\n", + "from cmath import *\n", + "\n", + "#Given data :\n", + "Za=complex(3,3);#in ohm\n", + "Zm=complex(6,3);#in ohm\n", + "\n", + "#Calculations\n", + "#Phase difference 90 degree : so angle between Im and V is 90-45 degree\n", + "theta=90-45;#in degree\n", + "f=50;#in Hz\n", + "#Formula : tand(theta)=(Xc-Xl)/R\n", + "Xl=3;#in ohm and Xc=1/(2*%pi*C)\n", + "R=6;#in ohm\n", + "C=1./(2*pi*f*(tan(theta*pi/180)*R+Xl));\n", + "\n", + "#Result\n", + "print \"Capacitance(in micro Farad) :\",(C*10**6);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance(in micro Farad) : (353.677651315+0j)\n" + ] + } + ], + "prompt_number": 7 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/chapter_3_1.ipynb b/Electrical_Machines_II/chapter_3_1.ipynb new file mode 100755 index 00000000..2c125106 --- /dev/null +++ b/Electrical_Machines_II/chapter_3_1.ipynb @@ -0,0 +1,68 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:d44762913f69b2b88f8f79d544b4c89ac298c3f68155a24b93b531ef54634cd3" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 3: Fractional K.W. \n", + "(FKW) Motors" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 3.1, Page 242" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from math import *\n", + "from cmath import *\n", + "\n", + "#Given data :\n", + "Za=complex(3,3);#in ohm\n", + "Zm=complex(6,3);#in ohm\n", + "\n", + "#Calculations\n", + "#Phase difference 90 degree : so angle between Im and V is 90-45 degree\n", + "theta=90-45;#in degree\n", + "f=50;#in Hz\n", + "#Formula : tand(theta)=(Xc-Xl)/R\n", + "Xl=3;#in ohm and Xc=1/(2*%pi*C)\n", + "R=6;#in ohm\n", + "C=1./(2*pi*f*(tan(theta*pi/180)*R+Xl));\n", + "\n", + "#Result\n", + "print \"Capacitance(in micro Farad) :\",(C*10**6);" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Capacitance(in micro Farad) : (353.677651315+0j)\n" + ] + } + ], + "prompt_number": 7 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/chapter_4.ipynb b/Electrical_Machines_II/chapter_4.ipynb new file mode 100755 index 00000000..632c52df --- /dev/null +++ b/Electrical_Machines_II/chapter_4.ipynb @@ -0,0 +1,183 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:9e27d5dc9ee72aa8af2a87478e959a0f346f3cc04c2f1acb86835886294034d4" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 4: Special Purpose Machines" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.1, Page 281" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=8;#no. of poles\n", + "T=5;#no. of teeteh in each pole\n", + "Nr=50.;#no. of teeth in rotor\n", + "\n", + "#Calculations\n", + "Ns=P*T;#no. of teeth\n", + "Beta=(Nr-Ns)*360/(Nr*Ns);#in degree\n", + "\n", + "#Result\n", + "print \"Stepping Angle(in degree) : = \",Beta;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stepping Angle(in degree) : = 1.8\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.2, Page 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Beta=2.5;#in degree\n", + "f=3600;#in PPs\n", + "Resolution=360;#in degree\n", + "\n", + "#Calculations&Results\n", + "print \"Resolution : \",Resolution;\n", + "Beta1=Resolution/Beta;#steps/revolution\n", + "print \"Steps per revolution : \",Beta1;\n", + "n=Beta*f/360;#in rps\n", + "print \"Steps required for making 25 revolution(in rps) : \",n;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Resolution : 360\n", + "Steps per revolution : 144.0\n", + "Steps required for making 25 revolution(in rps) : 25.0\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.3, Page 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Beta=15;#in degree\n", + "Phase=3;#no. of phase\n", + "\n", + "#Calculations&Results\n", + "#Formula : Beta*Nr=(360/phase)\n", + "Nr=(360./Phase)/Beta;#unitless\n", + "print \"No. of rotor :\",Nr;\n", + "#Case (i) : when Ns>Nr : Formula : Beta=(Ns-Nr)*360/(Ns*Nr)\n", + "Ns=Nr/(1-Beta*Nr/360);\n", + "print \"No. of stator poles if Ns>Nr :\",Ns;\n", + "#Case (ii) : when Nr>Ns : Formula : Beta=(Nr-Ns)*360/(Ns*Nr)\n", + "Ns=Nr/(1+Beta*Nr/360);\n", + "print \"No. of stator poles if Nr>Ns :\",Ns;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "No. of rotor : 8.0\n", + "No. of stator poles if Ns>Nr : 12.0\n", + "No. of stator poles if Nr>Ns : 6.0\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.4, Page 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Beta=1.8;#in degree\n", + "m=4;#no. of phases\n", + "\n", + "#Calculations\n", + "#Formula : Beta*Nr=(360/m)\n", + "Nr=(360/m)/Beta;#unitless\n", + "\n", + "#Results\n", + "print \"No. of rotor teeth : = \",Nr;\n", + "print \"In multi stack motor rotor teeth equals to stator teeth equals to : \",Nr;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "No. of rotor teeth : = 50.0\n", + "In multi stack motor rotor teeth equals to stator teeth equals to : 50.0\n" + ] + } + ], + "prompt_number": 5 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/chapter_4_1.ipynb b/Electrical_Machines_II/chapter_4_1.ipynb new file mode 100755 index 00000000..632c52df --- /dev/null +++ b/Electrical_Machines_II/chapter_4_1.ipynb @@ -0,0 +1,183 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:9e27d5dc9ee72aa8af2a87478e959a0f346f3cc04c2f1acb86835886294034d4" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 4: Special Purpose Machines" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.1, Page 281" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "P=8;#no. of poles\n", + "T=5;#no. of teeteh in each pole\n", + "Nr=50.;#no. of teeth in rotor\n", + "\n", + "#Calculations\n", + "Ns=P*T;#no. of teeth\n", + "Beta=(Nr-Ns)*360/(Nr*Ns);#in degree\n", + "\n", + "#Result\n", + "print \"Stepping Angle(in degree) : = \",Beta;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Stepping Angle(in degree) : = 1.8\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.2, Page 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Beta=2.5;#in degree\n", + "f=3600;#in PPs\n", + "Resolution=360;#in degree\n", + "\n", + "#Calculations&Results\n", + "print \"Resolution : \",Resolution;\n", + "Beta1=Resolution/Beta;#steps/revolution\n", + "print \"Steps per revolution : \",Beta1;\n", + "n=Beta*f/360;#in rps\n", + "print \"Steps required for making 25 revolution(in rps) : \",n;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Resolution : 360\n", + "Steps per revolution : 144.0\n", + "Steps required for making 25 revolution(in rps) : 25.0\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.3, Page 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Beta=15;#in degree\n", + "Phase=3;#no. of phase\n", + "\n", + "#Calculations&Results\n", + "#Formula : Beta*Nr=(360/phase)\n", + "Nr=(360./Phase)/Beta;#unitless\n", + "print \"No. of rotor :\",Nr;\n", + "#Case (i) : when Ns>Nr : Formula : Beta=(Ns-Nr)*360/(Ns*Nr)\n", + "Ns=Nr/(1-Beta*Nr/360);\n", + "print \"No. of stator poles if Ns>Nr :\",Ns;\n", + "#Case (ii) : when Nr>Ns : Formula : Beta=(Nr-Ns)*360/(Ns*Nr)\n", + "Ns=Nr/(1+Beta*Nr/360);\n", + "print \"No. of stator poles if Nr>Ns :\",Ns;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "No. of rotor : 8.0\n", + "No. of stator poles if Ns>Nr : 12.0\n", + "No. of stator poles if Nr>Ns : 6.0\n" + ] + } + ], + "prompt_number": 7 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 4.4, Page 282" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "#Given data :\n", + "Beta=1.8;#in degree\n", + "m=4;#no. of phases\n", + "\n", + "#Calculations\n", + "#Formula : Beta*Nr=(360/m)\n", + "Nr=(360/m)/Beta;#unitless\n", + "\n", + "#Results\n", + "print \"No. of rotor teeth : = \",Nr;\n", + "print \"In multi stack motor rotor teeth equals to stator teeth equals to : \",Nr;" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "No. of rotor teeth : = 50.0\n", + "In multi stack motor rotor teeth equals to stator teeth equals to : 50.0\n" + ] + } + ], + "prompt_number": 5 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file diff --git a/Electrical_Machines_II/screenshots/capacitance_3.png b/Electrical_Machines_II/screenshots/capacitance_3.png Binary files differnew file mode 100755 index 00000000..4bf742ff --- /dev/null +++ b/Electrical_Machines_II/screenshots/capacitance_3.png diff --git a/Electrical_Machines_II/screenshots/capacitance_3_1.png b/Electrical_Machines_II/screenshots/capacitance_3_1.png Binary files differnew file mode 100755 index 00000000..4bf742ff --- /dev/null +++ b/Electrical_Machines_II/screenshots/capacitance_3_1.png diff --git a/Electrical_Machines_II/screenshots/flux_per_pole.png b/Electrical_Machines_II/screenshots/flux_per_pole.png Binary files differnew file mode 100755 index 00000000..53a57349 --- /dev/null +++ b/Electrical_Machines_II/screenshots/flux_per_pole.png diff --git a/Electrical_Machines_II/screenshots/flux_per_pole_1.png b/Electrical_Machines_II/screenshots/flux_per_pole_1.png Binary files differnew file mode 100755 index 00000000..53a57349 --- /dev/null +++ b/Electrical_Machines_II/screenshots/flux_per_pole_1.png diff --git a/Electrical_Machines_II/screenshots/rotor_current_2.png b/Electrical_Machines_II/screenshots/rotor_current_2.png Binary files differnew file mode 100755 index 00000000..4b8c4b2b --- /dev/null +++ b/Electrical_Machines_II/screenshots/rotor_current_2.png diff --git a/Electrical_Machines_II/screenshots/rotor_current_2_1.png b/Electrical_Machines_II/screenshots/rotor_current_2_1.png Binary files differnew file mode 100755 index 00000000..4b8c4b2b --- /dev/null +++ b/Electrical_Machines_II/screenshots/rotor_current_2_1.png |
