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 /Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen | |
| 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 'Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen')
10 files changed, 3523 insertions, 0 deletions
diff --git a/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/1-magnetic_circuit.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/1-magnetic_circuit.ipynb new file mode 100644 index 0000000..5189dbf --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/1-magnetic_circuit.ipynb @@ -0,0 +1,235 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 1: magnetic circuit" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: current.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Bc=0.8;Hc=510;Bg=0.8;\n", +"A=4*%pi*10^-7;lg=1.5*10^-3;\n", +"lc=0.36;N=500;\n", +"\n", +"Fg=Bg/A*(2*lg)\n", +" \n", +"Fc=Hc*lc\n", +" \n", +"F=Fc+Fg\n", +"\n", +"i=F/N\n", +" \n", +"Pre=Bc/Hc\n", +" \n", +"RelPre=Pre/A\n", +" \n", +"F=Hc*lc\n", +"\n", +"i=F /N\n", +" \n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: fluxdensity.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"A=4*%pi*10^-7;lc=360;\n", +"N=500;i=4;lg=2*10^-3;\n", +"\n", +"m=-A*(lc/lg)\n", +" \n", +"c=(N*i*A)/(lg)\n", +" \n", +"Hc=(N*i)/(lc)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: airgapflux.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"N1=500;I1=10;N2=500;I2=10;\n", +"Ibafe=3*52*10^-2;A=4*%pi*10^-7;\n", +"b=1200;Ag=4*10^-4;Ac=4*10^-4;\n", +"lg=5*10^-3;Ibecore=51.5*10^-2;c=2.067*10^-4;\n", +"d=4.134*10^-4;\n", +"\n", +"F1=N1*I1\n", +" \n", +"F2=N2*I2\n", +" \n", +"Pre=1200*A\n", +" \n", +"Rbafe=(Ibafe)/(Pre*Ac)\n", +" \n", +"Rg=lg/(A*Ag)\n", +" \n", +"Rbecore=Ibecore/(Pre*Ac)\n", +" \n", +"Bg=d/(Ag)\n", +" \n", +"Hg=Bg/A\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: magneticflux.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Irad=20;Orad=25;Dia=22.5\n", +"N=250;i=2.5;\n", +"l=2*%pi*Dia*10^-2;\n", +"B=1.225;\n", +"\n", +"radius=1/2*(Irad+Orad)\n", +" \n", +"H=(N*i)/l\n", +" \n", +"A=%pi*((Orad-Irad)/2)^2*10^-4\n", +" \n", +"z=(1.225)*(%pi*6.25*10^-4)\n", +" \n", +"y=(N*z)\n", +" \n", +"L=(y/i)\n", +" \n", +"core=(B/H)\n", +"\n", +"l=(2*%pi*22.5*10^-2)\n", +" \n", +"Rcore=(l)/(core*A)\n", +"\n", +"L=(N^2)/(Rcore)\n", +" \n", +" \n", +" \n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5: fluxvoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"n=500;E=100;A=0.001;b=1/120;\n", +"f=1.2;\n", +" \n", +"max1=(E/1000)*(b)\n", +" \n", +"max2=(f*A)\n", +" \n", +"E=(120*n*max2*2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6: dimensions.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"lg=0.4*10^-2;Bg=0.8;\n", +"Hm=42*10^3;A=4*%pi*10^-7;Ag=2.5*10^-4;\n", +"Bm=0.95;\n", +"\n", +"Hg=Bg/A\n", +" \n", +"lm=(lg/Hm)*Hg\n", +" \n", +"Am=(Bg*Ag)/(Bm)" + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/10-POWER_SEMICONDUCTOR_CONVERTERS.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/10-POWER_SEMICONDUCTOR_CONVERTERS.ipynb new file mode 100644 index 0000000..8338459 --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/10-POWER_SEMICONDUCTOR_CONVERTERS.ipynb @@ -0,0 +1,356 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 10: POWER SEMICONDUCTOR CONVERTERS" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: firingangle.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Ka=0.09;N=1000;\n", +"Ia=30;Ra=0.4;V=120;\n", +"RevEa=-90;\n", +"\n", +"Ea=Ka*N\n", +" \n", +"Vo=Ea+(Ia*Ra)\n", +"\n", +"a=Vo*%pi\n", +" \n", +"b=2*sqrt(2)*V\n", +" \n", +"c=a/b\n", +" \n", +"angle=acosd(c)\n", +"\n", +"P=Vo*Ia\n", +" \n", +"S=V*Ia\n", +" \n", +"Pf=P/S\n", +"\n", +"Vo1=RevEa+(Ia*Ra)\n", +"\n", +"a=Vo1*%pi\n", +" \n", +"b=2*sqrt(2)*V\n", +" \n", +"c=a/b\n", +" \n", +"Angle=acosd(c)\n", +" \n", +"Pdc=Ea*Ia\n", +" \n", +"Pr=Ia^2*Ra\n", +" \n", +"Ps=Pdc-Pr" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: ouputvoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Vp=120;Angle=60;\n", +"\n", +"t0=%pi/2\n", +"\n", +"t1=t0:0.01:(210/360*2*%pi);\n", +"\n", +"integrate('2^.5*120*sin(t)','t',t0,t1)\n", +"\n", +"Vo=((3*sqrt(6))/(2*%pi))*120*cos(%pi*Angle/180)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: powerfactor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=480;Ka=0.3;N=1500;\n", +"Ia=130;Ra=0.1;No=1000;\n", +"\n", +"Vp=V/sqrt(3)\n", +" \n", +"Ea=Ka*N\n", +"\n", +"Vo=Ea+(Ia*Ra)\n", +"\n", +"a=Vo*%pi\n", +" \n", +"b=3*sqrt(6)*Vp\n", +" \n", +"c=a/b\n", +" \n", +"Angle=acosd(c)\n", +" \n", +"IA=sqrt(2/3)*Ia\n", +" \n", +"S=3*Vp*IA\n", +"\n", +"Ps=Vo*Ia\n", +" \n", +"Pf=Ps/S\n", +" \n", +"Ea1=Ka*No\n", +"\n", +"Vo1=-300+(Ia*Ra)\n", +"\n", +"a=Vo1*%pi\n", +"\n", +"b=3*sqrt(6)*Vp\n", +" \n", +"c=a/b\n", +" \n", +"Angle=acosd(c)\n", +"\n", +"Pdc=Ea1*Ia\n", +" \n", +"Pr=Ia^2*Ra\n", +" \n", +"Ps=Pdc-Pr" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: controlfiringangle.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"OutP=100*0.746;Eff=0.8;\n", +"Pf=0.85;V=460;\n", +"\n", +"S=OutP/(Eff*Pf)\n", +" \n", +"Il=S/(sqrt(3)*V)\n", +" \n", +"Ip=Il/sqrt(3)\n", +"\n", +"Is=Ip/sqrt(2)\n", +" \n", +"Vs=sqrt(2)*V\n", +" \n", +"Angle=acosd(Pf)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5: supplypower.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Ka=0.1;Ra=0.2;N=400;\n", +"Ia=100;V=120;N1=350;\n", +"Io=-100;\n", +"\n", +"Ea=Ka*N\n", +" \n", +"Vo=Ea+(Ia*Ra)\n", +"\n", +"Pmotor=Ea*Ia\n", +" \n", +"Pr=Ia^2*Ra\n", +" \n", +"Ps=V*Ia*0.5\n", +"\n", +"Vo=Ea1+(Ia*Ra)\n", +"\n", +"Vo=Ea1+(Io*Ra)\n", +" \n", +"Pmotor1=Ea1*Io\n", +" \n", +"Pr1=Ia^2*Ra\n", +" \n", +"Ps=V*Io*1/8" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6: outputpower.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=300;I=540;Angle=45;\n", +"\n", +"t0=0:0.1:%pi;\n", +" \n", +"t=0;\n", +"integrate('540*sin((x-45*%pi/180))','x',t,t0)/%pi\n", +"\n", +" Is=242.89;\n", +" \n", +"Ps=V*I \n", +" \n", +"Vo1=(4*V)/(%pi*sqrt(2))\n", +"\n", +"Pout=Vo1*Io/sqrt(2)*cos(%pi*Angle/180)\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7: phasevoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=600;Phase=3;\n", +"\n", +"Vl=sqrt(2/3)*V\n", +" \n", +"Vp=V*sqrt(2)/Phase" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8: angleshift.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=120;Vo=50;Vo1=100;\n", +" \n", +"Angleshift=(Vo^2/V^2)*180\n", +" \n", +"Angleshift=(Vo1^2/V^2)*180\n", +" \n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9: outputcurrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"I=100;\n", +"PulseScr=sqrt(1/3);\n", +"PulseOut=sqrt(2/3);\n", +"\n", +"a=PulseScr\n", +"\n", +"msScr=a*I\n", +" \n", +"b=PulseOut\n", +" \n", +"RmsOut=b*I" + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/2-TRANSFORMERS.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/2-TRANSFORMERS.ipynb new file mode 100644 index 0000000..67377a1 --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/2-TRANSFORMERS.ipynb @@ -0,0 +1,454 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 2: TRANSFORMERS" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: powertake.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"resimp=9;inte=1;V=10;A=1/3;\n", +"\n", +"I=V/(inte+resimp)\n", +" \n", +"P=I^2*resimp\n", +"\n", +"R=(A^2)*(resimp)\n", +" \n", +"I1=V/(inte+I)\n", +" \n", +"P1=I1^2*inte" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: parameters.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Vh=220;Ih=4.55;Wl=100;\n", +"Vl=150;Il=2.5;Wh=215;\n", +"Vhrated=2200;Vlrated=220;\n", +"Ihrated=4.55;Ilrated=45.5;\n", +"\n", +"Rcl=Vlrated^2/Wl\n", +" \n", +"Poc=(Vlrated^2/Rcl)\n", +" \n", +"Icl=Vlrated/Rcl\n", +"\n", +"Iml=(Il^2-Icl^2)^(1/2)\n", +" \n", +"Xml=Vlrated/Iml\n", +"\n", +"A=Vhrated/Vlrated\n", +" \n", +"Rch=A^2*Rcl\n", +"\n", +"Xmh=A^2*Xml\n", +"\n", +"Reqh=215/Ihrated^2\n", +" \n", +"Psc=Ihrated^2*Reqh\n", +"\n", +"Zeqh=Vl/Ihrated\n", +" \n", +"Xeqh=(Zeqh^2-Reqh^2)^(1/2)\n", +"\n", +"Reql=Reqh/A^2\n", +"\n", +"Xeql=Xeqh/A^2\n", +"\n", +"P=(Poc/(Vlrated*Il))\n", +"\n", +"Psh=(Psc/(Vl*Ihrated))\n", +"\n", +" \n", +" \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: voltageregulatio.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Fullload=75;Ia=4.55;Vl=2200;\n", +"\n", +"Fulload=Fullload/100\n", +" \n", +"Ih=Fullload*Ia\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[x1,y1]=polar2rect(Vl,0)\n", +" \n", +"[x2,y2]=polar2rect(35.46,-53.13)\n", +" \n", +"[x3,y3]=polar2rect(106.73,36.87)\n", +" \n", +"X1=x1+%i*y1\n", +"\n", +"X2=x2+%i*y2\n", +" \n", +"X3=x3+%i*y3\n", +" \n", +"X=X1+X2+X3\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[V,Angle]=rect2polar(2306.5,35.67)\n", +" \n", +"VolReg=(V-Vl)/Vl*100\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"[x1,y1]=polar2rect(Vl,0)\n", +" \n", +"[x2,y2]=polar2rect(35.46,53.13)\n", +"\n", +"[x3,y3]=polar2rect(106.73,143.13)\n", +" \n", +"X1=x1+%i*y1\n", +" \n", +"X2=x2+%i*y2\n", +" \n", +"X3=x3+%i*y3\n", +"\n", +"X=X1+X2+X3\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[V1,Angle1]=rect2polar(2135.89,92.4)\n", +" \n", +"VolReg=(V1-Vl)/Vl*100" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: efficiency.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=0.75;I=10000;A=0.6;\n", +"Pc=100;Reqh=10.4;Ih=(0.75*4.55)^2;\n", +"Reql=0.104;V2=220;B=1;\n", +"\n", +"Pout=V*I*A\n", +"\n", +"Pcu=(Ih*Reqh)\n", +"\n", +"Eff=Pout/(Pout+Pc+Pcu)\n", +" \n", +"Eff=Pout/(Pout+Pc+Pcu)*100\n", +"\n", +"I2=(100/0.104)^(1/2)\n", +" \n", +"Pout1=V2*I2*B\n", +"\n", +"Eff1=Pout1/(Pout1+Pc+Pcu)\n", +"\n", +"Eff1=Pout1/(Pout1+Pc+Pcu)*100" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5: efficiency.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Power=50;Lo1=0.5;\n", +"Lo2=0.75;Lo3=1;Lo4=1.1;\n", +"Pf1=1;Pf2=0.8;Pf3=0.9;\n", +"Pf4=1;Ho1=6;Ho2=6;Ho3=3;Ho4=3;Ho=6;Pc=200;\n", +"Pcu=500;\n", +"\n", +"EngOut=(Lo1*Power*Ho1*Pf1)+(Lo2*Power*Ho2*Pf2)+(Lo3*Power*Ho3*Pf3)+(Lo4*Power*Ho4*Pf4)\n", +" \n", +"A=Pc/1000\n", +"\n", +"TotalHour=Ho+Ho1+Ho2+Ho3+Ho4\n", +"\n", +"Coreloss=A*TotalHour\n", +" \n", +"B=Pcu/1000\n", +" \n", +"Copperloss=(Lo1^2*B*Ho1)+(Lo2^2*B*Ho2)+(Lo3^2*B*Ho3)+(Lo4^2*B*Ho4)\n", +" \n", +"Totalloss=Coreloss+Copperloss\n", +" \n", +"Eff=EngOut/(EngOut+Totalloss)*100" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6: kvarating.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"P=100000;Vs=2000;Vp=200;Ih=500;\n", +"Vl=2000;\n", +"\n", +"Iab=P/Vp\n", +"\n", +"Ibc=P/Vs\n", +" \n", +"Il=Ih+50\n", +" \n", +"Vh=Vl+200\n", +" \n", +"Kva1=(Vl*Il)/(1000)\n", +" \n", +"Kva2=(Vh*Ih)/(1000)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7: volcurregulatio.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Power=120000;;Phase=3;\n", +"V=230;Pri=2300;Sec=230;Z=0.012+%i*0.016;Pf=0.85;\n", +"\n", +"Is=Power/(sqrt(Phase)*V)\n", +" \n", +"I2=Is/sqrt(Phase)\n", +" \n", +"a=Pri/V\n", +" \n", +"I1=I2/a\n", +" \n", +" Zeq=(Z)*10^2\n", +" \n", +" a=acos(Pf)\n", +" \n", +" Deg=(a*180)/%pi\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[a,b]=polar2rect(Pri,0)\n", +" \n", +"A=a+%i*b\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[c,d]=polar2rect(I1,-Deg)\n", +" \n", +"A1=c+%i*d\n", +" \n", +"A2=A1*(Zeq)\n", +"\n", +"A3=A2+A\n", +" \n", +"V1=2332.4;\n", +"\n", +"PriVol=sqrt(Phase)*V1\n", +"\n", +"VR=(V1-Pri)/Pri*100" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8: supplyvoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Pri=1330;Sec=230;Zl=0.12+%i*0.25;Phase=3;V=230;\n", +"Z=0.8+%i*5;Power=27;\n", +"Zz=0.003+%i*0.015;Pf=0.9\n", +"\n", +"A=(Pri/Sec)^2*(Zl)\n", +"\n", +"Req=4.01;\n", +" \n", +"Xeqh=8.36;\n", +" \n", +"a=(sqrt(Phase)*Pri)/V\n", +" \n", +"Reql=0.8;\n", +" \n", +"Xeql=5;\n", +"Rr=0.003;\n", +" \n", +"Xx=0.015;\n", +" \n", +"R=(Reql+Req)*(1/10^2)+Rr\n", +"\n", +"X=(Xeql+Xeqh)*(1/10^2)+Xx\n", +" \n", +"Vl=V/sqrt(Phase)\n", +" \n", +"Il=(Power*10^3)/(Phase*133)\n", +" \n", +"Angle=-acos(%pi*Pf/180)\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9: impedances.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Vh=2200;Vl=220;Pb=10000;\n", +"I=0.25;a=10;Z=10.4+%i*31.3;\n", +"\n", +"Ib=Pb/Vh\n", +" \n", +"Il=Pb/Vl\n", +" \n", +"Zb=Vh/Ib\n", +" \n", +"Zl=Vl/Il\n", +" \n", +"Ih=I/Ib\n", +" \n", +"Zeq=Z/Zb\n", +" \n", +"Zeql=Z*(1/100)\n", +" \n", +"Zpu=Zeql/Zl\n", +" \n", +"Pcu=Ib^2*10.4\n", +" \n", +"Ppu=Pcu/Pb" + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/3-ELECTROMECHNICAL_ENERGY_CONVERSION.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/3-ELECTROMECHNICAL_ENERGY_CONVERSION.ipynb new file mode 100644 index 0000000..242c322 --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/3-ELECTROMECHNICAL_ENERGY_CONVERSION.ipynb @@ -0,0 +1,159 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 3: ELECTROMECHNICAL ENERGY CONVERSION" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: fieldenergy.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Hc=670;G=5;D=10;A=5;B=10;Bg=1;Z=4*%pi*10^-7;N=250;Area=700;\n", +"\n", +"Lc=2*(A+B)+2*(G+D)\n", +" \n", +"Hg=Bg/Z\n", +" \n", +"Lc=60/100\n", +" \n", +"Hg=Bg/Z\n", +" \n", +"Ni=(Hc*Lc)+(Hg*2*G*10^-3)\n", +" \n", +"I=Ni/N\n", +" \n", +"Vdc=I*G\n", +" \n", +"Wfc=Area/2\n", +" \n", +"Vc=2*(G*10^-2*D*10^-2*0.20)+2*(A*10^-2*B*10^-2*0.10)\n", +" \n", +"Wfc=Wfc*Vc\n", +" \n", +"Wfg=1.0/(2*Z)\n", +" \n", +"Vg=2*(G*10^-2*10*10^-2*0.005)\n", +" \n", +"Wfg=(Wfg*G*10^-2*10^-3)\n", +" \n", +"Wf=Wfc+Wfg" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: mechanicalforce.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"I=3;G=0.05;\n", +"\n", +"Lam=(0.09*I^(.5)/G)\n", +"\n", +"t1=0:0.1:3;\n", +" \n", +"t0=0;\n", +" \n", +"a=integrate('((0.09*2)/(G*I))*I^(0.5)','i',t0,t1)\n", +" \n", +"Wf=((0.09*2)/(G*I))*I^(1.5)\n", +" \n", +"Fm=-0.09*(2/3)*I^(1.5)*(1/G^2)\n", +" \n", +"Wf1=(G^2*Lam^3)/(0.09^2*I)\n", +"\n", +"Lam1=(0.09*I^(.5)/G)\n", +" \n", +"Fm=-((Lam1^3)*2*G)/(I*0.09^2)\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: liftingforce.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"A=4*%pi*10^-7;N=300;V=120;\n", +"R=6;G=5*10^-3;Ag=6*6*10^-4;\n", +"Lg=2*5*10^-3;\n", +"Vo=2*6*6*5*10^-7;\n", +"\n", +"I=V/R\n", +" \n", +"Bg=(A*N*I)/(2*G)\n", +"\n", +"Wf=(Bg^2)/(2*A)*(Vo)\n", +"\n", +"Fm=(Bg^2)/(2*A)*(2*Ag)\n", +"\n", +"L=(N^2*A*Ag)/(Lg)\n", +"\n", +"Irms=V/(sqrt(6^2+15.34^2))\n", +"\n", +"Brms=(A*N*Irms)/(2*G)\n", +"\n", +"Fm=(Brms^2)/(2*A)*(2*Ag)\n", +"\n", +" \n", +" \n", +" \n", +" " + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/4-DC_MACHINES.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/4-DC_MACHINES.ipynb new file mode 100644 index 0000000..b518dfb --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/4-DC_MACHINES.ipynb @@ -0,0 +1,512 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 4: DC MACHINES" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10: speedtorque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Ia=40;Ra=0.25;La=10;\n", +"ACv=265;Angle=30;\n", +"Ka=0.18;\n", +"\n", +"A=cos(%pi*Angle/180)\n", +" \n", +"Vt=(2*sqrt(2)*ACv*A)/(%pi)\n", +" \n", +"Ea=Vt-(Ia*Ra)\n", +" \n", +"N=Ea/Ka\n", +"\n", +"Sec=(Ka*60)/(2*%pi)\n", +" \n", +"T=Sec*Ia\n", +" \n", +"P=Vt*Ia\n", +"\n", +"P=(Ia^2*Ra)+Ea*Ia" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11: firingangle.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=480;RPM=1800;\n", +"Ia=16.5;Ra=0.0874;\n", +"Ka=0.33;\n", +"Angle=30;N=1800;\n", +"\n", +"Vp=(V/sqrt(3))\n", +" \n", +"Vt=(3*sqrt(6)*Vp)/(%pi)\n", +" \n", +"Ea=Vt-(Ia*Ra)\n", +" \n", +"No=(Ea/Ka)\n", +" \n", +"A=cos(%pi*Angle/180)\n", +"\n", +"Vt1=Vt*A\n", +" \n", +"Ea1=Vt1-(Ia*Ra)\n", +" \n", +"No1=Ea1/Ka\n", +"\n", +"Eaf=Ka*N\n", +" \n", +"Vtf=Eaf+(Ia*Ra)\n", +" \n", +" Angle=Vtf/Vt\n", +" \n", +" a=acos(Angle)\n", +" \n", +" Alfa=a*180/%pi\n", +"\n", +"Eas=Vtf-(Ia*Ra)\n", +" \n", +"Nos=Eas/Ka\n", +" \n", +"SpeedReg=(Nos-N)/N*100" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 12: speedcontrol.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Supv=120;Ia=20;Ra=0.5;\n", +"Ka=0.05;\n", +"\n", +"Vt=Ia*Ra\n", +" \n", +"A=(Vt/Supv)\n", +" \n", +"Ea=Supv-(Ia*Ra)\n", +" \n", +"N=Ea/Ka" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1.b: torque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"P=4;A=2;Z=462;Wn=(1000/60)*(2*%pi);z=0.0276;\n", +"Icoil=100;\n", +"\n", +"Ka=(Z*P)/(2*2*%pi)\n", +" \n", +"Wn\n", +" \n", +"Ea=(Ka*z*Wn)\n", +" \n", +"Icoil\n", +"\n", +"Ia=2*Icoil\n", +" \n", +"T=(Ka*z*Ia)\n", +" \n", +"Pa=Ea*Ia\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: terminalvoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Ia=120;Vt=100;Ea=100;\n", +"Ra=0.1;Ar=0.06;EA=98;\n", +"Ifeff=1.4;\n", +"\n", +"Vt1=Ea-(Ia*Ra)\n", +" \n", +"Ifeff1=1-Ar\n", +"\n", +"V=EA-(Ia*Ra)\n", +"\n", +"Ea=(Vt)+(Ia*Ra)\n", +"\n", +"Ifactual=Ifeff+Ar" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: fieldcurrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Eaa=111;Ia=120;Ra=0.1;\n", +"Rfw=80;Eac=85;If=0.5;\n", +"v=100;\n", +"\n", +"Vt=Eaa-(Ia*Ra)\n", +"\n", +"Rf=v/1\n", +" \n", +"Rfc=Rf-Rfw\n", +" \n", +"Rfcrit=Eac/If\n", +" \n", +"Rfc=Rfcrit-Rfw" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: fullloadvoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Ia=120;Ra=0.1;Vt=80;Vt1=75;AB=17;Ea=6;\n", +"\n", +"V=Ia*Ra\n", +"\n", +"Ia=AB/Ra\n", +" \n", +"Ia1=Ea/Ra" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5: seriesturns.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Ia=120;Vt=100;Rsr=0.01;\n", +"Ra=0.1;Ia=120;Rf=100;\n", +"Ifeff=1.45;If=1.01;\n", +"Nf=1200;Ifar=0.06;\n", +"\n", +"It=Ia-If\n", +" \n", +"If=(Vt+It*Rsr)/Rf\n", +" \n", +"Ea=Vt+(It*Rsr)+(Ia*Ra)\n", +" \n", +"Nsr=(Ifeff-If+Ifar)*(1200)/(It)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6: resistance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Vt=100;Ra=0.1;Ia=6;If=0.99;Rfw=80;\n", +"Ia1=5;Iarated=120;N=1000;\n", +"Afl=0.95;Prot=497.5;\n", +"\n", +"Eanl=Vt-(Ia*Ra)\n", +" \n", +"Rf=Vt/If\n", +" \n", +"Rfc=Rf-Rfw\n", +" \n", +"Prot=Ea*5\n", +" \n", +"Eanl=Vt-(Ia*Ra)\n", +" \n", +"Eafl=Vt-(Iarated*Ra)\n", +" \n", +"Wfl=(Eafl/Eanl)*N\n", +" \n", +"Wm=(Wfl/60)*2*%pi\n", +"\n", +"T=(Eafl*Iarated)/Wm\n", +"\n", +"Pout=(Eafl*Iarated)-(Prot)\n", +"\n", +"Pin=(Vt)*(Iarated+If)\n", +"\n", +"Eff=(Pout/Pin)*100\n", +"\n", +"Wfl1=(Eafl/Eanl)*(1/Afl)*N\n", +" \n", +"Wm1=(Wfl1/60)*(2*%pi)\n", +"\n", +"T=(Eafl*Iarated)/(Wm1)\n", +"\n", +"Eff1=(Pout/Pin)*100\n", +"\n", +"Wm=(1000/60)*(2*%pi)\n", +" \n", +"Ka=Eanl/Wm\n", +"\n", +"Ia=1.5*120\n", +" \n", +"Tstart=(Ka*Ia)\n", +"\n", +"Ifeff=If-Ifar\n", +"\n", +"Ea1=93.5\n", +" \n", +"Ka1=(Ea1/Wm)\n", +" \n", +"Tstart1=(Ka1*Ia)\n", +"\n", +"\n", +"\n", +" \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7: ampereturns.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"If=0.99;Vt=100;Ia=120;Ra=0.1;\n", +"Rpm1=932;Rpm2=1000;\n", +"Ifeff=0.86;Nf=1200;\n", +"rpm1=1000;\n", +"Ifeff1=1.32;Rpm3=800;\n", +"EA=65;\n", +"\n", +"Ea=Vt-(Ia*Ra)\n", +" \n", +"Ea1=Rpm2/Rpm1*Ea\n", +"\n", +"Ifar=If-Ifeff\n", +" \n", +"At=Nf*Ifar\n", +" \n", +"Ea2=Rpm2/Rpm3*Ea\n", +"\n", +"Nsr=(Ifeff1-If+Ifar)*(Nf)/(Ia+If)\n", +" \n", +"Ifeff2=If-(Nsr*(Ia+If))/Nf-Ifar\n", +"\n", +"N=(Ea/EA)*Rpm2" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 8: powerresistance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Vt=220;Ia=25;Ra=0.6;\n", +"Rsr=0.4;Rae=0;N=300;\n", +"N1=200;Wm=(2*%pi/60);\n", +"Hp=746;\n", +"\n", +"Ea=Vt-Ia*(Ra+Rsr+Rae)\n", +"\n", +"P=Ea*Ia\n", +" \n", +"a=P/Hp\n", +" \n", +"T=(Ea*Ia)/(N*Wm)\n", +"\n", +"Ksr=T/Ia^2\n", +"\n", +"T1=(N1/N)^2*T\n", +"\n", +"P=(T1*N1*Wm)\n", +"\n", +"a=T1/Ksr\n", +"\n", +"a=sqrt(Ia)\n", +"\n", +"Rae=(Vt-Ea-(Ia*Ra)+(Ia*Rsr))/Ia\n", +"\n", +"P=Ea*Ia\n", +" \n", +"\n", +" \n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 9: startingcurrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Iarated=100;Vt=100;\n", +"Ra=0.1;\n", +"Ia1=200;\n", +"\n", +"Iastart=Vt/Ra\n", +"\n", +"Rae=(Vt-20)/(200)\n", +" \n", +"Ea2=Vt-Iarated*(Ra+Rae)\n", +" \n", +"Rae2=(Vt-Ea2-20)/(200)\n", +"\n", +"Ea3=Vt-Ia*(Ra+Rae2)\n", +" \n", +"Ea3=Vt-Iarated*(Ra+Rae2)\n", +"\n", +"Rae3=(Vt-Ea3-20)/200\n", +" \n", +"Ea4=Vt-Iarated*(Ra+Rae3)\n", +" \n", +"Rae4=(Vt-Ea4-20)/200\n", +"\n", +"Ia=(Vt-Ea4)/Ra\n", +" \n", +"R1=Rae-Rae2\n", +" \n", +"R2=Rae2-Rae3\n", +" \n", +"R3=Rae3-Rae4" + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/5-INDUCTION_ASYNCHRONOUS_MACHINES.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/5-INDUCTION_ASYNCHRONOUS_MACHINES.ipynb new file mode 100644 index 0000000..2a92c3c --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/5-INDUCTION_ASYNCHRONOUS_MACHINES.ipynb @@ -0,0 +1,501 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 5: INDUCTION ASYNCHRONOUS MACHINES" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 10: parasitictorque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"F=60;P=4;N=1740;R1=0.5;\n", +"R2=0.5;X1=1;X2=1;Xm=35;\n", +"I=1.1;Peak=10;H=5;H1=7;\n", +"Phase=3;\n", +"\n", +"Ns=120*F/P\n", +" \n", +"S1=(Ns-N)/Ns\n", +" \n", +"A=R2/S1\n", +"\n", +"Z1=(Xm*%i)*(A+X1*%i)/(A+X2*%i+Xm*%i)\n", +" \n", +"Rth=12.08;\n", +"\n", +"Pg1=Phase*((I*Peak)/sqrt(2))^2*Rth\n", +"\n", +"Wsyn=(Ns/60)*2*%pi\n", +" \n", +"T1=Pg1/Wsyn\n", +"\n", +"Ns1=-(120*H*F)/P\n", +"\n", +"S2=(Ns1-N)/Ns1\n", +" \n", +"B=R2/S2\n", +"\n", +"Xm1=H*Xm\n", +" \n", +"hX2=H*X2\n", +"\n", +"Z2=Xm1*%i*(B+hX2*%i)/(B+hX2*%i+Xm1*%i)\n", +" \n", +"Rth2=0.39;\n", +"\n", +"Pg2=Phase*((0.22*Peak)/sqrt(2))^2*Rth2\n", +"\n", +"Wsyn2=(Ns1/60)*2*%pi\n", +" \n", +"T2=Pg2/Wsyn2\n", +"\n", +"Ns3=(120*H1*F)/P\n", +" \n", +"S3=(Ns3-N)/Ns3\n", +" \n", +"C=R2/S3\n", +"\n", +"Xm3=H1*Xm\n", +"\n", +"hX3=H1*X2\n", +" \n", +"Z3=Xm3*%i*(C+hX3*%i)/(C+hX3*%i+Xm3*%i)\n", +" \n", +"Rth3=0.54;\n", +"\n", +"Pg3=Phase*((0.16*Peak)/sqrt(2))^2*Rth3\n", +"\n", +"Wsyn3=(Ns3/60)*2*%pi\n", +"\n", +"T3=Pg3/Wsyn3" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 11: voltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"S=2;Pole=50;F=50;\n", +"Slip=0.25;\n", +"\n", +"Pole=Pole*10^-2\n", +" \n", +"Vs=S*Pole*F\n", +" \n", +"Vs=Vs*3600/1000\n", +" \n", +"V=(1-Slip)*(Vs)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: frequency.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"F=60;P=4;S=0.05;\n", +"Ns=1800;V=460;Tr=0.5;\n", +"\n", +"Ns=(120*F)/(P)\n", +" \n", +"N=(1-S)*Ns\n", +"\n", +"F2=S*F\n", +" \n", +"Sliprpm=S*Ns\n", +"\n", +"A=S*Tr*V/sqrt(3)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: powerloss.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Hp=15;Loss=750;\n", +"F=60;P=4;N=1728;\n", +"\n", +"Mecp=P+Loss\n", +"\n", +"Ns=120*F/P\n", +"\n", +"S=(Ns-N)/Ns\n", +" \n", +"Pag=Mecp/(1-S)\n", +" \n", +"P2=S*Pag" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: parameters.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Nof=60;NoV=2200;\n", +"NoI=4.5;NoP=1600;\n", +"BF=15;BV=270;BI=25;\n", +"BP=9000;\n", +"R1=2.8;V=2200;F=60;\n", +"Rbl=4.8;\n", +"\n", +"Pnl=1600\n", +"\n", +"Prot=Pnl-(3*NoI^2*R1)\n", +"\n", +"Vt=V/sqrt(3)\n", +" \n", +"Znl=(Vt/NoI)\n", +"\n", +"Rnl=(Pnl)/(3*NoI^2)\n", +"\n", +"Xnl=(Znl^2-Rnl^2)^(1/2)\n", +" \n", +"Rbl=BP/(3*BI^2)\n", +"\n", +"Zbl=(BV)/(sqrt(3)*BI)\n", +" \n", +"Xbl=(Zbl^2-Rbl^2)^(1/2)\n", +"\n", +"Xbl=Xbl*(F/BF)\n", +" \n", +"X1=Xbl/2\n", +" \n", +"Xm=(Xnl-X1)\n", +"\n", +"R=Rbl-R1\n", +" \n", +"R2=((X1+Xm)/(Xm))^2*2\n", +" \n", +"Vth=(Xm)/(X1+Xm)\n", +" \n", +"Rth=(Vth^2)*R1" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: sliptorque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"R1=0.25;X1=0.5;X2=0.5;\n", +"R2=0.2;Xm=30;V=460;\n", +"N=1740;F=60;P=4;Phase=3;No=1800;\n", +"\n", +"V1=V/sqrt(Phase)\n", +"\n", +"Z1=(R1+%i*X1)+(%i*Xm*(R2+%i*X2)/(R2+%i*30.5))\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Rth,Angle]=rect2polar(0.44,0.99)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[a,b]=polar2rect(Rth,Angle)\n", +" \n", +"X=a+%i*b\n", +"\n", +"Ist=(V1/X)\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ist,Angle]=rect2polar(99.56,-224.0)\n", +" \n", +"Wsy=(1800/F)*2*%pi\n", +"\n", +"Vth=(V1*%i*Xm)/(R1+%i*30.5)\n", +" \n", +"Vth=265.3;\n", +" \n", +"Zth=%i*Xm*(%i*R1+%i*X1)/(R1+%i*30.5)\n", +" \n", +"Zth=(%i*Xm*(R1+%i*X1))/(R1+%i*30.5)\n", +" \n", +"Rth=0.24;\n", +"\n", +"Xth=0.49;\n", +"\n", +"a=Vth^2/((Rth+R2)^2+(Xth+X1)^2)\n", +"\n", +"Tst=((Phase/Wsy)*a*R2)\n", +" \n", +"S=(No-N)/No\n", +" \n", +"FL=R2/S\n", +"\n", +"Z1=(R2+%i*X2)+((%i*Xm)*(FL+%i*X2))/(FL+%i*30.5)\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[c,d]=rect2polar(5.78,2.09)\n", +"\n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[z,y]=polar2rect(c,d)\n", +"\n", +"X=z+%i*y\n", +"\n", +"Ifl=V1/X\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ifl,Ang]=rect2polar(40.63,-14.6)\n", +"\n", +" \n", +"Ratio=Ist/Ifl\n", +"\n", +"Pf=cos(%pi*19.8/180)\n", +" \n", +"z=Vth^2/((Rth+FL)^2+(Xth+X2)^2)\n", +" \n", +"T=(Phase/Wsy)*z*FL\n", +"\n", +"Pag=T*Wsy\n", +" \n", +"P2=S*Pag\n", +" \n", +"Pmech=(1-S)*Pag\n", +" \n", +"Pout=Pmech-1700\n", +" \n", +"Pinp=Phase*V1*Ifl*Pf\n", +"\n", +"EFFMotor=Pout/Pinp*100\n", +" \n", +" EFFint=(1-S)*100\n", +" \n", +" Stmax=R2/(Rth^2+(Xth+X2)^2)^0.5\n", +" \n", +" X=Rth+(Rth^2+(Xth+X2)^2)^0.5\n", +"\n", +" \n", +" Tmax=(Phase/2*Wsy)*X\n", +" \n", +" T=Tmax/T\n", +" \n", +" \n", +" \n", +" \n", +" \n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5: resistance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"F=60;P=6;RPM=1140;\n", +"RPM1=1000;R=0.2;\n", +"\n", +"Ns=(120*F/P)\n", +"\n", +"S1=(Ns-RPM)/(Ns)\n", +" \n", +"S2=(Ns-RPM1)/(Ns)\n", +" \n", +"Rext=(R/S1)*S2-(R)\n", +"" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6: startingtorque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"F=60;P=4;N=1710;RI=6;\n", +"\n", +"Ns=(120*F)/P\n", +" \n", +"Sfl=(Ns-N)/Ns\n", +" \n", +"Tst=RI^2*Sfl\n", +" \n", +"Stmax=(0.0875/0.91)^(1/2)\n", +" \n", +"NMaxT=(1-Stmax)*Ns\n", +" \n", +"Tmax=(1+Stmax^2)/(2*Stmax)*Tst" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7: speedpower.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"F1=60;F2=15;F3=120;F=60;P=6;\n", +"V=240;A=2;\n", +"\n", +"S=F2/F1\n", +" \n", +"Ns=(120*F)/P\n", +" \n", +"N=(1+S)*Ns\n", +" \n", +"N1=(1-S)*Ns\n", +" \n", +"S1=F3/F1\n", +" \n", +"n=(1+S1)*Ns\n", +" \n", +"n1=(1-S1)*Ns\n", +" \n", +"Pac=1/S\n", +" \n", +"Pac1=-1/S\n", +"\n", +"Pdc=-(1-(S))/S\n", +" \n", +"Pdc=-(1+(S))/-S\n", +" \n", +"Pac1=1/S1\n", +" \n", +"Pac1=1/-S1\n", +"\n", +"Pdc=-(1-S1)/S1\n", +" \n", +"Pdc=-(1+S1)/-S1\n", +" " + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/6-SYNCHRONOUS_MACHINES.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/6-SYNCHRONOUS_MACHINES.ipynb new file mode 100644 index 0000000..75e221a --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/6-SYNCHRONOUS_MACHINES.ipynb @@ -0,0 +1,530 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 6: SYNCHRONOUS MACHINES" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: powerfactor.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"IMp=500;Pf=0.8;Pf1=0.6;\n", +"SMp=300;Pfs=1;P1=400;\n", +"P2=4;\n", +"\n", +"Power=IMp*Pf\n", +" \n", +"RecPower=IMp*Pf1\n", +"\n", +"SyPower=SMp*Pfs\n", +"\n", +"TotalPower=Power+SyPower\n", +" \n", +"RecPower\n", +" \n", +"ComplesPower=sqrt(TotalPower^2+RecPower^2)\n", +"\n", +"PowerFactor=(TotalPower/ComplesPower)\n", +" \n", +"Power1=sqrt(Power^2-RecPower^2)\n", +" \n", +"KVAR=RecPower-Power1\n", +" \n", +"NewKVA=sqrt(TotalPower^2+KVAR)\n", +"\n", +"PowerFactor1=(TotalPower/NewKVA)\n", +" \n", +"Ism=P1/(sqrt(3)*P2)\n", +"\n", +"PowerfactorSYS=SyPower/P1" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: unsaturatedvalues.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=14000;Phase=3;Ra=0.07;V1=10;Is=490;Pf=0.8;\n", +"If=200;Vl=18000;\n", +"\n", +"Vb=V/sqrt(Phase)\n", +" \n", +"Ib=(V1*10^6)/(sqrt(Phase)*V)\n", +" \n", +"Zb=Vb/Ib\n", +" \n", +"Zsun=(Vl/sqrt(3))/Is\n", +" \n", +"Xsun=sqrt(Zsun^2-Ra^2)\n", +" \n", +"Xsun=Xsun/Zb\n", +" \n", +"Zssa=(V/sqrt(3))/Is\n", +"\n", +"Xssa=Zssa/Zb\n", +"\n", +"a=cos(0.8)\n", +" \n", +"Deg=a*180/%pi\n", +" \n", +"Zs=Zssa/Ra\n", +" \n", +"Zs=atan(Zs)\n", +"\n", +"Angle=Zs*180/%pi\n", +"\n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[a,b]=polar2rect(1,0)\n", +" \n", +"X1=a+%i*b\n", +"\n", +"[c,d]=polar2rect(1,-Deg)\n", +" \n", +"X2=c+%i*d\n", +" \n", +"[e,f]=polar2rect(0.84,Angle)\n", +" \n", +"X3=e+%i*f\n", +" \n", +"X=X1+(X2*X3)\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +" \n", +"[I,Angle]=rect2polar(1.54,0.64)\n", +" \n", +"Ef=I*V\n", +" \n", +"If=I*If\n", +"\n", +" \n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: excitationvoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"V=208;Poles=4;F=60;Phase=3;Vol=5000;Xs=8;\n", +"\n", +"Vt=V/sqrt(Phase)\n", +" \n", +"Ia=Vol/(sqrt(Phase)*V)\n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"[x1,y1]=polar2rect(Vt,0)\n", +" \n", +"[x2,y2]=polar2rect(Ia,-36.9)\n", +" \n", +"[x3,y3]=polar2rect(8,90)\n", +" \n", +"X1=x1+%i*0\n", +" \n", +"X2=x2+%i*0\n", +" \n", +"X2=x2+%i*y2\n", +" \n", +"X3=x3+%i*y3\n", +" \n", +"X=X1+(X2*X3)\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ef,Angle]=rect2polar(186.7,88.7)\n", +"\n", +"Newvol=1.2*Ef\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[x1,y1]=polar2rect(Newvol,21)\n", +"\n", +"[x2,y2]=polar2rect(Vt,0)\n", +" \n", +"[x3,y3]=polar2rect(Xs,90)\n", +" \n", +"X1=x1+%i*y1\n", +" \n", +"X1=x2+%i*y2\n", +"\n", +"X1=x3+%i*y3\n", +"\n", +"X=(X1-X2)/X3\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ia,Angle1]=rect2polar(11.11,-13.93)\n", +"\n", +"Pf=cos(%pi*51.5/180)\n", +"\n", +"a=sin(%pi*51.5/180)\n", +" \n", +"ReactiveKva=(Phase*Vt*Ia*a*10^-3)\n", +" \n", +"Pmax=(Phase*Ef*Vt)/Xs\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[x1,y1]=polar2rect(206.9,90)\n", +" \n", +"[x2,y2]=polar2rect(120,0)\n", +" \n", +"[x3,y3]=polar2rect(8,90)\n", +" \n", +"X1=x1+%i*y1\n", +" \n", +"X2=x2+%i*y2\n", +"\n", +"X3=x3+%i*y3\n", +" \n", +"X=(X1-X2)/X3\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ia,Angle2]=rect2polar(25.8,15)\n", +"\n", +"tan=Vt/Ef\n", +"\n", +"Ang=tan*180/%pi\n", +"\n", +"Pf=cos(%pi*Ang/180)\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: powerangle.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=208;F=60;Phase=3;Power=3000;Xs=8;\n", +"\n", +"Vt=V/sqrt(Phase)\n", +" \n", +"Ia=Power/(Phase*Vt)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"[x1,y1]=polar2rect(120,0)\n", +"\n", +"[x2,y2]=polar2rect(8.33,0)\n", +" \n", +"[x3,y3]=polar2rect(8,90)\n", +" \n", +"X=X1-(X2*X3)\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ef,Angle]=rect2polar(120,-66.64)\n", +" \n", +"Pmax=(Phase*Ef*Vt)/Xs\n", +" \n", +"Ws=(1800/F)*2*%pi\n", +"\n", +"Tmax=Pmax/Ws" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5: fieldcurrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=460;Phase=3;N=1200;Hp=125;\n", +"Ra=0.078;Xal=0.15;Xar=1.85;Nre=28.2;Nse=28.2;\n", +"\n", +"Ia=sqrt(Phase)*V\n", +" \n", +"Ia=121.4\n", +" \n", +"Vt=V/sqrt(Phase)\n", +"\n", +"Ea=Vt-(Ia*Ra)\n", +" \n", +"Xs=Xal+Xar\n", +"\n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[x1,y1]=polar2rect(Ea,0)\n", +" \n", +"[x2,y2]=polar2rect(Xs,90)\n", +" \n", +"X1=x1+%i*y1\n", +" \n", +"X2=x2+%i*y2\n", +" \n", +"Im=X1/X2\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Im,Angle1]=rect2polar(7.84D-15,-128.0)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[x1,y1]=polar2rect(Im,Angle1)\n", +" \n", +"[x2,y2]=polar2rect(Ia,0)\n", +" \n", +"X1=x1+%i*y1\n", +" \n", +"X2=x2+%i*y2\n", +" \n", +"X=X1-X2\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[a,b]=rect2polar(-121.4,-128)\n", +" \n", +"n=sqrt(2)/Phase*Nre\n", +" \n", +"If=(a/n)*(Xs/Xar)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6: statorcurrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=11;Phase=3;F=60;IncExe=150;DecExe=50;Xs=10;\n", +"Power=80000;\n", +"\n", +"Vt=V/sqrt(Phase)\n", +" \n", +"Vt=Vt*1000\n", +" \n", +"Ef=IncExe/100\n", +"\n", +"Ef1=DecExe/100\n", +" \n", +"Ia=(Vt-(Ef*Vt))/(Xs)\n", +" \n", +"Pf=cos(90/90*%pi/2)\n", +" \n", +"Ia1=(Vt-(Ef1*Vt))/(Xs)\n", +" \n", +"Pf1=cos(90/90*%pi/2)\n", +" \n", +"Ia=Power/(Phase*Vt)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[x1,y1]=polar2rect(4.2,0)\n", +" \n", +"[x2,y2]=polar2rect(10,90)\n", +" \n", +"X1=x1+%i*y1\n", +"\n", +"X2=x2+%i*y2\n", +" \n", +"X=X1*X2\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[a,b]=rect2polar(Vt,42)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7: powercurrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Vt=1;Ia=1;Xd=0.8;Xq=0.4;Loss=0.15;Angle=36.9;\n", +"\n", +"a=cos(%pi*Angle/180)\n", +"\n", +"b=sin(%pi*Angle/180)\n", +"\n", +"TanDeg=(Vt*Xq*a)/(Vt+(Ia*Xq*b))\n", +" \n", +"z=atan(TanDeg)\n", +" \n", +"Deg=(z*%pi/180)\n", +" \n", +"Deg=(z*180/%pi)\n", +"\n", +"Angl=Angle-Deg\n", +" \n", +"Id=Ia*sin(%pi*Angl/180)\n", +" \n", +"Iq=Ia*cos(%pi*Angl/180)\n", +" \n", +"Ef=(Vt*cos(%pi*Deg/180))-(Id*Xq)\n", +" \n", +"Pf=((Vt*Ef)/Xd)*sin(%pi*Deg/180)\n", +" \n", +"Pr=(Vt^2*(Xd-Xq)/(2*Xd*Xq))*sin(%pi*45.6/180) \n", +" \n", +"Pout=Vt*Ia*a\n", +" \n", +"Prmax=(Vt^2*(Xd-Xq))/(2*Xd*Xq)\n", +" \n", +"v=asin(Loss/Prmax)/2\n", +" \n", +"Deg1=(v*180/%pi)\n", +" \n", +"Id=Vt*cos(%pi*Deg1/180)/Xd\n", +" \n", +"Iq=Vt*sin(%pi*Deg1/180)/Xq\n", +" \n", +"Ia=(Id^2+Iq^2)^(1/2)\n", +" \n", +" Ang=atan(Id/Iq)\n", +" \n", +" Ang=(Ang*180/%pi)\n", +" \n", +"Phi=Ang+Deg1\n", +" \n", +"Pf=cos(%pi*Phi/180)" + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/7-SINGLE_PHASE_MOTORS.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/7-SINGLE_PHASE_MOTORS.ipynb new file mode 100644 index 0000000..e24879d --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/7-SINGLE_PHASE_MOTORS.ipynb @@ -0,0 +1,504 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 7: SINGLE PHASE MOTORS" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 2: ouputpower.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"F=60;P=4;N=1730;Zb=27.86;\n", +"R=1.35;X=1.63;R1=2.9;X1=3.26;\n", +"V=120;Prot=81.2;\n", +"\n", +"Ns=(120*F)/P\n", +" \n", +"S=(Ns-N)/N\n", +" \n", +"a=%i*Zb*((R/S)+%i*X)\n", +" \n", +"b=(R/S)+%i*(Zb+X)\n", +" \n", +"Zf=a/b\n", +" \n", +"Rf=13.06;\n", +" \n", +"Xf=16.31;\n", +" \n", +"a=%i*Zb*(R/(2-S)+%i*X)\n", +" \n", +"b=R/(2-S)+%i*(Zb+X)\n", +"\n", +"Zb=a/b\n", +" \n", +"Rb=0.61;\n", +"\n", +"Xb=1.55;\n", +" \n", +"Zinput=(R1+Rf+Rb)+%i*(X1+Xf+Xb)\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[a,b]=rect2polar(16.57,21.12)\n", +"\n", +"Iinput=V/a\n", +" \n", +"cos(%pi*b/180)\n", +" \n", +"Pinput=V*Iinput*ans\n", +" \n", +"Wsy=Ns*2*(%pi/F)\n", +"\n", +"T=Iinput^2*(Rf-Rb)/Wsy\n", +" \n", +"Pmech=T*Wsy*(1-S)\n", +"\n", +"OutputPower=Pmech-Prot\n", +"\n", +"Eff=OutputPower/Pinput*100\n", +" \n", +"Pgf=Iinput^2*Rf\n", +" \n", +"Pgb=Iinput^2*Rb\n", +" \n", +"airgap=Pgf+Pgb\n", +" \n", +"P2=S*Pgf+(2-S)*Pgb" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: externalresistance.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=120;F=60;Pole=4;Zm=1.5+4.0;Za=3+6;\n", +"Xa=6;Xm=4;Rm=1.5;Ra=3;\n", +"\n", +"Ra=(Xa/Xm)*(Rm+sqrt(18.25))\n", +" \n", +"C=(2*%pi*F)*(Xa+(Ra*Rm)/(Xm+sqrt(18.25)))\n", +" \n", +"a=((-Xm*Ra)+(sqrt(18.25)*sqrt(13.2)))\n", +" \n", +"Xc=Xa+(a/Rm)\n", +" \n", +"Ia=V/(3+%i*6)\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is,Angle]=rect2polar(8,-16)\n", +"\n", +"Im=V/(1.5+%i*4)\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is1,Angle1]=rect2polar(9.86,-26.3)\n", +" \n", +"Alfa=Angle1-Angle\n", +" \n", +"Ts=Is*sin(%pi*6.01/180)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[a,b]=polar2rect(Is1,Angle1)\n", +" \n", +"X=a+%i*b\n", +"\n", +"C=1/C*10^6\n", +"\n", +"a=((-Xm*Ra)+(sqrt(18.25)*sqrt(13.2)))\n", +" \n", +"Xc=Xa+(a/Rm)\n", +" \n", +"C=10^6/(2*%pi*F*Xc)\n", +" \n", +"Ia=V/(3+%i*6)\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is,Angle]=rect2polar(8,-16)\n", +"\n", +"Im=V/(1.5+%i*4)\n", +" \n", +"[Is1,Angle1]=rect2polar(9.86,-26.3)\n", +"\n", +"Alfa=Angle1-Angle\n", +"\n", +"Ts=Is*sin(%pi*6.01/180)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[a,b]=polar2rect(Is1,Angle1)\n", +" \n", +"X=a+%i*b\n", +" \n", +"[c,d]=polar2rect(Is,Angle)\n", +"\n", +"X1=c+%i*d\n", +" \n", +"X2=X+X1\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[I,Angle]=rect2polar(17.86,-42.3)\n", +" \n", +"Ia=V/(Ra+%i*Xa)\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ia,Angle]=rect2polar(9.3,-6.4)\n", +" \n", +"Alfa=69.33-34.53\n", +"\n", +"Ts=Ia*sin(%pi*Alfa/180)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[Is,Angle]=polar2rect(Ia,Angle)\n", +" \n", +"[Is1,Angle1]=polar2rect(28.1,-69.44)\n", +" \n", +"X=Is+%i*angle\n", +"\n", +"X1=Is1+%i*Angle1\n", +" \n", +"X2=Is+%i*Angle\n", +"\n", +"X=X1+X2\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is,Angle]=rect2polar(19.1,-32.7)\n", +"\n", +"Xc=10^6/(2*%pi*F*405)\n", +"\n", +"Ia=V/(Ra+(%i*6+%i*6.55))\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is,Angle]=rect2polar(2.16,-9.04)\n", +" \n", +"\n", +"Ia=V/(Ra+(%i*6-%i*6.55))\n", +" \n", +"[Is,Angle]=rect2polar(38.6,7.09)\n", +"\n", +"Alfa=69.44+Angle\n", +"\n", +"Ts=Is*sin(%pi*Alfa/180)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[Is,Angle]=polar2rect(28.1,-69.44)\n", +" \n", +"[Is1,Angle1]=polar2rect(39.34,10.4)\n", +"\n", +"X1=Is+%i*Angle\n", +" \n", +"X2=Is1+%i*Angle1\n", +"\n", +"X=X1+X2\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is,Angle]=rect2polar(48.56,-19.20)\n", +" \n", +"Ia=V/(Ra+(%i*Xa-%i*Xc))\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[I,Angle]=rect2polar(23.9,19.6)\n", +"\n", +"Alfa=69.44+39.5\n", +"\n", +"Ts=I*sin(%pi*Alfa/180)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[Is,Angle]=polar2rect(28.1,-69.44)\n", +" \n", +"[Is1,Angle1]=polar2rect(I,39.35)\n", +" \n", +"X=Is+%i*Angle\n", +" \n", +"X1=Is1+%i*Angle1\n", +" \n", +"X2=X+X1\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[I,Angle]=rect2polar(33.7,-6.7)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 5: maximumtorque.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=120;F=60;X1m=2;R1m=1.5;R2=1.5;\n", +"X1a=2;R1a=1.5;X2=2;Xmag=48;C=30;a=1;\n", +"Z1m=1.5;Zb=0.69+%i*0.98;Z1a=2.5;\n", +"Xc=%i*2-%i*88.4;Ra=2.5;\n", +"\n", +"Xc=10^6/(2*%pi*F*C)\n", +" \n", +"Zb\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[x,y]=rect2polar(0.69,0.98)\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[a,b]=polar2rect(V,0)\n", +" \n", +"X=a+%i*b\n", +"\n", +"z=(Z1m+%i*2+2*(Zb))\n", +" \n", +"Im=X/z\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is,Angle]=rect2polar(14.41,-19.81)\n", +"\n", +"y=(Z1a+Xc+2*(Zb))\n", +" \n", +"Ia=X/y\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is1,Angle1]=rect2polar(0.065,1.41)\n", +" \n", +"Wsy=(1800*2*%pi)/F\n", +" \n", +"Ts=2*(Is)*(Is1)*2*0.69*sin(%pi*141.1/180)/Wsy\n", +" \n", +"Zm=Z1m+%i*2+2*(Zb)\n", +" \n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Ip1,Angle1]=rect2polar(2.88,3.96)\n", +" \n", +"Za=Ra+%i*2+2*(Zb)\n", +" \n", +"R=3.88;Im=3.96;\n", +"Xc=Im-((Im*R-4.9*sqrt(26.22))/2.88)\n", +"\n", +"c=10^6/(2767.34)\n", +"\n", +"Cs=c-C\n", +" \n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[v,a]=polar2rect(V,0)\n", +" \n", +"X=v+%i*a\n", +"\n", +"R=3.88;Im=3.96;Xc=7.34;\n", +"\n", +"a=R+(%i*Im-%i*Xc)\n", +"\n", +"z=X/a\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[Is1,Angle1]=rect2polar(17.5,15.3)\n", +" \n", +"Is=24.4;Angle=-53.4;\n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[a,b]=polar2rect(Is,Angle)\n", +" \n", +"X1=a+%i*b\n", +" \n", +"[c,d]=polar2rect(Is1,Angle1)\n", +"\n", +"X2=c+%i*d\n", +" \n", +"X=X1+X2\n", +"\n", +"function[r,theta]=rect2polar(x,y)\n", +"r=sqrt(x^2+y^2);\n", +"theta=atan(y/x)*180/%pi;\n", +"endfunction\n", +"\n", +"[z,y]=rect2polar(32.04,-4.28)\n", +" \n", +"a=sin(%pi*95/180)\n", +"\n", +"Ts=2*(Is1)*(Is)*2*0.69*a/Wsy\n", +"\n", +"T=Ts/z" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 6: torquedev.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"V=120;F=60;N=2000;A=0.6;\n", +"Ohm=20;L=0.25;\n", +"V=120;F=60;N=2000;Ia=0.6;\n", +"a=20;L=0.25;\n", +"\n", +"Edc=V-(Ia*Ra)\n", +" \n", +"X=2*%pi*F*L\n", +" \n", +"Eac=(-Ia*Ra)+sqrt(V^2-(Ia*X)^2)\n", +" \n", +"Nac=N*(Eac/Edc)\n", +" \n", +"Pf=(Eac+(Ia*Ra))/V\n", +"\n", +"Pmech=Eac*Ia\n", +" \n", +"Wm=(Nac*2*%pi)/F\n", +"\n", +"T=Pmech/Wm" + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/8-SPECIAL_MACHINES.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/8-SPECIAL_MACHINES.ipynb new file mode 100644 index 0000000..ef4799e --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/8-SPECIAL_MACHINES.ipynb @@ -0,0 +1,141 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 8: SPECIAL MACHINES" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: shaftpositio.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"T=0.2;V=115;N=3000;F=60;\n", +"J=10^-5;\n", +"\n", +"Km=T/V\n", +" \n", +"Wm=(N*2*%pi/F)\n", +" \n", +"Fm=T/Wm\n", +" \n", +"Tm=J/Fm\n", +" \n", +"A=Km/Fm\n", +" \n", +"Kmv=A*V\n", +"\n", +"KmvT=A*Tm" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: maximumvoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Rw=1;Lw=30;I=3;TimeOn=2;RF=0.0675;\n", +"StepRate=300;Turns=100;TimeOff=1;\n", +"PeakI=3;\n", +"\n", +"R=Lw/TimeOn\n", +" \n", +"Rext=R-TimeOff\n", +" \n", +"Prext=(I^2*Rext)\n", +" \n", +"Vs=I*R\n", +"\n", +"Rext=R-Rw\n", +" \n", +"R1=Lw/TimeOff\n", +" \n", +"Rf=R1-R\n", +"\n", +"Energy=(1/2*Lw*I^2)\n", +" \n", +"Power=Turns*Rf\n", +" \n", +"Power=Turns*RF\n", +" \n", +"Vc=V+(PeakI*R)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4: inductionenergy.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Lw=30;R=15;Ia=3;V=45;\n", +"\n", +"Tow=Lw/R\n", +"\n", +"t1=0.7*Tow\n", +"\n", +"t0=0:0.1:t1;\n", +" \n", +"t=0;\n", +" \n", +"a=integrate('45*(-3+6*%e^(-x/2))','x',t,t0)\n", +" \n", +" Energy=(1/2)*Lw*Ia^2\n", +" \n", +" ProEnergy=(a/Energy)*100" + ] + } +], +"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/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/9-TRANSIENTS_AND_DYNAMICS.ipynb b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/9-TRANSIENTS_AND_DYNAMICS.ipynb new file mode 100644 index 0000000..81f46ec --- /dev/null +++ b/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/9-TRANSIENTS_AND_DYNAMICS.ipynb @@ -0,0 +1,131 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 9: TRANSIENTS AND DYNAMICS" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 1: armaturevoltage.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Rf=100;Lf=25;Ra=0.25;Laq=0.02;\n", +"Kg=100;Ll=0.15;V=200;Rl=1;\n", +"\n", +"tow=Lf/Kg\n", +"\n", +"log(0.1)\n", +"\n", +"t=2.30/4\n", +"\n", +"Towat=(Ll+Laq)/(Rl+Ra)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 3: maximumcurrent.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"Xd=0.9;Vt=1;Ia=1;Xd1=0.4;Xd2=0.2;Ta=0.2;\n", +"Td1=4;Td2=0.6;t=0.1;\n", +"\n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[x,y]=polar2rect(Ia,-25.8)\n", +"\n", +"X=x+%i*y\n", +"\n", +"[x1,y1]=polar2rect(Xd,90)\n", +"\n", +"X1=x1+%i*y1\n", +"\n", +"A=Vt+(X*X1)\n", +"\n", +"function[x,y]=polar2rect(r,theta)\n", +"x=r*cos(theta*%pi/180);\n", +"y=r*sin(theta*%pi/180);\n", +"endfunction \n", +"\n", +"[Ei,Angle]=rect2polar(1.39,0.81)\n", +"\n", +"[x,y]=polar2rect(Ia,-25.8)\n", +"\n", +"X=x+%i*y\n", +"\n", +"[x2,y2]=polar2rect(Xd1,90)\n", +"\n", +"X2=x2+%i*y2\n", +"\n", +"[Ei2,Angle1]=rect2polar(1.17,0.36)\n", +"\n", +"[x,y]=polar2rect(Ia,-25.8)\n", +"\n", +"X=x+%i*y\n", +"\n", +"[x3,y3]=polar2rect(Xd2,90)\n", +"\n", +"X3=x3+%i*y3\n", +"\n", +"[Ei2,Angle2]=rect2polar(1.08,0.10)\n", +"\n", +"Idc=sqrt(2)*(Ei2/Xd2)\n", +"\n", +"Td1=(Xd1/Xd)*Td1\n", +"\n", +"Td2=(Xd2/Xd1)*Td2\n", +"\n", +"Isc=sqrt(2)*(Td1+1.29*%e^(-0.562*0.1)+2.42*%e^(-3.3*0.1))+7.78*%e^(-5*0.1)\n", +"\n", +" " + ] + } +], +"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 +} |