Shaft design

Members:
Andrés Hernández Ruiz
Mario Martinez Roqueme

Diseño de máquinas I – Machine design I

Ing. Valéry José Lancheros Suárez

Universidad de Córdoba
Faculty of Engineering

Mechanical Engineering
Montería – Córdoba
2021
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES

Exercise 3-81.A countershaft carrying two V-belt pulleys is shown in the figure. Pulley A receives
power from a motor through a belt with the belt tensions shown. The power is transmitted through
the shaft and delivered to the belt on pulley B. Assume the belt tension on the loose side at B is 15
percent of the tension on the tight side.

(a) Determine the tensions in the belt on pulley B, assuming the bearings at a constant speed.
(b) Find the magnitudes of the bearing reaction forces, assuming the bearings act as simple
supports.
(c) Draw shear-force and bending-moment diagrams for the shaft. If needed, make one set for
the horizontal plane and another set for the vertical plane.
(d) At the point of maximum bending moment, determine the bending stress and the torsional
shear stress.
(e) At the point of maximum bending moment, determine the principal stresses and the
maximum shear stress.
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES

7-12 to 7-21. For the problem specified in the table, build upon the results of the original problem
to obtain a preliminary design of the shaft by performing the following tasks.

(a) Sketch a general layout, including means to locate the components and to transmit the
torque. Estimates for the component widths are acceptable at this point.
(b) Specify a suitable material for the shaft.
(c) Determine critical diameters of the shaft based on infinite fatigue life with a design factor
of 1.5. Check for yielding.
(d) Make any other dimensional decisions necessary to specify all diameters and axial
dimensions. Sketch the shaft to scale. Showing all proposed dimensions.
(e) Check the deflections at the gears, and the slopes at the gears and the bearings for
satisfaction of the recommended limits in Table 7-2. Assume the deflections for any pulleys
are not likely to be critical. If any of the deflections exceed the recommended limits, make
appropriate changes to bring them all within the limits.

In this work, the Matlab code will be attached, where the exercise is explained in more detail.

First, the unknown reactions and stresses of the exercise are sought, assuming that T1 is bigger than
T2, the equilibrium conditions applied and, we find the necessary forces.

To calculate the diameters that make up the shaft, we have two ways, the first is to assume initial
diameters, determine the critical point and calculate a critical diameter with said initial diameter,
the idea is that the assumed initial diameter, is greater than the critical diameter, but not so larger
to avoid oversize of the shaft.

Another way is to assume one or more of the Marin factors to find a critical diameter and propose
a diameter similar to Example 7-3 in Shigley's mechanical engineering design book.

We assume a diameter since it is the fastest method and propose that 𝑑1 = 𝑑5 = 1,45

𝑑3 = 1,8
𝑑2 = 𝑑4 = 1,7

In a economic view, we will use a cheap Steel, Aisi 1020 CD, their properties are:
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES

As we noted from the total moment graph, the maximum moment occurs at B where 𝑀 = 2111 𝑙𝑏 ∗
𝑖𝑛 𝑦 𝑇 = 1000 𝑙𝑏
As there is no alternating torque Ta = 0, therefore Tm = 1000 lb, but if the alternating moment is
presented, thus Ma = 2111 lb * in, which is completely reversible, therefore 𝑀𝑚 = 0.

Using table 7-1 of the guide book, for a milled, keyway we have that Kt = 2.14 and Kts = 3 .in addition,
we have the ratio of the radius r in the lower part of the groove and the diameter of the axis d, is r
/ d = 0.02, also taken from Table 7-1.

Therefore r = 0.02 * d2, we will use d2 which is to the right of point B since it is the area near this
point with the highest bending moment.

Thus, 𝑟 = 0.034

We proceed to calculate the fatigue concentration factors we must calculate the notch sensitivity
factors.

From equation 6-35a and 6-35b we have that

Now using equation 6-32

