diff options
Diffstat (limited to '10/CH6')
-rwxr-xr-x | 10/CH6/EX1/cha6_1.sce | 29 | ||||
-rwxr-xr-x | 10/CH6/EX2/cha6_2.sce | 61 | ||||
-rwxr-xr-x | 10/CH6/EX3/cha6_3.sce | 104 | ||||
-rwxr-xr-x | 10/CH6/EX4/cha6_4.sce | 30 | ||||
-rwxr-xr-x | 10/CH6/EX5/cha6_5.sce | 60 | ||||
-rwxr-xr-x | 10/CH6/EX6/cha6_6.sce | 42 | ||||
-rwxr-xr-x | 10/CH6/EX7/cha6_7.sce | 47 |
7 files changed, 373 insertions, 0 deletions
diff --git a/10/CH6/EX1/cha6_1.sce b/10/CH6/EX1/cha6_1.sce new file mode 100755 index 000000000..4a0d5ce83 --- /dev/null +++ b/10/CH6/EX1/cha6_1.sce @@ -0,0 +1,29 @@ +IMp=500;Pf=0.8;Pf1=0.6;
+SMp=300;Pfs=1;P1=400;
+P2=4;
+
+Power=IMp*Pf
+
+RecPower=IMp*Pf1
+
+SyPower=SMp*Pfs
+
+TotalPower=Power+SyPower
+
+RecPower
+
+ComplesPower=sqrt(TotalPower^2+RecPower^2)
+
+PowerFactor=(TotalPower/ComplesPower)
+
+Power1=sqrt(Power^2-RecPower^2)
+
+KVAR=RecPower-Power1
+
+NewKVA=sqrt(TotalPower^2+KVAR)
+
+PowerFactor1=(TotalPower/NewKVA)
+
+Ism=P1/(sqrt(3)*P2)
+
+PowerfactorSYS=SyPower/P1
diff --git a/10/CH6/EX2/cha6_2.sce b/10/CH6/EX2/cha6_2.sce new file mode 100755 index 000000000..6fd11e4fa --- /dev/null +++ b/10/CH6/EX2/cha6_2.sce @@ -0,0 +1,61 @@ +V=14000;Phase=3;Ra=0.07;V1=10;Is=490;Pf=0.8;
+If=200;Vl=18000;
+
+Vb=V/sqrt(Phase)
+
+Ib=(V1*10^6)/(sqrt(Phase)*V)
+
+Zb=Vb/Ib
+
+Zsun=(Vl/sqrt(3))/Is
+
+Xsun=sqrt(Zsun^2-Ra^2)
+
+Xsun=Xsun/Zb
+
+Zssa=(V/sqrt(3))/Is
+
+Xssa=Zssa/Zb
+
+a=cos(0.8)
+
+Deg=a*180/%pi
+
+Zs=Zssa/Ra
+
+Zs=atan(Zs)
+
+Angle=Zs*180/%pi
+
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+
+[a,b]=polar2rect(1,0)
+
+X1=a+%i*b
+
+[c,d]=polar2rect(1,-Deg)
+
+X2=c+%i*d
+
+[e,f]=polar2rect(0.84,Angle)
+
+X3=e+%i*f
+
+X=X1+(X2*X3)
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[I,Angle]=rect2polar(1.54,0.64)
+
+Ef=I*V
+
+If=I*If
+
+
+
\ No newline at end of file diff --git a/10/CH6/EX3/cha6_3.sce b/10/CH6/EX3/cha6_3.sce new file mode 100755 index 000000000..2c6fbc4b6 --- /dev/null +++ b/10/CH6/EX3/cha6_3.sce @@ -0,0 +1,104 @@ +function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+V=208;Poles=4;F=60;Phase=3;Vol=5000;Xs=8;
+
+Vt=V/sqrt(Phase)
+
+Ia=Vol/(sqrt(Phase)*V)
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+[x1,y1]=polar2rect(Vt,0)
+
+[x2,y2]=polar2rect(Ia,-36.9)
+
+[x3,y3]=polar2rect(8,90)
+
+X1=x1+%i*0
+
+X2=x2+%i*0
+
+X2=x2+%i*y2
+
+X3=x3+%i*y3
+
+X=X1+(X2*X3)
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[Ef,Angle]=rect2polar(186.7,88.7)
+
+Newvol=1.2*Ef
+
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+
+[x1,y1]=polar2rect(Newvol,21)
+
+[x2,y2]=polar2rect(Vt,0)
+
+[x3,y3]=polar2rect(Xs,90)
+
+X1=x1+%i*y1
+
+X1=x2+%i*y2
+
+X1=x3+%i*y3
+
+X=(X1-X2)/X3
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[Ia,Angle1]=rect2polar(11.11,-13.93)
+
+Pf=cos(%pi*51.5/180)
+
+a=sin(%pi*51.5/180)
+
+ReactiveKva=(Phase*Vt*Ia*a*10^-3)
+
+Pmax=(Phase*Ef*Vt)/Xs
+
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+
+[x1,y1]=polar2rect(206.9,90)
+
+[x2,y2]=polar2rect(120,0)
+
+[x3,y3]=polar2rect(8,90)
+
+X1=x1+%i*y1
+
+X2=x2+%i*y2
+
+X3=x3+%i*y3
+
+X=(X1-X2)/X3
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[Ia,Angle2]=rect2polar(25.8,15)
+
+tan=Vt/Ef
+
+Ang=tan*180/%pi
+
+Pf=cos(%pi*Ang/180)
+
\ No newline at end of file diff --git a/10/CH6/EX4/cha6_4.sce b/10/CH6/EX4/cha6_4.sce new file mode 100755 index 000000000..1c792025f --- /dev/null +++ b/10/CH6/EX4/cha6_4.sce @@ -0,0 +1,30 @@ +V=208;F=60;Phase=3;Power=3000;Xs=8;
+
+Vt=V/sqrt(Phase)
+
+Ia=Power/(Phase*Vt)
+
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+[x1,y1]=polar2rect(120,0)
+
+[x2,y2]=polar2rect(8.33,0)
+
+[x3,y3]=polar2rect(8,90)
+
+X=X1-(X2*X3)
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[Ef,Angle]=rect2polar(120,-66.64)
+
+Pmax=(Phase*Ef*Vt)/Xs
+
+Ws=(1800/F)*2*%pi
+
+Tmax=Pmax/Ws
diff --git a/10/CH6/EX5/cha6_5.sce b/10/CH6/EX5/cha6_5.sce new file mode 100755 index 000000000..9e01b753f --- /dev/null +++ b/10/CH6/EX5/cha6_5.sce @@ -0,0 +1,60 @@ +V=460;Phase=3;N=1200;Hp=125;
+Ra=0.078;Xal=0.15;Xar=1.85;Nre=28.2;Nse=28.2;
+
+Ia=sqrt(Phase)*V
+
+Ia=121.4
+
+Vt=V/sqrt(Phase)
+
+Ea=Vt-(Ia*Ra)
+
+Xs=Xal+Xar
+
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+
+[x1,y1]=polar2rect(Ea,0)
+
+[x2,y2]=polar2rect(Xs,90)
+
+X1=x1+%i*y1
+
+X2=x2+%i*y2
+
+Im=X1/X2
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[Im,Angle1]=rect2polar(7.84D-15,-128.0)
+
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+
+[x1,y1]=polar2rect(Im,Angle1)
+
+[x2,y2]=polar2rect(Ia,0)
+
+X1=x1+%i*y1
+
+X2=x2+%i*y2
+
+X=X1-X2
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[a,b]=rect2polar(-121.4,-128)
+
+n=sqrt(2)/Phase*Nre
+
+If=(a/n)*(Xs/Xar)
diff --git a/10/CH6/EX6/cha6_6.sce b/10/CH6/EX6/cha6_6.sce new file mode 100755 index 000000000..b8dd10c2a --- /dev/null +++ b/10/CH6/EX6/cha6_6.sce @@ -0,0 +1,42 @@ +V=11;Phase=3;F=60;IncExe=150;DecExe=50;Xs=10;
+Power=80000;
+
+Vt=V/sqrt(Phase)
+
+Vt=Vt*1000
+
+Ef=IncExe/100
+
+Ef1=DecExe/100
+
+Ia=(Vt-(Ef*Vt))/(Xs)
+
+Pf=cos(90/90*%pi/2)
+
+Ia1=(Vt-(Ef1*Vt))/(Xs)
+
+Pf1=cos(90/90*%pi/2)
+
+Ia=Power/(Phase*Vt)
+
+function[x,y]=polar2rect(r,theta)
+x=r*cos(theta*%pi/180);
+y=r*sin(theta*%pi/180);
+endfunction
+
+[x1,y1]=polar2rect(4.2,0)
+
+[x2,y2]=polar2rect(10,90)
+
+X1=x1+%i*y1
+
+X2=x2+%i*y2
+
+X=X1*X2
+
+function[r,theta]=rect2polar(x,y)
+r=sqrt(x^2+y^2);
+theta=atan(y/x)*180/%pi;
+endfunction
+
+[a,b]=rect2polar(Vt,42)
diff --git a/10/CH6/EX7/cha6_7.sce b/10/CH6/EX7/cha6_7.sce new file mode 100755 index 000000000..b95fc4382 --- /dev/null +++ b/10/CH6/EX7/cha6_7.sce @@ -0,0 +1,47 @@ +Vt=1;Ia=1;Xd=0.8;Xq=0.4;Loss=0.15;Angle=36.9;
+
+a=cos(%pi*Angle/180)
+
+b=sin(%pi*Angle/180)
+
+TanDeg=(Vt*Xq*a)/(Vt+(Ia*Xq*b))
+
+z=atan(TanDeg)
+
+Deg=(z*%pi/180)
+
+Deg=(z*180/%pi)
+
+Angl=Angle-Deg
+
+Id=Ia*sin(%pi*Angl/180)
+
+Iq=Ia*cos(%pi*Angl/180)
+
+Ef=(Vt*cos(%pi*Deg/180))-(Id*Xq)
+
+Pf=((Vt*Ef)/Xd)*sin(%pi*Deg/180)
+
+Pr=(Vt^2*(Xd-Xq)/(2*Xd*Xq))*sin(%pi*45.6/180)
+
+Pout=Vt*Ia*a
+
+Prmax=(Vt^2*(Xd-Xq))/(2*Xd*Xq)
+
+v=asin(Loss/Prmax)/2
+
+Deg1=(v*180/%pi)
+
+Id=Vt*cos(%pi*Deg1/180)/Xd
+
+Iq=Vt*sin(%pi*Deg1/180)/Xq
+
+Ia=(Id^2+Iq^2)^(1/2)
+
+ Ang=atan(Id/Iq)
+
+ Ang=(Ang*180/%pi)
+
+Phi=Ang+Deg1
+
+Pf=cos(%pi*Phi/180)
|