Estudo de um modelo convectivo-difusivo-reativo em combustão no método de elementos finitos

Neste trabalho, apresenta-se um estudo numérico de um modelo convectivo-difusivo-reativo em combustão baseado no Método de Elementos Finitos. Primeiramente, apresenta-se o desenvolvimento das equações de balanço (quantidade de movimento, massa, espécie e energia) que modelam um processo de mistura molecular e reação química, irreversível, de passo único e exotérmica entre duas espécies químicas F (Combustível) e O (Oxidante). Tais espécies reagem e formam um produto P, conforme vFF +vOO ! vPP + calor, onde vF , vO e vP são os coeficientes estequiométricos molares. No modelo, considera-se que a reação é de primeira ordem com respeito a cada um dos reagentes e que a taxa de reação específica segue a cinética de Arrhenius. Em seguida, o modelo é estudado numericamente considerando-se um domínio retangular e condições de contorno do tipo Neumann. Tanto a Técnica das Diferenças Finitas como a Técnica de Elementos Finitos são utilizadas na discretização espacial das equações do modelo. Para a integração no tempo, utiliza-se a método de Runge-Kutta simplificado de três estágios. Os diferentes códigos computacionais obtidos, tanto pela Técnica de Diferenças Finitas como de Elementos Finitos, são comparados frente ao problema de interesse. Observa-se que ambas as técnicas apresentam resultados equivalentes. Além disso, os códigos desenvolvidos são robustos (capazes de lidar com vários conjuntos de parâmetros), de baixo custo e precisos. Por fim, apresenta-se uma revisão do trabalho de Zavaleta [48], no qual obtem-se uma estimativa local do erro na aproximação do problema estudado pela Técnica de Elementos Finitos.

Métodos numéricos para ecuaciones diferenciales rígidas. Aplicación a la semidiscretización del método de elementos finitos

295 p.

Análise de guias de ondas pelos métodos vetorial magnético e dos elementos-finitos

Neste trabalho, é apresentada uma formulação apropriada à análise de guias de ondas eletromagnéticos, cobrindo do espectro de microondas até o da óptica. Nas regiões a partir do ultravioleta, os comprimentos de onda são equivalentes às dimensões atômicas e a formulação necessita de uma abordagem quântica, que não é considerada neste estudo. A formulação é fundamentada nos métodos vetorial magnético e dos elementos finitos (MEF), em meios não homogêneos, anisotrópicos e não dissipativos, embora a dissipação possa ser facilmente introduzida na análise. Deu-se preferência à formulação com o campo magnético em vez do elétrico, pelo fato do campo magnético ignorar descontinuidades elétricas. Ele é contínuo em regiões de permeabilidade homogênea, propriedade dos meios dielétricos em geral ( = 0), independente da permissividade dos respectivos meios, conquanto os campos elétricos sejam descontínuos entre regiões de permissividades diferentes.

Análise da dissipação das tensões em dentes humanos restaurados com facetas laminadas de cerâmica, com três tipos de preparos, através do método dos elementos finitos

O objetivo deste trabalho é analisar in vitro a dissipação de tensões em incisivos centrais superiores humanos restaurados com facetas de cerâmica feldspática, através da análise do método dos elementos finitos, considerando cargas funcionais de mastigação e corte dos alimentos, em função de três tipos de preparos utilizados: sem proteção incisal; com proteção incisal em ângulo e com proteção incisal em degrau palatino. Foram utilizadas modelagens bidimensionais de um incisivo central superior e suas estruturas de suporte, simulando três situações: (Primeira modelagem) incisivo central superior com desgaste vestibular (em forma de janela); (Segunda modelagem) incisivo central superior com desgaste vestibular e proteção incisal em plano inclinado; (Terceira modelagem) incisivo central superior com desgaste vestibular, e proteção incisal com degrau palatino. Foi considerada uma carga (P=100N) com uma inclinação de 45 concentrada, simulando a região de contato do incisivo central inferior com o superior durante a mastigação e uma na região de contato topo a topo dos incisivos superior e inferior, simulando o corte dos alimentos. Após a análise dos dados obtidos pela distribuição de tensões, pode-se concluir que quanto à dissipação das tensões em todo o sistema proposto, com a aplicação de carga em 45, não foram observadas mudanças no estado tensional nos três diferentes preparos. Quando foi aplicada carga vertical, simulando o contato de topo, houve variação no estado tensional no sistema do dente com preparo em janela. Nas facetas, com a aplicação de carga em 45, nos preparos em janela e com proteção incisal em plano inclinado o resultado foi semelhante nos valores tensionais enquanto, nas facetas em dentes preparados com proteção incisal com degrau palatino, a distribuição foi mais homogênea tendo valores superiores, mostrando que o abraçamento do dente diminuiu a flexão.