Calculamos los factores de concentración por fatiga son:

𝐾𝑓 = 1,9006 𝑦 𝐾𝑓𝑠 = 2,4308

We proceed to calculate the Marín factors:

𝐾𝑎 = 𝑎𝑆𝑢𝑡 𝑏 where 𝑎 = 2, 𝑏 = −0,217 thus 𝐾𝑎 = 0,8826

𝐾𝑑 = 0,8306 because 0,11 < 𝑑 < 2

𝐾𝑐 = 1
𝐾𝑑 = 𝐾𝑓 = 𝐾𝑒 = 1
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES
Now, we proceed to calculate the diameter d2, using equation

7-8, with n = 1.5

We have d = 1.436 which is less than 1.7, so it is a good result.

We verify that the axis has not exceeded the yield point, for this we use Von Misses from equation
7-15.

Therefore, 𝜎´𝑚𝑎𝑥 =14588 psi

With which we obtain a fluency safety factor of:

𝑆𝑦
𝑛𝑦 = = 3,9 so the material does not flow.
𝜎´𝑚𝑎𝑥

We proceed to present the design proposal.

Finally, we look for the deflections in the bearings and gears using the Engilab Beam 2D program
with the previously mentioned data. The program shows us the following results:

For the XY plane

 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES
EngiLab Beam.2D 2018 Analysis Report
Model: "Model1.bea"
Report created: Tuesday, 27 April 2021, 12:48:00

Node Displacements

Node ID X-Displacement Y-Displacement Z-Rotation

1 0.000000000 0.000000000 0.0009068201
2 0.000000000 0.000339723 0.0009041471
3 0.000000000 0.005527839 0.0002616980
4 0.000000000 0.005732950 0.0000694731
5 0.000000000 0.004604634 -0.0004336501
6 0.000000000 0.004001191 -0.0005288644
7 0.000000000 0.000276011 -0.0007350122
8 0.000000000 0.000000000 -0.0007365396
Sign convention: Positive according to Global Axes
Displacements reported with 9 decimal places
Rotations (in Radians) reported with 10 decimal places
Background color: Free DOF Constrained DOF Spring DOF

Element End Forces

Element ID Axial Force Shear Force Bending Moment

1 Start 0.0000 -222.7273 0.000
1 End 0.0000 -222.7273 -83.523
2 Start 0.0000 -222.7273 -83.523
2 End 0.0000 -222.7273 -1781.818
3 Start 0.0000 127.2727 -1781.818
3 End 0.0000 127.2727 -1622.727
4 Start 0.0000 127.2727 -1622.727
4 End 0.0000 127.2727 -922.727
5 Start 0.0000 127.2727 -922.727
5 End 0.0000 127.2727 -763.636
6 Start 0.0000 127.2727 -763.636
6 End 0.0000 127.2727 -47.727
7 Start 0.0000 127.2727 -47.727
7 End 0.0000 127.2727 0.000
Sign convention: According to the diagrams' sign convention
Axial force: Positive when member is in tension
Shear force: Positive when it rotates the member clockwise
Bending moment: Positive when the "bottom" fiber is in tension ("top" fiber in compression)
Forces reported with 4 decimal places
Bending moments reported with 3 decimal places
Background color: Normal DOF Hinge DOF

Support Reactions

Node ID X-Force Y-Force Z-Moment

1 -222.7273
8 -127.2727
Sign convention: Positive according to Global Axes
Forces reported with 4 decimal places
Bending moments reported with 3 decimal places
Background color: Free DOF Constrained DOF Spring DOF
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES

For the XZ plane

EngiLab Beam.2D 2018 Analysis Report
Model: "XZ.bea"
Report created: Tuesday, 27 April 2021, 12:55:19

Node Displacements

Node ID X-Displacement Y-Displacement Z-Rotation

