MEF - Aplicando em 2D

Introduo, teoria, aplicaes e programao
usando Python
Prof. Marco Andr Argenta

Prof. Marcos Arndt
Desafios:
1. As interfaces so arestas ou faces;
2. A continuidade dos parmetros nodais no implica
necessariamente em continuidade das variveis

interpoladas nas interfaces;
3. Obteno mais elaborada das funes de interpolao;
4. Necessidade de maior flexibilidade na definio da
geometria dos elementos.
29/07/2012
Mtodo dos Elementos Finitos I
Arbtrio do campo de deslocamentos na

forma de interpolao de seus
deslocamentos nodais -> matriz N
Definio da geometria do elemento com as

mesmas funes de interpolao N do campo de
deslocamentos (definio isoparamtrica)
Formulao de
Deslocamentos
Obteno das relaes deformaodeslocamentos nodais -> matriz B
Especificao das relaes tensodeformaes -> matriz D
Soriano (2009)
29/07/2012
Integrao da matriz de rigidez e dos vetores de

foras nodais equivalentes, para obter o sistema
Keue = fe
Campo de deslocamentos:
u u v w em 3D ou u u v em 2D
T
Aproximao do campo de deslocamentos:

u Nu(e)
N : matriz de f. de forma
Relaes deformao-deslocamentos:
Lu
L : matriz de operadores diferenciais
Relaes tenso-deformaes:
D 0
29/07/2012
D : matriz constitutiva
s0 : vetor de tenses iniciais
Funcional Energia Potencial Total:

1 T
T
J (u) ( D 0 )dV pT udV q T udS
2
V
V
Sq
foras de
volume
foras de
superfcie
1
T
J (u) ( (Lu )T DLu 0 Lu )dV pT udV q T udS
2
V
V
Sq
Condio de estacionariedade:
J (u) ((Lu)T DLu (Lu)T 0 )dV uT pdV uT qdS 0
V
Sq
u Nu(e)
29/07/2012
Condio de estacionariedade:
J (u) ((LN u ( e) )T DLNu (e) (LN u (e) )T 0 )dV (Nu ( e) )T pdV (Nu (e) )T qdS 0
V
Sq
B LN
J (u) ((Bu ( e) )T DBu( e) (Bu ( e) )T 0 )dV (Nu (e ) )T pdV (Nu (e) )T qdS 0
V
( e)T
( B DBdV )u
T
(e)
( e )T
Sq
( BT 0 dV NT pdV NT q dS )
V
Sq
Vlida para quaisquer variaes:

T
(e)
B DBdV u BT 0 dV NT pdV NT q dS
V
V
Sq
V
29/07/2012
K ( e )u ( e ) f ( e )
Forma variacional fraca:
T
(e)
V
V
Sq
V
K ( e ) BT DBdVe
Matriz de rigidez
Ve
fs 0 BT 0 dVe
(e)
Tenses iniciais
Ve
(e)
fs 0 f p f q
(e)
(e)
(e)
f p NT pdVe
(e)
Foras de volume
Ve
f q NT q dSe
(e)
Foras de superfcie
Se
29/07/2012
Equao diferencial de equilbrio:

LT p 0
Mtodo dos resduos ponderados com as funes de

forma N:
T
T
(L p)dV 0
Aplicando o Teorema de Green ao 1o termo:

T
T
T
N
n
dS
(
LN
)
dV
pdV 0
S
normal externa ao contorno
n q
T
T
T
(
LN
)
dV
N
p
dV
qdS
V
29/07/2012
Como: D 0
Lu e u Nu(e)
T
(e)
T
T
T
(
LN
)
DLNu
dV
(
LN
)
s
dV
N
p
dV
N
0
qdS
V
B LN
T
(e)
V
V
S
V
K (e )
(e)
fs 0
(e)
fp
(e)
fq
K ( e )u ( e ) f ( e )
29/07/2012
Elemento finito triangular de trs pontos nodais: mais

verstil e simples em elasticidade bidimensional, embora
requeira grande refinamento de malha para resultados
acurados.
-Pontos nodais numerados

no sentido anti horrio;
- 2 deslocamentos por n
29/07/2012
10
Solues propositivas lineares em termos

generalizados:
1

u ( x, y ) 1 x y 2
v( x, y ) 1 x

3
de parmetros
4

y 5

6
Troca para parmetros nodais:

u1 1 x1

