Servicios Personalizados
Revista
Articulo
Indicadores
- Citado por SciELO
- Accesos
Links relacionados
- Similares en SciELO
Compartir
Ingeniería, investigación y tecnología
versión On-line ISSN 2594-0732versión impresa ISSN 1405-7743
Ing. invest. y tecnol. vol.9 no.1 Ciudad de México ene./mar. 2008
Educación en Ingeniería
Descomposición cruzada separable en la solución del problema de asignacióndistribución
Separable cross decomposition to solve the assignrouting problem
M. ElizondoCortés y R. AcevesGarcía
Posgrado de Ingeniería. Sección de Investigación de Operaciones.
Universidad Nacional Autónoma de México
Emails: mareli@avantel.net, aceves@servidor.unam.mx
Recibido: agosto de 2006
Aceptado: marzo de 2007
Resumen
El Problema de Inventario Ruteo (Inventory Routing Problem) surge en un contexto logístico que se presenta en las empresas y que pretende satisfacer las demandas de un conjunto de clientes distribuidos geográficamente, utilizando una flotilla de vehículos de capacidad limitada que se encuentran en un almacén central, al menor costo posible. El IRP es un problema NPduro que en aplicaciones reales suele ser de gran tamaño. Para su resolución se diseñó una estrategia que utiliza de forma conjunta, la descomposición cruzada y la relajación Lagrangena separable en la solución de la fase de asignacióndistribución, con lo que se obtienen un esquema tipo pingpong entre los dos subproblemas, que son del tipo transporte, para el cual se tiene un algoritmo de solución muy eficiente de orden O(n3) fácil de implementar para el problema completo.
Descriptores: Cadena de suministro, distribución, ruteo, descomposición cruzada separable.
Abstract
The InventoryRouting Problem emerges on a logistical context, that is presented into the companies and that it seeks to satisfy the demands of a group of clients distributed geographically, using a flotilla of vehicles of limited capacity, which are in a central ware house, at the small est possible cost. The IRP is a NPhard problem that is usually great size in real applications. For its solution was designed an strategy that uses of combined form, the crossed de composition and the separable Lagrangean relaxation in order to solve the assigndistribution phase, with what it is obtained a pingpong type scheme between two subproblems, which are from transport type, with which it is obtained a very efficient algorithm of order O(n3) and easy to implement for the complete problem.
Keywords: Supply chain, distribution, routing, separable cross decomposition.
Introducción
Uno de los obstáculos principales para la competitividad de las pequeñas y medianas empresas en México, es la falta de recursos y cultura para invertir en proyectos de desarrollo logístico. El desarrollo de estrategias poco costosas y fáciles de implementar que mejoren sus indicadores de desempeño, les proporcionarán instrumentos para sortear el inestable y cambiante entorno económico mexicano.
La investigación en el área logística de las últimas dos décadas, se ha enfocado a desarrollar varias estrategias encaminadas a coordinar la toma de decisiones en las actividades dentro de los elementos de las cadenas de suministro para mejorar su desempeño y efectividad en términos de costo, tiempos de respuesta, suministro a tiempo y servicio al cliente. Problemas como el de inventario manejado por el vendedor o problema de inventario ruteo (Elizondo, 2005) y el plan de reabastecimiento continuo, han sido algunas de las estrategias más promisorias para la coordinación de la cadena de suministro (Xu et al., 2001).
En particular, el problema de inventario ruteo (Inventory Routing Problem, IRP), modela una situación que se presenta comúnmente en las empresas e involucra en un solo modelo a las dos actividades más costosas de la cadena de suministro: el manejo de inventarios y la distribución física de productos. El IRP típico, considera que una compañía de distribución opera desde un almacén central y abastece a un gran número de clientes geográficamente distribuidos (Baita et al. 1998). En la aplicación práctica, el IRP resulta ser de gran tamaño, y adicionalmente se considera como un problema NPduro para un contexto de NPcompletez, debido a que el componente de ruteo de vehículos, además incluye restricciones de inventario. De tal forma que, la obtención óptima de la solución en tiempos razonables se torna prácticamente imposible.
El IRP tratado en Elizondo (2005) trabaja bajo dominio de tiempo (Baita et al. 1998), en intervalos discretos, en los cuales se deciden las cantidades a enviar y las rutas a seguir. En dicho trabajo, se analizaron estrategias de solución para problemas similares, tales como los de Federgruen y Zipkin (1984), Chien et al. (1989), Christiansen (1999) y Campbell et al. (2002), en los cuales se observó que si bien la mayoría ofrecen resultados satisfactorios para sus IRP específicos, un grave obstáculo común, fue el tratar con un problema entero mixto grande (problema de asignacióndistribución) para el cual no se han presentado alternativas de solución adecuadas. El presente artículo expone una estrategia de solución, a partir del modelo entero mixto denominado asignacióndistribución (AD), para el problema de decidir la asignación de clientes a vehículos y la cantidad de producto que le será entregada a cada uno de ellos. Dicha estrategia resultó ser muy eficiente y mejora los resultados obtenidos hasta ahora. Para su resolución, se utiliza la técnica de descomposición cruzada separable propuesta por Aceves (1996).
El problema de AsignaciónDistribución (AD)
Sea el problema AD, en el cual un número finito m de vehículos disponibles distribuyen un solo tipo de producto, para atender la demanda de una cierta población de usuarios concentrada en n puntos discretos, cada uno con demanda dj , cuando un vehículo en particular es seleccionado, se incurre en un costo fijo fi por la utilización del vehículo i y un costo variable cij xij que está en función del costo unitario de viaje del vehículo i al destino j atendiendo la fracción xij de la demanda del cliente j. Asumiendo que la capacidad de cada vehículo es limitada, el problema consiste en decidir cuáles de los posibles vehículos serán utilizados, de tal manera que sus capacidades no sean excedidas y las demandas satisfechas; así como qué patrón de distribución deberá utilizar, tal que el costo total de establecer las unidades en servicio conformado por costos fijos más costos variables, sea minimizado en un horizonte finito de planeación. El problema AD se puede formular como sigue:
donde
m: número de vehículos disponibles
n: número de clientes
dj: demanda del cliente j
fi: costo fijo para el vehículo i
ai: capacidad del vehículo i
cij: costo (en función de la distancia recorrida) de distribución al cliente j utilizando el vehículo i
xij: fracción de la demanda total atendida del cliente j utilizando el vehículo i.
yi = 1 si se usa el vehículo i,. 0 en otro caso.
La restricción (1) asegura la atención total de la demanda, (2) establece la distribución sólo con vehículos activos, (3) considera el uso de suficientes vehículos para atender la demanda y (4) considera no exceder la capacidad del vehículo.
El AD tienen dos decisiones inherentes: elegir los vehículos que se utilizarán y la forma de distribuir mejor la demanda para atender a los clientes. Esta complejidad lo hace un atractivo campo para el uso de técnicas de descomposición, ya que si la decisión discreta de elegir el vehículo se ha tomado, el problema continuo de distribución generalmente es más fácil de resolver. Además tiene una estructura especial que se puede explotar por las técnicas de descomposición.
La investigación desarrollada por Aceves (1996), incorpora en un mismo proceso la descomposición cruzada (Van Roy, 1983) y la relajación Lagrangeana separable (Guinard and Kim, 1983), logrando la ventaja de que ninguna de las restricciones originales desaparece, por lo que no es necesario elegir entre la calidad de la cota que se obtiene y el grado de dificultad del problema que queda. También se establece, que no es necesario utilizar el problema maestro en la solución, i.e., es posible resolver al problema completo iterando únicamente entre subproblemas, evitando por completo a los problemas maestros. A este procedimiento se le denominó descomposición cruzada separable y es el que se emplea para resolver el problema AD.
Aplicando descomposición de Benders al problema primal (AD), se obtiene:
Sujeto a:
que resulta en un problema del tipo transporte; y aplicando el esquema de relajación Lagrangena separable propuesto por Aceves (1996), el problema AD se puede establecer a través del siguiente proceso.
Copiar
duplicar la restricción
dualizar la restricción de igualdad y separar al problema de tal forma que se obtengan dos subproblemas con los siguientes conjuntos de restricciones:
A= {(1), (2), (3)}
y B {(2), (3), (4)}.
Del conjunto de restricciones (A) se obtiene un problema en variables yi que al aplicarle el principio de linearización entera y realizar la separación para cada vehículo i,
{yi = l, i=l,...,m}
con valores conocidos de λ, i = 1,...,m, se obtiene la solución óptima
vi= {λi ai con λi. < 0 y 0 en otro caso}.
Y minimizar la contribución total generada por cada vi, se obtiene un problema del tipo mochila 0 1 de la forma:
Si se considera ahora el subproblema originado por el conjunto de restricciones (B), y sólo se usa la parte correspondiente a la variable yi, resulta en un problema con el mismo conjunto de restricciones que el mochila anterior, y con valor de
µi = (fi vi) / 2
la función objetivo es la misma. Por lo cual, al realizar operaciones y reagrupar términos entre ambos problemas, nos queda
que resulta ser también un problema del tipo mochila. En consecuencia, el problema dual
Lagrangeano (PDL) se puede formular como los subproblemas:
y
Cuya solución estará dada por
v (PDL)= v(Q) + v(S).
Como Q y S son función de λi, es posible establecer para λi= fi / ai, o de forma más general
que el valor de la función objetivo de Q es cero, con vi =λi ai. Por lo cual, cualquier valor de y cumple para obtener el mínimo del problema.
Ahora, si se considera una solución factible que cumpla con la restricción (4), se tiene que el producto
para yi =1 y λi, < 0 i = 1,...,m
y de esta forma aún se cumple la propiedad de dualidad débil de programación lineal. Por lo que es posible reforzar al problema S con la restricción (4), con lo que se obtiene un problema del tipo transporte. Suceso muy interesante, ya que existen algoritmos de solución muy eficientes.
Sujeto a:
Con todos estos resultados (descomposición de Benders y Lagrangeana separable), fue posible desarrollar un algoritmo de solución más sencillo y rápido, que los obtenidos hasta el momento (Aceves, 1996).
Algoritmo de descomposición cruzada separable para el problema de AsignaciónDistribución
1. Iniciar
2. Resolver. Subproblema dual SDλK (problema de transporte), para obtener
2.1 Calcular
3. Probar
entonces terminar. De otra forma, identificar cuáles yik = 1 para i = 1,...,m.
4. Resolver. Subproblema primal
(problema de transporte), para obtener v(SPyk).
5. Probar. Si v{SPyk )= v{SDλK),
entonces terminar. De otra forma, regresar a la etapa 2 pero ahora con λik).
En iteraciones sucesivas de la descomposición cruzada separable, las rutas más económicas se van haciendo cada vez más baratas, y en consecuencia, se les asigna cada vez más flujo hasta que se saturan, ya sea agotando la capacidad del vehículo o satisfaciendo la demanda, es decir:
De esta forma, es posible usar los vehículos que satisfacen la demanda, obtenidos de la solución del subproblema dual Lagrangeano SDλ como vehículos activos; es decir, los yi =1, con i = 1,.., m, que se fijará n en el subproblema primal SPy. Así, el algoritmo de descomposición cruzada separable termina en un número finito de iteraciones y obtiene la solución.
Eficiencia del algoritmo
Como se mencionó, el algoritmo produce una solución sucesiva de dos subproblemas del tipo de transporte. Es importante resaltar la enorme ventaja de tener este tipo de problemas, ya que sus algoritmos de solución son de bajo orden de complejidad, O(n3), siendo además éste, el orden que acota superiormente a cada una de las fases de la estrategia propuesta, situación que la hace extremadamente eficiente de manera global.
Para evaluar la eficiencia del algoritmo se generaron 7 problemas: 4 problemas pequeñ os (menos de 10 clientes), 2 medianos (30 y 35 clientes) y uno grande (150 clientes). Las capacidades, demandas y costos fijos, fueron tomados de las instancias y ejemplos mostrados en Aceves (1996). En ellos, las capacidades de los vehículos variaron entre 4, 8, 16 y 24m3 y los costos fijos variaron entre 20, 25, 30, 35, 40 y 50 unidades monetarias. Los datos generados como variables aleatorias fueron las demandas de los clientes, producidas con distribución uniforme (1,10) y las matrices de distancia entre clientes y entre almacén central y clientes, las cuales emularon un área de operaciones semejante a la del Distrito Federal también con distribución uniforme, las matrices de distancias generadas fueron densas, de hecho, completamente llenas.
Cada problema se solucionó para obtener el óptimo utilizando Bifurcación y Acotación. Después, cada problema se resolvió con el algoritmo propuesto, por medio de un programa de cómputo desarrollado en Pascal.
Los resultados se muestran en la tabla 1. Entre los resultados obtenidos se encuentran:
1) Tiempo de solución extremadamente corto, ya que en la solución de los problemas se realizaron como máximo tres iteraciones y en todos los casos se obtuvieron los resultados en menos de 2 segundos
2) La calidad de las soluciones resultó excelente, ya que sólo un problema presentó brecha de optimación de 0.1%, en todos los demás problemas evaluados se obtuvo el óptimo.
Conclusiones
Se obtuvo una estrategia O(n3) con excelente calidad de resultados, bajo esfuerzo computacional, fácil de aplicarse e implementarse, la cual supera las técnicas conocidas hasta el momento para solucionar un problema logístico grande que es esencialmente NPduro.
Las técnicas de optimización deben utilizarse al afrontar los retos económicos, sociales, políticos y tecnológicos a los que se enfrentan los tomadores de decisiones logísticos al tratar de hacer más competitivas a sus organizaciones, de tal forma, que al dirigir los esfuerzos necesarios al auge de las pequeñas y medianas empresas mexicanas, se contribuye al desarrollo de la economía del país.
Referencias
Aceves García R. Un algoritmo para resolver el problema de localización de servicios con restricciones de demanda y adicionales. Tesis (Doctorado en Ingeniería). México. Universidad Nacional Autónoma de México, Ciudad Universitaria. 1996. [ Links ]
Baita F., Ukovich W., Pesenti R. and Favaretto D. Dynamic routingandinventory problems: A Review. Transportation Research. 32(8):585598. 1998 [ Links ]
Campbell A., Clarke L. and Savelsbergh M. Inventory routing in practice. In: The vehicle routing problem (Eds. P. Toth and D. Vigo), pp. 309330. Philadelphia: Society for Industrial and Applied Mathematics, 2002. [ Links ]
Chien W., Balakrishnan A. and Wong R. An integrated inventory allocation and vehicle routing problem. Transportation Science, (23):67. 1989. [ Links ]
Christiansen M. Decomposition of a combined inventory and time constrained ship routing problem. Transportation Science, (33): 316., 1999. [ Links ]
ElizondoCortés M. Una estrategia para resolver el problema de inventario distribución. Tesis (Doctorado en ingeniería). México. Universidad Nacional Autónoma de México, Ciudad Universitaria. 2005. [ Links ]
Federgruen A. and Zipkin P. A comprised vehicle routing and inventory allocation problem. Operations Research, (32): 1019 1037. 1984. [ Links ]
Guignard M. and Kim S. A strong Lagrangean relaxation for capacitated plant location problem. Technical report #56. Department of Statistics. The Wharton School, University of Pennsylvania, 1983. [ Links ]
Van Roy T.J. Cross decomposition for mixed integer programming. Mathematical Programming, (25): 4663. 1983. [ Links ]
Xu K., Dong Y. and Evers P. Towards better coordination of the supply chain. Transportation Research Part E, (37): 3554. 2001. [ Links ]
Semblanza de los autores
Mayra ElizondoCortés. Obtuvo el grado de maestra en ingeniería en investigación de operaciones por la Universidad Nacional Autónoma de México, más tarde el de doctora en investigación de operaciones en la misma institución. Su línea de investigación incluye optimización y simulación aplicados a logística y cadena de suministro. Es profesora de tiempo completo e investigadora en el Posgrado de Ingeniería de la UNAM, en donde ha dirigido tesis de maestría y ha participado en comités doctorales, colaborando en diversos proyectos. Ha participado en varios congresos y seminarios, así también ha publicado artículos en revistas internacionales.
Ricardo AcevesGarcía. Es ingeniero químico egresado de la Universidad Autónoma de Puebla. Obtuvo el grado de maestro y doctor en investigación de operaciones por la Facultad de Ingeniería de la UNAM. Ha trabajado en varios proyectos dentro del campo de transporte e investigación de operaciones para organismos públicos y privados. Actualmente es profesor de tiempo completo e investigador en la Facultad y el Posgrado de Ingeniería de la UNAM.