1 0.000000000 0.000000000 0.0015674992
2 0.000000000 0.000587628 0.0015660231
3 0.000000000 0.010748566 0.0008957138
4 0.000000000 0.011730167 0.0006693799
5 0.000000000 0.011892537 -0.0007161775
6 0.000000000 0.010748566 -0.0011196423
7 0.000000000 0.000797253 -0.0021233841
8 0.000000000 0.000000000 -0.0021273204
Sign convention: Positive according to Global Axes
Displacements reported with 9 decimal places
Rotations (in Radians) reported with 10 decimal places
Background color: Free DOF Constrained DOF Spring DOF

Element End Forces

Element ID Axial Force Shear Force Bending Moment

1 Start 0.0000 -122.9947 0.000
1 End 0.0000 -122.9947 -46.123
2 Start 0.0000 -122.9947 -46.123
2 End 0.0000 -122.9947 -983.957
3 Start 0.0000 -122.9947 -983.957
3 End 0.0000 -122.9947 -1137.701
4 Start 0.0000 -122.9947 -1137.701
4 End 0.0000 -122.9947 -1814.171
5 Start 0.0000 -122.9947 -1814.171
5 End 0.0000 -122.9947 -1967.914
6 Start 0.0000 327.9857 -1967.914
6 End 0.0000 327.9857 -122.995
7 Start 0.0000 327.9857 -122.995
7 End 0.0000 327.9857 0.000
Sign convention: According to the diagrams' sign convention
Axial force: Positive when member is in tension
Shear force: Positive when it rotates the member clockwise
Bending moment: Positive when the "bottom" fiber is in tension ("top" fiber in compression)
Forces reported with 4 decimal places
Bending moments reported with 3 decimal places
Background color: Normal DOF Hinge DOF

Support Reactions

Node ID X-Force Y-Force Z-Moment

1 -122.9947
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES
8 -327.9857
Sign convention: Positive according to Global Axes
Forces reported with 4 decimal places
Bending moments reported with 3 decimal places
Background color: Free DOF Constrained DOF Spring DOF
 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES

From what has been presented, we will note that the deflections in the bearings O and C and in the gears A
and B are respectively:
𝛿𝑂 = 0 𝑖𝑛
𝛿𝐶 = 0 𝑖𝑛
𝛿𝐴 = 0,01204 𝑖𝑛
𝛿𝐵 = 0,01149 𝑖𝑛

In the Engilab Beam 2D files, we will note that the deflections are within the range proposed in table 7-2.

Below is the Matlab code.

 UNIVERSIDAD DE CÓRDOBA
FACULTAD DE INGENIERÍAS
INGENIERÍA MECÁNICA
DISEÑO DE MAQUINAS I – DISEÑO DE EJES
 Matlab Code
Andrés Hernández Ruiz
Mario Martinez Roqueme
April 27, 2021

