initial commit / add all books

1 files changed, 45 insertions, 0 deletions

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)