## versión impresa ISSN 0103-6513

### Prod. vol.23 no.4 São Paulo oct./dic. 2013  Epub 14-Mayo-2013

#### http://dx.doi.org/10.1590/S0103-65132013005000023

Modelos lineares e não lineares inteiros para problemas da mochila bidimensional restrita a 2 estágios

Linear and nonlinear integer models for constrained two-stage two-dimensional knapsack problems

Horacio Hideki YanasseI, *; Reinaldo MorabitoII

Ihoracio@lac.inpe.br, INPE, Brasil
IImorabito@ufscar.br, UFSCar, Brasil

RESUMO

Palavras-chave: Problemas de corte e empacotamento. Mochila bidimensional. Corte guilhotinado-2 estágios. Modelos lineares e não lineares inteiros. Indústria de móveis.

ABSTRACT

In this work we review some linear and nonlinear integer models to generate two stage two-dimensional guillotine cutting patterns, including the constrained, non constrained, exact and non exact cases. These problems are particular cases of the two dimensional knapsack problems. We also present new models to generate these cutting patterns, based on adaptations and extensions of models that generate one-group constrained two dimensional cutting patterns. Two stage patterns arise in different cutting processes like, for instance, in the furniture industry and wooden hardboards. The models are useful for the research and development of more efficient methods, exploring particular structures, the model decomposition, model relaxations etc. They are also useful to evaluate the performance of heuristics, since they allow (at least for problems of moderate sizes) an estimative of the optimality gap of the solutions obtained by heuristics. To illustrate the application of the models we analyze the results of some computational experiments with instances of the literature and other generated randomly. The results were produced using a known commercial software and they show that the necessary computational effort to solve the models can be very different.

Keywords: Cutting and packing problems. Two-dimensional knapsack. Two-stage guillotine cut. Linear and nonlinear integer models. Furniture industry.

1. Introdução

Problemas de corte aparecem em muitos processos industriais em que rolos de papel e alumínio, chapas de vidro e fibra de vidro, barras e chapas de metais e madeira, pedaços de tecido ou de couro etc. são cortados para produzir uma certa quantidade de peças menores demandadas. O problema é determinar a melhor maneira de cortar os objetos maiores para produzir os itens demandados de modo a otimizar um ou mais objetivos, por exemplo, minimizar perdas de material. Para revisões e edições especiais em problemas de corte e empacotamento e suas aplicações industriais, os leitores podem consultar Dyckhoff e Waescher (1990), Lirov (1992), Dowsland e Dowsland (1992), Sweeney e Paternoster (1992), Dyckhoff e Finke (1992), Martello (1994a, b), Bischoff e Waescher (1995), Mukhacheva (1997), Dyckhoff, Scheithauer e Terno (1997), Arenales, Morabito e Yanasse (1999), Wang e Waescher (2002), Hifi (2002), Lodi, Martello e Monaci (2002), Waescher, Haussner e Schumann (2007), Morabito, Arenales e Yanasse (2009), Kendall, Daniels e Burke (2010) e ESICUP (EURO..., 2011).

Figura 1. Padrões de corte 2-estágios: (a) caso exato, (b) caso não exato.

Figura 2. Padrões de corte 1-grupo: (a) e (c) casos exatos, (b) e (d) casos não exatos.

Poucos trabalhos foram encontrados na literatura que apresentam formulações matemáticas para problemas de corte guilhotinados 2-estágios, vistos como casos particulares do problema da mochila bidimensional (WAESCHER; HAUSSNER; SCHUMANN, 2007). Modelos para tais problemas podem ser úteis para a pesquisa e desenvolvimento de métodos de solução mais efetivos, explorando características especiais e estruturas particulares, decomposição do modelo, relaxações do modelo etc. Esses modelos podem ser úteis também para a avaliação de desempenho de heurísticas, pois eles permitem (pelo menos para problemas de tamanho moderado) uma estimação do gap de otimalidade de suas soluções. Motivados por isso, neste estudo revemos modelos inteiros lineares e não lineares da literatura para gerar padrões 2-estágios e propomos novos modelos, incluindo os casos exato e não exato e restrito e irrestrito. Tais modelos também podem ser usados no procedimento de geração de colunas de Gilmore e Gomory, ou combinados com heurísticas de redução por repetição exaustiva (HINXMAN, 1980), para resolver problemas de corte de estoque. Além disso, eles podem ser combinados em formulações do problema bidimensional de empacotamento de bins (WAESCHER et al., 2007).

Este artigo está organizado da seguinte maneira: na seção 2, descrevemos brevemente os modelos lineares e não lineares existentes para o problema de corte 2-estágios. Na seção 3, propomos modelos lineares inteiros alternativos, alguns deles baseados na adaptação ou extensão de modelos de padrões de corte 1-grupo. Apesar da comparação efetiva dos modelos estar além do escopo deste trabalho, na seção 4, para efeito de ilustração, apresentamos os resultados de alguns experimentos computacionais da aplicação dos modelos para resolver alguns problemas exemplos da literatura e outros gerados de maneira aleatória. Utilizamos um software comercial bem conhecido de linguagem de modelagem GAMS e o resolvedor de programas lineares inteiros mistos CPLEX. Os resultados mostram que os esforços computacionais requeridos pelos modelos podem ser bastante diferentes, pelo menos para problemas de tamanho moderado. Finalmente, na seção 5, apresentamos as considerações finais e discutimos perspectivas de pesquisas futuras.

2. Revisão de modelos de corte 2-estágios

Nesta breve revisão admitimos, sem perda de generalidade, que os cortes do primeiro estágio são realizados paralelos ao comprimento da chapa, e os do segundo estágio são paralelos à largura da chapa. Admitimos também que a orientação das peças está fixada, ou seja, as peças não podem sofrer rotação. Sejam os seguintes parâmetros:

