initial commit / add all books

9 files changed, 260 insertions, 0 deletions

diff --git a/3793/CH5/EX5.1/exp_5_1.sce b/3793/CH5/EX5.1/exp_5_1.sce
new file mode 100644
index 000000000..b0ff219b9
--- /dev/null
+++ b/3793/CH5/EX5.1/exp_5_1.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+S=250;

+V=3330;

+r=.25;

+x=3.5;

+pf=0.8;

+i=(S*1000)/(sqrt(3)*V);

+Il=i;

+Vp=(V/sqrt(3));

+z=complex(r,x);

+Ia=complex(34.99,-26.24);

+E=Vp+(Ia*z);

+

+Vv=sqrt((real(E)^2)+(imag(E)^2));

+reg=((Vv-Vp)/Vp)*100;

+

+mprintf("regulation=%f percent",reg);

diff --git a/3793/CH5/EX5.2/exp_5_2.sce b/3793/CH5/EX5.2/exp_5_2.sce
new file mode 100644
index 000000000..64302f117
--- /dev/null
+++ b/3793/CH5/EX5.2/exp_5_2.sce
@@ -0,0 +1,42 @@
+clear;

+clc;

+V=33;

+S=45;

+pf=0.8;

+x=10;

+P=S*pf;

+z=complex(0,10);

+I=(P*1000)/(sqrt(3)*V*pf);

+Vp=V/sqrt(3);

+Ia=I*(complex(0.8,-0.6));

+E=(Vp*1000)+(z*Ia);

+mprintf("Excitation voltage per phase %.3f + %.3f\n",real(E),imag(E));

+Vv=sqrt((real(E)^2)+(imag(E)^2));

+del=asind((10*10)/((Vv*10^(-3))*Vp));

+x=Vv*cosd(del);

+y=Vv*sind(del);

+VV=complex(x,y);

+II=(VV-(Vp*1000))/complex(0,10);

+pf1=cosd(atand(imag(II)/real(II)));

+mprintf("Armature current  %.3f + %.3f\n",real(II),imag(II));

+mprintf("load angle %f\n",del);

+mprintf("power factor %f\n",pf1);

+maxdel=90;

+x1=Vv*cosd(maxdel);

+y1=Vv*sind(maxdel);

+VV1=complex(x1,y1);

+II1=(VV1-(Vp*1000))/complex(0,10);

+pf2=cosd(atand(imag(II1)/real(II1)));

+

+Pmax=(3*Vp*(Vv*10^(-3)))/10;

+mprintf("Maximum power %f",Pmax);

+p=P/3; // at minimum excitation pf is unity there fore power  is per phase

+

+E1=(p*10)/Vp;

+x2=E1*cosd(maxdel);

+y2=E1*sind(maxdel);

+VV2=complex(x2,y2);

+II2=((VV2-(Vp))*1000)/complex(0,10);

+pf3=cosd(atand(imag(II2)/real(II2)));

+

+

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');

+

diff --git a/3793/CH5/EX5.4/exp_5_4.sce b/3793/CH5/EX5.4/exp_5_4.sce
new file mode 100644
index 000000000..cf76dda7b
--- /dev/null
+++ b/3793/CH5/EX5.4/exp_5_4.sce
@@ -0,0 +1,39 @@
+clear;

+clc;

+E=1;

+Ig=.8;

+pf=.8;

+Xd=1.1;

+Xq=.8;

+p=E+complex(0,Xq)*Ig*complex(.8,-.6);

+angle=atand(imag(p)/real(p));

+dell=acosd(pf);

+Iq=Xq*cosd(dell+angle);

+Id=Xq*sind(dell+angle);

+function [A] = p2z(R,Theta)

+ if argn(2) <> 2 then

+    error("incorrect number of arguments.");

+ end

+ if ~and(size(R) == size(Theta)) then

+   error("arguments must be of the same dimension.");

+ end

+ A = R.*exp(%i*%pi*Theta/180.);

+endfunction

+Iqq=p2z(Iq,(20.3015));

+Idd=p2z(Id,(-69.685));

+

+EE=E+complex(0,Xd)*Idd+complex(0,Xq)*Iqq;

+mprintf("Excitation voltage and load angles are %.4f+%.4f pu and %.4f degree\n",real(EE),imag(EE),angle);

+step=2*%pi/100;

+delta=0:step:%pi;

+PP=(1.6673*E/Xd)*sin(delta);

+plot(delta,PP,'-k');

+set(gca(),"auto_clear","off")

+PP1=.5*(1/Xq-1/Xd)*sin(2*delta);

+plot(delta,PP1,'k.');

+PP2=PP+PP1;

+plot(delta,PP2,'k');

+xlabel('power angle in radians');

+ylabel('Power output in per unit');

+title('Plot of power generated vs power angle');

+set(gca(),"auto_clear","on")

diff --git a/3793/CH5/EX5.5/exp_5_5.sce b/3793/CH5/EX5.5/exp_5_5.sce
new file mode 100644
index 000000000..a57e5c163
--- /dev/null
+++ b/3793/CH5/EX5.5/exp_5_5.sce
@@ -0,0 +1,19 @@
+clear;

