diff options
Diffstat (limited to '1076/CH3')
35 files changed, 1184 insertions, 0 deletions
diff --git a/1076/CH3/EX3.1/3_1.sce b/1076/CH3/EX3.1/3_1.sce new file mode 100755 index 000000000..9ac356a32 --- /dev/null +++ b/1076/CH3/EX3.1/3_1.sce @@ -0,0 +1,27 @@ +clear;
+clc;
+
+
+sg1=10e6;
+vg1=13.2e3;
+sg2=15e6;
+vg2=13.2e3;
+sm1=8e6;
+vm1=12.5e3;
+sm2=12e6;
+vm2=12.5e3;
+Xg=15;
+Xm=20;
+
+sb=50e6;
+vb=13.8e3;
+
+xg1=Xg * (vg1/vb)^2 * (sb/sg1);
+xg2=Xg * (vg2/vb)^2 * (sb/sg2);
+xm1=Xm * (vm1/vb)^2 * (sb/sm1);
+xm2=Xm * (vm2/vb)^2 * (sb/sm2);
+ mprintf ("\nReactance of Generator 1= %.2f percent",xg1);
+ mprintf ("\nReactance of Generator 2= %.2f percent",xg2);
+ mprintf ("\nReactance of Motor 1= %.2f percent",xm1);
+ mprintf ("\nReactance of Motor 2= %.2f percent",xm2);
+
diff --git a/1076/CH3/EX3.10/3_10.sce b/1076/CH3/EX3.10/3_10.sce new file mode 100755 index 000000000..c6e8f15d5 --- /dev/null +++ b/1076/CH3/EX3.10/3_10.sce @@ -0,0 +1,45 @@ +clear;
+clc;
+
+d=100;
+f=50;
+r=.153;
+l=1.21e-3;
+c=.00958e-6;
+xl=2*%pi*f*l;
+xc=1/(2*%pi*f*c);
+Vr=110e3/sqrt(3);
+Pr=20e6;
+pf=.9;
+pfa=-1*acos(pf);
+Irm=Pr/(3*Vr*pf);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+Z=complex(r,xl) * d;
+Y=complex(0,1/xc) *d;
+
+//disp(abs(Z),Y,abs(Ir));
+
+Vs=(Vr*(1+((Z*Y)/2)))+ (Ir*Z);
+V=abs(Vs);
+vs=sqrt(3)*V*1e-3;
+mprintf("\nSending End Voltage = %.2f kV",vs);
+
+Is=(Vr*Y*(1+(Z*Y/4)))+(Ir*(1+(Y*Z/2)));
+I=abs(Is);
+mprintf("\nSending End Current = %.2f A",I);
+
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+phi2=atan(imag(Is)/real(Is))*(180/%pi);
+phi=phi1-phi2;
+pfs=cosd(phi);
+mprintf("\nSending End Power factor = %.3f ",pfs);
+
+
+Vr0=V/(1+(Y*Z/2));
+V0=abs(Vr0);
+VR=(V0-Vr)/Vr;
+mprintf("\nVoltage Regulation = %.2f ",VR*100);
+
+eff=Pr*100/(3*pfs*V*I)
+mprintf("\nEfficiency = %.0f percent ",eff);
diff --git a/1076/CH3/EX3.11/3_11.sce b/1076/CH3/EX3.11/3_11.sce new file mode 100755 index 000000000..11a33480a --- /dev/null +++ b/1076/CH3/EX3.11/3_11.sce @@ -0,0 +1,49 @@ +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;
+
+A=round(abs(A)*10000)/10000 * exp(%i * round(atan(imag(A)/real(A))*100)/100)
+B=round(abs(B)*1000)/1000 * exp(%i * round(atan(imag(B)/real(B))*100)/100)
+
+//see(B)
+
+
+Vr=220e3/sqrt(3); //incorrectly taken as 127021 in textbook.
+Vr=round(Vr)
+Pr=40e6;
+pf=.9;
+pfa=-1*acos(pf);
+Irm=Pr/(3*Vr);
+Ir=complex(Irm *pf, Irm * round(sin(pfa)*100)/100);
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs);
+vs=sqrt(3)*V*1e-3;
+phi1=atand(imag(Vs)/real(Vs));
+mprintf("\nSending End Voltage = %.2f kV",vs);
+
+Is=(C*Vr)+(D*Ir);
+I=abs(Is);
+phi2=atand(imag(Is)/real(Is))
+mprintf("\nSending End Current = %.2f A",I);
+
+phi=phi2-phi1;
+pfs=cosd(phi);
+mprintf("\nSending End Power factor = %.3f ",pfs);
+mprintf("\nSending End Power Angle = %.3f ",phi1);
+
+MVA=sqrt(3) * vs* I /1000;
+mprintf("\nSending End Power = %.3f ",MVA);
+
+disp("difference in results is due to taking Vr= 127021V instead of 127017V")
diff --git a/1076/CH3/EX3.12/3_12.sce b/1076/CH3/EX3.12/3_12.sce new file mode 100755 index 000000000..b19fee92b --- /dev/null +++ b/1076/CH3/EX3.12/3_12.sce @@ -0,0 +1,36 @@ +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;
+
+Vr=220e3/sqrt(3);
+Pr=40e6;
+pf=.9;
+pfa=-1*acos(pf);
+Ir=0;
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs);
+vs=sqrt(3)*V*1e-3;
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+mprintf("\nSending End Voltage = %.2f kV",vs);
+
+Is=(C*Vr)+(D*Ir);
+I=abs(Is);
+phi2=atan(imag(Is)/real(Is))*(180/%pi);
+mprintf("\nSending End Current = %.1f A",I);
+
+phi2=phi2+180;
+phi=phi1-phi2;
+pfs=cosd(phi);
+mprintf("\nSending End Power factor = %.4f ",pfs);
diff --git a/1076/CH3/EX3.13/3_13.sce b/1076/CH3/EX3.13/3_13.sce new file mode 100755 index 000000000..0aa12e3b8 --- /dev/null +++ b/1076/CH3/EX3.13/3_13.sce @@ -0,0 +1,18 @@ +clear;
+clc;
+
+Z=complex(14.1, 51.48);
+Y=complex(0, 1.194e-3);
+Zc=sqrt(Z/Y);
+g=sqrt(Z*Y);
+A=cosh(g);
+B=Zc* sinh(g);
+C=sinh(g)/Zc;
+D=A;
+
+mprintf("\nZc= %s", string(round(abs(Zc)*1000)/1000) +'/_'+ string(round(atand(imag(Zc)/real(Zc))*100)/100) )
+mprintf("\npropgation const= %s", string(round(abs(g)*1000)/1000) +'/_'+ string(round(atand(imag(g)/real(g))*100)/100) )
+mprintf("\nA= %s", string(round(abs(A)*1000)/1000) +'/_'+ string(round(atand(imag(A)/real(A))*100)/100) )
+mprintf("\nB= %s", string(round(abs(B)*1000)/1000) +'/_'+ string(round(atand(imag(B)/real(B))*100)/100) )
+mprintf("\nC= %s", string(round(abs(C)*1000)/1000) +'/_'+ string(round(atand(imag(C)/real(C))*100)/100 +180) )
+mprintf("\nD= %s", string(round(abs(D)*1000)/1000) +'/_'+ string(round(atand(imag(D)/real(D))*100)/100) )
diff --git a/1076/CH3/EX3.14/3_14.sce b/1076/CH3/EX3.14/3_14.sce new file mode 100755 index 000000000..a94aa6f5d --- /dev/null +++ b/1076/CH3/EX3.14/3_14.sce @@ -0,0 +1,29 @@ +clear;
+clc;
+
+Z=complex(200*cosd(80), 200*sind(80));
+Y=complex(.0013*cosd(90),.0013*sind(90));
+YZ=Y*Z;
+
+A=1+(YZ/2)+((YZ)^2/24);
+B=Z * (1+(YZ/6)+((YZ)^2/120));
+phiA=atan(imag(A)/real(A))*(180/%pi);
+phiB=atan(imag(B)/real(B))*(180/%pi);
+
+P=60e6
+pf=.8;
+Vs=round(220/sqrt(3))*1e3;
+VIrm=P/(3*pf)
+pfa=acos(pf);
+VIr=complex(VIrm *pf, VIrm * sin(pfa));
+
+pfa=pfa*(180/%pi);
+a=(round(abs(A)*1000)/1000)^2;
+b=round(((2*(abs(A)* (VIrm)*cosd(phiA)*abs(B)*cosd(phiB-pfa)) )+ (2*(abs(A)* (VIrm)*sind(phiA)*abs(B)*sind(phiB-pfa))) - (Vs^2))/1e7)*1e7;
+c=abs(B)^2 * (VIrm)^2;
+Vr=sqrt((-b+sqrt((b*b)-(4*a*c)))/(2*a));
+
+vr=sqrt(3) * Vr / 1000;
+Ir=VIr/(Vr*pf);
+mprintf("Receiveing End Line voltage= %.0f kV", fix(vr));
+mprintf("\n Receiveing End Line Current= %.0f A", Ir);
diff --git a/1076/CH3/EX3.15/3_15.sce b/1076/CH3/EX3.15/3_15.sce new file mode 100755 index 000000000..c7495c1a4 --- /dev/null +++ b/1076/CH3/EX3.15/3_15.sce @@ -0,0 +1,16 @@ +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);
+
+Z1=Zc * sinh(g);
+Y1=(1/Zc) *1e3*( (cosh(g)-1)/sinh(g));
+
+
+mprintf("\nZpi= %s", string(round(abs(Z1)*100)/100) +'/_'+ string(round(atand(imag(Z1)/real(Z1))*100)/100) )
+mprintf("\nYpi/2= %s *1E-3", string(round(abs(Y1)*10000)/10000) +'/_'+ string(round(atand(imag(Y1)/real(Y1))*10)/10) )
+
diff --git a/1076/CH3/EX3.16/3_16.sce b/1076/CH3/EX3.16/3_16.sce new file mode 100755 index 000000000..0bf417b48 --- /dev/null +++ b/1076/CH3/EX3.16/3_16.sce @@ -0,0 +1,63 @@ +clear;
+clc;
+
+Vr=132e3/sqrt(3);
+P=40e6;
+pf=.8;
+Irm=P/(3*Vr)
+pfa=-1* acos(pf);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+Z=complex(52, 200)
+Y=complex(0, 1.5e-3)
+YZ=Y*Z;
+
+A=1+ (YZ/2);
+D=A;
+B=Z;
+C=Y*(1+(YZ/4));
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs)
+vs=V*sqrt(3)*1e-3;
+mprintf("(a)\nSending End Voltage= %.0f kV", vs)
+
+Is=(C*Vr)+(D*Ir);
+I=abs(Is)
+mprintf("\nSending End Current= %.1f A", I)
+
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+phi2=atan(imag(Is)/real(Is))*(180/%pi);
+phi=phi1-phi2;
+pfs=cosd(phi);
+mprintf("\nSending End pf= %.3f ",pfs)
+
+Ps=sqrt(3)* vs * I * pfs /1000;
+mprintf("\nSending End Power= %.1f ",Ps)
+
+
+//(b)
+Zc=sqrt(Z/Y);
+g=sqrt(Z*Y);
+A=cosh(g);
+B=Zc* sinh(g);
+C=sinh(g)/Zc;
+D=A;
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs)
+vs=V*sqrt(3)*1e-3;
+mprintf("\n\n\n(b)\nSending End Voltage= %.1f kV", vs)
+
+Is=(C*Vr)+(D*Ir);
+I=abs(Is)
+mprintf("\nSending End Current= %.1f A", fix(I*10)/10)
+
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+phi2=atan(imag(Is)/real(Is))*(180/%pi);
+phi=phi1-phi2;
+pfs=cosd(phi);
+mprintf("\nSending End pf= %.3f ",pfs)
+
+Ps=sqrt(3)* vs * I * pfs /1000;
+mprintf("\nSending End Power= %.1f ",Ps)
diff --git a/1076/CH3/EX3.17/3_17.sce b/1076/CH3/EX3.17/3_17.sce new file mode 100755 index 000000000..ec32fdce1 --- /dev/null +++ b/1076/CH3/EX3.17/3_17.sce @@ -0,0 +1,13 @@ +clear;
+clc;
+
+Y1=500^-1;
+Y2=1000^-1;
+Z=100;
+
+A= 1+Y2 * Z;
+B=Z;
+C=Y1+Y2+(Y1*Y2*Z);
+D=1+Y1 * Z
+
+mprintf("A= %.1f ; B= %.1f ohm ; C=%.1f *1e-3seimens; D= %.1f", A, B, C*1e3, D);
diff --git a/1076/CH3/EX3.18/3_18.sce b/1076/CH3/EX3.18/3_18.sce new file mode 100755 index 000000000..9d3a2973c --- /dev/null +++ b/1076/CH3/EX3.18/3_18.sce @@ -0,0 +1,50 @@ +clear;
+clc;
+
+A1=complex(.98 * cosd(2), .98* sind(2));
+B1=complex(28 * cosd(69), 28* sind(69));
+C1=complex(.0002 * cosd(80), .0002* sind(80));
+D1=A1;
+
+A2=complex(.95 * cosd(3), .95* sind(3));
+B2=complex(40 * cosd(85), 40* sind(85));
+C2=complex(.0004 * cosd(90), .0004* sind(90));
+D2=A2;
+
+//(a)
+
+A= (A1* A2) + (B1* C2);
+B= (A1* B2) + (B1* D2);
+C= (C1* A2) + (D1* C2);
+D= (C1* B2) + (D1* D2);
+
+
+mprintf("(a)");
+mprintf("\nA= %s", string(round(abs(A)*1000)/1000) +'/_'+ string(round(atand(imag(A)/real(A))*10)/10) )
+mprintf("\nB= %s", string(round(abs(B)*100)/100) +'/_'+ string(round(atand(imag(B)/real(B))*100)/100) )
+mprintf("\nC= %s *1e-4", string(round(abs(C)*100000)/10) +'/_'+ string(round(atand(imag(C)/real(C))*10)/10 ) )
+mprintf("\nD= %s", string(round(abs(D)*1000)/1000) +'/_'+ string(round(atand(imag(D)/real(D))*10)/10) )
+
+
+//(b)
+
+Vr=110e3/sqrt(3);
+pf=.95;
+Irm=200
+pfa=-1* acos(pf);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs)
+vs=V*sqrt(3)*1e-3;
+mprintf("\n\n\n(b)\nSending End Voltage= %.2f kV", vs)
+
+Is=(C*Vr)+(D*Ir);
+I=abs(Is)
+mprintf("\nSending End Current= %.1f A", I)
+
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+phi2=atan(imag(Is)/real(Is))*(180/%pi);
+phi=phi1-phi2;
+pfs=cosd(phi);
+mprintf("\nSending End pf= %.2f ",pfs)
diff --git a/1076/CH3/EX3.19/3_19.sce b/1076/CH3/EX3.19/3_19.sce new file mode 100755 index 000000000..6ebf93de9 --- /dev/null +++ b/1076/CH3/EX3.19/3_19.sce @@ -0,0 +1,11 @@ +clear;
+clc;
+
+Vr=220e3/sqrt(3);
+d=300;
+f=50;
+
+V=-1*((Vr*(2*%pi*f)^2 * d*d *1e-10)/18);
+Vs=Vr+V;
+vs=sqrt(3) * Vs /1000;
+mprintf("Sending end voltage=%.2f kV, and voltage rise =%.0f V/phase", vs, -1*V)
diff --git a/1076/CH3/EX3.2/3_2.sce b/1076/CH3/EX3.2/3_2.sce new file mode 100755 index 000000000..351e90f80 --- /dev/null +++ b/1076/CH3/EX3.2/3_2.sce @@ -0,0 +1,50 @@ +clear;
+clc;
+
+vg=11e3;
+sg=90e6;
+xg=.25;
+
+st1=100e6;
+vt1a=10e3;
+vt1b=132e3;
+nt1=vt1a/vt1b;
+xt1=.06;
+
+st2=30e6 *3;
+vt2a=66e3 * sqrt(3);
+vt2b=10e3;
+nt2=vt2a/vt2b;
+xt2=.05;
+
+sm1=50e6;
+vm1=10e3;
+xm1=.2;
+
+sm2=40e6;
+vm2=10e3;
+xm2=.2;
+
+xl=100;
+
+Sb=sg;
+Vbg=vg;
+
+Xg=xg * (vg/Vbg)^2 * (Sb/sg);
+Xt1=xt1 * (vt1a/Vbg)^2 * (Sb/st1);
+
+Vbl=Vbg/nt1;
+Xl=xl * (Sb) / (Vbl)^2 ;
+
+Vbm=Vbl/nt2;
+
+Xt2=xt2 * (vt2b/Vbm)^2 * (Sb/st2);
+Xm1=xm1 * (vm1/Vbm)^2 * (Sb/sm1);
+Xm2=xm2 * (vm2/Vbm)^2 * (Sb/sm2);
+
+ mprintf ("\nReactance of Generator = %.2f ",Xg);
+ mprintf ("\nReactance of Transformer 1= %.4f ",Xt1);
+ mprintf ("\nReactance of Line = %.3f ",Xl);
+ mprintf ("\nReactance of Transformer 2= %.3f ",Xt2);
+ mprintf ("\nReactance of Motor 1= %.3f ",Xm1);
+ mprintf ("\nReactance of Motor 2= %.3f ",Xm2);
diff --git a/1076/CH3/EX3.20/3_20.sce b/1076/CH3/EX3.20/3_20.sce new file mode 100755 index 000000000..83b93c3c1 --- /dev/null +++ b/1076/CH3/EX3.20/3_20.sce @@ -0,0 +1,21 @@ +clear;
+clc;
+
+E=11e3/sqrt(3);
+E1=complex(E*cosd(0), E*sind(0));
+E2=complex(E*cosd(40), E*sind(40));
+Z=complex(0,15);
+I=(E1-E2)/Z;
+
+if (real(I)<0) then
+ mprintf("E1=generator, E2=motor");
+ else
+ mprintf("E2=generator, E1=motor");
+end
+
+S1=3 * E1 *1e-6* conj(I);
+S2=3 * E2 *1e-6* conj(I);
+
+mprintf("\n\n(a)Real Power consumed by E2= %.3f MW,\nPower delivered by E1= %.3f MW", -1*real(S2), -1*real(S1) );
+mprintf("\n\n(b) Reactive Power supplied by E1= %.3f MVAr,\nPower supplied by E2= %.3f MVAr", imag(S1),imag(S2)*-1 );
+mprintf("\n\n(C) Reactive Power absorbed by line= %.3f MVAr", 2* imag(S1) );
diff --git a/1076/CH3/EX3.21/3_21.sce b/1076/CH3/EX3.21/3_21.sce new file mode 100755 index 000000000..e71659c2f --- /dev/null +++ b/1076/CH3/EX3.21/3_21.sce @@ -0,0 +1,49 @@ +clear;
+clc;
+
+Vr=132e3/sqrt(3);
+P=50e6;
+pf=.8;
+Irm=P/(3*Vr)
+pfa=-1* acos(pf);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+A=complex(.98*cosd(3), .98*sind(3));
+B=complex(110*cosd(75), 110*sind(75));
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs)
+vs=V*sqrt(3)*1e-3;
+mprintf("\n\n\n(a)\nSending End Voltage= %.1f kV", vs)
+
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+mprintf("\nPower Angle= %.2f ",phi1)
+
+
+Ss=(((vs)^2 * conj(A/B))-((sqrt(3)*Vr/1000)*(sqrt(3)*Vs/1000)/conj(B)));
+Ps=real(Ss);
+Qs=imag(Ss);
+
+mprintf("\n\n\n(b)\nSending End Active Power= %.1f MW", Ps)
+mprintf("\nSending End Reactive Power= %.1f MVar lagging", Qs)
+
+Pl=Ps-(P * cos(pfa)*1e-6);
+Ql=Qs-(P* -1*sin(pfa)*1e-6);
+mprintf("\n\n\n(c)\nLine Loss= %.1f MW", Pl)
+mprintf("\nMVar absorbed by line= %.1f MVar", Ql)
+
+Pr=(P * cos(pfa)*1e-6);
+Qr=(P * sin(pfa)*1e-6);
+Vs1=140;
+Vr1=132;
+
+bd=acos(( Pr+real(((Vr1)^2 * conj(A/B))) ) * (abs(B)/(Vs1 *Vr1)));
+Qr1= (((Vs1 *Vr1)/abs(B))*sin(bd))- imag(((Vr1)^2 * conj(A/B)));
+Q=-Qr-Qr1;
+mprintf("\n\n\n(d)\nCapacity of static capacitor= %.1f MVar", Q)
+
+Vs2=132;
+Vr2=132;
+bd2=asin(( imag(((Vr2)^2 * conj(A/B))) ) * (abs(B)/(Vs2 *Vr2)));
+P2= (((Vs2 *Vr2)/abs(B))*cos(bd2))- real(((Vr2)^2 * conj(A/B)));
+mprintf("\n\n\n(e)\nPower supplied= %.2f MW",P2)
diff --git a/1076/CH3/EX3.22/3_22.sce b/1076/CH3/EX3.22/3_22.sce new file mode 100755 index 000000000..56fa15ab4 --- /dev/null +++ b/1076/CH3/EX3.22/3_22.sce @@ -0,0 +1,18 @@ +clear;
+clc;
+
+Pr=90e6
+pf=.9;
+S=Pr/(3*pf);
+
+P=Pr/3
+Q=sqrt(S^2 - P^2);
+
+V1=220e3/sqrt(3);
+V2=220e3/sqrt(3);
+
+R=15
+X=50;
+
+tr=sqrt(1-(((R*P)+(X*Q))/(V1^2)));
+mprintf("tap setting tr= %.4f, ts=%.3f",fix(1e4/tr)/1e4,tr)
diff --git a/1076/CH3/EX3.23/3_23.sce b/1076/CH3/EX3.23/3_23.sce new file mode 100755 index 000000000..b6da1da6c --- /dev/null +++ b/1076/CH3/EX3.23/3_23.sce @@ -0,0 +1,32 @@ +clear;
+clc;
+
+Vb=132
+Sb=100
+X=.15
+
+v1=125
+V1=v1/Vb;
+Q1=50;
+Qpu1=Q1/Sb;
+
+Vn1=(V1 + sqrt(V1^2 - (4*Qpu1*X)))/(2*1);
+vn1=Vn1 * Vb;
+Vo1=33;
+t1=vn1/Vo1;
+
+v2=140
+V2=v2/Vb;
+Q2=20;
+Qpu2=Q2/Sb;
+
+Vn2=(V2 + sqrt(V2^2 - (4*Qpu2*X)))/(2*1);
+vn2=Vn2 * Vb;
+Vo2=33;
+t2=vn2/Vo2;
+
+tm=(t1+t2)/2;
+dt=tm-t1;
+ts=dt*100/tm
+
+mprintf("tap setting = +- %.0f percent", ts);
diff --git a/1076/CH3/EX3.26/3_26.sce b/1076/CH3/EX3.26/3_26.sce new file mode 100755 index 000000000..72c176dbc --- /dev/null +++ b/1076/CH3/EX3.26/3_26.sce @@ -0,0 +1,28 @@ +clear;
+clc;
+
+Vr1=132
+Vs1=140
+VA=40;
+pf=.8;
+Pr=VA*pf;
+pfa=-1* acos(pf);
+Qr=(VA * sin(pfa));
+
+A=complex(.98*cosd(3), .98*sind(3));
+B=complex(110*cosd(75), 110*sind(75));
+
+bd1=acos(( Pr+real(((Vr1)^2 * conj(A/B))) ) * (abs(B)/(Vs1 *Vr1)));
+Qr1= (((Vs1 *Vr1)/abs(B))*sin(bd1))-imag(((Vr1)^2 * conj(A/B)));
+
+Q1=-Qr-Qr1;
+mprintf("\n\n\n(a)\nCapacity of static capacitor= %.2f MVar leading", fix(Q1*100)/100)
+
+Prn=0;
+Qrn=0;
+
+bd2=acos((Prn+real(((Vr1)^2 * conj(A/B))) ) * (abs(B)/(Vs1 *Vr1)));
+Qr2= (((Vs1 *Vr1)/abs(B))*sin(bd2))-imag(((Vr1)^2 * conj(A/B)));
+
+Q2=-Qrn-Qr2;
+mprintf("\n\n\n(b)\nCapacity of static capacitor= %.1f MVar lagging", -Q2)
diff --git a/1076/CH3/EX3.28/3_28.sce b/1076/CH3/EX3.28/3_28.sce new file mode 100755 index 000000000..ee35b6148 --- /dev/null +++ b/1076/CH3/EX3.28/3_28.sce @@ -0,0 +1,29 @@ +clear;
+clc;
+
+Vr=220
+Vs=240
+
+A=complex(.9*cosd(1), .9*sind(1));
+B=complex(140*cosd(84), 140*sind(84));
+
+b=atan(imag(B)/real(B));
+d=b
+vs=complex (Vs *cos(d), Vs *sin(d));
+P= (((Vs *Vr)/abs(B))*cos(b-d))- real(((Vr)^2 * conj(A/B)));
+mprintf("(a) Max Power Transmitted = %.2f MW", P);
+
+Pr=80;
+Sr=100;
+Qr=sqrt(Sr^2 - Pr^2);
+bd1=round(acos(( Pr+real(((Vr)^2 * conj(A/B))) ) * (abs(B)/(Vs *Vr)))*1000)/1000;
+Qr1= (((Vs *Vr)/abs(B))*sin(bd1))-imag(((Vr)^2 * conj(A/B)));
+
+Q=+Qr-Qr1;
+mprintf("\n\n\n(b)\nCapacity of static capacitor= %.2f MVar leading", Q)
+disp("There is a calculation error in the textbook. 40-49.37=10.63")
+
+d1=b-bd1;
+d1=d1*180/%pi
+mprintf("\n\n\n(c)\nLoad Angle factor= %.2f deg", d1)
+
diff --git a/1076/CH3/EX3.29/3_29.sce b/1076/CH3/EX3.29/3_29.sce new file mode 100755 index 000000000..37ea5c7fe --- /dev/null +++ b/1076/CH3/EX3.29/3_29.sce @@ -0,0 +1,26 @@ +clear;
+clc;
+
+A=complex(.93*cosd(3), .93*sind(3));
+B=complex(150*cosd(70), 150*sind(70));
+D=A;
+
+C=((A*D)-1)/B
+
+Z=complex(100*cosd(70), 100*sind(70));
+Y=complex(.00025*cosd(-75), .00025*sind(-75));
+
+m1=[A B; C D]
+m2=[1 0 ; Y 1 ]
+m3=[1 Z ; 0 1 ]
+
+A0=m1*m2*m3;
+A1=A0(1,1)
+B1=A0(1,2)
+C1=A0(2,1)
+D1=A0(2,2)
+
+mprintf("\nA0= %s", string(round(abs(A1)*1000)/1000) +'/_'+ string(round(atand(imag(A1)/real(A1))*10)/10) )
+mprintf("\nB0= %s", string(round(abs(B1)*1000)/1000) +'/_'+ string(round(atand(imag(B1)/real(B1))*100)/100) )
+mprintf("\nC0= %s", string(round(abs(C1)*100000)/100000) +'/_'+ string(round(atand(imag(C1)/real(C1))*1)/1 ) )
+mprintf("\nD0= %s", string(round(abs(D1)*1000)/1000) +'/_'+ string(round(atand(imag(D1)/real(D1))*10)/10) )
diff --git a/1076/CH3/EX3.3/3_3.sce b/1076/CH3/EX3.3/3_3.sce new file mode 100755 index 000000000..0c08e8582 --- /dev/null +++ b/1076/CH3/EX3.3/3_3.sce @@ -0,0 +1,34 @@ +clear;
+clc;
+
+s1=30e6;
+v1=132e3;
+s2=20e6;
+v2=11e3;
+v3=6.6e3;
+s3=10e3;
+
+xa=.07;
+xb=.09;
+xc=.04;
+va=v1;
+vb=v1;
+vc=v2;
+sa=s1;
+sb=s1;
+sc=s2;
+
+Sb=s1;
+Vb1=v1;
+Vb2=v2;
+Vb3=v3;
+
+Xa=xa * (va/Vb1)^2 * (Sb/sa);
+Xb=xb * (vb/Vb1)^2 * (Sb/sb);
+Xc=xc * (vc/Vb2)^2 * (Sb/sc);
+
+X = [1 1 0; 1 0 1; 0 1 1 ]^(-1) * [ Xa; Xb; Xc];
+
+mprintf("\nX1= %.2f pu",X(1,1));
+mprintf("\nX2= %.2f pu",X(2,1));
+mprintf("\nX3= %.2f pu",X(3,1));
diff --git a/1076/CH3/EX3.30/3_30.sce b/1076/CH3/EX3.30/3_30.sce new file mode 100755 index 000000000..ca3437ff8 --- /dev/null +++ b/1076/CH3/EX3.30/3_30.sce @@ -0,0 +1,14 @@ +clear;
+clc;
+
+f=50;
+l=200;
+Z=complex(14.1, 51.48);
+Y=complex(0, 1.194e-3);
+g=sqrt(Z*Y);
+
+b=imag(g)/l;
+wl=2*%pi/b;
+v=f*wl
+mprintf("\nwavelength = %.4f*1e3 km",wl*1e-3);
+mprintf("\nVelocity of Propagation = %.2f*1e5 km/sec",v*1e-5);
diff --git a/1076/CH3/EX3.31/3_31.sce b/1076/CH3/EX3.31/3_31.sce new file mode 100755 index 000000000..23e537197 --- /dev/null +++ b/1076/CH3/EX3.31/3_31.sce @@ -0,0 +1,55 @@ +clear;
+clc;
+
+clear;
+clc;
+
+
+vr=220e3/sqrt(3);
+Vb=vr
+Vr=vr/Vb
+Sr=40e6;
+Sb=40e6;
+Ib=Sb/(3*Vb)
+Zb=Vb/Ib;
+
+pf=.9;
+pfa=-1*acos(pf);
+Irm=(Sr/(3*vr))/Ib;
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+
+d=500;
+z=complex(.105, .3768);
+y=complex(0, 2.822e-6);
+Z1=z*d;
+Y1=y*d;
+Z=Z1/Zb;
+Y=Y1*Zb;
+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;
+
+
+Vs=(A*Vr)+(B*Ir);
+V=abs(Vs);
+vs=sqrt(3)*V*1e-3*Vb;
+phi1=atan(imag(Vs)/real(Vs))*(180/%pi);
+mprintf("\nSending End Voltage = %.2f kV",vs);
+
+Is=(C*Vr)+(D*Ir);
+I=abs(Is)*Ib;
+phi2=atan(imag(Is)/real(Is))*(180/%pi);
+mprintf("\nSending End Current = %.1f A",I);
+
+phi=phi2-phi1;
+pfs=cosd(phi);
+mprintf("\nSending End Power factor = %.3f ",pfs);
+
+MVA=sqrt(3) * vs* I /1000;
+mprintf("\nSending End Power = %.2f ",MVA);
+
+
diff --git a/1076/CH3/EX3.32/3_32.sce b/1076/CH3/EX3.32/3_32.sce new file mode 100755 index 000000000..9bc6b2e6c --- /dev/null +++ b/1076/CH3/EX3.32/3_32.sce @@ -0,0 +1,33 @@ +clear;
+clc;
+
+VAt1=10
+VAt2=10
+Xt1=.1;
+Xt2=.08
+Vt1a=13.8
+Vt1b=138
+Vt2a=138
+Vt2b=69
+n1=Vt1b/Vt1a
+n2=Vt2b/Vt2a
+
+
+Sb=10
+Vbb=138;
+Vba=Vbb / n1;
+Vbc=Vbb*n2;
+
+Zbc=Vbc^2/Sb;
+
+R=300;
+Rpu=R/Zbc;
+
+v=66
+V=v/Vbc;
+I=V/Rpu;
+Va=V+(I*complex(0,Xt1+Xt2));
+va=abs(Va)*Vba;
+
+mprintf("Voltage in ckt A = %.3f kV", va);
+
diff --git a/1076/CH3/EX3.33/3_33.sce b/1076/CH3/EX3.33/3_33.sce new file mode 100755 index 000000000..f589a8701 --- /dev/null +++ b/1076/CH3/EX3.33/3_33.sce @@ -0,0 +1,41 @@ +clear;
+clc;
+
+xg1=.2
+xg2=.3
+xt1=.2;
+xt2=.06
+zl=complex(40,150);
+
+Vg1=250
+Vg2=250;
+Vt1a=250
+Vt1b=800;
+Vt2a=1000;
+Vt2b=500
+nt1=Vt1b/Vt1a;
+nt2=Vt2b/Vt2a;
+
+sg1=2000;
+sg2=2000;
+st1=4000;
+st2=8000;
+
+Vb1=250;
+Vb2=Vb1*nt1;
+Vb3=Vb2*nt2;
+
+Sb=5000;
+Zbl=Vb2^2/Sb;
+
+Zl=zl/Zbl;
+Xt1=xt1 /( (Vb1/Vt1a)^2 * (st1/Sb));
+Xt2=xt2 / ((Vb2/Vt2a)^2 * (st2/Sb));
+Xg1=xg1 / (sg1/Sb);
+Xg2=xg2 / (sg2/Sb);
+
+ mprintf ("\nReactance of Generator 1= %.1f ",Xg1);
+ mprintf ("\nReactance of Generator 2= %.2f ",Xg2);
+ mprintf ("\nReactance of Transformer 1= %.2f ",Xt1);
+ disp (round(Zl*1e4)/1e4,"Impedance of Line = ");
+ mprintf ("\nReactance of Transformer 2= %.4f ",Xt2);
diff --git a/1076/CH3/EX3.34/3_34.sce b/1076/CH3/EX3.34/3_34.sce new file mode 100755 index 000000000..ea117c64e --- /dev/null +++ b/1076/CH3/EX3.34/3_34.sce @@ -0,0 +1,59 @@ +clear;
+clc;
+
+sg1=10
+sg2=20
+st1=10
+st2=10*3;
+
+vg1=6.6;
+vg2=11.5;
+vt1a=6.6
+vt1b=115
+vt2a=75*sqrt(3);
+vt2b=7.5*sqrt(3);
+nt1=vt1b/vt1a;
+nt2=vt2b/vt2a;
+xg1=.1
+xg2=.1
+xt1=.15
+xt2=.1
+
+
+Sb=20;
+Vb1=6.6;
+Vb2=Vb1*nt1;
+Vb3=Vb2*nt2;
+
+Xg1=xg1*Sb/sg1;
+Xg2=xg2*Sb/sg2;
+Xt1=xt1 * Sb/st1;
+Xt2= xt2 * Sb/st2 * (vt2a/Vb2)^2
+
+ mprintf ("\nReactance of Generator 1= %.1f",Xg1);
+ mprintf ("\nReactance of Generator 2= %.1f",Xg2);
+ mprintf ("\nReactance of Transformer 1= %.1f ",Xt1);
+ mprintf ("\nReactance of Transformer 2= %.5f ",fix(Xt2*1e5)/1e5);
+
+
+Zbl=Vb2^2/Sb;
+
+xab=100;
+xad=100
+xbc=150;
+xbd=200;
+xdc=50;
+
+Xab=xab/Zbl;
+Xad=xad/Zbl;
+Xbc=xbc/Zbl;
+Xbd=xbd/Zbl;
+Xdc=xdc/Zbl;
+
+printf("\nXab=%.4f\t", Xab);
+printf("Xad=%.4f\t", Xad);
+printf("Xbc=%.4f\t", Xbc);
+printf("Xbd=%.4f\t", Xbd);
+printf("Xdc=%.4f\t", Xdc);
+
+
diff --git a/1076/CH3/EX3.35/3_35.sce b/1076/CH3/EX3.35/3_35.sce new file mode 100755 index 000000000..225b577ac --- /dev/null +++ b/1076/CH3/EX3.35/3_35.sce @@ -0,0 +1,61 @@ +clear;
+clc;
+
+sg1=25
+sg2=15
+sg3=30
+st1=30
+st2=15
+st3=10*3;
+
+vg1=6.6;
+vg2=6.6;
+vg3=13.2;
+vt1a=6.6
+vt1b=115
+vt2a=6.6
+vt2b=115
+vt3a=69*sqrt(3);
+vt3b=6.9*sqrt(3);
+nt1=vt1b/vt1a;
+nt2=vt2b/vt2a;
+nt3=vt3b/vt3a;
+
+xg1=.2
+xg2=.15
+xg3=.15
+xt1=.1
+xt2=.1
+xt3=.1
+
+Sb=30;
+Vb1=6.6;
+Vb2=Vb1*nt1;
+Vb3=Vb2*nt2;
+Vb4=Vb2*nt3;
+
+Xg1=xg1*Sb/sg1;
+Xg2=xg2*Sb/sg2;
+Xg3=xg3*Sb/sg3* (vg3/Vb4)^2;
+Xt1=xt1 * Sb/st1;
+Xt2=xt2 * Sb/st2;
+Xt3= xt2 * Sb/st3 * (vt3b/Vb4)^2
+
+ mprintf ("\nReactance of Generator 1= %.2f pu",Xg1);
+ mprintf ("\nReactance of Generator 2= %.1f pu",Xg2);
+ mprintf ("\nReactance of Generator 3= %.4f pu",Xg3);
+ mprintf ("\nReactance of Transformer 1= %.1f pu",Xt1);
+ mprintf ("\nReactance of Transformer 2= %.1f pu",Xt2);
+ mprintf ("\nReactance of Transformer 3= %.3f pu",Xt3);
+
+
+Zbl=Vb2^2/Sb;
+
+xl1=120;
+xl2=90
+
+Xl1=xl1/Zbl;
+Xl2=xl2/Zbl;
+
+ mprintf ("\nReactance of Line 1 = %.4f pu",Xl1);
+ mprintf ("\nReactance of line 2 = %.3f pu",Xl2);
diff --git a/1076/CH3/EX3.36/3_36.sce b/1076/CH3/EX3.36/3_36.sce new file mode 100755 index 000000000..9e527e814 --- /dev/null +++ b/1076/CH3/EX3.36/3_36.sce @@ -0,0 +1,53 @@ +clear;
+clc;
+
+vg=11e3;
+sg=80e6;
+xg=.25;
+
+st1=100e6;
+vt1a=11e3;
+vt1b=220e3;
+nt1=vt1a/vt1b;
+xt1=.05;
+
+st2=150e6;
+vt2a=230e3 ;
+vt2b=33e3;
+nt2=vt2a/vt2b;
+xt2=.04;
+
+r=250;
+
+zl=complex(5,100);
+
+Sb=100e6;
+Vbg=11e3;
+
+Xg=xg * (vg/Vbg)^2 * (Sb/sg);
+Xt1=xt1 * (vt1a/Vbg)^2 * (Sb/st1);
+
+Vbl=Vbg/nt1;
+Zl=zl * (Sb) / (Vbl)^2 ;
+
+Vbm=Vbl/nt2;
+Xt2=xt2 * (vt2a/Vbl)^2 * (Sb/st2);
+R=r * (Sb) / (Vbm)^2 ;
+
+
+Vg=vg/Vbg;
+i=Vg/(Zl+complex(R,Xt1+Xt2+Xg))
+I=round(abs(i)*1e4)/1e4;
+
+Ic=I*Sb/(sqrt(3)*Vbg);
+Il=I*Sb/(sqrt(3)*Vbl);
+Ir=I*Sb/(sqrt(3)*Vbm);
+Vload=Ir*r/1000;
+Vloadll=sqrt(3) * Vload;
+Pr=3*Ir*Ir*r/1e6;
+mprintf("\n Generator Current = %.1f A",Ic);
+mprintf("\n Line Current = %.3f A",Il);
+mprintf("\n Load Current = %.1f A",Ir);
+mprintf("\n Load Voltage = %.2f kV",Vloadll);
+mprintf("\n Load Power = %.3f MW",Pr);
+
diff --git a/1076/CH3/EX3.38/3_38.sce b/1076/CH3/EX3.38/3_38.sce new file mode 100755 index 000000000..02bdaea57 --- /dev/null +++ b/1076/CH3/EX3.38/3_38.sce @@ -0,0 +1,30 @@ +clear;
+clc;
+
+z1=complex(4,6)
+z2=complex(3,2)
+
+Vs=3.3e3/sqrt(3)
+Is=250
+pf=.8
+pfa=acos(pf)
+I=Is *(exp(%i * -pfa))
+
+I1=I * z2/(z1+z2)
+pfa1=atan(imag(I1)/real(I1))
+pf1=cos(pfa1)
+mprintf("\n(a) Current in OH line = %.1f A pf= %.3f", abs(I1), pf1)
+
+I2=I * z1/(z1+z2)
+pfa2=atan(imag(I2)/real(I2))
+pf2=cos(pfa2)
+mprintf("\n(b) Current in cable = %.2f A pf= %.2f", abs(I2), pf2)
+
+vr=sqrt((Vs)^2-imag(I1*z1)^2)- real(I1*z1)
+Vr=vr*sqrt(3)/1000;
+mprintf("\n(c) Receiving end voltage = %.3f KV", Vr)
+
+d=atan(imag(I1*z1)/(Vr+real(I1*z1)))
+phi=pfa-d;
+pfr=cos(phi)
+mprintf("\n(d) Receiving end pf = %.1f lagging", pfr)
diff --git a/1076/CH3/EX3.39/3_39.sce b/1076/CH3/EX3.39/3_39.sce new file mode 100755 index 000000000..8666498a7 --- /dev/null +++ b/1076/CH3/EX3.39/3_39.sce @@ -0,0 +1,23 @@ +clear;
+clc;
+
+l=300
+R=.4 *3
+X=.8*3
+Vs=11e3/sqrt(3);
+P=3000;
+pf=.8
+pfa=acos(pf)
+VIr=P/(3*pf)
+
+a=1;
+b=-Vs
+c=VIr * 1e3 * ((R*cos(pfa))+(X*sin(pfa)))
+vr=(-b+sqrt((b*b )- (4*a*c)))/(2*a)
+Ir=VIr*1e3/vr;
+Vr=vr*sqrt(3)/1000;
+mprintf("\nReceiving End Voltage = %.2f kV",Vr)
+
+Pl=3* (Ir)^2 * R/ 1000;
+eff=P*100/(P+Pl)
+mprintf("\nefficiency = %.2f percent",eff)
diff --git a/1076/CH3/EX3.4/3_4.sce b/1076/CH3/EX3.4/3_4.sce new file mode 100755 index 000000000..8ee6cf655 --- /dev/null +++ b/1076/CH3/EX3.4/3_4.sce @@ -0,0 +1,34 @@ +clear;
+clc;
+
+d=15;
+Vr=11e3/sqrt(3);
+pfr=.8;
+Pd=5e6;
+Pl=.12*Pd;
+l=1.1e-3;
+L=l*d;
+
+I= Pd/(3*pfr*Vr);
+R=Pl/(3*I*I);
+X=2 * %pi* 50 *L;
+pfa=acos(pfr);
+Vs=Vr + (I * R * pfr) + (I * X * sin(pfa));
+vs=sqrt(3)*Vs;
+
+VR=(Vs-Vr)/Vr;
+
+mprintf("\n(a) Voltage Regulation = %.2f percent ", VR*100);
+
+pfa0=atan(R/X);
+pf0=cos(pfa0);
+
+mprintf("\n(b) pf at VR=0 = %.3f ", pf0);
+
+I0= (I* pfr)/pf0;
+Ic= (I * sin(pfa))+(I0*sin(pfa0));
+Xc=Vr/Ic
+C=1/(100*%pi*Xc);
+
+mprintf("\n(c) C = %.1f e-6 F", C*1e6);
+
diff --git a/1076/CH3/EX3.5/3_5.sce b/1076/CH3/EX3.5/3_5.sce new file mode 100755 index 000000000..64ef78abb --- /dev/null +++ b/1076/CH3/EX3.5/3_5.sce @@ -0,0 +1,22 @@ +clear;
+clc;
+
+Vs=11e3/sqrt(3);
+Pd=1200e3;
+R=5.31;
+X=5.54;
+pfr=.8;
+
+
+VIr= Pd/(3*pfr);
+a=1;
+b=-1*Vs;
+pfa=acos(pfr)
+c=(VIr * R * pfr) + (VIr * X * sin(pfa))
+
+Vr=(-b+sqrt(b^2 - (4*a*c)))/(2*a);
+
+I=VIr/Vr;
+vr=sqrt(3)*Vr;
+
+mprintf("Recieving end Voltage = %.3f KV and Current = %.2f A ",vr/1000,I);
diff --git a/1076/CH3/EX3.6/3_6.sce b/1076/CH3/EX3.6/3_6.sce new file mode 100755 index 000000000..cde795ff0 --- /dev/null +++ b/1076/CH3/EX3.6/3_6.sce @@ -0,0 +1,30 @@ +clear;
+clc;
+
+Pd=1e6;
+pf=.8;
+v1=30e3;
+v2=10e3;
+Rl=25;
+Xl=12;
+rt=.8;
+xt=2.5;
+n=v1/v2;
+
+Rt=rt*(n^2);
+Xt=xt*(n^2);
+
+R=Rt+Rl;
+X=Xt+Xl;
+
+Vr=v1/sqrt(3);
+I=Pd/(3*Vr*.8);
+
+pfa=acos(pf)
+Vs=Vr + (I * R * pf) + (I * X * sin(pfa));
+vs=sqrt(3)*Vs*1e-3;
+
+VR=(Vs-Vr)/Vr;
+
+mprintf("\nSending End Voltage = %.2f KV",vs);
+mprintf("\nVoltage Regulation= %.2f",VR*100);
diff --git a/1076/CH3/EX3.7/3_7.sce b/1076/CH3/EX3.7/3_7.sce new file mode 100755 index 000000000..95ad2abe5 --- /dev/null +++ b/1076/CH3/EX3.7/3_7.sce @@ -0,0 +1,19 @@ +clear;
+clc;
+
+Vs=33e3/sqrt(3);
+Vr=30e3/sqrt(3);
+Pr=10e6;
+pf=.8;
+eff=.96;
+
+I=Pr/(3*Vr*pf);
+
+Ps=Pr/eff;
+Pl=Ps-Pr;
+
+R=Pl/(3*I*I);
+pfa=acos(pf);
+X=((Vs-Vr)-(I*R*pf))/(I*sin(pfa));
+
+mprintf("R= %.1f ohm per phase, X= %.1f ohm per phase",R,X);
diff --git a/1076/CH3/EX3.8/3_8.sce b/1076/CH3/EX3.8/3_8.sce new file mode 100755 index 000000000..9dc8a3004 --- /dev/null +++ b/1076/CH3/EX3.8/3_8.sce @@ -0,0 +1,33 @@ +clear;
+clc;
+
+R=48.7;
+X=80.2;
+Z=complex(R,X);
+c=8.42e-9;
+l=200;
+C=c*l;
+Y=complex(0,(C* 100*%pi));
+
+Vr=88e3/sqrt(3);
+Pr=13.5e6;
+pf=.9;
+pfa=-1* acos(pf);
+Irm=Pr/(3*Vr*pf);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+Vs=(Vr*(1+((Z*Y)/2)))+ (Ir*Z*(1+(Z*Y/4)));
+V=abs(Vs);
+vs=sqrt(3)*V*1e-3;
+phi=atan(imag(Vs)/real(Vs))*(180/%pi);
+mprintf("\nSending End Voltage = %.2f kV",vs);
+mprintf("\nSending End Power Angle = %.1f deg ",phi);
+
+Is=(Vr*Y)+(Ir*(1+(Y*Z/2)));
+I=abs(Is);
+mprintf("\nSending End Current = %.2f A",I);
+
+Vr0=V/(1+(Y*Z/2));
+V0=abs(Vr0);
+VR=(V0-Vr)/Vr;
+mprintf("\nVoltage Regulation = %.1f ",VR*100);
diff --git a/1076/CH3/EX3.9/3_9.sce b/1076/CH3/EX3.9/3_9.sce new file mode 100755 index 000000000..098790ab8 --- /dev/null +++ b/1076/CH3/EX3.9/3_9.sce @@ -0,0 +1,33 @@ +clear;
+clc;
+
+R=48.7;
+X=80.2;
+Z=complex(R,X);
+c=8.42e-9;
+l=200;
+C=c*l;
+Y=complex(0,(C* 100*%pi));
+
+Vr=88e3/sqrt(3);
+Pr=13.5e6;
+pf=.9;
+pfa=-1* acos(pf);
+Irm=Pr/(3*Vr*pf);
+Ir=complex(Irm *pf, Irm * sin(pfa));
+
+Vs=(Vr*(1+((Z*Y)/2)))+ (Ir*Z);
+V=abs(Vs);
+vs=sqrt(3)*V*1e-3;
+phi=atan(imag(Vs)/real(Vs))*(180/%pi);
+mprintf("\nSending End Voltage = %.2f kV",vs);
+mprintf("\nSending End Power Angle = %.2f deg ",phi);
+
+Is=(Vr*Y*(1+(Z*Y/4)))+(Ir*(1+(Y*Z/2)));
+I=abs(Is);
+mprintf("\nSending End Current = %.2f A",I);
+
+Vr0=V/(1+(Y*Z/2));
+V0=abs(Vr0);
+VR=(V0-Vr)/Vr;
+mprintf("\nVoltage Regulation = %.2f ",VR*100);
|