Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos baseado em volumes de controle em malhas não estruturadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação numérica de escoamentos de fluídos pelo método de elementos finitos baseado em volumes de controle com a técnica de passo fracionado

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação de grandes escalas de escoamentos incompressíveis com transferência de calor e massa por um método de elementos finitos de subdomínio

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Comparação de desempenho entre a formulação direta do Método dos Elementos de Contorno com Funções Radiais e o Método dos Elementos Finitos em problemas de Poisson e Helmholtz

O presente trabalho objetiva avaliar o desempenho do MECID (Método dos Elementos de Contorno com Interpolação Direta) para resolver o termo integral referente à inércia na Equação de Helmholtz e, deste modo, permitir a modelagem do Problema de Autovalor assim como calcular as frequências naturais, comparando-o com os resultados obtidos pelo MEF (Método dos Elementos Finitos), gerado pela Formulação Clássica de Galerkin. Em primeira instância, serão abordados alguns problemas governados pela equação de Poisson, possibilitando iniciar a comparação de desempenho entre os métodos numéricos aqui abordados. Os problemas resolvidos se aplicam em diferentes e importantes áreas da engenharia, como na transmissão de calor, no eletromagnetismo e em problemas elásticos particulares. Em termos numéricos, sabe-se das dificuldades existentes na aproximação precisa de distribuições mais complexas de cargas, fontes ou sorvedouros no interior do domínio para qualquer técnica de contorno. No entanto, este trabalho mostra que, apesar de tais dificuldades, o desempenho do Método dos Elementos de Contorno é superior, tanto no cálculo da variável básica, quanto na sua derivada. Para tanto, são resolvidos problemas bidimensionais referentes a membranas elásticas, esforços em barras devido ao peso próprio e problemas de determinação de frequências naturais em problemas acústicos em domínios fechados, dentre outros apresentados, utilizando malhas com diferentes graus de refinamento, além de elementos lineares com funções de bases radiais para o MECID e funções base de interpolação polinomial de grau (um) para o MEF. São geradas curvas de desempenho através do cálculo do erro médio percentual para cada malha, demonstrando a convergência e a precisão de cada método. Os resultados também são comparados com as soluções analíticas, quando disponíveis, para cada exemplo resolvido neste trabalho.

Cálculo automático de estruturas: análise estrutural de lajes através do método dos elementos finitos

O principal objectivo deste trabalho foi o desenvolvimento de um programa de cálculo automático de lajes baseado no método dos elementos finitos. Como introdução ao tema é efectuada uma análise da evolução das estruturas ao longo dos tempos, diferenciando os vários tipos de estruturas que existem. São explicitados os conceitos fundamentais da mecânica dos sólidos, bem como os vários tipos de análises estruturais, os tipos de elementos finitos mais utilizados, bem como a classificação dos diversos tipos de estruturas. Devido à importância que o método dos elementos finitos tem na engenharia actual, é apresentada uma breve descrição da sua evolução ao longo dos tempos. Como abordagem mais abrangente são apresentados os conceitos gerais do método, evoluindo depois para uma descrição mais pormenorizada, aplicado aos elementos de laje adoptados no programa de cálculo desenvolvido. As lajes podem ser modeladas através da teoria das lajes finas ou espessas, sendo apresentados os fundamentos teóricos e simplificações que sustentam estas duas abordagens. A validação do programa de cálculo elaborado nesta dissertação, é efectuado através de dois exemplos, um modelo simples e outro mais complexo onde se demonstram todas as potencialidades do programa. A validação é efectuada através da comparação dos resultados, obtidos pelo programa e por um programa de referência no cálculo estrutural, o SAP2000. O desenvolvimento deste trabalho tem um objectivo mais abrangente de no futuro, este módulo, poder ser incluído num pacote de cálculo estrutural alargado a outro tipo de estruturas. Com este pressuposto, foi efectuada uma descrição pormenorizada da organização do programa e das suas capacidades e desenvolvido um manual de utilização.

Simulação numérica do comportamento de lajes: aplicação do método dos elementos finitos a problemas tridimensionais de elasticidade linear usando o freeFEM++

O presente trabalho entende a dedução das Equações de Compatibilidade, Equilíbrio, bem como as Leis Constitutivas, de modo a resolver numericamente problemas de Elasticidade Linear formulados para Lajes de Betão. A solução numérica é obtida usando o Método dos Elementos Finitos, implementado em freeFEM++, e comparada com os resultados conseguidos no SAP2000 e com os dados tabelados nas Tabelas de Barès, validando-se desta forma os modelos implementados através do freeFEM++.

Cálculo automático de estruturas: análise de placas através do método dos elementos finitos

