1 files changed, 54 insertions, 0 deletions
diff --git a/3793/CH5/EX5.3/exp_5_3.sce b/3793/CH5/EX5.3/exp_5_3.sce
new file mode 100644
index 000000000..c59f88723
--- /dev/null
+++ b/3793/CH5/EX5.3/exp_5_3.sce
@@ -0,0 +1,54 @@
+clear;
+clc;
+function [Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
+    index=1;
+    if index==1;
+        X= input('percentage load');
+        Pg=X*Prtd;
+        Vt=Vt/sqrt(3);
+        pfang=Phi*%pi/180;
+        x=pfang:-.005:-pfang;
+        Pg=Pg*ones(1,length(x));
+        Ia=Pg./(3*abs(Vt)*cos(x));
+        Iacom=Ia.*(cos(x)+%i*sin(x));
+        Ecom=Vt+%i*Xs*Iacom;
+        Emag=abs(Ecom);
+        Ifld=Emag*1000/K;
+        plot(Ifld,Ia,'k');
+        
+        set(gca(),"auto_clear","off");
+        return;
+        
+        plot(Ifld,Ia,'k-');
+        return;
+       
+        plot(Ifld,Ia,'k-.');
+        return;
+        
+        plot(Ifld,Ia,'k.');
+        return;
+        
+        set(gca(),"auto_clear","on");
+        
+       
+        
+        else
+     
+    end
+endfunction
+Prtd=36;
+Vt=33;
+Xs=10;
+Phi=70;
+K=1500;
+
+[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
+[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
+[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
+[Ecom,Emag]=vcurves(Prtd,Xs,Vt,Phi,K);
+index=0;
+
+xlabel('Field current in amperes');
+ylabel('Armature current in amp');
+title('Plot of V-curves of a Synchronous machine');
+