Modelagem e simulação de sistema logístico de distribuição de carne de frango

O sistema logístico para distribuição de produtos acabados caracteriza-se pela integração dos serviços de comunicação, transporte e financeiros com a finalidade de atender às demandas do consumidor final. Estima-se que no estado do Espírito Santo, o consumo de carne de frango seja de 44,4 quilos per capita por ano. Para atender a esta demanda, o estado conta com matadouros-frigoríficos distribuídos pelo seu território, bem como, com a participação de outras empresas localizadas no país. Em sistemas de transportes, são característicos Problemas de Roteamento de Veículos (VRP), que precisam ser estudados, caracterizados e otimizados, normalmente, através de rotinas computacionais, que permitem avaliar maior quantidade de variáveis. O presente trabalho teve por objetivo caracterizar um VRP de um matadouro-frigorífico da região do Sul do Espírito Santo e desenvolver um aplicativo computacional que seja suporte para os gestores de logística, servindo para avaliar e propor rotas, e analisar parâmetros logísticos do processo de distribuição de produtos. No desenvolvimento do aplicativo computacional foi necessário caracterizar o sistema logístico da empresa, coletar e analisar os dados das operações logísticas, desenvolver as rotinas computacionais que representassem o sistema em estudo, verificar a confiabilidade dos resultados fornecidos pelo aplicativo, validá-lo e então, poder realizar as experimentações. O aplicativo desenvolvido permitiu reproduzir dados do sistema estudado e avaliar rotas segundo parâmetros logísticos. Pode-se concluir que o aplicativo computacional desenvolvido é útil aos gestores de logística, permitindo a avaliação das rotas praticadas e de novas configurações de rotas.

Modelagem e simulação de um secador intermitente de fluxos contracorrentes para frutos do cafeeiro

Este trabalho foi realizado com o objetivo de desenvolver um modelo computacional para simular a secagem de frutos café em um secador intermitente de fluxos contracorrente, empregando a linguagem de simulação EXTEND™ e o Modelo de Thompson (THOMPSON; PEART; FOSTER, 1968). Para validação do modelo desenvolvido foram utilizados dados experimentais obtidos por Silva (1991), em que foram empregados três níveis de temperatura do ar de secagem de 60, 80 e 100 °C. O modelo desenvolvido foi validado, sendo constatados desvios absolutos de 1,8% b.u e 1,1 kg e erros relativos de 11% e 1,6% na previsão dos parâmetros teor de água final e consumo de lenha, respectivamente. O modelo validado foi empregado na condução de experimentos tipo comparação de cenários. O primeiro experimento refere a alterações do ciclo operacional em que foram alterados os tempos de movimentação e de parada do fluxo da massa de grãos. E o segundo refere à alteração da configuração do secador quanto às alturas das câmaras de secagem e descanso. O ciclo operacional com os tempos de movimentação de um minuto e de parada de dezesseis minutos, para a temperatura do ar de secagem de 100 °C, proporcionou o melhor desempenho, sendo constatado tempo secagem de 12,3 h, consumo de lenha de 109,5 kg, consumo específico de energia de 7660 kJ.kg-1 de água evaporada, e capacidade de secagem de 87,86 kg.h-1. Quanto à configuração do secador, o melhor desempenho ocorreu para altura da câmara de secagem de 2,3 m usando a temperatura do ar de secagem de 100 °C, em que foram simulados tempo de secagem de 12,0 h, consumo de lenha de 106,5 kg, consumo específico de energia, de 7433 kJ.kg-1 de água evaporada, e capacidade de secagem de 90 kg.h-1. Desse modo, na condução da secagem de frutos de café em um secador intermitente de fluxos contracorrentes é recomendado o ciclo operacional com tempos de movimentação de um minuto e o de parada de dezesseis minutos, e não empregar a câmara de descanso. Essa conclusão está fundamentada em índices de desempenho do secador. Ressalta-se que não foram simulados os impactos nos parâmetros de qualidade.

Modelagem Sísmica via métodos das diferenças finitas: caso da bacia do Amazonas

Este trabalho tem por objetivo apresentar os resultados da modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações e múltiplas reflexões, tendo a Bacia do Amazonas como referência à modelagem. As condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas foram analisadas numericamente pelo método das diferenças finitas, visando melhor compreensão e controle da interpretação de dados sísmicos. A geologia da Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o Recente que atingem espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançam espessuras de centenas de metros e perfazem um volume total de aproximadamente 90000 Km³. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica o que impossibilita melhor interpretação dos horizontes refletores que se encontram abaixo destas soleiras. Para representar situações geológicas desse tipo foram usados um modelo (sintético) acústico de velocidades e um código computacional elaborado via método das diferenças finitas com aproximação de quarta ordem no espaço e no tempo da equação da onda. A aplicação dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas melhorou a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observada a influência significativa das reflexões múltiplas devido à camada de alta velocidade, isto provocou maior perda de energia e dificultou a interpretação dos alvos. Por esta razão recomenda-se a integração de dados de superfície com os de poço, com o objetivo de obter melhor imagem dos alvos abaixo das soleiras de diabásio.

