Introducción al Método de los Elementos Finitos
2
Algunas aplicaciones del MEF a problemas elípticos
Problema de Elasticidad
Consideremos un cuerpo elástico, isótropo y homogéneo B,
que ocupa el dominio acotado O?R3, con frontera
G=Gt?Gu tal que GtnGu=Ø y área Gu>0.
El cuerpo B está sometido a una fuerza volumétrica f,
y a una fuerza superficial t aplicada sobre Gt.
Se supone B fijo a lo largo de Gu.
Se busca determinar
desplazamientos u.
deformaciones e, dependientes de los desplazamientos de acuerdo a la
cinemática de la deformación.
tensiones s, dependientes de las deformaciones de acuerdo a la ley
constitutiva del material.
O
O
Gu
Gt
?u =0
sobre Gu
1? ?ui ?u j ?
??x
?xi ? ?
= ??e d +µe
?,µ?
, ctes. de Lamé.
, µ =
?: coeficiente de Poisson
1+?
(1+?)(1-2?)
n =?s n
Introducción al Método de los Elementos Finitos
3
Problema de Elasticidad
Nota: en adelante, usaremos las siguientes convenciones de notación:
Ecuación de equilibrio
CB Dirichlet (despl. impuesto)
CB Neumann (tracción impuesta)
?
?divs + f =0 en O
?
?sn = t sobre Gt
O
O
Gu
Gt
n
Ecuaciones de clausura
Cinemáticas: asumiendo pequeñas deformaciones:
eij = ? + ?
2? j
Constitutivas: asumiendo comportamiento elástico lineal (ley de Hooke):
,
– derivada parcial: ui, j =
?ui
?xj
3
– sumatoria: sij j ij j
j=1
+
E E? E :módulo de Elasticidad
3
sij kk ij ij,
k=1
con ? =
{
con V= v:v?? ? ? y v = 0 en Gu
?H (O)
1
? ?
dx+? ?
-?
2
+v
vi, j j,i
?
– sij
de pequeñas
2
-? sij ij ? ivids+? fividx = 0
Introducción al Método de los Elementos Finitos
4
Forma variacional del problema de Elasticidad
Dado
se llega a
haciendo
}
3
(V)
Hallar u?V/a(u,v) = L(v), ?v?V
(D)
sij, j + fi = 0 en O, i, j =1,2,3
O O
G O
O Gt O
sij, jvidx+
fividx = 0
sijnjvids+
fividx = 0
-?Osijvi, j
e (v)dx+ t
d +µe t v ds+
O Gt i i O
T. de Green
Ley de Hooke
CB
dx+?Gt tivids+?O fividx = 0
dx+?Gt tivids+?O fividx = 0
sijvi, j +s jivj,i
O
Simetría de s
Cinemática
O
deformaciones
L(v)
a(u,v)
i.e.,L(v) = ? v V ,?v?V,?? +.
a(u,v) =? ??ui,i ij ij ? ij(v)dx =? ??divudivv +µeij ij(v)?dx
continua, i.e., a(u,v) =? u v ,?u,v?V,? ?
H (O) H (O)
V-elíptica, i.e., a(v,v) =a v ,?v?V,a ?
, con v
+
H (O)
Demo.: a(v,v) = ?? (divv) dx+µ? eij(v)eij(v) dx
H (O)
=? vi
H (O)
= µ? eij ij(v) dx = µc v
H (O)
Introducción al Método de los Elementos Finitos
5
.
1
1
2
2
.
3
i-1
1
1 1
2
es simétrica, i.e., a(u,v) = a(v,u),?u,v?V.
Forma variacional del problema de Elasticidad (cont.)
Se puede demostrar que la forma lineal L(v) =?O fividx+?G tivids es continua,
t
Se puede demostrar que la forma bilineal
d +µe (u)?e (u)e
O O
1
2
2
+
.
+
O O
O
, c?
(v)e
Desigualdad de Korn
{
}
Vh 1 (K)? ? ,?K?Th
= v:v?V y v K ?? ?P
H (O)
Introducción al Método de los Elementos Finitos
6
1
u-uh
= Ch u H2(O)
MEF aplicado al problema de Elasticidad
Consideremos el problema de Elasticidad en O?R3.
Sea Th={K} una malla de tetraedros de O. Definimos el espacio de EF
3
El MEF aplicado al problema de Elasticidad consiste en
Hallar uh ?Vh /a(uh,v) = L(v), ?v?Vh
La solución uh?Vh satisface
{
}
V = v:v?? ?H0 1(O)? ? y divv = 0 en O
Introducción al Método de los Elementos Finitos
7
Problema de Stokes
Consideremos las ecuaciones de Stokes para el flujo estacionario de un fluido
Newtoniano incompresible encerrado en un dominio O?R3, sometido a una
fuerza volumétrica f :
sij, j + fi = 0 en O, Balance de cant. de movto.
en O, Ley const. de fluido Newtoniano
ui,i = 0 en O, Condición de incompresibilidad
ui = 0 sobre G, CB Dirichlet
u: velocidad
sij = 2µeij(u)- pdij
s: tensión
p: presión
µ: viscosidad
-µ?ui + p,i = fi
en O, Balance de cant. de movto. p/fluido Newtoniano
Definimos el espacio de funciones de prueba
3
Luego, podemos llevar el problema de Stokes a la forma variacional
(V)
Hallar u?V/a(u,v) = L(v), ?v?V
fi i v dx = µ? ?ui i v dx-µ?
?ui
Introducción al Método de los Elementos Finitos
8
Dado que µ >0, se demuestra (ídem problema de Poiss
Página siguiente |