Cursofem 0

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

v1

5 2
4
1 v2 Introducción al Método
S
de los Elementos Finitos

Alberto Cardona, Víctor Fachinotti


Cimec (UNL/Conicet), Santa Fe, Argentina
19/08/2014
Introducción
Modelos matemáticos en ciencia e ingeniería Ecuaciones algebraicas,
diferenciales o integrales
El desarrollo de las computadoras permitió usar estos modelos para resolver
problemas prácticos. Se pueden simular y resolver sistemas altamente
complicados en ciencia e ingeniería.
Permiten:
1. Reducir la necesidad de experiencias con modelos y prototipos (caras y
lentas).
2. Comparar fácilmente distintas alternativas de diseño para llegar al óptimo
ingenieril.

Disciplinas relacionadas:
1. CAD: Computer Aided Design
2. CAE: Computer Aided Engineering
3. CAM: Computer Aided Manufacturing

Introducción al Método de los Elementos Finitos 2


Introducción (cont)

Soluciones analíticas
Sólo casos simples

Modelos matemáticos
Casos prácticos, en computadora

Soluciones numéricas

Método de los Elementos Finitos (MEF) : técnica general para hallar soluciones
numéricas de sistemas de ecuaciones diferenciales e integrales.

Origen: ingeniería estructural, años 50/60, para solución de ecuaciones


diferenciales en derivadas parciales en elasticidad.

Su aplicación se generalizó, integrado a sistemas de CAD/CAE.

Introducción al Método de los Elementos Finitos 3


Introducción (cont)
Aplicaciones del MEF:
• Ingeniería estructural
• Resistencia de materiales
• Mecánica de fluidos
• Ingeniería nuclear
• Electromagnetismo
• Campos eléctricos
• Propagación de ondas
• Conducción del calor
• Procesos de convección – difusión
• Ingeniería de petróleo
• Procesos de reacción – difusión
• ........

Introducción al Método de los Elementos Finitos 4


Método de diferencias finitas vs. MEF
•Idea básica de un método numérico para resolver ecuaciones diferenciales:

Continuo con Discretización Problema discreto


infinitos GDL c/ # GDL finito

•El problema discreto puede resolverse en una computadora !!!

Solución discreta,
aproximada, con N GDL

•Diferencias finitas: método numérico clásico para resolver ecuaciones


diferenciales donde el problema discreto se obtiene reemplazando:
u un1  un

x x
Introducción al Método de los Elementos Finitos 5
Discretización en el MEF
1. Reformulación de la ecuación diferencial en un problema variacional
equivalente.
Ej: en ecuaciones elípticas, en casos simples, toma la forma de problema de
minimización:

