diff options
author | Prashant S | 2020-04-14 10:25:32 +0530 |
---|---|---|
committer | GitHub | 2020-04-14 10:25:32 +0530 |
commit | 06b09e7d29d252fb2f5a056eeb8bd1264ff6a333 (patch) | |
tree | 2b1df110e24ff0174830d7f825f43ff1c134d1af /Electrical_Engineering_by_R_Anandanatarajan | |
parent | abb52650288b08a680335531742a7126ad0fb846 (diff) | |
parent | 476705d693c7122d34f9b049fa79b935405c9b49 (diff) | |
download | all-scilab-tbc-books-ipynb-master.tar.gz all-scilab-tbc-books-ipynb-master.tar.bz2 all-scilab-tbc-books-ipynb-master.zip |
Initial commit
Diffstat (limited to 'Electrical_Engineering_by_R_Anandanatarajan')
5 files changed, 3937 insertions, 0 deletions
diff --git a/Electrical_Engineering_by_R_Anandanatarajan/1-DC_Machines.ipynb b/Electrical_Engineering_by_R_Anandanatarajan/1-DC_Machines.ipynb new file mode 100644 index 0000000..b4eb286 --- /dev/null +++ b/Electrical_Engineering_by_R_Anandanatarajan/1-DC_Machines.ipynb @@ -0,0 +1,1615 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 1: DC Machines" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.10: To_calculate_the_properties_of_shunt_generator.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.10, Page 1.26\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"IL=180;//Load current in A\n", +"V=220;//Terminal voltage in V\n", +"Ra=0.01;//Armature resistance in ohm\n", +"Rsh=40;//Shield field resistance in ohm\n", +"Wc=1000;//Constant losses in W\n", +"x=185;//Load current in A\n", +"\n", +"//CALCULATIONS\n", +"Ia=(IL+(V/Rsh));//Armature current in A\n", +"Eg=(V+(Ia*Ra));//Generated emf in V\n", +"Pm=(V*x)+Wc+(Ia^2*Ra)+(V^2/Rsh);//Output of the prime mover in W\n", +"nm=((V*Ia)/Pm)*100;//Mechanical efficiency\n", +"ne=((V*IL)/(Eg*Ia))*100;//Electrical efficiency\n", +"no=((V*IL)/(Pm))*100;//Overall efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Generated emf is %3.3f V \n b)Output of the prime mover is %3.2f W \n c)Mechanical efficiency is %3.2f percent \n d)Electrical efficiency is %3.2f percent \n e)Overall efficiency is %3.2f percent',Eg,Pm,nm,ne,no)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.11: To_calculate_the_back_emf_and_total_mechanical_power.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.11, Page 1.33\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"\n", +"//INPUT DATA\n", +"IL=15;//Load current in A\n", +"V=220;//Terminal voltage in V\n", +"Rsh=180;//Field resistance in ohm\n", +"Ra=0.03;//Armature resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Ish=(V/Rsh);//Field current in A\n", +"Ia=(IL-Ish);//Armature current in A\n", +"Eb=(V-(Ia*Ra));//Back emf in V\n", +"Pm=(Eb*Ia)/1000;//Total mechanical power in kW\n", +"\n", +"//OUTPUT\n", +"mprintf('i)Back emf is %3.2f V \nii)Total mechanical power developed in the armature is %3.2f kW',Eb,Pm)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.12: To_find_the_change_in_back_emf.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.12, Page 1.34\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=220;//Terminal voltage in V\n", +"IaFL=25;//Full load armature current in A\n", +"IaNL=5;//No load armature current in A\n", +"Ra=0.5;//Armature resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"EbNL=(V-(IaNL*Ra));//Back emf at no load in V\n", +"Eb=(V-(IaFL*Ra));//Back emf at full load in V\n", +"E=(EbNL-Eb);//Change in back emf from no load to full load in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Change in back emf from no load to full load is %3.0f V',E)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.13: To_find_the_speed_of_the_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.13, Page 1.34\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=4;//Number of poles\n", +"V=500;//Terminal voltage in V\n", +"Ia=80;//Armature current in A\n", +"Ra=0.4;//Armature resistance in ohm\n", +"A=2;//Number of parallel paths\n", +"Z=522;//Number of conductors\n", +"q=0.025;//Useful flux per pole in Wb\n", +"\n", +"//CALCULATIONS\n", +"Eb=(V-(Ia*Ra));//Back emf in V\n", +"N=(Eb*60*A)/(P*q*Z);//Speed of the motor in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Speed of the motor is %3.1f rpm',N)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.14: To_find_the_armature_resistance_and_maximum_armature_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.14, Page 1.35\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Eb=225;//Back emf in V\n", +"IL=40;//Line current in A\n", +"Rsh=150;//Field resistance in ohm\n", +"Ish=1.67;//Field current in A\n", +"\n", +"//CALCULATIONS\n", +"V=(Ish*Rsh);//Terminal applied voltage in V\n", +"Ia=(IL-Ish);//Armature current in A\n", +"Ra=(V-Eb)/Ia;//Armature resistance in ohm\n", +"Ia=(V/Ra);//Maximum armature current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('i)Armature resistance is %3.2f ohm \nii)Armature current will be maximum at the moment of start up and it is %3.2f A',Ra,Ia)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.15: To_find_the_back_emf.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.15, Page 1.36\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=400;//Terminal voltage in V\n", +"P=8000;//Motor output power in W\n", +"n=0.9;//Motor efficiency\n", +"Rsh=180;//Field resistance in ohm\n", +"Ra=0.6;//Armature resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"If=(V/Rsh);//Field current in A\n", +"Pi=(P/n);//Input power in W\n", +"IL=(Pi/V);//Load current in A\n", +"Ia=(IL-If);//Armature current in A\n", +"Eb=(V-(Ia*Ra));//Back emf in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Back emf is %3.0f V',Eb)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.16: To_find_the_total_power_developed.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.16, Page 1.37\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=30000;//Power rating of the dc machine in W\n", +"V=300;//Terminal voltage in V\n", +"Ra=0.04;//Armature resistance in ohm\n", +"Rsh=120;//Shunt field resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"IL=(P/V);//Load current in A\n", +"Ia=(IL+(V/Rsh));//Armature current in A\n", +"Eg=(V+(Ia*Ra));//Generated emf in V\n", +"P=(Eg*Ia);//Power developed in the armature in W\n", +"Ish=(V/Rsh);//Field current in A\n", +"Ia2=(IL-Ish);//Armature current in motor in A\n", +"Eb=(V-(Ia2*Ra));//Back emf in V\n", +"P1=(Eb*Ia2);//Power developed in the armature in W\n", +"\n", +"//OUTPUT\n", +"mprintf('Total power developed in the armature when \ni)the dc machine is operated as a generator is %3.0f W \nii)when the dc machine is operated as a motor is %3.1f W',P,P1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.17: To_find_the_armature_torque_and_armature_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.17, Page 1.43\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=4;//Number of poles\n", +"Z=726;//Number of conductors \n", +"A=2;//Number of parallel paths\n", +"q=(30*10^-3);//Flux per pole in Wb\n", +"Ia=45;//Total armature current in A\n", +"\n", +"//CALCULATIONS\n", +"Ta=(0.159*Z*q*Ia*P)/A;//Armature toque in N.m\n", +"\n", +"//OUTPUT\n", +"mprintf('Armature torque is %3.2f N.m',Ta)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.18: To_calculate_the_current_taken_and_diameter_of_the_motor_pulley.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.18, Page 1.43\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N=(1800/60);//Speed of the motor in rps\n", +"V=200;//Terminal voltage in V\n", +"Nl=(900/60);//Lathe speed in rps\n", +"F=300;//Force exerted in N\n", +"r=0.2;//Radius of the shaft in m\n", +"n=0.9;//Efficiency of the motor\n", +"Dp=0.3;//Diamter of the Lathe pulley in m\n", +"\n", +"//CALCULATIONS\n", +"Tsh=(F*r);//Shaft torque in N.m\n", +"Psh=(Tsh*2*3.14*Nl);//Shaft power in W\n", +"Pi=(Psh/n);//Input power in W\n", +"I=(Pi/V);//Current taken by the motor in A\n", +"Dm=((Nl*Dp)/N)*100;//Diameter of the motor pulley in cm\n", +"\n", +"//OUTPUT\n", +"mprintf('Current taken by the motor is %3.1f A \nDiameter of the motor pulley is %3.0f cm',I,Dm)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.19: To_calculate_the_armature_torque_and_horse_power.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.19, Page 1.45\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N=(300/60);//Speed of the motor in rps\n", +"P=4;//Number of poles\n", +"Z=732;//Number of conductors \n", +"I=80;//Current through each conductor in A\n", +"l=0.35;//Length of the conductor in m\n", +"n=0.8;//Efficiency of flux distribution\n", +"B=0.8;//Flux densty in Wb/m^2\n", +"D=0.8;//Diameter of the armature in m\n", +"\n", +"//CALCULATIONS\n", +"Ze=(Z*n);//Number of effective conductors\n", +"q=(B*l*2*3.14*(D/2))/4;//Flux per pole in Wb\n", +"Ta=(0.159*q*Ze*I*P);//Armature torque in N.m\n", +"F=(B*I*l);//Force on each conductor in N\n", +"T=(F*(D/2));//Torque due to on econductor in N.m\n", +"T1=(T*Ze);//Torque due to all conductors in N.m\n", +"Br=(2*3.14*N*60*T1)/(60*746);//Brake Horse power in HP\n", +"\n", +"//OUTPUT\n", +"mprintf('Armature torque is %3.0f N.m \n Horse power output is %3.1f HP',Ta,Br)\n", +"\n", +"//=================================END OF PROGRAM==============================\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.1: To_find_generated_emf.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.1, Page 1.14\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"I=450;//Load current in A\n", +"V=250;//Terminal voltage in V\n", +"Ra=0.04;//Armature resistance in ohm\n", +"Rf=50;//Shunt field resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Ia=(V/Rf)+I;//Total current in A\n", +"Eg=(V+(Ia*Ra));//Generated emf in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Generated emf is %3.1f V',Eg)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.20: To_find_the_armature_torque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.20, Page 1.46\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"\n", +"//INPUT DATA\n", +"IL=50;//Load current in A\n", +"V=220;//Terminal voltage in V\n", +"Ra=0.3;//Armature resistance in ohm\n", +"Rsh=220;//Field resistance in ohm\n", +"N=1200;//Speed of the motor in rpm\n", +"\n", +"//CALCULATIONS\n", +"Ish=(V/Rsh);//Field current in A\n", +"Ia=(IL-Ish);//Armature current in A\n", +"Eb=(V-(Ia*Ra));//Back emf in V\n", +"Ta=(9.55*Eb*Ia)/N;//Armature torque in N.m\n", +"\n", +"//OUTPUT\n", +"mprintf('Armature torque is %3.0f N.m',Ta)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.21: To_determine_the_speed_and_electro_magnetic_torque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=220;//Terminal voltage in V\n", +"P=(10*746);//Rating of the motor in W\n", +"Iao=5;//No load armature current in A\n", +"No=1200;//No load speed in rpm\n", +"Ra=0.3;//Armature resistance in ohm\n", +"Ial=35;//Armature load current in A\n", +"\n", +"//CALCULATIONS\n", +"Nl=(No*((V-(Ial*Ra))/(V-(Iao*Ra))));//Speed at load in rpm\n", +"Ebo=218.5;//Back emf at no load in V\n", +"EbL=209.5;//Back emf at full load in V\n", +"Tao=(9.55*Ebo*Iao)/No;//No load torque in N.m\n", +"TaL=(9.55*EbL*Ial)/Nl;//Load torque in N.m\n", +"\n", +"//OUTPUT\n", +"mprintf('Load speed is %3.0f rpm \n Load torque is %3.2f N.m',Nl,TaL)\n", +"\n", +"//=================================END OF PROGRAM==============================\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.22: To_calculate_the_speed_of_the_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=220;//Terminal voltage in V\n", +"Io=4;//No load current in A\n", +"No=800;//No load speed in rpm\n", +"IL=24;//Load current in A\n", +"Ra=0.25;//Armature resistance in ohm\n", +"Rsh=220;//Shunt field resistance in ohm\n", +"No=800;//No load speed in rpm\n", +"\n", +"//CALCULATIONS\n", +"Ish=(V/Rsh);//Field current in A\n", +"Iao=Io-Ish;//Armature current at no load in A\n", +"IaL=IL-Ish;//Armature current at load in A\n", +"Nl=(No*((V-(IaL*Ra))/(V-(Iao*Ra))));//Speed at load in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Speed of the motor at load is %3.0f rpm',Nl)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.23: To_calculate_the_speed_of_a_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=6;//Number of poles\n", +"A=6;//Number of parallel paths for lap wound\n", +"Z=600;//Number of conductors\n", +"IL=100;//Load current in A\n", +"V=120;//Terminal voltage in V\n", +"Ra=30;//Armature resistance in ohm\n", +"Rsh=0.06;//Shunt field resistance in ohm\n", +"q=(30*10^-3);//Flux per pole in Wb\n", +"\n", +"//CALCULATIONS\n", +"Ish=(V/Ra);//Field current in A\n", +"Ia=(IL-Ish);//Armature current in A\n", +"Eb=(V-(Ia*Rsh));//Back emf in V\n", +"N=(60*Eb*A)/(q*Z*P);//Speed of the motor in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Speed of the lap wound shunt motor is %3.0f rpm',N)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.24: To_calculate_the_speed.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"//INPUT DATA\n", +"Pg=120000;//Power delivered when generator in W\n", +"Ng=1000;//Prime mover speed in rpm\n", +"Vg=600;//Terminal voltage given by the generator in V dc\n", +"Pm=120000;//Power taken as motor in W\n", +"Vm=600;//Terminal voltage when motor in V dc\n", +"Ra=0.05;//Armature resistance in ohm\n", +"Rsh=200;//Field resistance in ohm\n", +"Vb=1;//Brush drop in V\n", +"Ng=1000;//Speed of the generator in rpm\n", +"//CALCULATIONS\n", +"//When operated as a generator\n", +"IL1=(Pg/Vg);//Load current in A\n", +"If1=(Vg/Rsh);//Filed current in A\n", +"Ia1=(IL1+If1);//Armature current in A\n", +"Eg=(Vg+(Ia1*Ra)+Vb);//Generated emf in V\n", +"//When operated as a motor\n", +"IL2=(Pm/Vm);//Load current in A\n", +"If2=(Vm/Rsh);//Field current in A\n", +"Ia2=(IL2-If2);//Armature current in A\n", +"Eb=(Vm-(Ia2*Ra)-Vb);//Back emf in V\n", +"Nm=(Ng*Eb)/Eg;//Speed of the motor in rpm\n", +"//OUTPUT\n", +"mprintf('Speed of the dc machine when operated as a motor is %3.0f rpm',Nm)\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.25: To_find_the_speed_and_torque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"R=0.05;//Total resistance of the motor in ohm\n", +"IL1=120;//Load current in A\n", +"V=220;//Terminal voltage in V\n", +"N=1200;//Speed in rpm\n", +"IL2=60;//Half load current in A\n", +"\n", +"//CALCULATIONS\n", +"//Tnew=0.25*Told\n", +"//Hence percentage change in torque is 75% since it is (Told-Tnew)/Told*100\n", +"Ebnew=(V-(IL1*R));//New back emf in V\n", +"Ebold=(V-(IL2*R));//Old back emf in V\n", +"Nnew=(N*Ebnew*IL1)/(Ebold*IL2);//New speed in rpm\n", +"Pspeed=(Nnew/N)*100;//Percentage change in speed in %\n", +"//Ianew=(Iaold/sqrt(2))\n", +"I=sqrt(2)*100;//Percentage in current \n", +"N1new=(sqrt(2)*Ebnew*N)/Ebold;//New speed in rpm\n", +"P1speed=(N1new/N)*100;//Percentage change in speed in %\n", +"\n", +"//OUTPUT\n", +"mprintf('i)Percentage in speed is %3.2f and Percentage in torque is 75\nii)New speed is %3.0f rpm and new current is (1/sqrt(2)) times old current',Pspeed,N1new)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.26: To_find_the_efficiency_of_the_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=220;//Terminal voltage in V\n", +"ILo=5;//No load current in A\n", +"Ra=0.3;//Armature resistance in ohm\n", +"Rsh=220;//Field resistance in ohm\n", +"IL=50;//Load current in A\n", +"\n", +"//CALCULATIONS\n", +"Lo=(ILo*V);//No load losses in W\n", +"Ish=(V/Rsh);//Shunt current in A\n", +"Iao=(ILo-Ish);//No load armature current in A\n", +"Lco=((Iao^2*Ra)+(Ish^2*Rsh));//No load copper losses in W\n", +"Ifl=(Lo-Lco);//Iron and friction losses in W\n", +"Ia=(IL-Ish);//Armature current in A\n", +"Vl=(Ia^2*Ra);//Variable losses in W\n", +"Tl=(Vl+Lco+Ifl);//Total losses in W\n", +"P=(V*IL);//Input power in W\n", +"n=((P-Tl)/P)*100;//Efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('Efficiency of the motor is %3.1f percent',n)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.27: To_calculate_the_properties_of_shunt_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=250;//Terminal voltage in V\n", +"IL=50;//Load current in A\n", +"N=1000;//Speed in rpm\n", +"Wi=1200;//Iron and friction losses in W\n", +"Ra=0.05;//Armature resistance in ohm\n", +"Rsh=125;//Field resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Ish=(V/Rsh);//Field current in A\n", +"Ia=(IL-Ish);//Armature current in A\n", +"Eb=(V-(Ia*Ra));//Back emf in V\n", +"Cu=((V*IL)-(Eb*Ia));//Copper losses in W\n", +"Ta=(9.55*Eb*Ia)/N;//Armature torque in N.m\n", +"Ts=(9.55*((Eb*Ia)-Wi))/N;//Shaft torque in N.m\n", +"n=(((Eb*Ia)-Wi)/(V*IL))*100;//Efficiency of the motor\n", +"\n", +"//OPUTPUT\n", +"mprintf('(i)Copper loss is %3.1f W\n(ii)Armature torque is %3.1f N.m\n(iii)Shaft torque is %3.2f N.m\n(iv)Efficiency is %3.1f percent',Cu,Ta,Ts,n)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.28: To_find_the_speed_and_load_current_and_speed_regulation.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"A=2;//Number of parallel paths \n", +"Z=926;//Number of conductors\n", +"P=4;//Nmber of poles\n", +"V=220;//Line voltage in V\n", +"Io=3;//No load ine current in A\n", +"If=0.8;//No load field current in A\n", +"q=(6*10^-3);//No load field flux in Wb\n", +"Ra=0.9;//Armature resistance in ohm\n", +"T=30;//Load torque in N.m\n", +"\n", +"//CALCULATIONS\n", +"Ebo=(V-((Io-If)*Ra));//No load back emf in V\n", +"No=(Ebo*60*A)/(q*Z*P);//No load speed in rpm\n", +"Ia=(A*T)/(0.159*q*Z*P);//Armature current in A\n", +"IL=(Ia+If);//Load current in A\n", +"Eb=(V-(Ia*Ra));//Back emf in V\n", +"N=(Eb*60*A)/(q*Z*P);//Speed at load in rpm\n", +"R=((No-N)/No)*100;//Speed regulation in percent\n", +"\n", +"//OUTPUT\n", +"mprintf('No load speed is %3.0f rpm\nSpeed at load is %3.1f rpm\nSpeed regulation is %3.2f percent',No,N,R)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.29: To_find_the_change_in_speed.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=250;//Terminal voltage in V dc\n", +"N1=800;//Existing speed in rpm\n", +"Ra=0.05;//Armature resistance in ohm\n", +"Ia1=40;//Existing armature current in A\n", +"R=0.1;//Reduction in field flux\n", +"\n", +"//CALCULATIONS\n", +"Ia2=(Ia1/(1-R));//New armature current in A\n", +"N2=(N1*(V-(Ia1*Ra)))/((V-(Ia2*Ra))*(1-R));//New speed in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('New speed is %3.0f rpm',N2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.2: To_calculate_induced_emf_and_armature_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.2, Page 1.15\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Il=40;//Load current in A\n", +"V=400;//Terminal voltage in V\n", +"Ra=0.04;//Armature resistance in ohm\n", +"Rse=0.02;//Series field resistance in ohm\n", +"Rsh=300;//Shunt field resistance in ohm\n", +"V1=2;//Voltage drop across the brushes in V\n", +"\n", +"//CALCULATIONS\n", +" Ia=Il+(V/Rsh);//Armature current in A\n", +" Eg=V+(Ia*Ra)+(Ia*Rse)+V1;//Generated emf in V\n", +" \n", +"//OUTPUT\n", +"mprintf('Induced emf is %3.3f V \n Armature current is %3.2f A',Eg,Ia)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.30: To_find_the_resistance_to_be_included.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=300;//Terminal voltage in V\n", +"N1=600;//Existing speed in rpm\n", +"IL=30;//Load current in A\n", +"N2=800;//New speed in rpm\n", +"Ra=0.5;//Armature resistance in ohm\n", +"Rsh=125;//Field resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Ish1=(V/Rsh);//Field current in A\n", +"Ia1=(IL-Ish1);//Armature current in A\n", +"Ia2=(V-sqrt((V^2)-(4*Ra*(V-(Ia1*Ra))*Ia1*(N2/N1))));//New armature current in A\n", +"Ish2=(Ish1*Ia1)/Ia2;//New field current in A\n", +"Rsh2=(V/Ish2);//New field resistance in ohm\n", +"FR=(Rsh2-Rsh);//Field rheostat in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('The value of resistance to be included in the field is %3.2f ohm',FR)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.31: To_find_the_resistance_required.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=1000;//Initial speed in rpm\n", +"N2=600;//Final speed in rpm\n", +"Ia1=40;//Initial armature current in A\n", +"Ia2=30;//Final armature current in A\n", +"V=250;//Terminal voltage in V\n", +"Ra=0.5;//Armature resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"R=(V-((N2/N1)*(V-(Ia1*Ra))))/30;//Total resistance in ohm\n", +"Rs=(R-Ra);//Series resistance in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('Resistance required in series is %3.2f ohm',Rs)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.32: To_find_the_additional_field_resistance_to_be_included.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=1200;//Initial speed in rpm\n", +"N2=1500;//Final speed in rpm\n", +"Ia1=80;//Initial armature current in A\n", +"Ia2=100;//Final armature current in A\n", +"V=220;//Terminal voltage in V\n", +"Ra=0.05;//Armature resistance in ohm\n", +"Rsh1=220;//Initial shunt resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Rsh2=((N2/N1)*(V-(Ia1*Ra))*Rsh1)/(V-(Ia2*Ra));//New shunt resistance in ohm\n", +"Rs=(Rsh2-V);//Field resistance in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('Additional field resistance to be included in the field is %3.2f ohm',Rs)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.33: To_find_the_speed.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"clc\n", +"clear\n", +"//INPUT DATA\n", +"N1=1500;//Initial speed in rpm\n", +"N2=1200;//Final speed in rpm\n", +"Ia1=30;//Initial armature current in A\n", +"V=300;//Terminal voltage in V\n", +"Ra1=0.5;//Initial armature resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"R=(V-((N2/N1)*(V-(Ia1*Ra1))))/Ia1;//Total resistance in ohm\n", +"Rs=(R-Ra1);//Resistance to be added in ohm\n", +"n=((V-(Ia1*R))/V)*100;//Armature circuit efficiency\n", +"Nn2=(N2*(V-((Ia1/2)*R)))/(V-(Ia1*R));//New speed at half of the full load torque in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Resistance to be added to the existing armature resistance is %3.1f ohm \n Speed at half of the full load torque is %3.1f rpm',Rs,Nn2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.34: To_find_the_properties_of_shunt_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Pi=8800;//Input power in W\n", +"Ra=0.5;//Armature resistance in ohm\n", +"No=1260;//Speed of the motor at no load in rpm\n", +"V=240;//Line voltage in V\n", +"Pm=18800;//Gross mechanical power in W\n", +"V=240;//treminal voltage in V\n", +"\n", +"//CALCULATIONS\n", +"K=(V/No);//Constant of proportionality \n", +"Eb1=(240-sqrt((V^2)-(4*(Pi/2))))/2;//Back emf in V\n", +"Eb2=(240+sqrt((V^2)-(4*(Pi/2))))/2;//Back emf in V\n", +"I=(Pi/V);//Rated current in A\n", +"Ia=(V-Eb1)/Ra;//Armature current in A\n", +"Nn=(Eb2/K);//New speed in rpm\n", +"Ia2=(V-Eb2)/Ra;//Armature current in A\n", +"T=(60*Pi)/(2*3.14*Nn);//Torque developed in N.m\n", +"K2=(T/Ia2);//New constant\n", +"//TN=(0.5*10^-4)Nn^2\n", +"Nnn1=(-((K2*K)/Ra)+sqrt((((K2*K)/Ra))^2+(4*(0.5*10^-4)*((V*K2)/Ra))))/(2*0.5*10^-4);//New speed in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Armature current is %3.0f A \n Torque developed is %3.2f N.m \n New speed of the motor is %3.0f rpm',Ia2,T,Nnn1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.35: To_find_the_torque_and_power_and_speed.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=1500;//Initial speed in rpm\n", +"V1=270;//Terminal voltage in V\n", +"T=300;//Full load torque in N.m\n", +"N2=1200;//New speed in rpm\n", +"V2=(2*V1);//New terminal voltage in V\n", +"Ra=0.31;//Armature resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Ia=(T*2*3.14*N1)/(V1*60);//Full load current in A\n", +"Eb=(V1*(N2/N1));//Back emf in V\n", +"Pm=(Eb*Ia)/1000;//Mechanical power developed in kW\n", +"Eb2=(V2-(Ia*Ra));//Back emf at new terminal volatge in V\n", +"N=(Eb2*Ia*60)/(2*3.14*T);//New speed in rpm\n", +"Pm2=(Eb2*Ia/1000);//Mechanical power in kW\n", +"\n", +"//OUTPUT\n", +"mprintf('i)Full load current is %3.1f A, Full load power is %3.1f kW, Armature resistance is %3.2f ohm\nii)New motor torque is %3.0f N.m, Motor power is %3.1f kW, Motor speed is %3.0f rpm',Ia,Pm,Ra,T,Pm2,N)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.36: To_find_the_series_resistance_required.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=200;//Terminal voltage in V\n", +"Ra=0.05;//Armature resistance in ohm\n", +"Rse=0.03;//Field resistance in ohm\n", +"N1=1000;//Present speed in rpm\n", +"N2=800;//Required speed in rpm\n", +"Ia=40;//Armature current in A\n", +"\n", +"//CALCULATIONS\n", +"R=(V-((N2/N1)*(V-(Ia*(Ra+Rse)))))/Ia;//Total resistance in ohm\n", +"R1=(R-Ra-Rse);//Series resistance required to be connected in series with armature and field resistance in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('Series resistance required to be connected in series with armature and field resistance is %3.3f ohm',R1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.37: To_find_the_series_resistance_to_be_added.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=500;//Terminal voltage in V dc\n", +"I=30;//Line current in A\n", +"N1=600;//Initial speed in rpm\n", +"N2=500;//Required speed in rpm\n", +"R=0.5;//Total resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Eb1=(V-(I*R));//Back emf in V\n", +"Ka=(Eb1*60)/(I*N1);//Proportionality constant\n", +"T1=(Ka*I^2)/(2*3.14);//Torque developed at speed 600 rpm\n", +"T2=(T1*(N2/N1)^2);//Torque developed at speed 500 rpm\n", +"I2=sqrt((2*3.14*T2)/Ka);//New lin ecurrent in A\n", +"Eb2=(Ka*I2*N2)/60;//New back emf in V\n", +"R1=(V-Eb2)/I2;//Required series resistance in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('Series resistance to be added to armature field circuit is %3.1f ohm',R1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.38: To_calculate_the_Resistance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.21, Page 1.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=200;//Terminal voltage in V dc\n", +"I1=25;//Line current in A\n", +"Ra=0.5;//Armature resistance in ohm\n", +"Rse=0.3;//Field resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"//N2=0.75*N1\n", +"I2=sqrt((I1^2*(0.75)^3));//New line current in A\n", +"Eb1=(V-(I1*(Ra+Rse)));//Back emf in V\n", +"X=(V*I2);//X value for Resistance\n", +"R=(X-(0.75*Eb1))/I2^2;//Total resistance in ohm\n", +"Rs=(R-Ra-Rse);//Resistance to be connected in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('Resistance to be connected is %3.1f ohm',Rs)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.3: To_calculate_induced_emf_and_armature_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.3, Page 1.15\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V=250;//Terminal voltage in V\n", +"IL=40;//Load current in A\n", +"Ra=0.04;//Armature resistance in ohm\n", +"Rse=0.03;//Series field resistance in ohm\n", +"Rsh=100;//Shunt field resistance in ohm\n", +"Vbr=2;//Voltage drop across brushes in V\n", +"\n", +"//CALCULATIONS\n", +"Vsh=(V+(IL*Rse));//Voltage across shunt field in V\n", +"Ia=(IL+(Vsh/Rsh));//Armature current in A\n", +"Eg=(V+(IL*Rse)+(Ia*Ra)+Vbr);//Generated emf in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Induced emf is %3.1f V \nArmature current is %3.3f A',Eg,Ia)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.4: To_calculate_the_total_power_delivered.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.4, Page 1.16\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Ed=25000;//Power delivered by the generator in W\n", +"V=250;//Terminal voltage in V\n", +"Rsh=75;//Shunt field resistance in ohm\n", +"Ra=0.03;//Armature resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"IL=(Ed/V);//Load current in A\n", +"If=(V/Rsh);//Field current in A\n", +"Ia=(IL+If);//Armature current in A\n", +"Eg=(V+(Ia*Ra));//Generated emf in V\n", +"Pg=(Eg*Ia)/1000;//Generated power in kW\n", +"\n", +"//OUTPUT\n", +"mprintf('Total power delivered by the armature is %3.2f kW',Pg)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.5: To_calculate_the_emf_generated.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.5, Page 1.17\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"n=48;//Number of slots\n", +"z=16;//Number of conductors per slot\n", +"q=0.018;//Flux per pole in Wb\n", +"P=4;//Number of poles\n", +"N=1000;//Speed of armature in rpm\n", +"A=2;//Number of parallel paths \n", +"\n", +"//CALCULATIONS\n", +"Z=(n*z);//Number of conductors\n", +"Eg=(q*Z*N*P)/(60*A);//Generated emf in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Generated emf is %3.1f V',Eg)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.6: To_calculate_the_speed.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.6, Page 1.18\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=4;//Number of poles\n", +"Z=400;//Number of conductors\n", +"q=0.03;//Flux per pole in Wb\n", +"Eg=250;//Generated emf in V\n", +"A1=4;//Number of parallel paths in lap wound\n", +"A2=2;//Number of parallel paths in wave wound\n", +"\n", +"//CALCULATIONS\n", +"N1=(60*Eg*A1)/(q*Z*P);//Speed reqired in lap wound in rpm\n", +"N2=(60*Eg*A2)/(q*Z*P);//Speed reqired in wave wound in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Speed reqired in lap wound is %3.0f rpm \nSpeed reqired in wave wound is %3.0f rpm',N1,N2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.7: To_calculate_the_flux_per_pole.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.7, Page 1.18\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Eg1=250;//Existing generated emf in V\n", +"N1=800;//Existed rated speed in rpm\n", +"q1=0.03;//Existing flux in Wb\n", +"Eg2=300;//New generated emf in V\n", +"N2=1000;//New rated speed in rpm\n", +"\n", +"//CALCULATIONS\n", +"q2=(q1*N1*Eg2)/(Eg1*N2);//New flux per pole in Wb\n", +"\n", +"//OUTPUT\n", +"mprintf('New flux per pole is %3.4f Wb',q2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.8: To_calculate_the_terminal_voltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-1, Example 1.8, Page 1.19\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//NPUT DATA\n", +"n=200;//Number of turns\n", +"P=6;//Number of poles\n", +"A=P;//Since lap wound turns\n", +"Ra=0.0112;//Armature resistance in ohm\n", +"Ia=40;//Armature current in A\n", +"N=1000;//Armature speed in rpm\n", +"q=0.03;//Flux per pole in Wb\n", +"\n", +"//CALCULATIONS\n", +"Z=(n*2);//Total number of conductors\n", +"Eg=(q*Z*N*P)/(60*A);//Generated emf in V\n", +"IaRa=(Ia*Ra);//Armature drop in VI\n", +"V=(Eg-IaRa);//Terminal voltage in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Terminal voltage is %3.3f V',V)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.9: To_calculate_the_speed.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-9, Example 1.9, Page 1.20\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=4;//Number of poles\n", +"A=2;//Number of parallel paths for wave wound\n", +"Z=400;//Number of conductors\n", +"q=(20*10^-3);//Flux per pole in Wb\n", +"Ra=0.04;//Armature resistance in ohm\n", +"Rsh=75;//Shunt field resistance in ohm\n", +"V=250;//Terminal voltage in V\n", +"PL=(600*100);//Total load on the generator in W\n", +"Vld=10;//Line drop in V\n", +"\n", +"//CALCULATIONS\n", +"IL=(PL/V);//Load current in A\n", +"Ish=(V/Rsh);//Shunt field current in A\n", +"Ia=(IL+Ish);//Armature current in A\n", +"Eg=(V+(Ia*Ra)+Vld);//Generated emf in V\n", +"N=(60*Eg*A)/(q*Z*P);//Speed at which the generator should be driven in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Speed at which the generator should be driven is %i rpm',N)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Electrical_Engineering_by_R_Anandanatarajan/2-Transformers.ipynb b/Electrical_Engineering_by_R_Anandanatarajan/2-Transformers.ipynb new file mode 100644 index 0000000..3031c63 --- /dev/null +++ b/Electrical_Engineering_by_R_Anandanatarajan/2-Transformers.ipynb @@ -0,0 +1,1796 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 2: Transformers" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.10: Primary_current_and_power_factor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.10, Page 2.17\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Io=12;//Primary no load current in A\n", +"cosqo=0.25;//No load power factor \n", +"I2=220;//Secondary load current in A\n", +"cosq2=0.8;//Secondary power factor\n", +"K=(1/5);//Turn ratio\n", +"\n", +"//CALCULATIONS\n", +"qo=acosd(cosqo);//phase angle in degree\n", +"q2=acosd(cosq2);//Phase angle in degree\n", +"Ioc=complex((Io*cosqo),(Io*sind(qo)));//Io value in complex form\n", +"I2i=complex((K*I2*cosq2),(K*I2*sind(q2)));//I2i value in complex form\n", +"I1=(Ioc+I2i);//Primary current in A\n", +"X=sqrt((real(I1))^2+(imag(I1))^2);//Primary current in A\n", +"Y=cosd(atand(imag(I1)/real(I1)));//Power factor\n", +"\n", +"//OUTPUT\n", +"mprintf('Primary current is %3.1f A and Primary power factor is %3.2f',X,Y)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.11: Properties_of_an_ideal_step_up_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.11, Page 2.18\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"K=(330/110);//Turn ratio\n", +"N1=110;//Number of turns in the primary\n", +"N2=330;//Number of turns in the secandary\n", +"V1=4000;//Primary voltage in V\n", +"f=50;//Supply frquency in Hz\n", +"Z2=complex(120,40);//Secondary load\n", +"\n", +"//CALCULATIONS\n", +"q=(V1/(4.44*N1*f));//Flux in Wb\n", +"V2=(K*V1);//Secondary voltage in V\n", +"I2=(V2/Z2);//Secondary current in A\n", +"I1=K*I2;//Primary current in A\n", +"S=(V1*I1)/1000;//Transformer rating\n", +"P1=(V1*sqrt((real(I1))^2+(imag(I1))^2)*cosd((atand(imag(I1)/real(I1)))))/1000;//Real power in kW\n", +"R1=(V1*sqrt((real(I1))^2+(imag(I1))^2)*sind(-(atand(imag(I1)/real(I1)))))/1000;//Reactive power in KVAR\n", +"Zeq=(V1/I1);//Transformer equivalent impedence \n", +"a1=sqrt((real(I1))^2+(imag(I1))^2);\n", +"a2=sqrt((real(I2))^2+(imag(I2))^2);\n", +"b1=real(Zeq);\n", +"b2=imag(Zeq);\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Maximum flux in the core is %3.3f Wb\n(b)Primary current is %3.2f A and Secondary current is %3.2f A\n(c)Real power is %3.0f KW and Reactive power is %3.0f KVAR\n(d)Value of impedence consumed is %3.1f+j%3.1f',q,a1,a2,P1,R1,b1,b2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.12: Primary_current_and_peak_value_of_flux.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.12, Page 2.20\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=100;//Number of turns in the primary\n", +"N2=2000;//Number of turns in the secondary\n", +"V1=220;//Primary volatge in V\n", +"f=50;//Supply frequency in Hz\n", +"I2=3;//Secondary current in A\n", +"\n", +"//CALCULATIONS\n", +"K=(N2/N1);//Turn ratio\n", +"I1=(K*I2);//Primary current in A\n", +"q=(V1/(4.44*f*N1))*1000;//Peak vaue of flux linked with the secondary in m.Wb\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)The value of primary current is %3.0f A \n(b)The peak value of flux linked with the secondary is %3.1f m.Wb',I1,q)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.13: Secondary_voltage_and_primary_and_secondary_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.13, Page 2.21\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=1100;//Number of turns in the primary\n", +"N2=550;//Number of turns in the secandary\n", +"V1=200;//Primary voltage in V\n", +"R2=5;//Resistance in the secondary in ohm\n", +"\n", +"//CALCULATIONS\n", +"K=(N2/N1);//Turn ratio\n", +"V2=(K*V1);//Secondary voltage in V\n", +"I2=(V2/R2);//Current in the secondary in A\n", +"I1=(K*I2);//Current in the primary in A\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)Secondary voltage is %3.0f V\n(b)Primary current is %3.0f A\n(c)Secondary current is %3.0f A',V2,I2,I1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.14: Primary_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.16, Page 2.30\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=400;//Number of turns in the primary\n", +"N2=100;//Number of turns in the secondary\n", +"Io=4;//No load current in A\n", +"qo=0.3;//No load current power factor \n", +"I2=120;//Secondary current in A\n", +"q2=0.8;//Secondary current power factor\n", +"\n", +"//CALCULATIONS\n", +"K=(N2/N1);//Turn ratio\n", +"I2i=(K*I2);//Secondary current in A\n", +"I2ic=complex((I2i*q2),(I2i*sind(acosd(q2))));//Ixi in complex form\n", +"Ioc=complex((Io*qo),(Io*sind(acosd(qo))));//Ixi in complex form\n", +"Iic=(I2ic+Ioc);//Primary current in complex form\n", +"a1=sqrt((real(Iic))^2+(imag(Iic))^2);\n", +"a2=atand(imag(Iic)/real(Iic));\n", +"q=cosd(a2);//Phase angle in degree\n", +"\n", +"//OUTPUT\n", +"mprintf('Current taken by the primary is %3.2f A and power factor is %3.2f',a1,q)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.15: Power_delivered_and_current_taken.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.15, Page 2.23\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=6000;//Primary volatge in V\n", +"V2=500;//Secondary voltage in V\n", +"Z2=complex(4,3)\n", +"\n", +"//CALCULATIONS\n", +"K=(V2/V1);//Voltage transformation ratio\n", +"I2=(V2/Z2);//Secondary current in A\n", +"a1=sqrt((real(I2))^2+(imag(I2))^2);\n", +"a2=atand(imag(I2)/real(I2));\n", +"q=cosd(a2);//Phase angle in degree\n", +"P2=(V2*a1*q)/1000;//Power delivered in kW\n", +"I1=(K*a1);//Primary current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('Power delivered is %3.0f kW \nCurrent taken by an ideal transformer is %3.2f A',P2,I1)\n", +" \n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.16: Parameters_of_a_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.1, Page 2.4\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Q=25;//Rating of a transformer in KVA\n", +"V1=2000;//Primary voltage in V\n", +"V2=200;//Secondary volatge in V\n", +"R1=0.15;//Primary winding resistance in ohm\n", +"X1=0.25;//Primary leakage reactance in ohm\n", +"R2=0.04;//Secondary winding resistance in ohm\n", +"X2=0.015;//Secondary leakage reactance in ohm\n", +"\n", +"//CALCULATIONS\n", +"K=(V2/V1);//Voltage transformation ratio\n", +"Ro1=(R1+(R2/K^2));//Resistance referred to primary in ohm\n", +"Xo1=(X1+(X2/K^2));//Reactance referred to primary in ohm\n", +"Zo1=sqrt(Ro1^2+Xo1^2);//Impedence referred to primary in ohm\n", +"Ro2=(R2+(R1*K^2));//Resistance referred to secndary in ohm\n", +"Xo2=(X2+(X1*K^2));//Reactance referred to secondary in ohm\n", +"Zo2=sqrt(Ro2^2+Xo2^2);//Impedence referred to secondary in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)Resistance referred to primary is %3.2f ohm \n Reactance referred to primary is %3.2f ohm \n Impedence referred to primary is %3.1f ohm \n\n(b)Resistance referred to secndary is %3.4f ohm \n Reactance referred to secondary is %3.4f ohm \n Impedence referred to secondary is %3.3f ohm',Ro1,Xo1,Zo1,Ro2,Xo2,Zo2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.17: Parameters_of_a_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.17, Page 2.32\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"R1=3.5;//Primary Resistance in ohm\n", +"X1=5.2;//Primary reactance in ohm\n", +"R2=0.01;//Secondary Resistance in ohm\n", +"X2=0.02;//Secondary reactance in ohm\n", +"Q=40000;//Rating of the transformer in VA\n", +"V1=4000;//Primary voltage in V\n", +"V2=200;//Secondary voltage in V\n", +"\n", +"//CALCULATIONS\n", +"Z1=complex(R1,X1);//Primary impedence\n", +"Z2=complex(R2,X2);//Secondary impedence\n", +"I1=(Q/V1);//Primary current in A\n", +"I2=(Q/V2);//Secondary current in A\n", +"K=(I1/I2);//Current ratio\n", +"Ro1=(R1+(R2/K^2));//Resistance referred to primary in ohm\n", +"Xo1=(X1+(X2/K^2));//Reactance referred to primary in ohm\n", +"Zo1=(Z1+(Z2/K^2));//Impedence referred to primary in ohm\n", +"Ro2=(R2+(R1*K^2));//Resistance referred to secndary in ohm\n", +"Xo2=(X2+(X1*K^2));//Reactance referred to secondary in ohm\n", +"Zo2=(Z2+(Z1*K^2));//Impedence referred to secondary in ohm\n", +"a1=real(Zo1);\n", +"a2=imag(Zo1);\n", +"a3=real(Zo2);\n", +"a4=imag(Zo2);\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)Resistance referred to primary is %3.1f ohm \n Reactance referred to primary is %3.1f ohm \n Impedence referred to primary is %3.1f+j%3.1f ohm \n\n(b)Resistance referred to secndary is %3.5f ohm \n Reactance referred to secondary is %3.3f ohm \n Impedence referred to secondary is %3.5f+j%3.3f ohm',Ro1,Xo1,a1,a2,Ro2,Xo2,a3,a4)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.18: Parameters_of_a_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.18, Page 2.34\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Q=(40*1000);//Transformer rating in VA\n", +"V1=1600;//Primary voltage in V\n", +"V2=160;//Secondary voltage in V\n", +"f=50;//Frequency in Hz\n", +"R=10;//Turn ratio\n", +"\n", +"//CALCULATIONS\n", +"K=0.1;//Turn ratio\n", +"I2=(Q/V2);//Full load secondary current in A\n", +"Z2=(V2/I2);//Load impedence in ohm\n", +"Zo1=(Z2/K^2);//Impedence referred to high tension side in ohm\n", +"I2i=(K*I2);//Value of current referred to high tension side in A\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)Load impedence required for full load current is %3.2f ohm \n(b)Impedence referred to high tension side is %3.0f ohm\n(c)Value of current referred to high tension side is %3.0f A',Z2,Zo1,I2i)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.19: Priamary_and_secondary_resistance_and_reactance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.19, Page 2.35\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Q=80;//Transformer rating in KVA\n", +"V1=11000;//Primart voltage in V\n", +"V2=440;//Secondary voltage in V\n", +"Pcu=0.75;//Primary copper loss in kW\n", +"Scu=0.5;//Secondary copper loss in kW\n", +"\n", +"//CALCULATIONS\n", +"I2=(Q*1000)/V2;//Full load secondary current in A\n", +"I1=(Q*1000)/V1;//Full load primary current in A\n", +"R1=((Pcu)/I1^2)*1000;//Primary resistance in ohm\n", +"R2=(Scu*1000)/I2^2;//Secondary resistance in ohm\n", +"Xo1=(0.04*V1)/I1;// Equivalent primary reactance in ohm\n", +"K=(I1/I2);//Current ratio\n", +"X1=(Xo1/((R1+(R2/K))/(R2/K)));//Primary reactance in ohm\n", +"X2i=(Xo1/X1);//Equivalent secondary reactance in ohm\n", +"X2=(X2i*K);//Secondary reactance in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Primary resistance is %3.2f ohm and Secondary resistance in is %3.3f ohm\nb)Primary reactance is %3.2f ohm and Secondary reactance is %3.3f ohm',R1,R2,X1,X2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.1: Number_of_turns_and_full_load_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.1, Page 2.4\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=1000;//Voltage in primary circuit in V\n", +"V2=100;//Voltage in secondary ciricut in V\n", +"N2=60;//Number of turns in secondary \n", +"R=10000;//Rating of transformer in VA\n", +"\n", +"//CALCULATIONS\n", +"K=(V2/V1);//Voltage transformation ratio\n", +"N1=(N2/K);//Number of turns in primary\n", +"I1=(R/V1);//Current in the primary in A\n", +"I2=(R/V2);//Current in the secondary in A\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Number of turns in the primary is %3.0f turns \nb)Current in the primary is %3.0f A and Current in the secondary is %3.0f A',N1,I1,I2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.20: Primary_voltage_and_power_factor_and_efficiency.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.20, Page 2.37\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"K=(1/20);//Turn ratio\n", +"R1=30;//Primary resistance in ohm\n", +"R2=0.08;//Secondary resistance in ohm\n", +"X1=80;//Primary reactance in ohm\n", +"X2=0.3;//Secondary reactance in ohm\n", +"I=1.5;//No load current in A\n", +"cosqo=0.5;//Power factor\n", +"I2=200;//Load current in A\n", +"V2=500;//Secondary terminal voltage in V\n", +"cosq2=0.8;//Load power factor\n", +"q3=60;//Phase angle in degree\n", +"\n", +"//CALCULATIONS\n", +"q2=(acosd(cosq2));//Phase angle in degree\n", +"I2i=complex((I2*cosd(q2)),(I2*sind(-q2)));//Load current in complex form\n", +"V2i=complex(V2,0);//Secondary terminal voltage in complex form\n", +"Z2=complex(R2,X2);//Impedence in complex form\n", +"E2=(V2i+(I2i*Z2));//Terminal voltage in V\n", +"E1=(sqrt((real(E2))^2+(imag(E2))^2)/K);//Primary voltage in V\n", +"I2c=(K*I2);//Secondary current in A\n", +"I21c=complex((I2c*cosd(q2)),(I2c*sind(-q2)));//Load current in complex form\n", +"Io=complex((I*cosd(-q3)),(I*sind(-q3)));//No load current in A\n", +"I1c=(Io+I21c);//Total current\n", +"Z1=complex(R1,X1);//Primary impedence\n", +"V1=(E1+(I1c*Z1));//Primary applied voltage \n", +"V1i=(sqrt((real(V1))^2+(imag(V1))^2));//Primary applied voltage in V\n", +"A=((atand(imag(V1)/real(V1)))-((atand(imag(I1c)/real(I1c)))));//Angle between V1 and I1 in degree\n", +"p=cosd(A);//Power factor\n", +"Cu=(I2^2*(R2+(K^2*R1)));//Copper losses in W\n", +"C=(V1i*sqrt((real(Io))^2+(imag(Io))^2)*cosqo);//Constant losses in W\n", +"P=(V2*I2*cosq2);//Output power in W\n", +"n=(P/(P+Cu+C))*100;//Efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('Primary applied voltage is %3.2f V\nPrimary power factor is %3.2f \nEfficiency is %3.2f percent',V1i,p,n);\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.21: Primary_induced_emf_and_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.21, Page 2.39\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=1000;//Primary voltage in V\n", +"V2=300;//Secondary voltage in V\n", +"R1=0.2;//Primary resistance in ohm\n", +"X1=0.75;//Primary reactance in ohm\n", +"I1=50;//Primary current in A\n", +"cosq1=0.8;//Power factor\n", +"\n", +"//CALCULATIONS\n", +"E1=(V1-(I1*sqrt(R1^2+X1^2)));//Primary induced emf in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Primary induced emf is %3.1f V',E1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.22: Induced_emf_in_the_secondary.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.22, Page 2.40\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"L2=7500;//Load on secondary in W\n", +"V2=220;//Secondary voltage in V\n", +"cosq=0.8;//Power factor\n", +"R2=0.05;//Secondary resistance in ohm\n", +"X2=0.75;//Secondary reactance in ohm\n", +"V2i=200;//Secondary voltage in V\n", +"\n", +"//CALCULATIONS\n", +"I2=(L2/(V2*cosq));//Secondary current in A\n", +"q=acosd(cosq);//Phase angle in degree\n", +"I2c=complex((I2*cosd(q)),(I2*sind(-q)));//I2 in complex form\n", +"Z2=complex(R2,X2);\n", +"E2=(V2i+(I2c*Z2));//Induced imf in V\n", +"a1=real(E2);\n", +"a2=imag(E2);\n", +"\n", +"//OUTPUT\n", +"mprintf('Induced emf in the secondary is %3.2f+j%3.2f',a1,a2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.23: Primary_currrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.23, Page 2.40\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"K=(1000/200);//Voltage transformation ratio\n", +"R1=2;//Primary resistance in ohm\n", +"R2=200;//Secondary resistance in ohm\n", +"Vo=360;//Volts in V\n", +"\n", +"//CALCULATIONS\n", +"Z2i=(R2/K^2);//Equivalent secondary impedence in ohm\n", +"Zo1=(Z2i+R1);//Equivalent primary impedence in ohm\n", +"I1=(Vo/Zo1);//Primary current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('Primary current is %3.0f A',I1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.24: Secondary_voltage_and_primary_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.24, Page 2.41\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"K=(500/10);//Turn ratio in step up transformer \n", +"Z1=complex(0,6);//Primary reactance in ohm\n", +"Z2=complex(20000,-10000);//Secondary impedence in ohm\n", +"V1=100;//Primary voltage in V\n", +"\n", +"//CALCULATIONS\n", +"Z2i=(Z2/K^2);//Equivalent secondary impedence in ohm\n", +"Zo1=(Z1+Z2i);//Equivalent primary impedence in ohm\n", +"I1=(V1/Zo1);//Primary current in A\n", +"V2i=(I1*Z2i);//Equivalent secondary voltage in V\n", +"V2=(K*V2i);//Secondary voltage in V\n", +"X=sqrt((real(V2))^2+(imag(V2))^2);//X value for secondary voltage\n", +"Y=-(45+atand(imag(V2)/real(V2)));//Phase angle in degree\n", +"\n", +"//OUTPUT\n", +"mprintf('Secondary volatge is %3.0f V,%3.1f degree',X,Y)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.25: Efficiency_and_regulation.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.25, Page 2.44\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=2200;//Primary volatge in V\n", +"V2=220;//Secondary voltage in V\n", +"K=(V2/V1);//Voltage transformation ratio\n", +"R1=0.3;//Primary resistance in ohm\n", +"R2i=0.24;//Equivalent secondary resistance in ohm\n", +"Ro=300;//No load resistance in ohm\n", +"RL=0.4;//Load resistance in ohm\n", +"X1=0.8;//Primary reactance in ohm\n", +"X2i=0.9;//Equivalent secondary reactance in ohm\n", +"Xo=1100;//No load reactance in ohm\n", +"XL=0.3;//Load reactance in ohm\n", +"\n", +"//CALCULATIONS\n", +"ZLi=(complex(RL,R1)/K^2);//Equivalent load impedence in ohm\n", +"Z1=complex(R1,X1);//Primary impedence \n", +"Z2i=complex(R2i,X2i);//Equivalent secondary impedence\n", +"Zo=complex(Ro,Xo);//No load impedence \n", +"Zeq=((Zo*(Z1+Z2i+ZLi))/(Zo+Z1+Z2i+ZLi));//Equivalent impedence\n", +"I1=(V1/Zeq);//Primary current in A\n", +"I2i=((I1*Zo)/(Zo+Z1+Z2i+ZLi));//Equivalent secondary current in A\n", +"Io=((I1*((Z1+Z2i+ZLi)/(Zo+Z1+Z2i+ZLi))));//No load current in A\n", +"Pi=(V1*sqrt((real(I1))^2+(imag(I1))^2)*cosd(atand(imag(I1)/real(I1))))/1000;//Input power in kW\n", +"Pcu1=(((real(I1))^2+(imag(I1))^2)*R1);//Primary copper losses in W\n", +"Pcu2=(((real(I1))^2+(imag(I1))^2)*R2i);//Primary copper losses in W\n", +"C=(((real(Io))^2+(imag(Io))^2)*Ro);//Constant losses in W\n", +"n=(((Pi*1000)-Pcu2-C)/(Pi*1000))*100;//Efficiency\n", +"R=((V1-(sqrt((real(I2i))^2+(imag(I2i))^2)*sqrt((real(ZLi))^2+(imag(ZLi))^2)))/((sqrt((real(I2i))^2+(imag(I2i))^2)*sqrt((real(ZLi))^2+(imag(ZLi))^2))))*100;//Percentage Regulation\n", +"\n", +"//OUTPUT\n", +"mprintf('Efficiency is %3.1f percent \nRegulation is %3.2f percent',n,R)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.26: Secondary_terminal_voltage_at_full_load.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.26, Page 2.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"R1=6;//Primary resistance in ohm\n", +"R2=0.3;//Secondary resistance in ohm\n", +"X1=10;//Primary reactance in ohm\n", +"X2=0.5;//Secondary reactance in ohm\n", +"E1=2220;//primary induced emf in V\n", +"E2=220;//Secondary induced resistance in V\n", +"V1=2220;//Primary voltage drop in V\n", +"R=8;//Rate of transformer in KVA\n", +"K=E2/E1;//Transformer voltage ratio\n", +"cosQ=0.8;//Power factor\n", +"sinQ=0.6;//sine of Q\n", +"\n", +"//CALCULATIONS\n", +"R02=R2+(K^2*R1);//Resistance refered to the secondary in ohms\n", +"X02=X2+(K^2*X1);//Reactance refered to the secondary in ohms\n", +"I2=((R*1000)/E2);//Secondary full load current in A\n", +"V02=(I2*R02*cosQ)+(I2*X02*sinQ);//Secondary voltage drop in V\n", +"V2=E2-V02;//Secondary terminal voltage in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Secondary terminal voltage at full load is %3.1f V',V2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.27: Voltage_regulatio.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.27, Page 2.49\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Tr=36;//Transformer rating in terms of KVA\n", +"E1=5000;//Primary induced emf in V\n", +"E2=500;//Secondary induced emf in V\n", +"R01=22;//Winding resistance referred to the primary in ohm\n", +"X01=36;//Winding reactance referred to primary in ohm\n", +"cosQ1=0.8;//Primary power factor\n", +"cosQ2=0.8;//Secondary power factor\n", +"sinQ1=0.6;//sine of Q1\n", +"\n", +"//CALCULATIONS\n", +"I1=((X01*1000)/E1);//Full load primary current in A\n", +"Vd=(I1*R01*cosQ1)+(I1*X01*sinQ1);//Secondary voltage drop in V\n", +"V=(Vd/E1)*100;//Percentage voltage regulation in %\n", +"\n", +"//OUTPUT\n", +"mprintf('Percentage voltage regulation is %3.1f percent',V)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.28: Voltage_regulatio.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.28, Page 2.51\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Rp=1;//Percentage resistance drop in percentage\n", +"Xp=4;//Percentage reactance drop in percentage\n", +"cosQ1=0.8;//Lagging power factor\n", +"sinQ1=0.6;//Sine of Q1\n", +"cosQ2=1;//Power factor\n", +"sinQ2=0;//Sine of Q2\n", +"cosQ3=0.8;//Leading power factor\n", +"sinQ3=0.6;//Sine of Q3\n", +"\n", +"//CALCULATIONS\n", +"Vla=(Rp*cosQ1)+(Xp*sinQ1);//Percentage secondary voltage drop for lagging power factor in percentage\n", +"V=(Rp*cosQ2)+(Xp*sinQ2);//Percentage secondary voltage drop for unity power factor in percentage\n", +"Vle=(Rp*cosQ3)-(Xp*sinQ3);//Percentage secondary voltage drop for leading power factor in percentage\n", +"\n", +"//OUTPUT\n", +"mprintf('Secondary voltage drop for lagging power factor is %3.1f percent\nSecondary voltage drop for unity power factor is %3.1f percent\nsecondary voltage drop for leading power factor is %3.1f percent',Vla,V,Vle )\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.29: Power_factor_and_regulation.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.29, Page 2.52\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Resistance=3;//% Resistance drop\n", +"Reactance=6;//% Reactance drop\n", +"\n", +"//CALCULATIONS\n", +"q=atand(Reactance/Resistance);//Phase angle in degree\n", +"cosq=cosd(q);//Power factor\n", +"Regulation=((Resistance*cosq)+(Reactance*sind(q)));//% Regulation at the power factor\n", +"\n", +"//OUTPUT\n", +"mprintf('Power factor is %3.2f \nPercentage regulation at this power factor is %3.1f percent',cosq,Regulation)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.2: Properties_of_an_ideal_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.2, Page 2.5\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=400;//Number of turns in the primary\n", +"N2=30;//Number of turns in the secondary\n", +"Q=20000;//Rating of the transformer in VA\n", +"V1=2000;//Primary voltage in V\n", +"f=50;//Power supply frequency in Hz\n", +"\n", +"//CALCULATIONS\n", +"K=(N2/N1);//Voltage transformation ratio\n", +"I1=(Q/V1);//Current in the primary in A\n", +"I2=(I1/K);//Current in the secondary in A\n", +"V2=(K*V1);//Secondary voltage in V\n", +"q=(V1/(4.44*f*N1));//Maximum flux in the core in Wb\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)Full load primary current is %3.0f A and secondary current is %3.2f A \n(b)Induced emf in the secondary is %3.0f V \n(c)Maximum flux in the core is %3.3f Wb',I1,I2,V2,q)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.30: Parameters_of_a_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.30, Page 2.57\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=250;//Primary voltage in V\n", +"V2=100;//Secondary voltage in V\n", +"I1=0.4;//Primary current in A\n", +"Wo=36;//No load power input in W\n", +"\n", +"//CALCULATIONS\n", +"K=(V2/V1);//Voltage transformation ratio \n", +"q=acosd(Wo/(V1*I1));//Phase angle in degree\n", +"Im=(I1*sind(q));//Magnetising current in A\n", +"Iw=(I1*cosd(q));//Working current in A\n", +"I=(I1*V1*cosd(q));//Iron loss in W\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)Turns ratio is %3.1f \n(b)Magnetising current is %3.3f A \n(c)Working current is %3.3f A \n(d)Iron loss is %3.0f W',K,Im,Iw,I)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.31: Primary_voltage_and_power_factor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.31, Page 2.58\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"I2=400;//Full load secondary current in A\n", +"I1=(I2*0.2);//Full load secondary current in A\n", +"K=(I1/I2);//Turns ratio\n", +"Z1=complex(0.5,1.5);//Transformer parameter\n", +"Z2=complex(0.02,0.05);//Transformer parameter\n", +"\n", +"//CALCULATIONS\n", +"Zo1=Z1+(Z2/K^2);//Transformer parameter\n", +"Vsc=(I1*Zo1);//Primary voltage under short circuit test in V\n", +"[A B]=polar(Vsc);//Primary voltage under short circuit test in V in polar form\n", +"B=atand(imag(Zo1)/real(Zo1));//Phase angle in degree\n", +"\n", +"//OUTPUT\n", +"mprintf('Primary voltage under short circuit test is %3.1f and %3.2f degree V (polar form)',A,B)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.32: Equivalent_resistance_and_leakage_reactance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.32, Page 2.58\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Q=250*1000;//Rating of a transformer in VA\n", +"V1=11000;//Rated primary voltage in V\n", +"V2=2200;//Rated secondary voltage in V\n", +"N1=1000;//Number of turns in the primary\n", +"N2=200;//Number of turns in the secondary\n", +"R1=1.5;//Primary resistance in ohm\n", +"R2=0.05;//Secondary resistance in ohm\n", +"Vsc=600;//Primary voltage when secondary is short circuited in V\n", +"n=0.99;//Efficiency of the transformer\n", +"\n", +"//CALCULATIONS\n", +"K=(N2/N1);//Turn ratio\n", +"I1=(Q/(V1*n));//Full load primary current in A\n", +"Zo1=(Vsc/I1);//Equivalent reactance in ohm\n", +"R2i=(R2/K^2);//Equivalent secondary resistance in ohm\n", +"Ro1=(R1+R2i);//Equivalent primary resistance in ohm\n", +"Xo1=sqrt(Zo1^2-Ro1^2);//Equivalent ractance in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('Equivalent resistance referred to primary is %3.2f ohm \nEquivalent reactance referred to primary is %3.2f ohm',Ro1,Xo1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.33: Efficiency_of_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.33, Page 2.64\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"L=400;//Constant or Iron losses in W\n", +"C=700;//Full load copper loss in W\n", +"Q=40000;//Rating of transformer in VA\n", +"cosq=0.85;//Load power factor\n", +"\n", +"//CALCULATIONS\n", +"P=(Q*cosq);//Full load output in W\n", +"LC=(L+C);//Total full load losses in W\n", +"IP=(P+LC);//Full load input in W\n", +"n=(P/IP)*100;//Full load efficiency\n", +"P2=(0.5*Q*cosq);//Half load output in W\n", +"LC2=(L+(0.5^2*C));//Total losses at half loads in W\n", +"IP2=(P2+LC2);//Half load input in W\n", +"n2=(P2/IP2)*100;//Half load efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('Efficiency of the transformer at full load is %3.2f percent \nEfficiency of the transformer at half load is %3.2f percent',n,n2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.34: Parameters_of_a_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.34, Page 2.65\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Q=50000;//Rating of the transformer in VA\n", +"Pi=500;//Constat losses in W\n", +"Pcu=900;//Full load variable losses in W\n", +"cosq=0.8;//Power factor\n", +"\n", +"//CALCULATIONS\n", +"nFL=((Q*cosq)/((Q*cosq)+Pi+Pcu))*100;//Full load efficiency\n", +"L=(Q*sqrt(Pi/Pcu))/1000;//Load at which transformer operates at maximum efficiency in KVA\n", +"n=((L*1000)/((L*1000)+Pi+Pi))*100;//Maximum efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Full load efficiency is %3.2f percent \nb)Load at which transformer operates at maximum efficiency is %3.2f KVA \nc)Maximum efficiency is %3.2f percent',nFL,L,n)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.35: Efficiency_at_full_load.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.35, Page 2.66\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=5000;//Primary voltage in V\n", +"V2=200;//Secondary voltage in V\n", +"Q=60000;//Rating of transformer in VA\n", +"R1=8;//Primary resistance in ohm\n", +"R2=0.009;//Secondary resistance in ohm\n", +"Io=0.4;//No load primary current in A\n", +"cosq=0.29;//Power factor\n", +"\n", +"//CALCULATIONS\n", +"K=(V2/V1);//Turn ratio\n", +"Cu=((Q/V1)^2*(R1+(R2/K^2)));//Full load copper losses in W\n", +"C=(V1*Io*cosq);//Constant losses in W\n", +"I1=(Q/V1);//Primary current in A\n", +"nFL=((V1*I1*0.8)/((V1*I1*0.8)+(Cu+C)))*100;//Full load efficiency of the transformer\n", +"\n", +"//CALCULATIONS\n", +"mprintf('Full load efficiency of the transformer is %3.2f percent',nFL)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.36: Secondary_current_and_maximum_efficiency.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.36, Page 2.67\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=500;//Primary voltage in V\n", +"V2=100;//Secondary voltage in V\n", +"K=(V2/V1);//Turn ratio\n", +"R1=0.04;//Primary resistance in ohm\n", +"R2=0.03;//Secondary resistance in ohm\n", +"Pi=200;//Iron or constant lossses in W\n", +"\n", +"//CALCULATIONS\n", +"I2=sqrt(Pi/(R2+(R1*K^2)));//Secondary current in A\n", +"nmax1=((V2*I2)/((V2*I2)+Pi+Pi))*100;//Maximum efficiency at unit power factor\n", +"nmax8=((V2*I2*0.8)/((V2*I2*0.8)+Pi+Pi))*100;//Maximum efficiency at 0.8 power factor\n", +"\n", +"//OUTPUT\n", +"mprintf('Maximum efficiency at unit power factor is %3.2f percent \nMaximum efficiency at 0.8 power factor is %3.2f percent',nmax1,nmax8)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.37: Constant_losses_and_full_load_copper_losses.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.37, Page 2.68\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"nFL=0.98;//Efficiency of transformer at full load 0.8 power factor\n", +"upf=0.99;//Efficiency of the transformer at half load\n", +"Q=500;//Transformer rating in KVA\n", +"cosq=0.8;//Power facotor\n", +"\n", +"//CALCULATIONS\n", +"L=((Q*1000*cosq)/nFL)-(Q*1000*cosq);//Full load losses in W\n", +"L2=((0.5*Q*1000*100)/99)-(0.5*Q*1000);//Half load losses in W\n", +"A=[0.25,0.25;\n", +" 1,0.25]\n", +"B=[(0.25*L);\n", +" L2]\n", +"A=inv(A)*B;//Soving for Pi and Pc\n", +"\n", +"\n", +"//OUTPUT\n", +"mprintf('Constant losses are %3.2f W\nFull load copper losses are %3.2f W',A(1),A(2))\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.38: All_day_efficiency.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.38, Page 2.71\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"T=4;//Total loss in kW\n", +"Q=120;//Rating of transformer in KVA\n", +"DF=4;//Duration of operation at full load in h\n", +"DH=4;//Duration of operation at half load in h\n", +"DN=16;//Duration of operation at no load in h\n", +"\n", +"//CALCULATIONS\n", +"EF=(Q*1*T);//Energy delivered for 4 hours full load in kWh\n", +"EH=(0.5*Q*1*T);//Energy delivered for 4 hours half load in kWh\n", +"EN=0;//Energy delivered for 16 hours\n", +"E24=(EH+EF+EN);//Total energy deliverd for 24 hours in kWh\n", +"C=(1.5*24);//Constant losses for 24 hours in kWh\n", +"C4=(1.5*4);//full load copper losses for 4 hours in kWh\n", +"Ch4=(0.5^2*1.5*4);//Half load copper losses for 4 hours in kWh\n", +"CN=0;//No load copper loss for 16 hours\n", +"TE=(C+C4+Ch4+CN);//Total energy losses for 24 hours\n", +"n=(E24/(E24+TE))*100;//All day efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('All day efficiency is %3.1f percent',n)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.39: All_day_efficiency.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.39, Page 2.72\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Q=10;//Rating of transformer in KVA\n", +"n=0.96;//Full load efficiency\n", +"DN=12;//Duration of no load in h\n", +"DH=6;//Duration of half load in h\n", +"D4=4;//Duration of 1/4th load in h\n", +"DF=2;//Duration of full load in h\n", +"\n", +"//CALCULATIONS\n", +"O=(Q*1);//Full load output in kW\n", +"L=((O/n)-O)*1000;//Full load total losses in W\n", +"Fcu=(L/2);//Full load copper ;osses in W\n", +"Fc=Fcu;//Constant losses\n", +"LN=0;//No load energy delivered for 12 h\n", +"LF=(DF*O);//Full load energy delivered for 2 hours\n", +"L6=(DH*O*0.5);//Half load energy delivered for 6 hours\n", +"L4=(D4*O*0.25);//1/4th load energy delivered for 4 hours\n", +"TE=(LN+LF+L6+L4);//Total energy delivered for 24 hours in kWh\n", +"LLC=(Fc*24);//Constant losses for 24 h\n", +"LLF=(DF*Fc);//Full load copper losses delivered for 2 hours\n", +"LL6=(DH*Fc*0.5^2);//Half load copper losses delivered for 6 hours\n", +"LL4=(D4*Fc*0.25^2);//1/4th load copper losses delivered for 4 hours\n", +"LTE=(LLC+LLF+LL6+LL4)/1000;//Total copper losses delivered for 24 hours in kWh\n", +"nall=((TE/(TE+LTE))*100);//All day efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('All day efficiency is %3.1f percent',nall)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.3: Number_of_turns_and_induced_emf.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.3, Page 2.6\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"A=(40*10^-4);//Area of cross section of the core A in m^2\n", +"B=8;//Maximum flux density in the core B in Wb/m^2\n", +"V1=2000;//Primary voltage in V\n", +"V2=200;//Secondary voltage in V\n", +"f=50;//Frequency in Hz\n", +"\n", +"//CALCULATIONS\n", +"N1=(V1/(4.44*B*A*f));//Number of turns in the primary\n", +"N2=(V2/(4.44*f*A*B));//Number of turns in the secondary\n", +"\n", +"//OUTPUT\n", +"mprintf('Number of turns in the primary is %3.0f \nNumber of turns in the secondary is%3.0f',N1,N2)\n", +" \n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.40: Current_and_output_of_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.40, Page 2.75\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"VLP=11000;//Primary line voltage in V\n", +"VLS=440;//Secondary line voltage in V\n", +"Vphp=11000;//Primary phase voltage in V\n", +"Vphs=(440/sqrt(3));//Secondary phase voltage in V\n", +"ILP=4;//Primary line current in A\n", +"q=0.8;//Power factor\n", +"\n", +"//CALCULATIONS\n", +"Iphp=(ILP/sqrt(3));//Primary phase current in A\n", +"K=(Vphs/VLP);//Turn ratio\n", +"I2ph=(Iphp/K);//Secondary phase current in A\n", +"P=(sqrt(3)*VLS*VLP*q)/10^5;//Output of the transformer in kW\n", +"\n", +"//OUTPUT\n", +"mprintf('Primary phase current is %3.2f A and Secondary phase current is %3.0f A \nOutput of the transformer is%3.0f kW',Iphp,I2ph,P)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.41: Parameters_of_an_ideal_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.41, Page 2.77\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"VLP=2200;//Primary line voltage in V\n", +"Vphp=VLP;//Primary phase voltage in V\n", +"VLS=440;//Secondary line voltage in V\n", +"Vload=440;//Load line phase voltage in V\n", +"Z=complex(8,6);//Load impedence in complex form\n", +"\n", +"//CALCULATIONS\n", +"X=sqrt((real(Z))^2+(imag(Z))^2);//X value for load current\n", +"Y=atand(imag(Z)/real(Z));//Phase angle in degree\n", +"ILS=(VLS/X);//Load current in A\n", +"PS=(sqrt(3)*VLS*ILS*cosd(Y))/1000;//Power delivered by secondary in kW\n", +"K=((Vload/sqrt(3))/VLP);//Turn ratio\n", +"IPS=(sqrt(3)*ILS);//Secondary phase current in A\n", +"IPP=(K*IPS);//Primary phase current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Load delivered by the secondary is %3.1f kW \nb)Current in primary is %3.1f A and Current in secondary is %3.2f A',PS,IPP,IPS)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.4: Number_of_turns.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.4, Page 2.7\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=2500;//primary voltage in V\n", +"V2=200;//Secondary voltage in V\n", +"e=(30*0.9);//Effective side of magnetic core in cm\n", +"A=(30*30*0.9*0.9*10^-4);//Area of cross section of the limb in m^2\n", +"B=1;//Maximum flux density in Wb/m^2\n", +"q=(B*A);//Maximum flux in Wb\n", +"f=50;//Frequency of power supply in Hz\n", +"\n", +"//CALCULATIONS\n", +"N1=(V1/(4.44*f*q));//Number of turns in the primary\n", +"N2=(V2/(4.44*f*q));//Number of turns in the secondary\n", +"\n", +"//OUTPUT\n", +"mprintf('Number of turns in the primary are %3.0f turns and Number of turns in the secondary are %3.0f turns',N1,N2)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.5: Magnetising_and_iron_loss_components.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.5, Page 2.10\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Io=0.8;//No load primary current in A\n", +"Wo=75;//No load primary poewr in W\n", +"V1=220;//Primary voltage in V\n", +"f=50;//Supply frequency in Hz\n", +"\n", +"//CALCULATIONS\n", +"Iw=(Wo/V1);//Iron loss component in A\n", +"Im=sqrt(Io^2-Iw^2);//Magnetising component in A\n", +"\n", +"//OUTPUT\n", +"mprintf('Iron loss component is %3.2f A \nMagnetising component is %3.3f A',Iw,Im)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.6: Core_and_iron_loss_and_magnetising_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.6, Page 2.11\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Io=6;//No load primary current in A\n", +"cosq=0.3;//Power factor\n", +"V1=220;//Primary voltage in V\n", +"V2=2200;//Secondary voltage in V\n", +"\n", +"//CALCULATIONS\n", +"Wo=(V1*Io*cosq);//Core loss in W\n", +"Iw=(Io*cosq);//Iron loss current in A\n", +"Im=sqrt(Io^2-Iw^2);//Magnetsising current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)Core loss is %3.0f W\n(b)Iron loss current is %3.1f A\n(c)Magnetising current is %3.2f A',Wo,Iw,Im)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.7: Properties_of_a_transformer.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.7, Page 2.12\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"V1=200;//Primary voltage in V\n", +"V2=2000;//Secondary voltage in V\n", +"Io=7;//Primay no load current in A\n", +"Wo=180;//Primary no load power in W\n", +"R1=0.05;//Primary winding resistance in ohm\n", +"\n", +"//CALCULATIONS\n", +"Fe=(Wo-(Io^2*R1));//Iron loss or core loss alone in W\n", +"cosq=(Wo/(V1*Io));//No load power factor\n", +"Iw=(Wo/V1);//Working component of current in A\n", +"Im=sqrt(Io^2-Iw^2);//Magnetising current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('(a)The core loss is %3.2f W\n(b)No load power factor is %3.3f lagging\n(c)Working component of current is %3.1f A\n(d)Magnetising current is %3.2f A',Fe,cosq,Iw,Im)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.8: Primary_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.8, Page 2.14\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Io=6;//No load prmary current in A\n", +"cosqo=0.2;//Primary no load power factor \n", +"I2=125;//Secondary load current in A\n", +"cosq2=0.8;//Secondary load power factor \n", +"V1=400;//Primary terminal voltage in V\n", +"V2=100;//Secondary terminal voltage in V\n", +"\n", +"//CALCULATIONS\n", +"K=(V2/V1);//Voltage transformation ratio\n", +"I2i=(K*I2);//Secondary current in A\n", +"q=(acos(cosqo)-acos(cosq2));//Value of angle in degees\n", +"I1=sqrt((Io^2)+(I2i^2)+(2*Io*I2i*cos(q)));//Primary current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('Primary current is %3.2f A',I1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2.9: No_load_current_and_phase_angle.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-2, Example 2.9, Page 2.16\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N1=760;//Number of turns in the primary\n", +"N2=180;//Number of turns in the secondary\n", +"I2=70;//Secondary load current in A\n", +"cosq=0.8;//Secondary load power factor\n", +"I1=30;//Primary current in A\n", +"cosq1=0.71;//Primary current power factor\n", +"\n", +"//CALCULATIONS\n", +"K=(N2/N1);//Ratio of turns\n", +"I2i=(K*I2);//Secondary current in A\n", +"I1i=complex((I1*cosq1),(I1*sind(acosd(cosq1))));//Primary current in A\n", +"I2c=complex((I2i*cosq),(I2i*sind(acosd(cosq))));//Secondary current in A\n", +"A1=sqrt((real(I1i))^2+(imag(I1i))^2);\n", +"A2=(atand(imag(I1i)/real(I1i)));\n", +"B1=sqrt((real(I2c))^2+(imag(I2c))^2);\n", +"B2=(atand(imag(I2c)/real(I2c)));\n", +"C=(A1*cosd(A2))-(B1*cosd(B2));\n", +"D=(A1*sind(A2))-(B1*sind(B2));\n", +"q=atand(D/C);//Phase angle in degree\n", +"p=cosd(q);//Power factor\n", +"Io=(D/sind(q));//No load current in A\n", +"\n", +"//OUTPUT\n", +"mprintf('No load current of the transformer is %3.2f A and its phase angle is %3.2f degree',Io,q)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Electrical_Engineering_by_R_Anandanatarajan/3-Three_Phase_Induction_Motor.ipynb b/Electrical_Engineering_by_R_Anandanatarajan/3-Three_Phase_Induction_Motor.ipynb new file mode 100644 index 0000000..be4a6d2 --- /dev/null +++ b/Electrical_Engineering_by_R_Anandanatarajan/3-Three_Phase_Induction_Motor.ipynb @@ -0,0 +1,312 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 3: Three Phase Induction Motor" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.1: Frequency_of_rotor_current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-3, Example 3.1, Page 3.6\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N=900;//Rotor speed in rpm\n", +"f=50;//Power supply frequency in Hz\n", +"P=6;//No. of poles\n", +"\n", +"//CALCULATIONS\n", +"Ns=(120*f)/P;//Synchronous speed in rpm\n", +"s=((Ns-N)/Ns)*100;//%slip \n", +"f1=(s*f)/100;//Frequency of rotor current in Hz\n", +"\n", +"//OUTPUT\n", +"mprintf('Slip of a 3 phase motor is %i percent\nFrequency of rotor current is %i Hz',s,f1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.2: Full_load_speed_of_the_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-3, Example 3.2, Page 3.6\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"N=600;//Speed of 12 pole 3 phase alternator in rpm\n", +"P=12;//No. of poles of alternator\n", +"n=6;//No. of poles in induction motor\n", +"s=2.5;//slip of the motor in %\n", +"\n", +"//CALCULATIONS\n", +"f=(N*P)/120;//Alternator supply frequency in Hz\n", +"Ns=(120*f)/n;//Synchronous speed in rpm\n", +"N1=(Ns-((s*Ns)/100));//Full load speed of the motor when the slip is 2.5%\n", +"\n", +"//OUTPUT\n", +"mprintf('Full load speed of the motor when the slip is 2.5 percent = %irpm',N1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.3: Slip_and_speed_of_rotor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-3, Example 3.3, Page 3.7\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=6;//Number of poles\n", +"f=50;//Supply frequency in Hz\n", +"f1=3;//Rotor current frequency in Hz\n", +"\n", +"//CALCULATIONS\n", +"s=(f1/f)*100;//Slip of the motor in %\n", +"Ns=(120*f)/P;//Synchronous speed in rpm\n", +"N=(Ns-((s*Ns)/100));//Speed of the motor in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Slip of the motor is %i percent\nSpeed of the motor is %i rpm',s,N)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.4: Shaft_output_and_torque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-3, Example 3.4, Page 3.12\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"VL=440;//Supply line voltage in V\n", +"P=4;//Number of poles\n", +"IL=75;//Line current in A\n", +"cosx=0.8;//Power factor\n", +"n=0.8;//Efficiency of the motor\n", +"s=0.03;//slip of the motor\n", +"f=50;//Frequency in Hz\n", +"\n", +"//CALCULATIONS\n", +"Pm=(sqrt(3)*VL*IL*cosx*n);//Output power in W\n", +"Ns=(120*f)/P;//Synchronous speed in rpm\n", +"N=(1-s)*Ns;//Actual speed in rpm\n", +"\n", +"//OUTPUT\n", +"mprintf('Shaft output power is %3.0f W\nActual speed is %i rpm',Pm,N)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.5: Parameters_of_induction_motor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-3, Example 3.5, Page 3.13\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=6;//Number of poles\n", +"f=50;//Supply frequency in Hz\n", +"Tm=120;//Shaft torque in N.m\n", +"f1=2;//Rotor current frequency in Hz\n", +"L=5;//Amount of constant losses in N.m\n", +"C=500;//Amount of core losses in W\n", +"\n", +"//CALCULATIONS\n", +"Ns=(120*f)/P;//Synchronous speed in rpm\n", +"s=(f1/f);//Slip of the motor \n", +"N=(1-s)*Ns;//Actual speed in rpm\n", +"P=(2*3.14*N*Tm)/60;//Shaft power in W\n", +"Pm=(2*3.14*N*(Tm+L))/60000;//Mechanical power output in kW\n", +"R=(s*Pm)/(1-s);//Rotor copper losses in kW\n", +"I=(Pm+R+(L/10));//Motor input in kW\n", +"n=(Pm/I)*100;//Machine efficiency\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Mechanical power output is %3.3f kW\nb)Rotor copper losses is %3.2fkW\nc)Motor input is %3.3f kW\nd)Machine efficiency is %3.1f percent',Pm,R,I,n)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.6: Slip_and_torque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-3, Example 3.6, Page 3.17\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"VL=11000;//Supply line voltage in V\n", +"P=12;//Number of poles\n", +"f=50;//Supply frequency in Hz\n", +"R2=0.2;//Rotor resistance in ohm\n", +"X2=1.2;//Rotor reactance at stand still in ohm\n", +"N=480;//Full load speed in rpm\n", +"\n", +"//CALCULATIONS\n", +"s=(R2/X2);//Slip at maximum torque\n", +"Ns=(120*f)/P;//Synchronous speed in rpm\n", +"s1=(Ns-N)/Ns;//Slip at full load\n", +"T=((R2^2+(s1^2*X2^2))/((2*X2)*(s1*R2)));//Ratio of maximum and full load torque\n", +"T1=((R2^2+X2^2)/(2*X2*R2));//Ratio of maximum and starting torque\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Slip at maximum torque is %3.2f \nb)Ratio of maximum and full load torque is %3.2f \nc)Ratio of maximum and starting torque is %3.2f',s,T,T1)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3.7: Maximum_torque_and_starting_torque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-3, Example 3.7, Page 3.18\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=6;//Number of poles\n", +"f=50;//Supply frequency in Hz\n", +"R2=0.4;//Rotor reisitance in ohm\n", +"X2=4;//Rotor standstill reactance in ohm\n", +"T1=2;//Ratio of maximum torque to starting torque\n", +"\n", +"//CALCULATIONS\n", +"Ns=(120*f)/P;//Synchronous speed in rpm\n", +"Sm=(R2/X2);//Slip at maximum torque\n", +"NTM=(Ns*(1-Sm));//Speed of the motor at maximum torque in rpm\n", +"T=((R2^2+X2^2)/(2*R2*X2));//Ratio of maximum torque to starting torque\n", +"Rext=(sqrt(X2^2/((2*T1)-1))-R2);//Additional resistance required for the ratio of maximum torque to the statring torque to be 2 in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('a)Speed of the motor at maximum torque is %i rpm \n b)Ratio of maximum torque to starting torque is %3.2f \n c)Additional resistance required for the ratio of maximum torque to the starting torque to be 2 is %3.1f ohm',NTM,T,Rext)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Electrical_Engineering_by_R_Anandanatarajan/5-Synchronous_and_Special_Machines.ipynb b/Electrical_Engineering_by_R_Anandanatarajan/5-Synchronous_and_Special_Machines.ipynb new file mode 100644 index 0000000..2c15556 --- /dev/null +++ b/Electrical_Engineering_by_R_Anandanatarajan/5-Synchronous_and_Special_Machines.ipynb @@ -0,0 +1,182 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 5: Synchronous and Special Machines" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.1: Emf_generated_and_line_voltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-5, Example 5.1, Page 5.6\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"f=50;//Frequency in Hz\n", +"Z=200;//Number of conductors\n", +"kp=1;//Pitch factor\n", +"kd=0.96;//Distribution factor\n", +"q=0.05;//Flux in Wb\n", +"\n", +"//CALCULATIONS\n", +"Eph=(2.22*kp*kd*f*q*Z);//EMF generated per phase in V\n", +"LV=(Eph*sqrt(3));//Line voltage in V\n", +"\n", +"//OUTPUT\n", +"mprintf('(i)Emf generated per phase is %3.1f V \n(ii)Line voltage is %3.1f V',Eph,LV)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.2: Induced_emf_per_phase.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-5, Example 5.2, Page 5.7\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=8;//Number of poles\n", +"f=50;//Frequency in Hz\n", +"Z=(36*8);//Number of conductors\n", +"q=0.04;//Flux in Wb\n", +"kp=1;//Pitch factor\n", +"kd=1;//Distribution factor\n", +"\n", +"//CALCULATIONS\n", +"Eph=(2.22*kp*kd*f*q*Z);//EMF generated per phase in V\n", +"\n", +"//OUTPUT\n", +"mprintf('Induced emf per phase is %3.1f V',Eph)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.3: Number_of_conductors.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-5, Example 5.3, Page 5.7\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"P=8;//Number of poles\n", +"EL=11000;//Line voltage of the alternator in kV\n", +"Eph=(EL/sqrt(3));//Phase voltage per pole in V\n", +"kp=1;//Pitch factor\n", +"kd=0.98;//Distribution factor\n", +"q=0.17;//Flux in Wb\n", +"f=50;//Frequency in Hz\n", +"\n", +"//CALCULATIONS\n", +"Z=(Eph/(2.22*kp*kd*f*q));//Number of conductors per phase\n", +"\n", +"//OUTPUT\n", +"mprintf('Number of conductors per phase is %3.0f',Z)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5.4: Synchronous_reactance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"//Chapter-5, Example 5.3, Page 5.7\n", +"//=============================================================================\n", +"clc\n", +"clear\n", +"\n", +"//INPUT DATA\n", +"Eph=(6.6*10^3)/sqrt(3);//Phase voltage in V\n", +"Isc=145;//Short circuit current in A\n", +"Ra=1;//Resistance of stator winding in ohm\n", +"\n", +"//CALCULATIONS\n", +"Zs=(Eph/Isc);//Synchronous impedence in ohm\n", +"Xs=sqrt(Zs^2-Ra^2);//Synchronous reactance in ohm\n", +"\n", +"//OUTPUT\n", +"mprintf('Synchronous reactance is %3.2f ohm',Xs)\n", +"\n", +"//=================================END OF PROGRAM==============================" + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/Electrical_Engineering_by_R_Anandanatarajan/6-Transmission_and_Distribution.ipynb b/Electrical_Engineering_by_R_Anandanatarajan/6-Transmission_and_Distribution.ipynb new file mode 100644 index 0000000..c3b782f --- /dev/null +++ b/Electrical_Engineering_by_R_Anandanatarajan/6-Transmission_and_Distribution.ipynb @@ -0,0 +1,32 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 6: Transmission and Distribution" + ] + }, +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} |