1 clc
2 % Determinacion de las tensiones desconocidas , supondremos que T1 > T2 por lo
3 % que T2 =0 ,15 T1 segun las condiciones del problema
4 % Datos
5
6 % Determination of the unknown stresses , we will assume that T1 > T2 so that T2 =0.15 T1 according
to problem conditions
7 % Data
8 d_A =8; % Diametro de la polea A Pulley Diameter A
9 d_B =6; % Diametro de la polea B Pulley Diameter B
10 d_s =1; % Diametro del eje Shaft diameter
11 % Sumatoria de momentos en el eje x Sum of moments in the x - axis
12 TA =1000;
13 T1 = TA /(0.85* d_B /2) ;
14 T2 =0.15* T1 ;
15 % Reacciones Reactions
16 Cy = -350*8/22;
17 Cz = -1.15*16* T1 /22;
18 Oz = - Cz -1.15* T1 ;
19 Oy = -350 - Cy ;
20
21
22 % Diagramas de fuerza cortante y momento flector en el pano xy
23 % Desde O hasta A
24 % Shear force and bending moment diagrams in the xy plane
25 % From O to A
26 x_OA =0:1:8;
27 Vy_xy_OA = Oy ;
28 Mz_xy_OA = Oy * x_OA ;
29 % Desde A hasta C
30 % From A to C
31 x_AC =8:1:22;
32 Vy_xy_AC = Oy +350;
33 Mz_xy_AC =350*( x_AC -8) + Oy * x_AC ;
34 % Graficos de las ecuaciones
35 % Plano xy
36 % Fuerza cortante
37
38 % Equation graphs
39 % Plane xy
40 % Shear force
41 XOA_xy =[0 8];
42 VOA_xy =[ Vy_xy_OA Vy_xy_OA ];
43 XAC_xy =[8 22];
44 VAC_xy =[ Vy_xy_AC Vy_xy_AC ];
45 subplot (3 ,2 ,1) ;
46 plot ( XOA_xy , VOA_xy , ’b ’) ; hold on
47 plot ( XAC_xy , VAC_xy , ’b ’) ;
48 grid on
49 grid minor
50 xlabel ( ’x [ in ] ’) ;
51 ylabel ( ’V [ lb ] ’) ;
52 title ( ’ Fuerza cortante plano xy Shear force plane xy ’) ;
53 % Momento flector Bending moment
54 subplot (3 ,2 ,2) ;
55 plot ( x_OA , Mz_xy_OA , ’r ’) ; hold on
56 plot ( x_AC , Mz_xy_AC , ’r ’) ;
57 grid on
58 grid minor
59 xlabel ( ’x [ in ] ’) ;
60 ylabel ( ’M [ lb * in ] ’) ;

1
61 title ( ’ Momento flector plano xy Bending moment plane xy ’) ;
62
63
64 % Diagramas de fuerza cortante y momento flector en el pano xz
65 % Desde O hasta B
66 % Shear force and bending moment diagrams in the plane xz
67 % From O to B
68
69 x_OB =0:1:16;
70 Vz_xz_OB = - Oz ;
71 My_xz_OB = - Oz * x_OB ;
72 % Desde B hasta C From B to C
73 x_BC =16:1:22;
74 Vz_xz_BC = - Oz +1.15* T1 ;
75 My_xz_BC = - Oz * x_BC -1.15* T1 *( x_BC -16) ;
76 % Graficos de las ecuaciones
77 % Plano xz
78 % Fuerza cortante
79
80 % Equation graphs
81 % Plane xy
82 % Shear force
83 XOB_xz =[0 8];
84 VOB_xz =[ Vz_xz_OB Vz_xz_OB ];
85 XBC_xz =[8 22];
86 VBC_xz =[ Vz_xz_BC Vz_xz_BC ];
87 subplot (3 ,2 ,3) ;
88 plot ( XOB_xz , VOB_xz , ’b ’) ; hold on
89 plot ( XBC_xz , VBC_xz , ’b ’) ;
90 grid on
91 grid minor
92 xlabel ( ’x [ in ] ’) ;
93 ylabel ( ’V [ lb ] ’) ;
94 title ( ’ Fuerza cortante plano xz Shear force plane xz ’) ;
95 % Momento flector Bending moment
96 subplot (3 ,2 ,4) ;
97 plot ( x_OB , My_xz_OB , ’r ’) ; hold on
98 plot ( x_BC , My_xz_BC , ’r ’) ;
99 grid on
100 grid minor
101 xlabel ( ’x [ in ] ’) ;
102 ylabel ( ’M [ lb * in ] ’) ;
103 title ( ’ Momento flector plano xz Bending moment plane xz ’) ;
104
105 % Procedemos a graficar el momento total
106 % We proceed to plot the total momentum
107 L_OA =[0 8];
108 L_AB =[8 16];
109 L_BC =[16 22];
110 MO =0;
111 MA = sqrt (( Mz_xy_AC (1) ) ^2+( My_xz_OB (9) ) ^2) ;
112 MB = sqrt (( My_xz_OB ( end ) ) ^2+( Mz_xy_AC (9) ) ^2) ;
113 MC = sqrt (( Mz_xy_AC ( end ) ) ^2+( My_xz_BC ( end ) ) ^2) ;
114 MAT =[ MO MA ];
115 MBT =[ MA MB ];
116 MCT =[ MB MC ];
117 subplot (3 ,2 ,5) ;
118 plot ( L_OA , MAT , ’g ’) ; hold on
119 plot ( L_AB , MBT , ’g ’) ; hold on
120 plot ( L_BC , MCT , ’g ’) ;
121 grid minor ;
122 title ( ’ Momento total ’) ;
123 xlabel ( ’x [ in ] ’) ;
124 ylabel ( ’M [ lb * in ] ’) ;
125
126 % Para conocer el momento flector maximo , vemos la grafica , el cual es M =2111 lb * in
127 % To know the maximum bending moment , we see the graph , which is M =2111 lb * in .
128 M =2111;
129 Sig =32* M /( pi * d_s ^3) ;
130 Tau =16* TA /( pi * d_s ^3) ;
131
132 % Procedemos a hallar los esfuerzos principales con las ecuaciones del
133 % circulo de Mohr
134
135 % We proceed to find the principal stresses with the equations of the Mohr ’ s circle
136

