Plan Maestro de Producción Basado en Programación Lineal Entera para Una Empresa de Productos Químicos
Plan Maestro de Producción Basado en Programación Lineal Entera para Una Empresa de Productos Químicos
Plan Maestro de Producción Basado en Programación Lineal Entera para Una Empresa de Productos Químicos
Eduardo Gutiérrez
Article
Plan maestro de producción basado en
programación lineal entera para una empresa de
productos químicos
Suggested Citation: Zotelo, Yunuem Reyes; Mula, Josefa; Díaz-Madroñero, Manuel; González,
Eduardo Gutiérrez (2017) : Plan maestro de producción basado en programación lineal entera
para una empresa de productos químicos, Revista de Métodos Cuantitativos para la Economía
y la Empresa, ISSN 1886-516X, Universidad Pablo de Olavide, Sevilla, Vol. 24, pp. 147-168
Die Dokumente auf EconStor dürfen zu eigenen wissenschaftlichen Documents in EconStor may be saved and copied for your
Zwecken und zum Privatgebrauch gespeichert und kopiert werden. personal and scholarly purposes.
Sie dürfen die Dokumente nicht für öffentliche oder kommerzielle You are not to copy documents for public or commercial
Zwecke vervielfältigen, öffentlich ausstellen, öffentlich zugänglich purposes, to exhibit the documents publicly, to make them
machen, vertreiben oder anderweitig nutzen. publicly available on the internet, or to distribute or otherwise
use the documents in public.
Sofern die Verfasser die Dokumente unter Open-Content-Lizenzen
(insbesondere CC-Lizenzen) zur Verfügung gestellt haben sollten, If the documents have been made available under an Open
gelten abweichend von diesen Nutzungsbedingungen die in der dort Content Licence (especially Creative Commons Licences), you
genannten Lizenz gewährten Nutzungsrechte. may exercise further usage rights as specified in the indicated
licence.
https://creativecommons.org/licenses/by-sa/4.0/
REVISTA DE MÉTODOS CUANTITATIVOS PARA LA
ECONOMÍA Y LA EMPRESA (24). Páginas 147–168.
Diciembre de 2017. ISSN: 1886-516X. D.L: SE-2927-06.
www.upo.es/revistas/index.php/RevMetCuant/article/view/2885
147
Master Production Scheduling Based on Integer
Linear Programming for a Chemical Company
ABSTRACT
148
1. Introducción
La planificación de la producción es un proceso continuo y complejo que consiste en
determinar anticipadamente decisiones que permitan optimizar el uso de los recursos
productivos. Si se conoce a intervalos regulares de tiempo la cantidad de producción
que será demandada, es posible desarrollar programas para la producción, distribución
del trabajo y compra de materiales, logrando que la compañía alcance el rendimiento
pleno, sin que sea necesario el almacenamiento de grandes stocks por un tiempo
excesivo, economizando el capital e interés de la empresa (Boiteux et al., 2007).
149
producción en la empresa objeto de estudio. En la Sección 4, se describe la metodología
de trabajo y la aplicación de los diferentes enfoques para la planificación revisados en la
Sección 2. En la Sección 5, se propone el modelo PMP utilizando programación lineal
entera. Seguidamente, en la Sección 6, se evalúan los resultados obtenidos por los
diferentes escenarios, derivados de la modificación de algunas características del
sistema, lo cual es posible debido a la flexibilidad del modelo. Finalmente, se presentan
las conclusiones y las líneas futuras de investigación en la Sección 7.
2. Revisión de la literatura
De manera jerárquica, la planificación de producción abarca desde las instalaciones de
producción, incluyendo estrategias de localización de plantas y sus capacidades,
métodos de pronósticos, hasta llegar al nivel de planta donde los temas incluyen
planificación táctica y operativa, planificación de materiales y gestión de inventarios.
150
considerando restricciones para la producción mínima de acuerdo a un tamaño de lote
para la industria dedicada a la fabricación de pinturas.
151
Actualmente, el sistema reactivo de producción en la empresa objeto de estudio
está basado en cubrir los límites de inventario preestablecidos de forma empírica como
mínimos y máximos; estos niveles establecen, además, el espacio en almacén para cada
producto. A través de un informe de existencias generado por el sistema ERP, se
observa que, cuando las reservas son iguales o menores al mínimo, se genera una orden
de producción por la cantidad necesaria para alcanzar el nivel máximo. Esto ha
generado un conflicto, principalmente, respecto a la optimización de los recursos y el
nivel de servicio, al realizar pedidos por lotes incompletos, rebasar por mucho el
espacio de almacén o desabastecerse de un momento a otro.
Esta situación también ha generado conflictos internos, entre los que destacan la
desconfianza de los trabajadores, dados los cambios súbitos en la planificación diaria
para cubrir el desabastecimiento de productos, situación que implica en algunos casos el
pago de tiempo extra; y, por otro lado, la incertidumbre a la que se enfrenta el área de
compras al no contar con información confiable en la cual basar sus negociaciones de
compra de materia prima con los proveedores.
4. Metodología de trabajo
Para el diseño del modelo propuesto, se tomaron como base diferentes trabajos descritos
en la literatura, específicamente aquellos relacionados con la programación lineal
entera, MRP, PMP y cadena de suministro; adecuándolo a un modelo mononivel con
demanda discreta (Mula et al., 2006a, 2006b; Mula et al., 2014; Venkataraman y
Nathan, 1994). En Mula et al. (2006a, 2006b), los autores proponen un modelo de
programación lineal a medio plazo con restricciones de capacidad en un entorno
multiproducto, multiperiodo y multinivel para determinar el PMP y, con ello, los niveles
de inventario, la demanda diferida y los niveles de capacidad utilizada en un horizonte
específico de planificación. Estos objetivos son compartidos con el modelo propuesto,
sin considerar la información referente a las materias primas, que es tomada en cuenta
para la formulación del MRP, pero que no aplican en este trabajo; por lo que
únicamente se consideran conceptos como la minimización de costes por tiempo ocioso,
tiempo extra y los relacionados con la producción, el almacenamiento y la demanda
diferida de los artículos. Del modelo de Mula et al. (2014) se rescata la sincronización
de la información de la demanda a través de las múltiples etapas de la cadena de
suministro, concretamente en cuanto a los límites de inventario de cada una de ellas. El
modelo considera niveles iniciales de inventario, stock de seguridad y tamaños de lotes
predeterminados, introduciendo además la política de demanda diferida nula. Por su
parte, Venkataraman y Nathan (1994), a través del desarrollo de un PMP multiobjetivo
utilizando un horizonte rodante en un entorno de demanda determinista y múltiples
líneas de producción para una industria química dedicada a la fabricación de pinturas,
aporta al modelo propuesto el concepto de “tamaño de lote mínimo”, holgura respecto a
los niveles de producción e inventario en comparación a un nivel deseado, exceso de
inventario, tiempo ocioso y tiempo extra.
152
planificación (Venkataraman y Nathan, 1994). La propuesta presentada está diseñada en
base a la demanda en el almacén de producto terminado, que funciona a su vez como
centro de distribución, y está limitada por los recursos y tiempos implicados para la
producción. Los resultados del modelo brindarán información valiosa para el
requerimiento de materiales y la programación detallada, con la posibilidad de
ampliarse para considerar tiempos de entrega al consumidor final en tienda (véase
Figura 1).
Aprovisionamiento
Producción: Información para el modelo Distribución
Almacén
Almacén
de
de Planta de Producto
Materia producción Terminad
Prima
o CEDIS
Para realizar el modelo, se requirió información del proceso, los productos y los
recursos implicados. Por tanto, el modelo se define así:
Dados:
− Datos de productos: demanda, tiempo de espera para la entrega, tamaño de
lote de pedido, tanque utilizado, demanda diferida inicial, inventario inicial,
inventario mínimo y máximo permitido en almacén, coste de producción,
inventario y demanda diferida por periodo.
− Datos de los recursos productivos: tiempo requerido del recurso para la
producción, capacidad disponible, coste del tiempo extra y ocioso de cada
recurso.
− Datos para la planificación: número de periodos temporales en el horizonte
de planificación.
Determinar:
− La cantidad a producir y el número de lotes de cada producto.
− Inventario del producto al final de cada periodo.
− Demanda diferida por producto.
− Desviación con respecto al mínimo y máximo establecido en almacén.
− Tiempo ocioso y extra en cada recurso por periodo.
Para cumplir los objetivos principales:
− Minimizar los costes de producción, inventario y demanda diferida.
− Minimizar el tiempo ocioso y extra de cada recurso productivo.
153
− Minimizar la desviación con respecto a los límites establecidos en inventario.
Asumiendo que:
− Solo se puede realizar la fabricación de lotes completos.
− El coste por demanda diferida es significativamente alto.
− Se elaborarán modificaciones al modelo según el nivel de servicio requerido
por la empresa para los productos seleccionados.
154
La formulación del modelo de programación lineal entera es la siguiente:
Función objetivo:
ூ ்
Restricciones:
a) Balance de inventario:
ܸܰܫ௧ = ܸܰܫ,௧ିଵ + ݇,௧ି்ௌ ݈ݐ − ݀௧ − ܴ݀,௧ିଵ + ܴ݀௧ ∀݅, ݐ (2)
b) Balance de capacidad:
ூ
c) No negatividad:
(7)
ܸܰܫ௧ , ݇௧ , ݈݊݉ܫ௧ , ݔ݉ܫℎ௧ ≥ 0
d) Variables enteras:
(8)
݇௧ , ܸܰܫ௧ , ܴ݀௧ , ݈݊݉ܫ௧ , ݔ݉ܫℎ௧ ∈ ℤ
155
ser utilizado. Las ecuaciones (7) y (8) establecen que, dadas las características del
problema, las variables de decisión deberán ser enteras y positivas.
Escenario 1: lotes fijos y niveles de inventario mínimo y máximo flexibles, sin demanda
diferida.
156
preestablecidos para el inventario, de acuerdo a la capacidad del almacén, pueden
desviarse del valor establecido, pero esta desviación será minimizada por la función
objetivo. En este caso, la formulación del modelo a utilizar es la misma que la del
modelo inicial, expresado en las ecuaciones (1) a (8).
Escenario 2: lotes fijos sin niveles de inventario mínimo y máximo, sin demanda
diferida.
Función objetivo:
ூ ்
Restricciones:
a) Balance de inventario:
Ecuaciones (2) y (5).
b) Balance de capacidad:
Ecuación (6).
157
c) No negatividad:
Escenario 3: niveles de inventario mínimo y máximo flexibles sin lotes fijos, sin
demanda diferida.
Función objetivo:
ூ ்
Restricciones:
a) Balance de inventario:
ܸܰܫ௧ = ܸܰܫ,௧ିଵ + ܲ,௧ି்ௌ − ݀௧ − ܴ݀,௧ିଵ + ܴ݀௧ ∀݅, ݐ (14)
Ecuaciones (3) a (5).
b) Balance de capacidad:
c) No negatividad:
(16)
ܸܰܫ௧ , ܲ௧ , ݈݊݉ܫ௧ , ݔ݉ܫℎ௧ ≥ 0
158
d) Variables enteras:
donde Pit corresponde a la cantidad a producir del producto i en el periodo t sin tener en
cuenta la multiplicidad de tamaño de lote.
Escenario 4: sin lotes fijos ni niveles de inventario mínimo y máximo, sin demanda
diferida
Función objetivo:
ூ ்
Restricciones:
a) Balance de inventario:
Ecuación (14).
b) Balance de capacidad:
Ecuación (15).
c) No negatividad:
Ecuación (16).
d) Variables enteras:
Ecuación (17).
159
Tabla 5. Casos destacados del escenario 4
Tiempos de entrega de producto
=0 >0
A4 ●
B4 ●
Escenario 5: lotes fijos y niveles de inventario mínimo y máximo flexibles con un nivel
de servicio preestablecido.
161
Tabla 7. Resultados de la optimización para los escenarios propuestos en una política de demanda
diferida nula
Uno de los resultados arrojados por el modelo es el tiempo ocioso y extra por
recurso (véanse Tablas 8 y 13), lo cual funciona como una herramienta para equilibrar
el uso de las instalaciones y evitar cuellos de botella. La Tabla 8 muestra, por ejemplo,
el tiempo ocioso y extra para el recurso 4 dedicado a la fabricación de tres productos en
el escenario A1. Los elevados tiempos extra en algunos periodos temporales sugieren
que debería considerarse la producción utilizando otros recursos, principalmente,
durante los períodos 2, 7, y 12.
162
Tabla 8. Tiempo ocioso y extra para el recurso 4 en el escenario A1
163
El escenario 5, por su parte, brinda a la empresa la opción de elegir la proporción
de demanda que está dispuesta a satisfacer de acuerdo al coste, mientras que la
desviación respecto a los niveles de inventario se optimizará en base al coste de las
unidades almacenadas. En la Tabla 10 se observan los resultados en costes para la
optimización del escenario 5 en los casos señalados en la Tabla 6. Para los escenarios
B5, D5 y F5, la única variable que cambia es la proporción satisfecha de demanda
teniendo como premisa un tiempo de suministro nulo. Es posible observar que la
diferencia en los costes entre un 95 y 90% de proporción satisfecha es mínima, aunque
un 8% mayor en cuanto a tiempo ocioso. En cambio, comparando el 95% con el 80%,
se observa una reducción de costes, no siendo así para tiempo ocioso y tiempo extra,
que aumentan en un 13,7% y un 0,22% en comparación con el primer caso. En estas
condiciones, es mucho mejor inclinarse por un 80% en la proporción de demanda
satisfecha que por un 90%. Este comportamiento se observa también en los casos donde
TS>0, como lo señala la Tabla 11.
Tabla 10. Resultados de la optimización para una política con demanda diferida
Demanda diferida de acuerdo a la proporción que la empresa está dispuesta a satisfacer
ESCENARIO 5
Lotes fijos y niveles de inventario mínimo y máximo flexibles
A5 (95%) B5 (95%) C5 (95%)
Valor objetivo z= 21.207.417,69 22.460.133,29 23.709.535,00
Coste total MXN$ 21.172.285,69 MXN$ 22.460.133,29 MXN$ 23.709.535,00
Coste de producción cpt= MXN$ 14.720.551,00 MXN$ 14.767.613,73 MXN$ 15.003.386,92
Coste de inventario cit= MXN$ 6.428.883,00 MXN$ 5.574.411,00 MXN$ 6.161.961,00
Coste por demanda diferida crdt= - - -
Coste por tiempo ocioso ctoc= MXN$ 1.688,45 MXN$ 1.429,72 MXN$ 1.316,56
Coste por tiempo extra ctex= MXN$ 21.163,24 MXN$ 20.653,84 MXN$ 20.695,52
Inventario < mínimo itnl= 4.027 MXN$ $0,00 MXN$ 0,00
Inventario > máximo itxh= 31.105 MXN$ $2.096.025,00 MXN$ 2.522.175,00
D5 (90%) E5 (90%) F5 (80%)
Valor objetivo z= 22.697.033,15 22.904.506,23 21.917.484,33
Coste total MXN$ 22.697.033,15 MXN$ 22.904.506,23 MXN$ 21.917.484,33
Coste de producción cpt= MXN$ 14.728.142,42 MXN$ 14.728.142,42 MXN$ 146.565.871,35
Coste de inventario cit= MXN$ 5.742.096,00 MXN$ 5.873.721,00 MXN$ 5.337.609,00
Coste por demanda diferida crdt= - - -
Coste por tiempo ocioso ctoc= MXN$ 1.544,01 MXN$ 1.330,55 MXN$ 1.625,66
Coste por tiempo extra ctex= MXN$ 20.790,71 MXN$ 20.577,26 MXN$ 20.698,31
Inventario < mínimo itnl= MXN$ 0,00 MXN$ 0,00 MXN$ 0,00
Inventario > máximo itxh= MXN$ $2.204.460,00 MXN$ $2.280.735,00 MX:$ 1.900.680,00
G5 (80%)
Valor objetivo z= 22.968.593,14
Coste total MXN$ 22.968.593,14
Coste de producción cpt= MXN$ 14.837.831,47
Coste de inventario cit= MXN$ 5.834.766,00
Coste por demanda diferida crdt= -
Coste por tiempo ocioso ctoc= MXN$ 1.197,32
Coste por tiempo extra ctex= MXN$ 20.543,35
Inventario < mínimo itnl= MXN$ 0,00
Inventario > máximo itxh= MXN$ 2.274.255,00
164
tiempos ociosos y tiempos extra para los recursos implicados (véase Tabla 13). A partir
del análisis de los resultados, la empresa podrá tomar la mejor decisión en base a sus
intereses.
Tabla 11. Desviación de las diferentes proporciones de demanda satisfecha respecto al 95%
TS=0 TS>0
Tabla 12. Variables a considerar para la planificación de producción según el escenario F5.
Producto Periodos Lotes Producción Inventario Demanda diferida < Mínimo > Máximo
7 1 1 40 122 54 0 70
7 2 0 0 48 39 0 0
7 3 1 40 49 61 0 0
7 4 1 40 32 43 0 0
7 5 1 40 40 60 0 0
7 6 3 120 49 0 0 0
7 7 1 40 45 12 0 0
7 8 1 40 33 1 0 0
7 9 2 80 52 4 0 0
7 10 2 80 85 0 0 33
7 11 1 40 70 0 0 18
7 12 1 40 69 0 0 17
12 1 0 0 0 81 0 41
12 2 0 0 11 28 0 0
12 3 1 40 12 32 0 0
12 4 2 80 65 54 0 14
12 5 0 0 73 28 0 0
12 6 1 40 27 37 0 0
12 7 1 40 47 36 0 0
12 8 0 0 24 28 0 0
12 9 1 40 0 29 0 0
12 10 1 40 0 32 0 0
12 11 1 40 0 28 0 0
12 12 1 40 1 30 0 0
34 1 0 0 115 67 0 37
34 2 0 0 48 55 0 0
34 3 0 0 48 103 0 0
34 4 0 0 48 137 0 0
34 5 0 0 48 192 0 0
34 6 0 0 48 238 0 0
34 7 2 480 257 0 0 179
34 8 0 0 213 0 0 135
34 9 0 0 158 0 0 80
34 10 0 0 114 0 0 36
34 11 0 0 61 0 0 0
34 12 1 240 260 0 0 182
165
En el caso A5, el resultado de la función objetivo indica desviación absoluta, en
cuanto a unidades producidas, sin incluir costes.
7. Conclusiones
166
probado con datos históricos de demanda, por lo que en trabajos futuros se propone
incluir estos métodos y, además, la posibilidad de incorporar a este modelo un MRP
agregando información respecto a las materias primas y los envases y embalajes.
Referencias
Boiteux, O.D.; Corominas, A. y Lusa, A. (2007): “Estado del arte sobre planificación
agregada de la producción”. Enginyeria d’Organització i Logística Industrial, 4(2),
1–39.
Dzielinski, B.P. y Gomory, R.E. (1965): “Optimal programming of lot sizes, inventory
and labor allocations”. Management Science, 11(9), 874–890.
Flores, B.E. y Whybark, D.C. (1986): “Multiple criteria ABC analysis”. International
Journal of Operations and Production Management, 6(3), 38–46.
Hax, A.C. y Meal, H.C. (1975): “Hierarchical integration of production planning and
scheduling”. En M.A. Geisler (ed.): TIMS Studies in Management Sciences, Volume
1: Logistics. Nueva York: North Holland/American Elsevier, pp. 53–69.
Lasdon, L.S. y Terjung, R.C. (1971): “An efficient algorithm for multi-item
scheduling”. Operations Research, 19, 946–969.
Mangiameli, P.M. (1979): The effects of managerial policies on aggregate plans, the
master production schedule, and departmental plans. Tesis doctoral, Ohio: Ohio
State University.
167
Moniz, S.; Barbosa-Póvoa, A.P. y de Sousa, J.P. (2014): “Simultaneous regular and
non-regular production scheduling of multipurpose batch plants: A real chemical-
pharmaceutical case study”. Computers and Chemical Engineering, 67, 83–102.
Mula, J.; Lyons, A.C.; Hernández, J.E. y Poler, R. (2014): “An integer linear
programming model to support customer-driven material planning in synchronised,
multi-tier supply chains”. International Journal of Production Research, 52(14),
4267–4278.
Mula, J.; Poler, R. y García, J.P. (2006a): “MRP with flexible constraints: A fuzzy
mathematical programming approach”. Fuzzy Sets and Systems, 157(1), 74–97.
Mula, J.; Poler, R.; García-Sabater, G.S. y Lario, F.C. (2006b): “Models for production
planning under uncertainty: A review”. International Journal of Production
Economics, 103(1), 271–285.
Till, J.; Sand, G.; Urselmann, M. y Engell, S. (2007): “A hybrid evolutionary algorithm
for solving two-stage stochastic integer programs in chemical batch scheduling”.
Computers and Chemical Engineering, 31(5-6), 630–647.
168