(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 varían en forma
continua (ej.. desplazamiento de un cuerpo elástico, temperatura, etc.)
F(v) es la energía total asociada a v.
(M) es una caracterización equivalente de la solución de la ecuación diferencial
como aquella función en V que minimiza la energía total del sistema
considerado.

Introducción al Método de los Elementos Finitos 6


Discretización en el MEF (cont)
2. En general, la dimensión de V es infinita (las funciones de V no pueden
expresarse a través de un número finito de parámetros). Luego, (M) no puede
resolverse en forma analítica.

Para hallar una solución, la idea del MEF es reemplazar V por un conjunto Vh
de funciones simples que dependen de un número finito de parámetros:
(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 aproximación suficientemente buena de u , la solución de (M).
3. Usualmente elegimos:
Vh  V / vh Vh  vh V
En este caso, (M h ) corresponde al método clásico de Ritz-Galerkin
Introducción al Método de los Elementos Finitos 7
Discretización en el MEF (cont)
4. Característica particular del MEF: funciones de Vh son funciones polinomiales por tramos

Veremos más adelante que se pueden


hacer formulaciones variacionales más
generales que el problema de
minimización; ej.: métodos de Galerkin

Pasos para resolver por MEF:


1. Formulación variacional del problema
2. Discretización por MEF: construcción del espacio dimensional finito Vh
3. Solución del problema discreto
4. Implementación del método en computadora: programación

Introducción al Método de los Elementos Finitos 8


Comentarios
• Existen distintas formulaciones variacionales dependiendo, por ejemplo, de
la elección de las variables independientes
• La elección del subespacio de aproximación de dimensión finita Vh , o sea la
elección del tipo de elemento finito, está influenciada por:
– Formulación variacional
– Requerimientos de precisión
– Propiedades de regularidad de la solución exacta
– ...
• Para resolver el problema discreto necesitamos algoritmos de optimización
y/o algoritmos de solución de grandes sistemas de ecs. algebraicas lineales o
no lineales

Introducción al Método de los Elementos Finitos 9


Ventajas del MEF respecto de DF
– Tratamiento de geometrías complicadas
– Condiciones de borde generales
– Propiedades materiales no lineales o variables
– La estructura clara del método permite crear códigos multipropósito
generales
– Fundamentos teóricos sólidos
– Confiabilidad
– Posibilidad de estimación de error

Introducción al Método de los Elementos Finitos 10


Modelo para cálculo de tensiones en una brida

11
Verificación estructural de contenedor de material radiactivo

Introducción al Método de los Elementos Finitos 12


Verificación térmica de contenedor de material radiactivo

Nodal temperature [K]


T=1800 s
1100.

Verificación de daños por incendio 1000.

900.
Pt. A
800.
Pt. A
Pt. B
700.
Pt. C

600.
Pt. B
500.

Pt. C
400.

1000. 2000. 3000. 4000. 5000. 6000. 7000. 8000.


Time [s]

Evolución de
temperaturas
durante un incendio

Objetivo: Verificar el cumplimiento de la


reglamentación para contenedores de transporte.

Introducción al Método de los Elementos Finitos


Simulación de SMD (shaped metal deposition)
Temperaturas

Modelo CAD
Macrografía

Modelo FEM

Proyecto RAPOLAC Introducción al Método de los Elementos Finitos 14


(Unión Europea)
Simulation of shaped metal deposition process

Model t = 500 s
Validation

t = 55 s

Proyecto RAPOLAC
t = 15 s
(Unión Europea)

Introducción al Método de los Elementos Finitos 15


Simulación de SMD (shaped metal deposition)

Calculo tensiones RAPOLAC

16
Introducción al Método de los Elementos Finitos
Análisis de colada continua

Crack micrography

Objetivo: análisis termomecánico acoplado de la solidificación para determinación de


temperaturas, deformaciones, tensiones, macrosegregación

Comitente: Instituto Argentino de Siderurgia

Introducción al Método de los Elementos Finitos 17


Análisis de colada continua
Two spraying regions (original) Three spraying regions (proposed)

Temperature [ºC]

Plastic strain

Cross sections 10.5 m below meniscus


Introducción al Método de los Elementos Finitos 18
Simulación de tratamiento térmico de rodillos de laminación

Análisis de temperaturas, tensiones y deformaciones plásticas

Objetivo: Análisis fisuración de rodillos de laminación tricapa de fundición blanca de alto Cr.
Comitente: Fundicion San Cayetano – Burzaco – Prov. Buenos Aires
Introducción al Método de los Elementos Finitos 19
Reducción del depósito de polvo en una instalación de sinterizado

Precipitador electrostático y conducto de


entrada planta SINTER
Patrón de flujo en la manga de entrada
al precipitador
Objetivo: Optimizar el patrón de flujo en el
conducto principal del sistema, para
minimizar la deposición de partículas sólidas. Deflector diseñado
para reducir la
Comitente: Ternium Siderar, San Nicolás
(Bs.As.) deposición

Introducción al Método de los Elementos Finitos 20


Sistema de aireación en standpipe de equipo de cracking

La aireación es necesaria para


mantener los sólidos fluidizados.
Una aireación insuficiente produce
un flujo irregular y poca diferencia
de presión, en tanto una
sobreaireación puede conducir a
un flujo intermitente y reducida
diferencia de presión.

Objetivo: Determinar la influencia de los


diferentes puntos de aireación en la fluidización de
la mezcla gas-catalizador, para mejorar la
eficiencia del equipo de cracking.
Comitente: Repsol YPF, Ensenada (Bs.As.)

Introducción al Método de los Elementos Finitos 21


Análisis de interacción fluido estructura en régimen compresible

Encendido de un
motor cohete.
Deformación de la
estructura debido a
Propagación de ondas de choque en un
la onda de choque
contenedor de gas detenido súbitamente

Determinación del número de Mach crítico para


flujo sobre una placa plana

Objetivo: Desarrollar un código con


acoplamiento bidireccional entre el análisis
estructural y el análisis fluidodinámico.
Comitente: Open Engineering - ESTEC/ESA
Introducción al Método de los Elementos Finitos 22
Programas

• Usaremos un programa en Matlab / Octave, que distribuimos a través de la


página del curso, y sobre el cual realizarán programación

• Además, haremos algunos prácticos 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)

Introducción al Método de los Elementos Finitos 23

También podría gustarte