Introduo, teoria, aplicaes e programao
usando Python
Prof. Marco Andr Argenta
Prof. Marcos Arndt
Refinamentos do MEF:
1. Refinamento h: aumento do nmero de elementos que
compem a malha;
2. Refinamento p: aumento do grau e/ou do nmero das
funes de forma no elemento sem alterar a malha. No
MEF convencional, o refinamento p corresponde ao
aumento do grau do polinmio interpolador da soluo;
3. Refinamento hp: combinao do refinamento da malha
(refino h) simultaneamente com a variao na ordem
do polinmio aproximador (refino p).
25/07/2012
Mtodo dos Elementos Finitos I
Consistem na associao de funes vrtice e funes
bolha.
Funo vrtice: No nula em uma das extremidades do
intervalo (elemento) e nula na outra extremidade.
Relacionada aos ns da malha. Seu suporte corresponde
a 2 elementos vizinhos.
(Soln et al., 2004)
25/07/2012
Mtodo dos Elementos Finitos I
Funo bolha: Se anula em ambas as extremidades do
intervalo (elemento).
(Soln et al., 2004)
25/07/2012
Mtodo dos Elementos Finitos I
Definio:
Seja u uma funo desconhecida e sua aproximao dada
n
por
u u N i ai
i 1
A aproximao hierrquica se um aumento de n para
n+1 no altera as funes de forma Ni (i = 1 a n).
As funes de forma hierrquicas foram introduzidas por
Zienkiewicz et al por volta de 1971.
(Zienkiewicz, Gago e Kelly, 1983)
25/07/2012
Mtodo dos Elementos Finitos I
Em um refinamento p, se o conjunto de funes de forma
de uma aproximao de ordem p constitui um subconjunto
do conjunto de funes de forma de uma aproximao de
ordem p+1, este refinamento denominado hierrquico.
25/07/2012
Mtodo dos Elementos Finitos I
Vantagens: (Campion e Jarvis, 1996)
1. Reteno dos coeficientes da matriz de rigidez quando
a ordem da interpolao aumenta;
2. Altas taxas de convergncia sem necessidade de
refinar a malha;
3. Melhora do condicionamento das equaes.
Desvantagens: (Leung e Chan (1998), e Ribeiro (2001))
1. Novos g.l. no tem significado fsico direto;
2. Polinmios de alta ordem podem conduzir a mal
condicionamento do sistema de equaes.
25/07/2012
Mtodo dos Elementos Finitos I
2 ns: polinmio de ordem 1
N2
N1 1
3 ns: polinmio de ordem 2
N1 2 3 1
2
25/07/2012
N 2 4 4
Mtodo dos Elementos Finitos I
N3 2 2
3 ns: polinmio de ordem 2 hierrquico
N2
N1 1
N3 4 2 4
25/07/2012
Mtodo dos Elementos Finitos I
Podem ser obtidos a partir da integrao dos polinmios
de Legendre. No caso de um domnio x = [1, 1], os
polinmios de Lobatto tem a seguinte forma:
2k 1
k 1
L0 ( x) 1 ; L1 ( x) x ; Lk ( x)
xLk 1 ( x)
Lk 2 ( x) p/ k 2
k
k
Lk 1 ( x) 2 L2k 1 ( x)dx 2 2k 1
1
(Solin et al., 2004)
25/07/2012
Mtodo dos Elementos Finitos I
10
25/07/2012
Mtodo dos Elementos Finitos I
11
25/07/2012
Mtodo dos Elementos Finitos I
12
Polinmios K ortogonais:
Simple hierrquico:
Nodal : Lagrange
linear
(1 ) k 1 (1 )
k
4
Condicionamento da matriz de
rigidez do elemento mestre
25/07/2012
Mtodo dos Elementos Finitos I
13
Conjunto de funes de forma composto pelos 4 primeiros
polinmios de Hermite (MEF clssico) e por polinmios
especiais de Legendre na forma de Rodrigues (Peano,
1976):
f r ( )
( r 1) / 2
n 0
(1) n (2r 2n 7)!! ( r 2 n 1)
n
2 n!(r 2n 1)!
onde r!! r (r 2)(r 4)(2 ou 1) ; 0!! (1)!! 1
(r 1) / 2 denota a parteinteira desta operao
Estas funes e suas derivadas primeiras se anulam nas
extremidades do elemento (f.f. tipo bolha). (Bardell, 1996)
25/07/2012
Mtodo dos Elementos Finitos I
14
Bardell, N. S.; An engineering application of the h-p version of the finite element
method to the static analysis of a Euler-Bernoulli beam. Computers & Structures,
v. 59, n. 2, p. 195-211, 1996.
25/07/2012
Mtodo dos Elementos Finitos I
15
25/07/2012
Mtodo dos Elementos Finitos I
16
Polinmios K ortogonais:
25/07/2012
Mtodo dos Elementos Finitos I
17
01/07/2013
Mtodo dos Elementos Finitos I
18
01/07/2013
Mtodo dos Elementos Finitos I
19
01/07/2013
Mtodo dos Elementos Finitos I
20