Teórico Método de Elementos Finitos
Teórico Método de Elementos Finitos
Teórico Método de Elementos Finitos
Introduccin al
Mtodo de los Elementos Finitos
(FEM)
Clculo Estructural II
FCEFyN Dpto. Estructuras
Juan Gir - 2014
Estructuras aeronuticas
Diseo de elementos de mquinas
Estructuras civiles y mecnicas
Transferencia de calor
Aerodinmica y mecnica de los fluidos
Geomecnica
Hidrulica
Mquinas elctricas y electromagnetismo.
Bioingeniera.
UNC FCEFyN
CE II - 2014
Sistema
Fsico
Modelo
Conceptual
Discretizacin
Modelo
Matemtico
Solucin
Modelo
Discreto
Solucin
Discreta
Errores de la
Solucin
Errores de la Discretizacin y
Solucin
Errores de la Formulacin Matemtica, Discretizacin
y Solucin
Errores de la Idealizacin, Formulacin Matemtica, Discretizacin y
Solucin
Verificacin
Validacin
UNC FCEFyN
CE II - 2014
UNC FCEFyN
CE II - 2014
Desarrollo de un elemento
UNC FCEFyN
CE II - 2014
Funciones de Aproximacin
Las funciones de aproximacin de los
desplazamientos son claves en la definicin de
los elementos, ya que son las que permiten
representar el comportamiento en el interior del
dominio en funcin del desplazamiento de unos
pocos nudos. Estas deben ser:
1) Continuas y derivables
2) Integrables
3) Representar correctamente los
desplazamientos en todo el dominio
4) Representar el movimiento de cuerpo rgido.
UNC FCEFyN
CE II - 2014
Coordenadas Triangulares
a1i
x2
K
(i)
(j)
AJ
P
AK
a2 i
Ai
J
(k)
x1
Ai ( x1 , x2 )=
1
1
( JK JP )=
2
2
A j ( x1 , x2 )=
1
1
( KI KP )=
2
2
Ak ( x1 , x2 )=
UNC FCEFyN
CE II - 2014
1
( IJ IP )=
2
1
2
t1
t2
a1i
x1 x1j
t1
a1j
x1 x1k
t1
a1k
x1 x1i
a2i
x2 x2j
t2
a2j
x2 x2k
t2
a2k
x2 x2i
t3
0=
0
t3
0=
0
t3
0=
0
1
( x1j x1 ) a2i ( x2j x2 ) a1i
2
1
( x1k x1 ) a2j ( x2k x2 ) a1j
2
1
( x1i x1 ) a2k ( x2i x2 ) a1k
2
10
Coordenadas Triangulares
UNC FCEFyN
CE II - 2014
11
UNC FCEFyN
CE II - 2014
a1i
x2
K
(i)
(j)
AJ
P
AK
a2 i
Ai
J
(k)
x1
J
I
J
I
12
Ecuaciones Fundamentales
Ecuaciones de Equilibrio
11 21
+
+ F1 =
0
x1
x2
12
22
+
+ F2 =
0
x1
x2
Relaciones constitutivas
1
11
1
E
22 =
2
(1 ) 0 0
12
12
u1
u2
=
; 22
x1
x2
u2
1 u
= 1 +
x2
2 x1
11
22
12
Deformaciones en funcin de
los desplazamientos
Ecuaciones Cinemticas
11
=
0
0
1
2
=
11
u1 u1 i u1 j u1 k
=
+
+
x1 i x1 j x1 k x1
u1m a2m
1 i i
j j
k k
11 =u1 a2 + u1 a2 + u1 a2 =
2A
2A
1
k
u2i a1i + u2j a1j + u2k a1=
=
22
2A
u2m a1m
2A
13
2A
E
2
2
2
2
h.dA
=
+
+
+
Wn
2
2(1
)
11
22
12
12
12
2
2(1 ) A
=
W
n
Wn
=
h.E. A
2
2
2
2
2
2(1
)
(
)
11
22
12
12
12
2(1 2 )
Un =
F1 ( u1i i + u1j j + u1k k ) dA F2 ( u2i i + u2j j + u2k k ) dA
A
1dA
=
A
2 dA
=
A
3 dA
=
A
A
3
FA
FA
Un =
1 ( u1i + u1j + u1k ) 2 ( u2i + u2j + u2k )
3
3
UNC FCEFyN
CE II - 2014
14
0
;
0
;
0
=
=
=
u1i
u1j
u1k
n = Wn + U n
= W + U=
W + U
n
=
n 1=
n 1
u1i
=
=
=
0
;
0
;
0
u2i
u2j
u2k
22
22
11
12 U
h.E. A 11
2
+
+
+
+
2(1
)
11 i
+
22
11
22
12
u1
u1i
u1i
u1i
u1i u1i
2(1 2 )
22
22
11
12 U
h.E. A 11
2
=
+
+
+
+
2(1
)
11
+
22
11
22
12
u2i 2(1 2 ) u2i
u2i
u2i
u2i
u2i u2i
AFi
Eh
1 m m
m m
i
m m
i
m m
i
a
u
a
a
u
a
a
u
a
u
a
=
)
+
)
+
(
1 1
2
1 2
2
1 1
2 2
1
u1i 4 A(1 2 )
2
3
Eh
1 m i AFi
m i m 1 m i
m
m i
u
a
a
a
a
u
a
a
a1 a2 )
=
(
)
+
)
+
2
1
1
2
2
1
2
1
2
i
u2 4 A(1 )
2
2
3
UNC FCEFyN
CE II - 2014
15
kii12
kij11
kij12
kik11
kii22
k 12
ji
kij21
k 11
jj
kij22
k 12
jj
kik21
k 11
jk
k 22
ji
kki12
k 21
jj
kkj11
k 22
jj
kkj12
k 21
jk
11
kkk
kki22
kkj21
kkj22
kkk21
11
kim
=
CE II - 2014
Eh
m i 1 m i
a a +
a1 a1
2 2 2
4 A(1 )
2
12
kim
=
Eh
1 m i
m i
a
a
a2 a1
+
1 2
2
4 A(1 )
2
kim21
=
Eh
m i 1 m i
a1 a2
a2 a1 +
2
4 A(1 )
2
kim22
=
UNC FCEFyN
j j
k 22
jk u2
F2 3
12 k
F1k
kkk
u1
k
22 k
kkk
u2
F2
Eh
m i 1 m i
a a +
a2 a2
2 1 1
4 A(1 )
2
16
k41=14
12
k33
k3322
12
k43
22
k43
12
k53
k5322
11
k34
k3421
11
k44
21
k44
11
k54
k5421
12
k34
k3422
12
k44
22
k44
12
k54
k5422
11
k35
k3521
11
k45
21
k45
11
k55
k5521
12
k35
k3522
12
k45
22
k45
12
k55
22
k55
90
A
10
3
20
4
80
Eh
4 4 1 4 4
(80) 2 + (60) 2
aa +
a1 a1 =
( x23 x25 ) 2 + ( x13 x15 ) 2 =
2 2 2
4 A(1 )
2
Eh
4 4 1 4 4
1 2
k44
a2 a1 =
=
a1 a2 +
[ (60).(80). + (80).(60). ]
4 A(1 2 )
2
k442 1
=
k42=24
Eh
4 4 1 4 4
a1=
a2 [ (80).(60). + (60).(80). ]
a2 a1 +
4 A(1 2 )
2
Eh
4 4 1 4 4
3
5 2
3
5 2
2
2
aa +
a2 a2 =
(
x
x
)
+
(
x
x
)
=
(
60)
+
(
80)
1
1
2
2
2 1 1
4 A(1 )
2
UNC FCEFyN
CE II - 2014
17
Otros elementos
UNC FCEFyN
CE II - 2014
UNC FCEFyN
CE II - 2014
19
K
M
O
N
N
J
L
I
Grado de la funcin de
aproximacin
Grado de la funcin de
deformacin
Cantidad de Nudos
Grados de
libertad
Lineal
Constante
Cuadrtica
Lineal
12
Cbica
Cuadrtica
18
Otros elementos
I
L
L
UNC FCEFyN
CE II - 2014
K
I
20
Axil simtricos
UNC FCEFyN
CE II - 2014
Bandas Finitas
21
UNC FCEFyN
CE II - 2014
22
Tipos de elementos
Grados de libertad
UNC FCEFyN
CE II - 2014
Funciones de
aproximacin
de los
desplazamientos
23
Ensamble de elementos
UNC FCEFyN
CE II - 2014
Armado de la
Matriz de Rigidez
global
K1112
K1122
K1211
K1221
K1212
K1222
0 0 K1411
0 0 K1421
K1412
K1422
12
K 21
K 2122
0
0
11
K 22
K 2221
0
0
12
K 22
K 2222
0
0
0
0
0
0
11
0 K 24
0 K 2421
0 0
0 0
12
K 24
K 2422
0
0
12
K 41
K 4122
0
0
11
K 42
K 4221
0
0
12
K 42
K 4222
0
0
0
0
0
0
11
0 K 44
0 K 4421
0 0
0 0
12
K 44
K 4422
0
0
0
0
0
0
0
0
0 0
0 0
K1111
21
K11
11
K 21
21
K 21
0
0
K = 11
K
41
K 4121
0
0
0
0
0
0
0
0
C
4
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
UNC FCEFyN
CE II - 2014
Armado de la
Matriz de Rigidez
global
C
4
K1111 + K1111
21
21
K
+
K
11
11
11
K 21
21
K 21
K 11
31
21
K 31
K = 11
11
K 41 + K 41
K 4121 + K 4121
K1112 + K1112
K1122 + K1122
K1211
K1221
K1212
K1222
K1311
K1321
K1312
K1322
K1411 + K1411
K1421 + K1421
K1412 + K1412
K1422 + K1422
12
K 21
K 2122
12
K 31
K 3122
11
K 22
K 2221
0
0
12
K 22
K 2222
0
0
0
0
11
K 33
K 3321
0
0
12
K 33
K 3322
11
K 24
K 2421
11
K 34
K 3421
12
K 24
K 2422
12
K 34
K 3422
12
12
K 41
+ K 41
K 4122 + K 4122
11
K 42
K 4221
12
K 42
K 4222
11
K 43
K 4321
12
K 43
K 4322
11
11
K 44
+ K 44
K 4421 + K 4214
12
12
K 44
+ K 44
K 4422 + K 4422
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 026
Armado de la
Matriz de Rigidez
global
K11A + K11B
A
K
21
K 31B
K = A
B
+
K
K
41
41
UNC FCEFyN
CE II - 2014
C
4
K12A
K 22A
K13B
0
K 42A
0
K 33B + K 33C + K 33D
C
K 43B + K 43
K14A + K14B
K 24A
K 34B + K 34C
C
K 44A + K 44B + K 44
0
0
K 35D
0
0
K 53D
C
+ K 63D
K 63
0
C
K 64
0
K 55D
K 65D
1
0
0
2
K 36C + K 36D 3
C
K 46 4
K 56D 5
C
D 6
K 66 + K 66
27
UNC FCEFyN
CE II - 2014
Objetivos
Facilidad para la definicin del modelo y sus
datos.
Adecuada representacin de las caractersticas
elsticas del objeto estudiado.
Ausencia de problemas numricos
Reduccin del esfuerzo de clculo (tiempo de
proceso)
Facilidad para la interpretacin de los
resultados.
UNC FCEFyN
CE II - 2014
29
Recomendaciones
a) Utilizar elementos de forma regular.
b) La malla debe respetar fielmente los contornos
del objeto representado.
c) Las mallas deben densificarse en las zonas
donde se espera el mayor gradiente de
tensiones.
d) Las mallas deben densificarse gradualmente,
y no en forma brusca, evitndose que
elementos finitos de tamaos muy diferentes
compartan un mismo nudo.
UNC FCEFyN
CE II - 2014
30
Recomendaciones
e) Las mallas deberan ser regulares en el
sentido de que cada nudo sea compartido por
una cantidad similar de elementos.
f) En un proceso de refinamiento es
recomendable que las mallas mas densas
estn incluidas en las anteriores.
g) Las mallas densas son costosas y deben
evitarse. Solo densificar las mallas de manera
localizada y no general.
UNC FCEFyN
CE II - 2014
31
UNC FCEFyN
CE II - 2014
32
UNC FCEFyN
CE II - 2014
33
UNC FCEFyN
CE II - 2014
34
Malla 1
Malla 2
Malla 3
Malla 4
Malla 5
UNC FCEFyN
CE II - 2014
35
Convergencia
UNC FCEFyN
CE II - 2014
Condiciones de Consistencia
Completitud
Las funciones de desplazamiento deben ser
capaces de representar los desplazamientos del
cuerpo rgido y los estados de deformacin
constante.
Compatibilidad
Dentro de los elementos, y a travs de sus bordes,
las funciones de desplazamiento deben ser
continuas. No deben presentarse aberturas entre
elementos bajo condiciones de carga.
UNC FCEFyN
CE II - 2014
37
Condiciones de Convergencia
Consistencia
Debe haber consistencia entre el modelo fsico y
el modelo discreto que lo representa.
Estabilidad
A medida que se afina la malla la solucin tiende
en forma asinttica a un cierto valor, a condicin
de que no haya excesiva distorsin en los
elementos.
Exactitud
La aproximacin asinttica tiene a valores
correctos de la solucin del problema
UNC FCEFyN
CE II - 2014
38
Convergencia
Completitud
Estabilidad
Compatibilidad
Consistencia
Convergencia
Exactitud
UNC FCEFyN
CE II - 2014
39
Convergencia
f
UNC FCEFyN
CE II - 2014
40
Convergencia
UNC FCEFyN
CE II - 2014
Modelo
Elementos
Nudos
Grados de
Libertad
14.d
96
65
120
14.e
150
96
180
14.f
600
341
660
Convergencia
UNC FCEFyN
CE II - 2014
UNC FCEFyN
CE II - 2014
43
UNC FCEFyN
CE II - 2014
44
45
Elementos .....:
550
Nodos..........:
193
Grados Libertad: 1.000
Elementos: Vigas y TTC
UNC FCEFyN
CE II - 2014
46
47
Elementos .....:
234
Nodos..........: 1.488
Grados Libertad: 4.464
Elementos: Thick Shell 16 nodos
UNC FCEFyN
CE II - 2014
48
Elementos .....:
192
Nodos..........: 1.308
Grados Libertad: 7.800
Elementos: Thick Shell 16 nodos,
slidos de 20 nodos y
general shell de 8 nodos
UNC FCEFyN
CE II - 2014
49
UNC FCEFyN
CE II - 2014
50