diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /2123/CH5 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '2123/CH5')
41 files changed, 1007 insertions, 0 deletions
diff --git a/2123/CH5/EX5.1/Exa_5_1.sce b/2123/CH5/EX5.1/Exa_5_1.sce new file mode 100755 index 000000000..34af6f7b3 --- /dev/null +++ b/2123/CH5/EX5.1/Exa_5_1.sce @@ -0,0 +1,24 @@ +//Example No. 5.1
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+T=10;//turns
+Coil=144;//no. of coils
+R=0.011;//ohm
+fi=0.05;//Wb(flux per pole)
+N=200;//rpm
+par_paths=2;//for wave winding
+T_path=Coil*T/par_paths;//no. of turns in each parallel path
+R_path=R*T_path;//ohm
+Ra=R_path/par_paths;//ohm(armature resistance)
+disp(Ra,"Armature resistance in ohm : ");
+p=12;//poles
+emf=par_paths*Coil*T*p*fi*N/60/2;//V
+R1=1000;//ohm
+IL=emf/R1;//A
+Ia=IL;//A
+T=par_paths*Coil*T*p*fi*Ia/2/%pi/par_paths;//N-m
+disp(T,"Torque in N-m : ");
diff --git a/2123/CH5/EX5.10/Exa_5_10.sce b/2123/CH5/EX5.10/Exa_5_10.sce new file mode 100755 index 000000000..1f957a7a8 --- /dev/null +++ b/2123/CH5/EX5.10/Exa_5_10.sce @@ -0,0 +1,30 @@ +//Example No. 5.10
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+V=240;//V
+Ra=0.4;//ohm
+N1=600;//rpm
+Ifl=25;//A
+Radd=1;//ohm
+//If1=If2
+//T1=T2 leads to If1*Ia1=If2*Ia2: Ia1=Ia2
+Ia1=25;//A
+Ia2=25;//A
+Eb1=V-Ia1*Ra;//V
+Eb2=V-Ia2*(Ra+Radd);//V
+N2=N1*Eb2/Eb1;//rpm
+disp(N2,"Speed at full load torque in rpm : ");
+//T3=2*T1
+//If3=If1
+Ia3=2*Ia1;//A
+Eb3=V-Ia3*(Ra+Radd);//V
+N3=N1*Eb3/Eb1;//rpm
+disp(N3,"Speed at twice the full load torque in rpm : ");
+Eb4=0;//V(at speed zero Eb=0)
+Ia4=V/(Ra+Radd);//V
+T4ByT1=Ia4/Ia1;//(field constant)
+disp("Stalling torque is "+string(T4ByT1)+" times of full load torque.");
diff --git a/2123/CH5/EX5.11/Exa_5_11.sce b/2123/CH5/EX5.11/Exa_5_11.sce new file mode 100755 index 000000000..5dbefc396 --- /dev/null +++ b/2123/CH5/EX5.11/Exa_5_11.sce @@ -0,0 +1,22 @@ +//Example No. 5.11
+clc;
+clear;
+close;
+format('v',5);
+
+//Given Data :
+V=250;//V
+Ra=1;//ohm
+Ia1=25;//A
+N1=900;//rpm
+If=2;//A
+N2=1100;//rpm
+Eb1=V-Ia1*Ra;//V
+//If1*Ia1=If2*Ia2
+//Eb2=V-Ia2*Ra;//V
+//-Ia2^2*Ra+Ia2*V-Eb1*Ia1*N2/N1=0;
+polynomial=[-Ra V -Eb1*Ia1*N2/N1];
+Ia2=roots(polynomial);//A
+Ia2=Ia2(2);//A(wide range not allowed)
+If2=Ia1/Ia2*If;//A
+disp(If2,"New value of field current in A : ");
diff --git a/2123/CH5/EX5.12/Exa_5_12.sce b/2123/CH5/EX5.12/Exa_5_12.sce new file mode 100755 index 000000000..a96b0170d --- /dev/null +++ b/2123/CH5/EX5.12/Exa_5_12.sce @@ -0,0 +1,32 @@ +//Example No. 5.12
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+V=230;//V
+f=50;//Hz
+Rf=200;//ohm
+Ra=0.3;//ohm
+T=50;//N-m
+N=900;//rpm
+Kv=0.8;//V/A-rad/s
+Kt=0.8;//N-m/A^2
+Vm=V*sqrt(2);//V
+Vf=2*Vm/%pi;//V
+If=Vf/Rf;//A
+disp(If,"Field current in A :");
+//T=Kt*If*Ia
+Ia=T/Kt/If;//A
+omega=N*2*%pi/60;//rad/s
+Eb=Kv*omega*If;//V
+Va=Eb+Ia*Ra;//V
+//Va=Vm/%pi*(1+cosd(alfa_a))
+alfa_a=acosd(Va/Vm*%pi-1);//degree
+disp(alfa_a,"Fringe angle of converter in degree : ");
+Pout=Ia*Va;//W
+Iin=sqrt(2/2/180*Ia^2*integrate('1','omega',alfa_a,180));
+VAin=V*Iin;//VA
+pf_in=Pout/VAin;//lagging
+disp(pf_in,"Power factor of convertyer(lagging) : ")
diff --git a/2123/CH5/EX5.13/Exa_5_13.sce b/2123/CH5/EX5.13/Exa_5_13.sce new file mode 100755 index 000000000..88e84131a --- /dev/null +++ b/2123/CH5/EX5.13/Exa_5_13.sce @@ -0,0 +1,25 @@ +//Example No. 5.13
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+V=230;//V
+f=50;//Hz
+Rf=200;//ohm
+Ra=0.25//ohm
+Kv=1.1;//V/A-rad/s
+Kt=1.1;//N-m/A^2
+alfa_a=45;//degree
+Ia=50;//A
+alfa_f=0;
+Vf=2*V*sqrt(2)/%pi*cosd(alfa_f);//V
+Va=2*V*sqrt(2)/%pi*cosd(alfa_a);//V
+If=Vf/Rf;//A
+T=Kt*Ia*If;//N-m
+disp(T,"Torque developed in N-m : ");
+Eb=Va-Ia*Ra-2;//V
+omega=Eb/Kv/If;//rad/s
+N=omega*60/2/%pi;//rpm
+disp(N,"Motor speed in rpm : ");
diff --git a/2123/CH5/EX5.14/Exa_5_14.sce b/2123/CH5/EX5.14/Exa_5_14.sce new file mode 100755 index 000000000..3f5d301b4 --- /dev/null +++ b/2123/CH5/EX5.14/Exa_5_14.sce @@ -0,0 +1,20 @@ +//Example No. 5.14
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+V=400;//V
+Ra=0.3//ohm
+Rf=250;//ohm
+Ia=50;//A
+Kv=1.3;//V/A-rad/s
+N=1200;//rpm
+alfa_f=0;
+Vf=3*sqrt(3)*V*sqrt(2)/sqrt(3)/%pi*cosd(alfa_f);//V
+If=Vf/Rf;//A
+Eb=Kv*If*2*%pi*N/60;//V
+Va=Eb+Ia*Ra;//V
+alfa_a=acosd(Va/3/sqrt(3)/V/sqrt(2)*sqrt(3)*%pi);//degree
+disp(alfa_a,"Fringe angle of converter in degree : ");
diff --git a/2123/CH5/EX5.15/Exa_5_15.sce b/2123/CH5/EX5.15/Exa_5_15.sce new file mode 100755 index 000000000..35863bcf2 --- /dev/null +++ b/2123/CH5/EX5.15/Exa_5_15.sce @@ -0,0 +1,23 @@ +//Example No. 5.15
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+V=500;//V
+Ia=200;//A
+Ra=0.1//ohm
+Kv=1.4;//V/A-rad/s
+Kt=1.4;//N-m/A^2
+If=2;//A
+cycle=0.5;//sec
+Pin=cycle*V*Ia/1000;//KW
+disp(Pin,"Input power in KW : ");
+Va=cycle*V;//V
+Eb=Va-Ia*Ra;//V
+omega=Eb/Kv/2;//rad/s
+N=omega*60/2/%pi;//rpm
+disp(N,"Speed in rpm : ");
+T=Kt*2*Ia;//N-m
+disp(T,"Torque in N-m : ");
diff --git a/2123/CH5/EX5.16/Exa_5_16.sce b/2123/CH5/EX5.16/Exa_5_16.sce new file mode 100755 index 000000000..25d0b04a1 --- /dev/null +++ b/2123/CH5/EX5.16/Exa_5_16.sce @@ -0,0 +1,23 @@ +//Example No. 5.16
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+Ra=0.1//ohm
+Rb=7.5//ohm
+Kv=1.4;//V/A-rad/s
+Ia=120;//A
+If=1.6;//A
+cycle=0.35;//sec
+
+Vavg=Rb*Ia*(1-cycle);//V
+disp(Vavg,"Average voltage across chopper in volt : ");
+Pb=Ia^2*Rb*(1-cycle)^2;//W
+disp(Pb,"Power dissipated in watts : ");
+emf=Vavg+Ra*Ia;//V
+omega=emf/Kv/If;//rad/s
+N=omega*60/2/%pi;//rpm
+disp(N,"Speed in rpm : ");
+//Answer of Pb & speed is wrong in the book.
diff --git a/2123/CH5/EX5.17/Exa_5_17.sce b/2123/CH5/EX5.17/Exa_5_17.sce new file mode 100755 index 000000000..19dc43bd9 --- /dev/null +++ b/2123/CH5/EX5.17/Exa_5_17.sce @@ -0,0 +1,52 @@ +//Example No. 5.17
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+V=220;//V
+f=50;//Hz
+L=0.012;//H
+Ra=0.72;//ohm
+K=2;//V/rad/s
+T=60;//N-m
+alfa=90;//degree
+Va=3*sqrt(3)*V*sqrt(2)/2/%pi*(1+cosd(alfa));//V
+Ia=5;//A
+disp(Ia,"Armature Current in A : ");
+T1=Ia*K;//N-m
+disp(T1,"Torque in N-m : ");
+Eb=Va-Ia*Ra;//V
+omega=Eb/K;//rad/s
+N1=omega*60/2/%pi;//rpm
+disp(N1,"Speed in rpm : ");
+disp("");
+Ia=10;//A
+disp(Ia,"Armature Current in A : ");
+T2=Ia*K;//N-m
+disp(T2,"Torque in N-m : ");
+Eb=Va-Ia*Ra;//V
+omega=Eb/K;//rad/s
+N2=omega*60/2/%pi;//rpm
+disp(N2,"Speed in rpm : ");
+Ia=20;//A
+disp(Ia,"Armature Current in A : ");
+T3=Ia*K;//N-m
+disp(T3,"Torque in N-m : ");
+Eb=Va-Ia*Ra;//V
+omega=Eb/K;//rad/s
+N3=omega*60/2/%pi;//rpm
+disp(N3,"Speed in rpm : ");
+Ia=30;//A
+disp(Ia,"Armature Current in A : ");
+T4=Ia*K;//N-m
+disp(T4,"Torque in N-m : ");
+Eb=Va-Ia*Ra;//V
+omega=Eb/K;//rad/s
+N4=omega*60/2/%pi;//rpm
+disp(N4,"Speed in rpm : ");
+plot([T1 T2 T3 T4],[N1 N2 N3 N4]);
+title('Speed Torque Characteristics');
+xlabel('Torque(N-m)');
+ylabel('speed(RPM)');
diff --git a/2123/CH5/EX5.18/Exa_5_18.sce b/2123/CH5/EX5.18/Exa_5_18.sce new file mode 100755 index 000000000..00d28591c --- /dev/null +++ b/2123/CH5/EX5.18/Exa_5_18.sce @@ -0,0 +1,23 @@ +//Example No. 5.18
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+V=400;//V
+f=50;//Hz
+I=50;//A
+Ra=0.1;//ohm
+K=0.3;//V/rpm
+Ia=5;//A
+alfa=30;//degree
+Vavg=3*sqrt(3)*V*sqrt(2)/sqrt(3)/2/%pi*(1+cosd(alfa));//V
+Eb=Vavg-Ia*Ra;//V
+N=Eb/K;//rpm
+disp(N,"No load speed in rpm : ");
+Speed=1600;//rpm
+Eb=Speed*K;//V
+Vin=Eb+I*Ra;//V
+alfa=acosd(Vin/3/sqrt(3)/V/sqrt(2)*sqrt(3)*2*%pi-1);//degree
+disp(alfa,"Fringe angle in degree : ");
diff --git a/2123/CH5/EX5.19/Exa_5_19.sce b/2123/CH5/EX5.19/Exa_5_19.sce new file mode 100755 index 000000000..7c8f7dafc --- /dev/null +++ b/2123/CH5/EX5.19/Exa_5_19.sce @@ -0,0 +1,29 @@ +//Example No. 5.19
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+V=230;//V
+f=50;//Hz
+Rf=200;//ohm
+Ra=0.25//ohm
+Kv=1.1;//V/A-rad/s
+Kt=1.1;//N-m/A^2
+alfa_a=45;//degree
+Ia=50;//A
+alfa_f=0;
+Vf=2*V*sqrt(2)/%pi*cosd(alfa_f);//V
+Va=2*V*sqrt(2)/%pi*cosd(alfa_a);//V
+If=Vf/Rf;//A
+T=Kt*Ia*If;//N-m
+Eb=Va-Ia*Ra-2;//V
+omega=Eb/Kv/If;//rad/s
+Eg=-Eb;//V
+Va=Eg+Ia*Ra+2;//V
+alfa=acosd(Va/2/V/sqrt(2)*%pi);//degree
+disp(alfa,"Fringe angle to converter in degree : ");
+P=abs(Va)*Ia;//W(power fed back to source)
+disp(P,"Power fed back to source in Watts : ");
+//Answer wrong in the book.
diff --git a/2123/CH5/EX5.2/Exa_5_2.sce b/2123/CH5/EX5.2/Exa_5_2.sce new file mode 100755 index 000000000..e3feeaa1e --- /dev/null +++ b/2123/CH5/EX5.2/Exa_5_2.sce @@ -0,0 +1,18 @@ +//Example No. 5.2
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+Ia=110;//A
+V=480;//volt
+Ra=0.2;//ohm
+p=6;//poles
+C=864;//conductors
+fi=0.05;//Wb(flux per pole)
+back_emf=V-(Ia*Ra);//Volt
+N=back_emf*60*p/C/p/fi;//rpm
+disp(N,"Speed in rpm : ");
+T=C*p*fi*Ia/2/%pi/p;//N-m
+disp(T,"Torque in N-m : ");
diff --git a/2123/CH5/EX5.20/Exa_5_20.sce b/2123/CH5/EX5.20/Exa_5_20.sce new file mode 100755 index 000000000..a48f18551 --- /dev/null +++ b/2123/CH5/EX5.20/Exa_5_20.sce @@ -0,0 +1,15 @@ +//Example No. 5.20
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+V=240;//V
+alfa=100;//degree
+Ra=6//ohm
+Ia=1.8;//A
+Vm=V*sqrt(2);//V
+Vdc=Vm/%pi*(1+cosd(alfa));//Volt
+Eb=Vdc-Ia*Ra;//V
+disp(Eb,"Back emf in volt : ");
diff --git a/2123/CH5/EX5.21/Exa_5_21.sce b/2123/CH5/EX5.21/Exa_5_21.sce new file mode 100755 index 000000000..039f26e01 --- /dev/null +++ b/2123/CH5/EX5.21/Exa_5_21.sce @@ -0,0 +1,27 @@ +0//Example No. 5.21
+clc;
+clear;
+close;
+format('v',5);
+
+//Given Data :
+V1=230;//V
+N1=1500;//rpm
+Ra=1;//ohm
+Ia=10;//A
+T=5;//N-m
+//V=K*omega+Ia*Ra
+K=V1/(N1*2*%pi/60+Ia*Ra);//V-s/rad or N-m/A
+Ia=T/K;//A
+alfa1=30;//degree
+V=2*V1*sqrt(2)/%pi*cosd(alfa1);//Volt
+omega=(V-Ia*Ra)/K;//rad/s
+N=omega*60/2/%pi;//rpm
+disp(N,"Parrt(a) Speed in rpm : ");
+alfa=45;//degree
+N=950;//rpm
+V=2*V1*sqrt(2)/%pi*cosd(alfa);//Volt
+Ia=(V-K*2*%pi/60*N)/Ra;//A
+T=K*Ia;//N-m
+disp(T,"Part(b) Torque in N-m : ");
+//Answer is wrong in the book.
diff --git a/2123/CH5/EX5.22/Exa_5_22.sce b/2123/CH5/EX5.22/Exa_5_22.sce new file mode 100755 index 000000000..1d7ca2728 --- /dev/null +++ b/2123/CH5/EX5.22/Exa_5_22.sce @@ -0,0 +1,26 @@ +//Example No. 5.22
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+V1=500;//V
+N1=1500;//rpm
+Ia=100;//A
+V2=350;//V
+Ra=1.1;//ohm
+alfa=45;//degree
+N2=1200;//rpm
+//V=K*omega+Ia*Ra
+K=(V1-Ia*Ra)/(N1*2*%pi/60);//V-s/rad or N-m/A
+V=3*sqrt(3)*V2*sqrt(2)/2/%pi/sqrt(3)*(1+cosd(alfa));//Volt
+Ia=(V-K*N2*2*%pi/60)/Ra;//A
+disp(Ia,"RMS soirce current in A : ");
+Vin_rms=Ia*sqrt(120/180);//V
+Iavg=Ia/3;//A
+disp(Iavg,"Average thyristor current in A : ");
+Irms=Ia/sqrt(3);//A
+disp(Irms,"RMS thyristor current in A : ");
+pf_in=V*Ia/sqrt(3)/V2/Vin_rms;//lagging
+disp(pf_in,"Input power factor)lagging : ");
diff --git a/2123/CH5/EX5.23/Exa_5_23.sce b/2123/CH5/EX5.23/Exa_5_23.sce new file mode 100755 index 000000000..43b1bbe30 --- /dev/null +++ b/2123/CH5/EX5.23/Exa_5_23.sce @@ -0,0 +1,21 @@ +//Example No. 5.23
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+T1=40;//N-m
+N1=500;//rpm
+J=0.01;//N-m_sec^2/rad
+T2=100;//N-m
+N2=1000;//rpm
+disp("Te=J*d(omega)/dt+D*omega+TL");
+d_omegaBYdt=(T2-T1)/J;//
+//t=omega/d_omegaBYdt+A;
+omega1=N1*2*%pi/60;//rad/s
+t=0;//s(initial time)
+A=t-omega1/d_omegaBYdt;//
+omega2=N2*2*%pi/60;//rad/s
+t=omega2/d_omegaBYdt+A;//s
+disp(t,"Time taken by the motor in sec : ");
diff --git a/2123/CH5/EX5.24/Exa_5_24.sce b/2123/CH5/EX5.24/Exa_5_24.sce new file mode 100755 index 000000000..b2d92a87a --- /dev/null +++ b/2123/CH5/EX5.24/Exa_5_24.sce @@ -0,0 +1,28 @@ +//Example No. 5.24
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+f=400;//Hz
+V=200;//V
+T=30;//N-m
+N=1000;//rpm
+R=0.2;//ohm
+L=2;//mH
+Kv=1.5;//V-sec/rad
+Kt=1.5;//N-m/A
+Ia=T/Kt;//A
+omega=N*2*%pi/60;//rad/s
+Eb=Kv*omega;//V
+alfa=(Eb+Ia*R)/V;
+T=1/f*1000;//ms
+Ton=alfa*T;//ms
+Toff=T-Ton;//ms
+Imax=V/R*[(1-exp(-alfa*T*10^-3*R/(L*10^-3)))/(1-exp(-T*10^-3*R/(L*10^-3)))]-Eb/R;//A
+disp(Imax,"(a) Maximum motor armature current in A : ");
+Imin=V/R*[(exp(alfa*T*R/L)-1)/(exp(T*R/L)-1)]-Eb/R;//A
+disp(round(Imin),"(a) Minimum motor armature current in A : ");
+Iexc=Imax;//A
+disp(Iexc,"(b) Excursion of armature current in A : ");
diff --git a/2123/CH5/EX5.25/Exa_5_25.sce b/2123/CH5/EX5.25/Exa_5_25.sce new file mode 100755 index 000000000..2c62785a7 --- /dev/null +++ b/2123/CH5/EX5.25/Exa_5_25.sce @@ -0,0 +1,20 @@ +//Example No. 5.25
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V=230;//V
+f=50;//Hz
+Rf=1.5;//ohm
+Kt=0.25;//N-m/A
+T=25;//N-m
+Kv=0.25;//V-sec/rad
+Vdc=2*sqrt(2)*V/%pi;//V
+Em=Vdc;//V
+Ia=sqrt(T/Kt);//A
+disp(Ia,"Average motor current in A : ")
+omega_m=(Em-Ia*Rf)/Kv/Ia;//rad/s
+N=omega_m*60/2/%pi;//RPM
+disp(N,"Motor speed in RPM : ");
diff --git a/2123/CH5/EX5.26/Exa_5_26.sce b/2123/CH5/EX5.26/Exa_5_26.sce new file mode 100755 index 000000000..9b4041366 --- /dev/null +++ b/2123/CH5/EX5.26/Exa_5_26.sce @@ -0,0 +1,23 @@ +//Example No. 5.26
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=675;//V
+alfa1=90.5;//degree
+N1=350;//rpm
+Ia1=30;//A
+N2=500;//rpm
+Rf=0.22;//ohm
+Ra=0.22;//ohm
+Ia2=Ia1*N2/N1;//A
+disp(Ia2,"Armature current of converter in A : ");
+Va1=V1*sqrt(2)/%pi*(1+cosd(alfa1));//V
+Eb1=Va1-Ia1*(Ra+Rf);//V
+//Eb1/Eb2=Ia1*N1/(Ia2*N2)
+//Eb2=Va2-Ia2*(Ra+Rf)
+Va2=Eb1*Ia2*N2/(Ia1*N1)+Ia2*(Ra+Rf);//V
+alfa2=acosd(Va2/V1/sqrt(2)*%pi-1);//degree
+disp(alfa2,"Fringe angle of converter in degree : ");
diff --git a/2123/CH5/EX5.27/Exa_5_27.sce b/2123/CH5/EX5.27/Exa_5_27.sce new file mode 100755 index 000000000..400c156df --- /dev/null +++ b/2123/CH5/EX5.27/Exa_5_27.sce @@ -0,0 +1,26 @@ +//Example No. 5.27
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=230;//V
+P=15;//hp
+N=1500;//rpm
+V2=220;//V
+Ke=0.03;//V/A-s
+Kt=0.03;//N-m/A^2
+alfa=45;//degree
+Vm=V1*sqrt(2);//V
+omega=N*2*%pi/60;//rad/s
+T=4*Kt*Vm^2*cosd(alfa)^2/(%pi^2*(Ke*omega)^2);//N-m
+Ia=sqrt(T/Kt);//A
+disp("part (a) : ");
+disp(T,"Torque in N-m : ");
+disp(Ia,"Armature current in A : ");
+disp("part (b) : ");
+Ia=Vm*(1+cosd(alfa))/(%pi*(Ke*omega));//A
+T=Kt*Ia^2;//N-m
+disp(Ia,"Armature current in A : ");
+disp(T,"Torque in N-m : ");
diff --git a/2123/CH5/EX5.28/Exa_5_28.sce b/2123/CH5/EX5.28/Exa_5_28.sce new file mode 100755 index 000000000..c8f1bd542 --- /dev/null +++ b/2123/CH5/EX5.28/Exa_5_28.sce @@ -0,0 +1,22 @@ +//Example No. 5.28
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=230;//V
+N=1000;//rpm
+P=15;//hp
+Rt=0.2;//ohm
+Ke=0.03;//V/A-s
+Kt=0.03;//N-m/A^2
+alfa=30;//degree
+Vm=V1*sqrt(2);//V
+omega=N*2*%pi/60;//rad/s
+V=Vm/%pi*(1+cosd(alfa));//V
+//V=Ke*Ia*omega+Ia*Rt
+Ia=V/(Ke*omega+Rt);//A
+disp(Ia,"Motor current in A : ");
+T=Kt*Ia^2;//N-m
+disp(T,"Torque in N-m : ");
diff --git a/2123/CH5/EX5.29/Exa_5_29.sce b/2123/CH5/EX5.29/Exa_5_29.sce new file mode 100755 index 000000000..45da75754 --- /dev/null +++ b/2123/CH5/EX5.29/Exa_5_29.sce @@ -0,0 +1,19 @@ +//Example No. 5.29
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V=220;//V
+Vin=230;//V
+N1=1500;//rpm
+Ia1=10;//A
+Ra=3;//ohm
+N2=600;//rpm
+E1=V-Ia1*Ra;//V
+E2=E1*N2/N1;//V
+Ia2=Ia1/2;//A(because of Tnew=T/2)
+Vapp=E2+Ia2*Ra;//V
+alfa=acosd(Vapp*%pi/2/sqrt(2)/Vin);//degree
+disp(alfa,"Firing angle of converter in degree : ");
diff --git a/2123/CH5/EX5.3/Exa_5_3.sce b/2123/CH5/EX5.3/Exa_5_3.sce new file mode 100755 index 000000000..379207416 --- /dev/null +++ b/2123/CH5/EX5.3/Exa_5_3.sce @@ -0,0 +1,18 @@ +//Example No. 5.3
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+Ia=100;//A
+V=200;//volt
+N=600;//rpm
+Ra=0.05;//ohm
+Eff=85/100;//
+Ia1=Ia*Eff;//armature current in separately excited dc motor
+emf=V-Ia*Ra;//V(motoring mode induced emf)
+N1=500;//rpm(generating mode speed)
+Gen_emf=emf*N1/N;//V
+Vo=Gen_emf-Ia1*Ra;//V
+disp(Vo,"Voltage of source in Volt : ");
diff --git a/2123/CH5/EX5.30/Exa_5_30.sce b/2123/CH5/EX5.30/Exa_5_30.sce new file mode 100755 index 000000000..7ceddf6b9 --- /dev/null +++ b/2123/CH5/EX5.30/Exa_5_30.sce @@ -0,0 +1,16 @@ +//Example No. 5.30
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V=230;//V
+N=870;//rpm
+Ia=100;//A
+Ra=0.05;//ohm
+T=400;//N-m
+E=V-Ia*Ra;//V
+Vgen=V+Ia*Ra;//V
+N2=N*Vgen/E;//rpm
+disp(N2,"Motor speed in rpm : ");
diff --git a/2123/CH5/EX5.31/Exa_5_31.sce b/2123/CH5/EX5.31/Exa_5_31.sce new file mode 100755 index 000000000..bcec24489 --- /dev/null +++ b/2123/CH5/EX5.31/Exa_5_31.sce @@ -0,0 +1,23 @@ +//Example No. 5.31
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V=220;//V
+P=2.2;//KW
+N1=1000;//rpm
+Ra=2;//ohm
+f=250;//Hz
+alfa=0.9;//cycle
+N2=1200;//rpm
+N3=800;//rpm
+Ia1=P*1000/V;//A
+Ia2=Ia1*N2/N1;//A
+Eb2=alfa*V-Ia2*Ra;//V
+Eb3=Eb2*N3/N2;//V
+Ia3=Ia1*N3/N1;//A
+alfa3=(Eb3+Ia3*Ra)/V;//cycle
+ton=alfa3/f;//sec
+disp(ton,'On time of chopper in sec : ' );
diff --git a/2123/CH5/EX5.32/Exa_5_32.sce b/2123/CH5/EX5.32/Exa_5_32.sce new file mode 100755 index 000000000..8dd67e5b1 --- /dev/null +++ b/2123/CH5/EX5.32/Exa_5_32.sce @@ -0,0 +1,22 @@ +//Example No. 5.32
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V=230;//V
+N1=1000;//rpm
+Ia1=100;//A
+Ra=0.1;//ohm
+Rf=0.1;//ohm
+N2=800;//rpm
+Ia2=sqrt(2)*Ia1;//A(As T2=2*T1 & T proportional to Ia^2)
+Eb1=V-Ia1*(Ra+Rf);//V
+Eb2=N2*Ia2/(N1*Ia1)*Eb1;//V
+//Eb2=Ia2*(Ra+Rf+Rbraking)
+Rbraking=Eb2/Ia2-Ra-Rf;//ohm
+disp(Rbraking,'Braking resistance in ohm : ' );
+Ibraking=Eb2/Rbraking;//A
+disp(Ibraking,'Braking current in A : ' );
+//Braking current is not calculated in the textbook but asked in the example.
diff --git a/2123/CH5/EX5.33/Exa_5_33.sce b/2123/CH5/EX5.33/Exa_5_33.sce new file mode 100755 index 000000000..e297f4d3c --- /dev/null +++ b/2123/CH5/EX5.33/Exa_5_33.sce @@ -0,0 +1,29 @@ +//Example No. 5.33
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+P=6;//poles
+V=220;//V
+f=50;//Hz
+Ra=0.2;//ohm
+Rf=150;//ohm
+Z=150;//no. of conductors
+fi=0.02027;//Wb(flux)
+alfa=0;//degree
+alfa_a=45;//degree
+Ia=25;//A
+A=2;//
+T=Z*P*fi*Ia/(2*%pi*A);//N-m
+disp(T,"Totque in N-m : ");
+Vm=V*sqrt(2);//V
+Vdc=2*Vm/%pi*cosd(alfa_a);//V
+Eb=Vdc-Ia*Ra;//V
+N=Eb*60*A/(Z*P*fi);//rpm
+disp(N,"Speed in rpm : ");
+Pout=Vdc*Ia;//W
+pf=Pout/V/Ia;//lagging
+disp(pf,'Lagging power factor : ' );
+
diff --git a/2123/CH5/EX5.34/Exa_5_35.sce b/2123/CH5/EX5.34/Exa_5_35.sce new file mode 100755 index 000000000..dfe975ba7 --- /dev/null +++ b/2123/CH5/EX5.34/Exa_5_35.sce @@ -0,0 +1,18 @@ +//Example No. 5.35
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=200;//V
+N1=940;//rpm
+Ra=0.02;//ohm
+Ia=100;//A
+N2=500;//rpm
+Eb1=V1-Ia*Ra;//V
+//Eb1/Eb2=N1/N2
+//Eb2=V2-Ia*Ra;//V
+V2=Eb1*N2/N1+Ia*Ra;//V
+cycle=V2/V1;
+disp(cycle,"Duty cycle : ");
diff --git a/2123/CH5/EX5.35/Exa_5_35.sce b/2123/CH5/EX5.35/Exa_5_35.sce new file mode 100755 index 000000000..dfe975ba7 --- /dev/null +++ b/2123/CH5/EX5.35/Exa_5_35.sce @@ -0,0 +1,18 @@ +//Example No. 5.35
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=200;//V
+N1=940;//rpm
+Ra=0.02;//ohm
+Ia=100;//A
+N2=500;//rpm
+Eb1=V1-Ia*Ra;//V
+//Eb1/Eb2=N1/N2
+//Eb2=V2-Ia*Ra;//V
+V2=Eb1*N2/N1+Ia*Ra;//V
+cycle=V2/V1;
+disp(cycle,"Duty cycle : ");
diff --git a/2123/CH5/EX5.36/Exa_5_36.sce b/2123/CH5/EX5.36/Exa_5_36.sce new file mode 100755 index 000000000..611f6799e --- /dev/null +++ b/2123/CH5/EX5.36/Exa_5_36.sce @@ -0,0 +1,22 @@ +//Example No. 5.36
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=220;//V
+Ra=0.05;//ohm
+N1=1000;//rpm
+Ia=100;//A
+N2=500;//rpm
+Eb=V1-Ia*Ra;//V
+Ib=2*Ia;//A
+Rb=(V1+Eb)/Ib-Ra;//ohm
+disp(Rb,"Resistance to be added in ohm : ");
+Tb=Eb/(N1*2*%pi/60)*Ib;//N-m
+disp(Tb,"Initial braking torque in N-m : ");
+Eb2=Eb*N2/N1;//V
+Ib2=(V1+Eb2)/(Ra+Rb);//A
+Tb2=Eb2/(N2*2*%pi/60)*Ib2;//N-m
+disp(Tb2,"Initial braking torque in N-m : ");
diff --git a/2123/CH5/EX5.37/Exa_5_37.sce b/2123/CH5/EX5.37/Exa_5_37.sce new file mode 100755 index 000000000..34501a77f --- /dev/null +++ b/2123/CH5/EX5.37/Exa_5_37.sce @@ -0,0 +1,16 @@ +//Example No. 5.37
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=230;//V
+N1=870;//rpm
+Ia=100;//A
+Ra=0.05;//ohm
+T=400;//N-m
+Eb=V1-Ia*Ra;//V
+Vgen=V1+Ia*Ra;//V
+N2=N1*Vgen/Eb;//rpm
+disp(N2,"Speed in rpm : ");
diff --git a/2123/CH5/EX5.38/Exa_5_38.sce b/2123/CH5/EX5.38/Exa_5_38.sce new file mode 100755 index 000000000..70e82a537 --- /dev/null +++ b/2123/CH5/EX5.38/Exa_5_38.sce @@ -0,0 +1,26 @@ +//Example No. 5.38
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+P=10;//KW
+V1=230;//V
+N1=1200;//rpm
+Ra=0.5;//ohm
+Ke=0.182;//V/rpm
+V2=260;//V
+alfa=30;//degree
+Ia=30;//A
+Vm=V2*sqrt(2);//V
+Vdc=2*Vm/%pi*cosd(alfa);//V
+Eb=Vdc-Ia*Ra;//V
+Kt=Ke*60/2/%pi;//N-m/A
+T=Kt*Ia;//N-m
+disp(T,"Torque in N-m : ");
+N2=Eb/Ke;//rpm
+disp(N2,"Speed in rpm : ");
+Pout=Vdc*Ia;//W
+pf=Pout/V2/Ia;//lagging power factor
+disp(pf,"Lagging power factor : ");
diff --git a/2123/CH5/EX5.39/Exa_5_39.sce b/2123/CH5/EX5.39/Exa_5_39.sce new file mode 100755 index 000000000..fb31c3fe8 --- /dev/null +++ b/2123/CH5/EX5.39/Exa_5_39.sce @@ -0,0 +1,23 @@ +//Example No. 5.39
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+P=2.2;//KW
+V=220;//V
+N1=1000;//rpm
+Ra=2;//ohm
+f=250;//Hz
+alfa=0.9;//duty cycle
+N2=1200;//rpm
+N3=800;//rpm
+Ia1=P*1000/V;//A
+Ia2=Ia1*N2/N1;//A
+Eb1=alfa*V-Ia2*Ra;//V
+Eb2=Eb1*N3/N2;//V
+Ia3=Ia1*N3/N1;//A
+alfa3=(Eb2+Ia3*Ra)/V;//cycle
+ton=alfa3/f;//sec
+disp(ton*1000,'On time of chopper in milli seconds : ' );
diff --git a/2123/CH5/EX5.4/Exa_5_4.sce b/2123/CH5/EX5.4/Exa_5_4.sce new file mode 100755 index 000000000..965acc4d9 --- /dev/null +++ b/2123/CH5/EX5.4/Exa_5_4.sce @@ -0,0 +1,28 @@ +//Example No. 5.4
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+Ia1=10;//A
+V1=200;//volt
+N1=1800;//rpm
+Ra=0.6;//ohm
+Rfield=360;//ohm
+V2=180;//volt
+I_line=20;//A
+//fi2=V2/V1*fi1
+fi2BYfi1=V2/V1;
+//Ia1*fi1=Ia2*fi2
+Ia2=Ia1/fi2BYfi1;//A
+Eb1=V1-Ia1*Ra;//V
+Eb2=V2-Ia2*Ra;//V
+//Eb1/Eb2=fi1*N1/fi2/N2
+N2=N1/(Eb1/Eb2*fi2BYfi1);//rpm
+disp(N2,"Motor speed after supply voltage decreases in rpm : ");
+Ifield=V2/Rfield;//A
+Ia=I_line-Ifield
+//V2=Ia*(R+Ra)
+R=V2/Ia-Ra;//ohm
+disp(R,"Additional resistance in ohm : ");
diff --git a/2123/CH5/EX5.40/Exa_5_40.sce b/2123/CH5/EX5.40/Exa_5_40.sce new file mode 100755 index 000000000..a0212f0bb --- /dev/null +++ b/2123/CH5/EX5.40/Exa_5_40.sce @@ -0,0 +1,20 @@ +//Example No. 5.40
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V=220;//V
+Eff1=85/100;//Efficiency
+Eff2=80/100;//Efficiency
+Load=400;//Kg
+t=2.5;//ms
+Ra=0.1;//ohm
+g=9.81;//constant for gravity acceleration
+Pout=Load*g*t;//W
+IL=Pout/V/Eff1/Eff2;//A
+disp(IL,"Current drawn in ohm : ");
+Eb=V-IL*Ra;//V
+R=(V+Eb)/IL-Ra;//ohm
+disp(R,"Resistance to be added in ohm : ");
diff --git a/2123/CH5/EX5.41/Exa_5_41.sce b/2123/CH5/EX5.41/Exa_5_41.sce new file mode 100755 index 000000000..cbc0a1931 --- /dev/null +++ b/2123/CH5/EX5.41/Exa_5_41.sce @@ -0,0 +1,19 @@ +//Example No. 5.41
+clc;
+clear;
+close;
+format('v',9);
+
+//Given Data :
+V1=220;//V
+N1=1500;//rpm
+I=10;//A
+Ra=3;//ohm
+V2=230;//V
+N2=600;//rpm
+Eb1=V1-I*Ra;//V
+Eb2=Eb1*N2/N1;//V
+Ia=I/2;//A(at half rated torque)
+Vm=V1*sqrt(2);//V
+alfa=acosd((Eb2+Ia*Ra)*%pi/2/Vm);//degree
+disp(alfa,"Firing angle in degree : ");
diff --git a/2123/CH5/EX5.5/Exa_5_5.sce b/2123/CH5/EX5.5/Exa_5_5.sce new file mode 100755 index 000000000..c5c118afd --- /dev/null +++ b/2123/CH5/EX5.5/Exa_5_5.sce @@ -0,0 +1,21 @@ +//Example No. 5.5
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+Ia1=10;//A
+V1=200;//volt
+N1=1800;//rpm
+Ra=0.6;//ohm
+Rfield=360;//ohm
+V2=180;//volt
+I_line=20;//A
+
+Ia=Ia1-V1/Rfield;//A(At changeover time)
+emf=V1-Ia*Ra;//volt
+Ifield=emf/Rfield;//A(At changeover time)
+Iout=Ia1-Ifield;//A
+Rbraking=emf/Iout;//ohm(Braking Resistance)
+disp(Rbraking,"Braking resistance in ohm : ");
diff --git a/2123/CH5/EX5.6/Exa_5_6.sce b/2123/CH5/EX5.6/Exa_5_6.sce new file mode 100755 index 000000000..754b1e2f3 --- /dev/null +++ b/2123/CH5/EX5.6/Exa_5_6.sce @@ -0,0 +1,32 @@ +//Example No. 5.6
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+Ia1=10;//A
+V1=200;//volt
+N1=1800;//rpm
+Ra=0.6;//ohm
+Rfield=360;//ohm
+V2=180;//volt
+I_line=20;//A
+//Part (a)
+Ia=Ia1-V1/Rfield;//A(At changeover time)
+emf=V1-Ia*Ra;//volt
+Ifield=emf/Rfield;//A(At changeover time)
+Iout=Ia1-Ifield;//A
+Rbraking=emf/Iout;//ohm(Braking Resistance)
+
+I_initial=Iout;//A(Inotial current)
+t=30;//sec(time taken to stop)
+I_change_rate=I_initial/t;//A/s
+//i=I_initial-I_change_rate*t , for 0<t<30(during braking time)
+E_dissipated=integrate('(I_initial^2+(I_initial/30)^2/3*t^2-2*I_initial*I_initial/30*t)*Rbraking','t',0,t);//W-s
+disp(E_dissipated,"Part(a) Energy dissipated in watts-sec : ");
+//Part (b)
+//Rbraking=Rbraking-Rbraking/30*t;//ohm
+E=integrate('I_initial^2*(Rbraking-Rbraking/30*t)','t',0,30);//Watt-sec
+disp(E,"Part(b) Energy dissipated in watts-sec : ");
+//calculation of first part is not accurate in the book.
diff --git a/2123/CH5/EX5.7/Exa_5_7.sce b/2123/CH5/EX5.7/Exa_5_7.sce new file mode 100755 index 000000000..47d876791 --- /dev/null +++ b/2123/CH5/EX5.7/Exa_5_7.sce @@ -0,0 +1,22 @@ +//Example No. 5.7
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+I=50;//A
+V=200;//volt
+N=1000;//rpm
+Ra=0.2;//ohm
+Eb=V-I*Ra;//V
+Rt=(V+Eb)/2/I;//ohm(Total resistance required)
+disp(Rt-0.5,"Additional resistance required to limit the current in ohm : ");
+omega_m=N/60*2*%pi;//rad/s
+T=Eb*2*I/omega_m;//N-m
+disp(T,"Braking torque in N-m : ");
+Eb=0;//for speed=0
+I=V/Rt;//A
+//T proportional to I(for separately excited motor)
+T=T*(I/100);//N-m
+disp(T,"Torque when speed decreased to zero in N-m : ");
diff --git a/2123/CH5/EX5.8/Exa_5_8.sce b/2123/CH5/EX5.8/Exa_5_8.sce new file mode 100755 index 000000000..b251dadd9 --- /dev/null +++ b/2123/CH5/EX5.8/Exa_5_8.sce @@ -0,0 +1,34 @@ +//Example No. 5.8
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+Ra=0.2;//ohm
+Rf=100;//ohm
+N=500;//rpm
+Rb=2;//ohm
+E1=100;//V
+If1=2;//A
+If2=2.5;//A
+If3=3;//A
+E2=125;//V
+E3=150;//V
+//Ib=Rf*If1/2
+//Ia=If+Ib;//A
+omega_m=N/60*2*%pi;//rad/s
+Kefi1=E1/omega_m;
+Kefi2=E2/omega_m;
+Kefi3=E3/omega_m;
+T1=E1/omega_m*51*If1;//N-m
+T2=E2/omega_m*51*If2;//N-m
+T3=E3/omega_m*51*If3;//N-m
+Tload=300;//N-m
+Kefi=2.36;
+If=2.482;//A
+Ia=51*If;//A
+E=If*Rf/2+Ia*Ra;//V
+N=E/Kefi;//rad/s
+N=N*60/2/%pi;//rpm
+disp(N,"Speed of motor in rpm : ");
diff --git a/2123/CH5/EX5.9/Exa_5_9.sce b/2123/CH5/EX5.9/Exa_5_9.sce new file mode 100755 index 000000000..348a250d8 --- /dev/null +++ b/2123/CH5/EX5.9/Exa_5_9.sce @@ -0,0 +1,52 @@ +//Example No. 5.9
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+E1=200;//V
+E2=300;//V
+E3=400;//V
+E4=500;//V
+E5=600;//V
+E6=700;//V
+Ia1=20;//A
+Ia2=30;//A
+Ia3=40;//A
+Ia4=50;//A
+Ia5=60;//A
+Ia6=70;//A
+Rt=0.6;//ohm
+Tload=600;//N-m
+omega_m=Tload*2*%pi/60;//rad/s
+Kefi1=E1/omega_m;
+Kefi2=E2/omega_m;
+Kefi3=E3/omega_m;
+Kefi4=E4/omega_m;
+Kefi5=E5/omega_m;
+Kefi6=E6/omega_m;
+T1=E1/omega_m*Ia1;//N-m
+T2=E2/omega_m*Ia2;//N-m
+T3=E3/omega_m*Ia3;//N-m
+T4=E4/omega_m*Ia4;//N-m
+T5=E5/omega_m*Ia5;//N-m
+T6=E6/omega_m*Ia6;//N-m
+subplot(1,2,1);
+plot([Ia1 Ia2 Ia3 Ia4 Ia5 Ia6],[Kefi1 Kefi2 Kefi3 Kefi4 Kefi5 Kefi6])
+title('Ia Vs Kefi');
+xlabel("Ia(A)");
+ylabel("Kefi");
+subplot(1,2,2);
+plot([Ia1 Ia2 Ia3 Ia4 Ia5 Ia6],[T1 T2 T3 T4 T5 T6])
+title('Ia Vs T');
+xlabel("Ia(A)");
+ylabel("T(N-m)");
+//From the graph :
+T=600;//N-m
+Ia=63;//A
+Kefi=9.8;
+E=Kefi*omega_m;//V
+R=E/Ia;//ohm
+Rdb=R-Rt
+disp(Rdb,"Resistance for dynamic braking in ohm : ");
|