Estimativa do escoamento superficial em uma bacia hidrográfica com base em modelagem dinâmica e distribuída

A estimativa do escoamento superficial em bacias hidrográficas é de suma importância para conservação dos recursos naturais; entretanto, esse é um processo complexo e dinâmico, principalmente no contexto de sua variabilidade espacial. Dessa forma, torna-se adequada a aplicação dos Sistemas de Informações Geográficas (SIG) usando pequenas células de informação, pois assim é possível considerar o comportamento espacial das variáveis associadas à origem do escoamento superficial. Este trabalho teve como objetivo implementar os modelos hidrológicos Curva Número (CN-SCS) e Curva-Número Modificado (CN-MMS), com base na linguagem de programação do SIG PCRaster e em uma base de dados reduzida, de forma distribuída e dinâmica, com o intuito de estimar as lâminas de escoamento superficial geradas numa bacia hidrográfica de Latossolos, localizada no município de Nazareno, região dos Campos das Vertentes, Minas Gerais. Para aplicação do modelo CN-SCS foi preciso desenvolver um mapa com valores de CN no formato do PCRaster, enquanto para o modelo CN-MMS foram necessários os seguintes mapas: umidade volumétrica de saturação do solo, umidade volumétrica inicial do solo e profundidade de solo. Para simulação e avaliação de ambos os modelos, foram aplicados 18 eventos de chuva natural que provocaram escoamento superficial, durante o ano hidrológico 2004-2005, e suas respectivas lâminas de escoamento observadas. A análise do desempenho dos modelos foi feita aplicando-se análise de sensibilidade baseada no erro médio e na Raiz do Erro Quadrático (REQ). Tendo-se como referência essas estatísticas de precisão, pôde-se constatar que o modelo CN-MMS apresentou melhor calibração quando comparado ao modelo CN-SCS, devido à consideração direta da umidade inicial do solo. Contudo, a estruturação dos modelos no SIG PCRaster possibilitou o desenvolvimento de uma ferramenta computacional eficaz e útil para simulação do escoamento superficial, visto que propicia estruturação de rotinas computacionais considerando os problemas associados à variabilidade espacial dos dados de entrada dos modelos.

Modelagem da distribuição diamétrica em povoamentos de eucalipto submetidos a desbaste utilizando autômatos celulares

Este estudo apresenta um modelo de distribuição diamétrica baseado em um modelo de Autômatos celulares (AC) unidimensionais e redes neurais artificiais (RNA) para a simulação de desbaste. Cada célula do AC proposto representa uma classe de dap, sendo o estado futuro previsto em função do estado atual dessa célula, do estado de suas quatro células vizinhas e de sua idade atual e futura. Como regra de evolução, utilizou-se uma RNA. A exatidão foi avaliada empregando-se o procedimento estatístico L relação entre frequências observada e estimada e realismo biológico do modelo construído. Entre as redes treinadas, foram selecionadas as 10 que representavam a evolução da distribuição diamétrica com maior exatidão. Entre essas 10 RNA, sete apresentaram valores estimados estatisticamente iguais aos observados (p>0,01). O enfoque de modelagem proposto permite estimar distribuições diamétricas futuras com exatidão.

Rotina computacional para a determinação da velocidade de sedimentação das partículas do solo em suspensão no escoamento superficial

O presente trabalho teve como objetivo desenvolver uma rotina computacional para a determinação da velocidade de deposição de partículas em suspensão no escoamento superficial, verificar sua aplicação por intermédio de modelo de transporte de sedimentos e comparar os resultados obtidos com dados experimentais. Empregou-se na rotina o processo iterativo de Newton-Rapshon para a solução das equações empregadas na determinação da velocidade de deposição de partículas em suspensão no escoamento superficial, e na solução da equação do transporte de sedimentos empregou-se a técnica das diferenças finitas. Essas rotinas foram empregadas na implementação do modelo MTSES (Modelo para Transporte de Solutos no Solo e no Escoamento Superficial). As velocidades de queda das partículas obtidas pela rotina desenvolvida, em média, foram superestimadas, com erro relativo médio de 0,63%, o que possibilitou a utilização da rotina desenvolvida no MTSES. O modelo MTSES superestimou o total de sedimentos transportados pelo escoamento superficial para todas as intensidades de precipitação empregadas neste estudo, com variações porcentuais de 15,6 a 58,3%.

