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.7/Ex23_7.sce | 45 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 45 insertions(+)
 create mode 100755 2252/CH23/EX23.7/Ex23_7.sce

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

diff --git a/2252/CH23/EX23.7/Ex23_7.sce b/2252/CH23/EX23.7/Ex23_7.sce
new file mode 100755
index 000000000..c9c73c3ab
--- /dev/null
+++ b/2252/CH23/EX23.7/Ex23_7.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=3/100//slip
+Rl=R2_dash*(1-s)/s//load resistance
+Z1=.15+.4*j//stator impedance
+Z2_dash=.16+.4*j//equivalent rotor impedance
+Z=Rl+Z1+Z2_dash//total impedance per phase
+Vph=400/sqrt(3)//applied voltage per phase
+I2_dash=Vph/Z
+R0=200
+Xm=20
+Iw=Vph/R0
+Im=Vph/Xm
+I0=Iw-Im*j
+I1=I0+I2_dash
+[I1 theta]=rect2pol(I1)
+pf=cos(theta*%pi/180)
+p=mag(I2_dash)^2*Rl//output power per phase
+Pout=3*p
+Pin=3*Vph*I1*pf
+mprintf("By using approximate equivalent circuit, the values of different parameters are as under\nPer phase stator current=%f A\nPer phase rotor current=%f A\nOperating power factor=%f\nInput power=%f kW\nMechanical output power=%f kW\n",I1,I2_dash,pf,Pin/10^3,Pout/10^3)
+//refer Fig. 23.16 in the textbook
+Z2_dash=5.3+j*0.4
+Z=(j*Xm)*Z2_dash/(j*Xm+Z2_dash)//equivalent impedance
+Zin=Z1+Z//total input impedance
+I1=Vph/Zin
+I2_dash=I1*j*Xm/(j*Xm+Z2_dash)
+[I1 theta]=rect2pol(I1)
+pf=cos(theta*%pi/180)
+Pout=3*(mag(I2_dash))^2*R2_dash*(1-s)/s
+Pin=3*Vph*I1*pf
+mprintf("By solving the problem using the circuit in Fig. 23.16 in the book, the values of different parameters are as under\nPer phase stator current=%f A\nPer phase rotor current=%f A\nOperating power factor=%f\nInput power=%f kW\nMechanical output power=%f kW\n",mag(I1),I2_dash,pf,Pin/10^3,Pout/10^3)
-- 
cgit