1 files changed, 37 insertions, 0 deletions
diff --git a/3760/CH5/EX5.35/Ex5_35.sce b/3760/CH5/EX5.35/Ex5_35.sce
new file mode 100644
index 000000000..592ace7a2
--- /dev/null
+++ b/3760/CH5/EX5.35/Ex5_35.sce
@@ -0,0 +1,37 @@
+clc;
+v=11000; // rated voltage of motor
+P=20*10^6; // rated power of motor
+p=12; // number of poles
+f=50; // frequency 
+xd=5; // d-axis synchronous reactance
+xq=3; // q-axis synchronous reactance
+vt=v/sqrt(3); // per phase rated voltage
+ia=P/(sqrt(3)*v); // per phase armature current
+disp('case a');
+// from phasor diagram
+Ef=vt-%i*(ia*xq);
+de=atand((ia*xq),(vt)); // load angle
+id=ia*sind(de); // d-axis component of armature current
+Ef=ceil(abs(Ef)+id*(xd-xq));
+printf('Excitation voltage is %f V to neutral\n',Ef);
+disp('case b');
+po=(((Ef*vt*sind(de))/xd)+((vt^2*(xd-xq)*sind(2*de))/(2*xd*xq)))/1000;
+printf('Power is %f KW per phase\n',po);
+po=3*po*1000; 
+printf('Calculated power %f W is almost equal to given rated power i.e %f W\n',po,P);
+disp('case c'); 
+ps=((Ef*vt*cosd(de))/xd)+((vt^2*(xd-xq)*cosd(2*de))/(xd*xq));
+printf('Synchronozing power per electrical degree is %f KW\n',(3*ps*%pi)/(180*1000));
+ws=(4*%pi*f)/p; // synchronous speed
+t=(3*ps*%pi)/(180*ws);
+printf('Torque corresponding to synchronous power is %f N-m\n',t);
+disp('case d');
+printf('Synchronozing power per mechanical degree is %f KW\n',(3*ps*%pi*p)/(2*180*1000));
+t=(3*ps*%pi*p)/(2*180); 
+printf('Torque corresponding to synchronous power is %f N-m\n',t);
+disp('case e');
+de=acosd((-Ef*xq)/(4*Ef*(xd-xq))+(sqrt(0.5+((Ef*xq)/(4*vt*(xd-xq)))^2)));
+printf('Maximum value of power angle is %f degrees\n',de);
+pmax=(((Ef*vt*sind(de))/xd)+((vt^2*(xd-xq)*sind(2*de))/(2*xd*xq)))/1000;
+printf('Maximum power is %f KW per phase\n',pmax);
+printf('Maximum power is %f KW for 3 phase\n',pmax*3);