• L, W comprimento e largura da chapa
• m número total de tipos de peças
• li, wi comprimento e largura da peça tipo i (i = 1, ..., m)
• lmin, wmin menor comprimento e menor largura, respectivamente, das peças
• vi, bi valor (e.g., área) e demanda da peça tipo i, respectivamente.

Também sem perda de generalidade, admitimos que bi < , para i = 1, ..., m, uma vez que é um limitante superior para o número máximo de peças tipo i que podem ser cortados da chapa por um padrão 2-estágios guilhotinado. Dizemos que o problema é irrestrito se bi = para todo i; caso contrário é restrito.

2.1. Modelos inteiros não lineares

Vianna, Arenales e Gramani (2003) apresentaram o seguinte modelo inteiro não linear (VAG1) para o caso não exato. Sejam os parâmetros:

• K número de larguras diferentes wi.
• Pk número máximo possível de tiras L × wk (observe que Pk < )
• Ik = {i | wi < wk}

e as variáveis de decisão:

• número de vezes que a peça tipo i é cortada do p-ésimo padrão da tira L × wk
• µkp número de vezes que o p-ésimo padrão da tira L × wk é cortada ao longo da largura W da chapa

Sujeito a

com

A função objetivo (1) maximiza o valor total das peças cortadas no padrão. As restrições (2) impõem que o comprimento total das peças arranjadas nas tiras não exceda o comprimento da chapa. As restrições (3) impõem que a largura total das tiras arranjadas na placa não exceda a largura da chapa. As restrições (4) limitam o corte de peças a no máximo suas demandas. As restrições (5) referem-se à não negatividade e integralidade das variáveis. Note que sem as restrições (4), o modelo anterior corresponde ao caso do problema de corte 2-estágios irrestrito e ele, então, pode ser decomposto nos subproblemas explorados no método de duas fases de Gilmore e Gomory (1965), como mostrado em Vianna, Arenales e Gramani (2003).

O número de peças tipo i cortadas no padrão pode ser definido como:

Substituindo essa equação no modelo (1)-(5), obtemos o modelo VAG2:

Modelo VAG2

Sujeito a

com

Essa substituição tem o inconveniente de aumentar o número de variáveis e restrições, mas simplifica a apresentação da linearização apresentada a seguir. Observe que o caso de corte 2-estágios exato pode ser representado pelos modelos anteriores simplesmente redefinindo Ik = {i | wi = wk}. Como mostrado na seção 3, esses modelos inteiros não lineares podem ser reescritos como modelos lineares inteiros.

2.2. Modelos lineares inteiros

Lodi e Monaci (2003) apresentaram dois outros modelos lineares para o problema de corte 2-estágios não exato, baseados na restrição de empacotamento das peças em prateleiras (i.e., linhas formando níveis). Cada solução viável do problema é composta de prateleiras. Uma prateleira é uma tira (faixa) da chapa de comprimento L, com largura coincidente com a peça de maior largura cortada dela. A peça de maior largura de cada prateleira é dita ser a iniciante da prateleira. Numa solução viável do problema composta de prateleiras, cada item empacotado em uma prateleira pode ser cortado em no máximo 2 estágios (mais aparas).

No primeiro modelo (LM1) de Lodi e Monaci (2003), cada peça é considerada de maneira distinta. Assim, para cada peça tipo i, i = 1, ..., m, definimos bi itens distintos j, com lj = li, wj = wi, e vj = vi. Denotamos o número total de peças do problema por n, ou seja, e sem perda de generalidade, admitimos que w1 w2 ...wn. As variáveis de decisão do modelo são:

Note na definição anterior que j k porque a peça j só pode ser cortada (empacotada) de uma prateleira k cuja largura (wk) seja no mínimo igual à largura da peça j (wj), i.e., wk wj, j k. No modelo, admite-se que n prateleiras podem ser inicializadas em potencial, cada uma delas tendo como peça iniciante a peça k, k = 1, 2, ..., n. Assim, xkk = 1 (k = 1, ..., n) implica que a peça k é cortada da prateleira k, ou seja, a prateleira k é usada e a peça k é a peça iniciante dessa prateleira. O modelo LM1 então é definido por:

Modelo LM1

Sujeito a

com

A função objetivo (12) maximiza a soma dos valores das peças cortadas. As restrições (13) impõem que cada peça pode ser cortada no máximo uma vez, e somente de prateleiras cujas larguras sejam no mínimo iguais à sua largura. As restrições (14) impõem que as peças incluídas em cada prateleira não ultrapassem o seu comprimento. A restrição (15) impõe que as prateleiras usadas caibam na largura da chapa. As restrições (16) impõem que as variáveis de decisão sejam binárias.

No segundo modelo (LM2) de Lodi e Monaci (2003), as peças iguais são consideradas individualmente apenas nas prateleiras, para fins de sua inicialização. Define-se um mapeamento entre tipo de peças i, i = 1, ..., m, e prateleiras em potencial, k, k = 1, 2, ..., n, em que n é igual a como no modelo LM1. Sejam i = 1, ..., m, e α0 = 0; e βk, k = 1, 2, ..., n, o tipo de peça que inicializa a prateleira k, ou seja, para k = 1, 2, ..., n, βk = i sempre que αi - 1 + 1 < k < αi. Observe que se Wk é a largura da prateleira k, então para k = 1, 2, ..., n, Wk = wi sempre que αi - 1 + 1 < k < αi.

Admitindo-se novamente, sem perda de generalidade, que w1 w2 ... wm, qualquer peça do tipo i pode ser cortada de uma prateleira k, com k inteiro, 1 < k < αi, e qualquer prateleira k pode ser utilizada para cortar peças do tipo j, j inteiro, i < j < m, se αi - 1+1 < k < αi.

As variáveis de decisão (inteiras) do modelo LM2 são:

com i = 1, 2, ..., m; k = 1, 2, ..., αi.

O termo adicional da definição anterior indica que a peça tipo i que inicializa a prateleira é considerada à parte (se a prateleira é inicializada com esse tipo de item). Além dessas variáveis, temos as seguintes variáveis binárias:

