summaryrefslogtreecommitdiff
path: root/Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen
diff options
context:
space:
mode:
Diffstat (limited to 'Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen')
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/1-magnetic_circuit.ipynb235
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/10-POWER_SEMICONDUCTOR_CONVERTERS.ipynb356
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/2-TRANSFORMERS.ipynb454
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/3-ELECTROMECHNICAL_ENERGY_CONVERSION.ipynb159
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/4-DC_MACHINES.ipynb512
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/5-INDUCTION_ASYNCHRONOUS_MACHINES.ipynb501
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/6-SYNCHRONOUS_MACHINES.ipynb530
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/7-SINGLE_PHASE_MOTORS.ipynb504
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/8-SPECIAL_MACHINES.ipynb141
-rw-r--r--Principles_Of_Electric_Machines_And_Power_Electronics_by_P_C_Sen/9-TRANSIENTS_AND_DYNAMICS.ipynb131
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
+}