v1
4
1
v2
Introduccin al Mtodo
de los Elementos Finitos
Alberto Cardona, Vctor Fachinotti
Cimec (UNL/Conicet), Santa Fe, Argentina
19/08/2014
Introduccin
Modelos matemticos en ciencia e ingeniera
diferenciales o integrales
Ecuaciones algebraicas,
El desarrollo de las computadoras permiti usar estos modelos para resolver
problemas prcticos. Se pueden simular y resolver sistemas altamente
complicados en ciencia e ingeniera.
Permiten:
1. Reducir la necesidad de experiencias con modelos y prototipos (caras y
lentas).
2. Comparar fcilmente distintas alternativas de diseo para llegar al ptimo
ingenieril.
Disciplinas relacionadas:
1. CAD: Computer Aided Design
2. CAE: Computer Aided Engineering
3. CAM: Computer Aided Manufacturing
Introduccin al Mtodo de los Elementos Finitos
Introduccin (cont)
Soluciones analticas
Slo casos simples
Modelos matemticos
Casos prcticos, en computadora
Soluciones numricas
Mtodo de los Elementos Finitos (MEF) : tcnica general para hallar soluciones
numricas de sistemas de ecuaciones diferenciales e integrales.
Origen: ingeniera estructural, aos 50/60, para solucin de ecuaciones
diferenciales en derivadas parciales en elasticidad.
Su aplicacin se generaliz, integrado a sistemas de CAD/CAE.
Introduccin al Mtodo de los Elementos Finitos
Introduccin (cont)
Aplicaciones del MEF:
Ingeniera estructural
Resistencia de materiales
Mecnica de fluidos
Ingeniera nuclear
Electromagnetismo
Campos elctricos
Propagacin de ondas
Conduccin del calor
Procesos de conveccin difusin
Ingeniera de petrleo
Procesos de reaccin difusin
........
Introduccin al Mtodo de los Elementos Finitos
Mtodo de diferencias finitas vs. MEF
Idea bsica de un mtodo numrico para resolver ecuaciones diferenciales:
Continuo con
infinitos GDL
Discretizacin
Problema discreto
c/ # GDL finito
El problema discreto puede resolverse en una computadora !!!
Solucin discreta,
aproximada, con N GDL
Diferencias finitas: mtodo numrico clsico para resolver ecuaciones
diferenciales donde el problema discreto se obtiene reemplazando:
un1 un
u
x
x
Introduccin al Mtodo de los Elementos Finitos
Discretizacin en el MEF
1. Reformulacin de la ecuacin diferencial en un problema variacional
equivalente.
Ej: en ecuaciones elpticas, en casos simples, toma la forma de problema de
minimizacin:
(M)
Hallar u V
F (u) F (v) v V
donde:
V: conjunto de funciones admisibles
F :V
es un funcional
F (v)
v V
las funciones v V a menudo representan cantidades que varan en forma
continua (ej.. desplazamiento de un cuerpo elstico, temperatura, etc.)
F(v) es la energa total asociada a v.
(M) es una caracterizacin equivalente de la solucin de la ecuacin diferencial
como aquella funcin en V que minimiza la energa total del sistema
considerado.
Introduccin al Mtodo de los Elementos Finitos
Discretizacin en el MEF (cont)
2. En general, la dimensin de V es infinita (las funciones de V no pueden
expresarse a travs de un nmero finito de parmetros). Luego, (M) no puede
resolverse en forma analtica.
Para hallar una solucin, la idea del MEF es reemplazar V por un conjunto Vh
de funciones simples que dependen de un nmero finito de parmetros:
(Mh )
Hallar uh Vh / F (uh ) F (v) v Vh
Este problema es equivalente a un sistema de ecs. algebraicas. Se espera que uh
sea una aproximacin suficientemente buena de u , la solucin de (M).
3. Usualmente elegimos:
Vh V
/ vh Vh vh V
En este caso, (M h ) corresponde al mtodo clsico de Ritz-Galerkin
Introduccin al Mtodo de los Elementos Finitos
Discretizacin en el MEF (cont)
4. Caracterstica particular del MEF: funciones de Vh son funciones polinomiales por tramos
Veremos ms adelante que se pueden
hacer formulaciones variacionales ms
generales que el problema de
minimizacin; ej.: mtodos de Galerkin
Pasos para resolver por MEF:
1.
Formulacin variacional del problema
2.
Discretizacin por MEF: construccin del espacio dimensional finito Vh
3.
Solucin del problema discreto
4.
Implementacin del mtodo en computadora: programacin
Introduccin al Mtodo de los Elementos Finitos
Comentarios
Existen distintas formulaciones variacionales dependiendo, por ejemplo, de
la eleccin de las variables independientes
La eleccin del subespacio de aproximacin de dimensin finita Vh , o sea la
eleccin del tipo de elemento finito, est influenciada por:
Formulacin variacional
Requerimientos de precisin
Propiedades de regularidad de la solucin exacta
...
Para resolver el problema discreto necesitamos algoritmos de optimizacin
y/o algoritmos de solucin de grandes sistemas de ecs. algebraicas lineales o
no lineales
Introduccin al Mtodo de los Elementos Finitos
Ventajas del MEF respecto de DF
Tratamiento de geometras complicadas
Condiciones de borde generales
Propiedades materiales no lineales o variables
La estructura clara del mtodo permite crear cdigos multipropsito
generales
Fundamentos tericos slidos
Confiabilidad
Posibilidad de estimacin de error
Introduccin al Mtodo de los Elementos Finitos
10
Modelo para clculo de tensiones en una brida
11
Verificacin estructural de contenedor de material radiactivo
Introduccin al Mtodo de los Elementos Finitos
12
Verificacin trmica de contenedor de material radiactivo
Nodal temperature [K]
T=1800 s
1100.
Verificacin de daos por incendio
1000.
900.
Pt. A
800.
Pt. A
Pt. B
700.
Pt. C
600.
Pt. B
500.
Pt. C
400.
1000.
2000.
3000.
Evolucin de
temperaturas
durante un incendio
Objetivo: Verificar el cumplimiento de la
reglamentacin para contenedores de transporte.
Introduccin al Mtodo de los Elementos Finitos
4000.
5000.
6000.
7000.
8000.
Time [s]
Simulacin de SMD (shaped metal deposition)
Temperaturas
Modelo CAD
Macrografa
Modelo FEM
Proyecto RAPOLAC
(Unin Europea)
Introduccin al Mtodo de los Elementos Finitos
14
Simulation of shaped metal deposition process
Model
Validation
t = 500 s
t = 55 s
Proyecto RAPOLAC
(Unin Europea)
t = 15 s
Introduccin al Mtodo de los Elementos Finitos
15
Simulacin de SMD (shaped metal deposition)
Calculo tensiones RAPOLAC
Introduccin al Mtodo de los Elementos Finitos
16
Anlisis de colada continua
Crack micrography
Objetivo: anlisis termomecnico acoplado de la solidificacin para determinacin de
temperaturas, deformaciones, tensiones, macrosegregacin
Comitente: Instituto Argentino de Siderurgia
Introduccin al Mtodo de los Elementos Finitos
17
Anlisis de colada continua
Two spraying regions (original)
Three spraying regions (proposed)
Temperature [C]
Plastic strain
Cross sections 10.5 m below meniscus
Introduccin al Mtodo de los Elementos Finitos
18
Simulacin de tratamiento trmico de rodillos de laminacin
Anlisis de temperaturas, tensiones y deformaciones plsticas
Objetivo: Anlisis fisuracin de rodillos de laminacin tricapa de fundicin blanca de alto Cr.
Comitente: Fundicion San Cayetano Burzaco Prov. Buenos Aires
Introduccin al Mtodo de los Elementos Finitos
19
Reduccin del depsito de polvo en una instalacin de sinterizado
Precipitador electrosttico y conducto de
entrada planta SINTER
Patrn de flujo en la manga de entrada
al precipitador
Objetivo: Optimizar el patrn de flujo en el
conducto principal del sistema, para
minimizar la deposicin de partculas slidas.
Comitente: Ternium Siderar, San Nicols
([Link].)
Introduccin al Mtodo de los Elementos Finitos
Deflector diseado
para reducir la
deposicin
20
Sistema de aireacin en standpipe de equipo de cracking
La aireacin es necesaria para
mantener los slidos fluidizados.
Una aireacin insuficiente produce
un flujo irregular y poca diferencia
de presin, en tanto una
sobreaireacin puede conducir a
un flujo intermitente y reducida
diferencia de presin.
Objetivo: Determinar la influencia de los
diferentes puntos de aireacin en la fluidizacin de
la mezcla gas-catalizador, para mejorar la
eficiencia del equipo de cracking.
Comitente: Repsol YPF, Ensenada ([Link].)
Introduccin al Mtodo de los Elementos Finitos
21
Anlisis de interaccin fluido estructura en rgimen compresible
Encendido de un
motor cohete.
Deformacin de la
estructura debido a
la onda de choque
Propagacin de ondas de choque en un
contenedor de gas detenido sbitamente
Determinacin del nmero de Mach crtico para
flujo sobre una placa plana
Objetivo: Desarrollar un cdigo con
acoplamiento bidireccional entre el anlisis
estructural y el anlisis fluidodinmico.
Comitente: Open Engineering - ESTEC/ESA
Introduccin al Mtodo de los Elementos Finitos
22
Programas
Usaremos un programa en Matlab / Octave, que distribuimos a travs de la
pgina del curso, y sobre el cual realizarn programacin
Adems, haremos algunos prcticos con un programa comercial. Si disponen
de un programa, pueden usarlo. De lo contrario, pueden procurarse el
siguiente:
Samcef-Student (versiones disponibles en internet o de lo contrario, para
uso en red CIMEC)
Introduccin al Mtodo de los Elementos Finitos
23
