6 files changed, 172 insertions, 0 deletions
diff --git a/3793/CH11/EX11.1/exp_11_1.sce b/3793/CH11/EX11.1/exp_11_1.sce
new file mode 100644
index 000000000..628114475
--- /dev/null
+++ b/3793/CH11/EX11.1/exp_11_1.sce
@@ -0,0 +1,25 @@
+clear;
+clc;
+f=50;
+s=100;
+H=3.5;
+p=.16;
+sb=500;
+ip=.18;
+Pole=4;
+K=H*s;
+mprintf("Kinetic energy stored is %.3f MJ\n",K);
+Pa=(ip-p)*sb;
+A=(Pa*360*f)/(2*H*s);
+mprintf("acceleration of generator is %.3f degree electrical per second sqr\n",A);
+a=7.5;
+accp=a/f;
+Ns=120*f/Pole;
+rotora=(2.996*accp)^2; //change in rotor angle equation obtained with the help of eq11.1.1 and integrating it
+
+mprintf("change in rotor angle is %.3f rad \n",rotora);
+del=.202;
+rv=2.996*sqrt(del);
+vel=(rv*120)/(3.14*Pole);
+NN=Ns+vel;
+mprintf("Speed is %.3f rpm",NN);
diff --git a/3793/CH11/EX11.2/exp_11_2.sce b/3793/CH11/EX11.2/exp_11_2.sce
new file mode 100644
index 000000000..5aa5422d2
--- /dev/null
+++ b/3793/CH11/EX11.2/exp_11_2.sce
@@ -0,0 +1,17 @@
+clear;
+clc;
+pole=2;
+s=50;
+v=11;
+pf=.8;
+H=6;
+Ns=3000;
+inpinc=62000;
+K=H*s;
+mprintf("Kinetic energy stored is %.3f MJ\n",K);
+po=s*pf;
+pi=inpinc*735.5/10^6;
+f=50;
+Ap=pi-po;
+A=(Ap*360*f)/(2*H*s);
+mprintf("acceleration of generator is %.3f degree electrical per second sqr\n",A);
diff --git a/3793/CH11/EX11.3/exp_11_3.sce b/3793/CH11/EX11.3/exp_11_3.sce
new file mode 100644
index 000000000..8c608bc74
--- /dev/null
+++ b/3793/CH11/EX11.3/exp_11_3.sce
@@ -0,0 +1,47 @@
+clear;
+clc;
+function [S,pin,deltam,itr,delta1]=pinstab (V,X1,Xt,Xd,pf,typee,tolr);
+    Xtot=X1+Xt+Xd;
+    Pu=input('Generator  output power');
+    phi=acosd(pf);
+    Qu=Pu*tand(phi);
+    if typee==0 
+        S=Pu+%i*Qu;
+    else
+        S=Pu-%i*Qu;
+    end
+    I=conj(S)/conj(V);
+    Edash=V+I*(%i*Xtot);
+    Edash=abs(Edash);
+    delta0=asin(Pu*Xtot/(Edash*V));
+    itr=0;
+    deltam=input('initial estimate of deltam');
+    ndeltam=0;
+    difff=abs(ndeltam-deltam);
+    while difff > tolr;
+        itr=itr+1;
+        fdeltam=cos(delta0)-(deltam-delta0)*sin(deltam)-cos(deltam);
+        dfdeltam=(deltam-delta0)*cos(deltam);
+        ndeltam=deltam+fdeltam/dfdeltam;
+        difff=abs(ndeltam-deltam);
+        deltam=ndeltam;
+        
+    end
+    delta1=%pi-deltam;
+    pin=(Edash*V/Xtot)*sin(delta1);
+    deltam=deltam*180/%pi;
+    delta1=delta1*180/%pi;
+    mprintf("Magnitude of power input without loosing synchronism is %.4f",pin);
+    
+    
+    
+    
+endfunction
+V=1;
+X1=.15;
+Xt=.2;
+Xd=.15;
+pf=.8;
+typee=0;
+tolr=.001;
+[S,pin,deltam,itr,delta1]=pinstab (V,X1,Xt,Xd,pf,typee,tolr);
diff --git a/3793/CH11/EX11.4/exp_11_4.sce b/3793/CH11/EX11.4/exp_11_4.sce
new file mode 100644
index 000000000..d3d458f24
--- /dev/null
+++ b/3793/CH11/EX11.4/exp_11_4.sce
@@ -0,0 +1,38 @@
+clear;
+clc;
+function [deltam]=stabnr (tolr);
+    itr=0;
+    ndeltam=0;
+    deltam=input('Initial estimate of deltam');
+    difff=abs(ndeltam-deltam);
+    while difff>tolr;
+        itr=itr+1;
+        fdeltam=2.7202*cos(deltam)+.8*deltam-2.8247;
+        dfdeltam=2.7202*sin(deltam);
+        ndeltam=deltam+fdeltam/dfdeltam;
+        difff=abs(ndeltam-deltam);
+        deltam=ndeltam;
+        
+    end
+    deltam=deltam*180/%pi;
+    mprintf("Maximum swing of the rotor angle is %.4f degree, since it is less than(pi-delta0) therefore system will remain stable\n",deltam);
+endfunction
+phi=-acosd(.8);
+S=complex(.8,.6);
+pu=.8;
+V=1;
+I=conj(S)/conj(V);
+Xtot=.5;
+E=V+%i*I*Xtot;
+E=abs(E);
+delta0=asin(pu*Xtot/(E*V));
+mprintf("Deta0=%.4f radian\n",delta0);
+tolr=.001;
+[deltam]=stabnr (tolr);
+deltac=acos(-.0866/(2.7202));
+H=6;
+M=H/(%pi*50);
+pi=.8;
+t=sqrt(2*M*(deltac-delta0)/pi);
+mprintf("Critical angle is %.4f radian and time is %.4f seconds",deltac,t);
+
diff --git a/3793/CH11/EX11.5/exp_11_5.sce b/3793/CH11/EX11.5/exp_11_5.sce
new file mode 100644
index 000000000..eeda6a7fc
--- /dev/null
+++ b/3793/CH11/EX11.5/exp_11_5.sce
@@ -0,0 +1,21 @@
+clear;
+clc;
+Xtf=.2+.2+(.3*.6/0.9);
+pi=0.9;
+po=pi;
+del1=asin(Xtf*pi/(1.2*1));
+
+Pm=1.2*1/Xtf;
+//fault condition
+Xtf1=(.4*.3+.3*.3+.3*.4)/.3;
+Pm1=1.2*1/Xtf1;
+//post fault condition
+Xtf2=.2+.2+.3;
+Pm2=1.2*1/Xtf2;
+delm=(%pi-(asin(pi/Pm2)));
+
+delc=acos((pi*(delm-del1)+Pm2*cos(delm)-Pm1*cos(del1))/(Pm2-Pm1));
+
+
+mprintf("rotor angle is %.3f radian \n",del1);
+mprintf("Critical clearing angle is %.3f radian",delc);
diff --git a/3793/CH11/EX11.9/exp_11_9.sce b/3793/CH11/EX11.9/exp_11_9.sce
new file mode 100644
index 000000000..86d29a0e0
--- /dev/null
+++ b/3793/CH11/EX11.9/exp_11_9.sce
@@ -0,0 +1,24 @@
+clear;
+clc;
+pf=.8;
+f=50;
+rp=.8;
+X=.4;
+Xd=.2;
+H=10;
+v=1;
+Xeq=Xd+X;
+Ig=rp/(v*pf);
+angle=acosd(pf);
+E=sqrt(((v+Ig*Xd*Xeq)^2)+((Ig*Xd*pf)^2));
+del=atand((Ig*Xd*pf)/(v+Ig*Xd*Xeq));
+P=(E*v)/Xeq;
+mprintf("Steady state power limit is %.3f pu\n",P);
+Pc=cosd(del)*P;
+mprintf("Synchronizing power coefficient is %.3f pu\n",Pc);
+M=H/(3.14*f);
+gaama=sqrt(Pc/M);
+fre=gaama/(2*%pi);
+mprintf("frequency of free oscillation  is %.3f Hz\n",fre);
+T=1/fre;
+mprintf("time period of free oscillation is %.3f s\n",T);