Sejam Lk o comprimento da peça βk, ou seja, o comprimento da peça que inicializa a faixa k, e Wk a largura da faixa k que é igual à largura da peça βk, ou seja, a largura da peça que inicializa a faixa k, k = 1, ..., n. O modelo LM2 é definido por:

Modelo LM2

Sujeito a

com

A função objetivo (17) maximiza a soma dos valores das peças cortadas. As restrições (18) impõem que no máximo pode-se cortar a demanda de cada tipo de peça, e que essas peças somente podem ser cortadas de prateleiras cujas larguras sejam no mínimo iguais à largura da peça. As restrições (19) impõem que as peças incluídas em cada prateleira não ultrapassem o comprimento da chapa. A restrição (20) impõe que as prateleiras usadas caibam na largura da chapa. As restrições (21) impõem integralidade das variáveis xik, e, além disso, impõem que nas prateleiras de largura maior ou igual à largura de uma peça, o número de vezes que essa peça pode ser cortada é no máximo a sua demanda. As restrições (22) impõem que as variáveis qk sejam binárias.

Em adição às restrições (18)-(22), Lodi e Monaci (2003) incluem as seguintes restrições redundantes para fortalecer o modelo LM2:

Lodi e Monaci mostraram que os modelos LM1 e LM2 são competitivos comparados com o algoritmo exato de Hifi e Roucairol (2001), utilizando um método de resolução do tipo branch-and-bound do CPLEX, com a adição de algumas desigualdades lineares para evitar simetrias.

Para o modelo LM1, as seguintes desigualdades evitam algumas simetrias:

De maneira equivalente, para o modelo LM2, as seguintes desigualdades evitam algumas simetrias:

Os modelos LM1 e LM2 podem ser estendidos facilmente para o caso exato impondo (ou eliminando) que xjk = 0 se wk wj no modelo LM1, e Wk wj no modelo LM2.

3. Modelos lineares inteiros alternativos

3.1. Modelo linear inteiro 1

O modelo VAG2 (6)-(11) apresentado na seção 2.1 pode ser linearizado da seguinte forma (HARJUNKOSKI et al., 1997). Seja em que βkps-{0,1} e sk é tal que: ou seja, sk é o número mínimo de bits necessário para a representação binária de µkp (o mesmo poderia ser feito escolhendo-se ao invés de µkp; entretanto, o número de variáveis seria maior, na maioria dos casos). Com isso, a restrição não linear (9) pode ser reescrita como:

que, por sua vez, pode ser substituída pelo seguinte conjunto de restrições lineares:

em que M é um número suficientemente grande Observe que se βkps=1, então corretamente segue que por outro lado, se βkps = 0, então corretamente segue que O modelo (6)-(11) é reescrito como:

Modelo 1:

Sujeito a

com

Similarmente ao comentado na seção 2.1, o caso exato pode ser tratado pelo modelo 1 simplesmente redefinindo Ik = {i | wi = wk}.

É possível adicionar ao modelo 1 algumas desigualdades lineares para evitar simetrias:

Essas desigualdades apenas ordenam as tiras de mesma largura em um padrão, de forma não crescente de sua contribuição para a função objetivo.

A linearização anterior foi também explorada em Yanasse e Morabito (2006, 2008) para desenvolver modelos lineares inteiros (de formulações inteiras não lineares), para gerar padrões de corte guilhotinados p-grupo. Observamos que outras transformações 0-1 para as variáveis λj poderiam ser consideradas. Por exemplo, a transformação sugerida em Johnston e Sadinlija (2004) é interessante no caso em que existem valores específicos de λj que podem ser descartados a priori e o número de valores que λj pode assumir não é grande. Isto porque a transformação é pseudopolinomial, portanto, o número de variáveis 0-1 cresce mais rapidamente que a transformação apresentada anteriormente.

3.2. Modelo linear inteiro 2

Scheithauer (2002) apresentou um modelo linear inteiro para o problema de geração de padrões 1-grupo não exatos, que pode ser estendido para gerar padrões guilhotinados 2-estágios não exatos. Por conveniência, apresentamos a seguir o modelo mencionado de Scheithauer. Sejam os parâmetros:

• P,Q número máximo de tiras da esquerda para a direita e do fundo para o topo, respectivamente, no padrão

e as variáveis de decisão:

• Lj comprimento da j-ésima tira (da esquerda para a direita) no padrão (j = 1, ..., P)
• Wk largura da k-ésima tira (do fundo ao topo) no padrão (k = 1, ..., Q)

Sujeito a

com

O modelo (38)-(47) pode ser estendido para modelar o caso de padrões 2-estágios não exatos, redefinindo-se as variáveis xijk, que agora indicam se a peça tipo i é colocada na j-ésima posição da tira k (ao invés do retângulo Lj × Wk).

Modelo 2:

com

Para adaptar o modelo 2 para o caso exato, poderíamos pensar que seria suficiente substituir a desigualdade (52) pela igualdade:

No entanto, essa restrição é válida apenas no caso em que existe alguma peça do tipo i colocada na j-ésima posição da tira k. Se não existe nenhuma peça nessa posição, então essa restrição não precisa ser satisfeita. Para impor a restrição (56) apenas quando existe alguma peça do tipo i colocada na j-ésima posição da tira k, introduzimos as seguintes restrições:

em que M é um número suficientemente grande (e.g., M = W). Dessa forma, o modelo 2-estágios para o caso exato é obtido simplesmente substituindo as restrições (52) pelas restrições (57) e (58) no modelo 2.

3.3. Modelo linear inteiro 3

Um outro modelo linear pode ser proposto a partir do modelo 2, definindo-se variáveis inteiras xik para o número de peças tipo i na tira k, baseadas nas variáveis xijk do modelo 2, assim como variáveis binárias yik indicando se uma peça tipo i é alocada na tira k.

