10 files changed, 283 insertions, 0 deletions
diff --git a/1076/CH19/EX19.1/19_1.sce b/1076/CH19/EX19.1/19_1.sce
new file mode 100755
index 000000000..5caa4a327
--- /dev/null
+++ b/1076/CH19/EX19.1/19_1.sce
@@ -0,0 +1,38 @@
+clear;
+clc;
+
+V=400e3;
+X=96;
+Y=.001 * exp(%i * 90/180*%pi);
+Sb=500e6;
+Pr=2;
+
+Z=complex(0,X);
+Zb=V*V/Sb;
+A=1+(Z*Y/2);
+B=Z/Zb;
+
+
+Vs=1;
+Vr=.95;
+d1=asin(Pr*abs(B)/(Vs*Vr));
+Qr1=((Vs*Vr)/abs(B)) * cos(d1) - (abs(A)* Vr*Vr/abs(B))
+Qs1= (abs(A)* Vs*Vs/abs(B)) - ((Vs*Vr)/abs(B)) * cos(d1)
+
+mprintf("\n(a) Qr= %.3f MVAR, Qs=%.3f MVAR", Qr1, Qs1);
+
+Vs=1;
+Vr=.9;
+d2=asin(Pr*abs(B)/(Vs*Vr));
+Qr2=((Vs*Vr)/abs(B)) * cos(d2) - (abs(A)* Vr*Vr/abs(B))
+Qs2= (abs(A)* Vs*Vs/abs(B)) - ((Vs*Vr)/abs(B)) * cos(d2)
+
+mprintf("\n(b) Qr= %.3f MVAR, Qs=%.3f MVAR", Qr2, Qs2);
+
+Vs=1;
+Vr=1;
+d3=asin(Pr*abs(B)/(Vs*Vr));
+Qr3=((Vs*Vr)/abs(B)) * cos(d3) - (abs(A)* Vr*Vr/abs(B))
+Qs3= (abs(A)* Vs*Vs/abs(B)) - ((Vs*Vr)/abs(B)) * cos(d3)
+
+mprintf("\n(c) Qr= %.3f MVAR, Qs=%.3f MVAR", fix(Qr3*1000)/1000, fix(Qs3*1000)/1000);
diff --git a/1076/CH19/EX19.10/19_10.sce b/1076/CH19/EX19.10/19_10.sce
new file mode 100755
index 000000000..e97b2b153
--- /dev/null
+++ b/1076/CH19/EX19.10/19_10.sce
@@ -0,0 +1,42 @@
+clear;
+clc;
+
+d=500;
+z=complex(.105, .3768);
+y=complex(0, 2.822e-6);
+Z=z*d;
+Y=y*d;
+YZ=Y*Z;
+
+A=1+(YZ/2)+((YZ)^2/24);
+B=Z * (1+(YZ/6)+((YZ)^2/120));
+C=Y * (1+(YZ/6)+((YZ)^2/120));
+D=A;
+
+Ys=-.6*Y;
+
+As=1;
+Ds=1;
+Bs=0;
+Cs=Ys;
+
+Anew=A + (B*Cs);
+
+Vr=round(220e3/sqrt(3));
+Pr=40e6;
+pf=.9;
+pfa=-1*acos(pf);
+Irm=Pr/(3*Vr);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+Vs=(A*Vr)+(B*Ir);
+Vr0=abs(Vs)/(round(abs(Anew)*100)/100);
+
+VRc= (Vr0-Vr)*100/Vr
+mprintf("\n(a) Voltage Regulation = %.2f percent",VRc);
+
+Vr02=abs(Vs)/abs(A);
+
+VRc2= (Vr02-Vr)*100/Vr
+mprintf("\n(b) Voltage Regulation(uncompensated) = %.2f percent",VRc2);
+
diff --git a/1076/CH19/EX19.11/19_11.sce b/1076/CH19/EX19.11/19_11.sce
new file mode 100755
index 000000000..dcd88574b
--- /dev/null
+++ b/1076/CH19/EX19.11/19_11.sce
@@ -0,0 +1,37 @@
+clear;
+clc;
+
+v1=220;
+v2=132
+vb1=220;
+n=132/220
+vb2=vb1*n
+
+Sb=200;
+Zb= vb2 *vb2/Sb;
+
+x1=75;
+x2=70;
+x3=90
+
+X1=x1/Zb;
+X2=x2/Zb;
+X3=x3/Zb;
+X1=fix(X1*100)/100
+X2=fix(X2*100)/100
+X3=fix(X3*100)/100
+
+
+
+Xt1=.08
+Xt2=.08
+
+X=((Xt1+X1)^-1 +(Xt2+X2)^-1 +(X3)^-1 )^-1;
+X=fix(X*1000)/1000
+
+dV=4
+dS= Sb/X;
+dQ=round(dS*1000/vb2)*1e-3 * dV
+
+mprintf("\n(a)\n X1= %.2f pu; X2= %.2f pu; X3= %.2f pu", X1, X2, X3)
+mprintf("\n(b)\n VAR injected = %.2f MVAR",dQ)
diff --git a/1076/CH19/EX19.2/19_2.sce b/1076/CH19/EX19.2/19_2.sce
new file mode 100755
index 000000000..40be50051
--- /dev/null
+++ b/1076/CH19/EX19.2/19_2.sce
@@ -0,0 +1,18 @@
+clear;
+clc;
+
+V=10000;
+P=12.5e6;
+f=50;
+Xl=4;
+
+
+d=asin(P*Xl/(V*V));
+
+VL=2*V*sin(d/2);
+QL=VL^2/4;
+Qc=QL/2;
+
+C=Qc/(2*%pi*f*V*V);
+
+mprintf("C=%.1f e-6F",C*1e6);
diff --git a/1076/CH19/EX19.3/19_3.sce b/1076/CH19/EX19.3/19_3.sce
new file mode 100755
index 000000000..36da75e83
--- /dev/null
+++ b/1076/CH19/EX19.3/19_3.sce
@@ -0,0 +1,20 @@
+clear;
+clc;
+
+xs=.22;
+xl=.15;
+Sb=1000;
+Vr=1;
+
+X=xl+xs;
+
+Pr=1;
+
+pf=.8;
+pfa=acos(pf);
+Qr=Pr*tan(pfa);
+
+Vs=complex(Vr + (X * Qr / Vr) , (X * Pr / Vr));
+V=abs(Vs);
+
+mprintf("Vr = %.2f ang (%.1f) deg pu",V, fix(atand(imag(Vs)/real(Vs))*10)/10)
diff --git a/1076/CH19/EX19.5/19_5.sce b/1076/CH19/EX19.5/19_5.sce
new file mode 100755
index 000000000..a0da9823c
--- /dev/null
+++ b/1076/CH19/EX19.5/19_5.sce
@@ -0,0 +1,9 @@
+clear;
+clc;
+
+vf=3;
+VA=8000;
+
+dQ=vf/100 * VA
+
+mprintf("capacity of SVS= +- %d MVAR", dQ);
diff --git a/1076/CH19/EX19.6/19_6.sce b/1076/CH19/EX19.6/19_6.sce
new file mode 100755
index 000000000..e8d98364c
--- /dev/null
+++ b/1076/CH19/EX19.6/19_6.sce
@@ -0,0 +1,26 @@
+clear;
+clc;
+
+V=220e3;
+Zl=complex(.8, .2);
+
+Xline=.2;
+Xt=.05;
+Sb=100e6;
+Vb=220e3;
+v=V/Vb;
+
+X=Xline+ Xt;
+I=conj(Zl/v)
+phi1=atand(imag(I)/real(I))
+
+Vbus=1+ I * X*exp(%i * %pi/2)
+
+phi2=atand(imag(Vbus)/real(Vbus))
+
+vbus=abs(Vbus)
+vbus=round(vbus *1000)/1000
+vbus=vbus*Vb*1e-3;
+pf=cosd(-phi1+phi2)
+
+mprintf("Voltage at bus = %.2f Kv, pf= %.3f lagging", vbus, pf)
diff --git a/1076/CH19/EX19.7/19_7.sce b/1076/CH19/EX19.7/19_7.sce
new file mode 100755
index 000000000..cdac2ad30
--- /dev/null
+++ b/1076/CH19/EX19.7/19_7.sce
@@ -0,0 +1,26 @@
+clear;
+clc;
+//the data used is from Ex 19.6, not 19.5 as incorrectly mentioned in statement
+V=220e3;
+Zl=complex(.8,0);
+
+Xline=.2;
+Xt=.05;
+Sb=100e6;
+Vb=220e3;
+v=V/Vb;
+
+X=Xline+ Xt;
+I=conj(Zl/v)
+phi1=atand(imag(I)/real(I))
+
+Vbus=1+( I * X*exp(%i * %pi/2))
+
+phi2=atand(imag(Vbus)/real(Vbus))
+
+vbus=abs(Vbus)*Vb*1e-3;
+pf=cosd(-phi1+phi2)
+
+
+mprintf("Voltage at bus = %.1f Kv, pf= %.2f", vbus, pf)
+disp("the data used is from Ex 19.6, not 19.5 as incorrectly mentioned in statement")
diff --git a/1076/CH19/EX19.8/19_8.sce b/1076/CH19/EX19.8/19_8.sce
new file mode 100755
index 000000000..595d82490
--- /dev/null
+++ b/1076/CH19/EX19.8/19_8.sce
@@ -0,0 +1,36 @@
+clear;
+clc;
+
+Z=complex(180*cosd(75) , 180*sind(75));
+Y=complex(1e-3*cosd(90) , 1e-3*sind(90));
+
+g=sqrt(Y*Z);
+Zc=sqrt(Z/Y);
+Zc=round(abs(Zc)) * exp (%i * atan(imag(Zc)/real(Zc)))
+
+Z1=Zc * sinh(g);
+Y1=(1/Zc) *( (cosh(g)-1)/sinh(g));
+A=cosh(g);
+B=round(abs(sinh(g) )*1000) * exp (%i * atan(imag(sinh(g) )/real(sinh(g) ))) * Zc/1000
+C=B / (Zc *Zc)
+D=A;
+
+
+mprintf("\n(a)");
+
+mprintf("\nA =D= %.4f ang (%.2f) deg",abs(A), round(atand(imag(A)/real(A))*100)/100)
+mprintf("\nB = %.2f ang (%.2f) deg ohm",abs(B), round(atand(imag(B)/real(B))*100)/100)
+mprintf("\nC = %.5f ang (%.2f) deg ohm",abs(C), round(atand(imag(C)/real(C))*100)/100 +180)
+
+
+B2=B-(0.6 * imag(Z) *%i);
+Y2=(cosh(g)-1)/(sinh(g)*Zc);
+
+A2=1+(B2*Y2)
+D2=A2;
+C2=(2* Y2 )+ (B2*Y2*Y2);
+mprintf("\n(b)");
+
+mprintf("\nA =D= %.4f ang (%.2f) deg",abs(A2), round(atand(imag(A2)/real(A2))*100)/100)
+mprintf("\nB = %.0f ang (%.2f) deg ohm",abs(B2), round(atand(imag(B2)/real(B2))*100)/100)
+mprintf("\nC = %.3f ang (%.1f) deg ohm",abs(C2), round(atand(imag(C2)/real(C2))*100)/100 +180)
diff --git a/1076/CH19/EX19.9/19_9.sce b/1076/CH19/EX19.9/19_9.sce
new file mode 100755
index 000000000..b3d3e4ce9
--- /dev/null
+++ b/1076/CH19/EX19.9/19_9.sce
@@ -0,0 +1,31 @@
+clear;
+clc;
+
+Z=complex(180*cosd(75) , 180*sind(75));
+Y=complex(1e-3*cosd(90) , 1e-3*sind(90));
+YZ=Z*Y;
+
+
+B=Z;
+A=1+(YZ/2)
+D=A
+C= Y* (1+(YZ/4))
+
+mprintf("\n(a)");
+mprintf("\nA =D= %.3f ang (%.2f) deg",abs(A), round(atand(imag(A)/real(A))*100)/100)
+mprintf("\nB = %.0f ang (%.0f) deg ohm",abs(B), round(atand(imag(B)/real(B))*100)/100)
+mprintf("\nC = %.6f ang (%.1f) deg ohm",abs(C), round(atand(imag(C)/real(C))*100)/100 +180)
+
+
+B2=B-(0.6 * imag(Z) *%i);
+YZ2=B2*Y;
+A2=1+(YZ2/2)
+D2=A2
+C2= Y* (1+(YZ2/4))
+
+
+mprintf("\n(b)");
+
+mprintf("\nA =D= %.3f ang (%.2f) deg",fix(abs(A2)*1000)/1000, round(atand(imag(A2)/real(A2))*100)/100)
+mprintf("\nB = %.2f ang (%.2f) deg ohm",abs(B2), round(atand(imag(B2)/real(B2))*100)/100)
+mprintf("\nC = %.6f ang (%.1f) deg ohm",abs(C2), round(atand(imag(C2)/real(C2))*100)/100 +180)