u2 1 x2
u 1 x
3
3
29/07/2012
y1 1

y2 2
y3 3
u ( x, y ) 1 x
1 x1
y 1 x2
1 x3
y1
y2
y3
u1

u2
u
3
11
u ( x, y ) 1 x
x2 y3 y2 x3
1
y
y2 y3
2 Ae
x2 x3
x1 y3 y1 x3
y1 y3
x1 x3
x1 y2 y1 x2 u1

y1 y2 u2
x1 x2
u3
onde Ae a rea do tringulo expressa por

1 x1
1
Ae 1 x2
2
1 x3
29/07/2012
y1
1
y2 ( x2 y3 x1 y2 x3 y1 x2 y1 x1 y3 x3 y2 )
2
y3
12
u ( x, y ) N1
N2
u1

N 3 u2
u
3
v( x, y ) N1
N2
v1

N 3 v2
v
3
1
Ni
(ai bi x ci y )
2 Ae
a1 x2 y3 y2 x3 a2 x1 y3 y1 x3
a3 x1 y2 y1 x2
y
,
b
y
e
1
2
b3 y1 y2
2
3
1
3
c x x
c x x
c x x
2
3
2
1
3
1
2
1
3
29/07/2012
13
1
Ni
(ai bi x ci y )
2 Ae
Tm as propriedades do delta de Kronecker:

N i ( x j , y j ) ij com i, j 13
3
3
3
N 1,
N i , x 0, N i , y 0
i 1
i 1
i 1
29/07/2012
14
10/07/2013
15
Armazenando todos os deslocamentos nodais em um

T
(e)
nico vetor:
u u1 v1 u2 v2 u3 v3
u N1
u
v 0
03/08/2012
0
N1
N2
0
0
N2
N3
0
u1
v
1
0 u2
(e)
Nu

N 3 v2
u3

v3
16
Como o campo de deslocamentos deste elemento linear

e em cada um dos lados do elemento h dois parmetros
nodais (gl) para definir o deslocamento u e dois
parmetros nodais para definir o deslocamento v, a
interface entre dois elementos permanece retilnea aps a
deformao. Ou seja, no h descontinuidade de
deslocamentos entre elementos.
Como ocorre continuidade de deslocamentos atravs da
interface, diz-se que o elemento compatvel ou
conforme.
aconselhvel utilizar ngulos entre 30 e 150 e razes
entre segmentos de lados menores do que quatro para
adequada discretizao.
03/08/2012
17
Estado plano de tenso: Slido de espessura pequena,

superfcie mdia plana e aes externas atuantes apenas
nesse plano (plano xy).
s x 0, s y 0, s z 0
xy 0, xz 0, yz 0
Estado plano de deformao: Slidos longos, de apoio e

de sees transversais e aes externas constantes na direo
longitudinal z, alm de extremidades supostas indeslocveis
nessa direo. Pode-se estudar apenas o comportamento de
um segmento de espessura unitria.
x 0, y 0, z 0
xy 0, xz 0, yz 0
03/08/2012
18
Equaes de Compatibilidade (relaes deformaodeslocamento):

x u, x x
0
u

Lu y v, y 0
y
u v y x v
xy , y , x
03/08/2012
19
Equaes Constitutivas (relaes tenso-deformao):

D 0
Estado plano de tenso:

s x
E

D s y
2
1

xy
0
1
1
0
0 0 (1 )
x

y
2 xy
(s x s y )
Estado plano de deformao: s z (s x s y )

s x
0
(1 )
E

D s y
(1 )
0
(1 )(1 2 ) 0
0
(1 2 )
xy
03/08/2012
x

y
2 xy
20
B LN
0
x
N1
B 0
y
0
y x
1
Ni
(ai bi x ci y )
2 Ae
N2
N3
N1
N2
b1 0 b2
1
B
0 c1 0
2 Ae
c1 b1 c2
0
c2
b3
0
b2
c3
0
N 3
0
c3
b3
A matriz B tem coeficientes constantes. Isto implica em

estado de deformao constante e consequentemente em
tenso constante. Por este motivo o elemente conhecido
como CST (constant strain triangle).
03/08/2012
21
Matriz de Rigidez Elementar:

K ( e ) BT DBdVe
Ve
K ( e ) BT DB tdAe tBT DB dAe K ( e ) tAe BT DB

Ae
Ae
Estado plano de tenso: para espessura t constante