Modelagem de máquina pneumática recolhedora de frutos de café em terreiro usando análise por elementos finitos

Teve-se o objetivo de desenvolver um modelo matemático por meio de análise de elementos finitos, utilizando o programa computacional ANSYS®, versão 5.7, para otimizar o projeto de máquina recolhedora de frutos de café no terreiro. A modelagem da máquina foi realizada com base no levantamento das características aerodinâmicas dos frutos de café e da vazão de ar necessária para o transporte pneumático dos frutos. Foram obtidas, experimentalmente, as pressões estáticas nos dutos da máquina, sendo esses valores comparados com os resultados determinados pelo programa ANSYS, no intuito de validar o modelo. Com base nos resultados numéricos obtidos, concluiu-se que a modelagem desenvolvida apresentou resultados próximos aos determinados experimentalmente, obtendo erro relativo médio nos valores simulados de pressão de 9,2%. Por meio da modelagem, identificaram-se faixas de pressão que dificultariam o transporte pneumático dos frutos de café em alguns pontos da máquina. Esses problemas foram corrigidos e, com isso, o fluxo de ar proporcionado pelo ventilador foi suficiente para succionar os frutos de café no terreiro e transportá-los para dentro do reservatório da máquina. A modelagem desenvolvida atendeu às necessidades propostas no trabalho para o recolhimento dos frutos de café utilizando transporte pneumático eficientemente.

Rotina computacional e equação simplificada para modelar o transporte de sedimentos num Latossolo Vermelho Distrófico

Os objetivos do presente trabalho foram desenvolver rotina computacional para a solução da equação de Yalin e do diagrama de Shields e avaliar uma equação simplificada para modelar a capacidade de transporte de sedimento num Latossolo Vermelho Distrófico que possa ser utilizada no Water Erosion Prediction Project - WEPP, assim como em outros modelos de predição da erosão do solo. A capacidade de transporte de sedimento para o fluxo superficial foi representada como função-potência da tensão cisalhante, a qual revelou ser aproximação da equação de Yalin. Essa equação simplificada pôde ser aplicada em resultados experimentais oriundos de topografia complexa. A equação simplificada demonstrou acuracidade em relação à equação de Yalin, quando calibrada utilizando-se da tensão média cisalhante. Testes de validação com dados independentes demonstraram que a equação simplificada foi eficiente para estimar a capacidade de transporte de sedimento.

Metodologia para obtenção do hidrograma de escoamento superficial em encostas e canais. Parte II: modelo computacional e análise de sensibilidade

Estimativas de vazão máxima de escoamento superficial são necessárias para o projeto de obras hidráulicas em bacias urbanas e rurais. A dificuldade em aplicar os procedimentos disponíveis para calcular a variação do escoamento superficial com o tempo e de seu valor máximo deve-se à inexatidão dos métodos usados para esse objetivo e à variabilidade nos resultados que podem ser obtidos por profissionais que usem o mesmo procedimento. Dessa forma, a investigação de um método que produza estimativas confiáveis da vazão máxima e do hidrograma de escoamento superficial é de grande interesse. Neste trabalho, desenvolveu-se e avaliou-se a sensibilidade de um software (HIDROGRAMA 2.1) que permite a obtenção do hidrograma de escoamento superficial, da vazão máxima e seu tempo de ocorrência, da altura e da velocidade máximas do escoamento, do volume e da lâmina de escoamento superficial em encosta e em canais. O modelo apresentou grande sensibilidade ao período de retorno, à taxa de infiltração estável e ao comprimento da encosta e do canal.

Modelagem de um sistema de informação para rastreabilidade na indústria do vinho baseado em uma arquitetura orientada a serviços

Apresenta-se neste trabalho pesquisa científica que levou à modelagem de um sistema de informação para o registro de dados de rastreabilidade destinado à indústria do vinho, segundo os conceitos de uma arquitetura computacional orientada a serviços (service oriented architecture, SOA). Desde 2005, a manutenção de tais registros é obrigatória para todos os produtores que pretendem exportar para países da União Europeia. Ademais, os consumidores finais, inclusive brasileiros, preocupados com a segurança alimentar, apresentam demanda crescente por informações. Buscou-se uma solução que contemple a indústria coletivamente, por meio de consórcios ou associações de produtores, diluindo custos e compartilhando benefícios. A partir do levantamento bibliográfico realizado e de contatos com pesquisadores e representantes do setor produtivo vinícola brasileiro, em Bento Gonçalves -RS, pesquisaram-se tópicos de tecnologia da informação ligados ao tema, além de detalhes sobre o processo produtivo do vinho. A partir de um modelo do processo produtivo do vinho sugerido pelos autores, empregou-se a Unified Modeling Language (UML), para a modelagem do sistema, que utiliza a tecnologia de web services para o tráfego de dados. Com um protótipo funcional, constataram-se a viabilidade do modelo e a possibilidade de sua utilização em outros domínios.

Difusão do cloreto de sódio no processo de salga de queijos: modelagem matemática com o emprego do método de elementos finitos

A transferência de um soluto (cloreto de sódio), através de uma matriz sólida tridimensional (queijo) foi estudada aplicando-se o método de elementos finitos. A formulação variacional (Galerkin) do problema diferencial (modelo de difusão) teve como base teórica a 2ª lei de Fick. Os procedimentos para integração no tempo foram o de Crank-Nicolson e o de Euler-modificado, que foram escolhidos por apresentarem estabilidade incondicional. O programa computacional desenvolvido mostrou-se versátil para resolver situações de amostragem em condições mais realistas e pode ser aplicado para geometrias complexas. O modelo proposto permitiu uma boa estimativa do ganho de sal no queijo, usando um coeficiente de difusão cujo valor pode ser obtido por extrapolação de dados experimentais. A aplicação do método numérico (MEF), com o esquema de Crank-Nicolson, na simulação da difusão do cloreto de sódio na salga de queijos, mostrou boa aproximação quando os resultados foram comparados com os valores experimentais encontrados na literatura especializada.

O uso da simulação computacional como ferramenta de análise : um estudo de caso em empresa de câmaras de bronzeamento

A concorrência nos mercados locais e internacionais tem levado as empresas a tomarem-se cada vez mais competitivas e atentas as mud3nças no mercado, possibilitando-Ihes melhorias em seus processos produtivos. Este trabalho apresenta a modelagem e análise de melhorias no processo produtivo de uma empresa de montagem de câmaras de bronzeamento. A simulação computacional foi utilizada como ferramenta de análise das possíveis oportunidades de melhorias a serem introduzidas no processo produtivo. O processo original foi modelado com o objetivo de identificar e caracterizar a distribuição de tempos e utilização dos equipamentos e rerorsos presentes no processo produtivo da empresa. Diferentes propostas de melhorias foram investigadas e uma nova modelagem do processo foi sugerida, considerando exigências do mercado e características do processo produtivo empresarial. Comparações matemáticas entre os resultados foram elaboradas e as conclusões sobre o estudo são apresentadas ao final do trabalho. Os resultados da simulação estão sendo usados para avaliar possíveis alterações no processo produtivo da empresa, bem como para gerar documentos que registram conhecimento corporativo sobre possíveis arranjos na produção.

Modelagem mecânica e aproximação por métodos estabilizados de escoamentos multicomponentes

Este estudo foi motivado pela possibilidade de se empregar os conhecimentos da engenharia mecânica na solução de problemas de engenharia de alimentos por métodos numéricos, assim como pela utilização da dinâmica dos fluidos computacional (CFD) em mais um campo de pesquisa. A idéia básica foi a aplicação do método de elementos finitos na solução de problemas de escoamentos envolvendo mistura de diferentes componentes. Muitos alimentos apresentam-se como fluidos, e seu comportamento material pode ser newtoniano ou não newtoniano, às vezes descrito por relações constitutivas bastante complexas. Utilizou-se uma teoria de misturas apoiada nos conceitos de mecânica do contínuo para a modelagem mecânica do que se passou a considerar como um sistema multicomponente. Necessitou-se de uma detalhada revisão sobre os postulados clássicos da mecânica para que se pudesse recolocá-los, com alguma segurança e embasamento teórico, para sistemas multicomponentes. Tendo em mãos a modelagem do balanço de momentum e massa em sistemas multicomponentes, pôde-se aproximar estas equações através do método de elementos finitos. A literatura aponta que o método clássico de Galerkin não possui a eficiência necessária para a solução das equações de escoamento, que envolvem uma formulação mista onde se faz necessário tomar compatíveis os subespaços de velocidade e pressão, e também devido à natureza assimétrica da aceleração advectiva, o que também aparece como uma dificuldade na solução de problemas de advecçãodifusão, nos casos de advecção dominante. Assim, fez-se uso do método estabilizado tipo GLS, o qual supera as dificuldades enftentadas pelo método de Galerkin clássico em altos números de Reynolds, adicionando termos dependentes da malha, construídos de forma a aumentar a estabilidade da formulação de Galerkin original sem prejudicar sua consistência. Os resultados numéricos dividem-se em três categorias: problemas de transferência de quantidade de movimento para fluidos newtonianos, problemas de transferência de quantidade de movimento para fluidos com não linearidade material e problemas de advecção e difusão de massa em misturas. A comparação de algumas aproximações obtidas com as de outros autores se mostraram concordantes. A aproximação de problemas de fluidos segundo os modelos Carreau e Casson geraram os resultados esperados. A aproximação de um problema de injeção axial com mistura de dois fluidos produziu resultados coerentes, motivando a aplicação prática da aproximação por métodos estabilizados de problemas de misturas.

Modelagem de situações de punção em lajes de concreto armado através do método dos elementos finitos

O objetivo deste trabalho é o desenvolvimento de um programa computacional, elaborado através do método dos elementos finitos, para análise tridimensional de situações de punção em lajes de concreto armado com e sem armadura de cisalhamento. São utilizados modelos constitutivos elasto-viscoplásticos para representar o comportamento do concreto. Para a simulação do concreto submetido a um carregamento de longa duração, empregou-se o modelo de camadas superpostas. Este modelo admite que o concreto é constituído por um número conveniente de camadas com a mesma deformação, porém características materiais diferentes. A tensão total é obtida pela soma das contribuições de cada camada. A armadura é introduzida no modelo como uma linha de material mais rígido no interior do elemento de concreto. Considera-se aderência perfeita entre o concreto e o aço. Assim, os deslocamentos ao longo das barras de aço são determinados a partir dos deslocamentos nodais dos elementos de concreto. Desta forma, a matriz de rigidez da armadura tem a mesma dimensão da matriz de rigidez do elemento de concreto, as quais são somadas para a obtenção da matriz de rigidez total Para representar o comportamento após à fissuração do concreto, é utilizado o modelo de fissuras distribuídas, onde apenas as propriedades dos materiais são modificadas e a malha de elementos finitos não necessita de atualização. Sugerem-se duas formulações para a consideração da contribuição do concreto entre fissuras e da redução na capacidade de transferência de corte devido ao aumento da abertura da fissura. Este último é incluído através de uma simplificação, que consiste em atribuir ao módulo de elasticidade transversal, correspondente ao plano fissurado, um valor reduzido. Para comprovar a validade do programa computacional desenvolvido, comparam-se os resultados numéricos com valores experimentais de ensaios realizados por Coelho (1999), Trautwein (2001) e Fusco (1988). Obtiveram-se boas aproximações para lajes com e sem armadura de cisalhamento. Desta forma, o programa computacional desenvolvido pode ser empregado para generalizar resultados experimentais em lajes cogumelo com variação na resistência do concreto, na bitola das barras de armadura, no espaçamento entre as barras e no número de camadas de armadura de cisalhamento.

Modelagem estocástica de sistemas hierárquicos de estoques

Esta dissertação apresenta a Modelagem Estocástica de Sistemas Hierárquicos de Estoques (SHEs) como um instrumento de uso gerencial, orientado para a maximização do nível de serviço oferecido aos consumidores finais e a minimização do nível global de estoques em uma cadeia de abastecimento. Tal modelagem envolve otimização de sistemas não-lineares com integrais indefinidas, normalmente tratados através de simulação, dada a impossibilidade de estabelecer uma solução analítica para o problema. Neste trabalho, diferentemente, a solução é obtida através da otimização do sistema, a partir de uma modelagem que estabelece os níveis de estoque-alvo, os pontos de reposição de estoques, os níveis de serviço oferecidos aos consumidores finais e o nível ótimo global de estoques. A sua operacionalização utiliza um aplicativo computacional, projetado e desenvolvido especificamente para este fim. Estruturas hierárquicas, quando otimizadas localmente, em cada um de seus níveis, podem apresentar resultados subótimos em nível global. Assim, não é suficiente identificar os ótimos locais para compor uma solução otimizadora para o sistema. Para resolver esse problema, a maior parte dos pesquisadores utiliza simulação. A qualidade das soluções assim obtidas geralmente é inferior à qualidade da solução obtida mediante otimização global do sistema A presente modelagem tem como objetivo gerar soluções mais qualificadas para o problema da otimização de SHEs. Nesta dissertação, as principais abordagens para otimização de SHEs são apresentadas e analisadas, com destaque para o trabalho de Rosenbaum (1981), cujos resultados, obtidos através de simulação, são comparados com os resultados deste trabalho.