Cálculo de Carga en Molinos de Bolas
Cálculo de Carga en Molinos de Bolas
PROYECTO DE TESIS
Huacho - Perú
2012
pág. 1
ÍNDICE GENERAL
Carátula. 01
Índice. 02
Índice de figura. 05
Índice de cuadros 06
Glosario abreviatura. 07
Índice de anexo. 08
I. PLANTEAMIENTO DEL PROBLEMA.
1.1. Antecedentes. 09
1.2. Formulación del Problema. 10
1.2.1. Problema General. 10
1.2.2. Problema Específico. 10
1.3. Objetivos. 10
1.3.1. Objetivos Generales. 10
1.3.2. Objetivos Específicos. 10
1.4. Justificación. 11
1.4.1. Justificación Práctica. 11
1.4.2. Justificación Metodológica. 11
1.4.3. Justificación Técnica. 11
1.4.4. Justificación Social. 11
1.4.5. Justificación Económica. 12
1.5. Delimitación. 12
1.5.1. Delimitación Territorial. 12
1.5.2. Delimitación Tiempo y Espacio. 12
1.5.3. Delimitación de Recursos. 12
1.6. Viabilidad de Estudio. 13
II. MARCO TEÓRICO.
II.1. Antecedentes dela Investigación. 14
II.1.1. Investigación Relacionada con el Estudio. 14
II.1.2. Otras Publicaciones. 19
II.2. Bases Teóricas. 22
II.2.1. Molienda. 22
II.2.2. Tipos de Molinos. 23
II.2.3. Clases de Molienda. 23
II.2.4. Variables en el procesamiento de Minerales. 24
pág. 2
II.2.5. Variable de diseño. 25
II.2.6. Variable de Operación. 25
II.2.7. Relación energía – tamaño. 27
II.2.8. Antigua leyes de la conminución. 27
II.2.9. Postulado de Bond. 28
II.2.10.Diseño Factorial Simple. 30
II.2.11.Cálculo de Efectos. 31
II.2.12.Análisis de Variancia. 32
II.2.13.Modelo Matemático a Escala Codificada. 33
II.2.14.Análisis de Residuos. 34
II.2.15.Decodificación del Modelo Matemático a Escala Natural. 35
II.3. Definición de Términos Básicos. 35
II.4. Formulación de Hipótesis. 36
II.4.1. Hipótesis General. 36
II.4.2. Hipótesis Específicas. 36
II.5. Operacionalización de Variables. 37
III. METODOLOGÍA.
III.1. Diseño Metodológico. 37
III.1.1. Tipo de Investigación. 37
III.1.2. Método de Investigación. 37
III.1.3. Estrategias o Procedimientos de Contrastación de Hipótesis o
Cumplimiento de Objetivos Técnicos. 37
III.2. Población y Muestra. 38
III.2.1. Población. 38
III.2.2. Muestra. 38
III.3. Técnica de Recolección de Datos. 38
III.3.1. Técnicas. 38
III.3.2. Instrumentos. 39
III.4. Técnica de procesamiento de la Información. 39
III.5. Aspectos Éticos. 39
III.6. Matriz de Consistencia. 40
IV. RECURSOS, PRESUPUESTOS Y CRONOGRAMA.
IV.1. Recursos. 42
IV.1.1. Recursos Humanos. 42
IV.1.2. Recursos Institucionales. 42
pág. 3
IV.1.3. Recursos Financieros. 42
IV.2. Presupuestos. 42
IV.3. Cronograma. 43
V. BIBLIOGRÁFIA Y OTRAS FUENTES.
Anexos.
pág. 4
ÍNDICE FIGURA
Pág.
Fig. Nº 2-1 : Variación de la velocidad critica del molino 26
Fig. Nº 2-2 : Representación geométrica de un diseño, para k = 3 replicado
en el centro. 31
pág. 5
ÍNDICE CUADROS
Pág.
Cuadro Nº 2-1 : Variables de operación 37
Cuadro Nº 2-2 : Matriz de consistencia 1. 40
Cuadro Nº 2-3 : Matriz de consistencia 2. 41
Tabla Nº 2-4 : Presupuestos. 42
Cuadro Nº 2-5 : Cronograma 43
pág. 6
GLOSARIO ABREVIATURA
% : Porcentaje
UNJFSC : Universidad Nacional José Faustino Sánchez Carrión.
Wi : Word índex (índice de trabajo )
S.A : Sociedad anónima
Kg/Ton : Kilogramo por tonelada
TMSD : Toneladas métricas secar por día
Kw-h/T : Kilo watts hora por tonelada corta
Nc : Velocidad crítica
Kw : Kilo watts
W : watts
m : Micrones
P80 : 80 por ciento pasante del producto
F80 : 80 por ciento pasante en la alimentación
HPGR : Molienda de Rollos a presión alta
Pág. : Pagina
Fig. : Figura
A- : Anexo
EAPIM : Escuela académica profesional ingeniería metalúrgica.
pág. 7
ÍNDICE ANEXO
pág. 8
ÍNDICE DE TABLA
Tabla Nº 2-1 : 100
pág. 9
ÍNDICE DE FIGURA.
Fig. Nº 2-1 : 101
Fig. N 3-10 :
Fig. N 3-11 :
pág. 10
SÍMBOLOS
F80 : Ochenta porciento pasante alimentación.
HPGR : High PressureGrindingRolls (Molienda de Rollos a presión alta).
Kg/Ton : Kilogramos por toneladas.
Kw h/T : Kilo watt hora por tonelada.
P80 : Ochenta porciento pasante producto.
Pag. : Página.
S.A : Sociedad Anónima
Wi : Índice de trabajo
Vc : Velocidad crítica tangencial.
UNJFSC : Universidad Nacional José Faustino Sánchez Carrión.
TMSD : Toneladas métricas secas día.
SAG : Molienda semiautógena.
MINEX : Minería y Exportaciones.
pág. 11
CAPITULO I
ASPECTOS PRELIMINARES
pág. 12
Las principales pérdidas de energía están dadas, enforma de calor y ruido, en las
transmisiones delaccionamiento del molino.
pág. 13
Realizando la disgregación de minerales sulfurado de cobre experimentalmente a nivel
laboratorio nos permitirá optimizar la molienda usando el diseño factorial en la
Universidad Nacional José Faustino Sánchez Carrión 2012.
Para el tratamiento de los datos se aplicará la técnica del diseño factorial simple a fin
de procesar e identificar la influencia de las variables en estudio.
En el análisis de los datos se aplicará el software estadístico SPSS versión 19, a fin de
evaluar, interpretar y optimizar las variables en estudio.
pág. 14
A fin de desarrollar el presente trabajo de investigación se han estudiado las siguientes
variables:
Todos los pequeños mineros que no cuente con planta de beneficio en el distrito de Nazca
provincia de Nazca departamento de Ica.
Con esta técnica será posible conocer la forma como se desarrollan el análisis y
disgregación de minerales sulfurado de cobre a nivel laboratorio para la optimización de la
molienda
pág. 15
Se usará el análisis estadístico, usando programas de cálculo como Excel, SPSS, para
luego mostrar la información, mediante tablas, registros, figuras, promedios, medianas,
desviación estándar y otros.
pág. 16
CAPITULO II
MARCO TEÓRICO
2.2.1. Molienda.[2,5,21]
La molienda por atrición depende del grado de resbalamiento que existe entre las
zonas de desprendimiento y la base de carga. Cuando el diseño de los forros es bueno, la
transferencia de energía del casco a la carga de bolas, es igualmente buena; como si fuera un
conjunto de engranajes internos, por lo que el efecto de impacto es máximo. Si el grado de
resbalamiento de bolas sobre el casco es alto; se incrementa el efecto de molienda por
abrasión o atrición y se originará una mayor cantidad de lamas, así como un producto más
fino; pero es a expensas de la potencia y desgaste de forros.
pág. 17
a) Liberación de las partículas valiosas del mineral.
b) Llevar la granulometría de las partículas a un tamaño adecuado para la siguiente etapa.
c) Homogenizar la pulpa para la siguiente etapa.
Cuando se trabaja en circuito cerrado quiere decir que se trabaja con un clasificador para
incrementar la eficiencia de la molienda.
o Molinos de barras.
o Molinos de bolas (Molino de cónicos, Horizontales y Verticales).
o Molinos autógenos
o Molino semi autógenos.
o Molinos de guijarros
pág. 18
d. Molienda Diferencial: Denominada selectiva, aplicable a productos estructuralmente
heterogéneos, permite con ciertas precauciones aumentar la reducción de los
componentes fiables, algo mas que la de los componentes mas duros. De esta manera
es posible una selección más o menos grosera mediante clasificación volumétrica. Los
resultados de la molienda diferencial son esencialmente función de la naturaleza física
del material a tratar.
e. Molienda Formal: En donde las operaciones de fragmentación en las cuales la forma de
los fragmentos obtenidos influye en, el valor del producto.
f. Automolienda: Trata de operaciones de reducción volumétrica en las cuales los
elementos a moler se fragmentan por percusión y fricción recíprocas.
g. Molienda Criogena: Fragmentación de sólidos previamente enfriados a muy bajas
temperaturas (en general mediante nitrógeno liquido).
o Carga de mineral.
o Suministro de agua.
o Medios de molienda.
o Condición de los blindajes.
o Tiempo de molienda.
o Velocidad de operación de los molinos.
o Carga Circulante.
pág. 19
eficiencia. En estacondición se impiden los efectos de cataratas ycascadas de los medios
moledores, que son losefectos de los que depende la molienda.
b) Volumen de la carga: El volumen de la carga en un molino de bolasse expresa como el
porcentaje del volumen entrelas corazas que es llenado con bolas y mineral.Cuando el
molino está detenido, el volumen decarga puede obtenerse en forma rápida, midiendoel
diámetro interno entre corazas y la distanciadesde la parte superior del molino y la
carga.
c) Potencia versus carga en el molino. La máxima potencia se consume en el caso deque
la carga ocupe aproximadamente el 50 % del volumen.
d) Tamaño del molino: El tamaño del molino se determinará sólo enbase a la potencia
requerida para moler. Laestimación del tamaño puede realizarse en base afórmulas
empíricas.
pág. 20
b) Porcentaje de sólido en el molino.
c) Tamaño de Bolas.
Cuando una máquina opera sobre un mineral consume una cierta potencia, la que puede
variar en el tiempo dependiendo del esfuerzo que ella deba realizar en cada instante.
La energía unitaria utilizada por esa máquina para efectuar su trabajo se puede expresar
de la siguiente forma:
pág. 21
Si estamos en el estado estacionario: E=Pt ( KWH )
t
t
P (ξ )
E =∫
Por regla general especificamos como energía especifica: ^ d (ξ )
0 W (ξ )
Donde: P ( t ) =potencia en KW
W ( t )=carga de material en el equipo(TM )
t = tiempo en horas
E=Pt /W
^ con ( x ), de tal manera que a cierta granulometría, cuanta
La idea es relacionar E
energía se ha consumido.
pág. 22
- Vencer la energía de las piezas de las máquinas.
- Deformaciones elásticas de la máquina.
- Producción de ruido, calor y vibraciones de la instalación.
- Generación de electricidad.
- Roce entre partículas y piezas de la máquina.
- Pérdidas de eficiencia en la transmisión de energía eléctrica y mecánica.
Bond basó su Tercera Ley de la Conminución en tres principios de los que a su vez se
basan en mecanismos observados durante la reducción de tamaño de las partículas.
pág. 23
concluir que la energía consumida es proporcional a la diferencia entre la raíz
cuadrada de la superficie específica obtenida después y antes de la conminución. Esto
es:
1 1
[
E B=C ( S^ 1 ) 2 −( S^ 0 ) 2
^ ]
Reemplazando la superficie especifica en términos del tamaño promedio, resulta:
1 1
1 1
E B=K B
^
[( ) ( ) ]
x 1/122 1
2
−
x 1/122 0
2
1
α 1
Donde:K B=C S
αV ρ [ ] 2
1 1
E B=K B
^
[ −
√ P80 √ F 80 ]
F. Bond definió el parámetro K B en función del Work Index. W i (Índice de Trabajo del
material), que corresponde al trabajo total (expresado en KWh/Tc) necesario para
reducir una tonelada corta de material desde un tamaño teóricamente infinito ( F 80 → ∞ )
hasta partículas que en un 80% sean inferiores a 100 micrones
( P¿¿ 80=100 μm , aproximadamente67 %<200 mallas)¿;). Entonces:
1 1 KB
W i =K B [√ − =
100 √ ∞ 10 ]
De donde: K B=10 W iy finalmente, al reemplazar, tenemos:
pág. 24
1 1
E B=10 W i
^
[ −
√ P80 √ F 80 ]
Donde:
P80 = tamaño 80% pasante del producto (μm).
F80 = tamaño 80% pasante de la alimentación (μm).
WI = Índice de Trabajo del Material (KWh/Tc).
En el diseño factorial simple se asume que la relación de las variables con la respuesta
está representada por un modelo matemático lineal, el cual debe de predecir adecuadamente
todos los puntos del diseño, inclusive el punto central del diseño, por seguridad de que no se
cumpla esta hipótesis podemos afirmar que el modelo matemático lineal no es suficiente para
explicar las respuestas especialmente en el punto central del diseño y podemos asumir una
posible existencia de curvatura. Para ello la adición de las pruebas en el punto central del
diseño nos sirve para dar estimación promedia del efecto curvatura.
pág. 25
Fig. 2-2: Representación geométrica de un diseño, para k = 3 replicado en el centro.
El cálculo de efectos nos determinara cuando afecta las variables sobre la respuesta o
llamado también criterio de optimización de un proceso. El cual se define como la diferencia
entre los valores medios de las respuestas, cuando en los tratamientos respectivos la variable
se encuentra en su nivel inferior y cuando se encuentra en su nivel superior, definido por:
Efecto=
∑ Y O −∑ Y
N
2
Donde:
pág. 26
- SS Total = suma total de cuadrados corregidos.
- SS Efectos = suma total de cuadrados debido a los efectos o tratamientos.
- SS Error = suma total de cuadrados debido al error.
N 2
SS . . Efecto=
(∑ )
i=1
XijYi
b. La suma total de cuadrados para el error (SS Error) está definido por:
no 2
SS . . Error=
(∑ )
i=1
Yº −Yº
Donde:
Yº= Réplicas en el punto central del diseño.
no= numero de réplicas en el centro del diseño.
Yº= promedio de todas la réplicas.
2
N∗no ( Y −Yº )
SS . . Curvatura=
( N +no)
pág. 27
La manera precisa para determinar estadísticamente las variables o factores que
pueden ser significativos es por medio del teorema de Cochran el cual se resume en la
siguiente expresión:
SS .. Efecto
f1
Fo=
SS .. Error
f2
- f1= grados de libertad de los efectos e interacciones generalmente iguales a 1.
- f2= grados de libertad de la suma de cuadrados del error, lo cual está calculado
como la diferencia de los grados de libertad de la suma total de cuadrados menos
la suma de los efectos.
- f3= grados de libertad de la suma de cuadrados definido como el número de
pruebas experimentales realizadas (incluyendo las réplicas) menos 1.
[ B ] =([ X ]T [ X ])−1 ([ X ]T [Y ])
Con estos elementos podemos calcular la suma de cuadrados de residual del modelo
mediante la siguiente ecuación:
pág. 28
Nr
(Y −Yº )2
SS . . Re siduales=∑
l =1 Nr−l
Donde:
Y = Y observado.
Yº = Y calculado.
SS ..Re siduales
Fo=
Media. . decuadrados. . del. . error
pág. 29
Para el término lineal:
k =1
bj dij
ajZj = −∑
Δ Zj ji=1 Δ Zj
pág. 30
pág. 31
CAPITULO V
PARTE EXPERIMENTAL
a) Chancadora de quijada
b) 01 molino de bolas para efectuar la molienda de minerales de capacidad de 01 Kg de
muestra 8”X8”.
c) Estufa para secar muestra.
d) Balanza digital hasta 05 Kg de capacidad.
e) Cuarteador jones.
f) Rote up.
g) Filtro prensa.
pág. 32
X1 Densidad de pulpa (gr/mL).
X2 Tiempo de molienda (minutos).
X3 Porcentaje de carga de molienda.
A base de los datos anteriores se construye la matriz de diseño a escala codificada y a
escala natural para el diseño factorial simple:
Tabla Nº 5-x: Matriz del diseño factorial Escala codificada y natural molienda de mineral.
N Y
1 69.62
2 69.89
3 80.81
Tabla Nº 5-x: Matriz del diseño factorial con réplicas en el centro y respuesta de molienda.
N Y
1 69.62
pág. 33
Es el procedimiento apropiado para probar la significancia de los efectos. La siguiente
tabla muestra el análisis de varianza para los datos.
N Y
1 69.62 94.00
2 69.89 94.00
3 80.81 94.60
Ecuación general:
La ecuación codificada:
X0 X1 X2 X3 X1X2 X1X3 X2X3 X1X2X3
65,30 5,511 -0,509
pág. 34
Para ello multiplicamos la matriz de variables independientes con las columnas X2, X1X3
que resulto no significativa, y Y obsson los datos observado de las pruebas realizadas y se calcula
con el modelo matemático para producir Ycorrcalculamos los valores predichos por el modelo.
SSMR= 0,736152932
Fo= 11,45708926
pág. 35
5.7 DESCODIFICACIÓN DEL MODELO MATEMÁTICO A ESCALA NATURAL Y SU
ANÁLISIS.
Como se ha estado trabajando las variables de X1, X2, X3, representada por los valores de
(-1) y (+1), para que en el modelo se deba remplazar los valores reales se debe descodificar.
a0= 51,29
a2= 2,204
a3= -0,058
Y
= a0 a1Z1 a2Z2 a3Z3 a4Z1Z2 a5Z1Z3 a6Z2Z3 a7Z1Z2Z3
Y 2,204
= 51,29 5 -0,0576
pág. 36
De la ecuación podemos se observa que el tiempo tiene una influencia positiva en la
liberación del mineral, mientras que la relación de densidad y porcentaje de carga moledora
tiene una influencia negativa en la recuperación.
Mínimo Máximo
Porcentaje liberación (%) 51,46 57,59
Tiempo (min.) 10,00 15,00
58
% Libración -200m
56
54
Tiempo (min)
52
50
48
9.0 10.0 11.0 12.0 13.0 14.0 15.0 16.0
Tiempo (min)
pág. 37
58
58
57
% Libración -200m 57
56
56 mineral +bolas
55
55
54
300 320 340 360 380 400 420 440 460 480
Mineral + bolas
0 2 4 6 8 10 12
pág. 38
pág. 39
3. BIBLIOGRAFIA Y OTRAS FUENTES:
pág. 40
15. Manual de minería (2010). Estudios Mineros del Perú S.A.C. extraído el 05 de octubre del
2011de http:// www.estudiosmineros.com.
16. Manzaneda Cábala, J. (1990). "Procesamiento de Minerales " Lima: Cáceres.
17. Menéndez Aguado, J.M. (2001).Aplicación de la Simulación Matemática a la
Determinación de Consumos Energéticos en Fragmentación. Tesis Doctoral, Dpto.
Explotación y Prospección de Minas, Universidad de Oviedo España.
18. Misra, V. N., Das S. &Subbaiah, T. (2005) Emerging trends in mineral processing and
extractive metallurgy. ISBN 81-7764-858-6 :Alliedpublishers PVT. LTD.
19. Montgomery, D. (1991). Diseño y Análisis de Experimentos. México: Iberoamericana.
20. Mosquera Santos, E. V. (1983). Ampliación y optimización Metalúrgica de la planta
concentradora en la minera el Brocal S.A. Tesis Ingeniero Metalúrgico, UNMSM,Lima
Perú.
21. Porras Castillo, D. L. (1995). Procesamiento de Minerales. Lima: San maros.
22. Quiroz Núñez, I. (1986). Ingeniería Metalúrgica Operaciones Unitarias en
Procesamiento de Minerales. Lima: San Marcos.
23. Rao G. V. &MisraVibhuti, N. (2004). Mineral Processing Technology. : Allied publishers
PVT. LTD.
24. Rodríguez Avello, Á. (1981). El Beneficio de los Minerales. Madrid: Rocas y Minerales.
25. SEPÚLVEDA Jaime E. (2006). Criterios para la selección, aplicación y evaluación de
medios de molienda. I Encuentro Internacional de Procesamiento de
Minerales,Antofagasta Chile. En: www.molycop.cl.
26. Silva Díaz, E. (1988). Evaluación y optimización de las operaciones en la concentradora
de la minera Raura. Tesis Ingeniero Metalúrgico, UNMSM. Perú.
27. Torres Rubilar, M. E. (2010). Modelación y simulación de molienda con rodillos a alta
presión – HPGR. Tesis Magister en Ingeniería metalurgia Extractiva,Facultad de Ciencias
Físicas y Matemáticas, Universidad de Chile,Santiago de Chile. En:
www.cybertesis.uchile.cl/tesis/uchile/2010.
28. Vargas Gonzales, J. W. (2010). Evaluación de un circuito de molienda y clasificación.
Tesis Ingeniero de minas. Facultad de ingeniería en ciencias de la tierra. Escuela Superior
Politécnica del Litoral. Guayaquil Ecuador. En: www.dspace.espol.edu.ec.
pág. 41
30. Venugopal, R., Sharma, T., Saxena, V.K., &Mandre, N.R.(2005). Mineral Processing
Technology. India: McGraw-Hill.
31. Wills, B. A. & Napier Munn, T. J. (2006). Mineral processing technology. Queensland:
ISBN-10: 0-7506-4450-8. En www.elsevier.com.
32. YoplacCastromonte, E. (2008).Modelos matemáticos para reducción de minerales y
flotación. Lima: www.ciplima.org.pe.
33.
pág. 42
ANEXOS
pág. 43
Fig. A-1: Mapa de localización dela UNJFSC EAPIM.
pág. 44
Fig. A-2: Mapa Nacional de localización Huacho.
pág. 45
Fig. A-3: Mapa localización Nazca - Planta Concentradora MINEX.
pág. 46