2
137 Sig_1 = Sig /2+ sqrt (( Sig /2) ^2+( Tau ) ^2) ;
138 Sig_2 = Sig /2 - sqrt (( Sig /2) ^2+( Tau ) ^2) ;
139 Tau_max = sqrt (( Sig /2) ^2+( Tau ) ^2) ;
140
141
142 % Suponiendo los diametros que conforman el eje tenemos
143 % Assuming the diameters that make up the axis we have
144
145 % d1 = d5 =1 ,45
146 % d3 =1 ,8
147 % d2 = d4 =1 ,7
148 d1 =1.45; % in
149 d5 = d1 ; % in
150 d3 =1.8; % in
151 d2 =1.7; % in
152 d4 = d2 ; % in
153 % Como vemos en la grafica , el momento maximo se da en B , por lo que tenemos
154 % un momoento alternante pero un torque medio
155 % As we can see in the graph , the maximum moment occurs at B , so we have an
156 % alternating moment and a mean torque equals to :
157 Ma = M ;
158 Tm = TA ;
159 Ta =0;
160 Mm =0;
161 % El acero a usar es acero poco caro The steel to be used is a cheap steel
162 Sy =57; % kpsi
163 S_ut =68; % kpsi
164 E =27*10^6; % psi
165 % De la tabla 7 -1 tenemos para una chavetera fresada con Kt =2.14 y Kts =3 y que r =0 ,02 d , donde d
= d2 debido a que alli se
166 % presenta el mayor momento
167 Kt =2.14;
168 Kts =3;
169 d = d2 ;
170 r =0.02* d2 ;
171 % Procedemos a calcular los factores de sensibilidad a la muesca
172 % Now we will calculate the notch sensitive factors
173 raiz_a =0.246 -3.08*(10^ -3) * S_ut +1.51*(10^ -5) * S_ut ^2 -2.67*(10^ -8) * S_ut ^3;
174 raiz_as =0.19 -2.51*(10^ -3) * S_ut +1.35*(10^ -5) * S_ut ^2 -2.67*(10^ -8) * S_ut ^3;
175 q =1/(1+( raiz_a / sqrt ( r ) ) ) ;
176 qs =1/(1+( raiz_as / sqrt ( r ) ) ) ;
177 % Ahora calculando los concent radores de esfuerzo por fatiga
178 % Now we calculate the
179 Kf =1+ q *( Kt -1) ;
180 Kfs =1+ qs *( Kts -1) ;
181 % Ya teniendo los con centrado res de fatiga calculamos el limite de
182 % resistencia a la fatiga
183 % Buscamos los factores de Marin
184 % We determine the Marin factors
185 if S_ut <=200
186 S_exp =0.5* S_ut ;
187 elseif S_ut >200
188 S_exp =100;
189 end
190 % Procedemos a buscar los factores de Marin
191 % We determine the Marin factors
192 disp ( ’ Elija un acabado superficial para determinar ka ’) ;
193 disp ( ’1 - - - Esmerilado ( Ground ) ’) ;
194 disp ( ’2 - - - Maquinado o laminado en frio ( Machined or cold drawn ) ’) ;
195 disp ( ’3 - - - Laminado en caliente ( Hot rolled ) ’) ;
196 disp ( ’4 - - - Como sale de la forja ( As forged ) ’) ;
197 op = input ( ’ Digite su eleccion ( Choose your option ) = ’) ;
198 switch ( op )
199 case 1
200 a =1.21;
201 b = -0.067;
202 ka = a *( S_ut ) ^( b ) ;
203 case 2
204 a =2.0;
205 b = -0.217;
206 ka = a *( S_ut ) ^( b ) ;
207 case 3
208 a =11;
209 b = -0.65;
210 ka = a *( S_ut ) ^( b ) ;
211 case 4

