Teórico Método de Elementos Finitos

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 50

Clculo Estructural II

Dpto. Estructuras FCEFyN UNC


2014

Introduccin al
Mtodo de los Elementos Finitos
(FEM)
Clculo Estructural II
FCEFyN Dpto. Estructuras
Juan Gir - 2014

Orgenes del Mtodo de los Elementos Finitos


Collar, Duncan y Frazer plantearon a la
aeroelasticidad en forma matricial (1934-1938).
Propusieron la formulacin y terminologa para la
discretizacin del continuo utilizada actualmente.
Argyris sistematiz el concepto del ensamble de
sistemas estructurales a partir de componentes
elementales y algebra matricial (1954).
Turner desarroll un modelo de un ala delta para
su estudio areoelastico con elementos tringulos
(1953 a 1956).
Clough denomin FEM al nuevo mtodo (1960).
UNC FCEFyN
CE II - 2014

Orgenes del Mtodo de los Elementos Finitos


Zienkiewicz impuls la difusin del mtodo a
travs de la publicacin de su reconocido libro
(1960).
Wilson y Bathe contribuyeron a la utilizacin
practica del mtodo con la difusin, sin costo, de
excelentes programas de calculo, tales como el
SAP Structural Analysis Program (1965).
La NASA inici el proyecto que condujo al
NASTRAN (1965).
Westinghouse envi al mercado el programa
ANSYS (1969).
UNC FCEFyN
CE II - 2014

Campos de aplicacin de los Elementos Finitos

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

Desarrollo de modelos discretos


Idealizacin

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

Desarrollo de modelos discretos

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

Ai ( x1j x1 )a2i ( x2j x2 )a1i


i ( x1 , x2=
) =
2A
A
Aj ( x1k x1 )a2j ( x2k x2 )a1j
j ( x1 , x=
=
2)
2A
A
Ak ( x1i x1 )a2k ( x2i x2 )a1k
k ( x1 , x=
=
2)
2A
A

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

Funcin de aproximacin lineal

CE II - 2014

a1i

x2

K
(i)
(j)

AJ

P
AK

a2 i

Ai
J

(k)

x1

u1 ( x1 , x2 ) =u1i i ( x1 , x2 ) + u1j j ( x1 , x2 ) + u1k k ( x1 , x2 )


u2 ( x1 , x2 ) =u2i i ( x1 , x2 ) + u2j j ( x1 , x2 ) + u2k k ( x1 , x2 )
K

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

u1m a1m u2m a2m


12 =
4A
UNC FCEFyN
CE II - 2014

13

Energa de Deformacin de un elemento


1
( 1111 + 22 22 + 2 1212 ) h.dA

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 )

Potencial de las cargas exteriores


Un =
( F1 u1 + F2 u2 ) dA +
( f1 u1 + f 2 u2 ) dS
A

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

Mnima Energa Potencial Total

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

Matriz de Rigidez del Elemento Tringulo


kii11
21
kii
k 11
ji
21
k ji
k 11
ki21
kki

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

kik12 u1i F1i



kik22 u2i F2i
j
j
k 12
F1
A
jk u1
=

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

Ejemplo del clculo de la Matriz de Rigidez


11
k33
21
k33
11
k43
A
K = 21
k43
k 11
53
21
k53

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

Tipos de elementos de dos y tres dimensiones

Tringulos de Tensin Lineal y Cuadrtica


Cuadrilteros
Tetraedros de Tensin Constante
Prismas rectangulares
Elementos placa y slidos regulares
Elementos isoparamtricos
Elementos axil simtricos.
Bandas finitas
Elementos de compuestos
Elementos prismticos rectos y curvos

UNC FCEFyN
CE II - 2014

19

Complejidad de elementos y grados de libertad


K

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

Otros elementos finitos

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

Criterios para el ensamble


de elementos generacin
de mallas

UNC FCEFyN
CE II - 2014

Criterios para la generacin de mallas

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

Criterios para la generacin de mallas

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

Criterios para la generacin de mallas

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

Criterios para la generacin de mallas

UNC FCEFyN
CE II - 2014

32

Criterios para la generacin de mallas

UNC FCEFyN
CE II - 2014

33

Criterios para la generacin de mallas

UNC FCEFyN
CE II - 2014

34

Criterios para la generacin de mallas

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

Implementacin: Contexto Sistema MEF

UNC FCEFyN
CE II - 2014

43

Implementacin: Ncleo de Calculo

UNC FCEFyN
CE II - 2014

44

Sistemas de clculo ms difundidos


SAP (Structural Analysis Program), SAP4, SAP 6,
NONSAP, SAP7 y SAP 2000, desde 1970.
STRUDL (STRUctural Design Language), 1967.
NASTRAN (Nasa STRuctural ANalysis), 1964.
ANSYS (Analysis System), 1969.
ASKA, 1970
NISA (Numerically Integrated elements of
System Analysis), 1978.
ABAQUS, 1978.
CATIA (Computer Aided Three-Dimensional
Interactive), 1980.
SOLIDWORKS, 1995.
UNC FCEFyN
CE II - 2014

45

Ejemplo: chasis de auto de frmula 2

Elementos .....:
550
Nodos..........:
193
Grados Libertad: 1.000
Elementos: Vigas y TTC
UNC FCEFyN
CE II - 2014

46

Ejemplo: chasis de un camin

Elementos .....: 1.400


Nodos..........: 4.400
Grados Libertad: 24.000
Elementos: General shell 8 nodos
UNC FCEFyN
CE II - 2014

47

Ejemplo: anlisis de mltiple de escape

Elementos .....:
234
Nodos..........: 1.488
Grados Libertad: 4.464
Elementos: Thick Shell 16 nodos
UNC FCEFyN
CE II - 2014

48

Ejemplo: modelo de carcaza de diferencial

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

También podría gustarte