+clc;

+//case1

+V2=2300;

+S=150;

+V1=11500+V2;

+a=(V1-V2)/V2; //two winding transformer ratio

+

+aa=a+1; //autotransformer ratio

+output=((1+a)/a)*S;

+mprintf("output for 1 case %f KVA\n",output);

+//case 2

+V11=13.8;

+V22=11.5;

+a1=(V11-V22)/V22;

+

+output1=((1+a1)/a1)*S;

+mprintf("output for 2 case %f KVA\n",output1);

+

diff --git a/3793/CH5/EX5.6/exp_5_6.sce b/3793/CH5/EX5.6/exp_5_6.sce
new file mode 100644
index 000000000..78ff89fb0
--- /dev/null
+++ b/3793/CH5/EX5.6/exp_5_6.sce
@@ -0,0 +1,20 @@
+clear;

+clc;

+v1=66;

+v2=11;

+v3=6.6;

+s1=20;

+s2=10;

+s3=5;

+Xps=0.1;

+Xpt=0.12;

+Xst=.08;

+//now these rectance in pu and converted into 50 MVA base

+xps=Xps*(50/s1);

+xpt=Xpt*(50/s1);

+xst=Xst*(50/s2);

+Xp=complex(0,((xps+xpt-xst)/2));

+Xs=complex(0,((xps-xpt+xst)/2));

+Xs1=complex(0,((-xps+xpt+xst)/2));

+mprintf(" pu leakage reactances are %f, %f and %f",imag(Xp),imag(Xs),imag(Xs1));

+

diff --git a/3793/CH5/EX5.7/exp_5_7.sce b/3793/CH5/EX5.7/exp_5_7.sce
new file mode 100644
index 000000000..0be71feab
--- /dev/null
+++ b/3793/CH5/EX5.7/exp_5_7.sce
@@ -0,0 +1,17 @@
+clear;

+clc;

+v=220;

+s=5;

+z=4.5;

+Vb=11;

+sb=50;

+Zb=(Vb^2)/s;

+Zpu=z/Zb;

+mprintf("pu leakage reactance is %f\n",Zpu);

+a=Vb/v;

+Zs=z/(a^2);

+//case2

+vb1=220;

+Zb1=(vb1^2)/s;

+Zpu1=Zs/Zb1;

+mprintf("Ratio of pu leakage reactances are %f",Zpu1);

diff --git a/3793/CH5/EX5.8/exp_5_8.sce b/3793/CH5/EX5.8/exp_5_8.sce
new file mode 100644
index 000000000..194e3d102
--- /dev/null
+++ b/3793/CH5/EX5.8/exp_5_8.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+s=5;

+v1=11;

+v2=66;

+X1=.08;

+xm=75;

+Z1=(v1^2)/s;

+X11=X1*Z1;

+Xmm=xm*Z1;

+mprintf(" Actual reactance for primary X1=%f ohm and Xm=%f ohm\n",X11,Xmm);

+//case2

+Z2=(v2^2)/s;

+X2=X1*Z2;

+X2m=xm*Z2;

+mprintf(" Actual reactance for secondary X1=%f ohm and Xm=%f ohm\n",X2,X2m);

+mprintf("The pu values are independent of the side to which they are refeered. Therefore the pu values of X1 and Xm remain unchanged for all types of 3 phase transformer connections.")

+

diff --git a/3793/CH5/EX5.9/exp_5_9.sce b/3793/CH5/EX5.9/exp_5_9.sce
new file mode 100644
index 000000000..6d6734f8d
--- /dev/null
+++ b/3793/CH5/EX5.9/exp_5_9.sce
@@ -0,0 +1,33 @@
+clear;

+clc;

+s=50;

+vt=150;

+Sg1=50;

+vg1=11;

+Xg1=0.1;

+Sg2=40;

+vg2=6.6;

+xg2=.12;

+St1=100;

+Xt1=.15;

+St2=50;

+xt2=.1;

+vt1=220;

+l=75;

+pf=0.8;

+Z34=complex(30,150);

+Z35=complex(20,40);

+Z45=complex(25,60);

+Zb=(vt1^2)/100;

+z34=Z34/Zb;

+z35=Z35/Zb;

+z45=Z45/Zb;

+mprintf("reactances in pu are Z34=%f+j%f pu, Z35=%f+j%f pu and Z45=%f+j%f pu\n",real(z34),imag(z34),real(z35),imag(z35),real(z45),imag(z45));

+vbg1=11;

+vbg2=6.313;

+Xg11=Xg1*(vt/Sg1)*((vbg1/vbg2)^2);

+Xt11=Xt1*(vt/St1)*((vbg1/vbg2)^2);

+Xg22=xg2*(vt/Sg2)*((vg2/vbg2)^2);

+Xt22=xt2*(vt/St2)*((vg2/vbg2)^2);

+l1=(l/vt)*complex(.8,.6);

+mprintf("reactance of transformer and generator in pu are Xg1=%f pu, Xg2=%f pu, Xt1=%f pu, Xt2=%f pu and load=%f+j%f pu ",Xg11,Xg22,Xt11,Xt22,real(l1),imag(l1));