diff options
Diffstat (limited to '1076/CH16')
28 files changed, 890 insertions, 0 deletions
diff --git a/1076/CH16/EX16.1/16_1.sce b/1076/CH16/EX16.1/16_1.sce new file mode 100755 index 000000000..199a2ade1 --- /dev/null +++ b/1076/CH16/EX16.1/16_1.sce @@ -0,0 +1,28 @@ +clear
+clc
+
+I1=100
+I2=150
+I3=200
+
+l1=150
+l2=100
+l3=100
+
+r=.1/1000
+
+Va=200
+
+rac=l1*r;
+rcd=l2*r
+rbd=l3*r
+
+dvc=(I1+I2+I3)*rac;
+dvd=(I1+I2+I3)*rac + (I2+I3)*rcd ;
+dvb=(I1+I2+I3)*rac + (I2+I3)*rcd + (I3*rbd);
+
+Vc=Va-dvc
+Vd=Va-dvd
+Vb=Va-dvb
+
+mprintf("voltage at B= %.2f V C= %.2f V D= %.2f V ", Vb, Vc, Vd)
diff --git a/1076/CH16/EX16.10/16_10.sce b/1076/CH16/EX16.10/16_10.sce new file mode 100755 index 000000000..07bd6e9f3 --- /dev/null +++ b/1076/CH16/EX16.10/16_10.sce @@ -0,0 +1,52 @@ +clear
+clc
+
+I=200
+
+r1=.05
+r2=.06
+r3=.02
+r4=.04
+r5=.03
+r6=.01
+ra=.02
+rb=.03
+
+I1=100
+I3=30
+I4=50
+I5=20
+//(a)
+dv=0
+A=[
+(ra) (-rb) (r6+r5 +r4+r3)
+1 1 0
+(ra+r1+r2) (-rb) -(r1+r2)
+]
+
+B=[
+dv+ (r5*(I5) +r4*(I5+I4)+r3*(I5+I4+I3))
+I
+dv+(r2*I1)
+]
+i=inv(A)*B
+
+mprintf("\n(a) Ia = %.0f A, Ib= %.0f A", i(1), i(2))
+
+//(b)
+dv=-5
+A=[
+(ra) (-rb) (r6+r5 +r4+r3)
+1 1 0
+(ra+r1+r2) (-rb) -(r1+r2)
+]
+
+B=[
+dv+(r5*(I5) +r4*(I5+I4)+r3*(I5+I4+I3))
+I
+dv+(r2*I1)
+]
+i=inv(A)*B
+
+mprintf("\n(b) Ia = %d A, Ib= %d A", i(1), i(2))
+
diff --git a/1076/CH16/EX16.11/16_11.sce b/1076/CH16/EX16.11/16_11.sce new file mode 100755 index 000000000..cebc32df5 --- /dev/null +++ b/1076/CH16/EX16.11/16_11.sce @@ -0,0 +1,13 @@ +clear
+clc
+
+L=1200
+L1=900
+L2=600
+
+r=1.5
+
+x=(L1*L + L*L*r/2)/(L1 + L2 + (r*2*L/2))
+y=L-x
+mprintf("\nCurrent in CB= %.0f a A", x)
+mprintf("\nCurrent in CA= %.0f a A", y)
diff --git a/1076/CH16/EX16.12/16_12.sce b/1076/CH16/EX16.12/16_12.sce new file mode 100755 index 000000000..8d744c5c3 --- /dev/null +++ b/1076/CH16/EX16.12/16_12.sce @@ -0,0 +1,20 @@ +clear
+clc
+
+V=250
+
+L=[
+150 200 280 320 390 450 500
+20 40 35 25 10 20 30
+]
+
+r=.1/(2*500)
+
+D=L(1,:)
+I=L(2,:)'
+
+dv=2*r*D*I
+
+Ve=V-dv
+
+mprintf("Voltage at far end is %.2f V", Ve)
diff --git a/1076/CH16/EX16.13/16_13.sce b/1076/CH16/EX16.13/16_13.sce new file mode 100755 index 000000000..d1a09e4f5 --- /dev/null +++ b/1076/CH16/EX16.13/16_13.sce @@ -0,0 +1,24 @@ +clear
+clc
+
+V=250
+
+L=[
+150 50 80 40 70 60 50 150
+0 20 60 95 120 130 150 180
+]
+
+r=.1/(2*500)
+
+D=L(1,:)
+I1=L(2,:)'
+I2=ones(8,1)
+
+dv1=2*r*D*I1
+dv2=2*r*D*I2
+
+Ia=dv1/dv2
+Ib=L(length(L))-Ia
+
+Vc=V-(2*r*((Ia*D(1))+((Ia-I1(2))*D(2))+((Ia-I1(3))*D(3))))
+mprintf("Ia= %.2f A, Ib= %.2f A, Vmin at C = %.3f V",Ia, Ib, Vc)
diff --git a/1076/CH16/EX16.14/16_14.sce b/1076/CH16/EX16.14/16_14.sce new file mode 100755 index 000000000..19149594d --- /dev/null +++ b/1076/CH16/EX16.14/16_14.sce @@ -0,0 +1,35 @@ +clear
+clc
+
+
+r1=.03
+r2=.02
+r3=.03
+r4=.04
+r5=.04
+r6=.01
+r7=.02
+
+I1=20
+I2=30
+I3=25
+I4=30
+I5=125
+I6=20
+
+dv=0
+A=[
+(r5+r1+r4+r3+r2) (r3+r2)
+(r3+r2) (r3+r2+r6+r7)
+]
+
+
+B=[
+((r4*I5)+(r3*(I5-I4))+ (r2*(I5-I4-I3)) + (r1*I1))
+((r3*(I5-I4))+ (r2*(I5-I4-I3))+(r7*I6))
+]
+i=inv(A)*B
+x=i(1)
+y=i(2)
+
+mprintf("\nEA = %.2f A, AB= %.2f A, ED= %.2f A, DF= %.2f A, DC= %.2f A, BC= %.2f A, FB= %.2f A", x,x-I1, I5-x,y, I5-I4-x-y, I5-I4-I3-x-y, y-I6)
diff --git a/1076/CH16/EX16.15/16_15.sce b/1076/CH16/EX16.15/16_15.sce new file mode 100755 index 000000000..ec1c8c624 --- /dev/null +++ b/1076/CH16/EX16.15/16_15.sce @@ -0,0 +1,20 @@ +clear
+clc
+
+z1=complex(2,1)
+z2=complex(2,3)
+z3=complex(1,2)
+
+ib= 40 * exp (%i * -1 * acos (.8))
+ic= 60 * exp (%i * -1 * acos (.6))
+
+i1=((ib*z2)+((ib+ic)*z3))/(z1+z2+z3)
+i2=i1-ib
+i3=i2-ic
+
+mprintf("\n Current A to B in Ampere")
+disp(round(i1*100)/100)
+mprintf("\n Current B to C in Ampere")
+disp(round(i2*100)/100)
+mprintf("\n Current A to C in Ampere")
+disp(round(-i3*100)/100)
diff --git a/1076/CH16/EX16.16/16_16.sce b/1076/CH16/EX16.16/16_16.sce new file mode 100755 index 000000000..09440c0b8 --- /dev/null +++ b/1076/CH16/EX16.16/16_16.sce @@ -0,0 +1,42 @@ +clear
+clc
+
+//(a)
+V=230
+df=5/100
+
+I1=20
+I2=10
+L1=300
+L2=200
+L=500
+I=I1+I2
+
+rho=.0286
+T=20
+dT=30
+a=.004
+rho1=rho*(1+(dT*a))
+rho1=round(rho1*1000)/1000
+
+A=1
+A1=I1/I
+A2=I2/I
+
+A=rho1*2*((L*round(I/A))+(L1*round(I1/A1)))/(df * V)
+A=round(A*100)/100
+
+A1=A*A1
+A2=A*A2
+
+mprintf("\n(a)CROSSECTIONS: SC= %.2f mm sq, CA= %.2f mm sq, CB= %.2f mm sq",A, A1, A2 )
+
+R=2*rho1 * (L/A)
+R1=2*rho1 * (L1/A1)
+R2=2*rho1 * (L2/A2)
+
+P=(I*I*R)+(I1*I1*R1)+(I2*I2*R2)
+
+PL=P*100/((V*I1)+(V*I2))
+mprintf("\n(b)Percentage power loss = %.1f percent",PL )
+
diff --git a/1076/CH16/EX16.17/16_17.sce b/1076/CH16/EX16.17/16_17.sce new file mode 100755 index 000000000..c7a005d27 --- /dev/null +++ b/1076/CH16/EX16.17/16_17.sce @@ -0,0 +1,49 @@ +clear
+clc
+
+I1=100
+I2=50
+I3=50
+I4=100
+I5=0
+
+L=500
+l1=50
+l2=100
+l3=100
+l4=150
+l5=L-l1-l2-l3-l4
+
+R=.075
+r=R/L
+U=1
+
+Va=230
+
+r1=l1*r;
+r2=l2*r
+r3=l3*r
+r4=l4*r
+r5=l5*r
+
+Ia= I1+I2+I3+I4+I5 + U*(l1+l2+l3+l4+l5)
+dvab=(Ia - (.5* U * l1))*r1;
+Vb=Va-dvab
+
+Ib= Ia- (U*(l1)) - I1
+dvbc=(Ib - (.5* U * l2))*r2;
+Vc=Vb-dvbc
+
+Ic= Ib- (U*(l2)) - I2
+dvcd=(Ic - (.5* U * l3))*r3;
+Vd=Vc-dvcd
+
+Id= Ic- (U*(+l3)) - I3
+dvde=(Id - (.5* U * l4))*r4;
+Ve=Vd-dvde
+
+Ie= Id- (U*(l4)) - I4
+dvef=(Ie - (.5* U * l5))*r5;
+Vf=Ve-dvef
+
+mprintf("voltage at B= %.4f V C= %.4f V D= %.4f V E= %.2f V F= %.1f V ", Vb, Vc, Vd,Ve, Vf)
diff --git a/1076/CH16/EX16.18/16_18.sce b/1076/CH16/EX16.18/16_18.sce new file mode 100755 index 000000000..2aea77dfa --- /dev/null +++ b/1076/CH16/EX16.18/16_18.sce @@ -0,0 +1,25 @@ +clear
+clc
+
+Vc=220
+
+Z=2*complex(.15, .2)
+Z1=.5*Z
+Z2=.5*Z
+
+Ib= 60 * exp (%i * -acos (.9))
+Ic= 100 * exp (%i * -acos (.8))
+
+Vb=Vc + (Ic * Z2)
+Iab=Ib + Ic
+
+Va=Vb + (Iab * Z1)
+Vs=abs(Va)
+
+Is= abs(Iab)
+pf=cosd(abs(atand(imag(Va)/real(Va))+abs(atand(imag(Iab)/real(Iab)))))
+
+mprintf("\nsending end voltage = %.1f V", Vs)
+mprintf("\nsending end current = %.2f A", Is)
+mprintf("\nsending end pf = %.3f", pf)
+
diff --git a/1076/CH16/EX16.19/16_19.sce b/1076/CH16/EX16.19/16_19.sce new file mode 100755 index 000000000..af0c830d0 --- /dev/null +++ b/1076/CH16/EX16.19/16_19.sce @@ -0,0 +1,19 @@ +clear
+clc
+
+ab=1000
+bc=900
+ac=600
+
+A=500
+B=750
+C=450
+
+theta=acos((ab^2 + ac^2 - bc^2)/(2*ac*ab))
+am=ac * cos(theta)
+cm=ac * sin(theta)
+
+x= ((B*ab)+(C*am))/(A+B+C)
+y= ((C*cm))/(A+B+C)
+
+mprintf("Coordinates wrt A are (%.2f, %.0f)", x,y)
diff --git a/1076/CH16/EX16.2/16_2.sce b/1076/CH16/EX16.2/16_2.sce new file mode 100755 index 000000000..b3c725194 --- /dev/null +++ b/1076/CH16/EX16.2/16_2.sce @@ -0,0 +1,33 @@ +clear
+clc
+
+I1=100
+I2=150
+I3=200
+
+l1=150
+l2=100
+l3=100
+
+r=.1/1000
+U=1
+
+Va=200
+
+rac=l1*r;
+rcd=l2*r
+rbd=l3*r
+
+Iac= I1+I2+I3 + U*(l1+l2+l3)
+dvc=(Iac - (.5* U * l1))*rac;
+Vc=Va-dvc
+
+Icd= I2+I3 + U*(l2+l3)
+dvd=(Icd - (.5* U * l2))*rcd ;
+Vd=Vc-dvd
+
+Idb= I3 + U*(l3)
+dvb=(Idb - (.5* U * l3))*rbd ;
+Vb=Vd-dvb
+
+mprintf("voltage at B= %.3f V C= %.3f V D= %.3f V ", Vb, Vc, Vd)
diff --git a/1076/CH16/EX16.20/16_20.sce b/1076/CH16/EX16.20/16_20.sce new file mode 100755 index 000000000..959b47f55 --- /dev/null +++ b/1076/CH16/EX16.20/16_20.sce @@ -0,0 +1,35 @@ +clear
+clc
+
+rho=1/58
+
+l1=90
+l2=90
+l3=100
+l4=80
+l5=90
+l6=90
+l7=90
+
+I1=50
+I2=40
+I3=20
+I4=30
+I5=45
+I6=50
+
+I=I1+I2+I3+I4+I5+I6
+L=l1+l2+l3+l4+l5+l6+l7
+Ia=((I1*l2)+((I1+I2)*l3)+((I1+I2+I3)*l4)+((I1+I2+I3+I4)*l5)+((I1+I2+I3+I4+I5)*l6)+((I1+I2+I3+I4+I5+I6)*l7))/L
+
+Ib=I-Ia
+V=230*2
+V1=440
+
+dv=V-V1
+
+r=dv/((Ia*(l1+l2+l3+l4))-((I1*l2)+((I1+I2)*l3)+((I1+I2+I3)*l4)))
+
+a=rho*2/r
+
+mprintf("crossection area (in mm sq)= %.2f", a)
diff --git a/1076/CH16/EX16.21/16_21.sce b/1076/CH16/EX16.21/16_21.sce new file mode 100755 index 000000000..8cba6a37c --- /dev/null +++ b/1076/CH16/EX16.21/16_21.sce @@ -0,0 +1,81 @@ +clear
+clc
+
+r1=.08
+r2=.1
+r3=.12
+r4=.14
+r5=.09
+r6=.16
+
+I1=20
+I2=50
+I3=25
+I4=40
+I5=30
+
+x=((r2*I1)+(r3*(I1+I2))+(r4*(I1+I2+I3))+(r5*(I1+I2+I3+I4))+(r6*(I1+I2+I3+I4+I5)))/(r1+r2+r3+r4+r5+r6)
+
+iab=x
+ibc=iab-I1
+icd=ibc-I2
+ide=icd-I3
+ief=ide-I4
+ifa=ief-I5
+
+iab=round(iab*10)/10
+ibc=round(ibc*10)/10
+icd=round(icd*10)/10
+ide=round(ide*10)/10
+ief=round(ief*10)/10
+ifa=round(ifa*10)/10
+
+
+mprintf("\n(a) Current A to B in Ampere = %.1f", iab)
+mprintf("\n Current B to C in Ampere = %.1f", ibc)
+mprintf("\n Current C to D in Ampere = %.1f", icd)
+mprintf("\n Current D to E in Ampere = %.1f", ide)
+mprintf("\n Current E to F in Ampere = %.1f", ief)
+mprintf("\n Current F to A in Ampere = %.1f", ifa)
+
+
+
+Vb=230+((r6*ifa)+(r5*ief)+(r4*ide))
+mprintf("\n Minimum voltage at B, in Volts = %.3f\n\n", Vb)
+
+//(b)
+r7=.1
+
+
+B=[
+(r4+r5+r6) -r7
+(r1+r2+r3) (r7+r1+r3+r2)
+]
+
+C=[
+((r5*I5)+(r4*(I5+I4)))
+((r2*(I5+I2+I3+I4))+(r1*(I1+I2+I3+I4+I5))+(r3*(I1+I2+I3)))
+]
+
+A= inv(B)*C
+x=A(1)
+y=A(2)
+iab=I1+I2+I3+I4+I5 - x-y
+ibc=iab-I1
+icd=ibc-I2
+iaf=x
+ife=iaf-I5
+ied=ife-I4
+iad=y
+
+mprintf("\n(b) Current A to B in Ampere = %.1f", iab)
+mprintf("\n Current B to C in Ampere = %.1f", ibc)
+mprintf("\n Current C to D in Ampere = %.1f", icd)
+mprintf("\n Current E to D in Ampere = %.1f", ied)
+mprintf("\n Current F to E in Ampere = %.1f", ife)
+mprintf("\n Current A to F in Ampere = %.1f", iaf)
+mprintf("\n Current A to D in Ampere = %.1f", iad)
+
+Ve=230-((r6*iaf)+(r5*ife))
+mprintf("\n Minimum voltage at E, in Volts = %.3f", Ve)
+
diff --git a/1076/CH16/EX16.22/16_22.sce b/1076/CH16/EX16.22/16_22.sce new file mode 100755 index 000000000..06861754c --- /dev/null +++ b/1076/CH16/EX16.22/16_22.sce @@ -0,0 +1,24 @@ +clear
+clc
+
+V=33e3
+L=10
+Ce=.8
+data=[
+4e6 .8 10
+2e6 .8 6
+1e6 .8 8
+]
+rho=.0286*1e3
+P=2000
+q=0
+for i=1:3
+ I(i,1)=data (i,1)/(sqrt(3) * V * data(i,2))
+ E(i,1)= I(i)^2 * rho*3* data (i,3) * 1e-3
+ q=q+E(i)
+end
+Q=q*365*Ce
+
+a=sqrt(Q/P)
+
+mprintf("a=%.2f sq mm", fix(a*100)/100)
diff --git a/1076/CH16/EX16.23/16_23.sce b/1076/CH16/EX16.23/16_23.sce new file mode 100755 index 000000000..48ee89056 --- /dev/null +++ b/1076/CH16/EX16.23/16_23.sce @@ -0,0 +1,56 @@ +clear
+clc
+
+V=220
+
+rho1=.0286
+T=20
+dT=30
+a=.004
+rho=rho1*(1+(dT*a))
+rho=round(rho*1000)/1000
+
+P1=5
+ef1=.81
+pf1=.77
+I1= P1*735.5/(V*pf1*ef1)
+I1=I1 * exp(%i *-1 * acos (pf1))
+I1=round(I1 *1000)/1000
+Ic=I1
+
+
+P2=26
+ef2=.87
+pf2=.85
+I2= P2*735.5/(V*pf2*ef2)
+I2=round(I2*100)/100
+I2=I2 * exp(%i *-1 * acos (pf2))
+Ib=fix((I1+I2)*100)/100
+
+P3=10
+ef3=.83
+pf3=.82
+I3= P3*735.5/(V*pf3*ef3)
+I3=I3 * exp(%i *-1 * acos (pf3))
+Ia=I1+I2+I3
+
+P=real (Ia * V )
+PLa=round(.05*P)
+
+l1=40
+l2=35
+l3=30
+
+PL1=2 * (round(abs(Ia)*100)/100)^2 * rho * l1
+PL2=2 * (fix(abs(Ib)*100)/100)^2 * rho * l2
+PL3=2 * (round(abs(Ic)*100)/100)^2 * rho * l3
+
+PL1=fix(PL1/100)*100
+PL2=round(PL2/100)*100
+PL3=round(PL3/10)*10
+
+
+a=(PLa/(PL1+PL2+PL3))^-1
+
+mprintf("a=%.2f sq mm", a)
+disp("the difference in the result is due to error in calculation of abs(I2)=141.12 instead of 144.12")
diff --git a/1076/CH16/EX16.24/16_24.sce b/1076/CH16/EX16.24/16_24.sce new file mode 100755 index 000000000..7a1bf6a6f --- /dev/null +++ b/1076/CH16/EX16.24/16_24.sce @@ -0,0 +1,22 @@ +clear
+clc
+
+Pp=3e3;
+Et=1e7
+Plp=220
+C=2.1
+
+lf=(Et/8760)/Pp
+lf=round(lf*100)/100
+mprintf("\nLoad Factor = %.2f ", lf)
+
+L=(.3*lf) + (.7*lf*lf)
+L=round(L*1000)/1000
+mprintf("\nLoss Factor = %.3f ", L)
+
+Loss=L*Plp
+mprintf("\nLoss = %.1f kW", Loss)
+
+CLY=Loss*8760*C
+mprintf("\nCost of Energy Loss = Rs %.1f ",CLY)
+
diff --git a/1076/CH16/EX16.25/16_25.sce b/1076/CH16/EX16.25/16_25.sce new file mode 100755 index 000000000..456709c42 --- /dev/null +++ b/1076/CH16/EX16.25/16_25.sce @@ -0,0 +1,53 @@ +clear
+clc
+
+xl=.1
+
+st1=10
+vt1a=66
+vt1b=11
+xt1=.1
+n1=vt1b/vt1a
+
+st2=5
+vt2a=11
+vt2b=3.3
+xt2=.08
+n2=vt2b/vt2a
+
+st3=05
+vt3a=11
+vt3b=.415
+xt3=.06
+n3=vt3b/vt3a
+
+Sm=5
+pfm=.8
+
+Sl=1
+pfl=.8
+
+Sb=10
+Vb1=66
+Vb2=Vb1*n1
+Vb3=Vb2*n2
+Vb4=Vb2*n3
+
+Xt2=xt2* (Sb/st2)
+Xt3=xt3* (Sb/st3)
+
+Il=(Sl/Sb) * exp(%i * -acos(pfl))
+Im=(Sm/Sb) * exp(%i * -acos(pfm))
+
+It=Im+Il
+Vt= 1- (It*(%i*(xt1+xl)))- (Il *%i* Xt3)
+
+mprintf("\n(a) Voltage at domestic load is %.3f pu", abs(Vt))
+
+Im=-1 * %i * 5 * Sm / Sb
+It=Il +Im
+Vt= 1- (It*(%i*(xt1+xl)))- (Il *%i* Xt3)
+
+mprintf("\n(b) Voltage at domestic load is %.3f pu", abs(Vt))
+
+
diff --git a/1076/CH16/EX16.26/16_26.sce b/1076/CH16/EX16.26/16_26.sce new file mode 100755 index 000000000..81804b030 --- /dev/null +++ b/1076/CH16/EX16.26/16_26.sce @@ -0,0 +1,28 @@ +clear
+clc
+
+V=400
+Im=40
+pfm=.75
+pfd=.95
+
+Pm=sqrt(3)* V * Im * pfm * 1e-3
+phi1=acos(pfm)
+phi2=acos(pfd)
+kvar1=Pm * tan ( phi1)
+kvar2=Pm * tan ( phi2)
+kvarC=kvar2-kvar1
+KVarC=(abs(kvarC)*10)/10
+Ic= KVarC* 1e3 /(3 * V)
+Ic=round(Ic*1000)/1000
+Xc=round(V*100/Ic)/100
+C=1/(2*%pi*50 * Xc)
+
+mprintf("\nkVAR of Capacitor Bank = %.1f KVAR, Phase Current Of Capacitor Bank=%.2f A, C= %.2f e-6 F/phase",KVarC, Ic, C*1e6)
+
+
+
+
+
+
+
diff --git a/1076/CH16/EX16.27/16_27.sce b/1076/CH16/EX16.27/16_27.sce new file mode 100755 index 000000000..3b1e25e20 --- /dev/null +++ b/1076/CH16/EX16.27/16_27.sce @@ -0,0 +1,23 @@ +clear
+clc
+
+V=400
+Im=40
+pfm=.75
+pfd=.95
+
+Pm=sqrt(3)* V * Im * pfm * 1e-3
+phi1=acos(pfm)
+phi2=acos(pfd)
+kvar1=Pm * tan ( phi1)
+kvar2=Pm * tan ( phi2)
+kvarC=kvar2-kvar1
+Ic= abs(kvarC)* 1e3 /(3 * V)
+Xc=V/Ic
+C=1/(2*%pi*50 * Xc)
+
+
+Imx=Im * sin (acos(pfm))
+Iline= (Im * pfm) - (%i * (Imx - (Ic*sqrt(3))))
+dL=(1-(abs(Iline)/Im)^2)*100
+mprintf("\npercentage reduction in power loss=%.2f",dL)
diff --git a/1076/CH16/EX16.28/16_28.sce b/1076/CH16/EX16.28/16_28.sce new file mode 100755 index 000000000..723b2c50f --- /dev/null +++ b/1076/CH16/EX16.28/16_28.sce @@ -0,0 +1,25 @@ +clear
+clc
+
+t1=15
+P1=1.3
+pf1=.75
+
+t2=9
+P2=.4
+pf2=.8
+
+pfd=.95
+
+kvar1i=P1*1e3*tan (acos(pf1))
+kvar1f=P1*1e3*tan (acos(pfd))
+kvarr1=kvar1i-kvar1f
+
+kvar2i=P2*1e3*tan (acos(pf2))
+kvar2f=P2*1e3*tan (acos(pfd))
+kvarr2=kvar2i-kvar2f
+
+SBC=abs(kvarr2-kvarr1)
+FBC=min(kvarr2,kvarr1)
+
+mprintf("Switch Bank Capacity: %.2f KVAR, Fixed Bank Capacity: %.2f KVAR",SBC, FBC)
diff --git a/1076/CH16/EX16.3/16_3.sce b/1076/CH16/EX16.3/16_3.sce new file mode 100755 index 000000000..33236cf34 --- /dev/null +++ b/1076/CH16/EX16.3/16_3.sce @@ -0,0 +1,22 @@ +clear;
+clc
+
+r1=.02
+r2=.05
+r3=.03
+
+Ic1=100
+Id1=180
+
+Va=255;
+Vb=250
+
+dV=abs(Va-Vb)
+Ia=(dV+(r1*0)+(r2*Ic1)+(r3*(Id1+Ic1)))/(r1+r2+r3)
+
+Ib=-(Ia-(Ic1+Id1))
+
+Vc=Va-Ia*r1
+Vd=Vc-((Ia-Ic1)*r2)
+
+mprintf("IA= %.0f A, IB=%.0fA, Vc=%.2f V, Vd=%.2f V", Ia, Ib, Vc,Vd)
diff --git a/1076/CH16/EX16.4/16_4.sce b/1076/CH16/EX16.4/16_4.sce new file mode 100755 index 000000000..f1fa4e611 --- /dev/null +++ b/1076/CH16/EX16.4/16_4.sce @@ -0,0 +1,29 @@ +clear
+clc
+
+l=300
+l1=120
+l3=120
+l2=l-l1-l3
+U=.25
+Ic1=40
+Id1=60
+
+r=.1/100
+r1=l1*r;
+r2=l2*r
+r3=l3*r
+
+Va=300
+Vb=300
+
+dV=abs(Va-Vb)
+Ia=(dV+(r1*.5*U*l1)+(r2*.5*U*l2)+(r3*.5*U*l3)+(r2*(Ic1+U*l1))+(r3*(Ic1+U*l1 +Id1+U*l2)))/(r1+r2+r3)
+I=Ic1+Id1+(U*l)
+Ib=I-Ia
+
+Vc=Va-(Ia-.5*U*l1)*r1
+Vd=Vb-((Ib-.5*U*l3)*r3)
+
+mprintf("IA= %.1f A, IB=%.1fA, Vc=%.2f V, Vd=%.2f V", Ia, Ib, Vc,Vd)
+
diff --git a/1076/CH16/EX16.5/16_5.sce b/1076/CH16/EX16.5/16_5.sce new file mode 100755 index 000000000..41d04e27f --- /dev/null +++ b/1076/CH16/EX16.5/16_5.sce @@ -0,0 +1,19 @@ +clear
+clc
+
+l=400
+l1=100
+l2=250
+l3=400
+r=.25/1000
+
+I1=100
+I2=120
+I3=80
+
+V=240
+
+dv=r*((l1*I1)+(l2*I2)+(l3*I3))
+Ve=V-dv
+
+mprintf("\nVoltage at end is %.0f V, drop = %.0f V", Ve, dv)
diff --git a/1076/CH16/EX16.6/16_6.sce b/1076/CH16/EX16.6/16_6.sce new file mode 100755 index 000000000..243119169 --- /dev/null +++ b/1076/CH16/EX16.6/16_6.sce @@ -0,0 +1,31 @@ +clear
+clc
+
+l1=100
+l2=250
+l3=400
+r=.25/1000
+l=.125/1000
+
+I1=100
+I2=120
+I3=80
+
+pf1=.7
+pf2=1
+pf3=.8
+
+phi1=acos(pf1)
+phi2=acos(pf2)
+phi3=acos(pf3)
+
+Z1=l1 * ((r*cos(phi1))+(l*sin(phi1)))
+Z2=l2 * ((r*cos(phi2))+(l*sin(phi2)))
+Z3=l3 * ((r*cos(phi3))+(l*sin(phi3)))
+
+V=240
+
+dv=(Z1*I1)+(Z2*I2)+(Z3*I3)
+Ve=V-dv
+
+mprintf("\nVoltage at end is %.2f V, drop = %.2f V", Ve, dv)
diff --git a/1076/CH16/EX16.7/16_7.sce b/1076/CH16/EX16.7/16_7.sce new file mode 100755 index 000000000..ea9a2cf29 --- /dev/null +++ b/1076/CH16/EX16.7/16_7.sce @@ -0,0 +1,30 @@ +clear
+clc
+
+V=240
+
+P1=50e3
+P2=50e3
+P3=200e3
+Pm=500e3
+pfm=.8
+
+
+I1=P1/V
+I2=P2/V
+I3=P3/V
+Im=Pm/(3*V*pfm)
+
+i1=fix(I1*exp(%i *0*%pi/3)*100)/100
+i2=fix((fix(I2*100)/100)*exp(%i *-2*%pi/3)*1000)/1000
+i3=fix(I3*exp(%i *2*%pi/3)*100)/100
+in=abs(round((i1+i2+i3)*100)/100)
+
+iR=sqrt((I1)^2+(Im)^2+(2*I1*Im*pfm))
+iY=sqrt((I2)^2+(Im)^2+(2*I2*Im*pfm))
+iB=sqrt((I3)^2+(Im)^2+(2*I3*Im*pfm))
+
+mprintf("\nIR= %.0f A", iR)
+mprintf("\nIY= %.0f A", iY)
+mprintf("\nIB= %.0f A", iB)
+mprintf("\nIn= %.2f A", in)
diff --git a/1076/CH16/EX16.8/16_8.sce b/1076/CH16/EX16.8/16_8.sce new file mode 100755 index 000000000..da61cd2db --- /dev/null +++ b/1076/CH16/EX16.8/16_8.sce @@ -0,0 +1,32 @@ +clear
+clc
+
+V=230
+
+I1=80
+I2=70
+I3=50
+
+pf1=.8
+pf2=.9
+pf3=1
+
+phi1=acos(pf1)
+phi2=acos(pf2)
+phi3=acos(pf3)
+
+i1=I1*exp(%i *0*%pi/3) * exp(%i * -phi1)
+i2=I2*exp(%i*-2*%pi/3) * exp(%i * -phi2)
+i3=I3*exp(%i *2*%pi/3) * exp(%i * -phi3)
+
+r=.2
+
+in=i1+i2+i3
+
+dvR=r*i1
+dvN=r*in
+
+VR=V+dvR+dvN
+
+mprintf("VR= %.1f V ang (%.1f)deg V", abs(VR), atand(imag(VR)/real(VR)))
+
diff --git a/1076/CH16/EX16.9/16_9.sce b/1076/CH16/EX16.9/16_9.sce new file mode 100755 index 000000000..a47357512 --- /dev/null +++ b/1076/CH16/EX16.9/16_9.sce @@ -0,0 +1,20 @@ +clear
+clc
+
+Vb=240
+
+Ib=100* exp(%i *-1* acos (.8))
+Ia=100* exp(%i *-1* acos (.6))
+
+z=complex(.2, .3)
+
+Va=round((Vb + (Ib * z/2)))
+Isa = (Ia * exp(%i * 1 * atan(imag(Va)/real(Va)))) + Ib
+//Isa=Isa*%i
+
+Vs=Va + (Isa * z/2)
+Vs=round(Vs*100)/100
+vs=abs(Vs)
+phi=atand(imag(Vs)/real(Vs))
+
+mprintf("|Vs| = %.2f A, phase angle between Vs & Vb = %.2f deg", fix(vs*100)/100, phi)
|