3
212 a =12.7;
213 b = -0.758;
214 ka = a *( S_ut ) ^( b ) ;
215 otherwise
216 disp ( ’ No es una opcion ( Is not an option ) ’) ;
217 end
218 disp ( ’ El elemento esta girando ? esto para determinar kb ( The element is rotating ? this is to
determine kb ) ’) ;
219 disp ( ’1 - - - Si esta girando , ( is spinning ) ’) ;
220 disp ( ’2 - - - No esta girando ,( is not spinning ) ’) ;
221 disp ( ’3 - - - Solo carga axial ’) ;
222 op2 = input ( ’ Digite su opcion ( Choose your option ) = ’) ;
223 switch op2
224 case 1
225 if d >=0.11 && d <=2
226 kb =( d /0.3) ^( -0.107) ;
227 elseif d >2 && d <10
228 kb =0.91* d ^( -0.157) ;
229 end
230 case 2
231 de =0.37* d ;
232 if de >=0.11 && de <=2
233 kb =( de /0.3) ^( -0.107) ;
234 elseif de >2 && de <10
235 kb =0.91*( de ) ^( -0.157) ;
236 end
237 case 3
238 kb =1;
239 end
240 disp ( ’ El elemento esta en tension , carga axial o torsion ? , esto para determinar kc ( Is the
element in tension , axial load or torsion ? this is to determine kc ) ’) ;
241 disp ( ’1 - - - Flexion ’) ;
242 disp ( ’2 - - - Axial ’) ;
243 disp ( ’3 - - - Torsion ’) ;
244 disp ( ’4 - - - Carga combinada , ( combined load ) ’) ;
245 op3 = input ( ’ Escoja la opcion ( Choose an option ) = ’) ;
246 switch op3
247 case 1
248 kc =1;
249 case 2
250 kc =0.85;
251 case 3
252 kc =0.59;
253 case 4
254 kc =1;
255 end
256 % Como desconocemos datos de temperatura , confiabilida y efectos varios
257 % Since we do not know data on temperature , reliability and various effects
258 kd =1;
259 ke =1;
260 kf =1;
261 S_e = ka * kb * kc * kd * ke * kf * S_exp ;
262
263 % Procedemos a calcular el diametro con un factor de seguridad n =1 ,5
264 % Now we can calculate the diameter with a factor of safety n =1 ,5
265 n =1.5;
266 A = sqrt (4*( Kf * Ma ) ^2+3*( Kfs * Ta ) ^2) ;
267 B = sqrt (4*( Kf * Mm ) ^2+3*( Kfs * Tm ) ^2) ;
268 d_cr =((16* n / pi ) *( A /( S_e *1000) + B /( S_ut *1000) ) ) ^(1/3) ;
269 if d_cr < d2
270 fprintf ( ’ El diametro critico del eje es = %1.2 f \ n ’ , d_cr ) ;
271 else
272 disp ( ’ Proponga otro diametro ’) ;
273 end
274
275 % Ahora bien , calculamos el esfuerzo maximo con Von Misses ecuacion 7 -15 para verificar
276 % que no haya falla debido a la fluencia
277
278 % Now , we calculate the maximum stress with Von Misses equation 7 -15 to verify that there is no
failure due to yelding
279 Sigma_pmax = sqrt ((32* Kf *( Mm + Ma ) /( pi *( d_cr ^3) ) ) ^2+3*(16* Kfs *( Ta + Tm ) /(( pi * d_cr ^3) ) ) ^2) ;
280 % Por lo que usando la ecuacion 7 -16
281 % By using equation 7 -16
282 ny = Sy *1000/ Sigma_pmax ;
283 if ny >1
284 fprintf ( ’ El factor de seguridad a fluencia es %1.2 f por lo tanto se encuentra en regimen

4
 elastico \ n ’ , ny ) ;
