Implementação do método das características na modelagem de problemas de convecção natural em cavidades cilíndricas

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

Implementação de um gerador tridimensional de malhas de elementos finitos, com aplicações à simulação computacional em odontologia

Pós-graduação em Matemática - IBILCE

Metodologia dos observadores de estado para diagnose de falhas em sistemas contendo elementos finitos de placas de Kirchoff

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Desenvolvimento de um programa paramétrico para simulação do processo de estampagem utilizando elementos finitos

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

Estudo dos critérios de falhas em compósitos PPS/fibras de carbono utilizando elementos finitos

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

Influência da espessura da camada de cimento e da variação da temperatura no comportamento mecânico de fragmentos cerâmicos utilizados em restaurações estéticas

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

Análise do comportamento experimental e numérico de prismas de alvenaria estrutural submetidos a ações verticais utilizando elementos finitos volumétricos

Pós-graduação em Engenharia Civil - FEIS

Caracterização físico-química, morfológica, análise mecânica e de elementos finitos 3D, de diferentes placas e parafusos metálicos e técnicas de fixação interna, empregadas em fraturas de ângulo mandibular

Pós-graduação em Odontologia - FOA

Análise termoelástica experimental e computacional por elementos finitos da parede externa em concreto armado de um forno de cozimento de anodo

Nesta dissertação foi desenvolvido um estudo sobre o comportamento de uma estrutura em concreto armado submetida a elevadas variações térmicas. A estrutura analisada consiste na parede externa de um forno de cozimento de anodo utilizado na produção de Alumínio, e está submetida a variações térmicas de cerca de 125C devido à própria operação do forno. As motivações principais para o estudo foram a ocorrência de grandes deformações e o surgimento de fissuras na estrutura, o que poderia provocar a inutilização do forno. O objetivo do trabalho foi a investigação das causas destas patologias, assim como o estabelecimento de propostas para reforço da estrutura. Para isso, foram realizados estudos experimentais e computacionais do comportamento da estrutura. Inicialmente, a estrutura foi monitorada utilizando-se transdutores de deslocamento e sensores de temperatura (termopares), conectados a um sistema de aquisição de dados para obtenção e armazenamento automática das amostras ao longo do tempo. Em seguida, foram desenvolvidos modelos computacionais em Elementos Finitos com auxilio do programa computacional Algor, para determinação da distribuição de temperatura e as correspondentes tensões e deformações de origem térmica na estrutura. Nestes estudos, foram realizadas análises estacionárias e transientes de condução de calor, seguidas de análises de tensões de origem térmica. Como conclusão do estudo, tem-se que a metodologia proposta para a solução do problema foi bastante satisfatória, solucionando o problema de forma precisa, porém econômica.

Comportamento mecânico de sistemas prótese/implante em região dos elementos finitos e interferometria holográfica

Pós-graduação em Odontologia - FOAR

Modelagem 1D e 2,5D de dados do método CSEM marinho em meios com anisotropia transversal inclinada

Neste trabalho apresentamos a solução do campo eletromagnético gerado por um dipolo elétrico horizontal em meios transversalmente isotrópicos com eixo de simetria vertical (TIV) e com eixo de simetria inclinado (TII). Para modelos unidimensionais, o campo eletromagnético foi obtido por duas metodologias distintas: (1) solução semi-analítica das equações de Maxwell com auxílio de potenciais vetores no caso TIV e (2) em modelos com anisotropia transversal inclinada o campo eletromagnético foi separado em primário e secundário, e então, o campo secundário foi calculado pelo método de elementos finitos no domínio (k_x, k_y, z) da transformada de Fourier. Para estruturas bidimensionais, foi aplicada a mesma metodologia usado nos modelos TII unidimensionais, onde o campo secundário foi calculado pelo método de elementos finitos no domínio (x, k_y, z), da transformada de Fourier, com a utilização de malhas não estruturadas para discretização dos modelos. Estas respostas foram usados para avaliar os efeitos da anisotropia elétrica nos dados CSEM marinho 1D e 2,5D.

Modelagem 2,5D dos campos usados no Método Eletromagnético a Multi-Frequência - EMMF

Esta tese mostra a modelagem 2,5D de dados sintéticos do Método Eletromagnético a Multi-frequência (EMMF). O trabalho é apresentado em duas partes: a primeira apresenta os detalhes dos métodos usados nos cálculos dos campos gerados por uma bobina horizontal de corrente colocada sobre a superfície de modelos bidimensionais; e a segunda, usa os resultados obtidos para simular os dados medidos no método EMMF, que são as partes real e imaginária da componente radial do campo magnético gerado pela bobina. Nesta segunda parte, observamos o comportamento do campo calculado em diversos modelos, incluindo variações nas propriedades físicas e na geometria dos mesmos, com o intuito de verificar a sensibilidade do campo observado com relação às estruturas presentes em uma bacia sedimentar. Com esta modelagem, podemos observar as características dos dados e como as duas partes, real e imaginária, contribuem com informações distintas e complementares. Os resultados mostram que os dados da componente radial do campo magnético apresentam muito boa resolução lateral, mesmo estando a fonte fixa em uma única posição. A capacidade desses dados em distinguir e resolver estruturas alvo será fundamental para o trabalho futuro de inversão, bem como para a construção de seções de resistividade aparente.