Nesta dissertação desenvolveu-se um programa de cálculo automático que efectua a análise estrutural de placas, sujeitas a acções estáticas e considerando um comportamento isotrópico elástico linear dos materiais. Este programa foi elaborado em linguagem Fortran. O método que se utilizou para o cálculo das estruturas é baseado na teoria dos elementos finitos, tendo sido adoptados elementos finitos com dois graus de liberdade por nó (duas translações), associados ao método dos deslocamentos para a sua resolução. Para permitir a representação dos resultados obtidos pelo programa, utilizou-se um módulo gráfico comercial: GiD. A validação do programa desenvolvido foi feita a partir da comparação dos resultados calculados de alguns exemplos com os obtidos através de um programa de cálculo estrutural de referência: SAP2000. Esta dissertação contempla ainda um manual, que tem como finalidade o fornecimento de indicações para a correcta utilização do software desenvolvido, o qual inclui orientações para a introdução de dados e respectiva interpretação de resultados.

Cálculo automático de estruturas: análise estrutural de pórticos planos utilizando o método dos elementos finitos

O principal objectivo da presente dissertação foi o desenvolvimento de um programa de cálculo automático de estruturas reticulares planas. A execução deste trabalho insere-se num programa mais vasto e ambicioso, em que se pretende desenvolver um pacote de cálculo automático de estruturas para utilização académica. Este pacote incluirá o tipo de elementos mais utilizado na Engenharia Civil, podendo efectuar análises lineares e não-lineares, considerando acções estáticas e dinâmicas. Pretende-se também dotar este pacote com o dimensionamento e verificação da segurança dos elementos estruturais. O programa foi desenvolvido recorrendo ao Método dos Elementos Finitos. Em face da aplicação deste método é efectuada uma breve resenha histórica, incluindo o seu desenvolvimento, campos de aplicação e importância na Engenharia Civil para o cálculo de estruturas. Foram implementados elementos finitos do tipo barra 2D de secção constante, em que se admitiu um material isotrópico com comportamento elástico linear. O programa permite a aplicação de cargas estáticas nos nós e elementos de barra, assentamentos de apoio, e fornece os resultados de forma gráfica e numérica. A validação do programa de cálculo foi efectuada através da comparação dos resultados obtidos pelo programa com um outro Programa de cálculo estrutural comercial (SAP2000) e com valores apresentados na bibliografia [LNEC, 1984]. Com o objectivo de incluir este trabalho num pacote de programas mais abrangentes, foi elaborado um fluxograma onde são descritos os blocos principais que compõem o programa desenvolvido e um manual de utilização devidamente pormenorizado.

Análise da estrutura de um edifício em betão armado usando um método dos elementos finitos: aplicação do método dos elementos finitos a problemas tridimensionais de elasticidade linear usando o freeFEM++

Pretende-se, utilizando o Modelo da Elasticidade Linear em freeFEM++, determinar os esforços e deslocamentos de um edifício alto submetido apenas à acção do peso próprio da estrutura e, efectuar estudos comparativos dos resultados obtidos com o SAP2000. O trabalho inicia-se com a introdução da teoria da elasticidade linear, onde são feitas as deduções das Equações de Compatibilidade, Equilíbrio e as Leis Constitutivas, de modo a resolver numericamente o problema de Elasticidade Linear mencionado. O método de elementos finitos será implementado em freeFEM++ com auxílio do GMSH que é uma poderosa ferramenta com capacidade de gerar automaticamente malhas tridimensionais de elementos finitos e com relativa facilidade de pré e pós-processamento.

Modelização de um dispositivo de travagem automático de via, através do método dos elementos finitos

Apresenta-se a modelização, através do Método dos Elementos Finitos, de um sistema de proteção automática de comboios, denominado DTAV – Dispositivo de Travagem Automático de Via. A modelização é realizada com o auxílio do software Ansys 13.0. Descreve-se a evolução histórica dos sistemas de proteção automática de comboios, desde a época vitoriana até à atualidade. São categorizados os vários tipos de sistemas de proteção existentes, consoante as suas funcionalidades e os seus princípios físicos de funcionamento. É apresentada uma breve descrição da normalização em vigor, aplicada aos sistemas de proteção automáticos de comboios. Descreve-se em pormenor o sistema DTAV, nomeadamente a sua funcionalidade, o seu princípio físico de funcionamento e o conjunto de equipamentos de que é constituído. Apresenta-se uma breve introdução sobre o método dos elementos finitos, enquadrando o modelo criado com os princípios físicos em que se fundamenta, nomeadamente através da descrição das equações de eletromagnetismo. Desenvolve-se modelo do sistema DTAV, evidenciado as etapas da construção e os resultados obtidos, validados por um conjunto de especificações e ensaios prévios realizados em laboratório.

Desenvolvimento de programa em linguagem Matlab® para cálculo pelo método de elementos finitos

