From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2252/CH23/EX23.6/Ex23_6.sce | 45 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 45 insertions(+)
 create mode 100755 2252/CH23/EX23.6/Ex23_6.sce

(limited to '2252/CH23/EX23.6')

diff --git a/2252/CH23/EX23.6/Ex23_6.sce b/2252/CH23/EX23.6/Ex23_6.sce
new file mode 100755
index 000000000..7f6e6c72f
--- /dev/null
+++ b/2252/CH23/EX23.6/Ex23_6.sce
@@ -0,0 +1,45 @@
+
+function[r,theta]=rect2pol(A)
+    x=real(A)
+    y=imag(A)
+    r=sqrt(x^2+y^2)
+    theta=atand(y/x)
+endfunction
+function[r]=mag(A)
+    x=real(A)
+    y=imag(A)
+    r=sqrt(x^2+y^2)
+endfunction
+j=%i
+R2_dash=.16
+s=.03
+X2=.4
+Z2_dash=R2_dash/s+X2*j//effective rotor impedance referred to stator
+R0=200
+Xm=20*j
+Z=1/(1/R0+1/Xm+1/Z2_dash)//equivalent impedance
+Z1=.15+.4*j//stator impedance
+Zin=Z1+Z//total input impedance
+V=400//applied voltage
+//calculating stator current
+V1=V/sqrt(3)//per phase stator voltage
+I1=V1/Zin
+[I1 theta1]=rect2pol(I1)
+mprintf("Stator current=%f A at %f power factor lagging\n",I1,cos(theta1*%pi/180))
+//calculating rotor current
+I1=V1/Zin
+E1=V1-I1*Z1
+Iw=E1/R0//per phase core loss component of no load current
+Im=E1/Xm//per phase magnetising current
+I0=Iw+Im//no load current
+I2_dash=I1-I0
+[I2_dash theta2]=rect2pol(I2_dash)
+mprintf("Per phase rotor current=%f A lagging by %f degrees\n",I2_dash,-theta2)
+//calculating mechanical output power
+P=I2_dash^2*R2_dash*(1-s)/s//mechanical power output per phase
+Pout=3*P
+mprintf("Total mechanical power output=%f kW\n",Pout/1000)
+//calculating input power drawn by the motor
+Pin=3*V1*mag(I1)*cos(theta1*%pi/180)
+mprintf("Total input power drawn by the motor=%f kW\n",Pin/1000)
+//answers vary from the textbook due to round off error
-- 
cgit