initial commit / add all books

16 files changed, 336 insertions, 0 deletions

diff --git a/503/CH5/EX5.1/ch5_1.sci b/503/CH5/EX5.1/ch5_1.sci
new file mode 100755
index 000000000..f50bde3c2
--- /dev/null
+++ b/503/CH5/EX5.1/ch5_1.sci
@@ -0,0 +1,20 @@
+// To calculate harmanic factor for stator

+

+clc;

+S=36;    //no of slots

+q=3;    //no of phases

+p=4;    //no of poles

+m=S/(q*p);    //slots/pole/phase

+g=180*p/S;    //gamma elec

+function [k]=bfctr(n)

+    k=sind(m*n*g/2)/(m*sind(n*g/2));

+endfunction

+

+K_b=bfctr(1);

+disp(K_b,'K_b(fundamental)');

+

+K_b=bfctr(3);

+disp(K_b,'K_b(third harmonic)');

+

+K_b=bfctr(5);

+disp(K_b,'K_b(fifth harmonic)');
\ No newline at end of file
diff --git a/503/CH5/EX5.10/ch5_10.sci b/503/CH5/EX5.10/ch5_10.sci
new file mode 100755
index 000000000..9b8c4fb4d
--- /dev/null
+++ b/503/CH5/EX5.10/ch5_10.sci
@@ -0,0 +1,27 @@
+// (a)to calculate field current and flux/pole(b)to calculate open ckt ph and line voltages

+// (c)to caculate field current

+

+clc;

+B_peak=1.65;

+g=.008;

+u_o=4*%pi*10^-7;

+P=4;

+K_b=.957;

+N_field=364/2;

+I_f=B_peak*%pi*g*P/((4*u_o)*(K_b*N_field));

+disp(I_f,'field current(A)');

+l=1.02;    //rotor length

+r=.41/2;    //rotor radius

+phi=(4/P)*B_peak*l*r;

+disp(phi,'flux/pole(Wb)');

+N_ph=3*11*P/2;

+ga=60/3;    //slot angle

+m=3;

+f=50;

+K_b=sind(m*ga/2)/(m*sind(ga/2));    //breadth factor

+E_ph=sqrt(2)*%pi*K_b*f*N_ph*phi;

+disp(E_ph,'E_ph(V)');

+E_line=sqrt(3)*E_ph;

+disp(E_line,'E_line(V)');

+I_fnew=.75*I_f;

+disp(I_fnew,'I_f(new)(A)');
\ No newline at end of file
diff --git a/503/CH5/EX5.11/ch5_11.sci b/503/CH5/EX5.11/ch5_11.sci
new file mode 100755
index 000000000..5042cb8f8
--- /dev/null
+++ b/503/CH5/EX5.11/ch5_11.sci
@@ -0,0 +1,20 @@
+// to find fundamental mmf wave,speed and its peak value

+

+clc;

+p=4;

+S=60;

+g=180*p/S;

+ph=3;

+m=S/(p*ph);    //slots/pole/phase

+K_b=sind(m*g/2)/(m*sind(g/2));    //breadth factor

+I_L=48;

+I_P=I_L/sqrt(3);

+I_Pmax=I_P*sqrt(2);

+c=24;        //conductors

+N_ph=S*c/(ph*2);      //turns/phase

+F_m=(4/%pi)*K_b*(N_ph/p)*I_Pmax;

+disp(F_m,'F_m(AT/pole)');

+F_peak=(3/2)*F_m;

+disp(F_peak,'F_peak(AT/pole)');

+n=120*f/P;

+disp(n,'speed(rpm)');

diff --git a/503/CH5/EX5.12/ch5_12.sci b/503/CH5/EX5.12/ch5_12.sci
new file mode 100755
index 000000000..0eaf9c428
--- /dev/null
+++ b/503/CH5/EX5.12/ch5_12.sci
@@ -0,0 +1,10 @@
+// to calculate resultant air gap flux/pole

+

+clc;

+F1=400;

+F2=850;

+a=123.6;

+Fr=sqrt(F1^2+F2^2+2*F1*F2*cosd(a));

+P=1.408*10^-4;        //permeance/pole

+phi_r=P*Fr;