Análisis del comportamiento de un rodamiento de vuelco mediante elementos finitos y técnicas de subestructuras

[ES]estudio del comportamiento de un rodamiento de vuelco situado en la unión buje-pala de un aerogenerador mediante diferentes técnicas de modelización por elementos finitos y su comparación. Los análisis realizados son: estático con contacto Rough, estático con contacto Bonded y mediante técnicas de subestructuras. En el estudio se ha utilizado ANSYS Mechanical APDL.

Distribuição de tensões em próteses parciais fixas posteriores livres de metal com retentores intracoronários: análise em elementos finitos

Para reabilitar a ausência de um elemento dentário posterior, as próteses parciais fixas (PPF) com retentores intracoronários são uma alternativa aos implantes osseointegrados. O objetivo deste estudo foi avaliar a distribuição de tensões nessas próteses com três combinações de materiais: cerâmica de zircônia parcialmente estabilizada por ítria (ZPEI) revestida por cerâmica de fluorapatita (α), cerâmica de dissilicato de lítio (β) ou compósito fibrorreforçado (γ). Na composição α, foram analisadas a presença ou ausência da cerâmica de revestimento na parede cervical das caixas proximais e três variações na área total da seção transversal dos conectores (4 mm de largura x 3,2, 4,2 ou 5,2 mm de altura). Em 8 modelos bidimensionais de elementos finitos, uma carga vertical de 500 N foi aplicada na fossa central do pôntico e as tensões principais máximas (tração) e mínimas (compressão) foram apontadas em MPa. Inicialmente foram avaliados os 6 modelos com PPF de ZPEI e suas variações. Os maiores valores das tensões de tração foram encontrados no terço cervical dos conectores. Quando presente nestas regiões, a cerâmica de revestimento recebeu tensões acima do limite de sua resistência à flexão. Na comparação entre os modelos sem cerâmica de revestimento na parede cervical das caixas proximais, mesmo aquele com conectores de 3,2 x 4 mm, cuja infraestrutura apresentava 2,5 x 3 mm, poderia ser recomendado para uso clínico. Altos valores de tensões de compressão foram registrados entre o terço oclusal e médio dos conectores, correspondente à união entre as cerâmicas, o que poderia ocasionar, devido à flexão, falhas adesivas. Posteriormente, o modelo de ZPEI com a cerâmica de fluorapatita ausente da parede cervical das caixas proximais e área total dos conectores de 4,2 x 4 mm foi comparado aos dois outros materiais com conectores de mesma área. Na PPF de dissilicato de lítio, os valores representaram uma provável violação do limite de sua resistência à flexão. A PPF de compósito fibrorreforçado apresentou tensões bem abaixo do limite de resistência à flexão de sua infraestrutura, mas, como no modelo de ZPEI, tensões compressivas se concentraram com alto valor entre o terço oclusal e médio dos conectores, local de união entre a resina composta e a infraestrutura de fibras. Os resultados mostraram que a cerâmica de dissilicato de lítio e a presença da cerâmica de fluorapatita na parede cervical das caixas proximais deveriam ser contraindicadas para a condição proposta. Parece viável uma área de conectores na infraestrutura de ZPEI com no mínimo 2,5 x 3 mm. A PPF de compósito fibrorreforçado apresenta resistência estrutural para a situação estudada, mas, como também aquelas compostas de ZPEI, aparenta ter como pontos fracos a adesão entre a infraestrutura e o material de cobertura e a própria resistência deste último.

Caracterización a impacto de caucho reciclado mediante elementos finitos

Análisis de caucho reciclado de manera hiperelástica mediante métodos de ajuste de Mínimos Cuadrados con programa MATLAB y Curve fitting mediante ANSYS. Para la parte viscoelástica se usa Algoritmo de Optimicación con MATLAB. Comprobación de resultados y fiabilidad.

