Trabajo Final
Trabajo Final
Trabajo Final
RESOLUCION DE PROBLEMAS DE
INGENIERA METALURGICA
1
INDICE
1.- Sumario 3
2.-Introduccion 5
3.1 .- Problema 1 15
3.1.1 .-Método de Newton Raphson multivariable 15
3.1.2.- Método del punto fijo sucesivo 18
3.1.3.- Método del punto fijo simultáneo 19
3.2 .- Problema 2 21
3.3 .- Problema 3 25
3.3.1 .- Buscar objetivo 25
3.3.2 .- Solver sin restricciones 26
3.3.3 .- Solver con restricciones 27
3.4.- Problema 4 29
3.5.- Problema 5 37
5.- Bibliografía 47
2
1. SUMARIO
Este trabajo tiene como finalidad resolver problemas de ingeniera principalmente asociados a
la metalurgia, aplicando diversos métodos de resolución y así comparando los resultados
conseguidos predecir que método es más eficaz.
Los problemas a solucionar en este trabajo son cinco. El primero trata de un molino de bolas
donde se requiere determinar el diámetro del molino y la fracción de llenado de las bolas
mediante un sistema de ecuaciones no lineales, aquí el objetivo es demostrar los
conocimientos que adquirimos de los distintos métodos numéricos para la resolución de
sistemas de ecuaciones tantos lineales y no lineales. Para ello utilizamos los métodos
siguientes:
Donde se obtuvo que el método que convergía con mayor rapidez es el de MNRMV.
El segundo problema trata del traslado de los ánodos de cobre y nos requiere determinar la
potencia suministrada por la fuente de circuito. Para ellos aplicamos leyes de Kirchhoff
obteniendo un sistema de ecuaciones lineales y utilizamos el siguiente método
El tercer problema nos pide determinar la distancia entre las mandíbulas de una trituradora
para ello utilizamos las herramientas de Excel:
Buscar Objetivo
Solver
Acá el solver es un método más eficiente ya que se pueden añadir restricciones los cuales
logran hacer más precisa la solución.
El cuarto problema tiene relación a la termodinámica, se requiere determinar que metal logra
reducir a un oxido de níquel que se forma en un reactor, para ello debemos calcular su
variación de energía estándar utilizamos los siguientes métodos:
Regla trapezoidal
Regla de Simpson 1/3
Regla de Simpson 3/8
Si bien todos los métodos se acercaron bastante a la solución real no poseemos de un criterio
para determinar cual método es más eficaz debido a que estos métodos son exactos para
polinomios a lo más de tercer grado.
3
Y finalmente el quinto problema es una EDO que consiste en calcular la cantidad de
contaminante en un tanque de relave a un determinado tiempo. Para ello se utilizo:
Para los métodos anteriormente nombrados se debió desarrollar los teoremas y algoritmos
concernientes y con el apoyo del programa Excel y Visual Basic que nos permiten programar y
crear plantillas que resuelvan la ecuación y los sistemas de ecuaciones.
4
2. INTRODUCCION
Desplazamientos simultáneos
Sumando x a ambos lados de cada una de las ecuaciones, se forma un nuevo sistema de
ecuaciones:
Desplazamientos sucesivos
Criterio de convergencia:
5
, i= 1,…., n. Se puede generar nuevos esquemas numéricos
a través de nuevos despejes algebraicos. Sin embargo, no siempre se obtiene una sucesión
convergente del MPFMV. Aún para una sucesión convergente, es posible que la convergencia
dependa de valores iniciales de las iteraciones.
En general, no existe un criterio de convergencia que involucre todos los casos posibles
de esquemas numéricos del MPFMV. Pero con la base de obtener una sucesión convergente de
iteraciones, se puede decir el siguiente criterio que entrega condiciones suficientes para la
convergencia del MPFMV;
6
Toda ecuación con , puede ser no lineal o lineal, sin embargo se debe
tener al menos una ecuación no lineal. Luego, el sistema no lineal mostrado anteriormente, se
resuelve de la siguiente forma:
Algoritmo:
Teorema:
La cual posee solución única si y solo si es una matriz NO SINGULAR, en donde una matriz no
singular debe cumplir las siguientes condiciones:
7
es invertible:
todas las filas (y columnas) de son L.I.:
Entonces, si
Así,
Algoritmo:
Herramientas de Excel
Excel para realizar búsquedas consta de dos herramientas las cuales son:
Buscar Objetivo:
Esta es una herramienta rápida y fácil de usar. La cual nos permite mediante un algoritmo de
iteración interna del Excel, definir una formula con un valor deseado (función objetivo), a
través de cambios en otra celda (parámetro), que interviene en la fórmula. Por ejemplo se
quiere resolver la ecuación: por ejemplo
8
iii) Se cambia el valor de la celda, método por el cual se lograra mediante la definición
de un valor inicial de x, cuyo valor debe ser lo mas cercano posible a la solución
Solver:
Esta herramienta al igual que la búsqueda de objetivo es muy fácil y rápida de usar, pero la
diferencia es que esta herramienta nos permite agregar condiciones o restricciones a la
búsqueda.
Con esta herramienta se puede obtener el valor óptimo para una formula y además ajusta los
valores en las celdas cambiantes que se especifiquen (celdas ajustables), para así generar el
resultado especificado en la fórmula de la celda objetivo. Más aún se pueden aplicar
restricciones para la restricción de los valores que pueden utilizar Solver y estas restricciones
pueden hacer referencia a otras celdas a las que afecte la fórmula de la celda objetivo.
Integración numérica
9
Ésta, puede aproximar el integrando por el polinomio que interpola a en
Integrando se obtiene:
Luego, haciendo un análisis geométrico, nos damos cuenta que la regla aproxima el área bajo la
curva f(x) mediante un trapecio inscrito.
Cabe destacar que la regla del trapecio se puede aplicar de dos maneras: un solo trapecio o
varios trapecios, siendo para las dos la misma aplicación, sólo lo que varía es el número de
10
intervalos, de lo que se puede observar que a mayor número de intervalos, más precisión tiene
el método.
Algoritmo:
Regla de simpson
Y obtenemos,
Algoritmo:
11
ii).- Aplicamos la aproximación como corresponde y
Algoritmo:
Una ecuación diferencial ordinaria es una ecuación que contiene las derivadas de una o
más variables dependientes con respecto a una o más variables independientes en donde se
clasifican de acuerdo a su tipo, orden y linealidad.
12
El cual supondremos que tiene solución única, donde , la cual es
Algo muy importante de los PVI es que se basan en tomar una partición en N subintervalos
[a,b] ( ) y obtener sucesivamente N números que
Para obtener las soluciones del PVI existen dos formas habituales que son los métodos
analíticos y los métodos numéricos:
Los métodos analíticos generan soluciones exactas que se resuelven por edos exactas, factor
integrante, Bernoulli entre otras que son vistos en cursos anteriores. Nuestro objetivo es
abordar los métodos numéricos que aunque se obtienen soluciones aproximadas permiten
entregar soluciones rápidas y con un grado de aproximación requerido.
obtener es decir, se necesitan conocer los puntos anteriores para obtener los
siguientes puntos. Considerando el PVI consiste en transformar la derivada y(x) ’en una
pendiente lineal a través de la aproximación:
13
: De donde,
Una vez calculada esta aproximación, se puede utilizar para obtener la aproximación
de , a saber,
,…, .
Método de Heun
Llamado método predictor corrector, este es un método de paso múltiple que sugiere
una mejora al método de Euler, es decir, combina un método implícito con un método explícito
donde en la primera etapa se obtiene una estimación del algoritmo explicito y luego se
14
Donde el valor de lo podemos aproximar a través el método de Euler, así:
Entonces:
Método de Runge-Kutta 44
Este método de cuarto orden utiliza una serie de Taylor de orden N=4 y consiste
principalmente en calcular la aproximación de de la siguiente manera:
Método de Ralston
15
Este método es similar al método de runge-kutta y Euler explicito, es decir, un método de paso
simple que usan la información de la solución calculada en para obtener
es decir, se necesitan conocer los puntos anteriores para obtener los siguientes puntos y tener
la solución de la EDO.
PROBLEMA 1
En un molino de bolas utilizado para reducir a polvo minerales de cobre mediante la rotación
de un tambor que está girando a 72% de la velocidad critica, contiene bolas de acero y tienen
una carga de 289,47 [ton], L= 6 [m], la potencia consumida es P=3389[kW] y la densidad de la
bola es de ρ=7,75 [ton/m3].
Determinar el diámetro del molino (D) y la fracción de llenado de las bolas de acero (J)
Ecuaciones:
Solución:
A simple vista notamos que ambas ecuaciones son no lineales generando así un sistema de
ecuaciones no lineales.
Por lo tanto para resolver este problema utilizaremos métodos de Newton-Raphson
multivariable y punto fijo con iteraciones simultáneas y sucesivas.
(1)
(2)
16
(x)=0
Jacobiano de la
Luego el Jacobiano de la en
17
Ingresándolo estos datos a Excel se tiene:
Iteramos hasta que se cumple el criterio de error, de la siguiente imagen adjunta se observa
que el método converge en la tercera iteración
18
Método del Punto Fijo con iteraciones sucesivas:
En la iteración inicial nos damos los valores de J=0.35 obteniendo el valor de D=5.949868428
como el valor está alejado de la solución real seguimos iterando hasta la iteración 40 donde el
valor de J es de 0.389 y obtenemos el valor de D=5.822586148, el que claramente es un valor
mucho más cercano a la solución y con un error de iteración menor al 0.1%. a l ingresar esto a
19
las casillas de Excel se obtiene:
Este grafico muestra la diversificación del error mediante este método en función del número
de iteraciones:
20
Método del punto fijo con iteraciones simultáneas
Definiendo un vector (J , D)=(0.4 , 6); este se evalúa en las ecuaciones iniciales, donde los
resultados de esta se utilizan para formar otro vector que será utilizado nuevamente en las
ecuaciones y así hasta que las incógnitas sean constantes.
21
En la tabla se observa que en la iteración numero 66 las incógnitas comienzan a no variar.
Conclusión: a partir de los resultados obtenidos con los métodos se concluye que el MNRMV
converge a la solución con mayor rapidez. Mientras que los métodos de punto fijo con
iteraciones sucesivos o simultáneos convergen con mayor precisión.
Discusión: el MNRMV por converger cuadráticamente hace que se reduzca el error mas
velozmente.
PROBLEMA 2
En una minera se analiza el circuito de una maquina eléctrica que se utiliza para el traslado de
los ánodos de cobre desde los moldes hacia una piscina con agua para disminuir la
temperatura. Un operador necesita calcular la potencia suministrada por la fuente del circuito.
22
Para calcular las corrientes que gobiernan la potencia suministrada por la fuente del circuito se
necesitan usar las 2 leyes de kirchoff.
Ecuaciones de mallas
Ecuaciones de nodos
0=I1-I2-I4
0=I2-I5-I3
Y “b” es el vector:
23
A simple vista podemos observar que nuestro sistema tiene solución única puesto que
tenemos 5 ecuaciones y 5 incógnitas (I 1,I2,I3,I4,I5). Sin embargo también debe cumplirse que la
matriz A sea no singular y una de las condiciones que debe tener es que el determinante de A
sea distinto de cero.
Luego para calcular el vector “x” hay que hacer una multiplicación de una matriz de 5X5 con un
vector 1X3:
Utilizando el comando MMULT tenemos la solución del sistema, es decir las corrientes del
circuito.
24
Ahora haciendo un análisis de sensibilidad para estudiar el comportamiento de la solución del
sistema y para ello cambiamos un poco el valor de la casilla que está en blanco.
25
Ahora calculamos la norma de frobenius de A y de A -1 ingresada en el Excel de circuito.
Discusión: al estar mal condicionada pueden haber varias razones: que el tamaño de la matriz
un poco grande, pocos dígitos (mala precisión). Además haciendo un análisis de sensibilidad de
la matriz A se infiere además a partir del grafico que la potencia entregada por el circuito es
directamente proporcional a la inversa del valor a11(la casilla en donde variamos los valores).
PROBLEMA 3
Calcular la distancia entre la mandíbula inferior abierta (oss) de una trituradora giratoria para
reducir 1800 Tph desde un tamaño de material de alimentación de F 80=32 [in] a un producto de
P80=3,8[in] de diámetro 3[m]. Las propiedades del material son un índice de trabajo de impacto
de Wi=14,5 [kWh/Ton] y densidad de 2700 Kg/m 3.Las constantes K1=0,9 y una porosidad es
igual a 0,4.
26
La resolución de este problema se hace mediante la siguiente formula
En donde,
Este problema lo podemos resolver mediante las planillas de Excel las cuales se mostraran a
continuación:
27
Solver sin restricción
la misma forma que en el caso anterior, situamos los datos dados en la planilla Excel y en la
herramienta “Solver” de la siguiente forma:
Luego, ejecutamos, y la herramienta “Solver” nos entrega la solución oss =1,00, la cual
podemos apreciar en la siguiente imagen:
28
Solver con restricción
Procedemos a colocar los datos en la planilla Excel de la misma forma en la cual hemos
trabajado hasta ahora y en la herramienta “Solver”, pero esta vez vamos a agregar restricciones
la cual consiste en que 0 ≤ oss ≤ 1
29
Conclusión: Las herramientas de Excel soy muy útiles ya que además de ser muy fácil de
usarlas son rápidas en la resolución de la función. Podemos encontrar el valor de la variable
buscada (en nuestro caso la distancia entre la mandíbula inferior abierta, oss), sin tener la
necesidad de despejar la ecuación. Además en el caso especifico del “Solver”, esta herramienta
en comparación con la “Búsqueda de objetivo” será más óptima ya que esta herramienta
posee la opción de restringir la función, la cual al hacer esto la solución de la ecuación será más
óptima en un rango determinado.
Discusión: Como estos dos métodos son herramientas de Excel con lleva a que Excel realice las
iteraciones en forma interna, la cual se demorará según el tipo de ecuación con la que uno
trabaje, si es muy compleja Excel se demorará un poco más en comparación con una ecuación
simple.
30