+disp(phi_r,'air gap flux/pole(Wb)');
\ No newline at end of file
diff --git a/503/CH5/EX5.13/ch5_13.sci b/503/CH5/EX5.13/ch5_13.sci
new file mode 100755
index 000000000..77fea79d9
--- /dev/null
+++ b/503/CH5/EX5.13/ch5_13.sci
@@ -0,0 +1,36 @@
+//To calculate resultant AT/pole and peak air gap flux density, rotor AT/pole, stator AT and its angle with the resultant AT, stator currrent

+

+clc;

+ph=3;

+S=36;

+c=8*2;

+p=4;

+f=50;

+N_ph=S*c/(ph*2);      //turns/phase

+ga=180*p/S;

+m=S/(p*ph);    //slots/pole/phase

+K_b=sind(m*ga/2)/(m*sind(ga/2));    //breadth factor

+V_L=400;

+V_ph=V_L/sqrt(3);

+phi_r=V_ph/(4.44*K_b*f*N_ph);

+disp(phi_r,'phi_r(Wb/pole)');

+D=.16;

+l=0.12;

+PA=%pi*l*D/4;        //pole area

+B_rav=phi_r/PA;

+B_rpeak=(%pi/2)*B_rav;

+g=2*10^-3;

+u_o=4*%pi*10^-7;

+F_r=g*B_rpeak/u_o;

+disp(F_r,'F_r(AT/pole)');

+T=60;        //torque(Nm)

+d=26;

+F2=T/((%pi/2)*(p/2)^2*phi_r*sind(d));

+disp(F2,'F2(AT/pole)');

+F1=sqrt(F2^2+F_r^2-2*F2*F_r*sind(d));

+disp(F1,'F1(AT/pole)');

+w=acosd((F1^2+F_r^2-F2^2)/(2*F1*F_r));

+disp(w,'angle(deg)');

+K_w=K_b;

+I_a=F1/((3/2)*(4*sqrt(2)/%pi)*K_w*(N_ph/p));

+disp(I_a,'I_a(A)');
\ No newline at end of file
diff --git a/503/CH5/EX5.14/ch5_14.sci b/503/CH5/EX5.14/ch5_14.sci
new file mode 100755
index 000000000..8d39a413f
--- /dev/null
+++ b/503/CH5/EX5.14/ch5_14.sci
@@ -0,0 +1,31 @@
+//to determine in F2,peak rotor AT, max torque, ele i/p at max torque(motoring mode),open ckt voltage(generating mode)

+

+clc;

+disp('motoring mode');

+K_w=.976;

+N_pole=746;

+p=4;

+I_f=20;

+F2=(4/%pi)*K_w*(N_pole/p)*I_f;

+disp(F2,'F2(AT)');

+B_r=1.6;

+D=.29;

+l=.35;

+T_max=(p/2)*(%pi*D*l/2)*F2*B_r;

+disp(T_max,'T_max');

+f=50;

+w_m=4*%pi*f/p;

+P_in=T_max*w_m;

+disp(P_in,'P_in(W)');

+

+disp('generating mode');

+m=S/(3*p);

+ga=180*p/S;

+K_b=sind(30)/(3*sind(15/2));

+K_w=K_b;

+u_o=4*%pi*10^-7;

+phi_r=((2*D*l/p)*(u_o/g))*F2;

+N_ph=20*p*4/2;

+E_ph=4.44*K_b*f*N_ph*phi_r;

+E_l=sqrt(3)*E_ph;

+disp(E_l,'E_l(V)');
\ No newline at end of file
diff --git a/503/CH5/EX5.15/ch5_15.sci b/503/CH5/EX5.15/ch5_15.sci
new file mode 100755
index 000000000..bbf6dee51
--- /dev/null
+++ b/503/CH5/EX5.15/ch5_15.sci
@@ -0,0 +1,12 @@
+// to find motor speed

+

+clc;

+n=1500;        //speed of sync generator

+p=4;

+f=n*p/120;

+

+p_im=6;

+n_s=120*f/p_im;

+s=0.05;        //slip

+n_im=(1-s)*n_s;

+disp(n_im,'speed of induction motor(rpm)');
\ No newline at end of file
diff --git a/503/CH5/EX5.16/ch5_16.sci b/503/CH5/EX5.16/ch5_16.sci
new file mode 100755
index 000000000..462bc65a8
--- /dev/null
+++ b/503/CH5/EX5.16/ch5_16.sci
@@ -0,0 +1,23 @@
+//to find voltage available b/w slip rings and its freq

+

+clc;

+disp('(a)');

+f=50;