Diseño, Modelado y Análisis por Elementos Finitos de un Cuadro de Bicicleta de Carretera.

[ES]El objetivo de este proyecto es llevar a cabo el diseño conceptual, modelado y análisis estructural de un cuadro de bicicleta de geometría de carretera, empleando herramientas de Diseño Asistido por Ordenador y software de Análisis con Elementos Finitos. El modelo inicial será analizado y optimizado de cara a que cumpla los requerimientos resitentes del escenario de cargas asumido. El propósito de este proyecto es completar la primera fase de un hipotético caso de producción.

Simulación por elementos finitos del conformado de una pieza para carrocería de un automóvil.

[ES]El presente proyecto tiene como objetivo optimizar los parámetros de conformado a partir de la geometría de los troqueles de una pieza de un automóvil, con el fin de asegurar los parámetros de calidad que se exigen a este tipo de procesos se cumplan. Para ello, se dispone de un programa de simulación en elementos finitos cuyo nombre es Pam-Stamp 2G, en el cual se va a llevar a cabo el diseño. Además, se dispone de la matriz de la pieza requerida para una mayor facilidad de resolución. El proyecto constará de una descripción detallada del proceso de diseño del conformado durante todo el curso académico, complementado con un diagrama Gantt, así como las decisiones adoptadas durante el trabajo para lograr un desarrollo óptimo del mismo. Por último, contendrá también un presupuesto en el que se detallará los gastos que ha supuesto el desarrollo completo del trabajo.

Avaliação de ligações soldadas de perfis tubulares em estruturas de aço através do método dos elementos finitos.

A necessidade de se obter soluções para os diversos problemas estruturais na engenharia, associado ao amplo crescimento e difusão da construção metálica fez surgir, por volta dos anos 60, os perfis tubulares. As ligações soldadas entre perfis tubulares são as ligações de maior difusão nas construções com este tipo de perfil, sendo necessário portanto, estudos mais detalhados que tornem a utilização de ligações soldadas mais seguras e otimizadas. O Eurocode 3 possui um item específico para as ligações entre perfis tubulares, em especial as ligações soldadas devido a suas particularidades. Infelizmente, a norma brasileira NBR8800 não aborda as ligações envolvendo perfis tubulares. Para a execução da análise numérica no presente trabalho, faz-se necessário uma revisão bibliográfica. Com estes resultados, foram feitas modelagens de ligações tubulares com o método dos elementos finitos de forma a otimizar os modelos a serem utilizados em uma análise paramétrica futura. Desta forma, foram desenvolvidos dois modelos numéricos, um considerando ligação tipo T entre perfis tubulares quadrados e outro para uma ligação tipo K entre perfis circulares. Estes modelos foram caracterizados através de elementos de casca com seis graus de liberdade por nó considerando-se adicionalmente o efeito de membrana. A análise não-linear realizada considerou a não-linearidade do material através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bi-linear e a não-linearidade geométrica foi atribuída através da formulação de Lagrange atualizado. Dentre as principais conclusões obtidas no presente trabalho, pode-se citar que os resultados para as ligações tipo T, o Eurocode 3, fornece resultados que precisam ser observados com cautela. Todavia, para as ligações do tipo K, os resultados numéricos mostraram-se sempre inferiores aos valores através do Eurocode 3, representando um dimensionamento a favor da segurança.

Implementação paralela de um código de elementos finitos em 2D para as Equações de Navier-Stokes para fluidos incompressíveis com transporte de escalares.

O estudo do fluxo de água e do transporte escalar em reservatórios hidrelétricos é importante para a determinação da qualidade da água durante as fases iniciais do enchimento e durante a vida útil do reservatório. Neste contexto, um código de elementos finitos paralelo 2D foi implementado para resolver as equações de Navier-Stokes para fluido incompressível acopladas a transporte escalar, utilizando o modelo de programação de troca de mensagens, a fim de realizar simulações em um ambiente de cluster de computadores. A discretização espacial é baseada no elemento MINI, que satisfaz as condições de Babuska-Brezzi (BB), que permite uma formulação mista estável. Todas as estruturas de dados distribuídos necessárias nas diferentes fases do código, como pré-processamento, solução e pós-processamento, foram implementadas usando a biblioteca PETSc. Os sistemas lineares resultantes foram resolvidos usando o método da projeção discreto com fatoração LU por blocos. Para aumentar o desempenho paralelo na solução dos sistemas lineares, foi empregado o método de condensação estática para resolver a velocidade intermediária nos vértices e no centróide do elemento MINI separadamente. Os resultados de desempenho do método de condensação estática com a abordagem da solução do sistema completo foram comparados. Os testes mostraram que o método de condensação estática apresenta melhor desempenho para grandes problemas, às custas de maior uso de memória. O desempenho de outras partes do código também são apresentados.