d
03/08/2012
1
2
22
Estado plano de deformao: para espessura unitria

K (e) Ae BT DB
1 2
e
2
03/08/2012
23
Vetor de foras equivalentes:

- Foras de volume uniformes:
fp
(e)
N pdVe t N pdAe
T
Ve
fp
03/08/2012
(e)
Ae
p x N1dAe

p y N1dAe
p
N
dA
tAe
x
2
e
t
p y N 2 dAe 3
p N dA
x 3 e
p y N 3 dAe
px
p
y
p x

py
px

p y
px
p
py
N dA
i
vol. pirmide 1 3 Ae h
24

- Foras de superfcie uniformes:
fq
(e)
N q dSe
T
Se
q x
q
q y
para o lado 1-3:
fq
03/08/2012
(e)
q x N1dS
q x

q
q y N1dS
y
q x N 2 dS l13 0

q y N 2 dS 2 0
q N dS
q x
x
3

q y
q y N 3 dS
N dS rea do tringulo 1 2lh

i
25
Obteno das tenses:

D
Lu
u Nu(e)
DLu DLNu (e) DBu(e)
u1
v
1
s x
u2

s y DB

v2
xy
u3

v3
03/08/2012
26
Elemento finito retangular de quatro pontos nodais:

utilizado apenas em modelos com contorno retangular,
mas com comportamento melhor que o CST por conter
mais termos no polinmio interpolador.
-Pontos nodais numerados

no sentido anti horrio;
- 2 deslocamentos por n
03/08/2012
27
Solues propositivas lineares em termos de parmetros

generalizados:
5
1

u ( x, y ) 1 x y xy 2
v( x, y ) 1 x y xy 6
3
7
4
8
Troca para parmetros nodais: usando um sistema local
de coordenadas com origem no n 4.
u1 1 a 0 0 1
u 1 a b ab
2
2

u3 1 0 b 0 3
u4 1 0 0 0 4
03/08/2012
28
Troca para parmetros nodais:

u ( x, y ) 1 x
u ( x, y ) 1 x
03/08/2012
1
1
xy
1
a
a
0
0
0 0 u1
b ab u2

b 0 u3

0 0 u4
0
0
1 u1
0
1a
u
0
0
1
a
2
xy

0
0
1b
1 b u3
1
ab
1
ab
1
ab
1
ab
u4
29
u ( x, y ) N1
N1
N2
u1
u

N 4 2
u3
u4
N2
N3
N3
x xy
N4
a ab
v( x, y ) N1
xy
ab
N2
N3
v1
v

N 4 2
v3
v4
y xy
x y xy
1
b ab
a b ab
Funes de forma bilineares: quadrticas no domnio do

elemento e lineares no seu contorno.
03/08/2012
30
N1
N3
03/08/2012
N2
N4
31
Como o campo de deslocamentos em cada um dos
lados do elemento linear e h dois parmetros nodais

(gl) para definir o deslocamento u e dois parmetros
nodais para definir o deslocamento v, a interface entre
dois elementos permanece retilnea aps a deformao.
Ou seja, no h descontinuidade de deslocamentos
entre
elementos,
satisfazendo
o
critrio
de
conformidade.
O campo de deslocamentos contm o termo xy, sendo
completo e mais refinado que o adotado no CST.
Conduz a melhores resultados do que dois elementos
CST que formem o mesmo retngulo.
03/08/2012
32
Armazenando todos os deslocamentos nodais em um

nico vetor: u(e) u1 v1 u2 v2 u3 v3 u4 v4 T
u N1
u
v 0
03/08/2012
N2
N3
N4
N1
N2
N3
u1
v
1
u2

0 v2
(e)
Nu

N 4 u3
v3

u4
v
4
33
Lu LNu (e)
B LN
0
x
N1
B 0
y
0
y x
1 y
a ab
B 0
x
ab
03/08/2012
y
ab
0
x
ab
1 y
a ab
0
x
ab
N2
N3
N4
N1
N2
N3
0
x
ab
y
ab
ab
0
1 x
b ab
0
1 x
b ab
y
ab
1 y

a ab
0
1 x

b ab
0
N 4
1 x

b ab
1 y

a ab
0
34
Matriz de Rigidez Elementar:

K ( e ) BT DBdVe
Ve
K ( e ) t BT DBdAe
Ae
a b
K ( e ) t BT DBdydx
0 0
03/08/2012
35
Estado plano de tenso: para espessura t constante
03/08/2012
36
Estado plano de deformao: para espessura unitria
03/08/2012
37