Variáveis:

• Wk largura da k-ésima tira (do fundo ao topo) no padrão (k = 1,..., Q)

ou seja, o número de peças tipo i na tira k,

Modelo 3:

com

em que M é um número suficientemente grande A descrição dessas restrições é similar à das restrições do modelo 2. Note que a função objetivo (59) e as restrições (60) e (63) são equivalentes a (48), (51) e (53) do modelo 2, com a substituição das variáveis As restrições (61) garantem que as larguras das peças em cada tira sejam menores ou iguais à largura da tira, similarmente às restrições (52) do modelo 2, e as restrições (62) e (66) são iguais às restrições (49) e (54) do modelo 2. As restrições (64) garantem que xik = 0 se yik = 0, e que yik = 1 se xik > 0. Já as restrições (65) não são necessárias no modelo; elas foram adicionadas apenas para garantir que yik = 0 se xik = 0. Finalmente as restrições (67) referem-se aos domínios das variáveis; note que Wk não precisam ser inteiras, assim como nas restrições (55) do modelo 2.

Similarmente ao modelo 2, para adaptar o modelo 3 para o caso exato, poder-se-ia pensar que seria suficiente substituir a desigualdade (61) pela igualdade:

Novamente, essa restrição é válida apenas se existe ao menos uma peça tipo i colocada na tira k. Para impor isso somente quando existe tal peça, introduzimos as seguintes restrições:

em que M é um número suficientemente grande (e.g., M = W). Portanto, o modelo para o caso exato 2-estágios é obtido substituindo as restrições (61) pelas restrições (69)-(70) no modelo 3.

Na Tabela 1 os modelos lineares inteiros LM1, LM2, 1, 2 e 3 são comparados em termos de número de variáveis e restrições (sem contabilizar as restrições para evitar simetrias).

4. Experimentos computacionais

Por ilustração, na Tabela 2 apresentam-se os dados de entrada de um problema simples bidimensional guilhotinado 2-estágios restrito analisado em Vianna et al. (2003). Na Tabela 3 resumimos os resultados computacionais obtidos com os modelos LM1, LM2, 1, 2 e 3 para ambos os casos, exato e não exato. Observe que todos os modelos encontram (e provam) a solução ótima do problema em tempos de execução (em segundos) relativamente pequenos. Na Figura 3 são apresentados os padrões ótimos encontrados pelos modelos para os casos não exato irrestrito (a), não exato restrito (b) e exato restrito (c).

Tabela 2. Exemplo restrito de Vianna et al. (2003) com (L, W) = (100, 100) com m = 5 tipos de peças e n = 22 peças.

Tabela 3. Resultados dos modelos LM1, LM2, 1, 2 e 3 (padrões 2-estágios) para o exemplo de Vianna, Arenales e Gramani (2003).

Figura 3. Padrões de corte do exemplo em Vianna et al. (2003): (a) caso não exato irrestrito (valor ótimo = 99,86), (b) caso não exato restrito (valor ótimo = 98,86), (c) caso exato restrito (valor ótimo = 91,26).

Para analisar o desempenho dos modelos em outros exemplos, geramos aleatoriamente 10 exemplares de problemas 2-estágios restritos, com chapa (L, W) = (100, 100) e m = 5 tipos de peças (li, wi), amostrados de distribuições uniformes no intervalo [0,1L, 0,5L] e [0,1W, 0,5W], respectivamente (depois de amostrado, li e wi foram simplesmente arredondados). As quantidades bi também foram amostradas (e depois arredondadas) de distribuições uniformes no intervalo [1, ï£°L/liï£»ï£°W/wiï£»]. Na Tabela 4 os resultados obtidos com os modelos LM1, LM2, 1, 2 e 3 para o caso não exato estão resumidos. Usando o GAMS/CPLEX, todos os modelos encontraram as soluções ótimas de todos os exemplares, mas o modelo LM1 não foi capaz de provar a otimalidade da solução do exemplo 2 (dentro de um limite de tempo imposto de 300 segundos - o número entre parênteses na linha do exemplo 2 corresponde ao gap entre o valor da melhor solução obtida e o valor do limitante superior neste tempo limite) e também teve dificuldades para resolver o exemplo 4. Note que os modelos LM2 e 3 tiveram um desempenho superior aos demais modelos.

