Método de Elementos Finitos

Em construção - as referencias e imagens ainda estão faltando

Grupo: Antônio Carlan, Gabriela Pereira, Renan Soares e Victor Gandara

Objetivo: Apresentar uma introdução do Método de Elementos Finitos, abordando seu funcionamento e demonstrando suas etapas através de um exemplo, demonstrando seu potencial e o que o diferencia de outros métodos estudados nas cadeiras de métodos computacionais até então vistos.

Introdução ao Método de Elementos Finitos(MEF)

Aplicação

Este é um método numérico de resolução aproximada para problemas de valores sobre contorno em equações diferenciais. É largamente utilizado na industrial para o desenvolvimento de novas tecnologias e produtos, pode ser utilizado tanto no processo de desenvolvimento quanto de otimização, através da simulação da física de situações em que se espera que o projeto atue.

[imagem1] Simulação de deformação de um veiculo https://en.wikipedia.org/wiki/Finite_element_method

O que torna este método diferente é a forma com que ele discretiza e faz suas aproximações. Sua discretização é feita na geometria, ou seja, repartindo o objeto em diversos elementos menores. Desta forma é possível que se resolva problemas em geometrias complexas que seriam muito difíceis ou até impraticáveis com outros métodos. Também permite trabalhar com materiais compostos, ou de diferentes densidades em um objeto.

[imagen2] Modelo normal e modelo discretizado e representado por uma malha. [imagem3] Simulação de estresse sofrido no objeto.

O Método

O conceito fundamental do MEF é a aproximação de uma quantidade continua(p. ex., tensão elétrica, temperatura, etc.) através de um modelo discreto composto por um conjunto de funções simples que são definidas em um número finito de subdomínios(elementos)

A resolução de um problema através do MEF envolve quatro etapas principais:

1. Discretização da região: Transformar o objeto de estudo em vários elementos menores que juntos formam o objeto original de forma aproximada e a isso chama-se Malha.
2. Aproximação da solução: Obtenção das equações correspondentes de cada elemento, através de interpolação polinomial.
3. Montagem do sistema global: Juntar os elementos e suas equações em uma matriz que fornecerá a solução global.
4. Resolução do sistema: Resolver a matriz gerada para obter a solução aproximada do problema.

Discretização: A Malha

A malha define o domínio do problema, ela é a base do método. Porém não é parte do algoritmo em si. Normalmente é criada automaticamente por um programa especifico para isso, como GMSH e NetGen. Apenas sendo fornecida as especificações que se deseja para a malha, vão influenciar na quantidade de nós, na dimensão do domínio e na precisão da aproximação.

Ela é composta por um conjunto de pontos interligados que descrevem o objeto. Estes pontos são chamados de nós e a união destes nós forma um elemento, e a união destes elementos formam o objeto que é nosso domínio.

[imagem nos.png]

1. Dimensão: A malha pode ser unidimensional, bidimensional ou tridimensional.

[imagem dimensão print slide 5]

1. Ordem: Quanto maior a ordem da malha mais precisa será a aproximação dos resultados entre os nós, porém aumenta o consumo computacional.

[imagem ordem print slide 5]

1. Densidade de elementos: Pode-se escolher por gerar uma malha mais ou menos detalhada, isso implica no consumo computacional. Quanto mais detalhada a malha, ou seja, quanto maior o numero de nós da malha, menor será o erro de aproximação do problema real.

[imagem densidade de malhas]

Funcionamento do método

Cada nó desta malha terá atribuída uma variável do problema físico que esta sendo implementado(o valor da grandeza contínua para aquele ponto), como Tensão Elétrica, Dilatação, Temperatura, etc. Esta variável é chamada Valor nodal, ela representa a solução da equação naquele ponto (Variáveis nodais = incógnitas ou graus de liberdade do problema).

[imagem print slide 7]

Já a grandeza contínua é aproximada em cada elemento por um polinômio que é definido usando os valores nodais da grandeza. Ou seja, nós iremos descobrir os valores nodais dos nós e entre eles,os valores da grandeza contínua serão dados por uma aproximação polinomial.

Lembrando que os elementos são formados pelos nós e que os valores nodais estão associado aos nós e não ao elemento diretamente.

Então, cada elemento terá um polinômio que descreve a grandeza contínua na região entre os nós que formam este elemento, no caso de uma matriz com elementos de primeira ordem, será uma reta entre os valores nodais. Os polinômios são escolhidos/deduzidos de forma que a continuidade seja mantida na fronteira entre os elementos, utilizando os métodos matemáticos adequados para isso.

"Este polinômio resulta de um ajuste dos valores nodais para que se tenha uma boa aproximação da grandeza real, este ajuste é feito através da minimização de alguma grandeza associada ao problema."

Interpolação

Os valores nodais, que são nossas incógnitas, precisam ser isolados e deixados em função das coordenadas dos nós, que possuímos através da malha de elementos finitos já gerada.

[imagem gráfico slide 9]

$V_{1}eV_{2}$ -> valores nodais $x_{1}ex_{2}$ -> coordenadas dos nós $L=x_{2}-x_{1}$

Aproximação polinomial no caso da malha de primeira ordem:

$V(x)=a_{1}+a_{2}x$

As variáveis a1 e a2 são coeficientes da interpolação, porém não as temos e não queremos deixa-las no calculo, para isso:

$V_{1}=a_{1}+a_{2}x_{1}$

$V_{2}=a_{1}+a_{2}x_{2}$

logo:

$<center><math>{\begin{bmatrix}V_{1}\\V_{2}\\\end{bmatrix}}={\begin{bmatrix}1&x_{1}\\1&x_{2}\\\end{bmatrix}}{\begin{bmatrix}a_{1}\\a_{2}\\\end{bmatrix}}$

</math>

Portanto:

$a_{1}={\frac {1}{L}}(V_{1}x_{2}-V_{2}x_{2})$

$a_{2}={\frac {1}{L}}(V_{1}-V_{2})$

Substituindo os valores de $a_{1}$ e $a_{2}$ na equação de interpolação anterior $V(x)=a_{1}+a_{2}x$ ficaremos com:

$V(x)={\frac {1}{L}}[V_{1}x_{2}-V_{2}x_{1}+(V_{2}-V_{1})x]$

Ou, como é definido a aproximação nodal:

$V(x)={\frac {(x_{2}-x)}{L}}V_{1}+{\frac {(x-x_{1})}{L}}V_{2}$

Esta fica conhecida por aproximação nodal pois a grandeza $V(x)$ pode ser obtida em qualquer ponto do intervalo $[x_{1},x_{2}]$ a partir dos valores $V_{1}$ e $V_{2}$.