INVESTIGAÇÃO OPERACIONAL 7ª Aula Método dos Transportes e de Distribuição O método dos transportes é um dos métodos de programação linear e deve o nome à sua aplicação em problemas que envolvem a optimização do transporte de bens. O método de distribuição é outro método de programação linear destinado à alocação (ou distribuição) de pessoas por tarefas, podendo-se considerar um tipo de problemas de transportes. As aplicações relativas a problemas de Transportes e Alocação envolvem inúmeras variáveis de decisão e restrições. No entanto, uma grande parte dos coeficientes das variáveis nas restrições são zero, o que permite que as simplificações introduzidas pelo método dos transportes levem a um menor volume de cálculo. Cecília Rocha # 1 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula Exercício- Exemplo Suponha que Inglaterra, França e Espanha produzem todo o trigo, cevada e aveia disponível no mundo. A procura mundial de trigo corresponde à produção de 125 milhões de acres de solo. Com o mesmo objectivo são necessários 60 milhões de acres para cevada e 75 milhões de acres para aveia. O total de solo agrícola disponível para este propósito, em Inglaterra, França e Espanha é de, respectivamente, 70 milhões, 110 milhões e 80 milhões de acres. O número de horas de trabalho necessárias para produzir 1 acre de trigo é de 18h em Inglaterra, 13 em França e 16 em Espanha. No caso do cevada são necessárias 15h em Inglaterra e 12h em França e em Espanha. Para o aveia são precisas 12h em Inglaterra, 10 em França e 16 em Espanha. O custo da hora de trabalho para produção de trigo é de 3 u.m., 2.4 u.m. e 3.3 u.m., respectivamente em Inglaterra, França e Espanha. Para a produção de cevada o custo da hora de trabalho será de 2.7 u.m., 3.0 u.m. e 2.8 u.m. em Inglaterra, França e Espanha. No caso da aveia haverá um custo da hora de trabalho de 2.3 u.m. em Inglaterra, 2.5 u.m. em França e 2.1 u.m. em Espanha. O problema é definir a melhor distribuição da produção em cada país, de forma a satisfazer as necessidades mundiais de trigo, cevada e aveia mas minimizando o custo de produção total. a) Formular este problema como um Problema de Transportes, construindo o quadro de custos e requisitos; b) Utilize uma rotina automática do SOLVER para encontrar um solução óptima para o problema; c) Utilize o método de Vogel e o método do Custo Mínimo para determinar uma solução básica admissível inicial; d) Resolva pelo Método dos Transportes. Cecília Rocha # 2 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) [70] [110] [80] Cecília Rocha # 3 Inglaterra França Espanha x11 18*3.0 u.m. = 54 x12 15*2.7 u.m. = 40.5 x13 12*2.3 u.m. = 27.6 x21 13*2.4 u.m. = 31.2 x22 12*3.0 u.m. = 36 x23 10*2.5 u.m. = 25 x31 16*3.3 u.m. = 52.8 x33 16*2.10 u.m. = 33.6 Trigo Cevada Aveia [125] [60] [75] 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Formulação do Problema Minimizar Z = 54 x11 + 40.5 x12 + 27.6 x13 + 31.2 x21 + 36 x22 + 25 x23 + 52.8 x31 + 33.6 x32 + 33.6 x33 s.a.: Restrições de Origem x11 + x12 + x13 = 70 x21 + x22 + x23 = 110 x31 + x32 + x33 = 80 Restrições de Destino x11 + x21 x12 + x31 + x22 x13 + x23 = 120 + x32 = 60 + x33 = 75 O facto de todos os coeficientes das variáveis terem o valor UM e estarem dispostos da forma evidenciada neste problema, é que permite distinguir este problema dos restantes e não a sua aplicabilidade a este tipo específico de problemas. Pode ser aplicado noutras situações em que esta estrutura se repita. Cecília Rocha # 4 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Modelo do Problema de Transportes Genericamente, o problema de transportes refere-se à distribuição de qualquer tipo de bem, proveniente de um conjunto de fornecedores – denominados Origens – para um conjunto de clientes – denominados Destinos. Assim, uma origem i (i = 1, 2, ..., m) pode fornecer si unidades de um dado bem por vários destinos e um dado destino j (j = 1, 2, ..., n) tem uma procura de dj unidades a receber das origens. Um pressuposto básico é que o custo de distribuição entre a origem i e o destino j é proporcional ao número de unidades transportadas, sendo ci j o custo de distribuição unitário. Para o exercício exemplo apresentado o Quadro de Custos e Requisitos é o seguinte: Destinos Cevada Aveia 18 * 3.0 = 54 15 * 2.7 = 40.5 12 * 2.3 = 27.6 70 França 13 * 2.4 = 31.2 12 * 3.0 = 36.0 10 * 2.5 = 25.0 110 Espanha 16 * 3.3 = 52.8 12 * 2.8 = 33.6 16 * 2.1 = 33.6 80 125 60 75 Inglaterra Origens Procura Cecília Rocha # 5 Oferta Trigo 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Modelo do Problema de Transportes Genericamente, o Quadro de Custos e Requisitos é o seguinte: Destinos D2 ... Dn S1 c11 c12 ... c1n s1 S2 c21 c22 ... c2n s2 Sm cm1 cm2 ... cmn sm Procura d1 d2 ... dn Origens Oferta D1 ... A formulação deste tipo de problema é a seguinte: m n Minimizar Sujeito a: cij xij i 1 j 1 n xij si j 1 Cecília Rocha # 6 m e xij d j i 1 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Propriedades das Soluções nos Problemas de Transportes Solução Inteira: Para os problemas de transportes, em que a procura e oferta têm valores inteiros, todas as variáveis básicas em todas as Soluções Básicas Admissíveis (incluindo a óptima) também têm valor inteiro. Solução Admissível: É condição necessária e suficiente para que um problema de transportes tenha alguma solução admissível que: m n i 1 j 1 si dj Esta propriedade pode ser analisada observando a seguinte situação m n m n i 1 j 1 i 1 j 1 si dj xi j Cecília Rocha # 7 Esta condição “A oferta total tem de igualar a procura total” implica que, nos casos em que tal não aconteça, seja necessário adoptar uma Origem Fictícia ou um Destino Fictício, consoante exista maior procura ou maior oferta, respectivamente. 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Propriedades das Soluções nos Problemas de Transportes – Exemplo 1 Solução Admissível: Mês Cecília Rocha # 8 Distribuição da Produção Uma empresa constrói aviões comerciais para diversas companhias de aviação. A última etapa de produção consiste na construção dos motores a jacto e sua instalação (operação muito rápida) na fuselagem já totalmente montada. Esta empresa tem estado a trabalhar para cumprir alguns contratos já adjudicados que envolvem o fornecimento de um grande número de aviões, pelo que a produção de motores a jacto terá de ser equacionada para os próximos 4 meses, de forma a minimizar o custo de construção e armazenamento. Para satisfazer os prazos de entrega estabelecidos nos contratos, a empresa deverá ter construídos no final do 1º, 2º, 3º e 4º meses, pelo menos, 10, 25, 50 e 70 motores a jacto. As instalações onde serão construídos os motores variam conforme seja necessária a sua afectação para outras actividades, pelo que o número de motores construídos será bastante distinto. A capacidade máxima de construção mensal está indicada no quadro seguinte. Dado que existe variação no custo de construção dos motores, pode compensar que alguns motores sejam construídos nos meses anteriores à sua colocação na fuselagem, o que implicará custo de armazenamento. No entanto, é uma opção que deverá ser analisada. Os custos envolvidos são indicados no quadro seguinte. Produção Mínima Capacidade Máxima Custo de Construção/unidade Custo de Armazenamento/unidade 1 10 (10) 25 1.08 0.015 2 15 (25) 35 1.11 0.015 3 25 (50) 30 1.10 0.015 4 20 (70) 10 1.13 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Preparação da resolução pelo Método dos Transportes Origem Destino Custo associado à construção do motor e seu eventual armazenamento Não se pode instalar um motor que ainda não foi construído cij = M, para obrigar a variável correspondente a ser zero. Deverá ser a mais favorável em termos de diminuição de custos. Como não se conhece esse valor vamos assumir a capacidade máxima de construção de motores em cada mês. Procura, dj Cecília Rocha # 9 Se i j, Se i > j, Oferta, si Número de motores a construir no mês i e a instalar no mês j Parâmetros, cij Colocação de motores na fuselagem nos meses j = 1, 2, 3 e 4 Variáveis de decisão, xij Produção de motores nos meses i = 1, 2, 3 e 4 Corresponde aos valores incluídos nos contratos assinados para fornecimento de aviões. No entanto, ao analisarmos a propriedade da existência de solução admissível verifica-se que o número de motores a produzir (100, no máximo) é superior ao número de motores necessários (70). Assim, será necessário adicionar um Destino Fictício (Df) cuja procura será igual ao “excesso” de oferta (30). O custo de construção e armazenamento de motores associado a este destino fictício tem de ser nulo (não existe !...) 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Preparação da resolução pelo Método dos Transportes Assim, o quadro de custos e requisitos deste problema será: Destinos D1 Origens D3 S1 1.080 S2 M 1.110 S3 M M 1.100 S4 M M 10 15 Procura Cecília Rocha # 10 D2 D4 1.08+1*0.015= 1.095 1.08+2*0.015=1.110 1.08+3*0.015=1.125 DF Oferta 0 25 0 35 1.10+1*0.015=1.115 0 30 M 1.130 0 10 25 20 (30) 1.11+1*0.015=1.125 1.11+2*0.015=1.140 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Propriedades das Soluções nos Problemas de Transportes – Exemplo 2 Solução Admissível: Distribuição de Recursos Uma empresa é responsável pela distribuição de água a 4 concelhos (Berdoo, Los Devils, San Go e Hollyglass), numa região muito árida, pelo que é necessário comprar e transportar a água de outras regiões para abastecer estes concelhos. As fontes de fornecimento de água possíveis têm captação em 3 linhas de água – Rio Colombo, Rio Sacron e Rio Calorie. Esta empresa compra a água a esses fornecedores e revende-a à população dos 4 concelhos. É possível abastecer todos os concelhos a partir de qualquer dos rios, excepto a localidade de Hollyglass que não pode ser abastecida pelo rio Calorie. Consoante o rio que abasteça estes concelhos assim variará o custo da água, devido ao percurso a percorrer e aos acidentes geográficos a ultrapassar. Independentemente desta situação o custo da água no cliente será sempre o mesmo. A administração da empresa enfrenta agora o problema de garantir o fornecimento mínimo necessário aos 4 concelhos e distribuir toda a água disponível proveniente dos 3 rios da forma mais económica possível. Os custos e necessidades envolvidos são indicados no quadro seguinte. Custo por distância entre os rios e os concelhos Cecília Rocha # 11 Berdoo Los Devils San Go Hollyglass Oferta de água disponível Rio Colombo 16 13 22 17 50 Rio Sacron 14 13 19 15 60 Rio Calorie 19 20 23 Ñ pode haver fornecimento 50 Necessidades mínimas 30 70 0 10 Fornecimento pretendido 50 70 30 Toda a que puder ser fornecida 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Preparação da resolução pelo Método dos Transportes Origem, i Destino, j É a máxima possível disponibilizada pelos rios. Procura, dj Cecília Rocha # 12 Custo associado ao transporte e distribuição de água entre os 3 rios e os 4 concelhos Dado que não existe possibilidade de fornecimento de água a Hollyglass a partir do rio Calorie, este parâmetro será M Oferta, si Quantidade de água a fornecer do rio i e a distribuir para o concelho j Parâmetros, cij Distribuição de água para os 4 concelhos Berdoo. Los Devils, San Go e Hollyglass Variáveis de decisão, xij Fornecimento de água a partir dos rios Colombo, Sacron e Calorie A quantidade de água a receber por cada concelho é uma variável de decisão com limite inferior (necessidades mínimas) e limite superior. Este limite superior será a quantidade pretendida, a não ser que o somatório das quantidades pretendidas exceda a água disponível após satisfação das necessidades mínimas. Neste contexto, o concelho de Hollyglass não poderá dispor de mais de 60 u.a. = (50+60+50) – (30+70+0). Como a procura tem de ser um valor único, temos um problema a resolver. Poderemos começar por supor que não existem necessidades mínimas. Assim, as quantidades pretendidas serão as únicas restrições em relação à quantidade de água a distribuir pelos concelhos. No entanto, teremos de fazer um ajuste ao problema dado que não cumprimos a propriedade da solução admissível, ou seja, precisaremos de uma Origem Fictícia (Of) que deverá garantir 50 u.a. = (50+70+30+60) - (50+60+50). O custo associado a esta Origem Fictícia será zero. 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Procura, dj (cont.) Para entrarmos em consideração com as necessidades mínimas de cada concelho, será preciso garantir que a proveniência da água não seja a Origem Fictícia. Assim: Para o concelho de San Go, não há qualquer problema dado que não tem necessidades mínimas; Para o concelho de Hollyglass, também não haverá nenhum inconveniente, dado que as suas necessidades mínimas de 10 u.a. estão garantidas (a quantidade pedida excede em 10 u.a. a oferta atribuída à Origem Fictícia); No caso do concelho de Los Devils, como a quantidade pretendida é igual às necessidades mínimas, toda a procura terá de ser satisfeita pelas origens reais, o que implica considerar um parâmetro M para esta situação na Origem Fictícia; Quanto ao concelho de Berdoo, como a Origem Fictícia tem capacidade para satisfazer a quantidade pretendida por este concelho, é preciso assegurar que as necessidades mínimas sejam satisfeitas pelos outros 3 rios. Assim, terá de ser feito o desdobramento do concelho de Berdoo em 2 fracções, uma correspondente às necessidades mínimas e outra ao excedente de água pretendido, respectivamente, com procura de 30 e 20. Será ainda indispensável garantir que as necessidades mínimas não sejam satisfeitas pela Origem Fictícia, atribuindo um parâmetro M ao respectivo custo. Destinos Berdoo* Los Devils San Go Hollyglass Rio Colombo 16 16 13 22 17 50 Rio Sacron 14 14 13 19 15 60 Rio Calorie 19 19 20 23 M 50 OF M 0 M 0 0 50 Procura 30 20 70 30 60 Origens Cecília Rocha # 13 Oferta Berdoo 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Método Simplex aplicado ao Problema de Transportes Como os Problemas de Transportes são um dos tipos de problemas de programação linear, por isso, é possível resolvê-los pelo método Simplex dado nas 2 aulas anteriores. No entanto, dada a especificidade destes problemas, o método simplex pode ser simplificado – Método Simplex dos Transportes. Preparação do Método Após construir o quadro dos coeficientes das restrições para o método simplex, converter a função objectivo para a forma de maximização e introduzir as variáveis artificiais z1, z2, ..., zm+n, obter-se-á o seguinte quadro simplex: Variáveis Equação Básicas Z (0) Coeficientes Z ... xij ... zi -1 cij M 0 1 1 0 1 ... zm+j M ... Lado Direito 0 (1) ... zi (i) si ... zm+j (m+j) 1 dj ... (m+n) Cecília Rocha # 14 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Preparação do Método Falta agora realizar algumas operações algébricas antes da 1ª iteração para eliminar os coeficientes das variáveis (artificiais) básicas iniciais da linha (0) que sejam diferentes de zero. Após essas operações, a nova linha (0) terá a seguinte forma: Variáveis Básicas Z Equação (0) Coeficientes Z -1 ... xij Cij – ui - vj ... zi M - ui ... zm+j M - vj ... Lado Direito m n i 1 j 1 si ui dj v j Onde: ui – múltiplo da linha original (i) que tem de ser subtraído (directa ou indirectamente) à linha original (0) no método simplex, durante todas as operações que levam ao quadro actual vj – múltiplo da linha original (m+j) que tem de ser subtraído (directa ou indirectamente) à linha original (0) no método simplex, durante todas as operações que levam ao quadro actual Se xij é uma variável não básica, então cij – ui – vj é interpretada como a taxa a que Z se irá alterar à medida que xij aumenta Cecília Rocha # 15 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Preparação do Método Em primeiro lugar, não são necessárias variáveis artificiais porque se pode obter uma solução básica inicial com métodos auxiliares simples A linha (0) pode ser obtida sem utilizar qualquer outra linha, calculando os valores actuais de ui e vj directamente. Dado que cada variável básica tem de ter coeficiente zero na linha (0), os valores de ui e vj podem ser obtidos pela resolução de um conjunto de equações: cij – ui – vj = 0 Cecília Rocha # 16 para cada i e j em que xij é variável básica A variável básica de saída pode ser identificada facilmente sem utilizar os coeficientes das variáveis básicas de entrada, assim como, a nova SBA pode ser detectada imediatamente sem se realizarem nenhumas operações algébricas. Deste modo podemos prescindir de quase todo o quadro do método simplex. Além dos dados de base (parâmetros cij, oferta si e procura dj), o método simplex dos transportes só precisa da SBA inicial, dos valores actuais de ui e vj e dos valores resultantes da operação cij – ui – vj para as variáveis não básicas xij. Estes dados podem ser organizados num quadro denominado – Quadro Simplex dos Transportes. 2001/2002 INVESTIGAÇÃO OPERACIONAL 7ª Aula (cont.) Preparação do Método Formato do quadro simplex dos transportes Iteração ? 1 2 1 2 Destino ... ... n Origem c11 c12 c1n c21 c22 c2n cm1 cm2 cmn Oferta s1 s2 ... m Procura d1 d2 sm dn Vj = cij - ui cij xij Cecília Rocha # 17 ui = cij - vj Variável Básica cij cij – ui - vj Z= Variável Não Básica 2001/2002