+p=6;

+n_s=120*f/p;

+n=-1000;

+s=(n_s-n)/n_s;

+f_s=f*s;

+disp(f_s,'slip freq(Hz)');

+v2=100;

+V2=s*v2;

+disp(V2,'slip ring voltage(V)');

+

+disp('(b)');

+n=1500;

+s=(n_s-n)/n_s;

+f_s=abs(f*s);

+disp(f_s,'slip freq(Hz)');

+v2=100;

+V2=s*v2;

+disp(V2,'slip ring voltage(V)');

diff --git a/503/CH5/EX5.18/ch5_18.sci b/503/CH5/EX5.18/ch5_18.sci
new file mode 100755
index 000000000..d8e275fcc
--- /dev/null
+++ b/503/CH5/EX5.18/ch5_18.sci
@@ -0,0 +1,21 @@
+//to find no of poles, slip and freq of rotor currents at full load, motor speed at twice of full load

+

+clc;

+n_s=600;

+f=50;

+P=120*f/n_s;

+disp(p,'no of poles');

+n=576;

+s=(n_s-n)/n_s;

+disp(s,'slip');

+f2=s*f;

+n_r=s*n_s;

+disp(n_r,'rotor speed wrt rotating field(rpm)');

+ss=f2*s;

+n=(1-ss)*n_s;

+disp(n,'motor speed(rpm)');

+nn=528;

+s_old=s;

+s_new=(n_s-nn)/n_s;

+fac=s_new/s_old;

+disp(fac,'factor is');
\ No newline at end of file
diff --git a/503/CH5/EX5.19/ch5_19.sci b/503/CH5/EX5.19/ch5_19.sci
new file mode 100755
index 000000000..0ddd488ca
--- /dev/null
+++ b/503/CH5/EX5.19/ch5_19.sci
@@ -0,0 +1,21 @@
+// to calculate amplitude of travelling wave mmf,peak value of air flux density, velocity of wave, current freq at some desired velocity

+

+clc;

+K_w=.925;

+N_ph=48;

+I=750/sqrt(2);

+wndnglgth=2;

+wavelgth=wndnglgth/0.5;

+p=2*wavelgth;

+F_peak=(3/2)*(4*sqrt(2)/%pi)*K_w*(N_ph/p)*I;

+disp(F_peak,'F_peak(A/m)');

+g=.01;

+u_o=4*%pi*10^-7;

+B_peak=u_o*F_peak/g;

+disp(B_peak,'B_peak(T)');

+f=25;

+B=.5;

+v=f*B;    disp(v,'velocity(m/s)');

+vv=72*10^3/3600;        //given velocity

+f=vv/0.5;

+disp(f,'freq(Hz)');
\ No newline at end of file
diff --git a/503/CH5/EX5.2/ch5_2.sci b/503/CH5/EX5.2/ch5_2.sci
new file mode 100755
index 000000000..45b0164d3
--- /dev/null
+++ b/503/CH5/EX5.2/ch5_2.sci
@@ -0,0 +1,20 @@
+// to find the frequency and phase and line voltages

+

+clc;

+n=375;        //speed in rpm

+p=16;        //no of poles

+f=n*p/120;

+disp(f,'freq(Hz)');

+S=144;        //no of slots

+c=10;        //no of conductors/slot

+t=S*c/2;        //no of turns

+ph=3;

+N_ph=t/ph;    //no of turns/ph

+g=180*p/S;    //slots angle

+m=S/(p*ph);    //slots/pole/phase

+K_b=sind(m*g/2)/(m*sind(g/2));    //breadth factor

+phi=0.04;    //flux per pole

+E_p=4.44*K_b*f*N_ph*phi;

+disp(E_p,'phase voltage(V)');

+E_l=sqrt(3)*E_p;

+disp(E_l,'line voltage(V)');

diff --git a/503/CH5/EX5.3/ch5_3.sci b/503/CH5/EX5.3/ch5_3.sci
new file mode 100755
index 000000000..64a8f2207
--- /dev/null
+++ b/503/CH5/EX5.3/ch5_3.sci
@@ -0,0 +1,22 @@
+// to find the phase and line voltages

+

+clc;

+f=50;    //freq

+n=600;    //speed in rpm

+p=120*f/n;

+ph=3;

+m=4;        //slots/pole/ph

+S=p*ph*m;    //slots

+t=12;        //turns per coil