A seguir, realizamos alguns experimentos com exemplares bem conhecidos da literatura de problemas restritos de corte e empacotamento com os modelos LM1, LM2, 1, 2 e 3, para obter suas soluções ótimas 2-estágios não exatos (Tabela 5). Os dados de entrada e as soluções (sem restrições de estágios) estão publicados em Christofides e Whitlock (1977) (exemplos CW1, CW2 e CW3 com m = 7, 10 e 20, respectivamente), Wang (1983) (exemplo W com m = 20) e Oliveira e Ferreira (1990) (exemplos OF1 e OF2 com m = 10). Esses exemplos também foram analisados em Viswanathan e Bagchi (1993), Christofides e Hadjiconstantinou (1995), Daza, Alvarenga e Diego (1995), Morabito e Arenales (1996) e Fayard, Hifi e Zissimopoulos (1998). Também utilizamos alguns exemplares irrestritos gcut* da OR-Library analisados em Beasley (1985) e Cintra et al. (2008) (os exemplos gcut1, gcut2 e gcut3 com (L, W) = (250, 250) e m = 10, 20, 30, respectivamente, o exemplo gcut12 com (L, W) = (1000, 1000) e m = 50, e o exemplo gcut13 com (L, W) = (3000, 3000) e (m = 32), para verificar o desempenho dos modelos em casos que a demanda bi de peças do tipo i é grande (i.e., bi = para todo i). Os demais exemplos gcut*w que aparecem na Tabela 5 correspondem aos problemas irrestritos gcut* com os respectivos comprimentos das peças li trocados pelas respectivas larguras wi, e vice-versa, as larguras wi pelos comprimentos li. Note que isso é equivalente a considerar nos problemas gcut* que os cortes do primeiro estágio são paralelos à largura W da chapa, e os cortes do segundo estágio são paralelos ao comprimento L da chapa. E o exemplo gcut13r é um problema irrestrito com (L, W) = (3000, 3000) e m = 64 tipos de peças (i.e., os 32 tipos de peças do exemplo gcut13, mais os 32 tipos de peças do exemplo gcut13w), ou seja, é equivalente a considerar no problema gcut13 que as peças podem sofrer rotações de 90º.

Observe na Tabela 5 que, usando o GAMS/CPLEX, os modelos LM2 e LM1 tiveram os melhores desempenhos nos problemas restritos, seguidos do modelo 3, modelo 2 e modelo 1. Nos problemas irrestritos, o modelo LM2 e o modelo 3 em geral tiveram os melhores desempenhos. No caso dos exemplos gcut13*, os modelos não foram capazes de encontrar (e provar) a solução ótima dentro do limite de tempo de 300 segundos (exceto para o modelo LM2 no exemplar gcut13w). Assim como na Tabela 4, o número entre parênteses na Tabela 5 corresponde ao gap entre o valor da melhor solução obtida (limitante inferior) e o valor do limitante superior neste tempo limite. Os valores das soluções ótimas dos exemplos gcut13 e gcut13r reportados na literatura são 8906216 e 8997780, respectivamente (BEASLEY, 1985; CINTRA ET AL., 2008); portanto, os gaps de otimalidade das soluções do modelo LM2 para esses problemas na Tabela 5 são bem pequenos (0,002% e 0,03%, respectivamente). Em particular, não foi possível resolver o modelo LM1 para o exemplo gcut13r por falta de memória computacional. Além disso, note que, ao contrário do que se poderia esperar, as soluções obtidas pelos modelos para o gcut13r no tempo limite são piores do que as soluções obtidas pelos mesmos modelos para o gcut13w no mesmo tempo limite. Verificou-se que os modelos precisam de mais tempo para obter soluções para o gcut13r tão boas ou melhores do que as obtidas para o gcut13w; por exemplo, o modelo LM2 precisa de mais 100 segundos para obter para o gcut13r a mesma solução obtida para o gcut13w.

As seis últimas linhas da Tabela 5 apresentam resultados adicionais obtidos com versões restritas dos problemas gcut1-3, gcut12-13, limitando-se nesses exemplos as quantidades de peças bi de cada tipo i = 1,.., m. Os exemplos gcut*ª correspondem aos respectivos exemplos gcut* com bi = 1, i = 1,.., m, enquanto o exemplo gcut13b corresponde ao exemplo gcut13 com bi = 10. Note que nesses problemas restritos os modelos LM2 e LM1 novamente tiveram os melhores desempenhos dentro de limite de tempo computacional, seguidos do modelo 3. Convém observar que ao se relaxar o limite de tempo computacional, os modelos LM1, LM2 e 3 não foram capazes de obter a solução ótima para o exemplo gcut13b devido à falta de memória computacional.

Finalmente, para ilustrar o desempenho dos modelos LM1, LM2, 1, 2 e 3 em um exemplo prático, consideramos um exemplo apresentado em Morabito e Arenales (2000) de uma indústria de móveis. Nas Tabelas 6 e 7 apresentam-se os dados de entrada e os resultados obtidos pelo GAMS/CPLEX para este exemplo para os casos exato e não exato (os valores das soluções na tabela correspondem às porcentagens das áreas totais das peças cortadas no padrão). Os gaps entre as melhores soluções obtidas (limitantes inferiores) e os limitantes superiores dos modelos no tempo limite de 300 segundos estão apresentados entre parênteses na Tabela 7. No caso não exato, a melhor solução foi obtida pelos modelos LM2, 1 e 3 dentro do limite de tempo (o gap de otimalidade do valor desta solução é menor que 0,3%). No caso exato, apenas o modelo 3 foi capaz de provar a otimalidade da solução obtida dentro do limite de tempo. Todos os outros modelos também encontraram essa mesma solução dentro do limite de tempo, exceto o modelo LM1 (que necessitou de mais de 600 segundos para encontrá-la).

5. Considerações finais

Neste estudo revemos modelos existentes e propomos novos modelos para gerar padrões guilhotinados restritos 2-estágios bidimensionais, incluindo os casos exato e não exato. Esses modelos são úteis para a pesquisa e desenvolvimento de métodos de solução mais eficientes, explorando características específicas e estruturas particulares, decomposição do modelo, relaxações do modelo etc. Tais modelos também são úteis para avaliar o desempenho de heurísticas, pois eles permitem (ao menos para problemas de tamanho moderado) estimar o gap de otimalidade de soluções heurísticas.

Uma comparação mais efetiva dos desempenhos dos modelos está além do escopo deste trabalho. No entanto, para ilustrar a aplicação dos modelos, apresentamos resultados de alguns experimentos computacionais usando um software comercial bem conhecido, a linguagem de modelagem GAMS e o otimizador CPLEX. Esses resultados mostraram que os esforços computacionais para solucionar os modelos podem ser bem diferentes. O modelo 3 (proposto neste estudo) obteve melhor desempenho que o modelo 1 (linearização proposta neste estudo do modelo não linear de Vianna, Arenales e Gramani (2003)) e o modelo 2 (extensão do modelo linear de SCHEITHAUER, 2002). Além disso, os modelos 1, 2 e 3 foram competitivos com os modelos LM1 e LM2 propostos em Lodi e Monaci (2003) nos experimentos aqui realizados, principalmente o modelo 3.

Cabe observar que para o caso em que a orientação das peças não está fixa é suficiente adicionar mais peças aos modelos 1, 2 e 3 correspondentes às peças originais rotacionadas em 90º. A restrição na limitação do número de peças é ajustada de acordo, isto é, limitamos a soma das peças rotacionadas mais as peças em sua orientação original. A extensão dos modelos 1, 2 e 3 para tratar problemas de corte bidimensionais guilhotinados restritos com mais do que 2 estágios é um tópico para pesquisa futura. Outra linha de pesquisa interessante seria estudar limitantes inferiores e superiores efetivos, restrições de simetria e desigualdades válidas específicas para os modelos 1, 2 e 3, com vistas a reduzir os tempos de execução necessários para resolvê-los.

Um comentário final é que para resolvermos o problema de corte de estoque com padrões bidimensionais 2-estágios, usando a abordagem de Gilmore e Gomory (1965) com geração de colunas, precisamos gerar um padrão de corte 2-estágios restrito a cada iteração do método simplex, isto é, resolver um problema da mochila bidimensional com 2 estágios. Assim, é grande o interesse em se ter um método eficiente de resolução do problema de geração destes padrões de cortes. Com os modelos 1, 2 e 3, pudemos resolver exemplares deste problema da mochila bidimensional com vários tipos de itens em poucos segundos. Além disso, notamos que em vários exemplares os tempos computacionais para se obter uma solução ótima são bem menores do que os tempos computacionais totais para se provar a otimalidade dessa solução. Uma perspectiva interessante para pesquisa futura seria explorar isso num procedimento de geração de colunas baseado nos modelos 1, 2 e 3 dentro da abordagem de Gilmore e Gomory, e comparar os resultados obtidos com outras abordagens conhecidas da literatura, como por exemplo as recentemente propostas em Macedo, Alves e Carvalho (2010) e Silva, Alvelos e Carvalho (2010).

Gostaríamos de agradecer a um dos revisores anônimos pelos úteis comentários e sugestões. Esta pesquisa foi parcialmente financiada pelo CNPq e pela Fapesp.

Referências

ALEM, D. J.; MORABITO, R. Production planning in furniture settings via robust optimization. Computers & Operations Research, v. 39, n. 2, p. 139-150, 2012. http://dx.doi.org/10.1016/j.cor.2011.02.022        [ Links ]

ARENALES, M.; MORABITO, R.; YANASSE, H. (Eds.). Special issue: Cutting and packing problems. Pesquisa Operacional, v. 19, n. 2, p. 107-299, 1999.         [ Links ]

BEASLEY, J. Algorithms for unconstrained two-dimensional guillotine cutting. Journal of the Operational Research Society, v. 36, p. 297-306, 1985.         [ Links ]

BELOV, G.; SCHEITHAUER, G. A branch-and-bound-and-price algorithm for one-dimensional stock cutting and two-dimensional two-stage cutting. European Journal of Operational Research, v. 171, p. 85-106, 2006. http://dx.doi.org/10.1016/j.ejor.2004.08.036        [ Links ]

BISCHOFF, E.; WAESCHER, G. (Eds.). Special issue: Cutting and packing. European Journal of Operational Research, v. 84, n. 3, 1995. http://dx.doi.org/10.1016/0377-2217(95)00018-L        [ Links ]

CARNIERI, C.; GUILLERMO, A.; GAVINHO, L. Solution procedures for cutting lumber into furniture parts. European Journal of Operational Research, v. 73, p. 495-501, 1994. http://dx.doi.org/10.1016/0377-2217(94)90244-5        [ Links ]

CHRISTOFIDES, N.; HADJICONSTANTINOU, E. An exact algorithm for orthogonal 2-D cutting problems using guillotine cuts. European Journal of Operational Research, v. 83, p. 21-38, 1995. http://dx.doi.org/10.1016/0377-2217(93)E0277-5        [ Links ]

CHRISTOFIDES, N.; WHITLOCK, C. An algorithm for two-dimensional cutting problems. Operations Research, v. 25, n. 1, p. 30-44, 1977. http://dx.doi.org/10.1287/opre.25.1.30        [ Links ]

CINTRA, G. F. et al. Algorithms for two-dimensional cutting stock and strip packing problems using dynamic programming and column generation. European Journal of Operational Research, v. 191, p. 61-85, 2008. http://dx.doi.org/10.1016/j.ejor.2007.08.007        [ Links ]

CUI, Y. An exact algorithm for generating homogeneous T-shape cutting patterns. Computers & Operations Research, v. 32, p. 143-152, 2005. http://dx.doi.org/10.1016/S0305-0548(03)00208-9        [ Links ]

DAZA, V. P.; ALVARENGA, A. G.; DIEGO, J. Exact solutions for constrained two-dimensional cutting problems. European Journal of Operational Research, v. 84, p. 633-644, 1995. http://dx.doi.org/10.1016/0377-2217(95)00028-O        [ Links ]

DOWSLAND, K.; DOWSLAND, W. Packing problems. European Journal of Operational Research, v. 56, p. 2-14, 1992 http://dx.doi.org/10.1016/0377-2217(92)90288-K        [ Links ]

DYCKHOFF, H. A new linear programming approach to the cutting stock problem. Operations Research, v. 29, n. 6, p. 1092-1104, 1981. http://dx.doi.org/10.1287/opre.29.6.1092        [ Links ]

DYCKHOFF, H.; FINKE, U. Cutting and packing in production and distribution: Typology and bibliography. Heidelberg: Springler-Verlag Co., 1992. http://dx.doi.org/10.1007/978-3-642-58165-6        [ Links ]

DYCKHOFF, H.; SCHEITHAUER, G.; TERNO, J. Cutting and packing. In: AMICO, M.; MAFFIOLI, F.; MARTELLO, S. (Eds.). Annotated bibliographies in combinatorial optimisation. New York: John Wiley & Sons, 1997. p. 393-414.         [ Links ]

DYCKHOFF, H.; WAESCHER, G. (Eds.). Special issue: Cutting and packing. European Journal of Operational Research, v. 44, n. 2, 1990.         [ Links ]

EURO SPECIAL INTEREST GROUP ON CUTTING AND PACKING - ESICUP. Apdio, 2011. Disponível em: <http://www.apdio.pt/sicup/>         [ Links ].

FARLEY, A. Practical adaptations of the Gilmore-Gomory approach to cutting stock problems. OR Spektrum, v. 10, p. 113-123, 1983. http://dx.doi.org/10.1007/BF01720210        [ Links ]

FAYARD, D.; HIFI, M.; ZISSIMOPOULOS, V. An efficient approach for large-scale two-dimensional guillotine cutting stock problems. Journal of the Operational Research Society, v. 49, p. 1270-1277, 1998.         [ Links ]

FORONDA, S.; CARINO, H. A heuristic approach to the lumber allocation and manufacturing in hardwood dimension and furniture manufacturing. European Journal of Operational Research, v. 54, p. 151-162, 1991. http://dx.doi.org/10.1016/0377-2217(91)90294-6        [ Links ]

GILMORE, P.; GOMORY, R. Multistage cutting stock problems of two and more dimensions. Operations Research, v. 14, p. 94-120, 1965. http://dx.doi.org/10.1287/opre.13.1.94        [ Links ]

GRAMANI, M.; FRANÇA, P. The combined cutting stock and lot-sizing problem in industrial processes. European Journal of Operational Research, v. 174, n. 1, p. 509-521, 2006. http://dx.doi.org/10.1016/j.ejor.2004.12.019        [ Links ]

GRAMANI, M. C. N. ; FRANÇA, P. M.; ARENALES, M. N. A lagrangean relaxation approach to a coupled lot-sizing and cutting stock problem. International Journal of Production Economics, v. 119, p. 219-227, 2009. http://dx.doi.org/10.1016/j.ijpe.2009.02.011        [ Links ]

HARJUNKOSKI, I. et al. Different strategies for solving bilinear integer non-linear programming problems with convex transformations. Computers & Chemical Engineering, v. 21, p. 487-492, 1997.         [ Links ]

HIFI, M. The DH/KD algorithm: A hybrid approach for unconstrained two-dimensional cutting problems. European Journal of Operational Research, v. 97, n. 1, p. 41-52, 1997. http://dx.doi.org/10.1016/S0377-2217(96)00060-4        [ Links ]

HIFI, M. (Ed.). Special issue on cutting and packing. Studia Informatica Universalis, v. 2, p. 1-161, 2002.         [ Links ]

HIFI, M.; ROUCAIROL, C. Approximate and exact algorithms for constrained (un)weighted two-dimensional two-staged cutting stock problems. Journal of Combinatorial Optimization, v. 5, p. 465-494, 2001. http://dx.doi.org/10.1023/A:1011628809603        [ Links ]

HINXMAN, A. The trim-loss and assortment problems: A survey. European Journal of Operational Research, v. 5, p. 8-18, 1980. http://dx.doi.org/10.1016/0377-2217(80)90068-5        [ Links ]

JOHNSTON, R. E.; SADINLIJA, E. A new model for complete solutions to one-dimensional cutting stock problems. European Journal of Operational Research, v. 153, p. 176-183, 2004. http://dx.doi.org/10.1016/S0377-2217(02)00704-X        [ Links ]

KENDALL, G.; DANIELS, K.; BURKE, E. K. Special issue: Cutting, packing, layout and space allocation. Annals of Operations Research, v. 179, n. 1, 2010.         [ Links ]

LIROV, Y. (Ed.). Special issue: Cutting stock: Geometric resource allocation. Mathematical and Computer Modelling, v. 16, n. 1, 1992.         [ Links ]

LODI, A.; MARTELLO, S.; MONACI, M. Two-dimensional packing problems: a survey. European Journal of Operational Research, v. 141, p. 241-252, 2002. http://dx.doi.org/10.1016/S0377-2217(02)00123-6        [ Links ]

LODI, A.; MONACI, M. Integer programming models for 2-staged two-dimensional knapack problems. Mathematical Programming, v. 94, p. 257-278, 2003. http://dx.doi.org/10.1007/s10107-002-0319-9        [ Links ]

MACEDO, R.; ALVES, C.; CARVALHO, J. M. V. Arc-flow model for two-dimensional guillotine cutting stock problem. Computers & Operations Research, v. 37, p. 991-1001, 2010. http://dx.doi.org/10.1016/j.cor.2009.08.005        [ Links ]

MARTELLO, S. (Ed.). Special issue: Knapsack, packing and cutting, Part I: One-dimensional knapsack problems. INFOR, v. 32, n. 3, 1994a.         [ Links ]

MARTELLO, S. (Ed.). Special issue: Knapsack, packing and cutting, Part II: Multidimensional knapsack and cutting stock problems. INFOR, v. 32, n. 4, 1994b.         [ Links ]

MORABITO, R.; ARENALES, M. Staged and constrained two-dimensional guillotine cutting problems: An and/or-graph approach. European Journal of Operational Research, v. 94, p. 548-560, 1996. http://dx.doi.org/10.1016/0377-2217(95)00128-X        [ Links ]

MORABITO, R.; ARENALES, M. Optimizing the cutting of stock plates in a furniture company. International Journal of Production Research, v. 38, n. 12, p. 2725-2742, 2000. http://dx.doi.org/10.1080/002075400411457        [ Links ]

MORABITO, R.; BELLUZZO, L. Optimising the cutting of wood fibre plates in the hardboard industry. European Journal of Operational Research, v. 183, p. 1405-1420, 2007. http://dx.doi.org/10.1016/j.ejor.2005.11.066        [ Links ]

MORABITO, R.; GARCIA, V. The cutting stock problem in a hardboard industry: A case study. Computers & Operations Research, v. 25, n. 6, p. 469-485, 1998. http://dx.doi.org/10.1016/S0305-0548(97)00087-7        [ Links ]

MORABITO, R.; ARENALES, M. N.; YANASSE, H. H. (Eds.). Special issue on Cutting, packing and related problems. International Transactions in Operations Research, v. 16, n. 6, 2009. http://dx.doi.org/10.1111/j.1475-3995.2009.00739.x        [ Links ]

MORABITO, R.; PUREZA, V. A heuristic approach based on dynamic programming and AND/OR-graph search for the constrained two-dimensional guillotine cutting problem. Annals of Operation Research, v. 179, p. 297-315, 2010. http://dx.doi.org/10.1007/s10479-008-0457-4        [ Links ]

MUKHACHEVA, E. A. (Ed.). Decision making under conditions of uncertainty: cutting -packing problems. Ufa: The International Scientific Collection, 1997.         [ Links ]

OLIVEIRA, J. F.; FERREIRA, J. S. An improved version of Wang's algorithm for two-dimensional cutting problems. European Journal of Operational Research, v. 44, p. 256-266, 1990. http://dx.doi.org/10.1016/0377-2217(90)90361-E        [ Links ]

POLDI, K. C.; ARENALES, M. N. Heurísticas para o problema de corte de estoque unidimensional inteiro. Pesquisa Operacional, v. 26, p. 473-492, 2006. http://dx.doi.org/10.1590/S0101-74382006000300003        [ Links ]

POLDI, K. C.; ARENALES, M. N. Heuristics for the one-dimensional cutting stock problem with limited multiple stock lengths. Computers & Operations Research, v. 36, p. 2074-2081, 2009. http://dx.doi.org/10.1016/j.cor.2008.07.001        [ Links ]

RANGEL, S.; FIGUEIREDO, A. G.; ALTAMIRO, G. O problema de corte de estoque em indústrias de móveis de pequeno e médio portes. Pesquisa Operacional, v. 28, p. 451-472, 2008. http://dx.doi.org/10.1590/S0101-74382008000300004        [ Links ]

RIEHME, J.; SCHEITHAUER, G.; TERNO, J. The solution of two-stage guillotine cutting stock problems having extremely varying order demands. European Journal of Operational Research, v. 91, p. 543-552, 1996. http://dx.doi.org/10.1016/0377-2217(95)00200-6        [ Links ]

SANTOS, S. G.; ARAUJO, S. A.; RANGEL, M. S. Integrated cutting machine programming and lot sizing in furniture industry. Pesquisa Operacional para o Desenvolvimento, v. 3, p. 249-166, 2011.         [ Links ]

SCHEITHAUER, G. On a two-dimensional guillotine cutting problem. In: INTERNATIONAL FEDERATION OF OPERATIONAL RESEARCH SOCIETIES - IFORS, 16., 2002, Edinburgh. Proceedings ... Edinburgh, 2002.         [ Links ]

SILVA, E.; ALVELOS, F.; CARVALHO, J. M. V. An integer programming model for two- and three-stage two-dimensional cutting stock problems. European Journal of Operational Research, v. 205, p. 699-708, 2010. http://dx.doi.org/10.1016/j.ejor.2010.01.039        [ Links ]

SWEENEY, P.; PATERNOSTER, E. Cutting and packing problems: A categorised, application-oriented research bibliography. Journal of the Operational Research Society, v. 43, p. 691-706, 1992.         [ Links ]

VIANNA, A. C. G.; ARENALES, M. N.; GRAMANI, M. C. N. Two-stage and constrained two-dimensional guillotine cutting problems. São Carlos: USP, 2003. p. 1-28. (Notas - Série Computação, n. 69).         [ Links ]

VISWANATHAN, K.; BAGCHI, A. Best-first search methods for constrained two-dimensional cutting stock problems. Operations Research, v. 41, n. 4, p. 768-776, 1993. http://dx.doi.org/10.1287/opre.41.4.768        [ Links ]

YANASSE, H. H.; KATSURAYAMA, D. M. Checkerboard patterns: proposals for its generation. International Transactions in Operational Research, v. 12, p. 21-45, 2005. http://dx.doi.org/10.1111/j.1475-3995.2005.00488.x        [ Links ]

YANASSE, H. H.; ZINOBER, A.; HARRIS, R. Two-dimensional cutting stock with multiple stock sizes. Journal of the Operational Research Society, v. 42, n. 8, p. 673-683. 1991.         [ Links ]

YANASSE, H. H.; MORABITO, R. Linear models for one-group two-dimensional guillotine cutting problems. International Journal of Production Research, v. 44, n. 17, p. 3471-3491, 2006. http://dx.doi.org/10.1080/00207540500478603        [ Links ]

YANASSE, H. H.; MORABITO, R. A note on linear models for two-group and three-group two-dimensional guillotine cutting problems. International Journal of Production Research, v. 46, n. 21, p. 6189-6206, 2008. http://dx.doi.org/10.1080/00207540601011543        [ Links ]

WAESCHER, G.; GAU, T. Heuristics for the integer one-dimensional cutting stock problem: a computational study. OR Specktrum, v. 18, p. 131-144, 1996. http://dx.doi.org/10.1007/BF01539705        [ Links ]

WAESCHER, G.; HAUSSNER, H.; SCHUMANN, H. An improved typology of cutting and packing problems. European Journal of Operational Research, v. 183, p. 1109-1130, 2007. http://dx.doi.org/10.1016/j.ejor.2005.12.047        [ Links ]

WANG, P. Two algorithms for constrained two-dimensional cutting stock problems. Operations Research, v. 31, p. 573-586, 1983. http://dx.doi.org/10.1287/opre.31.3.573        [ Links ]

WANG, P.; WAESCHER, G. (Eds.). Special issue on cutting and packing problems. European Journal of Operational Research, v. 141, p. 239-469, 2002. http://dx.doi.org/10.1016/S0377-2217(02)00122-4        [ Links ]

Recebido 01/06/2011
Aceito 27/09/2012

* INPE, São José dos Campos, SP, Brasil

Todo el contenido de esta revista, excepto dónde está identificado, está bajo una Licencia Creative Commons