Para o projeto de qualquer estrutura existente (edifícios, pontes, veículos, máquinas, etc.) é necessário conhecer as condições de carga, geometria e comportamento de todas as suas partes, assim como respeitar as normativas em vigor nos países nos quais a estrutura será aplicada. A primeira parte de qualquer projeto nesta área passa pela fase da análise estrutural, onde são calculadas todas as interações e efeitos de cargas sobre as estruturas físicas e os seus componentes de maneira a verificar a aptidão da estrutura para o seu uso. Inicialmente parte-se de uma estrutura de geometria simplificada, pondo de parte os elementos físicos irrelevantes (elementos de fixação, revestimentos, etc.) de maneira a simplificar o cálculo de estruturas complexas e, em função dos resultados obtidos da análise estrutural, melhorar a estrutura se necessário. A análise por elementos finitos é a ferramenta principal durante esta primeira fase do projeto. E atualmente, devido às exigências do mercado, é imprescindível o suporte computorizado de maneira a agilizar esta fase do projeto. Existe para esta finalidade uma vasta gama de programas que permitem realizar tarefas que passam pelo desenho de estruturas, análise estática de cargas, análise dinâmica e vibrações, visualização do comportamento físico (deformações) em tempo real, que permitem a otimização da estrutura em análise. Porém, estes programas demostram uma certa complexidade durante a introdução dos parâmetros, levando muitas vezes a resultados errados. Assim sendo, é essencial para o projetista ter uma ferramenta fiável e simples de usar que possa ser usada para fins de projeto de estruturas e otimização. Sobre esta base nasce este projeto tese onde se elaborou um programa com interface gráfica no ambiente Matlab® para a análise de estruturas por elementos finitos, com elementos do tipo Barra e Viga, quer em 2D ou 3D. Este programa permite definir a estrutura por meio de coordenadas, introdução de forma rápida e clara, propriedades mecânicas dos elementos, condições fronteira e cargas a aplicar. Como resultados devolve ao utilizador as reações, deformações e distribuição de tensões nos elementos quer em forma tabular quer em representação gráfica sobre a estrutura em análise. Existe ainda a possibilidade de importação de dados e exportação dos resultados em ficheiros XLS e XLSX, de maneira a facilitar a gestão de informação. Foram realizados diferentes testes e análises de estruturas de forma a validar os resultados do programa e a sua integridade. Os resultados foram todos satisfatórios e convergem para os resultados de outros programas, publicados em livros, e para cálculo a mão feitos pelo autor.

Projecto de equipamento de rotação de moldes de grande porte com recurso ao método de elementos finitos

A fabricação de moldes é uma actividade que se desenvolve em Portugal há várias décadas, com qualidade reconhecida mundialmente. A indústria automóvel é um dos principais destinos de grande parte dos moldes produzidos em Portugal. Algumas das peças produzidas nesses moldes possuem dimensões apreciáveis, tendo em conta o peso e dimensão média das peças injectadas em plástico. São exemplo disso os pára-choques dos automóveis e o tabliê. A produção dos moldes para estas peças é complexa, implicando um longo tempo de maquinagem e a manipulação dos moldes em diferentes fases, com vista ao acesso a todas as faces do molde. A IGM – Indústria Global de Moldes, SA. é a empresa responsável pela produção dos moldes para peças de média dimensão dentro do Grupo SIMOLDES. Atendendo à necessidade permanente de rodar a 90o moldes de grande porte, que podem apresentar pesos na ordem das 10 a 30 toneladas, e não existindo no mercado qualquer solução adaptável a esta necessidade, a empresa entendeu por bem levar a efeito o projecto desse equipamento, atendendo ao compromisso custo-benefício que torne viável a realização prática do mesmo. Após os esboços iniciais e uma discussão interactiva e iterativa com a empresa, foram analisadas as diferentes soluções entendidas como viáveis, tendo sido escolhido um dos anteprojectos realizados. Foram ainda discutidas as diversas alternativas de accionamento. Com base nesse anteprojecto, a estrutura foi optimizada e verificada através do Método de Elementos Finitos, tendo sido elaborado o projecto final, com o grau de detalhe necessário à sua fabricação e implementação na empresa.

Criação de programa informático em ambiente Matlab® para estudo de estruturas planas pelo Método de Elementos Finitos.

Os modelos a ser analisados pelo Método de Elementos Finitos são cada vez mais complexos e, nos tempos que correm, seria impensável realizar tais análises sem um apoio computorizado. Existe para esta finalidade uma vasta gama de programas que permitem realizar tarefas que passam pelo desenho de estruturas, análise estática de cargas, análise dinâmica e vibrações, visualização do comportamento físico (deformações) em tempo real, que permitem a otimização da estrutura. Sob o pretexto de permitir a qualquer utilizador uma análise de estruturas simples com o Método dos Elementos Finitos, surge esta tese, onde se irá criar de raiz um programa com interface gráfica no ambiente MATLAB® para análise de estruturas simples com dois tipos de elemento finito, triangular de deformação constante e quadrangular de deformação linear. O software desenvolvido, verificado por comparação com um software comercial dedicado para o efeito, efetua malhagem com elementos bidimensionais triangulares e quadriláteros e resolve modelos arbitrados pelo Método de Elementos Finitos, representando estes resultados visualmente e em formato tabular.

Determinação das tensões residuais e deformações resultantes do processo de soldadura TIG através do método dos elementos finitos

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica