diff options
Diffstat (limited to '1076/CH10')
30 files changed, 1229 insertions, 0 deletions
diff --git a/1076/CH10/EX10.1/10_1.sce b/1076/CH10/EX10.1/10_1.sce new file mode 100755 index 000000000..edd415b13 --- /dev/null +++ b/1076/CH10/EX10.1/10_1.sce @@ -0,0 +1,46 @@ +clear
+clc
+
+Sg1=100
+Vg1=11
+xg1=.15 *%i
+
+Sg2=50
+Vg2=11
+xg2=.1 *%i
+
+St1=100
+Vt1a=11
+Vt1b=132
+xt1=.1 *%i
+nt1=Vt1b/Vt1a
+
+St2=50
+Vt2a=11
+Vt2b=132
+xt2=.08 *%i
+nt2=Vt2b/Vt2a
+
+Sb=100
+Vb1=11
+Vb2=nt1*Vb1
+
+xl=.2 * 200 *%i
+Xl=xl/(Vb2*Vb2/Sb)
+
+Xg2=xg2*Sb/Sg2
+Xt2=xt2 * Sb/St2
+
+X=( ((xg1 +xt1) * (Xg2 +Xt2) )/ ((xg1 +xt1) + (Xg2 +Xt2) )) +Xl/2
+
+I=1/X
+
+Ib1=Sb*1e3/(Vb1*sqrt(3))
+Ib2=Sb*1e3/(Vb2*sqrt(3))
+If=abs(I*Ib2)
+Ifg=abs(I*Ib1)
+Ifg1=Ifg * (Xg2 +Xt2)/(xg1+xt1+Xt2+Xg2)
+Ifg2=Ifg * (xg1 +xt1)/(xg1+xt1+Xt2+Xg2)
+MVAf=abs(I*1*Sb)
+
+mprintf("Total fault current = %.2f A, Fault Level= %f MVA,\n Fault current supplied by generator 1=%f A, generator 2=%f A",If, MVAf, Ifg1, Ifg2)
diff --git a/1076/CH10/EX10.10/10_10.sce b/1076/CH10/EX10.10/10_10.sce new file mode 100755 index 000000000..c2c38f543 --- /dev/null +++ b/1076/CH10/EX10.10/10_10.sce @@ -0,0 +1,48 @@ +clear
+clc
+
+Z=[.2 .2 .05]' * %i
+S=30
+V=11
+I=S*1e6/(sqrt(3)*V*1e3)
+E=1
+
+Ia1=E/(Z(1)+Z(2)+Z(3))
+If= 3*abs(Ia1) * S*1e6/(sqrt(3)* V*1e3)
+Ia2=Ia1
+Ia0=Ia1
+
+a=exp(%i * 2 * %pi/3)
+A=[1 1 1
+1 a^2 a
+1 a a^2
+]
+
+Va1=E-Ia1*Z(1)
+Va2=0-Ia2*Z(2)
+Va0=0-Ia0*Z(3)
+
+Vp=[ Va0 Va1 Va2]'
+v=A*Vp
+
+vab=v(1)-v(2)
+vbc=v(2)-v(3)
+vca=v(3)-v(1)
+
+Vbll=V/sqrt(3)
+
+Vab=vab * Vbll
+Vbc=vbc * Vbll
+Vca=vca * Vbll
+Vll=[Vab Vbc Vca]
+mprintf("\n(a)\n line currents: If= %f A, Line voltages in KV ", If)
+mprintf("\nVab= %s", string(round(abs(Vll(1))*10)/10) +'/_'+ string(round(atand(imag(Vll(1))/real(Vll(1)))*10)/10) )
+mprintf("\nVbc= %s", string(round(abs(Vll(2))*10)/10) +'/_'+ string(round(atand(imag(Vll(2))/real(Vll(2)))*10)/10 +180) )
+mprintf("\nVca= %s", string(round(abs(Vll(3))*10)/10) +'/_'+ string(round(atand(imag(Vll(3))/real(Vll(3)))*10)/10 +180))
+
+
+
+If3=E/Z(1)
+if3=abs(If3)* S*1e6/(sqrt(3)* V*1e3)
+
+mprintf("\n(b)3 phase fault current is -%.0fi A", if3)
diff --git a/1076/CH10/EX10.11/10_11.sce b/1076/CH10/EX10.11/10_11.sce new file mode 100755 index 000000000..8860eb509 --- /dev/null +++ b/1076/CH10/EX10.11/10_11.sce @@ -0,0 +1,42 @@ +clear
+clc
+
+X=[
+ .25 .25 .05
+ .2 .2 .05
+ .06 .06 .06
+ .07 .07 .07
+ .1 .1 .3
+ .1 .1 .3
+]
+
+B=[
+ 100 11
+ 100 11
+ 100 11
+ 100 11
+ 100 220
+ 100 220
+ ]
+V1=11
+V2=220
+S=100
+Xe=3*.03
+//end 9
+X1=(((X(1,1)*B(1,1) *V1/(S*B(1,2))) +(X(3,1)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,1)*B(2,1) *V1/(S*B(2,2)))+(X(4,1)*B(4,1) *V1/(S*B(4,2))) +((X(5,1)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,1)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X2=(((X(1,2)*B(1,1) *V1/(S*B(1,2))) +(X(3,2)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,2)*B(2,1) *V1/(S*B(2,2)))+(X(4,2)*B(4,1) *V1/(S*B(4,2))) +((X(5,2)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,2)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X0=((X(3,3)*B(3,1) *V1/(S*B(3,2)))^-1 + ((Xe *B(4,1) *V1/(S*B(4,2))) + (X(2,3)*B(2,1) *V1/(S*B(2,2))) + (X(4,3)*B(4,1) *V1/(S*B(4,2))) +(((X(5,3)*B(5,1) *V2/(S*B(5,2))))^-1 +(X(6,3)*B(6,1) *V2/(S*B(6,2)))^-1)^-1 )^-1 )^-1
+
+Z1=%i * round(X1*1e3)/1e3
+Z2=%i * round(X2*1e3)/1e3
+Z0=%i * round(X0*1e3)/1e3
+
+Z=[Z1 Z2 Z0]'
+
+mprintf("\nZ1= %.3fj\n", imag(Z(1)))
+mprintf("Z2= %.3fj\n", imag(Z(2)))
+mprintf("Z0= %.3fj\n", imag(Z(3)))
+
+
diff --git a/1076/CH10/EX10.12/10_12.sce b/1076/CH10/EX10.12/10_12.sce new file mode 100755 index 000000000..c1214af7b --- /dev/null +++ b/1076/CH10/EX10.12/10_12.sce @@ -0,0 +1,139 @@ +clear
+clc
+
+X=[
+ .25 .25 .05
+ .2 .2 .05
+ .06 .06 .06
+ .07 .07 .07
+ .1 .1 .3
+ .1 .1 .3
+]
+
+B=[
+ 100 11
+ 100 11
+ 100 11
+ 100 11
+ 100 220
+ 100 220
+ ]
+V1=11
+V2=220
+S=100
+Xe=3*.03
+//end 9
+X1=(((X(1,1)*B(1,1) *V1/(S*B(1,2))) +(X(3,1)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,1)*B(2,1) *V1/(S*B(2,2)))+(X(4,1)*B(4,1) *V1/(S*B(4,2))) +((X(5,1)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,1)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X2=(((X(1,2)*B(1,1) *V1/(S*B(1,2))) +(X(3,2)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,2)*B(2,1) *V1/(S*B(2,2)))+(X(4,2)*B(4,1) *V1/(S*B(4,2))) +((X(5,2)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,2)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X0=((X(3,3)*B(3,1) *V1/(S*B(3,2)))^-1 + ((Xe *B(4,1) *V1/(S*B(4,2))) + (X(2,3)*B(2,1) *V1/(S*B(2,2))) + (X(4,3)*B(4,1) *V1/(S*B(4,2))) +(((X(5,3)*B(5,1) *V2/(S*B(5,2))))^-1 +(X(6,3)*B(6,1) *V2/(S*B(6,2)))^-1)^-1 )^-1 )^-1
+
+Z1=%i * round(X1*1e3)/1e3
+Z2=%i * round(X2*1e3)/1e3
+Z0=%i * round(X0*1e3)/1e3
+Z=[Z1 Z2 Z0]'
+//end ques 11
+
+a=exp(%i * 2 * %pi /3)
+A= [1 1 1 ; 1 a^2 a; 1 a a^2]
+
+//12(a)
+Ia1a=1/(Z1+Z2+Z0)
+Ia1a=round(Ia1a *1e2)/1e2
+mIa1a=[ Ia1a Ia1a Ia1a]'
+mIa=A*mIa1a
+Ia=round(mIa(1)*100)/100
+Iaa = round(abs(mIa(1))*100)/100
+Iba=round(S*1e7/(sqrt(3) * V2 * 1e3))/10
+IFa=round(Iba * Iaa *100)/100
+
+mprintf("\n(a) Fault current = %.2f A",IFa)
+
+//12(b)
+
+Va1=round((1- (Z1 * Ia1a))*100)/100
+Va2=round((0- (Z2 * Ia1a))*100)/100
+Va0=round((0- (Z0 * Ia1a))*100)/100
+mVa1=[ Va0 Va1 Va2]'
+mVa=A*mVa1
+v=mVa * V2 / sqrt(3)
+v=round(v *10000)/10000
+
+mprintf("\n\n(b)line to neutral voltages in KV ")
+mprintf("\nVa= %s", string(fix(abs(v(1))*100)/100) )
+mprintf("\nVb= %s", string(fix(abs(v(2))*100)/100) +'/_'+ string(round(atand(imag(v(2))/real(v(2)))*100)/100 +180 ) )
+mprintf("\nVc= %s", string(fix(abs(v(3))*100)/100) +'/_'+ string(round(atand(imag(v(3))/real(v(3)))*100)/100 +180))
+
+//12(c) --- g2
+
+Ia1g2= Ia1a * (((X(1,1)*B(1,1) *V1/(S*B(1,2))) +(X(3,1)*B(3,1) *V1/(S*B(3,2))) ))/((((X(1,1)*B(1,1) *V1/(S*B(1,2))) +(X(3,1)*B(3,1) *V1/(S*B(3,2))) )+((X(2,1)*B(2,1) *V1/(S*B(2,2)))+(X(4,1)*B(4,1) *V1/(S*B(4,2))) +((X(5,1)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,1)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)))
+Ia1g2=round(Ia1g2 *1e2)/1e2
+
+Ia2g2= Ia1a * (((X(1,2)*B(1,1) *V1/(S*B(1,2))) +(X(3,2)*B(3,1) *V1/(S*B(3,2))) ))/((((X(1,2)*B(1,1) *V1/(S*B(1,2))) +(X(3,2)*B(3,1) *V1/(S*B(3,2))) )+((X(2,2)*B(2,1) *V1/(S*B(2,2)))+(X(4,2)*B(4,1) *V1/(S*B(4,2))) +((X(5,2)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,2)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)))
+Ia2g2=round(Ia2g2 *1e2)/1e2
+
+Ia0g2= Ia1a * (((X(3,3)*B(3,1) *V1/(S*B(3,2))) ))/((((Xe *B(4,1) *V1/(S*B(1,2))) +(X(3,3)*B(3,1) *V1/(S*B(3,2))) )+(X(2,3)*B(2,1) *V1/(S*B(2,2)))+(X(4,3)*B(4,1) *V1/(S*B(4,2))) +((X(5,3)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,3)*B(6,1) *V2/(S*B(6,2)))^-1)^-1))
+Ia0g2=round(Ia0g2 *1e2)/1e2
+
+mIa1g2=[ Ia0g2 Ia1g2 Ia2g2]'
+mIag2=A*mIa1g2
+Ibc=round(S*1e7/(sqrt(3) * V1 * 1e3))/10
+Iag2=abs(mIag2) * Ibc
+
+
+mprintf("\n\n(c)line currents in A at generator 2 ")
+mprintf("\nIa= %s", string(round(abs(Iag2(1))*10)/10) )
+mprintf("\nIb= %s", string(round(abs(Iag2(2))*10)/10) )
+mprintf("\nIc= %s", string(round(abs(Iag2(3))*10)/10) )
+
+
+Va1g2=round((1- (X(2,1) *%i * Ia1g2))*100)/100
+Va2g2=round((0- (X(2,2) *%i * Ia2g2))*100)/100
+Va0g2=round((0- ((X(2,3) +Xe) *%i * Ia0g2))*10000)/10000
+mVa1g2=[ Va0g2 Va1g2 Va2g2]'
+mVag2=A*mVa1g2
+vg2=mVag2 * V1 / sqrt(3)
+
+mprintf("\n\nline to neutral voltages in KV at generator 2 ")
+mprintf("\nVa= %s", string(fix(abs(vg2(1))*100)/100) )
+mprintf("\nVb= %s", string(fix(abs(vg2(2))*100)/100) )
+mprintf("\nVc= %s", string(fix(abs(vg2(3))*100)/100) )
+
+//12(c) --- g1
+
+Ia1g1= (Ia1a-Ia1g2) * exp(%i * 1*%pi /6) *-1
+Ia1g1=round(Ia1g1 *1e2)/1e2
+
+Ia2g1= (Ia1a-Ia2g2)* exp(%i *-1*%pi /6) *-1
+Ia2g1=round(Ia2g1 *1e2)/1e2
+
+Ia0g1= 0
+
+mIa1g1=[ Ia0g1 Ia1g1 Ia2g1]'
+mIag1=A*mIa1g1
+mIag1=round(mIag1*1e1)/1e1
+Ibc=round(S*1e7/(sqrt(3) * V1 * 1e3))/10
+Iag1=abs(mIag1) * Ibc
+
+
+mprintf("\n\nline currents in A at generator 1")
+mprintf("\nIa= %s", string(round(abs(Iag1(1))*1)/1) )
+mprintf("\nIb= %s", string(round(abs(Iag1(2))*1)/1) )
+mprintf("\nIc= %s", string(round(abs(Iag1(3))*1)/1) )
+
+
+Va1g1=round((1- (X(1,1) *%i * (Ia1a-Ia1g2)))*10000)/10000
+Va1g1=Va1g1 * exp(%i * 1*%pi /6)
+Va2g1=round((0- (X(1,2) *%i * (Ia1a-Ia2g2)))*10000)/10000
+Va2g1=Va2g1 * exp(%i *- 1*%pi /6)
+Va0g1=0
+mVa1g1=[ Va0g1 Va1g1 Va2g1]'
+mVag1=A*mVa1g1
+vg1=mVag1 * V1 / sqrt(3)
+
+mprintf("\n\nline to neutral voltages in KV at generator 1 ")
+mprintf("\nVa= %s", string(fix(abs(vg1(1))*100)/100) )
+mprintf("\nVb= %s", string(fix(abs(vg1(2))*100)/100) )
+mprintf("\nVc= %s", string(fix(abs(vg1(3))*100)/100) )
+
diff --git a/1076/CH10/EX10.13/10_13.sce b/1076/CH10/EX10.13/10_13.sce new file mode 100755 index 000000000..a35cf9325 --- /dev/null +++ b/1076/CH10/EX10.13/10_13.sce @@ -0,0 +1,18 @@ +clear
+clc
+
+Sb=37.5
+Vb=33
+Zb=Vb*Vb/Sb
+Ib= Sb *1e6 / (sqrt(3) *Vb *1e3)
+x1=[.18 .12 .1]
+x2=[6.3 6.3 12.6]
+X2=x2/Zb
+
+X=x1+X2
+
+x=X(1)+X(2)+X(3)
+
+If=3*1/x
+IF=If*Ib
+mprintf("Fault current is %.1f A", IF)
diff --git a/1076/CH10/EX10.14/10_14.sce b/1076/CH10/EX10.14/10_14.sce new file mode 100755 index 000000000..878568275 --- /dev/null +++ b/1076/CH10/EX10.14/10_14.sce @@ -0,0 +1,41 @@ +clear
+clc
+
+Z=[.2 .2 .05]' * %i
+S=30
+V=11
+I=S*1e6/(sqrt(3)*V*1e3)
+E=1
+
+Ia1=E/(Z(1)+Z(2))
+Ia2=-Ia1
+Ia0=0
+Ibase=S*1e6/(sqrt(3)* V*1e3)
+
+Ia=0
+Ib=sqrt(3)*Ia1*Ibase
+Ic=-Ib
+
+mprintf("\nLine currents: (Ia Ib Ic) in Ampere")
+mprintf("\nIa= %d", Ia)
+mprintf("\nIb= -%d",abs(Ib))
+mprintf("\nIc= %d", abs(Ic))
+
+a=exp(%i * 2 * %pi/3)
+A=[1 1 1
+1 a^2 a
+1 a a^2
+]
+
+Va1=E-Ia1*Z(1)
+Va2=0-Ia2*Z(2)
+Va0=0-Ia0*Z(3)
+
+Vbll=V/sqrt(3)
+Vp=[ Va0 Va1 Va2]'
+v=A*Vp*Vbll
+
+mprintf("\nline to neutral voltages in KV ")
+mprintf("\nVa= %s", string(round(abs(v(1))*1000)/1000) +'/_'+ string(round(atand(imag(v(1))/real(v(1)))*10)/10) )
+mprintf("\nVb= %s", string(round(abs(v(2))*1000)/1000) +'/_'+ string(round(atand(imag(v(2))/real(v(2)))*10)/10 +180) )
+mprintf("\nVc= %s", string(round(abs(v(3))*1000)/1000) +'/_'+ string(round(atand(imag(v(3))/real(v(3)))*10)/10 +180))
diff --git a/1076/CH10/EX10.15/10_15.sce b/1076/CH10/EX10.15/10_15.sce new file mode 100755 index 000000000..dda24d54d --- /dev/null +++ b/1076/CH10/EX10.15/10_15.sce @@ -0,0 +1,45 @@ +clear
+clc
+
+clear
+clc
+
+Z=[.2 .2 .05]' * %i
+S=30
+V=11
+I=round(S*1e8/(sqrt(3)*V*1e3) )/1e2
+E=1
+
+Ia1=E/(Z(1)+(Z(2)*Z(3)/(Z(2)+Z(3))))
+Ia1=round(Ia1*1000)/1000
+Ia2=-Ia1 * Z(3)/(Z(2)+Z(3))
+Ia0=-Ia1 * Z(2)/(Z(2)+Z(3))
+
+a=exp(%i * 2 * %pi/3)
+A=[1 1 1
+1 a^2 a
+1 a a^2
+]
+
+Ia=A*[ Ia0 Ia1 Ia2]' * I
+
+mprintf("Line currents: ")
+mprintf("\nIa= %.2f ang(0) A",abs(Ia(1)))
+mprintf("\nIb= %.2f ang(%.2f) A",abs(Ia(2)),atand(imag(Ia(2))/real(Ia(2)))+180)
+mprintf("\nIc= %.2f ang(%.2f) A",abs(Ia(3)),atand(imag(Ia(3))/real(Ia(3))))
+
+
+If=Ia(2)+Ia(3)
+mprintf("\nFault current= %.0f ang(%.2f) A",abs(If),atand(imag(If)/real(If)))
+
+Va1=1-(Ia1*Z(1))
+Va2=Va1
+Va0=Va1
+
+Va=A*[Va0 Va1 Va2]' * V/sqrt(3)
+mprintf("\nLine to neutral voltages: ")
+mprintf("\nVa= %.3f ang(%.2f) kV",abs(Va(1)),atand(imag(Va(1))/real(Va(1))))
+mprintf("\nVb= %.3f KV",abs(Va(2)))
+mprintf("\nVc= %.3f KV",abs(Va(3)))
+
+
diff --git a/1076/CH10/EX10.16/10_16.sce b/1076/CH10/EX10.16/10_16.sce new file mode 100755 index 000000000..1ce736714 --- /dev/null +++ b/1076/CH10/EX10.16/10_16.sce @@ -0,0 +1,58 @@ +clear
+clc
+
+X=[
+ .25 .25 .05
+ .2 .2 .05
+ .06 .06 .06
+ .07 .07 .07
+ .1 .1 .3
+ .1 .1 .3
+]
+
+B=[
+ 100 11
+ 100 11
+ 100 11
+ 100 11
+ 100 220
+ 100 220
+ ]
+V1=11
+V2=220
+S=100
+Xe=3*.03
+//end 9
+X1=(((X(1,1)*B(1,1) *V1/(S*B(1,2))) +(X(3,1)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,1)*B(2,1) *V1/(S*B(2,2)))+(X(4,1)*B(4,1) *V1/(S*B(4,2))) +((X(5,1)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,1)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X2=(((X(1,2)*B(1,1) *V1/(S*B(1,2))) +(X(3,2)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,2)*B(2,1) *V1/(S*B(2,2)))+(X(4,2)*B(4,1) *V1/(S*B(4,2))) +((X(5,2)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,2)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X0=((X(3,3)*B(3,1) *V1/(S*B(3,2)))^-1 + ((Xe *B(4,1) *V1/(S*B(4,2))) + (X(2,3)*B(2,1) *V1/(S*B(2,2))) + (X(4,3)*B(4,1) *V1/(S*B(4,2))) +(((X(5,3)*B(5,1) *V2/(S*B(5,2))))^-1 +(X(6,3)*B(6,1) *V2/(S*B(6,2)))^-1)^-1 )^-1 )^-1
+
+Z1=%i * round(X1*1e3)/1e3
+Z2=%i * round(X2*1e3)/1e3
+Z0=%i * round(X0*1e3)/1e3
+
+//end q11
+
+a=exp(%i * 2 * %pi /3)
+A= [1 1 1 ; 1 a^2 a; 1 a a^2]
+
+
+Ia1=1/(Z1 + (Z2*Z0/(Z2+Z0)))
+Ia1=round(Ia1 *1e3)/1e3
+Ia2=(Z0/(Z2+Z0)) * Ia1 *-1
+Ia2=round(Ia2 *1e3)/1e3
+Ia0=(Z2/(Z2+Z0)) * Ia1 *-1
+Ia0=round(Ia0 *1e3)/1e3
+
+mIa1=[ Ia0 Ia1 Ia2]'
+mIa=A*mIa1
+Ib=round(S*1e7/(sqrt(3) * V2 * 1e3))/10
+Ia=round(abs(mIa)) * Ib
+
+
+mprintf("\nline currents are: in Amperes ")
+mprintf("\nIa= %s", string(round(abs(Ia(1))*10)/10) )
+mprintf("\nIb= %s", string(round(abs(Ia(2))*10)/10) +'/_'+ string(round(atand(imag(mIa(2)*Ib)/real(mIa(2)*Ib))*100)/100 +180 ))
+mprintf("\nIc= %s", string(round(abs(Ia(3))*10)/10) +'/_'+ string(round(atand(imag(mIa(3)*Ib)/real(mIa(3)*Ib))*100)/100 ))
diff --git a/1076/CH10/EX10.17/10_17.sce b/1076/CH10/EX10.17/10_17.sce new file mode 100755 index 000000000..712f09f74 --- /dev/null +++ b/1076/CH10/EX10.17/10_17.sce @@ -0,0 +1,19 @@ +clear
+clc
+
+S=50
+V=11
+Z=V*V/S
+
+If1=1870
+If2=2590
+If3=4130
+
+X1=V*1e3 /(sqrt(3)*If1)
+x1=X1/Z
+X2=((V*1e3*sqrt(3)/(sqrt(3)) )/(If2))- X1
+x2=X2/Z
+X0=((sqrt(3)* V *1e3)/If3) - X1-X2
+x0=X0/Z
+
+mprintf("x1= %.2f pu \n x2= %.2f pu \n x0= %.2f pu", x1, x2,x0)
diff --git a/1076/CH10/EX10.18/10_18.sce b/1076/CH10/EX10.18/10_18.sce new file mode 100755 index 000000000..9ada54a36 --- /dev/null +++ b/1076/CH10/EX10.18/10_18.sce @@ -0,0 +1,46 @@ +clear
+clc
+
+Vt=10.95
+V=11
+S=100
+I=round(S*1e7/(sqrt(3)*V*1e3))/10 //Error in evaluation of base current in textbook
+
+vt=Vt/V
+pf=.8
+P=40
+
+Il=round(P*1e7/(sqrt(3)*Vt*pf*1e3))/10
+Il=Il * exp (%i * acos(pf))/I
+
+x1g=.2
+x2g=.2
+x0g=.05
+
+Sm=50
+x1m=.2 * S/Sm
+x2m=.2 * S/Sm
+x0m=.05
+
+xt1=.05
+xt2=.05
+xt0=.15
+
+
+V=vt+ (Il * %i * x0m)
+x0m=.05* S/Sm
+Ia1=V/(%i*((x1g*(x1m+xt1)/(x1g+x1m+xt1))+(x2g*(x2m+xt2)/(x2g+x2m+xt2))+(x0g*(x0m+xt0)/(x0g+x0m+xt0))))
+Ia1=round(Ia1 * 1000)/1000
+Ig1=round((Ia1 * (x1m + xt1)/(x1m+x1g+xt1) + Il)*1000)/1000
+Im1=round((Ia1 * (x1g)/(x1m+x1g+xt1) - Il)*1000)/1000
+Ig2=round((Ia1 * (x2m + xt2)/(x2m+x2g+xt2))*1000)/1000
+Im2=round((Ia1 * (x2g)/(x2m+x2g+xt2))*1000)/1000
+Ig0=round((Ia1 * (x0m + xt0)/(x0m+x0g+xt0))*100)/100
+Im0=round((Ia1 * (x0g)/(x0m+x0g+xt0))*100)/100
+
+Im=round((Im1+Im2+Im0)*1000) *1e-3
+Ig=round((Ig1+Ig2+Ig0)*1000) *1e-3
+
+mprintf("\nCurrent Through motor = %.2f, ang (%.1f) deg", abs(Im)*I, atand(imag(Im)/real(Im))+180)
+mprintf("\nCurrent Through generator = %.2f, ang (%.1f) deg", abs(Ig)*I, atand(imag(Ig)/real(Ig)))
+disp("Error in evaluation of base current in textbook")
diff --git a/1076/CH10/EX10.19/10_19.sce b/1076/CH10/EX10.19/10_19.sce new file mode 100755 index 000000000..bc9ebdc57 --- /dev/null +++ b/1076/CH10/EX10.19/10_19.sce @@ -0,0 +1,18 @@ +clear
+clc
+
+S=100
+V=11
+Z=V*V/S
+
+SCA=1000
+SCB=650
+
+xa=S/SCA
+xb=S/SCB
+Xc=.5
+xc=Xc/Z
+
+X=round(((xa *(xb+xc))/(xa+xb+xc))*1000)/1000
+FS=S/X
+mprintf("Fault MVA= %.2f MVA",FS)
diff --git a/1076/CH10/EX10.2/10_2.sce b/1076/CH10/EX10.2/10_2.sce new file mode 100755 index 000000000..cfb8b26cb --- /dev/null +++ b/1076/CH10/EX10.2/10_2.sce @@ -0,0 +1,24 @@ +clear
+clc
+
+s=50
+v=11
+x=.15
+S=50
+V=11
+
+Xe=.15/4
+FMVA= round(10000/Xe)/10000
+Fault=FMVA * S
+mprintf("\n(a)fault level = %.3f pu", Fault)
+
+sga=.5*Fault
+sgb=800-sga
+
+xb=.15/2
+X=(S/sgb)-xb
+x=X*(V*V/S)
+mprintf("\n(b)X = %.3f ohms", x)
+
+
+
diff --git a/1076/CH10/EX10.20/10_20.sce b/1076/CH10/EX10.20/10_20.sce new file mode 100755 index 000000000..e8bda27c2 --- /dev/null +++ b/1076/CH10/EX10.20/10_20.sce @@ -0,0 +1,19 @@ +clear
+clc
+
+V=33
+S=75
+Z=V*V/S
+
+sg=15
+xg=.15
+Xg=xg *S/sg
+xt=.08
+x=Xg/3
+
+CB=750
+XF=S/CB
+xi=((x*xt)-(XF*(x+xt)))/(XF-x)
+xi=round(xi*10000)/10000
+X=xi*Z
+mprintf("reactance of reactor X= %.3f ohm",X)
diff --git a/1076/CH10/EX10.21/10_21.sce b/1076/CH10/EX10.21/10_21.sce new file mode 100755 index 000000000..78e4afecb --- /dev/null +++ b/1076/CH10/EX10.21/10_21.sce @@ -0,0 +1,43 @@ +clear
+clc
+
+V=6.6
+S=10
+Z=V*V/S
+I=S*1e6/(sqrt(3)*V*1e3)
+
+X1=.15
+X2=.75*X1
+X0=.3*X1
+
+x1=X1/3
+x2=X2/3
+x0=X0/3
+
+E=1
+
+x1=X1/3
+x2=X2/3
+x0=X0/3
+IFa=abs(3*E/(%i * (x1+x2+x0)))
+IFa=round(IFa*1000)/1000
+ifa=IFa * I
+mprintf("\n(a)fault current when all gen neutrals gounded= %.2f A",ifa)
+
+x1=X1/3
+x2=X2/3
+x0=X0
+IFb=abs(3*E/(%i * (x1+x2+x0)))
+IFb=round(IFb*1000)/1000
+ifb=IFb * I
+mprintf("\n(b)fault current when one gen neutral gounded= %.2f A",ifb)
+
+x1=X1/3
+x2=X2/3
+x0=X0
+R0=.3
+r0=round(3*R0/Z*1000)/1000
+IFc=abs(3*E/(complex(r0,(x1+x2+x0))))
+IFc=round(IFc*1000)/1000
+ifc=IFc * I
+mprintf("\n(c)fault current when one neutral gounded thru resistance= %.2f A",ifc)
diff --git a/1076/CH10/EX10.22/10_22.sce b/1076/CH10/EX10.22/10_22.sce new file mode 100755 index 000000000..1c214b13e --- /dev/null +++ b/1076/CH10/EX10.22/10_22.sce @@ -0,0 +1,79 @@ +clear
+clc
+
+x1g=.1
+x2g=.1
+x0g=.05
+
+x1t=.05
+x2t=.05
+x0t=.05
+
+x1l=.4
+x2l=.4
+x0l=.8
+
+x1lm=x1l/2
+x2lm=x2l/2
+x0lm=x0l/2
+
+X1=x1g+x1lm
+X2=x2g+x2lm
+X0=(x0g+x0lm)*(x0lm+x0t)/(x0g+x0lm+x0lm+x0t)
+
+X=X1+X2+X0
+Ia1=round((1/(%i*X))*1000)/1000
+Ia2=round((Ia1)*1000)/1000
+Ia0=round((Ia1)*1000)/1000
+Ia=Ia1+Ia2+Ia0
+IFa=abs(Ia)
+
+Va1=1-(Ia1 * X1 *%i)
+Va2=0-(Ia2 * X2 *%i)
+Va0=0-(Ia0 * X0 *%i)
+
+a=exp(%i * 2 * %pi/3)
+A=[1 1 1
+1 a^2 a
+1 a a^2
+]
+
+va1=[Va0 Va1 Va2]'
+Va=A*va1
+mprintf("\n(a)")
+mprintf("\If= %.3f ang(%.2f)",abs(Ia),270)
+mprintf("\nVa= %.3f ang(%.2f)",abs(Va(1)),atand(imag(Va(1))/real(Va(1))))
+mprintf("\nVb= %.3f ang(%.2f)",abs(Va(2)),atand(imag(Va(2))/real(Va(2)))+180)
+mprintf("\nVc= %.3f ang(%.2f)",abs(Va(3)),atand(imag(Va(3))/real(Va(3)))+180)
+
+
+X1=x1g+x1lm
+X2=x2g+x2lm
+X0=(x0g+x0lm)
+
+X=X1+X2+X0
+Ia1=round((1/(%i*X))*1000)/1000
+Ia2=round((Ia1)*1000)/1000
+Ia0=round((Ia1)*1000)/1000
+Ia=Ia1+Ia2+Ia0
+IFa=abs(Ia)
+
+Va1=1-(Ia1 * X1 *%i)
+Va2=0-(Ia2 * X2 *%i)
+Va0=0-(Ia0 * X0 *%i)
+
+a=exp(%i * 2 * %pi/3)
+A=[1 1 1
+1 a^2 a
+1 a a^2
+]
+
+va1=[Va0 Va1 Va2]'
+Va=A*va1
+mprintf("\n(b)")
+mprintf("\If= %.3f ang(%.2f)",abs(Ia),270)
+mprintf("\nVa= %.3f ang(%.2f)",abs(Va(1)),atand(imag(Va(1))/real(Va(1))))
+mprintf("\nVb= %.3f ang(%.2f)",abs(Va(2)),atand(imag(Va(2))/real(Va(2)))+180)
+mprintf("\nVc= %.3f ang(%.2f)",abs(Va(3)),atand(imag(Va(3))/real(Va(3)))+180)
+
+
diff --git a/1076/CH10/EX10.24/10_24.sce b/1076/CH10/EX10.24/10_24.sce new file mode 100755 index 000000000..073124204 --- /dev/null +++ b/1076/CH10/EX10.24/10_24.sce @@ -0,0 +1,22 @@ +clear
+clc
+
+X1=.4
+X2=.3
+X0=.05
+
+S=15
+V=13.2
+Z=V*V/S
+E=1
+If=1
+
+Xn=((3*E/(If))-(X1+X2+X0))/3
+xn=Xn*Z
+mprintf("\n(a)Xn= %.3f ohm ",xn)
+
+Rn=((3*E/(If))-((X1+X2+X0)*%i))/3
+rn=Rn*Z
+mprintf("\n(b)Rn= %.2f ohm ",rn)
+//the differnece in result is due to error in calculation in textbook
+disp("the differnece in result is due to error in calculation in textbook")
diff --git a/1076/CH10/EX10.27/10_27.sce b/1076/CH10/EX10.27/10_27.sce new file mode 100755 index 000000000..1774837a6 --- /dev/null +++ b/1076/CH10/EX10.27/10_27.sce @@ -0,0 +1,34 @@ +clear
+clc
+
+S=50
+data=[.05 20
+.08 50
+.06 30
+.08 50
+.04 30
+.05 40
+.05 50
+.05 40
+]
+for(i=1:8)
+ X(i)=round(data(i,1) * S/data(i,2)*10000)/10000
+end
+
+X1=round((((X(2)*X(8))+(X(2)*X(7))+(X(7)*X(8)))/X(2))*1000)/1000
+X2=round((((X(2)*X(8))+(X(2)*X(7))+(X(7)*X(8)))/X(7))*1000)/1000
+X3=round((((X(2)*X(8))+(X(2)*X(7))+(X(7)*X(8)))/X(8))*1000)/1000
+X4=round((((1/X(1)) + (1/X2))^-1)*1000)/1000
+X5=round((((1/X(4)) + (1/X3))^-1)*1000)/1000
+X6=round((X4*X5/(X4+X5+X1))*1000)/1000
+X7=round((X4*X1/(X4+X5+X1))*1000)/1000
+X8=round((X1*X5/(X4+X5+X1))*1000)/1000
+X9=round((X7+X(5))*1000)/1000
+X10=round((X8+X(6))*1000)/1000
+X11=round((((1/X10) + (1/X9))^-1)*1000)/1000
+X12=round((X11+X6)*1000)/1000
+X13=((1/X12) + (1/X(3)))^-1
+MVA=S/X13
+//Mismatch is due to error in calculation in the textbook
+mprintf("Fault MVA=%.2f MVA",MVA)
+disp("Mismatch is due to error in calculation in the textbook")
diff --git a/1076/CH10/EX10.28/10_28.sce b/1076/CH10/EX10.28/10_28.sce new file mode 100755 index 000000000..be728e206 --- /dev/null +++ b/1076/CH10/EX10.28/10_28.sce @@ -0,0 +1,38 @@ +clear
+clc
+
+S=30
+SF=1000
+V1=33
+V2=132
+Z1=V1*V1/S
+I2=S*1e6/(sqrt(3)*V2*1e3)
+
+
+x1g=%i*S/SF
+x2g=x1g*2/3
+x0g=x1g/3
+r0g=60/Z1
+z0g=r0g+x0g
+xt=%i *.1
+
+X1=xt+x1g
+X2=xt+x2g
+X0=xt
+
+FMVAa=round(abs(S/X1)*100)/100
+IFa=abs(1/X1)* I2
+mprintf("\n(a)3 phase fault Fault MVA=%.2f MVA, Fault Current=%.2f A",FMVAa,IFa)
+
+IFb=abs(3/(X1+X2+X0))* I2
+FMVAb=IFb * sqrt(3) * V2 * 1e-3
+mprintf("\n(b)single line to ground fault Fault MVA=%.2f MVA, Fault Current=%.2f A",FMVAb,IFb)
+
+IFc=abs(sqrt(3)/(X1+X2))* I2
+FMVAc=IFc * sqrt(3) * V2 * 1e-3
+mprintf("\n(c)L-L Fault MVA=%.2f MVA, Fault Current=%.1f A",FMVAc,IFc)
+
+IFd=abs(1/(X1+(X2*X0/(X2+X0))))* I2
+FMVAd=IFd * sqrt(3) * V2 * 1e-3
+mprintf("\n(d)L-L-G fault Fault MVA=%.2f MVA, Fault Current=%.2f A",FMVAd,IFd)
+
diff --git a/1076/CH10/EX10.29/10_29.sce b/1076/CH10/EX10.29/10_29.sce new file mode 100755 index 000000000..4efcb49c4 --- /dev/null +++ b/1076/CH10/EX10.29/10_29.sce @@ -0,0 +1,128 @@ +clear
+clc
+
+a=exp(%i * 2*%pi/3)
+
+V=33
+S=45
+SF=2000
+V2=132
+Z=V2*V2/S
+I=S*1e6 /(sqrt(3) * V2 * 1e3)
+
+X1=.4* 60
+X2=.4* 60
+X0=1 * 60
+
+x1=X1/Z
+x2=X2/Z
+x0=X0/Z
+
+xg=.0225
+xt=.1
+
+xf1=x1/2 + xg + xt
+xf2=x2/2 + xg + xt
+xf0=(x0/2 + xt)/2
+xf=xf1+xf2+xf0
+ia1=round(1000/xf)/1000
+iF1=3*ia1
+IF1=iF1*I
+mprintf("\n(a)Fault Current = %.2fA",IF1)
+
+IA1=ia1*-1*%i
+IA2=ia1*-1*%i
+IA0= ia1 * .5*-1*%i
+
+IA=round((IA1 +IA2 +IA0)*10000)/10000
+IB=round(((IA1*a*a) + (a*IA2) +IA0)*10000)/10000
+IC=round(((IA1*a) + (a*a*IA2) +IA0)*10000)/10000
+mprintf("\n(b)From T1 to P: (in order - IA,IB, IC in pu)")
+disp(IC, IB, IA)
+
+IA1b=0
+IA2b=0
+IA0b=ia1 * .5*-1*%i
+IAb=round((IA1b +IA2b +IA0b)*10000)/10000
+IBb=round(((IA1b*a*a) + (a*IA2b) +IA0b)*10000)/10000
+ICb=round(((IA1b*a) + (a*a*IA2b) +IA0b)*10000)/10000
+mprintf("\nFrom P to T2: (in order - IA,IB, IC in pu)")
+disp(ICb, IBb, IAb)
+
+
+
+//(c)
+
+
+Ia0=0
+Ia1=IA1 * exp(%i * -1*%pi/6)
+Ia2=IA2 * exp(%i * 1*%pi/6)
+
+Ia=round((Ia1 +Ia2 +Ia0)*1000)/1000
+Ib=round(((Ia1*a*a) + (a*Ia2) +Ia0)*1000)/1000
+Ic=round(((Ia1*a) + (a*a*Ia2) +Ia0)*1000)/1000
+mprintf("\n(c) Currents in lines connecting source to T1 (in order - Ia,Ib, Ic in pu)")
+disp(Ic, Ib, Ia)
+
+
+mprintf("\n(d) Currents in star wdg of T1")
+mprintf("\nIa= %.4fj, Ib= %.4fj, Ic= %.4fj", imag(IA), imag(IB), imag(IC))
+mprintf("\nIa1= %.4fj, Ia2= %.4fj, Ia0= %.4fj", imag(IA1), imag(IA2), imag(IA0))
+mprintf("\nIb1= %.4f /_%.3f, Ib2=%.4f/_%.3f, Ib0= %.4fj", abs(IA1 *a*a),atand(imag(IA1*a*a)/real(IA1*a*a))+180,abs(IA1 *a),atand(imag(IA2*a)/real(IA2*a)),imag(IA0))
+mprintf("\nIc1= %.4f /_%.3f, Ic2=%.4f/_%.3f, Ic0= %.4fj", abs(IA1 *a),atand(imag(IA1*a)/real(IA1*a)),abs(IA1 *a*a),atand(imag(IA2*a*a)/real(IA2*a*a))+180,imag(IA0))
+
+
+Iab1 = round(IA1*1e3 /sqrt(3))/1e3
+Iab2 = round(IA2*1e3 /sqrt(3))/1e3
+Iab0 = round(IA0*1e3 /sqrt(3))/1e3
+
+Ibc1=Iab1 * a*a
+Ibc2=Iab2 * a
+Ibc0=Iab0
+
+Ica1=Iab1 * a
+Ica2=Iab2 * a *a
+Ica0=Iab0
+
+Iab= Iab1 +Iab2 +Iab0
+Ibc= Ibc1 +Ibc2 +Ibc0
+Ica= Ica1 +Ica2 +Ica0
+
+mprintf("\n\nCurrents in delta wdg of T1")
+mprintf("\nIab= %.4fj, Ibc= %.4fj, Ica= %.4fj", imag(Iab), imag(Ibc), imag(Ica))
+mprintf("\nIab1= %.4fj, Iab2= %.4fj, Iab0= %.4fj", imag(Iab1), imag(Iab2), imag(Iab0))
+mprintf("\nIbc1= %.4f /_%.3f, Ibc2=%.4f/_%.3f, Ibc0= %.4fj", abs(Ibc1),atand(imag(Ibc1)/real(Ibc1))+180,abs(Ibc2),atand(imag(Ibc2)/real(Ibc2)),imag(Ibc0))
+mprintf("\nIca1= %.4f /_%.3f, Ica2=%.4f/_%.3f, Ica0= %.4fj", abs(Ica1),atand(imag(Ica1)/real(Ica1)),abs(Ica2),atand(imag(Ica2)/real(Ica2))+180,imag(Ica0))
+
+mprintf("\n\n Currents in star wdg of T2")
+mprintf("\nIa= %.4fj, Ib= %.4fj, Ic= %.4fj", imag(IAb), imag(IBb), imag(ICb))
+mprintf("\nIa1= %.4f, Ia2= %.4f, Ia0= %.4fj", imag(IA1b), imag(IA2b), imag(IA0b))
+mprintf("\nIb1= %.4f, Ib2= %.4f, Ib0= %.4fj", imag(IA1b *a*a), imag(IA1b *a),imag(IA0b))
+mprintf("\nIc1= %.4f, Ic2= %.4f, Ic0= %.4fj", imag(IA1b *a), imag(IA1b*a *a),imag(IA0b))
+
+
+Iab1b = round(IA1b*1e3 /sqrt(3))/1e3
+Iab2b = round(IA2b*1e3 /sqrt(3))/1e3
+Iab0b = round(IA0b*1e3 /sqrt(3))/1e3
+
+Ibc1b=Iab1b * a*a
+Ibc2b=Iab2b * a
+Ibc0b=Iab0b
+
+Ica1b=Iab1b * a
+Ica2b=Iab2b * a *a
+Ica0b=Iab0b
+
+Iabb= Iab1b +Iab2b +Iab0b
+Ibcb= Ibc1b +Ibc2b +Ibc0b
+Icab= Ica1b +Ica2b +Ica0b
+
+mprintf("\n\nCurrents in delta wdg of T2")
+mprintf("\nIab= %.3fj, Ibc= %.3fj, Ica= %.3fj", imag(Iabb), imag(Ibcb), imag(Icab))
+mprintf("\nIab1= %.3f, Iab2=%.3f, Iab0= %.3fj", imag(Iab1b), imag(Iab2b), imag(Iab0b))
+mprintf("\nIbc1= %.3f, Ibc2=%.3f, Ibc0= %.3fj", imag(Ibc1b),imag(Ibc2b),imag(Ibc0b))
+mprintf("\nIca1= %.3f, Ica2=%.3f, Ica0= %.3fj", imag(Ica1b),imag(Ica2b),imag(Ica0b))
+
+
+
+
diff --git a/1076/CH10/EX10.3/10_3.sce b/1076/CH10/EX10.3/10_3.sce new file mode 100755 index 000000000..9d7d4d878 --- /dev/null +++ b/1076/CH10/EX10.3/10_3.sce @@ -0,0 +1,28 @@ +clear
+clc
+
+G=[100 11 .2]
+T=[100 11 132 .05]
+L=[
+1.2e-3 100
+1e-3 50
+1e-3 50
+]
+B=[100 11]
+T(5)= T(3)/T(2)
+B(3)=B(2)* T(5)
+B(4)=B(3)^2/B(1)
+for(i=1:3)
+ L(i,3)= 2*%pi*50* L(i,1)* L(i,2);
+ L(i,4)= L(i,3)/B(4)
+end
+l=L(:,4)
+X1=l(1)*l(2)/(l(1)+l(2)+l(3))
+X2=l(3)*l(2)/(l(1)+l(2)+l(3))
+X3=l(1)*l(3)/(l(1)+l(2)+l(3))
+X=(((G(3)+T(4)+(X1))^-1) + ((G(3)+T(4)+(X3))^-1))^-1 +X2
+SF=round(10/X)/10
+S=SF* B(1)
+IF=round(10/X)/10
+I=IF*B(1)*1e6/(sqrt(3)*B(3)*1e3)
+mprintf("Fault level = %.0f MVA, Fault current %.1f A", S,I)
diff --git a/1076/CH10/EX10.30/10_30.sce b/1076/CH10/EX10.30/10_30.sce new file mode 100755 index 000000000..37f5dcf36 --- /dev/null +++ b/1076/CH10/EX10.30/10_30.sce @@ -0,0 +1,57 @@ +clear
+clc
+
+R1=4
+R2=2
+
+S=50
+V1=11
+V2=132
+V3=33
+Z1=V1*V1/S
+Z2=V2*V2/S
+Z3=V3*V3/S
+
+r1=3*R1/Z1
+r2=3*R2/Z3
+
+
+x1g=.4*%i
+x2g=.3*%i
+x0g=.1*%i
+
+x1t1=.08*%i
+x2t1=.08*%i
+x0t1=.08*%i
+
+x1t2=.05*%i
+x2t2=.05*%i
+x0t2=.05*%i
+
+x1t3=.04*%i
+x2t3=.04*%i
+x0t3=.04*%i
+
+x1t4=.06*%i
+x2t4=.06*%i
+x0t4=.06*%i
+
+
+X1l=20*%i
+X2l=20*%i
+X0l=50*%i
+
+x1l=X1l/Z2
+x2l=X2l/Z2
+x0l=X0l/Z2
+
+X1=x1g+x1t2+x1l+x1t1+x1t3
+X2=x2g+x2t2+x2l+x2t1+x2t3
+X0=r2+(((x0t2+x0l+x0t1)*x0t4/(x0t2+x0l+x0t1+x0t4))+x0t3)
+
+IF=abs(3*1/(X1+X2+X0))
+IB=S*1e6/(sqrt(3)*V3*1e3)
+If=IF*IB
+SF=IF*S
+mprintf("fault current= %.0fA, fault level=%.2f MVA",If,SF)
+
diff --git a/1076/CH10/EX10.31/10_31.sce b/1076/CH10/EX10.31/10_31.sce new file mode 100755 index 000000000..a1d27ac14 --- /dev/null +++ b/1076/CH10/EX10.31/10_31.sce @@ -0,0 +1,37 @@ +clear
+clc
+
+a=exp(%i *2*%pi/3)
+
+Z1=complex(2.8,1)
+Z2=complex(.1,.6)
+
+V=400
+E=V/sqrt(3)
+
+Ia1=E/(Z1+Z2)
+Ia2=-Ia1
+
+Ia=Ia1+Ia2
+Ib= (a^2-a)*Ia1
+Ic=-Ib
+
+disp( "Line Currents Ia, Ib, Ic, in amperes")
+mprintf("\nIa= %s", string(round(abs(Ia)*10)/10) +'/_'+ string(0) )
+mprintf("\nIb= %s", string(round(abs(Ib)*10)/10) +'/_'+ string(round(atand(imag(Ib)/real(Ib))*100)/100 -180) )
+mprintf("\nIc= %s", string(round(abs(Ic)*10)/10) +'/_'+ string(round(atand(imag(Ic)/real(Ic))*100)/100) )
+
+Va2=-Z2 * Ia2
+Vaa=3*Va2
+Van=(Z1*Ia1)+(Z2*Ia2)
+Vcn=(a*Z1*Ia1)+(a*a*Z2*Ia2)
+Vbn=(a*a*Z1*Ia1)+(a*Z2*Ia2)
+VNn=Va2
+
+
+mprintf("\n\n\nVaa= %s", string(round(abs(Vaa)*100)/100) +'/_'+ string(round(atand(imag(Vaa)/real(Vaa))*10)/10))
+mprintf("\nVan= %s", string(round(abs(Van)*100)/100) +'/_'+ string(round(atand(imag(Van)/real(Van))*10)/10))//error in value substitution in textbook
+mprintf("\nVbn= %s", string(round(abs(Vbn)*10)/10) +'/_'+ string(round(atand(imag(Vbn)/real(Vbn))*10)/10 -180))
+mprintf("\nVcn= %s", string(round(abs(Vcn)*10)/10) +'/_'+ string(round(atand(imag(Vcn)/real(Vcn))*10)/10 +180)) //error in value substitution in textbook
+mprintf("\nVNn= %s", string(round(abs(VNn)*100)/100) +'/_'+ string(round(atand(imag(VNn)/real(VNn))*10)/10))
+disp("error is due to mistake in value substitution in textbook")
diff --git a/1076/CH10/EX10.32/10_32.sce b/1076/CH10/EX10.32/10_32.sce new file mode 100755 index 000000000..1627a977c --- /dev/null +++ b/1076/CH10/EX10.32/10_32.sce @@ -0,0 +1,16 @@ +clear
+clc
+
+S=10
+xg=.1
+xe=.08
+
+X1= 1/((1/.1) + 1/(xe + ((xg+xe)/2)))
+FMVA1=S* (1/X1)
+mprintf("When reactors are used, fault level=%.2f MVA", FMVA1)
+
+X2= xg/3
+FMVA2=S* (1/X2)
+mprintf("\nWhen reactors are not used, fault level=%.0f MVA",FMVA2)
+
+
diff --git a/1076/CH10/EX10.33/10_33.sce b/1076/CH10/EX10.33/10_33.sce new file mode 100755 index 000000000..3e54418b9 --- /dev/null +++ b/1076/CH10/EX10.33/10_33.sce @@ -0,0 +1,32 @@ +clear
+clc
+
+S=25
+pf=.8
+P=15
+Vt=10.6
+V1=11
+V2=11 * 66/11
+I1=S*1e6/(sqrt(3)*V1*1e3)
+I2=S*1e6/(sqrt(3)*V2*1e3)
+Il=(P*1e6/(sqrt(3)*Vt*1e3*pf)) * exp (%i * acos(pf))
+vt=Vt/V1
+Z=V2*V2/S
+XL=10
+xl=XL/Z
+xt=.1
+xg=.15
+xm=.15
+
+Xth=xm * (xg+xt+xt+xl)/(xm + xg+xt+xt+xl)
+IF=vt/Xth
+If=IF*I2
+
+Ifg=IF*I1 *%i*-1* xm /(xm + xg+xt+xt+xl)
+ifg=abs(Ifg + Il)
+Ifm=IF*I1 *%i *-1* (xg+xt+xt+xl)/(xg+xt+xt+xl+xm)
+ifm=abs(Ifm - Il)
+
+mprintf("total fault current = %.0f A, current through generator=%.0f A, current through motor=%.0f A\n",If, ifg, ifm)
+//error in calculation of Ifm-I =-.623 - 6.891j instead of -.623-j5.96
+disp("error in calculation of Ifm-I =-.623 - 6.891j->(correct) instead of -.623-j5.96 -> incorrect")
diff --git a/1076/CH10/EX10.34/10_34.sce b/1076/CH10/EX10.34/10_34.sce new file mode 100755 index 000000000..28f6eac88 --- /dev/null +++ b/1076/CH10/EX10.34/10_34.sce @@ -0,0 +1,13 @@ +clear
+clc
+
+S=25
+V=11
+Z=V*V/S
+I=S*1e6/(sqrt(3)*V*1e3)
+Isc=6*I
+Xt=V*1e3/(sqrt(3)*Isc)
+Xi=.15*Z
+Xo=Xt-Xi
+x=Xo*100/Z
+mprintf("External reactance required is %.3f pu",x)
diff --git a/1076/CH10/EX10.35/10_35.sce b/1076/CH10/EX10.35/10_35.sce new file mode 100755 index 000000000..290403145 --- /dev/null +++ b/1076/CH10/EX10.35/10_35.sce @@ -0,0 +1,28 @@ +clear
+clc
+
+Ia=10*exp(%i *30 *%pi/180)
+Ib=15*exp(%i *-60*%pi/180)
+Ic=0-(Ia+Ib)
+
+Iac=(Ia +Ia +Ib)/3
+Icb=Ic +Iac
+Iba=Iac-Ia
+Ia0=(Ia +Ib +Ic)/3
+
+a=exp(%i * 2 * %pi/3)
+
+Ia1=(Ia + a*Ib + a*a*Ic)/3
+Ia2=(Ia + a*Ic + a*a*Ib)/3
+
+disp(round(Ia1*100)/100, "Ia1","(a)Symmetrical Line Components")
+disp(round(Ia2*100)/100, "Ia2")
+disp(round(Ia0*100)/100, "Ia0")
+
+Iac0=(Iac+Icb+Iba)/3
+Iac1=(Iac+ a*Icb+ a*a*Iba)/3
+Iac2=(Iac+a*a*Icb+a*Iba)/3
+
+disp(round(Iac1*100)/100, "Iac1","(b)Symmetrical Delta Components")
+disp(round(Iac2*100)/100, "Iac2")
+disp(round(Iac0*100)/100, "Iac0")
diff --git a/1076/CH10/EX10.4/10_4.sce b/1076/CH10/EX10.4/10_4.sce new file mode 100755 index 000000000..cdf01d30e --- /dev/null +++ b/1076/CH10/EX10.4/10_4.sce @@ -0,0 +1,32 @@ +clear
+clc
+
+T=[10 132 6.6 .15]
+M=[5 6.6 .3 .2 ]
+B=[10 6.6]
+T(5)= T(3)/T(2)
+B(3)=B(2)* T(5)
+B(4)= B(1)*1e6/(sqrt(3)*B(2)*1e3)
+M(5)=M(4) *B(1)/M(1)
+M(6)=M(3) *B(1)/M(1)
+
+X1=1/((1/M(5))+(1/M(5))+(1/T(4)))
+IF1=round(100/X1)/100
+I1=IF1*B(4)
+mprintf("\n(a) sub transient fault current=%.0f A", I1)
+
+It=round(100/T(4))/100
+Im=1/M(5)
+ID=It+Im
+iD=ID*B(4)
+mprintf("\n(b) current through D=%.0f A", iD)
+
+RD=iD*1.6
+mprintf("\n(c) current rating of D=%.0f A", RD)
+
+X2=1/((1/M(6))+(1/T(4)))
+IF2=round(100/X2)/100
+I2=IF2*round(B(4)*10)/10
+iCB=1.1 *I2
+mprintf("\n(d) current to be interrrupted by D=%.1f A", fix(iCB*10)/10)
+
diff --git a/1076/CH10/EX10.5/10_5.sce b/1076/CH10/EX10.5/10_5.sce new file mode 100755 index 000000000..d5a93953a --- /dev/null +++ b/1076/CH10/EX10.5/10_5.sce @@ -0,0 +1,33 @@ +clear
+clc
+
+G=[100 11 .25]
+M=[50 11 .2 40 .8]
+xl=.05
+vt=10.95
+
+B=[100 11]
+B(3)= B(1)*1e6/(sqrt(3)*B(2)*1e3)
+
+IL=M(4)*1e6 / (sqrt(3)*vt*1e3*M(5))
+Il=round(((IL/B(3) * exp(%i * acos(M(5)))))*1000)/1000
+
+Vt=round(1e3*vt/B(2))/1e3
+
+V=fix((Vt + Il*xl*%i)*1e4)*1e-4
+
+M(6)=M(3) *B(1)/M(1)
+
+xth= round( G(3) * (M(6) +xl)/(G(3) + (M(6) +xl)) *10000)/10000
+
+If=fix(V*1e3/(%i *xth))/1e3
+temp= fix(imag(If)*100)/100
+If=complex(real(If), temp)
+Ifg=fix((If * (M(6) +xl)/(G(3) + (M(6) +xl)))*100)/100
+Ifm=round((If * (G(3) )/(G(3) + (M(6) +xl)))*1000)/1000
+
+Ig=(Ifg + Il) * B(3)
+Im=(Ifm - Il) *B(3)
+
+mprintf("\ntotal generator current during fault= %s A", string(round(abs(Ig)*10)/10) +'/_'+ string(round(atand(imag(Ig)/real(Ig))*10)/10))
+mprintf("\ntotal motor current during fault= %s A", string(fix(abs(Im)*1)/1) +'/_'+ string(fix(atand(imag(Im)/real(Im))*10)/10 +180))
diff --git a/1076/CH10/EX10.6/10_6.sce b/1076/CH10/EX10.6/10_6.sce new file mode 100755 index 000000000..e69a143cc --- /dev/null +++ b/1076/CH10/EX10.6/10_6.sce @@ -0,0 +1,19 @@ +clear
+clc
+
+I=[
+5*exp(%i * %pi * 60 / 180)
+5*exp(%i * %pi * -60 / 180)
+0
+]
+a=exp(%i * 2 * %pi/3)
+A=[1 1 1
+1 a^2 a
+1 a a^2
+]
+
+Is=inv(A)*I
+
+mprintf("\nIa0= %s", string(round(abs(Is(1))*1000)/1000) +'/_'+ string(round(atand(imag(Is(1))/real(Is(1)))*100)/100) )
+mprintf("\nIa1= %s", string(round(abs(Is(2))*1000)/1000) +'/_'+ string(round(atand(imag(Is(2))/real(Is(2)))*100)/100) )
+mprintf("\nIa2= %s", string(round(abs(Is(3))*1000)/1000) +'/_'+ string(round(atand(imag(Is(3))/real(Is(3)))*100)/100 +180) )
diff --git a/1076/CH10/EX10.8/10_8.sce b/1076/CH10/EX10.8/10_8.sce new file mode 100755 index 000000000..b35e7a508 --- /dev/null +++ b/1076/CH10/EX10.8/10_8.sce @@ -0,0 +1,27 @@ +clear
+clc
+
+G=[50 11 0 .08
+30 11 0 .07]
+T=[50 11 220 .1
+30 220 11 .09]
+
+B=[50 11]
+
+T(1,5)= T(1,3)/T(1,2)
+T(2,5)= T(2,3)/T(2,2)
+B(3)=B(2)* T(1,5)
+B(4)=B(3)* T(2,5)
+
+B(5)= B(3)^2/B(1)
+
+Z=555.6
+z=Z/B(5)
+
+zt2=T(2,4) * B(1)/T(2,1)
+zg2=G(2,4) * B(1)/G(2,1)
+
+Zn=3
+zn=Zn *3 / ( B(4)^2/B(1))
+
+mprintf("zero seq netwk: xt1= %.1f, xt2=%.2f, Xg1=%.2f Xg2=%.3f, xl=%.3f, Zn=%.2fi",T(1,4), zt2, G(1,4), zg2, z,zn)
|