Análise da resolução do modelo de Hughes em sondagens magnetotelúricas

Na maioria dos métodos de exploração geofísica, a interpretação é feita assumindo-se um modelo da Terra uniformemente estratificado. Todos os métodos de inversão, inclusive o de dados eletromagnéticos, exigem técnica de modelamento teórico de modo a auxiliar a interpretação. Na literatura os dados são geralmente interpretados em termos de uma estrutura condutiva unidimensional; comumente a Terra é assumida ser horizontalmente uniforme de modo que a condutividade é função somente da profundidade. Neste trabalho uma técnica semi-analítica de modelagem desenvolvida por Hughes (1973) foi usada para modelar a resposta magnética de duas camadas na qual a interface separando as camadas pode ser representada por uma expansão em série de Fourier. A técnica envolve um método de perturbação para encontrar o efeito de um contorno senoidal com pequenas ondulações. Como a perturbação é de primeira ordem a solução obtida é linear, podemos então usar o princípio da superposição e combinar soluções para várias senoides de forma a obter a solução para qualquer dupla camada expandida em série de Fourier. Da comparação com a técnica de elementos finitos, as seguintes conclusões podem ser tiradas: • Para um modelo de dupla camada da Terra, as camadas separadas por uma interface cuja profundidade varia senoidalmente em uma direção, as respostas eletromagnética são muito mais fortes quando a espessura da primeira camada é da ordem do skin depth da onda eletromagnética no meio, e será tanto maior quanto maior for o contraste de condutividade entre as camadas; • Por outro lado, a resistividade aparente para este modelo não é afetada pela mudança na frequência espacial (v) do contorno; • Em caso do uso da solução geral para qualquer dupla camada na Terra cuja interface possa ser desenvolvida em série de Fourier, esta técnica produziu bons resultados quando comparado com a técnica de elementos finitos. A linerização restringe a aplicação da técnica para pequenas estruturas, apesar disso, uma grande quantidade de estruturas pode ser modelada de modo simples e com tempo computacional bastante rápido; • Quando a dimensão da primeira camada possui a mesma ordem de grandeza da estrutura, esta técnica não é recomendada, porque para algumas posições de sondagem, as curvas de resistividade aparente obtidas mostram um pequeno deslocamento quando comparadas com as curvas obtidas por elementos finitos.

Análise teórica do problema de Weaver da falha infinita, modo TE magnetotelúrico

É apresentada uma solução totalmente analítica do modelo da falha infinita para o modo TE magnetotelúrico, levando em conta a presença do ar, utilizando como base o trabalho de Sampaio apresentado em 1985, que apresenta uma solução parcialmente analítica e parcialmente numérica – solução híbrida. Naquela solução foram aplicadas oito condições de contorno, sendo que em quatro delas foram encontradas inconsistências matemáticas que foram dirimidas com alterações adequadas nas soluções propostas por Sampaio. Tais alterações propiciaram a chegarse à solução totalmente analítica aqui apresentada. A solução obtida foi comparada com a solução de Weaver, com a de Sampaio e com o resultado do método numérico dos elementos finitos para contrastes de resistividade iguais a 2, 10 e 50. A comparação da solução analítica, para o campo elétrico normalizado, com a solução de elementos finitos mostra que a solução analítica proporcionou resultados mais próximos, em comparação aos fornecidos por Weaver e por Sampaio. Este é um problema muito difícil, aberto para uma solução analítica definitiva. A solução apresentada aqui é, nesta direção, um grande passo.

Os efeitos Delaware e Groningen: um estudo quantitativo por elementos finitos

Os efeitos Delaware e Groningen são dois tipos de anomalia que afetam ferramentas de eletrodos para perfilagem de resistividade. Ambos os efeitos ocorrem quando há uma camada muito resistiva, como anidrita ou halita, acima do(s) reservatório(s), produzindo um gradiente de resistividade muito similar ao produzido por um contato óleo-água. Os erros de interpretação produzidos têm ocasionado prejuízos consideráveis à indústria de petróleo. A PETROBRÁS, em particular, tem enfrentado problemas ocasionados pelo efeito Groningen sobre perfis obtidos em bacias paleozóicas da região norte do Brasil. Neste trabalho adaptamos, com avanços, uma metodologia desenvolvida por LOVELL (1990), baseada na equação de Helmholtz para H_Φ, para modelagem dos efeitos Delaware e Groningen. Solucionamos esta equação por elementos finitos triangulares e retangulares. O sistema linear gerado pelo método de elementos finitos é resolvido por gradiente bi-conjugado pré-condicionado, sendo este pré-condicionador obtido por decomposição LU (Low Up) da matriz de stiffness. As voltagens são calculadas por um algoritmo, mais preciso, recentemente desenvolvido. Os perfis são gerados por um novo algoritmo envolvendo uma sucessiva troca de resistividade de subdomínios. Este procedimento permite obter cada nova matriz de stiffness a partir da anterior pelo cálculo, muito mais rápido, da variação dessa matriz. Este método permite ainda, acelerar a solução iterativa pelo uso da solução na posição anterior da ferramenta. Finalmente geramos perfis sintéticos afetados por cada um dos efeitos para um modelo da ferramenta Dual Laterolog.