fprintf('\n--------------Calculation for Bevel Gears Design--------------\n');
Zp = 24;
disp(['Zg(No. of Pinion Teeth) = ',num2str(Zp)]);
disp('-------------------------------------------------')
Zg = 32;
disp(['Zg(No. of Gear Teeth) = ',num2str(Zg)]);
disp('-------------------------------------------------')
Power = 12.5;
disp(['Power = ',num2str(Power),' KW']);
disp('-------------------------------------------------')
omg = 1440;
disp(['Angular velocity = ',num2str(omg),' rpm']);
disp('-------------------------------------------------')
Cs = 1.5;
disp(['Rated Torque rate = ',num2str(Cs),' ']);
Sut = 750;
disp(['Ultimate Strength = ',num2str(Sut),' N/mm2']);
%int m;
%b = 10*m;
v = 7.5;
disp([' V = ',num2str(v),' m/s ']);
FOS1 = 2;
disp([' FOS (initially assumed for "m" calculation = ',num2str(FOS1),' ']);
% b/A = 1/3
r = 2/3;
Cv1= 5.6/(5.6+sqrt(v));
disp([' Cv = ',num2str(Cv1),' ']);
Sigma_b = Sut/3;
fprintf('Calculation of Torque Mt\n');
Mt = 60*(10^6)*Power/(2*pi*omg);
disp([' Mt = ',num2str(Mt),' N-mm']);
disp('-------------------------------------------------')
gamma = rad2deg(atan(Zp/Zg));
disp([' gamma = ',num2str(gamma),' deg']);
fprintf('\n As both the gear and pinion material is same, then obviously pinion is weaker,
\n');
Zp1 = Zp/cos(deg2rad(gamma));
disp([' Zp1 = ',num2str(Zp1),' ']);
disp('-------------------------------------------------')
disp("Reading the Lewis Form Factor table :");
a = readmatrix('yfactor.txt')
for i = 1:1:34
z1 = a(i,end-1);
if z1 == Zp1
Y = a(i,end);
elseif Zp1 > z1 && Zp1 < a(i+1,end-1)
Y = ((a(i+1,end)-a(i,end))*(Zp1 - a(i,end-1))/(a(i+1,end-1) - z1) + a(i,end));
break;
else
continue;
end
end
disp(['Lewis Form Factor(Y) = ',num2str(Y)]);
disp('-------------------------------------------------')
%A = 3*b;
m = ((2*FOS1*Cs*Mt)/(Zp*10*Sigma_b*Cv1*Y*r))^(1/3);
disp(['Module(Calculated) = ',num2str(m),' mm ']);
% first preference value of module
for j = 0:0.25:10
if m == j
break;
elseif m > j && m < j+0.25
m = j+0.25;
else
continue;
end
end
disp(['Preferential Module = ',num2str(m)]);
disp('-------------------------------------------------')
disp("New Dimensions of Gear : ");
disp('-------------------------------------------------')
Dp = m*Zp;
disp([' New Dimension: Dp = ',num2str(Dp),' mm ']);
Dg = m*Zg;
disp([' New Dimension: Dg = ',num2str(Dg),' mm ']);
C = 11400;
disp([' C = ',num2str(C),' ']);
disp('-------------------------------------------------')
%New values of velocity and it's factors
disp("calculating new speed and Velocity factor");
disp('-------------------------------------------------')
v = pi*Dp*omg/(60*1000);
disp([' V = ',num2str(v),' m/s ']);
Cv = 5.6/(5.6+sqrt(v));
disp([' Cv = ',num2str(Cv),' ']);
Pt = 2*Mt/Dp;
disp([' Pt = ',num2str(Pt),' N']);
disp('-------------------------------------------------')
disp("Reading the bevel gear grades table :");
g = readmatrix("bevelgeargrades.txt") %importing the grade table
for k = 1:1:7
if(m <= g(k,1))
e = g(k,end);
break;
elseif (m > g(k,1) && m<g(k+1,1))
e = g(k+1,end);
break;
else
continue;
end
end
disp([' error = ',num2str(e),' mm ']);
b = 10*m;
disp([' New Value of Face Width = ',num2str(b),' mm ']);
disp('-------------------------------------------------')
A = sqrt((Dp/2)^2 + (Dg/2)^2);
disp([' Cone Distance (A) = ',num2str(A),' mm ']);
disp('-------------------------------------------------')
Pd = (21*v*(C*e*0.001*b + Pt)/(21*v + sqrt(C*e*0.001*b + Pt)));
disp([' Dynamic load (Pd) = ',num2str(Pd),' N']);
disp('-------------------------------------------------')
Peff = Cs*Pt + Pd;
disp([' Peff = ',num2str(Peff),' N']);
disp('-------------------------------------------------')
Sb = m*b*Sigma_b*Y*(1-(b/A));
disp([' S_b(New Modified Value ) = ',num2str(Sb),' N']);
disp('-------------------------------------------------')
FOS = Sb/Peff;
disp([' FOS (For Dynamic Load ) = ',num2str(FOS)]);
disp('-------------------------------------------------')
disp("Calculation of brinell Hardness Number");
disp('-------------------------------------------------')
FOS_w = 2;
disp("FOS (For wear assumed ) = 2 ");
Sw = Peff * FOS_w;
disp([' Sw = ',num2str(Sw),' N ']);
Q = 2*Zg/(Zg + Zp*tan(deg2rad(gamma)));
BHN = sqrt((Sw * cos(deg2rad(gamma)) * 100*100)/(0.75*b*Q*Dp*0.16));
disp([' BHN = ',num2str(BHN),' ']);
