INF Prac PDF
INF Prac PDF
INF Prac PDF
Siendo el día 10 de enero del año 2019, en las instalaciones de Makyl Engineering consulting
services E.I.R.L, a cargo del Gerente Ingeniero Lenin Chávez Callo y en coordinación con la
Universidad Nacional San Antonio Abad de Cusco, es que pude realizar mis practicas pre
profesionales 1 referidos a Ingeniería Mecánica, dichas prácticas tuvieron un periodo desde la
fecha de inicio indicada hasta el 8 de mayo en la ciudad de Lima y posteriormente finalizando
las practicas el 1 de junio del año en curso en la ciudad de Cusco, se realizó dicha práctica de
la siguiente manera:
1
INGENIERÍA MECÁNICA UNSAAC
ÍNDICE
2
INGENIERÍA MECÁNICA UNSAAC
3
INGENIERÍA MECÁNICA UNSAAC
CONCLUSIONES................................................................................................................................... 110
BIBLIOGRAFIA..................................................................................................................................... 112
LISTA DE FIGURAS
Figura 8 Viga en voladizo con una fuerza aplicada en el extremo. (Singiresu, 2010) ..................... 18
Figura 15. Diagrama de momento flector correspondiente a una fuerza unitaria en m1. (Elaboración
propia). ................................................................................................................................ 27
Figura 16. Diagrama de momento flector correspondiente a una fuerza unitaria en m2. (Elaboración
propia). ................................................................................................................................ 28
4
INGENIERÍA MECÁNICA UNSAAC
Figura 17. Diagrama de momento flector correspondiente a una fuerza unitaria en m3. (Elaboración
propia). ................................................................................................................................ 28
Figura 18. Diagrama de momento flector correspondiente a una fuerza unitaria en m4. (Elaboración
propia). ................................................................................................................................ 28
Figura 19. Tablas de integración para hallar la matriz de flexibilidad. (Alencastre Miranda, 2018) 30
Figura 26. Estudio en elementos finitos previo. (Elaboración Propia - SolidWorks) ..................... 34
Figura 27. DAC con los cables posicionados en cada canal. (Elaboración Propia) ...................... 35
Figura 28. Acelerómetro con su respectivo cable codificado. (Elaboración Propia) ..................... 35
Figura 29. Posicionamiento de los sensores en la estructura roja. (Elaboración Propia) .............. 35
Figura 30. Experimento del péndulo realizado con fines didácticos. (Elaboración Propia)............. 36
Figura 31. Entorno de programación gráfica de LabView para el experimento del péndulo.
(Elaboración Propia - LabView) ............................................................................................. 36
Figura 32. Resultados Medidos en el experimento del péndulo. (Elaboración Propia - Labview) ... 37
Figura 33. Programación grafica del programa de medición. (Elaboración Propia - Labview) ...... 37
Figura 34. Interfaz del Programa de Medición. (Elaboración Propia - Labview) .......................... 38
Figura 36. FFT de la Fuerza excitadora deseada. (Elaboración Propia - Labview) ....................... 38
Figura 38. Programación grafica del programa de procesamiento de datos. (Elaboración Propia -
Labview) .............................................................................................................................. 39
Figura 39. Interfaz del programa de procesamiento de datos. (Elaboración Propia - Labview)...... 40
5
INGENIERÍA MECÁNICA UNSAAC
Figura 42. Diagrama en el dominio de la frecuencia FFT – Módulo Amarillo. (Elaboración Propia)
........................................................................................................................................... 41
Figura 43. Densidad Espectral de Potencia – Módulo amarillo. (Elaboración Propia) .................. 42
Figura 44. Diagrama en el dominio de la frecuencia FFT– Módulo rojo. (Elaboración propia) ..... 43
Figura 49. Para facilitar este proceso los cables están etiquetados. (Elaboración Propia) ............. 45
Figura 52. Modo de vibración a 5.645 Hz. (Elaboración Propia-ARTeMIS Modal) ....................... 47
Figura 60. Evaluación Técnica Económica - Concepto Solución . (Elaboración Propia) ................ 59
6
INGENIERÍA MECÁNICA UNSAAC
Figura 72. Vista de ensamble del proyecto óptimo. (Elaboración Propia) .................................... 71
Figura 73. Faja Radial Stacker 230-CVB-002 – Antamina (Fuente: Informe interno Makyl) .......... 75
Figura 75. Ubicación geográfica de la mina Antamina y una vista del asentamiento minero. (Fuente:
Internet) ............................................................................................................................... 76
Figura 76. Modelo conceptual general de la Faja Transportadora Stacker Radial 230- CVB 002.
(Fuente: Informe interno Makyl) ............................................................................................. 77
Figura 77. Modelo conceptual general simplificado de la faja 230 – CVB 002. (Fuente: Informe
interno Makyl) ...................................................................................................................... 78
Figura 78. Modelo general simplificado de dos grados de libertad. (Fuente: Informe interno Makyl)
........................................................................................................................................... 79
Figura 79. Principio básico para hallar la rigidez idealizándolo como elemento viga. (Fuente:
Informe interno Makyl) .......................................................................................................... 81
Figura 80. Modelo Numérico creado para el cálculo de las rigideces. (Fuente: Informe interno
Makyl) ................................................................................................................................. 82
Figura 81. Vista en el plano XZ del Stacker Radial. (Fuente: Informe interno Makyl) ................... 84
Figura 82. Modelo conceptual discreto. (Fuente: Informe interno Makyl) .................................... 84
Figura 83. Momento de área equivalente I1. (Fuente: Informe interno Makyl) .............................. 86
Figura 84. Momento de área equivalente I2. (Fuente: Informe interno Makyl) .............................. 86
Figura 85. Momento de área equivalente I3. (Fuente: Informe interno Makyl) .............................. 86
Figura 86. Momento de área equivalente I4. (Fuente: Informe interno Makyl) .............................. 86
Figura 87. Cálculo de la constante de rigidez de un resorte helicoidal. (Meruane, 2012) ............... 88
Figura 91. Variación de la Frecuencia Natural en función del incremento de K1. (Fuente: Informe
interno Makyl) ...................................................................................................................... 94
7
INGENIERÍA MECÁNICA UNSAAC
Figura 92. Variación de la Frecuencia Natural en función del incremento de K2. (Fuente: Informe
interno Makyl) ...................................................................................................................... 94
Figura 93. Variación de la Frecuencia Natural en función del incremento de K3. (Fuente: Informe
interno Makyl) ...................................................................................................................... 95
Figura 94. Vista en el plano XZ del Stacker Radial. (Fuente: Informe interno Makyl) .................... 95
Figura 95. Modelo conceptual para el análisis. (Fuente: Informe interno Makyl) .......................... 96
Figura 97. Función de amplificación de M2. (Fuente: Informe interno Makyl) ........................... 100
Figura 98. Función de amplificación de M3. (Fuente: Informe interno Makyl)............................ 100
Figura 99. Función de amplificación de M4. (Fuente: Informe interno Makyl) ........................... 101
Figura 100. Idea conceptual de una posible solución mediante absorsor. (Fuente: Informe interno
Makyl) ............................................................................................................................... 101
Figura 101. Análisis modal del modelo discreto. (Fuente: Informe interno Makyl) ...................... 102
Figura 102. Frecuencia de Resonancia a 0.53Hz de la Estructura en el plano XZ. (Fuente: Informe
interno Makyl) .................................................................................................................... 103
Figura 103. Frecuencia de Resonancia a 1.68Hz de la Estructura en el plano XZ. (Fuente: Informe
interno Makyl) .................................................................................................................... 104
Figura 104. Frecuencia de Resonancia a 3.96Hz de la Estructura en el plano XZ. (Fuente: Informe
interno Makyl) .................................................................................................................... 104
Figura 105. Frecuencia de Resonancia a 3.96Hz de la Estructura en el plano XZ. (Fuente: Informe
interno Makyl) .................................................................................................................... 104
Figura 106. Vista del Stacker Radial con soluciones implementadas. (Fuente: Informe interno Makyl)
......................................................................................................................................... 105
Figura 107. Sistema de sujeción. (Fuente: Informe interno Makyl) ............................................ 106
Figura 108. Pórticos de anclaje superior e inferior. (Fuente: Informe interno Makyl)................. 106
Figura 109. Reforzamiento de los polines inferiores y superiores. (Fuente: Informe interno Makyl)
......................................................................................................................................... 107
Figura 110. Masas incorporadas entre el Ball Ring y el Travel Gear. (Fuente: Informe interno Makyl)
......................................................................................................................................... 107
Figura 111. Masas adicionales agregadas a la estructura. (Fuente: Informe interno Makyl) ........ 108
Figura 112. Reforzamiento del sistema motriz. (Fuente: Informe interno Makyl) ......................... 108
8
INGENIERÍA MECÁNICA UNSAAC
9
INGENIERÍA MECÁNICA UNSAAC
Capítulo 1
Introducción y generalidades
10
INGENIERÍA MECÁNICA UNSAAC
Periodo de prácticas
Las prácticas fueron de 90 días con 9 horas diarias de lunes a viernes y 5 horas los
días sábados.
11
INGENIERÍA MECÁNICA UNSAAC
Capítulo 2
Información general de la empresa
Empresa
Makyl Engineering Consulting & Services, se dedica principalmente a los rubros de:
Diseño y ejecución de proyectos en actividades mineras e industriales.
Objetivo
El objetivo principal de esta empresa es ser líder brindando los mejores servicios en
los rubros antes mencionado y capacitando constantemente a su personal que cuenta para
estos servicios.
Ubicación
Av. Canadá N° 3644 - Piso 3 – Lima, Perú
Infraestructura tecnológica
La Empresa Makyl Engineering Consulting & Services para los distintos servicios en
el rubro de la minería:
Para la simulación en elementos finitos, bien para el cálculo estructural y para un
análisis modal: SolidWorks, ANSYS y SAP 2000.
12
INGENIERÍA MECÁNICA UNSAAC
Organización de la empresa
13
INGENIERÍA MECÁNICA UNSAAC
Capítulo 3
Descripción de las principales actividades realizadas en la empresa
14
INGENIERÍA MECÁNICA UNSAAC
Capítulo 4
Marco teórico
Definición de vibración
Cualquier movimiento que se repite después de un intervalo de tiempo se llama
vibración u oscilación. El vaivén de un péndulo y el movimiento de una cuerda pulsada son
ejemplos comunes de vibración. La teoría de la vibración tiene que ver con el estudio de los
movimientos oscilatorios de los cuerpos y las fuerzas asociadas con ellos. (Meruane, 2012)
Conceptos básicos
Ciclo. Al movimiento de un cuerpo vibratorio desde su posición no perturbada o de equilibrio
hasta su posición en una dirección, y luego de vuelta a la posición de equilibrio, y luego a su
posición extrema en la otra dirección, y de vuelta a la posición de equilibrio se le llama ciclo
de vibración. (Meruane, 2012)
Frecuencia natural. Si se deja que un sistema vibre por sí mismo después de una
perturbación inicial, la frecuencia con la cual oscila sin la acción de fuerzas externas se
conoce como frecuencia natural. (Singiresu, 2010)
15
INGENIERÍA MECÁNICA UNSAAC
Sistemas de un grado de libertad. Se llama sistema de un solo grado de libertad, ya que una
coordenada (x) es suficiente para especificar la posición de la masa en cualquier momento.
No existe ninguna fuerza externa aplicada a la masa, de ahí que el movimiento resultante de
una perturbación inicial será una vibración libre.
16
INGENIERÍA MECÁNICA UNSAAC
17
INGENIERÍA MECÁNICA UNSAAC
Figura 8 Viga en voladizo con una fuerza aplicada en el extremo. (Singiresu, 2010)
18
INGENIERÍA MECÁNICA UNSAAC
19
INGENIERÍA MECÁNICA UNSAAC
20
INGENIERÍA MECÁNICA UNSAAC
(12)
21
INGENIERÍA MECÁNICA UNSAAC
Las ecuaciones (10), (11) y (13) son funciones de tiempo armónicas. El movimiento
es simétrico con respecto a la posición de equilibrio de la masa “m”. La velocidad es un
máximo y la aceleración es cero cada vez que la masa pasa por esta posición. En los
desplazamientos extremos, la velocidad es cero y la aceleración es un máximo. El sistema de
resorte-masa se conoce como oscilador armónico. (Singiresu, 2010)
La ecuación (11) se puede expresar en una forma diferente si introducimos la notación
Donde A y son las constantes nuevas, las cuales se pueden expresar en función de
y como:
22
INGENIERÍA MECÁNICA UNSAAC
23
INGENIERÍA MECÁNICA UNSAAC
ni a .(Meruane, 2012)
Pre-multiplicando por
24
INGENIERÍA MECÁNICA UNSAAC
Reordenando y definiendo :
escalar que satisface la anterior ecuación para un vector se denomina valor propio y
La primera ecuación nos dice que dos modos distintos son ortogonales con la matriz
de masa (también con la matriz de rigidez) y cuando se calcula el producto del mismo modo
con la matriz de masa, como en la segunda ecuación, el producto no es cero (Meruane, 2012).
25
INGENIERÍA MECÁNICA UNSAAC
y son matrices diagonales denominadas matriz de masa modal y matriz de rigidez modal
respectivamente. (Meruane, 2012)
Con el objetivo de desacoplar las ecuaciones de movimiento asumimos una solución,
26
INGENIERÍA MECÁNICA UNSAAC
Figura 15. Diagrama de momento flector correspondiente a una fuerza unitaria en m1.
(Elaboración propia).
27
INGENIERÍA MECÁNICA UNSAAC
Figura 16. Diagrama de momento flector correspondiente a una fuerza unitaria en m2.
(Elaboración propia).
Figura 18. Diagrama de momento flector correspondiente a una fuerza unitaria en m4.
(Elaboración propia).
28
INGENIERÍA MECÁNICA UNSAAC
En donde:
Corresponde a la deflexión
La energía elástica de deformación elástica varía en cada caso que se nos presente,
pero en el caso de flexión de vigas se puede demostrar que: (Alencastre Miranda, 2018)
29
INGENIERÍA MECÁNICA UNSAAC
Figura 19. Tablas de integración para hallar la matriz de flexibilidad. (Alencastre Miranda,
2018)
Por lo cual aplicamos el álgebra matricial para obtener una matriz de coeficientes
A través de la ley de Hooke, se obtiene que la matriz de rigidez, que viene a ser la
inversa de la matriz de rigidez. (Alencastre Miranda, 2018)
30
INGENIERÍA MECÁNICA UNSAAC
Capítulo 5
Pruebas de laboratorio
31
INGENIERÍA MECÁNICA UNSAAC
32
INGENIERÍA MECÁNICA UNSAAC
Cables conectores
En el laboratorio se disponen de cables desde 10 metros hasta 30 metros. En
nuestro caso solo será necesario usar los de 10 metros a menos.
Experimentación
33
INGENIERÍA MECÁNICA UNSAAC
34
INGENIERÍA MECÁNICA UNSAAC
Figura 27. DAC con los cables posicionados en cada canal. (Elaboración Propia)
Luego de haber realizado los pasos previos se procedió a programar en LabView para
el procesamiento y manejo de señales y así poder obtener nuestros resultados de medición.
35
INGENIERÍA MECÁNICA UNSAAC
Figura 30. Experimento del péndulo realizado con fines didácticos. (Elaboración Propia)
Figura 31. Entorno de programación gráfica de LabView para el experimento del péndulo.
(Elaboración Propia - LabView)
36
INGENIERÍA MECÁNICA UNSAAC
Figura 32. Resultados Medidos en el experimento del péndulo. (Elaboración Propia - Labview)
Figura 33. Programación grafica del programa de medición. (Elaboración Propia - Labview)
37
INGENIERÍA MECÁNICA UNSAAC
Una vez teniendo nuestro programa de mediciones listo para operar procedemos a
realizar impactos en la estructura con el martillo modal, procurando que no ocurran rebotes
en el golpe y utilizando la punta adecuada para la estructura o elemento a medir. A
continuación, se puede ver la entrada deseada y su respectiva FFT.
38
INGENIERÍA MECÁNICA UNSAAC
Los impactos deben ser en los puntos que previamente se establecieron, procurando
que sean en distintas direcciones de referencia en X, Y, Z. En nuestro caso se necesitaron 6
golpes, dos en cada eje ya que se tenía dos acelerómetros triaxiales y cuatro uniaxiales. Esto
en total suma 10 canales y más el martillo modal serian 11 canales.
Figura 38. Programación grafica del programa de procesamiento de datos. (Elaboración Propia -
Labview)
39
INGENIERÍA MECÁNICA UNSAAC
Figura 39. Interfaz del programa de procesamiento de datos. (Elaboración Propia - Labview)
40
INGENIERÍA MECÁNICA UNSAAC
Resultados
Figura 42. Diagrama en el dominio de la frecuencia FFT – Módulo Amarillo. (Elaboración Propia)
41
INGENIERÍA MECÁNICA UNSAAC
42
INGENIERÍA MECÁNICA UNSAAC
43
INGENIERÍA MECÁNICA UNSAAC
- 1 - 25.75
- 2 - 36
44
INGENIERÍA MECÁNICA UNSAAC
Figura 49. Para facilitar este proceso los cables están etiquetados. (Elaboración Propia)
Una vez grabada la data de aceleración de los canales disponibles se pasa a preparar la
data para el análisis.
45
INGENIERÍA MECÁNICA UNSAAC
46
INGENIERÍA MECÁNICA UNSAAC
47
INGENIERÍA MECÁNICA UNSAAC
48
INGENIERÍA MECÁNICA UNSAAC
Capítulo 6
Memoria descriptiva
Requerimientos generales
El proyecto tiene los siguientes requerimientos:
Diseñar un modelo zaranda a escala 1:15.
Realizar planos de ensamble y despice.
Construir el modelo zaranda.
Medir las vibraciones del modelo zaranda.
Ser transportable.
Limitaciones
Para la medición de vibraciones se utilizarán los acelerómetros proporcionados por la
empresa, debido a su elevado costo que imposibilita adquirir otro tipo de acelerómetros.
El presente capitulo solamente contemplará la fase de diseño del prototipo de
laboratorio Zaranda y no la fase de pruebas de vibraciones, debido a que como practicante
participe en la fase de diseño y luego me asignaron otra tarea que será el tema del capítulo 7
del presente informe.
Lista de exigencias
En la lista de exigencias se detalla las condiciones que debe cumplir el diseño para
satisfacción del cliente.
49
INGENIERÍA MECÁNICA UNSAAC
LISTA DE EXIGENCIAS
Categoría D/E Descripción
El Modelo a escala de la zaranda vibratoria tiene por
Función
E objetivo medir los modos de vibración a través del
Principal
análisis modal.
Energía E Baterías y/o DC a través de reguladores de voltaje.
Las dimensiones del modelo son de una relación de
Geometría E 1:15 y debe ser lo más aproximado posible al modelo
real.
El modelo contará con 6 grados de libertad:
Cinemática E Giro en el eje x, y y z
Adelante/atrás en el eje x,y y z
El equipo debe ser capaz de vibrar a distintos modos
Dinámica E
de vibración.
El diseño del sistema debe permitir que se pueda
Operatividad E
modificar la fuerza de excitación.
Los materiales deben encontrarse en el mercado
Fabricación E
nacional.
El modelo debe ser desarmable para facilitar el
Montaje E
traslado a largas distancias.
Material E Los materiales serán reusados.
Señal de entrada: Energía eléctrica del motor
Señales E
Señal de salida: Señal de control de respuesta.
Los costos de fabricación deberán ser lo mínimo
Costo E
posible.
50
INGENIERÍA MECÁNICA UNSAAC
Concepción de la solución
Estructura de funciones
Con las funciones obtenidas anteriormente se representa de manera esquemática la
estructura de funciones del Modelo Zaranda.
Concepto de solución
De acuerdo a la estructura de funciones definida se representa la matriz morfológica.
Matriz morfológica
En la matriz morfológica se coloca los portadores de función posibles.
51
INGENIERÍA MECÁNICA UNSAAC
MATRIZ MORFOLOGICA
Alternativa
FUNCIONES Alternativa 1 Alternativa 2 Alternativa 3 Alternativa 4
5
Motor de
Motor
inducción
DC
TRANSFORMACION
Solenoide
Masa
excéntrica
EXCITACIÓN
Resorte de Muelle
Resorte torsión
VIBRACION
52
INGENIERÍA MECÁNICA UNSAAC
Conceptos de soluciones
Se esquematiza el Concepto de Solución 1
53
INGENIERÍA MECÁNICA UNSAAC
54
INGENIERÍA MECÁNICA UNSAAC
55
INGENIERÍA MECÁNICA UNSAAC
56
INGENIERÍA MECÁNICA UNSAAC
Impor Soluciones
No. Criterios Técnicos tancia
(i) S1 S2 S3 S4
1 Función 12% 3 3 3 3
4 Cinemática 10% 3 2 2 3
6 Fabricación 10% 3 3 2 1
7 Complejidad 10% 3 3 3 1
8 Facilidad de manejo 8% 3 3 2 1
9 Transportabilidad 8% 4 4 3 1
10 Manipulación 8% 3 3 3 1
57
INGENIERÍA MECÁNICA UNSAAC
Import Soluciones
No. Criterios Económicos
ancia (i) S1 S2 S3 S4
58
INGENIERÍA MECÁNICA UNSAAC
DIAGRAMA DE EVALUACIÓN
100%
90%
S1
80%
S2
70%
Criterios Económicos
60%
S3
50%
S4
40%
30%
20%
10%
0%
0% 20% 40% 60% 80% 100%
Criterios TéCNICOS
59
INGENIERÍA MECÁNICA UNSAAC
60
INGENIERÍA MECÁNICA UNSAAC
viga, los soportes del puente y de los distintos perfiles que se utilizaran para crear la base y la
estructura de la criba vibratoria.
El perforado de los distintos agujeros necesario para unir distintas piezas se realizará
mediante el taladro de banco. El proceso de soldadura será por arco eléctrico con electrodos
revestidos (SMAW). Se realizará el esmerilado mediante la amoladora con disco de desbaste
para recortar los filos cortantes y dar mejor acabado superficial.
Se realizará la construcción del soporte de concreto en el cual están incluidos los
procesos de enrejado encofrado y vaciado. Se moldeará plástico a través de una matriz para
las piezas solicitadas como uniones o acoples.
61
INGENIERÍA MECÁNICA UNSAAC
62
INGENIERÍA MECÁNICA UNSAAC
63
INGENIERÍA MECÁNICA UNSAAC
64
INGENIERÍA MECÁNICA UNSAAC
65
INGENIERÍA MECÁNICA UNSAAC
66
INGENIERÍA MECÁNICA UNSAAC
Proyecto Preliminar 1:
- El acople del eje de la viga puente con el eje del motor es mediante una junta
mecánica de dos u.
Proyecto Preliminar 2:
- El acople del eje de la viga puente con el eje del motor es mediante
engranajes de reductor de velocidad.
Proyecto Preliminar 3:
- El acople del eje de la viga puente con el eje del motor es mediante bridas de
acople.
67
INGENIERÍA MECÁNICA UNSAAC
68
INGENIERÍA MECÁNICA UNSAAC
69
INGENIERÍA MECÁNICA UNSAAC
70
INGENIERÍA MECÁNICA UNSAAC
71
INGENIERÍA MECÁNICA UNSAAC
armónica. El concepto si satisface dicha función, pero aun las dimensiones de los resortes no
han sido establecidas.
La función de medición está claramente definida gracias a que el laboratorio cuenta
con acelerómetros imantados para que la maqueta pueda cumplir su misión principal la cual
es poder medir los distintos modos de vibración en una maquina como esta en escala real.
72
INGENIERÍA MECÁNICA UNSAAC
73
INGENIERÍA MECÁNICA UNSAAC
Elaboración del procedimiento para realizar un análisis dinámico a una faja tipo
Stacker radial
El presente documento se realizó por iniciativa de la empresa de obtener un
procedimiento a seguir para el análisis dinámico estructural de una faja apiladora radial o en
ingles un Radial Stacker, ya que este tipo de máquinas es usado intensamente en la industria
minera. Se recurrirá a un caso de referencia resuelto en la empresa anteriormente el cual fue
la faja Overland número 230-CVB-002 de la compañía Minera Antamina S.A. la cual
implementó el Debottlenecking Project 130K que tenía por objetivo incrementar la capacidad
de procesamiento y transporte de mineral a 180 KTPD, para lo cual se efectuó el up grade de
los sistemas motrices de la faja Overland 220-CVB-001 y 230-CVB-002. La producción se
sostuvo entre las 7000 tph y las 8000 tph, lo que deriva en el incremento de carga derivada a
la faja stacker 230-CVB-002 a un rango de 4000 a 6000 tph. Se tomaron medidas correctivas
para poder satisfacer la nueva demanda, pero no se consideró la influencia de estos cambios
en la dinámica estructural del sistema.
74
INGENIERÍA MECÁNICA UNSAAC
Figura 74. Faja Radial Stacker 230-CVB-002 – Antamina (Fuente: Informe interno Makyl)
Faja transportadora stacker radial. Según la clasificación por funcionalidad, dentro de las
fajas apiladoras, se encuentran los apiladores radiales o los Radial Stacker por su nombre en
inglés. Son capaces de girar sobre un apoyo con un gran rodamiento o Ball Ring que funciona
como pivote, permitiéndoles apilar distintas rumas de material gracias a su Angulo de giro
que logra a través de un Travel Gear.
Ubicación del caso del caso de estudio. Antamina está ubicada en Ubicada en el
departamento de Ancash a 270 Km. al noreste de Lima, Perú. Antamina se ha constituido en
uno de los mayores productores nacionales de concentrados de cobre y zinc del país.
75
INGENIERÍA MECÁNICA UNSAAC
Desarrollo del análisis del problema. Luego de realizar la visita y levantamiento de campo
el primer paso es configurar un modelo conceptual que nos ayude a representar el sistema y
replicar las condiciones de trabajo en la cual ocurren los problemas de vibración.
Para ello se recurre a crear un modelo conceptual general que servirá de punto de
partida para realizar distintos análisis. El proceso de obtener el modelo analítico requiere
tener conocimientos básicos sobre la teoría de vibraciones, así como también de estática y
resistencia de materiales. El resultado dependerá de la habilidad y practica del diseñador,
pero sobretodo paciencia y buen humor.
76
INGENIERÍA MECÁNICA UNSAAC
Figura 77. Modelo conceptual general de la Faja Transportadora Stacker Radial 230- CVB 002.
(Fuente: Informe interno Makyl)
77
INGENIERÍA MECÁNICA UNSAAC
Figura 78. Modelo conceptual general simplificado de la faja 230 – CVB 002. (Fuente: Informe
interno Makyl)
Enfoque del análisis. El enfoque de análisis que demos al problema dependerá de donde está
ocurriendo el problema, o mejor dicho en que plano está ocurriendo el problema. El primer
enfoque que se le dio al problema fue en el plano XY al cual el modelo se redujo a
simplemente un problema de 2 grados de libertad. Como se verá más adelante, este modelo
inicial pudo replicar ciertas frecuencias que afectaban a la estructura, pero para poder
entender el problema de manera más general y poder construir un modelo más fiable, se hizo
un modelo en el plano XZ en el cual se discretizó las masas y se hallaron momentos de
inercia equivalentes. Se mostrará cada coso por separado para poder observar el tipo de
análisis y cuáles son los criterios requeridos para cada uno.
78
INGENIERÍA MECÁNICA UNSAAC
Figura 79. Modelo general simplificado de dos grados de libertad. (Fuente: Informe interno
Makyl)
79
INGENIERÍA MECÁNICA UNSAAC
De aquí, podemos observar que se han hallado dos frecuencias naturales, lo cual es
obvio ya que estamos hablando de un sistema de dos grados de libertad. Con este resultado
nos aproximamos al caso real y podríamos plantear una solución intuitiva la cual se
desarrollaría a profundidad, pero la pregunta ahora es ¿Cómo hallo las masas y rigideces
equivalentes en el sistema?
Para cumplir nuestros objetivos tenemos que recurrir a las recomendaciones para
obtener las magnitudes vibratorias del sistema. En el caso de las masas equivalentes fue algo
más sencillo ya que se recurrió a las especificaciones técnicas ofrecidas por el fabricante de la
faja transportadora en los planos y a los pesos de cada componente que fue otorgado por el
personal técnico de la faja.
80
INGENIERÍA MECÁNICA UNSAAC
Tabla 10. Masa M1: Sistema Motriz. (Fuente: Informe interno Makyl)
Tabla 11. Masa M2: Estructura y polines. (Fuente: Informe interno Makyl)
Figura 80. Principio básico para hallar la rigidez idealizándolo como elemento viga. (Fuente: Informe interno
Makyl)
81
INGENIERÍA MECÁNICA UNSAAC
Figura 81. Modelo Numérico creado para el cálculo de las rigideces. (Fuente: Informe interno Makyl)
Tabla 12. Rigidez K1eq obtenida para configuración actual o antes de la solución. (Fuente: Informe interno
Makyl)
82
INGENIERÍA MECÁNICA UNSAAC
Tabla 13. Rigidez K2eq obtenida para configuración actual o antes de la solución. (Fuente: Informe interno
Makyl)
Tabla 14. Rigidez K3eq obtenida para configuración actual o antes de la solución. (Fuente: Informe interno
Makyl)
83
INGENIERÍA MECÁNICA UNSAAC
Más adelante se verán los resultados obtenidos y se comparara con la solución que se
propuso en ese momento, la cual ya se puede intuir, se basa en variar los parámetros
vibratorios, en este caso, solamente las rigideces.
Análisis en modelo discreto. Como bien se dijo, el modelo anterior de 2 grados de libertad
ayudó a replicar las frecuencias naturales e iterar las rigideces, pero no ayuda a replicar los
modos de vibración en la estructura principal de la faja Stacker Radial. Por lo tanto, en ese
momento se decidió optimizar el modelo y se optó por la discretización del sistema y así
calibrar los resultados a los sucedidos en el fenómeno real. Cabe resaltar que las frecuencias
que podrían inducir a la resonancia pueden ser del motor, engranes del reductor o de la polea
motriz, por lo cual se buscará si existe resonancia a través del siguiente análisis.
m1 m2 m3 m4
x
Apoyo Ball Ring Apoyo Travel Gear
z
L1 L2 L3 L4 L5
Figura 82. Vista en el plano XZ del Stacker Radial. (Fuente: Informe interno Makyl)
m1 m2 m3 m4
L1 L2 L3 L4 L5
84
INGENIERÍA MECÁNICA UNSAAC
estructura que soporta el moto-reductor (se considera la masa de la estructura hasta el ball
ring). La masa M2 representa la masa de una parte de la estructura (desde el ball ring hasta el
punto en que varía la sección) y su respectivo material a transportar.
La masa M3 representa la masa de una parte de la estructura (desde el punto donde la
sección de la estructura empieza a incrementar hasta el travel gear) y su respectivo material a
transportar. La masa M4 representa la masa de una parte de la estructura (desde el travel gear
hasta el final de la estructura) y su respectivo material a transportar.
L1 que es la longitud desde el centro de masa de la masa equivalente M1 hasta el Ball
Ring y tiene un momento de área I1. L2 que es la longitud desde el Ball Ring hasta centro de
masa de la masa equivalente M2, y tiene un momento de área I2. L3 que es la longitud desde
el centro de masa de la masa equivalente M2 hasta centro de masa de la masa equivalente M3
y tiene un momento de área I3. L4 que es la longitud desde el centro de masa de la masa
equivalente M3 hasta el travel gear centro de masa de la masa equivalente M4, y tiene un
momento de área I4. L5 que es la longitud desde el travel gear hasta el centro de masa de la
masa equivalente M5, y tiene un momento de área I5.
A continuación, podremos observar los parámetros del modelo discreto para los dos
casos en el cual este sin carga y con carga, antes de realizar alguna modificación de solución.
Tabla 15. Parámetros del modelo discreto. (Fuente: Informe interno Makyl)
Las siguientes figuras corresponden a las secciones de las cual se tomó el momento de
inercia para el cálculo analítico. Se utilizaron las secciones más críticas de cada sección.
85
INGENIERÍA MECÁNICA UNSAAC
Figura 84. Momento de área equivalente I1. (Fuente: Informe interno Makyl)
Figura 85. Momento de área equivalente I2. (Fuente: Informe interno Makyl)
Figura 86. Momento de área equivalente I3. (Fuente: Informe interno Makyl)
Figura 87. Momento de área equivalente I4. (Fuente: Informe interno Makyl)
86
INGENIERÍA MECÁNICA UNSAAC
0 x 1 K11
M 1 0 0
K12 K13 K14 x1
0 M2 0 0 x 2
K 22 K 23 K 24 x2
0
0 0 M3 0 x3 K 33 K 34 x3
0 0 0 M 4 Sim K 44 x4
x4
De forma simplificada se tiene:
M i x i K ij xi 0
Dónde: es la matriz de masa y es la matriz de rigidez.
87
INGENIERÍA MECÁNICA UNSAAC
Capítulo 7
Memoria de Cálculo
(1)
Datos Valor
G 8E+10 (Pa)
d 0.0018 (m)
R 0.007 (m)
n 8
88
INGENIERÍA MECÁNICA UNSAAC
(2)
(3)
Donde:
X: amplitud de vibración (m)
M: masa de la criba (kg)
m: Masa excéntrica (kg)
R: distancia excéntrica (m)
r: Relación de frecuencias w/wn
Para obtener el valor maximo se utiliza Sen(wt)=+1 o -1; valor máximo cuando wt:
90° o 270°. Sabemos que el motor adquirido y escogido tiene una velocidad nominal de 1200
RPM y se acoplara directamente.
89
INGENIERÍA MECÁNICA UNSAAC
(4)
Una vez verificado esto, remplazamos en la ecuación 3 los valores de M=23.55 kg;
m=0.250 kg y R=0.02 m y obtenemos que la máxima deflexión es aceptable.
Xmáx= 0.007 m
90
INGENIERÍA MECÁNICA UNSAAC
Sabiendo que K1 esté compuesto por 10 resortes y K2 por 14 resortes de 4782 N/m
cada uno, nuestros valores serian:
Tabla 17. Valores para el análisis vibracional a 2GDL (Elaboración Propia)
91
INGENIERÍA MECÁNICA UNSAAC
(9)
(10)
Factor Valor
a 114768
b 8033.76
c 8033.76
d 10604.5632
(11)
Tabla 19. Factores para hallar las frecuencias naturales (Elaboración Propia)
Fac
tor Valor
13.2087179
A 8
-
B 286937.851
115252321
C 0
Frecuencia
natural Valores
92
INGENIERÍA MECÁNICA UNSAAC
Elaboración del procedimiento para realizar un análisis dinámico a una faja tipo
Stacker radial
Análisis para el modelo de dos Grados de Libertad. Luego de haber replicado el problema
a través de nuestro modelo conceptual se podrán dar alternativas de solución las cuales estén
basadas en variar los parámetros de las magnitudes vibratorias, como la rigidez, la masa
incluso el amortiguamiento si es necesario.
En el caso del modelo conceptual de dos grados de libertad, ya tenemos todos
nuestros parámetros listos para iniciar el análisis y proponer alternativas de solución. Para
todo esto debemos auxiliarnos de algún software como Microsoft Excel que nos servirá como
hoja de cálculo e iterar las rigideces ya que las masas, para este caso se consideraran que es
un parámetro invariable.
93
INGENIERÍA MECÁNICA UNSAAC
Figura 92. Variación de la Frecuencia Natural en función del incremento de K1. (Fuente: Informe interno
Makyl)
Figura 93. Variación de la Frecuencia Natural en función del incremento de K2. (Fuente: Informe interno
Makyl)
94
INGENIERÍA MECÁNICA UNSAAC
Figura 94. Variación de la Frecuencia Natural en función del incremento de K3. (Fuente: Informe interno
Makyl)
Análisis del modelo discreto de 4 Grados de Libertad. Como se vio anteriormente se
obtuvo un procedimiento para hallar las frecuencias naturales del sistema en el plano XZ, de
una manera que se tiene un modelo representativo del problema de resonancia en la estructura
que soporta los polines de una faja transportadora Stacker Radial.
m1 m2 m3 m4
x
Apoyo Ball Ring Apoyo Travel Gear
z
L1 L2 L3 L4 L5
Figura 95. Vista en el plano XZ del Stacker Radial. (Fuente: Informe interno Makyl)
Entonces según lo mencionado en el capítulo anterior debemos aplicar el método de
flexibilidades para poder hallar la matriz de rigidez y aplicar el análisis modal y obtener las
frecuencias naturales del sistema.
Por lo tanto, luego de realizar todas las operaciones y obtener todos los elementos de
la matriz de flexibilidad, se obtiene la matriz de rigidez:
95
INGENIERÍA MECÁNICA UNSAAC
m1 m2 m3 m4
L1 L2 L3 L4 L5
Figura 96. Modelo conceptual para el análisis. (Fuente: Informe interno Makyl)
Se utilizaron las masas del sistema de las tablas que se vio en el capítulo anterior para
el caso del que el Stacker trabaje sin carga y con carga. Obviamente el más relevante es el
caso en que la faja transportadora se encuentre con carga. Luego se aplica el análisis modal
como vimos en las ecuaciones del capítulo anterior.
M i x i K ij xi 0
Como vemos en la ecuación característica del sistema, no podemos obtener las
frecuencias naturales directamente ya que no podemos dividir las matrices de rigidez y masas
por ser matrices. Por lo cual el espacio modal nos servirá de herramienta para poder hallar las
frecuencias naturales y los respectivos modos de vibración.
96
INGENIERÍA MECÁNICA UNSAAC
Con Sin
Modo
carga carga
Modo
0,72 0,78
1
Modo
1,94 2,15
2
Modo
4,00 4,42
3
Modo
9,83 9,99
4
De esta forma podemos presenciar el poder que tiene realizar el cálculo analítico
mediante un modelo conceptual ya que este nos permite replicar el problema y dar un primer
vistazo a cuál es el problema a enfrentar. Esto nos permite proyectar cuales podrían ser las
posibles soluciones que se deben implementar para que el sistema salga de la zona de
resonancia. Una vez obtenido el análisis modal, se procedió a realizar un análisis dinámico.
Para este análisis se aplica una fuerza puntual en la masa m1, la cual representa la masa de
desbalance producida por el desalineamiento angular entre el eje de salida del reductor y el
eje de entrada de la polea motriz. Esta fuerza fue calculada multiplicando la constante de
rigidez K1 por el desplazamiento de la masa M1 medido en campo. Mediante el análisis
matricial se logra conocer los desplazamientos, dado que tienen una naturaleza armónica,
solucionando el sistema definido por:
97
INGENIERÍA MECÁNICA UNSAAC
AMPLITUD DE VIBRACIÓN
(mm) A 1,913 Hz CASO ACTUAL
Sin Con
Masa
Carga Carga
M1 2.4 15.3
M2 4.8 36.2
M3 2.2 18.7
M4 1.9 13.4
Donde:
A través de este análisis se obtiene las amplitudes logradas tanto para el trabajo sin
carga como para el caso con carga completa cuando la faja opera a una frecuencia de trabajo
de 1.9133 Hz. Los resultados obtenidos en este análisis fueron:
Para poder obtener resultados se requiere tener la fuerza de excitación que está
induciendo al sistema a vibrar. En este caso se dedujo que era debido a un desbalance entre la
salida del reductor y el acople de la polea motriz. En otros casos puede haber diferentes
fuentes de excitación de las cuales el ingeniero debe deducir de donde viene el problema.
98
INGENIERÍA MECÁNICA UNSAAC
FUNCION DE AMPLIFICACION DE M1
20
18
16
14
12
AMPLITUD (mm)
10
4
CON CARGA
SIN CARGA
0
1 2 3 4 5 6 7 8 9 10 11 12
FRECUENCIA (Hz)
Figura 97. Función de amplificación de M1. (Fuente: Informe interno Makyl)
.
99
INGENIERÍA MECÁNICA UNSAAC
FUNCION DE AMPLIFICACION DE M2
20
18
16
14
12
AMPLITUD (mm)
10
CON CARGA
SIN CARGA
2
0
1 2 3 4 5 6 7 8 9 10 11 12
FRECUENCIA (Hz)
Figura 98. Función de amplificación de M2. (Fuente: Informe interno Makyl)
FUNCION DE AMPLIFICACION DE M3
20
18
16
14
12
AMPLITUD (mm)
10
CON CARGA
SIN CARGA
2
0
1 2 3 4 5 6 7 8 9 10 11 12
FRECUENCIA (Hz)
100
INGENIERÍA MECÁNICA UNSAAC
FUNCION DE AMPLIFICACION DE M4
20
18
16
14
12
AMPLITUD (mm)
10
2
CON CARGA
SIN CARGA
0
1 2 3 4 5 6 7 8 9 10 11 12
FRECUENCIA (Hz)
Figura 100. Función de amplificación de M4. (Fuente: Informe interno Makyl)
A partir de estos resultados se pueden proponer soluciones para atenuar los niveles de
vibración que están afectando a la estructura. Las ideas pueden variar en función de la
creatividad. Por ejemplo, una de las soluciones que propuso inicialmente Makyl pero que no
se implemento fue diseñar un absorsor de vibraciones como se ve en la figura.
ma
ka m2 m3 m4
x m1 I1 I1 I2 I2 I2
Figura 101. Idea conceptual de una posible solución mediante absorsor. (Fuente: Informe interno Makyl)
Análisis por elementos finitos. En esta etapa luego de haber entendido el problema a través
del modelo conceptual, se pasa a analizar la estructura a través de los elementos finitos.
Análisis modal FEM modelo discreto actual. Es necesario validar nuestro modelo
conceptual a través del análisis por elementos finitos y así verificar que no se haya cometido
101
INGENIERÍA MECÁNICA UNSAAC
algún error al momento de hallar las frecuencias naturales en nuestro modelo conceptual. A
continuación, se ven los resultados de los 2 primeros modos de una simulación FEM del
modelo discreto llevado a cabo en Femap.
Figura 102. Análisis modal del modelo discreto. (Fuente: Informe interno Makyl)
102
INGENIERÍA MECÁNICA UNSAAC
Figura 103. Frecuencia de Resonancia a 0.53Hz de la Estructura en el plano XZ. (Fuente: Informe interno
Makyl)
103
INGENIERÍA MECÁNICA UNSAAC
Figura 104. Frecuencia de Resonancia a 1.68Hz de la Estructura en el plano XZ. (Fuente: Informe interno
Makyl)
Figura 105. Frecuencia de Resonancia a 3.96Hz de la Estructura en el plano XZ. (Fuente: Informe interno
Makyl)
Figura 106. Frecuencia de Resonancia a 3.96Hz de la Estructura en el plano XZ. (Fuente: Informe interno
Makyl)
Implementación de la solución
104
INGENIERÍA MECÁNICA UNSAAC
Figura 107. Vista del Stacker Radial con soluciones implementadas. (Fuente: Informe interno Makyl)
105
INGENIERÍA MECÁNICA UNSAAC
Figura 109. Pórticos de anclaje superior e inferior. (Fuente: Informe interno Makyl)
Se propone el reforzamiento de los soportes de los polines superiores e inferiores
(bastidores) que otorga una mayor rigidez, de manera que las deflexiones de estos soportes
estén dentro de norma y que sus frecuencias naturales no se vean activadas por el ancho de
banda de la frecuencia de operación de los polines, evitando así el fenómeno de resonancia y
sus consecuencias estructurales como son solturas, roturas, desprendimientos y fatiga en
componentes.
106
INGENIERÍA MECÁNICA UNSAAC
Figura 110. Reforzamiento de los polines inferiores y superiores. (Fuente: Informe interno Makyl)
Figura 111. Masas incorporadas entre el Ball Ring y el Travel Gear. (Fuente: Informe interno
Makyl)
107
INGENIERÍA MECÁNICA UNSAAC
Figura 112. Masas adicionales agregadas a la estructura. (Fuente: Informe interno Makyl)
Figura 113. Reforzamiento del sistema motriz. (Fuente: Informe interno Makyl)
Análisis modal al sistema reforzado
Se han resaltado aquellos modos que tienen mayor masa modal, por lo que involucran
desplazamientos generalizados en toda la estructura.
108
INGENIERÍA MECÁNICA UNSAAC
109
INGENIERÍA MECÁNICA UNSAAC
CONCLUSIONES
Se estudió y comprendió la teoría que conlleva el análisis modal, ya sea el
experimental y operacional, llegando a concluir que son ensayos muy importantes
y requeridos en la industria y la minería, debido a su gran implicancia en
operación de grandes maquinarias.
SUGERENCIAS
Se recomienda implementar un laboratorio de vibraciones mecánicas y ensayos
dinámicos, en nuestra escuela profesional, debido a su creciente requerimiento en
el mercado y la falta de profesionales ilustrados en el tema.
110
INGENIERÍA MECÁNICA UNSAAC
PLANOS
111
INGENIERÍA MECÁNICA UNSAAC
BIBLIOGRAFIA
Alencastre Miranda, J. (2018). Material del curso de Vibraciones Mecánicas y
Dinámica Estructural. Lima.
Meruane, V. (2012). Vibraciones Mecánicas. Santiago, Chile.
Singiresu, R. (2010). Vibraciones Mecánicas. Prentice Hall.
112