EJERCICIO 1
Integrantes:
David Lema
Luis Pruna
Cristian Tutillo
Nivel: Sexto.
Materia: Diseño de Elementos de Máquinas.
Docente: Ing. Félix J. Manjarrés A. M. Sc.
Nrc: 3168.
Tema: Ejes
Diseñe el eje de salida III de un reductor de engranajes helicoidales de doble reducción de la
Figura 1. La potencia para transmitir por el eje de salida es P3 = 9,5 kW. El eje de salida gira a
una velocidad de rotación de w3 = 120 rpm. El engranaje helicoidal en el eje de salida tiene el
número de dientes de z = 108 con módulo mn = 3 mm, un ángulo de hélice de β = 12,4 ° y un
ancho de cara de B = 50 mm. El factor de seguridad se selecciona como 1,5. En la Figura 1, los
detalles de las dimensiones relevantes son A = 15 mm, = 40 mm, C = 20 mm, B2 = 50 mm y S =
8 mm.
Librerías a utilizar
In [10]: import os
os.system('cls')
import numpy as np
import matplotlib.pyplot as plt
import sympy
import sympy as sp
import pandas as pd
import math
from prettytable import PrettyTable
from matplotlib import rc
from IPython.display import Image
sympy.init_printing(use_unicode=True)
rc("mathtext", default = "regular")
#Colores
BLACK = '\033[30m'
RED = '\033[31m'
GREEN = '\033[32m'
YELLOW = '\033[33m'
BLUE = '\033[34m'
MAGENTA = '\033[35m'
CYAN = '\033[36m'
WHITE = '\033[37m'
RESET = '\033[39m'
Párametros de entrada
In [11]: σ_m1 = 172 # [MPa]
τ = 35 # [MPa]
w_3 = 120 # [rpm]
z = 108 # Número de dientes
β = 12.4 # [°]
α = 20 # [°]
# Datos adicionales
ka = 1.3
m = 3
pi = math.pi
cos = math.cos
tan = math.tan
SOLUCIÓN DEL EJERCICIO
1. Selección del material y el tratamiento térmico:
Material: Acero de medio carbono
Tratamiento térmico: Templado
Esfuerzo permisible: [σ_-1] = 172 MPa
2. Se determina el diámetro mínimo del eje:
Esfuerzo permisible:[20 ~ 40]
Se selecciona el valor de τ=35 Mpa
El mínimo diámetro aceptable puede se puede determinar:
6
9.55 ∗ 10 P
d ⩾ √ ∗
0, 2[τ ] n
3. Definir la disposición de los elementos montados en el eje:
Planteamiento de la posible disposición:
F igura 2
4. Determinar el diámetro para cada escalón del eje, considerando la caracteristica
fundamental de cada componente montado.
Se recomienda iniciar por el acoplador. La longitud y diámetro de este escalón debe ser
compatible con dicho acoplador.
El torque transmitido por el eje va a ser:
Con dicho parámetro del catálogo de acopladores, se selecciona:
LX3 → GB/T 5014 2003
El par nominal es: 1250000 Nm
Diámetro interior → d1 = 48 mm
Longitud total acoplamiento → L = 112 mm
Longitud sujeción → L_s = 84 mm
→ Para el escalón I:
dI = 48mm
LI = 82mm
→ Para el escalón II: Considerar un hombro de soporte para el acoplador de altura h =
(0.07 ~ 0.1)d
dI I = 58mm
LI I = 50mm
Se selecciona un tipo de rodamiento idóneo para la aplicación:
Radiales → Rodamientos de bolas
Axiales → Rodamientos de rodillos
Selección: 7313 C (tentativo)
dI I I = 65mm
LI I I = 33mm
→ Para el escalón VI (Se salta a este escalón por considerar todoa aquellos escalones que
poseen elementos montados sobre sí) se considera un cubo del engranaje helicoidal → B_2
= 50 mm, por tanto se deja en voladizo 2 mm para que contacte con el separador entre
engrane y cojintete:
LV I = 48mm
El engrane helicoidal seleccionado para este caso tiene un diámetro de cubo de 70 mmm.
→ Para el escalón V se utiliza lo denominado ANILLO DE HOMBRO
h > 0.07 ∗ d → Condici ón mínima
h = 5 mm
dV = 80 mm
b ≥ 1.4 ∗ h
LV = 10 mm
→ Para el segmento del eje IV Cuando se trata de rodamienros → Catálogo de los mismos,
se especifica "h"→ 6 mm
dI V = 77 mm
LI V = S + S + B1 + C − LV = 73 mm
→ Para el segmento del eje VII
Se considera:
dV I I = 65mm
LV I I = T + S + A + B2 − Lv
LV I I = 33 + 8 + 15 + 50 − 48 = 58 mm
In [13]: #Resultados obtenidos
di = 48 # [mm]
Li = 83 # [mm]
dii= 58 # [mm]
Lii= 50 # [mm]
diii= 65 # [mm]
Liii= 33 # [mm]
dvi = 48 # [mm]
Lvi = 70 # [mm]
dv = 80 # [mm]
Lv = 10 # [mm]
div = 77 # [mm]
Liv = 73 # [mm]
dvii = 77 # [mm]
Lvii = 73 # [mm]
5. Determinar elementos circunferenciales.
Constituyen ser los de movimiento circunferencial en el eje ---- caso chavetas
Se comprenden las cotas de la siguiente forma:
bhl = 201240 --- ENGRANE -- AJUSTE H7N6
bhl = 14970 --- ACOPLADOR-- AJUSTE H7K6
6. Fileteos, chaflanes y redondeos
Esta dimensión se escoge por criterio de diseño en base a la proporción de tamaño del eje y de los
componentes: -- 2mm
7. Analisis de fuerza
7.1.- Análisis de fuerza en el engrane (Fr,Fa,Ft)
m ∗ z
dz =
cos(β)
2 ∗ T
Ft =
dz
Ft ∗ tan(α)
Fr =
cos(β)
Fa = Ft ∗ tan(β)
Definición de vectores
In [14]: d = []
T = []
Tca = []
Ft = []
Fr = []
Fa = []
σca = []
τf = []
σa = []
Sσ = []
Sτ = []
Sca = []
Fs1 = []
In [15]: dz = ((m*z)/(cos(math.radians(β))))
dz = round (dz,2)
print("dz= ","{:.2f}".format(dz),'mm')
dz= 331.74 mm
Párametro de variación
In [16]: P_3 = np.arange(9.5,15.5,0.5) # Potencia [KW]
Bucle de solución
In [17]: for i in range(0, len(P_3)):
d.append((((9.55*(10**6))*(P_3[i]/w_3))/(0.2*τ))**(1/3))
T.append((9.55*10**6) * (P_3[i]/w_3))
Tca.append(T[i]*ka)
Ft.append(((2*T[i])/(dz)))
Fr.append(((Ft[i]*tan(math.radians(α)))/(cos(math.radians(β)))))
Fa.append((Ft[i]*tan(math.radians(β))))
In [18]: print(" Tabla de Iteración 1")
table= PrettyTable([RED+"P_3[KW]",BLACK+"d[mm]","T[N.mm]","Tca[N.mm]","Ft[N]","Fr[N]
for i in range(0,len(P_3)):
table.add_row(([RED+"{:.1f}".format(P_3[i]),BLACK+"{:.1f}".format(d[i]),round
round(Fr[i],2),round(Fa[i],2)]))
print(table)
Tabla de Iteración 1
+---------+-------+------------+------------+---------+---------+---------+
| P_3[KW] | d[mm] | T[N.mm] | Tca[N.mm] | Ft[N] | Fr[N] | Fa[N] |
+---------+-------+------------+------------+---------+---------+---------+
| 9.5 | 47.6 | 756041.67 | 982854.17 | 4558.04 | 1698.61 | 1002.15 |
| 10.0 | 48.4 | 795833.33 | 1034583.33 | 4797.93 | 1788.02 | 1054.89 |
| 10.5 | 49.2 | 835625.0 | 1086312.5 | 5037.83 | 1877.42 | 1107.64 |
| 11.0 | 50.0 | 875416.67 | 1138041.67 | 5277.73 | 1966.82 | 1160.38 |
| 11.5 | 50.8 | 915208.33 | 1189770.83 | 5517.62 | 2056.22 | 1213.13 |
| 12.0 | 51.5 | 955000.0 | 1241500.0 | 5757.52 | 2145.62 | 1265.87 |
| 12.5 | 52.2 | 994791.67 | 1293229.17 | 5997.42 | 2235.02 | 1318.62 |
| 13.0 | 52.9 | 1034583.33 | 1344958.33 | 6237.31 | 2324.42 | 1371.36 |
| 13.5 | 53.5 | 1074375.0 | 1396687.5 | 6477.21 | 2413.82 | 1424.11 |
| 14.0 | 54.2 | 1114166.67 | 1448416.67 | 6717.11 | 2503.22 | 1476.85 |
| 14.5 | 54.8 | 1153958.33 | 1500145.83 | 6957.0 | 2592.62 | 1529.6 |
| 15.0 | 55.5 | 1193750.0 | 1551875.0 | 7196.9 | 2682.02 | 1582.34 |
+---------+-------+------------+------------+---------+---------+---------+
8. Fuerzas en el eje; reacciones:
Fuerzas horizontales -> respecto a "x"
RH 1 = 1461N
RH 2 = 3097N
Fuerzas verticales -> respecto a "y"
RV 1 = 1539N
RV 2 = 160N
Decide como vanos (brazos) las distancias siguientes
L1 = 118.4mm
L2 = 113.6mm
L3 = 53.6mm
In [19]: # DATOS
RH1 = 1461 # N
RH2 = 3097 # N
RV1 = 1539 #N
RV2 = 160 # N
L1 = 118.4 # mm
L2 = 113.6 # mm
L3 = 53.6 # mm
9.Determinación de los diagramas de torque y momento flector
Los momentos flectores se van a dar con respecto a los planos: horizontal y vertical
M V 1 = RV 1. L2
M V 2 = RV 2. L3
M H = RH 1. L2
In [20]: MV1 = RV1 * L2
MV1 = round (MV1,0)
print("MV1= ","{:.0f}".format(MV1),'N.mm')
MV1= 174830 N.mm
In [21]: MV2 = RV2 * L3
MV2 = round (MV2,0)
print("MV2= ","{:.0f}".format(MV2),'N.mm')
MV2= 8576 N.mm
In [22]: MH = RH1 * L2
MH = round (MH,2)
print("MH= ","{:.1f}".format(MH),'N.mm')
MH= 165969.6 N.mm
Momentos totales:
2 2 1/2
M 1 = (M H + MV 1 )
2 2 1/2
M 2 = (M H + MV 2 )
In [23]: M1 = (MH**2 + MV1**2)**(1/2)
M1 = round (M1,0)
print("M1= ","{:.0f}".format(M1),'N.mm')
M1= 241063 N.mm
In [24]: M2 = (MH**2 + MV2**2)**(1/2)
M2 = round (M2,0)
print("M2= ","{:.0f}".format(M2),'N.mm')
M2= 166191 N.mm
10. Análisis de resistencia:
Se deberia verificar la resistencia del eje a combinacion (Torque y Flector)
Estática
2 2 1/2
(M 1 + (α ∗ T ) )
σca =
W
Torque:
T = 756000N . mm
Dinámica
El momento flector en el zona: a diametro mayor; concentracion de M, y accion T:
para este caso transision entre el escalon V y VI
L2 − 25
M = M1 ∗
L2
M
σf =
W
In [25]: M = (M1 * ((L2-25)/(L2)))
M = round (M,0)
print("M= ","{:.0f}".format(M),'Mpa')
M= 188012 Mpa
In [26]: W = 34300
σf = M/W
σf = round (σf,2)
print("σf= ","{:.2f}".format(σf),'Mpa')
σf= 5.48 Mpa
Torque:
T = 756000M pa
T
τf = = 11.0M pa
Wt
Fa 1002N
σa = = = 0.26M pa
pi
A 2
∗ 70
4
De tablas:
σb = 600M pa
σ1 = 240M pa
τ1 = 140M pa
Factores de diseño
Se determinan los factores de influencia
Eje de transmision
f iσ = 0.1
f iτ = 0.05
Referencial
Kσ = 2.95
Kτ = 2.33
V-VI
σ−1
Sσ =
Kσ ∗ σa + f iσ ∗ σm
τ1
Sτ =
Kτ ∗ τa + f iτ ∗ τm
Sσ ∗ Sτ
Sca = = 8.68
2 2 1/2
(Sσ + Sτ )
In [27]: # Datos
σb = 600 # Mpa
σ11 = 240 # Mpa
τ1 = 140 # Mpa
fiσ = 0.1
fiτ = 0.05
kσ = 2.95
kτ = 2.33
In [28]: # FACTOR DE CARGA a
a = 0.6
dVI = 70 # mm
In [29]: for i in range(0, len(P_3)):
T.append((9.55*10**6) * (P_3[i]/w_3))
σca.append((((M1**2)+(a*T[i])**2)**(1/2)/ W))
τf.append((16*T[i])/(pi * dVI**3))
Ft.append(((2*T[i])/(dz)))
Fa.append((Ft[i]*tan(math.radians(β))))
σa.append(Fa[i] / ((pi/4)*dVI**2))
Sσ.append(((σ11) / ((kσ*σf) + (fiσ*σa[i]))))
Sτ.append(((τ1) / ((kτ*(τf[i]/2)) + (fiτ*(τf[i]/2)))))
Sca.append(((Sσ[i]*Sτ[i]) / ((Sσ[i]**2 + Sτ[i]**2)**(1/2))))
In [ ]:
In [30]: print(" Tabla de Iteración 2")
table= PrettyTable([RED+"P_3[KW]",BLACK+"T[N.mm]","σca[MPa]","τf[MPa]","σa[MPa]","Sσ
for i in range(0,len(P_3)):
table.add_row(([RED+"{:.1f}".format(P_3[i]),BLACK+"{:.1f}".format(T[i]/1000),
round(Sσ[i],3),round(Sτ[i],3),round(Sca[i],3) ]))
print(table)
Tabla de Iteración 2
+---------+---------+----------+---------+---------+--------+-------+-------+
| P_3[KW] | T[N.mm] | σca[MPa] | τf[MPa] | σa[MPa] | Sσ | Sτ | Sca |
+---------+---------+----------+---------+---------+--------+-------+-------+
| 9.5 | 756.0 | 14.98 | 11.23 | 0.26 | 14.822 | 10.48 | 8.557 |
| 10.0 | 795.8 | 15.59 | 11.82 | 0.27 | 14.821 | 9.956 | 8.264 |
| 10.5 | 835.6 | 16.22 | 12.41 | 0.29 | 14.82 | 9.482 | 7.987 |
| 11.0 | 875.4 | 16.85 | 13.0 | 0.3 | 14.818 | 9.051 | 7.724 |
| 11.5 | 915.2 | 17.48 | 13.59 | 0.32 | 14.817 | 8.657 | 7.475 |
| 12.0 | 955.0 | 18.12 | 14.18 | 0.33 | 14.816 | 8.297 | 7.239 |
| 12.5 | 994.8 | 18.77 | 14.77 | 0.34 | 14.815 | 7.965 | 7.015 |
| 13.0 | 1034.6 | 19.41 | 15.36 | 0.36 | 14.813 | 7.658 | 6.803 |
| 13.5 | 1074.4 | 20.06 | 15.95 | 0.37 | 14.812 | 7.375 | 6.602 |
| 14.0 | 1114.2 | 20.72 | 16.54 | 0.38 | 14.811 | 7.111 | 6.411 |
| 14.5 | 1154.0 | 21.37 | 17.13 | 0.4 | 14.81 | 6.866 | 6.229 |
| 15.0 | 1193.8 | 22.03 | 17.73 | 0.41 | 14.808 | 6.637 | 6.057 |
+---------+---------+----------+---------+---------+--------+-------+-------+
In [58]: Fs = np.arange(1.5,3.8,0.2) # Factor de seguridad
In [59]: for i in range(0, len(Fs)):
Fs1.append((Fs[i]))
In [64]: print(" Tabla de Iteración 3")
table= PrettyTable([RED+"Fs",BLACK+"Sca"])
for i in range(0,len(Fs)):
table.add_row(([RED+"{:.1f}".format(Fs1[i]),BLACK+"{:.3f}".format(Sca[i]) ]))
print(table)
Tabla de Iteración 3
+-----+-------+
| Fs | Sca |
+-----+-------+
| 1.5 | 8.557 |
| 1.7 | 8.264 |
| 1.9 | 7.987 |
| 2.1 | 7.724 |
| 2.3 | 7.475 |
| 2.5 | 7.239 |
| 2.7 | 7.015 |
| 2.9 | 6.803 |
| 3.1 | 6.602 |
| 3.3 | 6.411 |
| 3.5 | 6.229 |
| 3.7 | 6.057 |
+-----+-------+
Gráficos iterativos
In [34]: fig, ax1 = plt.subplots(figsize=(8, 4))
ax1.set_title("Potencia transmitida vs Diámetro mínimo del eje", fontsize=16, color=
ax1.set_xlabel("Potencia transmitida (KW)", fontsize=16, color="blue")
ax1.plot(P_3,d, lw=2, color="green")
ax1.grid()
ax1.set_ylabel(r"Diámetro mínimo del eje (mm)", fontsize=16, color="blue")
for label in ax1.get_yticklabels():
label.set_color("black")
In [35]: fig, ax1 = plt.subplots(figsize=(8, 4))
ax1.set_title("Potencia transmitida vs Fuerza tangencial", fontsize=16, color="orang
ax1.set_xlabel("Potencia transmitida (KW)", fontsize=16, color="red")
ax1.plot(P_3,d, lw=2, color="purple")
ax1.grid()
ax1.set_ylabel(r"Fuerza tangencial (N)", fontsize=16, color="red")
for label in ax1.get_yticklabels():
label.set_color("black")
In [36]: fig, ax1 = plt.subplots(figsize=(8, 4))
ax1.set_title("Potencia transmitida vs Fuerza radial", fontsize=16, color="orange")
ax1.set_xlabel("Potencia transmitida (KW)", fontsize=16, color="blue")
ax1.plot(P_3,d, lw=2, color="red")
ax1.grid()
ax1.set_ylabel(r"Fuerza radial (N)", fontsize=16, color="blue")
for label in ax1.get_yticklabels():
label.set_color("black")
In [37]: fig, ax1 = plt.subplots(figsize=(8, 4))
ax1.set_title("Potencia transmitida vs Fuerza axial", fontsize=16, color="orange")
ax1.set_xlabel("Potencia transmitida (KW)", fontsize=16, color="purple")
ax1.plot(P_3,d, lw=2, color="blue")
ax1.grid()
ax1.set_ylabel(r"Fuerza axial (N)", fontsize=16, color="purple")
for label in ax1.get_yticklabels():
label.set_color("black")
In [38]: fig, ax1 = plt.subplots(figsize=(8, 4))
ax1.set_title("Esfuerzo axial vs Sca", fontsize=16, color="orange")
ax1.set_xlabel("Esfuerzo axial (MPa)", fontsize=16, color="purple")
ax1.plot(σa,Sca, lw=2, color="yellow")
ax1.grid()
ax1.set_ylabel(r"Sca", fontsize=16, color="purple")
for label in ax1.get_yticklabels():
label.set_color("black")
In [65]: fig, ax1 = plt.subplots(figsize=(8, 4))
ax1.set_title("Factor de seguridad vs Factor de seguridad calculado", fontsize=16, c
ax1.set_xlabel("Factor de seguridad", fontsize=16, color="blue")
ax1.plot(Fs,Sca, lw=2, color="red")
ax1.grid()
ax1.set_ylabel(r"Factor de seguridad calculado)", fontsize=16, color="blue")
for label in ax1.get_yticklabels():
label.set_color("black")
EJE EN ANSYS
Conclusiones
Se puede concluir que a medida que se aumente la potencia transmistida en el sistema
teniendo el mismo material del eje en todos los casos, este tiende a reducirse
considerablemente su factor de seguridad.
Se puede concluir que si se tiene una mayor potencia transmitida las fuerzas radiales y
axiales tienden a crecer de igual manera, es decir son directamente proporcianales entre
ellas.
Se puede concluir que si se deseara continuar con un factor de seguridad alto, se tendria
que ir cambiando el material que sea el más adeacudo para dicha aplicación.
Se puede concluir que si se deseara continuar con un factor de seguridad alto, se tendria
que ir cambiando el material que sea el más adeacudo para dicha aplicación.
Se concluye que según este calculo se pudo comparar si los valores de seguridad
calculados para el eje del ejercicio, cuando podía ser considerado como seguro o no.
Se logro variar los valores del factor de seguridad para ser comparados con el valor de
seguridad calculado , mediantes las respectivas ecuaciones.
