Análisis por Elementos Finitos
Clase 1
Ing. Eduardo Orcés P.
Mayo 18/2011
Ing. Eduardo Orcés
2011-I
TEMAS
Objetivos de aprendizaje
Métodos de Análisis en el Proceso de
Diseño de Ingeniería
Desarrollo Histórico del MEF
Introducción al Método de Elementos
Finitos (MEF)
Concepto de Rigidez Estructural
Ing. Eduardo Orcés 2011-I
Método de Elementos Finitos
El Método de Elementos Finitos es una técnica muy poderosa de
análisis matemático („discretización‟) usada en mecánica de sólidos y
en otras disciplinas de ingeniería.
Ing Eduardo Orcés 2011-I
Objetivos de aprendizaje
Los objetivos de aprendizaje de este curso son:
o Entender la perspectiva, tópicos importantes, terminología
y pasos en la solución de problemas usando el MEF.
o Conocer el comportamiento y uso de los distintos tipos de
elementos tratados en este curso.
o Estar en capacidad de crear modelos apropiados de EF
para problemas dados.
o Poder interpretar y evaluar la calidad de los resultados.
o Estar consciente de las limitaciones del MEF, el cual es
una herramienta matemática.
Ing Eduardo Orcés 2011-I
Conceptos Básicos
• El Método de Elementos Finitos (MEF)
está basado en la idea de armar objetos
grandes y complejos a partir de bloques
constitutivos simples. En la vida diaria
existen múltiples ejemplos, así como en
la ingeniería.
• Juguetes „Lego‟
• Edificios
• Máquinas, etc.
Ing Eduardo Orcés 2011-I
Conceptos Básicos
• El Análisis por medio de Elementos
Finitos consiste en subdividir todos los
sistemas en componentes individuales
o “elementos” cuyo comportamiento es
de fácil comprensión y, a continuación,
reconstruir el sistema original a partir de
estos componentes.
Ing Eduardo Orcés 2011-I
Conceptos Básicos
• El modelo matemático “discretizado”
que se crea permite simular lo que
sucede con el sistema real al aplicarle
cargas, e identificar las áreas críticas.
• En esta imagen, las
áreas de color rojo
indican una posible
falla.
Ing Eduardo Orcés 2011-I
Conceptos Básicos
• Una vez que se ha modelado un sistema en
la computadora, es posible modificar su
diseño original y “optimizarlo” para obtener el
mejor desempeño posible del mismo.
• En la figura, se muestra una
biela de conexión, antes y
después de su optimización.
Ing Eduardo Orcés 2011-I
Métodos de Análisis en el
Proceso de Diseño de Ingeniería
PROCESO GENERALIZADO DE DISEÑO
ESPECIFICACIÓN DEL CLIENTE
Función
Costo
Mantenibilidad
ESPECIFICACIÓN TÉCNICA DE DISEÑO
Cargas
Forma geométrica
Material
Manufactura
ETAPAS DE DISEÑO
Diseño preliminar
Diseño de factibilidad
Diseño de detalle
Ing. Eduardo Orcés 2011-I
Métodos de Análisis en el
Proceso de Diseño de Ingeniería
• DISEÑO PRELIMINAR
- Análisis limitado (“cálculos manuales”)
- Esbozo de amplio espectro de
soluciones alternativas que cumplen
las especificaciones
- Costeo preliminar
• DISEÑO DE FACTIBILIDAD
- Análisis más avanzado (“numérico”)
- Restricción de rango de soluciones a
aquellas más probables de cumplir las
especificaciones
- Costeo más preciso
• DISEÑO DETALLADO
- Análisis intensivo, numérico o de otro tipo
- Diseño para manufactura
- Lista de partes
- Costeo detallado
Ing. Eduardo Orcés 2011-I
Métodos de Análisis
• Cálculos manuales:
- Análisis de componentes individuales
- Condiciones de borde simplificadas
- Análisis elástico simplificado
- etc.
• Ventajas: - Rápidos y de bajo costo
- Bastante aproximados si se tiene experiencia
- Buenos para estimar costos
• Desventajas: - Pueden ser mal aplicados por ingenieros sin
experiencia
- Tienden a ser muy conservadores
- Fallan en situaciones de cargas muy dinámicas
Ing. Eduardo Orcés 2011-I
Métodos de Análisis
• Métodos experimentales:
- Extensómetros eléctricos
- Fotoelasticidad
- etc.
• Ventajas:
- Ideales para estados complejos de esfuerzos ([Link]. Concentración
de esfuerzos)
- Se incluyen condiciones prácticas de carga
• Desventajas:
- Se requiere experiencia para fabricación de modelos e
interpretación de resultados
- Requieren instrumentación costosa
- No conveniente para análisis de rutina
Ing. Eduardo Orcés 2011-I
Métodos de Análisis
• Métodos numéricos:
• Diferencias finitas:
- Basado en la discretización de las ecuaciones diferenciales
- Aplicaciones extensas en Mecánica de fluidos, Transferencia de
calor, Flujo de metales, Impacto, etc.
• Ventajas: - Historia extensa de desarrollo
- Conocido por los ingenieros
• Desventajas: - Poca aplicación en Mecánica de sólidos
- Dificultades con las condiciones de borde
Ing. Eduardo Orcés 2011-I
Métodos de Análisis
• Métodos numéricos:
• Elementos de frontera:
- Basado en el teorema de Green de integrar sobre la frontera las
ecuaciones diferenciales del campo
- Discretización hecha solo en la frontera
• Ventajas: - Más fácil de crear la malla
- Puede tratar con exactitud problemas de dominios infinitos
- Ha producido excelentes resultados en problemas de
concentración de esfuerzos, etc.
• Desventajas: - Poco conocido por los ingenieros
- Pocos programas comerciales
- No se ha aplicado todavía a problemas no-lineales y dinámicos
Ing. Eduardo Orcés 2011-I
Métodos de Análisis
• Métodos numéricos:
• Elementos Finitos:
- Basado en el Principio de la Energía Potencial Mínima
- Discretización hecha sobre todo el campo
- Amplia aplicación en Diseño
• Ventajas: - Relativamente fácil de aprender
- Se conocen bien las bases matemáticas del método, lo que permite
estimar las condiciones de precisión y convergencia
- Programas comerciales disponibles en un amplio rango de problemas
lineales, no-lineales, dinámicos, etc.
• Desventajas: - Debido a la facilidad de su uso, los analistas inexpertos pueden
cometer errores costosos en la creación del modelo y en la
interpretación de resultados
- En el campo no-lineal, los resultados pueden ser dudosos y difíciles
de interpretar
Ing. Eduardo Orcés 2011-I
Métodos de Análisis
• Otros métodos numéricos:
• Mecánica de Fluidos Computacional
(CFD)
• Método de Volúmenes Finitos (MVF)
• Métodos sin malla
Ing. Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Dominio del flujo Dominio estructural
Dominios multi-celulares para modelos fluido-estructurales
Ing Eduardo Orcés 2011-I
Desarrollo Histórico del MEF
Arquímedes (250 A.C.): Rectificación del
circulo como el límite de polígonos regulares inscritos
1 2
Área de un triángulo : Si R sen i
2
N
1 2 2
Área del círculo: S N Si R Nsen R2 , N
i 1 2 N
donde N número total de triángulos ("elementos")
Ing Eduardo Orcés 2011-I
Desarrollo Histórico del MEF
R. Courant (1943)
Método parecido al de EF. Utilizaba funciones de
interpolación sobre regiones triangulares
Turner, Clough, Martin, y Topp (1956):
Análisis de estructuras elásticas complejas en aeronáutica usando
elementos finitos
Zienkiewicz y Cheung (1967):
El Método del Elemento Finito en la Mecánica Estructural
y del Medio Contínuo
Ing Eduardo Orcés 2011-I
Desarrollo Histórico del MEF
Década de los 70
Aplicaciones en computadores grandes (‘mainframes’)
Década de los 80
Microcomputadores, pre- y postprocesadores
Década de los 90
Análisis de sistemas estructurales grandes
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Subestructuras
• Las subestructuras fueron inventadas en la
industria aeroespacial, en los años 60.
Ing Eduardo Orcés 2011-I
• Subestructuras de primer nivel
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Rib spacing controlled by
results from preliminary design model
Ing Eduardo Orcés 2011-I
Superelementos
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
La subestructuración de múltiples niveles fue inventada en
la industria naviera noruega en los años 60.
Ing Eduardo Orcés 2011-I
Los cohetes multi-etapa se adaptan de manera natural al
uso de subestructuras
Ing Eduardo Orcés 2011-I
Aplicaciones en Ingeniería del MEF
Problemas de Probs. de Valores Problemas de
Área
Equilibrio Característicos Propagación
Estructuras Análisis estático
Propagación de
cerchas,pórticos, Frecuencias naturales y
mecánicas y placas,cascarones, pandeo de estructuras
[Link] transientes
y aleatorias
civiles etc.
Estructuras Análisis estático [Link], pandeo de Respuesta a cargas
aeronáuticas alas,fuselajes, etc. estructuras aeronáuticas aleatorias y transientes
[Link],
recipientes a [Link].y estabilidad
Diseño presión,engra- Problemas de fractura
de bielas,engranajes y bajo cargas dinámicas
mecánico najes,pistones, máquinas herramientas
etc.
Distribuciones de Flujo transiente de calor
Transferencia
temperaturas en --- en toberas,álabes,
de calor sólidos y fluidos motores de CI, etc.
Hidráulica/Ing. Flujos de superficie Frecs. y modos naturales de Flujos inestables y
libre,capa
de recursos estuarios,lagos y puertos. propagación de ondas en
límite,transónicos.
hídricos/Hidro- Movimiento de lìquidos en [Link]ón
Estructuras
recipientes. transiente en acuíferos
dinámica hidráulicas y presas
Muros,túneles,presas Frecs.y modos naturales de
Propagación de ondas en
Geotecnia mecánica de rocas y [Link]ón
suelos y rocas
suelos,fundaciones suelo-estructura
Esfuerzos en Cargas de impacto en
Ingeniería
huesos,dientes,etc --- crá[Link]ámica de
Biomédica Prótesis,válvulas,etc. estructuras anatómicas
2011-I
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
Prueba de impacto de un recipiente de presión
Ing Eduardo Orcés 2011-I
Programas de Computadora
La mayoría de los programas incluyen
todas, o algunas de las siguientes
capacidades:
• Preprocesadores (construir modelo de
EF, cargas y restricciones)
• Resolvedores FEA (ensamblar y resolver
el sistema de ecuaciones)
• Postprocesadores (ordenar y mostrar los
resultados)
Ing Eduardo Orcés 2011-I
Programas comerciales de
Computadora
FEA de propósito general:
• ANSYS, COSMOS, NASTRAN, VisualFEA
Pre y Postprocesadores:
• PATRAN, HyperMesh
Análisis no-lineal, dinámico, impacto
• SAP, ABAQUS, Dyna-3D, MARK
Paquetes completos CAD / CAM /CAE
• SDRC/I-DEAS, SolidWorks, Inventor
Ing Eduardo Orcés 2011-I
Programas comerciales de
Computadora
Software de propósito general:
• MATLAB, Mathematica, C, C++
Análisis de problemas de Campos (Fluidos,
Ondas, etc.)
• Fluent, FEMLAB, Software de BEM
Ing Eduardo Orcés 2011-I
Ejemplo: Pórtico modelado con
Matlab
Ing Eduardo Orcés 2011-I
Ejemplo: Flujo de calor en placa,
modelado en Matlab
Ing Eduardo Orcés 2011-I
Ejemplo: Análisis dinámico de un
pórtico, en Matlab
Ing Eduardo Orcés 2011-I
Ejemplo: Gancho modelado con
SolidWorks
Ing Eduardo Orcés 2011-I
Proceso de Idealizacion
Idealizacion
Ing Eduardo Orcés 2011-I
Descomposicion
Ing Eduardo Orcés 2011-I
Ensamblaje y Solucion
Ing Eduardo Orcés 2011-I
Conceptos Básicos del MEF
1. Discretización de la estructura:
Dividir en “elementos finitos”, y generar “malla”.
Decidir sobre número, tipo, material, tamaño y
disposición de los elementos.
Ing Eduardo Orcés 2011-I
Ing Eduardo Orcés 2011-I
2011-I
Ing Eduardo Orcés 2011-I
2011-I
Modelo Modelo Sólido
Vigas 1D Shell 2D 3D
2011-I
2. Selección de la función de interpolación
más apropiada:
Describir el comportamiento de las cantidades físicas en
cada elemento mediante funciones asumidas lo más
sencillas posible (p. ej., lineal).
Debe satisfacer requerimientos de convergencia y
compatibilidad (polinomio, generalmente).
Funciones de interpolación φ(x,y):
* Triángulo (3 nodos): φ = a1 + a2 x + a3 y
* Cuadrilátero (4 nodos): φ = a1 + a2 x + a3 y + a4 x y
* Cuadrilátero (8 nodos): φ = a1 + a2 x + a3 y + a4 x y +…
+ a5 x2 + a6y2 + a7x2y + a8xy2
Ing. Eduardo Orcés 2011-I
Ing. Eduardo Orcés 2011-I
3. Obtención de las matrices de rigidez k(e)
y vectores de carga p(e) de los elementos:
Obtener las relaciones constitutivas para cada elemento
(Ley de Hooke, Fourier, etc.). En Mecánica de Sólidos,
implica determinar relaciones entre fuerzas nodales
equivalentes y desplazamientos nodales. En problemas
térmicos, determinar relaciones entre flujos de calor
nodales equivalentes y temperaturas nodales.
Se aplica, ya sea un método “directo” (equilibrio,
conservación de energía y masa), ó algún principio
apropiado de “minimización” (Energía Potencial Mínima,
Galerkin, etc.)
2011-I
4. Armar la matriz de rigidez global K de la
estructura:
Conectar (ensamblar) los elementos en los nodos para
sumar apropiadamente las contribuciones de cada
elemento, y así obtener un sistema aproximado de
ecuaciones para toda la estructura, que en forma
matricial es de la forma
KD=F
Esta ecuación matricial implica equilibrio (conservación)
global para toda la estructura.
K = Σ k (e) = Matriz de rigidez global
D = Vector de desplazamientos (temperaturas) nodales
F = Σ p (e) = Vector de fuerzas (flujos) nodales
Ing. Eduardo Orcés 2011-I
5. Aplicar las condiciones de borde
apropiadas:
En Mecánica de Sólidos, especificar C.B. esenciales o
de Dirichlet (desplazamientos conocidos, [Link]. cero),
y/o C.B. naturales o de Neumann (momentos
conocidos, [Link]. cero).
En Transferencia de Calor, temperaturas y/o flujos
conocidos.
Ing. Eduardo Orcés 2011-I
6. Resolver un sistema de ecuaciones
simultáneas para hallar los desplazamientos
(temperaturas):
En problemas lineales, aplicar Eliminación Gaussiana, o
equivalente. En problemas no-lineales, usar iteración,
modificando en cada paso la matriz K, y/o el vector F.
Ing. Eduardo Orcés 2011-I
7. Cálculo de esfuerzos y deformaciones
(flujos de calor) en los elementos:
Para cada elemento, calcular:
Deformación unitaria: ε = B d
Esfuerzos (Flujos): σ = E ε = E B d
donde:
d = Vector de desplazamientos (temperaturas)
nodales
B = Matriz de gradientes
E = Matriz de propiedades del material (E,μ; k)
Ing. Eduardo Orcés 2011-I
Matrices de Rigidez
• Matrices elementales:
• Resorte:
P1 k (u1 u 2 )
P2 k (u 2 u1 )
En forma matricial:
k k u1 P1
k (e) d P (e)
k k u2 P2
La matriz de rigidez k (e) es “singular” (no puede ser invertida) debido a que
el resorte “flota en el espacio”.
Ing. Eduardo Orcés 2011-I
• Barra en tensión:
Igual que el resorte con k = AE / L.
• Conducción de calor (Corriente eléctrica) en barra:
k 1 1 T1 q1
L 1 1 T2 q2
1 1 1 V1 I1
R 1 1 V2 I2
Ing. Eduardo Orcés 2011-I
• Matriz global (ensamblada):
• Resortes en serie:
• Equilibrio: F1 = P1 , F2 = P2 (1) + P2 (2) , F3 = P3 (A)
• Compatibilidad: d1 (1) = u1 , d2 (1) = u2
d1 (2) = u2 , d2 (2) = u3 (B)
• Relación Fuerza / Deformación:
P1 = - P2 (1) = k1 (u1 – u2 )
P3 = - P2 (2) = k2 (u3 – u2 ) (C)
Ing. Eduardo Orcés 2011-I
• (A), (B), (C) F1 = P1 = k1 (u1 – u2 ) = k1 u1 – k1 u2
F2 = P2 (1) + P2 (2) = k1 (u2 – u1 ) + k2 (u2 – u3 )
= - k1 u1 + (k1 + k2 )u2 – k2 u3
F3 = P3 = k2 (u3 – u2 ) = - k2 u2 + k2 u3
• En forma matricial:
F1 k1 k1 0 u1
F2 k1 k1 k 2 k2 u2 F K D
F3 0 k2 k2 u3
• Ensamblaje directo de la Matriz de Rigidez Global K :
Observando el resultado anterior, se concluye que se puede obtener
directamente la matriz global K si se superponen directamente las matrices
elementales k (e), tomando en cuenta a qué coordenada nodal corresponde
cada columna / fila de k (e) (conectividad ) .
Ing. Eduardo Orcés 2011-I
u1 u2 GDL globales u1 u2 u3
k1 k1 u1 k1 k1 0 u1
(1)
k Expansión al K (1) k1 k1 0 u2
k1 k1 u2 tamaño de K
0 0 0 u3
u2 u3 GDL globales u1 u2 u3
k2 k2 0 0 0
( 2) u2 u1
k K ( 2)
0 k2 k2 u2
k2 k2 u3 u3
0 k2 k2
u1 u2 u3
k1 k1 0 u1
(e)
K K k1 k1 k2 k2 u2
e
0 k2 k2 u3
Ing. Eduardo Orcés 2011-I
Observaciones:
• Las matrices K y k(e) son simétricas
• El tamaño de k(e) depende del número de grados de libertad por elemento
(2 para un resorte, 4 para una viga, etc.)
• El tamaño de K depende del número total de grados de libertad globales (3
en el ejemplo) puede ser muy grande!
Ing. Eduardo Orcés 2011-I
Tareas
Leer los siguientes capítulos del texto guía,
Introducción al estudio del Elemento Finito en
Ingeniería, Chandrupatla/Belegundu, 2ª ed., Edit.
Pearson:
- Cap. 1, Conceptos Fundamentales
- Cap. 2, Álgebra de Matrices y Eliminación de
Gauss
Ing. Eduardo Orcés 2011-I
