journal.gif (7158 bytes)

 

Evaluación del Desempeño por Simulación de un Sistemas Conducido por Eventos Discretos


ContenidoSiguiente artículo

Ramiro Fuentes Ayala

Facultad de Ciencia y Tecnología 
- Maestría en Control Moderno de Sistemas

Universidad Mayor de San Simón

Cochabamba- Bolivia

Resumen - Este trabajo tiene como objetivo desarrollar un modelo de simulación del sistema de transporte de cajas de la sección aparado de la fabrica de calzados MANACO y evaluar el desempeño de la banda transportadora (circular).

El desarrollo del proyecto esta dividido en las siguientes etapas: Definición del Problema y objetivos; Desarrollo del modelo conceptual y recolección de datos; Codificación y Verificación; Validación del modelo; Proyecto experimental; Ejecución y análisis de los resultados. Como resultado se obtuvo un modelo de simulación del sistema de transporte implementado en el lenguaje de inscripciones RENEW. 

Para evaluar el desempeño de la banda transportadora se evaluó varios índices de desempeño como; numero de cajas en los buffers, tiempo en que las maquinas esperan por materia prima, numero de cajas por unidad de tiempo que son transportados, influencia de la velocidad de la circular en su desempeño. 

 

 

Introducción

 

 

La complejidad de los sistemas de manufactura debido al avance tecnológico y a la demanda del mercado, requiere de un constante análisis con el objetivo de mejorar el desempeño del sistema de producción. Dicho análisis debe dar respuestas inmediatas, para lo cual el uso del computador y las técnicas de simulación son bastante eficientes en este propósito. Tomando en cuenta los riesgos que representa el realizar pruebas con el sistema real, las pruebas pasan a ser ejecutados en un modelo simulado por computador que debe representar adecuadamente el comportamiento del sistema en estudio. Mediante esta técnica nuevas alternativas pueden ser verificadas antes de ser implementadas en el sistema real. Los métodos de análisis pueden ser divididos en dos grandes clases: Métodos Analíticos y Simulación.

 El método adoptado fue el método por Simulación debido a la complejidad y aleatoriedades del sistema en estudio.

Problema - La línea de producción de la sección aparado de la fábrica de calzados Manaco, presenta fluctuaciones en el número de pares producidos. También se ha observado que los puestos de trabajo presentan periodos de espera por la materia prima que es transportado mediante una banda transportadora (circular).

Objeto - Proceso de transporte vía banda transportadora de la materia prima en la sección aparado de la fabrica de calzados Manaco.

Objetivos - Medir el desempeño de la banda transportadora por simulación para contribuir al mejoramiento del sistema de producción de la sección aparado.

Campo de acción - Modelación por simulación de sistemas de transporte de materia prima en la industria de manufactura.

Diagnostico - Actualmente la sección 422 (sección aparado del calzado) de la fabrica de calzados MANACO, tiene una producción de 1000 pares por día. Esta sección tiene una línea de producción que utiliza una banda transportadora (circular) para transportar cajas (conteniendo cortes de cuero para el aparado), hacia puestos de trabajo mediante compuertas controladas remotamente por un operador quien está encargado de distribuir dichas cajas. Dicha distribución se la realiza por demanda del puesto de trabajo, esto es, el buffer del puesto de trabajo (capacidad cuatro cajas) activa una lámpara en un tablero de mando cuando el buffer puede recibir otra caja (o dos). Esta lámpara es observada por el  encargado y procede al envío.

Existen varios factores que pueden influenciar en la sincronía de transporte y demanda de materia prima en los puestos de trabajo que influencia en la producción.

Tendencias: 

Métodos

 

 

El desarrollo del proyecto siguió las siguientes etapas: Definición del Problema y objetivos; Desarrollo del modelo conceptual y recolección de datos; Codificación y Verificación; Validación del modelo; Proyecto experimental; Ejecución y análisis de los resultados. Como resultado se obtuvo un modelo de simulación del sistema de transporte implementado en el lenguaje de inscripciones RENEW. Una vez validado el modelo se procedió a los experimentos para responder a los objetivos.

 

Diseño -