a b
(e)
T
T
f
N
p
dV
t
N
- Foras de volume uniformes: p
e
pdydx
Ve
px
p
py
03/08/2012
fp
(e)
p x N1dydx
Ae
p y N1dydx
Ae
px
p x N 2 dydx
p
Ae
y
p y N 2 dydx tab p x
t Ae

p x Ae N 3 dydx 4 p y
px
p
N
dydx
y Ae 3

p y
p x Ae N 4 dydx
p
N
dydx
y Ae 4
0 0
38

- Foras de superfcie uniformes:
para o lado 1-2:
q x
q
q y
03/08/2012
fq
(e)
fq
q x N1dS
q x

q
q y N1dS
y
q x
q x N 2 dS

q N dS
y 2 b q y

q x N 3 dS 2 0
q N dS
0
y
3

q x N 4 dS
0
0

q y N 4 dS
(e)
NT q dSe
Se
N dS rea do tringulo 1 2lh

i
39
Obteno das tenses:

D
Lu
DLu DLNu (e) DBu(e)
03/08/2012
u Nu(e)
u1
v
1
u2
s x

v2
s y DB

u3
xy
v3

u4
v
4
40
Montagem do Sistema Global de Equaes:
10/07/2013
41
Montagem do Sistema Global de Equaes:
10/07/2013
42

MEF - Aplicando em 2D

Enviado por

Direitos autorais:

Formatos disponíveis

MEF - Aplicando em 2D

Enviado por

Dados do documento

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

MEF - Aplicando em 2D

Enviado por

Direitos autorais:

Formatos disponíveis

Introduo, teoria, aplicaes e programao

Prof. Marco Andr Argenta

necessariamente em continuidade das variveis

Mtodo dos Elementos Finitos I

Arbtrio do campo de deslocamentos na

Definio da geometria do elemento com as

Obteno das relaes deformaodeslocamentos nodais -> matriz B

Especificao das relaes tensodeformaes -> matriz D

Integrao da matriz de rigidez e dos vetores de

Aproximao do campo de deslocamentos:

L : matriz de operadores diferenciais

Mtodo dos Elementos Finitos I

Funcional Energia Potencial Total:

Mtodo dos Elementos Finitos I

Vlida para quaisquer variaes:

Mtodo dos Elementos Finitos I

Forma variacional fraca:

Mtodo dos Elementos Finitos I

Equao diferencial de equilbrio:

Mtodo dos resduos ponderados com as funes de

Aplicando o Teorema de Green ao 1o termo:

normal externa ao contorno

Mtodo dos Elementos Finitos I

Mtodo dos Elementos Finitos I

Elemento finito triangular de trs pontos nodais: mais

-Pontos nodais numerados

Mtodo dos Elementos Finitos I

Solues propositivas lineares em termos

Troca para parmetros nodais:

Mtodo dos Elementos Finitos I

onde Ae a rea do tringulo expressa por

Mtodo dos Elementos Finitos I

Mtodo dos Elementos Finitos I

Tm as propriedades do delta de Kronecker:

Mtodo dos Elementos Finitos I

Mtodo dos Elementos Finitos I

Armazenando todos os deslocamentos nodais em um

Mtodo dos Elementos Finitos I

Como o campo de deslocamentos deste elemento linear

Mtodo dos Elementos Finitos I

Estado plano de tenso: Slido de espessura pequena,

Estado plano de deformao: Slidos longos, de apoio e

Mtodo dos Elementos Finitos I

Equaes de Compatibilidade (relaes deformaodeslocamento):

Mtodo dos Elementos Finitos I

Equaes Constitutivas (relaes tenso-deformao):

Estado plano de tenso:

Estado plano de deformao: s z (s x s y )

Mtodo dos Elementos Finitos I

A matriz B tem coeficientes constantes. Isto implica em

Mtodo dos Elementos Finitos I

Matriz de Rigidez Elementar:

K ( e ) BT DB tdAe tBT DB dAe K ( e ) tAe BT DB

Estado plano de tenso: para espessura t constante

Mtodo dos Elementos Finitos I

Estado plano de deformao: para espessura unitria

Mtodo dos Elementos Finitos I

Vetor de foras equivalentes:

Mtodo dos Elementos Finitos I

Vetor de foras equivalentes: