Mestodo Simplex para Una Mina Subterranea
Mestodo Simplex para Una Mina Subterranea
Mestodo Simplex para Una Mina Subterranea
SIMPLEX A MINERA
SUBTERRNEA
GRUPO:
G3 A1
I. INTRODUCCIN
El mtodo simplex es un mtodo muy prctico, ya que solo trabaja con los coeficientes de la funcin
objetivo y de las restricciones.
El mtodo Simplex es un procedimiento iterativo que permite mejorar la solucin de la funcin objetivo
en cada paso. El proceso concluye cuando no es posible continuar mejorando dicho valor, es decir, se
ha alcanzado la solucin ptima (el mayor o menor valor posible, segn el caso, para el que se
satisfacen todas las restricciones).
Partiendo del valor de la funcin objetivo en un punto cualquiera, el procedimiento consiste en buscar
otro punto que mejore el valor anterior. Como en el mtodo Grfico, dichos puntos son los vrtices del
polgono o (poliedro, si el nmero de variables es mayor de 2) que constituye la regin determinada
por las restricciones a las que se encuentra sujeto el problema (llamada regin factible). La bsqueda
se realiza mediante desplazamientos por las aristas del polgono, desde el vrtice actual hasta uno
adyacente que mejore el valor de la funcin objetivo. Siempre que exista regin factible, como su
nmero de vrtices y de aristas es finito, ser posible encontrar la solucin.
El mtodo Simplex se basa en la siguiente propiedad: si la funcin objetivo Z no toma su valor mximo
en el vrtice A, entonces existe una arista que parte de A y a lo largo de la cual el valor de Z aumenta.
Ser necesario tener en cuenta que el mtodo Simplex nicamente trabaja con restricciones del
problema cuyas inecuaciones sean del tipo "" (menor o igual) y sus coeficientes independientes sean
mayores o iguales a 0. Por tanto habr que estandarizar las restricciones para que cumplan estos
requisitos antes de iniciar el algoritmo del Simplex. En caso de que despus de ste proceso
aparezcan restricciones del tipo "" (mayor o igual) o "=" (igualdad), o no se puedan cambiar, ser
necesario emplear otros mtodos de resolucin, siendo el ms comn el mtodo de las Dos Fases.
El objetivo consistir en maximizar o minimizar el valor de la funcin objetivo
incrementar ganancias o reducir prdidas.
que genere,
El mtodo simplex es un procedimiento iterativo que permite tender progresivamente hacia la solucin
ptima. Es un procedimiento sistemtico y eficiente para encontrar y probar soluciones situadas en los
vrtices de optimalidad.
El mtodo requiere que las restricciones sean ecuaciones en lugar de inecuaciones, lo cual se aade
variables de holgura a cada inecuacin del modelo, variables que nunca pueden ser negativas y
tienen coeficiente 0 en la funcin objetiva
Conceptos importantes:
Solucin bsica: Valores de las variables que satisfacen las restricciones de igualdad de un
programa lineal en forma estndar, despus de que las variables no bsicas se toman como
cero.
Solucin bsica factible inicial: Valores de las variables que satisfacen las restricciones de
igualdad y de no negatividad de un programa lineal en forma estndar, despus de que las
variables no bsicas se toman como cero.
Variable de holgura: variable no negativa que se aade al lado izquierdo de una restriccin
menor o igual que, para obtener una restriccin de igualdad equivalente.
Variable artificial: variable no negativa que se aade al lado izquierdo de una restriccin
mayor o igual que, para obtener una restriccin de igualdad equivalente.
Iteracin: una serie de pasos de un algoritmo que se repiten.
Prueba de optimidad: Mtodo para determinar si la solucin obtenida es la ptima.
Mejora: proceso de encontrar soluciones factibles con valores de la funcin objetivo cada vez
mejores.
Preparando el modelo para adaptarlo al mtodo Simplex
La forma estndar del modelo de problema consta de una funcin objetivo sujeta a determinadas
restricciones:
Funcin objetivo: c1x1 + c2x2 +... + cnxn
Sujeto a:
a11x1 + a12x2 +
a21x1 + a22x2 +
...
am1x1 + am2x2 +
x1,..., xn 0
...
...
+
+
a1nxn =
a2nxn =
b1
b2
...
amnxn =
bm
(por ejemplo,
2.
3.
4.
5.
Tipo de optimizacin.
Como se ha comentado, el objetivo del mtodo consistir en optimizar el valor de la funcin objetivo.
Sin embargo se presentan dos opciones: obtener el valor ptimo mayor (maximizar) u obtener el valor
ptimo menor (minimizar).
Adems existen diferencias en el algoritmo entre el objetivo de maximizacin y el de minimizacin en
cuanto al criterio de condicin de parada para finalizar las iteraciones y a las condiciones de entrada y
salida de la base. As:
Objetivo de maximizacin
Objetivo de minimizacin
Solucin: Realmente no existe este problema dado que para que la solucin sea superior a 0 es
necesario que alguna restriccin tenga impuesta la condicin "" (y se tratara de un modelo para el
mtodo de las Dos Fases). En el caso planteado, la solucin real debe ser cero.
Cambio de signo de los trminos independientes.
Tambin se ha dicho que los trminos independientes (bi) de cada ecuacin deben ser no negativos
para poder emplear el mtodo Simplex. A tal fin, si alguna de las restricciones presenta un trmino
independiente menor que 0 habr que multiplicar por "-1" ambos lados de la inecuacin (teniendo en
cuenta que esta operacin tambin afecta al tipo de restriccin).
Ventajas: Con sta simple modificacin de signos en las restricciones correspondientes se posibilita la
aplicacin del mtodo Simplex al problema modelado.
Inconvenientes: Puede resultar que en las restricciones donde tengamos que modificar los signos de
las constantes, los tipos de desigualdad fueran "" (quedando tras la operacin del tipo "") siendo
necesario desarrollar el mtodo de las Dos Fases. Este inconveniente no es controlable, aunque
podra ocurrir el caso contrario y resultar beneficioso si los trminos independientes negativos se
presentan en todas aquellas restricciones con desigualdad de tipo "". Si existe alguna restriccin del
tipo "=" no supondra ninguna ventaja ni desventaja puesto que siempre sera de necesaria aplicacin
el mtodo de las Dos Fases.
Para normalizar una restriccin con una desigualdad del tipo "", hay que aadir una nueva variable,
llamada variable de holgura xs (con la condicin de no negatividad: xs 0). Esta nueva variable aparece
con coeficiente cero en la funcin objetivo, y sumando en la ecuacin correspondiente (que ahora s ser
una identidad matemtica o ecuacin de igualdad).
a11x1 + a12x2 b1
En caso de una desigualdad del tipo "", tambin hay que aadir una nueva variable llamada variable de
exceso xs (con la condicin de no negatividad: xs Esta nueva variable aparece con coeficiente cero en la
funcin objetivo, y restando en la ecuacin correspondiente.). 0
Surge ahora un problema con la condicin de no negatividad con esta nueva variable del problema. Las
inecuaciones que contengan una desigualdad de tipo "" quedaran:
------------------------------------------------------------------------------------------------------------------------Facultad de Ingeniera Geolgica, Minera y Metalrgica
a11x1 + a12x2 b1
a11x1 + a12x2 - 1xs = b1
Al realizar la primera iteracin con el mtodo Simplex, las variables bsicas no estarn en la base y
tomarn valor cero. En este caso la nueva variable xs, tras hacer cero a x1 y x2, tomar el valor -b1 y no
cumplira la condicin de no negatividad. Es necesario aadir otra nueva variable xr, llamada variable
artificial, que tambin aparecer con coeficiente cero en la funcin objetivo y sumando en la restriccin
correspondiente. Quedando entonces de la siguiente manera:
a11x1 + a12x2 b1
Al contrario de lo que cabra pensar, para las restricciones de tipo "=" (aunque ya son identidades) tambin
es necesario agregar variables artificiales xr. Como en el caso anterior, su coeficiente ser cero en la
funcin objetivo y aparecer sumando en la restriccin correspondiente.
a11x1 + a12x2 = b1
En el ltimo caso se hace patente que las variables artificiales suponen una violacin de las leyes del
lgebra, por lo que ser necesario asegurar que dichas variables artificiales tengan un valor 0 en la
solucin final. De esto se encarga el mtodo de las Dos Fases y por ello siempre que aparezcan este tipo
de variables habr que realizarlo.
En la siguiente tabla se resume segn la desigualdad el tipo de variable que aparece en la ecuacin
normalizada, as como su signo
Tabla
C1
C2
...
Cn
Base
Cb
P0
P1
P2
...
Pn
P1
Cb1
b1
a11
a12
...
a1n
P2
Cb2
b2
a21
a22
...
a2n
...
...
...
...
...
...
...
Pm
Cbm
bm
am1
am2
...
amn
Z0
Z1-C1
Z2-C2
...
Zn-Cn
Todos los valores incluidos en la tabla vendrn dados por el modelo del problema salvo los
valores de la fila Z (o fila indicadora). Estos se obtienen de la siguiente forma: Z j = (CbiPj) para i =
1..m, donde si j = 0, P0 = bi y C0 = 0, y en caso contrario Pj = aij.
Se observa, al realizar el mtodo Simplex, que en esta primera tabla ocupan la base todas las
variables de holgura y por ello (todos los coeficientes de las variables de holgura son 0 en la funcin
objetivo) el valor inicial de Z es cero.
Por este mismo motivo tampoco es necesario realizar los clculos de los costes reducidos en la
primera tabla, pudindose determinar directamente como el cambio de signo de los coeficientes de
cada variable en la funcin objetivo, esto es, -C j.
Condicin de parada:
Se cumple la condicin de parada cuando la fila indicadora no contiene ningn valor negativo entre
los costes reducidos (cuando el objetivo es la maximizacin), esto es, no existe posibilidad de
mejora.
Si no se cumple la condicin de parada es necesario realizar una iteracin ms del algoritmo, esto
es, determinar la variable que se vuelve bsica y la que deja de serlo, encontrar el elemento pivote,
actualizar los valores de la tabla y comprobar si se cumple nuevamente la condicin de parada.
Es tambin posible determinar que el problema no se encuentra acotado y su solucin siempre
resultar mejorable. En tal caso no es necesario continuar iterando indefinidamente y se puede
finalizar el algoritmo. Esta situacin ocurre cuando en la columna de la variable entrante a la base
todos los valores son negativos o nulos.
Elemento pivote:
El elemento pivote de la tabla queda marcado por la interseccin entre la columna de la variable
entrante y la fila de la variable saliente.
Actualizacin de la tabla:
Las filas correspondientes a la funcin objetivo y a los ttulos permanecern inalteradas en la nueva
tabla. El resto de valores debern calcularse como se explica a continuacin:
En la fila del elemento pivote cada nuevo elemento se calcula como:
Nuevo Elemento Fila Pivote = Anterior Elemento Fila Pivote / Pivote.
En el resto de las filas cada elemento se calcula:
Nuevo Elemento Fila = Anterior Elemento Fila - (Anterior Elemento Fila en Columna Pivote *
Nuevo Elemento Fila Pivote).
De esta forma se consigue que todos los elementos de la columna de la variable entrante sean
nulos salvo el de la fila de la variable saliente cuyo valor ser 1. (Es anlogo a utilizar el mtodo
de Gauss-Jordan para resolver sistemas de ecuaciones lineales).
Aplicaciones del algoritmo simplex
Para poder solucionar un problema mediante un algoritmo primero se debe extraer toda la
informacin que aporta el enunciado y preparar el problema de acuerdo a las necesidades del
mtodo resolutivo.
Los pasos para modelar un problema son los siguientes:
Paso 1: Determinar las variables de decisin y expresarlas algebraicamente.
o X1,..., Xn
Paso 2: Determinar las restricciones y se expresarlas como ecuaciones o inecuaciones
dependientes de las variables de decisin:
O A11X1 + A12X2 +... + A1nXn , , = b1
O A21X1 + A22X2 +... + A 2nXn , , = b2
o Am1X1 + Am2X2 + ... + AmnXn , , = bm
Paso 3: Expresar todas las condiciones implcitamente establecidas por la naturaleza de las
variables: que no puedan ser negativas, que sean enteras, que solo puedan tomar
determinados valores, ...
o X1,..., Xn 0
VENTAJAS
1. Se visualiza bien, como es el funcionamiento de los algoritmos.
2. Es una herramienta muy til que le permite al estudiante comprobar los resultados obtenidos en
un problema de P.L. y corregir errores que se puedan cometer.
3. Es de fcil accesibilidad ya que se encuentra de manera gratuita en la web.
En una empresa minera del sur del pas que opera por el mtodo de explotacin subterrnea,
actualmente desea implementar su extraccin y se ha propuesto iniciar con la explotacin de dos
unidades mineras, con las cuales se han propuesto extraer dos tipos de minerales, se sabe tambin
que la capacidad de la chancadora en la unidad A es de 45 tm/hr y en la unidad B de 80 tm/dia, y la
capacidad del molino para la unidad A es de 70tm/hr y en la unidad B de 50 tm/hr. En la unidad B la
mxima produccin de dicho mineral es de 800 tm/dia. El gerente de operaciones pide maximizar la
venta.
CAPACIDAD
DE CAPACIDAD
DE VALOR
CHANCADORA
MOLIENDA
MINERAL
GALERIA
(TM/HR)
45
(TM/HR)
70
(US$/TM)
25
A
GALERIA
80
50
35
22HR/DIA
22HR/DIA
DEL
IV. SOLUCIN
X2<=800.(3)
TABLA INICIAL
BASE
Variables
DECISIN
de VARIABLES
HOLGURA
DE SOLUCI
N
OPERACI
N
X1
X2
S1
S2
S3
S1
7700
1100
S2
16
15840
1760
S3
800
800
-25
-35
ITERACIN No 1
BASE
Variables
DECISIN
de VARIABLES
HOLGURA
DE SOLUCI
N
OPERACI
N
X1
X2
S1
S2
S3
S1
7700
1100
S2
16
15840
1760
S3
800
800
-25
-35
RESULTADO DE ITERACIN No
BASE
Variables
DECISIN
de VARIABLES
HOLGURA
DE SOLUCI
N
X1
X2
S1
S2
S3
X1
-7
2100
S2
16
-9
8640
S3
800
-25
35
28000
ITERACIN No
BASE
OPERACI
N
(6)
Variables
DECISIN
de VARIABLES
HOLGURA
DE SOLUCI
N
OPERACI
N
X1
X2
S1
S2
S3
X1
-7
2100
420
S2
16
-9
8640
540
S3
800
-25
35
28000
RESULTADO DE ITERACIN No
BASE
X1
Variables
DECISIN
de VARIABLES
HOLGURA
DE SOLUCI OPERACI
N
N
X1
X2
S1
S2
S3
X1
1/5
-7/5
420
S2
-16/5
67/5
1920
X2
800
385
00
TABLA FINAL
BASE
Variables
DECISIN
de VARIABLES
HOLGURA
DE SOLUCI OPERACI
N
N
X1
X2
S1
S2
S3
X1
1/5
-7/5
420
S2
-16/5
67/5
1920
X2
800
385
00
Como todos los coeficientes de la fila de la funcin objetivo son positivos, hemos llegado a
la solucin ptima. La solucin ptima viene dada por el valor de Z en la columna de los
valores solucin, en nuestro caso: 38500.
DONDE
X1 = 420
DONDE
X2 = 800
PROBLEMAS PROPUESTOS
Problema N2:
Impala Gold Company operaba una mina de oro en el Estado Libre de Orange, Sudfrica. La
operacin consista en minera subterrnea, a una profundidad de 1400 metros, para la extraccin de
mineral de oro. Las rocas se transportaban por los piques de la mina a un molino que trituraba la roca
y extraa el oro.
La mina Impala tena tres piques. La informacin acerca de estos piques se presenta en la tabla
adjunta. Observe que el mineral que se obtiene de cada una de las reas de los Piques tiene diferente
contenido de oro y distinto costo.
Las rocas que extraan de los tres piques subterrneos se enviaban al molino para ser trituradas y
refinadas. La capacidad del molino dependa de la finura del molido de las rocas. Si las rocas se
molan finamente, la capacidad del molino era de 240000 toneladas mensuales y se recuperaba el
95% del oro en la operacin. Las rocas de cada tnel se podan moler por separado. El costo de
molido fino de una tonelada de roca era $1.12 por tonelada. Si el molido de las rocas era grueso, la
capacidad del molino era de 250000 toneladas, pero la recuperacin de oro bajaba al 90%. El costo
de molido grueso de una tonelada de roca era de $0.85. La mina poda vender todo el oro que
produjera, a $0.80 el gramo.
El gerente de la mina estaba preocupado por la cantidad de mineral que deba extraer de cada una de
las reas de los Piques. Haba observado que la capacidad del molino no era suficiente para manejar
todos los Piques si stos operaban a toda su capacidad. El problema se complicaba an ms por el
requisito legal de que la mina no poda operar por encima del grado promedio de las reservas de
mineral. En la mina Impala este grado promedio era de 20 gramos por tonelada. Por lo tanto, exista la
restriccin legal de que la mezcla de rocas de los tres Piques no poda exceder un promedio de 20
gramos de oro por tonelada de mineral.
Formule un modelo de programacin lineal para maximizar los beneficios de la operacin de la mina.
Capacidad
de
transporte
del
Pique(tn/mes)
Grado
del
mineral(grAu/tn)
Costo
de
extraccin($/tn)
Pique N 1
85000
Pique N2
90000
Pique N3
95000
25
20
15
Solucin
Solucin del problema usando el tablero simplex:
X1F +X1G<=85000..(1)
X2F +X2G<=90000(2)
X3F +X3G<=95000(3)
Ecuacin de la recta:
F
G
+
=1
24000 25000
Inecuacin de la recta:
F
G
+
1
24000 25000
250000
F+G 25000
240000
25(X1F +X1G) +20(X2F +X2G) +15(X3F +X3G) <=20(X1F +X2F +X3F +X1G +X2G +X3G)
5 (X1F) -5 (X3F) +5 (X1G) -5(X3G) <=0
Funcin objetivo:
Maximizar:
Sujeto a:
X1F +X1G<=85000
------------------------------------------------------------------------------------------------------------------------Facultad de Ingeniera Geolgica, Minera y Metalrgica
X2F +X2G<=90000
X3F +X3G<=95000
1.04(X1F) + 1.04(X2F) + 1.04(X3F) + (X1G) + (X2G) + (X3G) <= 250000
5 (X1F) - 5(X3F) + 5(X1G) -5(X3G) <=0
Transformamos todas las inecuaciones a las ecuaciones:
Z - 11.88(X1F) - 9.08(X2F) - 6.28(X3F) - 11.15(X1G) - 8.55(X2G) - 5.95(X3G) = 0
X1F +X1G +S1 = 85000
X2F +X2G +S2 = 90000
X3F +X3G +S3 =95000
1.04(X1F) + 1.04(X2F) + 1.04(X3F) + (X1G) + (X2G) + (X3G) + S4 = 250000
5 (X1F) - 5(X3F) + 5(X1G) -5(X3G) + S5 = 0
Llenamos la tabla:
Primero la columna pivote en la fila objetivo el cual ser el elemento mas negativo de dicha fila
Luego dividimos la columna del extremo derecho con la columna pivote para ubicar cual ser nuestra fila pivote
el cual ser escogiendo al menor nmero positivo de dicha operacin:
Una vez ubicados la fila y columna pivote, la interseccin ser nuestro elemento pivote el cual tendremos que
convertirlo en la unidad, para hacer dicha operacin debemos dividir el nmero que convierte a la unidad, a
toda la fila.
Luego procedemos a iterar con nuestra fila pivote las dems filas convirtiendo a cero a todos los elementos que
estn arriba o abajo del elemento pivote, obtenindose la siguiente tabla:
Ubicamos nuestra fila y columna pivote y asi obtener nuestro elemento pivote
Ubicamos nuestra fila y columna pivote y asi obtener nuestro elemento pivote
Ubicamos nuestra fila y columna pivote y asi obtener nuestro elemento pivote
------------------------------------------------------------------------------------------------------------------------Facultad de Ingeniera Geolgica, Minera y Metalrgica
Iteramos la tabla hasta logra obtener que todos los coeficientes de la funcin objetivo sea positivo, obteniendo
la siguiente tabla:
Problema N3:
La empresa Minas Poracota S.A.Ubicado en la regin Arequipa, provincia de Condesuyo distrito de
Cayarani tiene asentada su unidad minera Poracota la cual se dedica a la explotacin y transporte de
mineral.Sierto da se le asigna al ingeniero de turno que se encargue de transportar los materiales
extrados de dos labores mineras una de chimenea y la otra de cruzero hacia tres plantas
concentradoras, sabiendo que la primera labor minera tiene una produccin de 10Tn y la segunda de
15Tn de material por turno ,adems las dos primeras plantas concentradoras pueden procesar 8Tn de
material ,mientras que la tercera necesita 9 toneladas por turno. El costo de transporte desde cada
labor minera a la planta concentradora viene dado por el siguiente cuadro:
Labor minera
Labor
minera(1)CH
Labor
minera(2)CZ
Planta
(1)
10
15
concentradora Planta
(2)
15
10
concentradora Planta
(3)
20
concentradora
10
Restricciones:
X1+x2+x3<=10 (1)
X2+x4+x6<=15.. (2)
X1+x2>=8. (3)
X3+x4>=8. (4)
X5+x6>=9. (5)
Planta
Utilizando el programa de
para la solucin del problema.
investigacin
de
operaciones
TORA
1 Ingresar al programa.
4 Se abrira la siguiente ventana en donde se debe poner titulo al problema ,luego se ingresara el
numero de variables y resticiones a continuacion se aciona Enter.
5Al abrirse la siguiente ventana se elegir la accin que se desee realizar sea esta
maximizar o minimizar, luego se llenara la tabla con las restricciones del problema.
6Luego de haber llenado la tabla por completo se hace click en el botn SOLVE Menu.
7Dar
click en SI, luego seguir el siguiente
orden: Solve Problem, Algebraic, Iterations, Bounded Simplex.
De la tabla siquiente se llega a la conclusion que si se quiere aminorar los costos de traslado del
material desde la laborminera a la planta consentradora se deberian repartir desde la labor de cruzero
------------------------------------------------------------------------------------------------------------------------Facultad de Ingeniera Geolgica, Minera y Metalrgica
8Tn a la concentradora (1) ,2Tn a la concentradora (2) y 0 Tn a la concentradora (3);De igual forma
para la labor minera de chimenea se deberian repartir 0Tn a la concentradora (1),6Tn ala
concentradora(2)y9Tn a la concentradora (3).
PROGRAMA LINDO
LINDO es una aplicacin para computadoras que se utiliza para resolver problemas de programacin
lineal, cuadrtica y entera.
Desde 1979 el programa LINDO ha sido una de las herramientas de optimizacin favoritas de las
comunidades Educativas y Empresariales. LINDO Systems se ha dedicado a proveer poderosas e
innovativas herramientas de optimizacin que tambin son flexibles y muy fciles de usar. LINDO tiene
una larga historia y es uno de los pioneros en crear poderosos programas de optimizacin.
En 1979 se vendi en Mxico la primer copia comercial de LINDO, la ayuda que este proporcion en
aquel momento, le hizo ganar popularidad muy rpidamente para luego ser utilizado en aplicaciones
industriales. En 1983 la versin LINDO/PC fue el primer paquete para programacin lineal, este
manejaba 60 restricciones y 120 variables. En 1996 apareci la versin 6.0 para WINDOWS.
Antes que aparecieran Lotus 1-2-3 o Excel, LINDO haba sido incorporado a la planilla de clculo
VisiCalc, el paquete se llam VINO y es el equivalente del SOLVER que viene con Excel. En estos
momentos existe una hoja de clculo llamada WHATS Best! la cual se integra a EXCEL o LOTUS 12-3; sta resuelve problemas de optimizacin con algoritmos de LINDO Systems.
Esta imagen es una pantalla obtenida del programa LINDO:
La
le
VII.
RECOMENDACIONES
------------------------------------------------------------------------------------------------------------------------Facultad de Ingeniera Geolgica, Minera y Metalrgica
Para una mejor compresin del mtodo del algoritmo del tablero simplex se recomienda tener
conocimientos previos de algebra vectorial o algebra lineal para poder comprender de mejor
manera algunas definiciones que estn involucradas en el algoritmo como definicin de rango
de una matriz, operaciones elementales fila columna, sistema de ecuaciones lineales entre
otros.
En lo posible se debera utilizar ms de un software en la solucin delos problemas de
programacin lineal porque de esa forma se puede comparar resultados y darse cuenta que
programa es ms idneo para cada tipo particular de problema.
VIII. CONCLUSIONES
IX.BIBLIOGRAFA