Para diseñar un modelo de simulación se desarrollo un modelo conceptual, a partir del análisis del sistema real tomando en cuenta los conceptos; modelado, simulación, eventos discretos y las redes de petri de alto nivel. El resultado obtenido es un modelo de simulación que fue codificado en un simulador de redes de petri de alto nivel, el cual usa un lenguaje de inscripciones. La figura 1.-  muestra la red que implementa la línea de producción de la sección aparado, que consta de un buffer de distribución (BD), circular_ida, circular_vuelta, maquinas (puestos de trabajo). Cada uno de estos recursos está asociado a lugares en la rede de petri. El lugar BD esta implementado por una red (que no se muestra)  y es utilizado en tiempo de ejecución, dicha red implementa el envió de cajas hacia la circular_ida  que lleva cajas hacia las maquinas que a su vez envía cajas hacia la circular_vuelta que devuelve cajas hacia el buffer de distribución. Esta red permite describir tanto la asociación de los recursos, la dinámica, y la funcionalidad del sistema. Todos estos recursos (buffer de distribución, circular, maquina) están implementados en otras redes que son invocadas en tiempo de ejecución.

 

Datos de entrada para la validación del modelo:

Frecuencia de envió a maquinas - Es la frecuencia con la que una maquina requiere cajas. Esto se muestra en la tabla 1.-

Para generar esta frecuencia se asigno un rango de valores a cada maquina destino, dichos valores sumados deben dar “1”. Luego se generan valores aleatorios entre 0 y 1 y se determina que maquina es destino. Y luego se calcula el tiempo que demora la caja en llegar a su destino, tomando en cuenta que son datos; la distancia de las maquinas y la velocidad de la circular.

Intervalo de tiempo de envío - Es el intervalo de tiempo en que se produce un envío (una caja o dos cajas apiladas). Se generan valores aleatorios siguiendo una distribución exponencial. Esto fue determinado después de un tratamiento estadístico de los datos recolectados. La prueba de Kolmogrov Smirnov (K-S)  mostró que esta distribución se aproxima mejor a una exponencial. Como se muestra en la tabla 2.- y la figura 2.-

 

Figura 1. Codificación del modelo de simulación

 

Velocidad de la banda transportadora (circular) - Es la velocidad tanto de ida como de vuelta, ya que por la parte superior se envían  y por la parte inferior retornan las cajas.

 

 Tiempo de procesamiento de las maquinas - Es el tiempo en que una maquina termina de procesar la materia prima de una caja. Estos datos fueron recolectados durante dos meses en dos turnos para cada maquina. Y según la prueba K-S estos datos se aproximaron a distribuciones normales mejor que otras. Como se puede observar para 10 maquinas en la tabla 3.- y en la figura 3.-

 

N

maq

frec.

fre_re

1

m4

6

1,8%

2

m6

36

10,8%

3

m8

20

6,0%

4

m12

12

3,6%

5

m14

17

5,1%

6

m16

4

1,2%

7

m18

15

4,5%

8

m20

16

4,8%

9

m26

15

4,5%

10

m27

3

0,9%

11

m28

13

3,9%

12

m29

5

1,5%

13

m30

13

3,9%

14

m32

6

1,8%

15

m34

6

1,8%

16

m36

23

6,9%

17

m37

3

0,9%

18

m38

22

6,6%

19

m41

4

1,2%

20

m42

6

1,8%

21

m43

9

2,7%

22

m44

5

1,5%

23

m45

4

1,2%

24

m46

4

1,2%

25

m47

21

6,3%

26

m49

5

1,5%

27

m51

3

0,9%

28

m53

16

4,8%

29

m54

22

6,6%

 

total

334

100,0%

Tabla 1 Frecuencia de envió de cajas

 

 

 

 

 

 

 Figura 2.- Distribución de frecuencias de envío de cajas

 

 

 

 

M27

M26

M20

M18

M14

M12

M16

M8

M6

M4

N

 

7

7

7

7

7

7

7

7

7

7

Normal Parameters

Mean

8,431

14,952

16,384

8,169

6,220

6,608

12,732

15,809

5,377

8,101

Std. Deviation  

1,538

12,901

6,484

2,575

1,600

3,526

4,865

4,584

1,467

6,200

