initial commit / add all books

1 files changed, 40 insertions, 0 deletions

diff --git a/3760/CH5/EX5.19/Ex5_19.sce b/3760/CH5/EX5.19/Ex5_19.sce
new file mode 100644
index 000000000..ca0b5df91
--- /dev/null
+++ b/3760/CH5/EX5.19/Ex5_19.sce
@@ -0,0 +1,40 @@
+clc;

+v=2000; // rated voltage of motor

+xsm=2; // synchronous reactance of motor

+xsg=3; // synchronous reactance of generator

+xt=1.5; // transmission line reactance

+ia=100; // current drawn by motor

+pf=1; // power factor

+disp('case a');

+vt=v/sqrt(3); // rated per phase voltage

+Efm=floor(sqrt(vt^2+(ia*xsm)^2)); // excitation EMF

+printf('Excitation EMF for motor is %f V\n',Efm);

+Efg=sqrt(vt^2+(ia*(xsg+xt))^2); // excitation EMF

+printf('Excitation EMF for alternator is %f V\n',Efg);

+disp('case b');

+de1=acosd(vt/Efm); // load angle for motor

+de2=acosd(vt/Efg); // load angle for alternator

+de=de1+de2; // power angle between Efm and Efg

+pt=(Efg*Efm*sind(de))/(xsm+xsg+xt);

+P=pt*3;

+printf('Per phase power transfer between alternator and motor is %f KW\n',pt/1000);

+printf('Net power transfer between alternator and motor is %f KW\n',P/1000);

+disp('case c');

+// from phasor diagram fig 5.54

+ia=sqrt(Efm^2+Efg^2)/(xsm+xsg+xt); 

+// for maximum transfer of power , power angle=90 degrees

+de=90

+pmax=(Efg*Efm*sind(de))/(xsm+xsg+xt);

+P=pmax*3;

+printf('Per phase maximum power transfer between alternator and motor is %f KW\n',pmax/1000);

+printf('Net maximum power transfer between alternator and motor is %f KW\n',P/1000);

+// from phasor diagrams determining various parameters needed to find power factor

+be=acosd(Efm/(ia*(xsm+xsg+xt)));

+Vp=sqrt((Efm-ia*xsm*cosd(be))^2+(ia*xsm*sind(be))^2); // phase voltage

+Vl=sqrt(3)*Vp; // line voltage

+printf('Armature current for given condition is %f A\n',ia);

+printf('Terminal voltage of synchronous motor is %f V\n',Vp);

+// from phasor diagram

+aoc=asind((ia*xsm*sind(be))/Vp);

+pf=cosd(90-be-aoc); 

+printf('Power factor angle of motor is %f leading',pf);