disp('DESIGN OF OVERALL DIMENSION OF TRANSFORMER\n' );

disp('Enter the following values: ' )

KVA=input('KVA rating: ');
disp('voltage per turn calculation give k :\n1 0.75-0.85 for 1phase core\n2 0.45 for
3phase core\n3 1-1.2 for 1phase shell\n4 1.3 for 3phase\n' );
k=input('k: ');
Et=k*sqrt(KVA);
f=input('Line frequency: ');
m=input('Number of phases: ');
Bm=input('Flux density: ');
Ki=input('Stacking factor: ');
disp('The values of peak flux per pole, net iron area and Gross iron area calculated
are');
PHlm=Et/(4.44*f);
Ai=PHlm/Bm%net iron area;
Agi=Ai/Ki %gross core section area;
c=input('Enter the type of core:\n1)Square\n2)Stepped\n3)3-Stepped\ n4)4-Stepped
\n');
switch c
case 1
ct=0.45;
d=sqrt(Ai/ct) % dia od circumsribing circle
a=sqrt(0.5)*d %width of the largest stamping
case 2
ct=0.56;
d=sqrt(Ai/ct)
a=0.85*d %dimensions
b=0.53*d
case 3
ct=0.6;
d=sqrt(Ai/ct)
a=0.42*d %dimensions
b=0.7*d
c=0.9*d
case 4
ct=0.62;
d=sqrt(Ai/ct)
a=0.36*d %dimensions
b=0.36*d
c=0.78*d
r=0.92*d
end;
disp('WINDOW DESIGN OF TRANSFORMER\n' )
KV=input('Primary winding voltage: ');
deltap=input('Primary current density: ');
Kw=input('Enter window space factor: ' )
Aw=(KVA*1000)/(3.33*Bm*f*Kw*deltap*Ai) %window area
ratiohw=input('Ratio - height to width of window in the range of 2-4: ' ); %Range 2-4
disp('The window width and window height are' );
Ww=sqrt(Aw/ratiohw) %window width
Hw=Ww*ratiohw%window height
disp('YOKE DESIGN OF TRANSFORMER\n' )
ratioyl=input('ratio - area of yoke to limbs: ' );
Dy=input('Depth of yoke: ');
disp('The flux density in yoke,yoke area, gross yoke area and height of yoke are
calculated as');
FDy=Bm/ratioyl%flux density in yoke
Ay=ratioyl*(PHlm/Bm) %yoke area
Agy=Ay/Ki %gross area of yoke
Hy=Agy/Dy%height of yoke
disp('OVERALL DIMENSION OF TRANSFORMER\n' )
disp('The distance betwn core centers,height, width and depth of transformers are
obtained as');
D=d+Ww%dist between adjacent core centres
H=Hw+2*Hy%height of frame
W=2*D+Dy%width of frame
Df=Dy%depth of frame
disp('Given data : 3ph Slip-ring motor, star connected, best power factor desired' )
KW=100
V=3300
f=50 %frequency
P=12 % number of poles
pf=0.9 %power factor
Bav=0.4 %Average gap density (Wb/m^2)
disp('Slot loading<=500A')
Kw=0.96 %winding factor
Ac=25000 %Electric loading (ac/m)
eff=0.9 %efficiency
Q=(KW/eff)/pf%input KVA
C0=11*Bav*Ac*Kw*0.001 %Output coefficient
ns=(2*f)/P %synchronous speed in rps
x=Q/(C0*ns)
disp('For best power T=sqrt(0.18L)' );
y=(0.18*P*P)/(pi*pi)
L=sqrt(x/y)
D=sqrt(y*L)
Et=V/(sqrt(3))
Fm=(Bav*pi*D*L)/P
Ts=round(Et/(4.44*f*Fm*Kw))
disp('Stator slot pitch should lie between 15-25mm' );
SS1=round((pi*D)/(0.015))
SS2=round((pi*D)/(0.025))
disp('stator slots Ss should lie between values of SS1 to SS2' )
disp('When qs=2')
Ss2=3*P*2
disp('When qs=3')
Ss3=3*P*3
disp('When qs=4')
Ss4=3*P*4
disp('When qs=5')
Ss5=3*P*5
disp('Value of Ss for q=3 and q=4 lie between specified limits' );
disp('check: ')
Iph=(Q*1000)/(sqrt(3)*V)
disp('For star connection Iph=IL');
disp('When Ss=%d,Ss3');
Zss3=round((6*Ts)/Ss3)
SL3=Zss3*Iph
disp('When Ss=%d,Ss4');
Zss4=round((6*Ts)/Ss4)
SL4=Zss4*Iph
disp('The slot loading for q=4 that is Ss=Ss3 exceeds 500 and is rejected' );
Sc=Zss4*Ss4
Tsm=round((Zss4*Ss4)/6)
%Main dimensions of dc machine
Pa=300; %KW rating
V=500; %voltage rating
P=6; %number of poles of the machine
Bav=0.67; %specific magnetic loading (Wb/m^2)
ac=25000; %specific electrical loading (A/m)
N=500; %speed of machine
Rap=0.75; %ratio of core length to pole pitch
ns=N/60; %speed in rps
Co=pi*pi*Bav*ac*(10^-3); %output coefficient
D2L=Pa/(Co*ns); % volume= D*D*L
RLD=pi*Rap/P; %Value of L/D
D=(D2L/RLD)^(1/3) %Value of D
L=RLD*D %Value of L
va=pi*D*ns %Peripheral velocity
if va<=30
disp('Peripheral velocity is within permissible limit of 30m/s' )
else
disp('Va exceeds limit')
end
Ia=Pa*(10^3)/V %armature current
Iz=Ia/2 %current per parallel path
if Iz<200
disp('Simplex winding can be used' )
else
disp('Current limit of 200 A/path was exceeded.Simplex winding cant be used' )
Iz=Ia/P %new current per path
if Iz<200
disp('Rather LAP winding can be used' )
end
end
Ez=Bav*L*va; %average induced emf/conductor
Zc=roundn(V/Ez,1) %armature conductors per path rounded to nearest
10
disp('Total no. of conductors')
Z=Zc*P %Total no. of conductors
disp('Condn 1: Slot pitch varies from 2.5 to 3.5 cm...' )
Smin1=round(pi*D/0.035)
Smax1=round(pi*D/0.025)
disp('Condn 2: ratio of no. of slots to pole varies from 9 to 16...' )
Smin2=P*9
Smax2=P*16
disp('Final Lower & Higher bounds of no. of slots...' )
Smin=max(Smin1,Smin2)
Smax=min(Smax1,Smax2)
disp('Condn 3: for lap winding no. of slots are a multiple of pole pair(3)...' )
Smin=ceil(Smin*2/P)*P/2; %recalibrating initial value to multiple of P/2
A=Smin:3:Smax
disp('Cond 4: to reduce flux pulsations Slots/pole must be equal to integer + 0.5' )
B=A(A/P==floor(A/P)+0.5)
disp('These are the possible values of no. of slots' )
disp('Slots per pole*Rap...')
B=B*Rap/P
disp('Choosing S which is nearly an integer' )
C=B-floor(B); %finding element with min fractional part
pos=find(C==min(C)); %position of optimal value in array
S= B(pos)*P/Rap
disp('No. of conductors per slot...' )
Zps=Z/S
disp('As Zs must be an even integer, the revised no. of conductors is...' )
Zps=round(Zps);
Zrev=S*Zps
Sload=Iz*Zps
if Sload<1500
disp('Slot loading is within the permissible limit' )
end
disp('No. of commutator segments...' )
NC=Zrev/2 %no. of coils
disp('Pitch of commutator segments...' )
PC=0.7*D*pi/NC*(10^3) %pitch of coils (in mm)
if 4<PC && PC<9
disp('Pitch of coils is well within the permissible limits of(4,9)' )
end
disp('Given data : Synchronous motor , salient pole rotor' );
Q=2000 %KW
N=200 %synchronous speed
V=2000 %voltage
f=50 %frequency
D=2 %stator bore
L=0.5 %stator length
Ssp=3.5 %slots per pole per phase
Zs=4 %conductors per slot
Cp=2 %circuits per phase
kw=1.2 %winding factor
Bav=1.5 %flux density
a=2
p=3 %phase
df=3 %winding depth
ATl=2 %ratio of full load mmf to armature mmf
kf=0.8 % space factor for field winding
loss=2000 %field losses
Eph=V/1.732; % phase voltage
nop=120*f/N; % number of poles
ts=Ssp*nop*p; %total slots
z= 314*Zs; %total conductors
tph= z/6; %turns per phase
disp('since there are two parallel paths , no of turns in series perphase will be
turns per phase divided by 2' );
notsp= tph/a;% no of turns in series per phase
fpp= Eph/(4.44*f*kw*(tph/a)); % flux per pole
disp(' flux per pole : ');
disp(fpp);
tfpp= 1.2*fpp; %total flux per pole
Ap= tfpp/Bav;%area of pole body
Lp=L;% taking length of pole body equal to length of the core
disp('length of the pole : ');
disp(Lp);
Wp= Ap/Lp; %width of pole body
disp(' width of the pole : ');
disp(Wp);
Iph= (Q*1000)/(3*Eph);%current
Iz=Iph/2;%conductoe current
ATa= (2.7*Iph*tph)/nop;%armature mmf per pole
ATf1= 2*ATa;%full load field mmf
hf= ATf1/((10^4)*((2000*0.5*0.03)^(1/2))); %height of field winding
disp(' height of the field winding : ' );
disp(hf);
hop=hf+5;%height of the pole
disp(' height of the pole : ');
disp(hop);