+N_ph=S*t/ph;

+g=180*p/S;

+K_b=sind(m*g/2)/(m*sind(g/2));    //breadth factor

+cp=10;        //coil pitch

+pp=S/cp;        //pole pitch

+theta_sp=(pp-cp)*g;    //short pitch angle

+K_p=cosd(theta_sp/2);

+phi=.035; 

+E_p=4.44*K_b*K_p*f*N_ph*phi;

+disp(E_p,'phase voltage(V)');

+E_l=sqrt(3)*E_p;

+disp(E_l,'line voltage(V)');

diff --git a/503/CH5/EX5.4/ch5_4.sci b/503/CH5/EX5.4/ch5_4.sci
new file mode 100755
index 000000000..cc120a890
--- /dev/null
+++ b/503/CH5/EX5.4/ch5_4.sci
@@ -0,0 +1,17 @@
+// to calculate flux/pole

+

+clc;

+S=42;

+p=2;

+ph=3;

+m=S/(p*ph);    //slots/pole/phase

+g=180*p/S;    //slots angle

+K_b=sind(m*g/2)/(m*sind(g/2));    //breadth factor

+cp=17;

+pp=S/p;

+theta_sp=(pp-cp)*g;    //short pitch angle

+K_p=cosd(theta_sp/2);

+N_ph=S*2/(ph*p*2);        //2 parallel paths

+E_p=2300/sqrt(3);

+phi=E_p/(4.44*K_b*K_p*f*N_ph);

+disp(phi,'flux/pole(Wb)');

diff --git a/503/CH5/EX5.5/ch5_5.sci b/503/CH5/EX5.5/ch5_5.sci
new file mode 100755
index 000000000..3507a6566
--- /dev/null
+++ b/503/CH5/EX5.5/ch5_5.sci
@@ -0,0 +1,18 @@
+// to calculate  useful flux/pole and ares of pole shoe

+

+clc;

+p=1500*1000;    //power

+v=600;

+I_a=p/v;

+cu=25*1000;    //copper losses

+R_a=cu/I_a^2;

+E_a=v+I_a*R_a;

+n=200;

+Z=2500;

+p=16;

+A=16;

+phi=E_a*60*A/(p*n*Z);

+disp(phi,'flux/pole(Wb)');

+fd=0.85;    //flux density

+a=phi/fd;

+disp(a,'area of pole shoe(m*m)');
\ No newline at end of file
diff --git a/503/CH5/EX5.6/ch5_6.sci b/503/CH5/EX5.6/ch5_6.sci
new file mode 100755
index 000000000..16bd6b9ef
--- /dev/null
+++ b/503/CH5/EX5.6/ch5_6.sci
@@ -0,0 +1,17 @@
+// To calculate em power developed,mech power fed, torque provided by primemover

+

+clc;

+phi=32*10^-3;    //flux/pole

+n=1600;        //speed in rpm

+Z=728;        //no of conductors

+p=4;

+A=4;

+E_a=phi*n*Z*(p/A)/60;

+I_a=100;

+P_e=E_a*I_a;

+disp(P_e,'electromagnetic power(W)');

+P_m=P_e;

+disp(P_m,'mechanical power(W) fed');

+w_m=2*%pi*n/60;

+T=P_m/w_m;

+disp(T,'primemover torque(Nm)');

diff --git a/503/CH5/EX5.9/ch5_9.sci b/503/CH5/EX5.9/ch5_9.sci
new file mode 100755
index 000000000..e5e0fd270
--- /dev/null
+++ b/503/CH5/EX5.9/ch5_9.sci
@@ -0,0 +1,21 @@
+// To determine peak value of fundamental mmf

+

+clc;

+f=50;

+n_s=300;

+p=120*f/n_s;

+P=400*1000;    //power

+V=3300;

+I_L=P/(sqrt(3)*V);

+I_P=I_L;

+I_m=sqrt(2)*I_P;    //max value of phase current

+S=180;

+g=180*p/S;

+ph=3;

+m=S/(p*ph);    //slots/pole/phase

+K_b=sind(m*g/2)/(m*sind(g/2));    //breadth factor

+c=8;        //conductors/1 coil side

+N_ph=S*c/(ph*2);      //turns/phase

+F_m=(4/%pi)*K_b*(N_ph/p)*I_m;

+F_peak=(3/2)*F_m;

+disp(F_peak,'peak mmf(AT/pole)');