285 else
286 disp ( ’ El material se deforma plasticamente ’) ;
287 end
288
289 % Se procede a calcular las deflexiones en el programa Engilab Beam 2D , el
290 % cual se adjunta en el informe
291
292 % The deflections are calculated in the Engilab Beam 2 D program , which is attached in the
report .
293
294 delta_O =0;
295 delta_C =0;
296 delta_A = sqrt ( 0 . 0 0 5 5 2 5 3 6 3 ^ 2 + 0 . 0 1 0 7 0 0 5 3 1 ^ 2 ) ;
297 delta_B = sqrt ( 0 . 0 1 0 7 6 2 8 6 0 ^ 2 + 0 . 0 0 4 0 1 1 1 5 7 ^ 2 ) ;
298 fprintf ( ’ Las deflexiones en los cojinetes O y C son %1.5 f in y %1.5 f in r es p ec ti va m en te \ n (
The deflections in bearings O and C are --- in and --- in respectively .) ’ , delta_O , delta_C )
;
299 fprintf ( ’ Las deflexiones en los engranes A y B son %1.5 f in y %1.5 f in r es pe c ti va me n te \ n ( The
deflections in gears A and B are --- in and --- in respectively .) ’ , delta_A , delta_B ) ;
300 % En el plano xy calcularemos la pendiente en O mediante la desviacion
301 % tangencial
302
303 % In the xy plane we will calculate the slope in O by means of the deviation tangential
304 I_d1 = pi *( d1 ^4) /64;
305 I_d2 = pi *( d2 ^4) /64;
306 I_d3 = pi *( d3 ^4) /64;
307 syms x
308 Mz_OA = Oy * x ;
309 Mz_AC =350*( x -8) + Oy * x ;
310 ThetaOC =(1/( E * I_d1 ) ) * int ( Mz_OA ,x ,0 ,0.375) +(1/( E * I_d2 ) ) * int ( Mz_OA ,x ,0.375 ,8) +(1/( E * I_d2 ) ) * int (
Mz_AC ,x ,8 ,9.25) +(1/( E * I_d3 ) ) * int ( Mz_AC ,x ,9.25 ,14.75) +(1/( E * I_d2 ) ) * int ( Mz_AC ,x
,14.75 ,21.625) +(1/( E * I_d1 ) ) * int ( Mz_AC ,x ,21.625 ,22) ;
311 ThetaOC_deg = ThetaOC *(180/ pi ) ;
312 fprintf ( ’ La pendiente de O con respecto a C es %1.4 f grados o %1.4 f radianes ( The slope of O
with respect to C is --- degrees or --- radians .) ’ , ThetaOC_deg , ThetaOC ) ;

5
Figure 1: Shear force and bending moment

Figure 2: shaft