diff options
Diffstat (limited to '3760/CH5/EX5.24/Ex5_24.sce')
| -rw-r--r-- | 3760/CH5/EX5.24/Ex5_24.sce | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/3760/CH5/EX5.24/Ex5_24.sce b/3760/CH5/EX5.24/Ex5_24.sce new file mode 100644 index 000000000..95035647e --- /dev/null +++ b/3760/CH5/EX5.24/Ex5_24.sce @@ -0,0 +1,26 @@ +clc;
+xd=1.2; // d axis synchronous reactance
+xq=0.8; // q axis synchronous reactance
+ra=0.025; // armature resistance
+vt=1; // pu rated per phase voltage
+disp('case a');
+disp('For lagging power factor')
+pf=0.8; // power factor
+ia=1*(pf-sqrt(1-pf^2)*%i); // armature current
+Ef1=vt+%i*(ia*xq)+ia*ra; // excitation EMF
+id=1*sind(atand(imag(Ef1),real(Ef1))+acosd(pf)); // d component of armature current
+iq=1*cosd(atand(imag(Ef1),real(Ef1))+acosd(pf)); // q component of armature current
+Ef=abs(Ef1)+id*(xd-xq); // excitation EMF
+printf('Excitation EMF is %f p.u. at a load angle of %f degrees\n',abs(Ef),atand(imag(Ef1),real(Ef1)));
+disp('case b');
+disp('For leading power factor')
+pf=0.8; // power factor
+ia=1*(pf+sqrt(1-pf^2)*%i); // armature current
+Ef1=vt+%i*(ia*xq)+ia*ra; // excitation EMF
+id=1*sind(atand(imag(Ef1),real(Ef1))-acosd(pf)); // d component of armature current
+iq=1*cosd(atand(imag(Ef1),real(Ef1))-acosd(pf)); // q component of armature current
+Ef=abs(Ef1)+id*(xd-xq); // excitation EMF
+printf('Excitation EMF is %f p.u.at a load angle of %f degrees\n',abs(Ef),atand(imag(Ef1),real(Ef1)));
+
+
+
|