summaryrefslogtreecommitdiff
path: root/3760/CH6/EX6.10/Ex6_10.sce
diff options
context:
space:
mode:
Diffstat (limited to '3760/CH6/EX6.10/Ex6_10.sce')
-rw-r--r--3760/CH6/EX6.10/Ex6_10.sce78
1 files changed, 78 insertions, 0 deletions
diff --git a/3760/CH6/EX6.10/Ex6_10.sce b/3760/CH6/EX6.10/Ex6_10.sce
new file mode 100644
index 000000000..634d1fc96
--- /dev/null
+++ b/3760/CH6/EX6.10/Ex6_10.sce
@@ -0,0 +1,78 @@
+
+
+clc;
+
+//from 6.9 problem
+P=4;
+r1=0.15;
+x1=0.45;
+r2=0.12;
+x2=0.45;
+Xm=28.5;
+s=0.04;
+V=400;
+f=50;
+Pfixed=400;
+t=1.2; // rotor effective turns ratio
+
+//for part a
+//According to the conditions and diagram
+t1=complex(r1,x1);
+t2=complex(0,Xm);
+t3=complex(r1,x2+Xm);
+Ze=(t1*t2)/(t3);
+Re=real(Ze);
+Xe=imag(Ze);
+t4=complex(Re,(x2+Xe));
+SmT=(r2)/(sqrt((Re*Re)+((x2+Xe)*(x2+Xe))));
+Ve=(V/sqrt(3))*(Xm/(x2+Xm));
+Ws=(4*%pi*f)/P;
+Tem=(3/Ws)*Ve^2*(1/2)*(1/(Re+sqrt(Re^2+(x2+Xe)^2)));
+Pm=Tem*(1-SmT)*Ws;
+Psh=Pm-Pfixed;
+Tsh=Psh/(Ws*(1-SmT));
+mprintf('for part a \n slip = %f \n maximun torque = %f Nm \n power output = %f KW \n',SmT,Tem, Psh/1000);
+
+
+//for part b
+s=1;
+I2st=(Ve)/(sqrt((r2+Re)*(r2+Re)+(x2+Xe)*(x2+Xe)));
+Test=(3/Ws)*I2st*I2st*(r2);
+mprintf(' for part b rotor current = %f A \n torque = %f Nm \n',I2st,Test);
+
+
+//for part c
+R=sqrt(Re^2+(x2+Xe)^2)-r2;
+Ra=R/(t^2);
+mprintf('for part c \n external resisitance value is = %f Ohm \n',Ra);
+
+//for part d
+s1=0.04;
+Pm=((3*(Ve)*(Ve))*r2*((1-s1)/s1))/(((Re+r2+((r2*(1-s1)/s1))))*((Re+r2+((r2*(1-s1)/s1))))+((x2+Xe)*(x2+Xe)));
+mprintf('for part d \n power developed is %f KW \n',Pm/1000);
+
+//for part e
+SmP=(r2)/(sqrt(((Re+r2)*(Re+r2))+((x2+Xe)*(x2+Xe)))+r2);
+Pmn=((3/2)*Ve*Ve)/(Re+r2+sqrt((r2+Re)*(r2+Re)+(x2+Xe)*(x2+Xe)));
+mprintf('for part e \n slip = %f \n power developed = %f KW',SmP,Pmn/1000);
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+