Desenvolvimento e implementação de algoritmo para análise de dutos fabricados ou reparados com material compósito carregados com pressão interna utilizando o método dos elementos finitos.

A integridade de dutos ganha importância à medida em que o desenvolvimento da indústria expande a malha para transportar líquidos e gases num contexto de aumento das exigências ambientais. Um importante aliado para manutenção da integridade de dutos são reparos de materiais compósitos. Estes materiais apresentam baixa densidade, capacidade de direcionar resistência de acordo com as direções das solicitações, execução de reparo a frio sem necessidade de interromper produção ou grande maquinário. Este trabalho tem como objetivo desenvolver e implementar um algoritmo de elementos finitos que permita avaliar os esforços e a resistência das paredes de um tubos fabricados ou reparados com laminados de material compósito carregados com pressão interna. Entre as vantagens de desenvolver um programa tem-se: agilidade de avaliação, menor custo com licença, menores exigências computacionais, possibilidade de desenvolver o programa e o melhor entendimento da modelagem dos fenômenos. Utiliza-se como entrada do programa o diâmetro do duto, pressão interna e parâmetros do laminado. A modelagem em elementos finitos é realizada a partir da teoria clássica de laminados. Aplicando o carregamento resultante da pressão interna, determina-se os deslocamentos e são calculadas as tensões e aplicado o critério de falha de Tsai-Hill em cada camada. Estudos experimentais e numéricos encontrados na literatura foram simulados com o programa gerado e os resultados para propriedades do laminado, tensões nos dutos e pressão de ruptura apresentam concordância com os resultados da literatura.O programa ainda tem sua estrutura modificada para encontrar a pressão de falha a partir dos dados do laminado. O programa implementado permite uma avaliação rápida de resistência do reparo e possibilita avaliar rapidamente a resposta a mudanças nos parâmetros de projeto do laminado.

Elementos finitos em fluidos dominados pelo fenômeno de advecção: um método semi-Lagrangeano.

Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção.

Estabilidade de talude com efeito sísmico a partir dos métodos de equilíbrio limite e de elementos finitos.

A Costa Peruana apresenta alta atividade sísmica, tornando imprescindível a execução de análises de estabilidade que considerem os eventos sísmicos. Com o desenvolvimento de novas ferramentas numéricas, a análise dinâmica está se tornando cada vez mais importante e usual na fase de projeto, não se justificando mais a execução de análises estáticas em locais tão vulneráveis a sismos. A presente dissertação tem como objetivo realizar a análise da estabilidade dos Penhascos de Lima, Peru, considerando quatro taludes distintos, em uma região bastante afetada por abalos sísmicos. As análises foram executadas pelos métodos de equilíbrio limite, elementos finitos e pseudo-estático, buscando-se uma comparação entre os diferentes métodos. O trabalho apresenta uma descrição completa dos taludes em estudo, uma avaliação da condição sismológica da região e finalmente a estabilidade dos Penhascos de Lima, fazendo uso dos programas computacionais Slide (método do equilíbrio limite) e Plaxis (método dos elementos finitos) em 2D. Os resultados mostraram que os três métodos adotados forneceram fatores de segurança compatíveis, principalmente quando se considera perfis menos estratificados. Para perfis homogêneos, as diferenças obtidas foram da ordem de 0,5 a 1,0 %. As análises ressaltaram a importância de considerar a condição dinâmica, e mostraram-se bastante sensíveis aos valores de carga sísmica adotado.

Observação e modelação do comportamento dinâmico do viaduto do Corgo: utilização integrada de modelos de identificação modal e modelos de elementos finitos de barra

Trabalho Final de Mestrado elaborado no Laboratório Nacional de Engenharia Civil (LNEC) para a obtenção do grau de Mestre em Engenharia Civil pelo Instituto Superior de Engenharia de Lisboa no âmbito do protocolo de Cooperação entre o ISEL e o LNEC