Most Extreme Differences

Absolute

,191

,342

,188

,195

,155

,277

,253

,310

,210

,305

 

Positive

,118

,342

,188

,195

,153

,277

,253

,212

,210

,305

 

Negative

-,191

-,188

-,156

-,178

-,155

-,164

-,178

-,310

-,141

-,176

Kolmogorov-

Smirnov Z

 

,506

,904

,497

,515

,409

,733

,668

,819

,554

,807

Asymp. Sig. (2-tailed)

 

,960

,387

,966

,953

,996

,655

,763

,513

,918

,533

 

Tabla 3.- Prueba de la Distribución normal de la tasa de producción algunas maquinas

 

 

Figura 3.- Histogramas de los tiempos de procesamiento la maquina m6

 

 

 

Salida de datos:

Tiempo de retorno de cajas al buffer de distribución - Es el intervalo de tiempo entre una caja y otra que retorna. La tabla 4.- describe los datos estadísticos del sistema  real (real1, real2, real3) y el simulado (sim1,sim2, sim3).

 

 

 

Tabla 4.- Promedio y desviación de las muestras del sistema real y simulado

 

 

Paired Samples Test

 

 

Paired Differences

 

 

 

 

t

df

Sig. (2-tailed)

 

 

Mean

Std. Deviation

Std. Error Mean

99% Confidence Interval of the Difference

 

 

 

 

 

 

 

 

 

Lower

Upper

 

 

 

Pair 1

REAL1 - SIM1

-,35167

1,29355

,28925

-1,17918

,47585

-1,216

19

,239

Pair 2

REAL2 - SIM2

,41683

1,02566

,22934

-,23931

1,07297

1,817

19

,085

Pair 3

REAL3 - SIM3

,35317

1,25933

,28159

-,45245

1,15879

1,254

19

,225

 

 

Tabla 5.- Comparación de los intervalos de tiempo de retorno de cajas, para el sistema real y el simulado.

 

La tabla 5.- muestra resultado de la comparación de las dos muestras. Como se puede observar los valores críticos “t” están dentro del rango –2.43 y +2.43 los cuales fueron extraídos de la tabla estadística.

Por lo tanto, el modelo simulado representa al sistema real y se acepta la hipótesis H0 lo que significa que las dos muestras pertenecen a una misma población.

Una vez validado el modelo se procedió a los experimentos para responder a los objetivos del proyecto.

Datos de Entrada del Modelo de Simulación - El modelo de simulación del sistema permite que el proyectista establezca los datos (parámetros de entrada) de interés con el objetivo de evaluar el sistema de acuerdo a las condiciones dadas. La entrada de datos principales del modelo pueden ser ingresados abriendo el archivo “enviacajas.rnw” que corresponde a la rede mostrada en la figura 4.-. el ingreso de datos se lo realiza modificando los parámetros del lugar maquinas [maq, dist, med, sigma, lim1, lim2]

 

 

 

 

Figura 4.- Modelo de envío de cajas codificado en el simulador RENEW.

 

Significado de los parámetros:

 “maq”  permite ingresar el número asignado a la máquina,

“dist” permite ingresar la distancia a la que se encuentra la máquina respecto del buffer de distribución,

“media” permite ingresar el valor medio de producción de la máquina, en cajas por minuto.

“sigma” permite ingresar el valor de desviación estándar del valor medio de producción de la maquina.

“lim1” permite ingresar el limite inferior del rango de clase de la máquina

“lim2” permite ingresar el limite superior del rango de clase de la máquina

Si el usuario no modifica los valores anteriores se usan los valores por defecto. El tiempo de simulación es mostrado continuamente por el simulador, esto no es necesario especificarlo, solo se observa que haya llegado al tiempo requerido. La velocidad de la banda transportadora puede ser modificada abriendo los archivos de redes “circularida.rnw” y “circularvuelta.rnw” . Por defecto se encuentra con velocidad nominal de 70.2 m/min (velocidad actual).

Datos para el experimento:

 

Numero de máquinas: 29

Velocidad de la banda transportadora (circular): 70.2m/min

Tiempo de procesamiento de las maquinas: Fueron utilizados los mismos datos que fueron utilizados para la validación del modelo.

 

Resultados y Discusión

 

 

Interpretación de los resultados usando RENEW - A cada instante se muestran los

estados en que se encuentran la banda transportadora, los buffers de las maquinas, las maquinas, numero de cajas  transportadas por la  banda transportadora, las cajas y el    tiempo que permanecen en el sistema.

El simulador no cuenta con el recurso que permita recolectar datos en un archivo, pero es posible implementarlo llamando métodos de java. Esto se implementó para el caso de la recolección de los datos de salida de la banda transportadora de vuelta. 

Para medir el tiempo en que la banda transportadora  (BT) permanece con cajas, se procedió a implementar una rede, que acumula el tiempo durante la simulación, conociendo este tiempo mas el tiempo de simulación, la resta de estos valores determina el tiempo en que la circular (arriba y abajo) permanece sin cajas.

Es posible observar el numero de cajas en cada buffer de maquina, simplemente abriendo el lugar buffer_maquina, donde se encuentran las fichas; [maq, dist, med, sigma] @tiempo ej:2’[6,2.40,5.377,1.467]@ 453.645 especifica dos cajas en la maquina m6, que se encuentra a una distancia de 2.4 metros, con media de producción de 5.377 y desviación estándar 1.467 con un tiempo 453.645 minutos en el sistema.

Ejecución y análisis de los resultados del modelo de simulación -  Como el objetivo

principal es medir el desempeño de la BT, condiciones iniciales fueron establecidos. Se

ejecuto el modelo simulado (siguiendo la metodología de sistema no terminales) con los

mismo valores utilizados para la validación, con el objetivo de determinar el porcentaje

de uso de la banda transportadora respecto de su capacidad total, ya que estas son

las condiciones reales de la banda transportadora. Esto se calculo de la siguiente

manera: (tiempo en que la banda esta vacía)/ (tiempo de simulación). Esto dio un

resultado de 56,78%, lo que quiere decir que la banda esta siendo utilizado el 56,78 %

de su capacidad nominal. 

Para verificar que la banda transportador es capaz de atender a los 29 puestos de trabajo, asumiendo que nunca falta material en el buffer, se utilizo las siguientes condiciones iniciales:  29 puestos de trabajo, con envío de cajas controlado por la demanda de una caja en el buffer (esto es, si la maquina termina de procesar una caja, se envía otra caja). Para lo cual se implementó una red para estas nuevas condiciones usando las redes para cada recurso anteriormente codificadas.

El experimento de simulación fue ejecutado tomando en cuenta datos de entrada como el promedio, desviación de los tiempos de procesamiento de las 29 maquinas.

Adoptando la metodología de sistemas no terminales  se siguió los siguientes pasos:

1.- Se asumió que el sistema es no terminal.

2.- El tiempo de “warmup”  fue determinado efectuando una simulación con un tiempo de 480 min, esto es una jornada de trabajo y observando el instante en que los buffers de las maquinas tengan por lo menos una caja en promedio. Dado un tiempo de warmup igual a 20min

3.- Como el simulador permite observar en tiempo de ejecución el numero de maquinas que están en operación y cuantas cajas se encuentran en buffer, fue posible determinar el instante en el que los buffer tienen por lo menos una caja y a partir de esto se consideró el estado de régimen.

4.- También se observó la correlación de los datos de salida para garantizar la independencia de los mismos.

Media de la muestra        3,7 min

Varianza                         32,7min       

Correlación                      0.0067

Num. de observaciones  40

 

5.- Se seleccionó el tamaño de los grupos igual a 10 veces el numero de observaciones, lo cual dio 400 muestras, entonces cada grupo formado tuvo 400 observaciones.

El numero de grupos formados fue de 20, representando 20 muestras.

Tomando en cuenta que en promedio de cajas de salida es  1.33 min, el tiempo final de simulación fue definido como:

Tiempo de simulación = “warmup”+(10grupos)*(400observaciones)*(1.33min) = 5340min (3.6 días de producción).

Interpretación y análisis de los Resultados de la Simulación - El numero de cajas transportadas de ida y vuelta  por minuto, en promedio es mostrado en el lugar cajas_ida (usando el  tiempo de simulación se calcula 3,05 cajas/min) y cajas_vuelta(2,4 cajas/min). Tomando en cuenta que la banda puede enviar dos cajas apiladas, la esperanza matemática nos dice que el numero de cajas transportadas de ida es:  3,05* (porcentaje de envió de una caja)+3,05* (porcentaje de envió de dos cajas) = 3,05*0.41+6,1*0.59 = 4,34 cajas/min. Esto significa que es la capacidad de transporte de la circular para atender la demanda de los puestos de trabajo. Numero de cajas promedio acumulados en los buffer de la maquinas, en promedio es de 3 lo cual nos indica que la banda transportadora atiende satisfactoriamente la demanda de los puestos de trabajo. Entonces la circular no es un factor para la ocurrencia de perjuicios en los puestos de trabajo. Para una velocidad de 40,2 m/min de la banda transportadora es posible atender satisfactoriamente, puesto que aun así es  posible entender a las maquinas. Debajo de esta velocidad aumenta el promedio de máquinas en espera.  Este resultado puede darnos nuevas alternativas de transporte mas lentos.

También se ha observado que la banda transportador esta siendo utilizada el 56,78% de su capacidad nominal.

 

Conclusión - Tomando en cuenta las grandes ventajas de utilizar la simulación como herramienta de proyecto, se puede aplicar en el desenvolvimiento de proyectos industriales para determinar estrategias de funcionamiento. Este trabajo tuvo como objetivo construir un modelo de simulación para evaluar el desempeño de la banda transportador de la sección aparado de la Fábrica de Calzados MANACO, lo cual fue cumplido con éxito.

Para lograr este objetivo el proyecto tuvo las siguientes etapas:

 

 

Contribuciones de este proyecto - Modelo del sistema de transporte de la sección aparado.

Conocimiento de todo el funcionamiento del sistema, Conocimiento de la influencia de la velocidad de la circular en la producción, Posibilidad de ejecutar nuevos proyectos a partir del modelo de simulación.

Perspectivas - Si bien el presente proyecto no soluciona el problema de tiempos de espera en los puestos de trabajo en su totalidad, ya que esto esta ligado a otros factores, pero ofrece una gran contribución a solucionar este problema. Por que basado en el modelo de simulación obtenido, es posible implementar experimentos que incluyan detalles como la secuencias de operaciones para cada tipo de modelo de calzado y de esta manera conocer cual es la secuencia de envíos y cuantas maquinas se asignaría a cada proceso del modelo de calzado, con el objetivo de determinar la mejor secuencia de envíos que no  produce tiempos de espera en las maquinas. Es posible realizar experimentos; como ser nuevas alternativas de transporte; Integración del sistema con un sistema real ;Todas estas alternativas pueden ser respondidas a través de un experimento de simulación.

 

Referencias Bibliográficas

 

[1]      N.Viswanadham, Y. Narahari, by Prentice Hall, Performance Modeling of Automated Manufacturing Systems, 1992

[2]      Agencia Española  de Cooperación  internacional, “I Curso Iberoamericano de Automática, modulo II (Las Redes de Petri)”, 4 al 8 de diciembre de 1995

[3]      G.W. Brams edición en español por Ma Pilar Nivela Alos, “Las redes de Petri Teoría y practica, tomo I y II ”, 1986

[4]      Jerry Banks y Johs S. Carson, II, “Discrete-Event System Simulation”, 1994

[5]      Paul A. Fishwick, “ The Role of Process Abastraction in Simulation” IEEE Transactions on Systems, Man And Cybernetics,Volume 18, Enero/Febrero de 1998, paginas 18  a 39

[6]      Prentice Hall  James Rumbouh, Michel Blaha, Willam Premerlani, Fredeerick Eddy, Modelamiento y Diseño Orientado a Objetos.1991

[7]      Osborne/Mc Graw-Hill Patrick Naughton “Manual JAVA”

[8]      Miller Irwin, Freund “Probabilidad y Estadistica para Ingenieros”

[9]      Osborne/McGraw-Hill Edicion en español Patrick Nughton “Manual de Java  1996

[10]    Editora da DAUFSC Janette Cardoso, Robert Valente”Redes de Petri” 1997