Tema 1.

Introducción

UNIVERSIDAD SANTO TOMAS- SECCIONAL

BUCARAMANGA
FACULTAD DE INGENIERÍA MECATRÓNICA

MÉTODO DE LOS ELEMENTOS FINITOS

Y APLICACIONES

TEMA
TEMA 11

INTRODUCCIÓN
 Tema 1. Introducción

ÍNDICE

Ejemplo de método numérico.

Ejemplo y definición de EDP.
Ejemplos y clasificación de las EDP.
1.- Métodos numéricos para
resolución de EDP
2.- Diferentes enfoques para plantear
el MEF
3.- Cronología
 Tema 1. Introducción

Ejemplo de método numérico

Aproximar la función cos(x) en el punto x=0 a una función

polinomial de la forma P(x)=C_0+C_1*x+C_2*x^2.
1) Suponer que en en cero la aproximación tiene el mismo valor que
la función.

2) Suponer que en en cero la derivada de la aproximación tiene el

mismo valor que la función.

3) Suponer que en en cero la curvatura (senguda derivada) tiene el

mismo valor que la función.
 Tema 1. Introducción

En Matlab…
 Tema 1. Introducción

En Matlab…
x=linspace(-5,5,100); %El dominio
y=cos(x); %La función real
p=1+0*x-0.5*x.^2;
p2=1+0*x-0.5*x.^2+(1/24)*x.^4;%función
aproximada
%Gráfica
f=figure
plot(x,y)
title("Función coseno y su
aproximación")
hold on
plot(x,p)
hold on
plot(x,p2)
legend({"Cos(x)","Función
aproximada en x=0", "p2"},
'location','southwest')
ylim([-2 2])
grid on
hold off
 Tema 1. Introducción

Expansión de Taylor y notación O(h^p).

P(x) se denomina polinomio de Taylor.

Ojo:
• Aparece el factorial
• Al agregar más términos no se afectan los anteriores.
• Si la aproximación se hace en otro lugar “a” diferente a cero
entonces en el polinomio x se transforma en (x-a)
Ejercicio: expresar la segunda derivada f´´(x) en términos de
sus derivadas anteriores O(h^p). Escribir las series de Taylor
en x+h y en x-h
 Tema 1. Introducción

Ecuación de calor en una barra

 Tema 1. Introducción

Introducción a las EDP

En la E.D.P. lineal si u representa la

Clasificación de variable independiente y X, Y las
las Ecuaciones
diferenciales variables independientes, entonces la
forma general es
Según la
Lineales y no
condición o tipo lineales
Según el orden
de problema  

Valor Inicial Primer orden

Valores en la
Segundo orden
En donde A,B,C,D,E son constantes
frontera
reales, si la E.D.P es homogénea
(G(x,y)=0) entonces:
Tercer orden,
Orden n, etc…
Hiperbólica  

Parabólica  

Elíptica  
 Tema 1. Introducción

Ejercicio: clasificar las siguientes ecuaciones diferenciales

 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

Problema de contorno:
Es un problema gobernado por:
•ecuaciones diferenciales o integrales en un dominio, y
•condiciones de contorno en la frontera del dominio

Ejemplo: problema elástico z, w

 x  xy  zx
   bx  0 P
x y z
dominio

t
 xy  y  yz
   by  0 σ  Dε ε  Lu
x y z
 zx  yz  z
   bz  0 y, v
x y z b

t x  l x  m xy  n zx
frontera

x, u
t y  l  xy  m y  n yz ~
uu
t z  l  zx  m yz  n z
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

Problema de contorno:
Es un problema gobernado por:
•ecuaciones diferenciales o integrales en un dominio, y
•condiciones de contorno en la frontera del dominio

Generalmente no se dispone de una solución analítica exacta

Soluciones numéricas aproximadas (errores)
•Diferencias finitas
•Funciones de prueba
•FEM
•BEM
•Métodos sin malla
•Volúmenes finitos…

Problemas físicos
•Discretos: Sol. (analítica o numérica) sencilla
•Continuos: Sol. más compleja, analítica o numérica (mediante
discretización)
 Tema 1. Introducción

Métodos numéricos para resolución de E.D.P

 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

1.1.- Diferencias finitas

f f f
 i 1, j i , j
x ij h

f f f
 i , j 1 i , j
y ij k

2 f f i 1, j 2 f i , j  f i 1, j

x 2 ij h2

•Es difícil tratar con geometrías o condiciones de contorno complejas.

•Posible mal condicionamiento numérico.
•El algoritmo de resolución depende de la ecuación (es difícil de generalizar).
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

1.1.- Diferencias finitas

Problema: Resolver la ecuación de calor de una región cuadrada limitada por
x=2 y y=2. con las condiciones de contorno dadas por: (con tamaño de paso
h=2/3)
 

 
 

Solución: sumar la expresión de la expansión de Taylor de la segunda

derivada hallada en el punto anterior.

 
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

1.1.- Diferencias finitas

RECORDAR TÉRMINO Oh
SERIES DE TAYLOR (término de
error)
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

1.2.- Método de las funciones de prueba

Aproximación para u ( x, y ) uˆ ( x, y ) a1  a2 xa3 y a4 x 2  a5 xy ...

todo el dominio:

Minimización de un funcional a1 , a2 , a3 , a4 , a5 , ...

o residuos ponderados:

DESVENTAJAS:
•Solución única para todo el subdominio
 puede requerir un número excesivo de términos
 grado polinómico elevado,
 posible mal condicionamiento numérico.
•Coeficientes ai no son fácilmente interpretables.
•Difícil satisfacer condiciones de contorno generales.
•Matriz de coeficientes llena.

ex, series de Taylor

 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

1.2.2.- Residuos ponderados

u campo incógnita exacto
Ecuación diferencial L(u ) f 0 L operador diferencial
A Dominio
f función definida en A
Residuo R(u ) L(u ) f

Forma integral W (u )   R (u )dA Y función de ponderación

A

Residuo y forma R(uˆ ) L(uˆ ) f û campo incógnita aproximado

integral para û W (uˆ )  R(uˆ )dA
A

Forma integral para W j (uˆ )  j R (uˆ )dAW j (a1 ,..,am ) a1,…am coeficientes incógnita
selección de Y’s A
j 1,..,m Yj funciones de ponderación

Anulación de forma W j ( a1 ,..,am )0

j 1,..,m
integral Yj
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

1.2.2.- Residuos ponderados

Colocación  ( x)  x  xi

Subdominio  ( x)  1,  i  

Galerkin  ( x )  h( x )

Mínimos cuadrados  ( x)  R( x)
H
 ( x)   [ R( x)]2 dx
0
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

1.2.1 - Planteamiento variacional

Formulación en la que una ecuación diferencial se recalcula en una integral

equivalente ponderando la ecuación diferencial de la variable dependiente y
una función de prueba (buscando un máximo o un mínimo).

Evaluación de un funcional en
ˆ (a1 ,..,am )
 j 1,..,m
función de los coeficientes incógnita
ˆ

Condición de estacionariedad del funcional 0 j 1,..,m
ai

Enfoque variacional en el problema elástico. Opciones:

 Por el Principio de Trabajos Virtuales

Obtención de  Minimizando la Energía Potencial Total
desplazamientos
 Otros
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

PROBLEMA
En la Figura se muestra una viga simplemente apoyada de longitud H y
momentos concentrados en los extremos Mo, con módulo de elasticidad E ,
sección transversal de área A y momento de inercia I .

Las condiciones de frontera son y (0) = 0, y (H ) = 0, donde y es la deflexión

de la viga. La ecuación diferencial que gobierna el fenómeno físico es:

d2y
EI 2  M ( x)  0
dx
 Tema 1. Introducción

1. Métodos numéricos para resolución de EDP

Una solución aproximada para la Ecuación podría ser:

x 
y ( x)  A sin  
H 

donde A es un coeficiente a ser determinado. La ecuación satisface las

condiciones de frontera de manera que y (0) = 0;0 y y (H ) = 0;0 y su forma
puede ser similar a la curva real de deflexión.
La solución exacta está dada por:

Mox
y ( x)  xH
2 EI

Aplicar cada uno de los métodos mostrados arriba para encontrar la

solución al problema planteado.
 Tema 1. Introducción

1.- Métodos numéricos para resolución de EDP

1.3.- Método de los elementos finitos (FEM)

Subdivisión en
elementos:

Aproximación de la solución en cada elemento mediante funciones de prueba

Aplicación de (ensamblado)
•principio variacional, o
•método de los residuos ponderados

Resolución del sistema de ecuaciones algebraicas

Postproceso y validación de los resultados numéricos

 Tema 1. Introducción

Método de elementos Finitos

En general los métodos empleados

son de dos tipos:
• Residuos ponderados: se aplican
• Principios variacionales: Asociados a directamente sobre la ecuación
la minimización de algún funcional diferencial y sus condiciones de
contorno, no precisando de ningún
funcional asociado

W( u )    R( u ) dA  0
A
 Tema 1. Introducción

Método de elementos Finitos

 Tema 1. Introducción

1.- Métodos numéricos para resolución de EDP

1.3.- Método de los elementos de contorno (BEM)

Transformación: ecuaciones diferenciales de dominio (D)

integrales de frontera (D-1)

Subdivisión del contorno en elementos

Resolución del sistema de ecuaciones algebraicas

 Tema 1. Introducción

2.- Diferentes enfoques del planteamiento del MEF

Ecs. diferenciales

Variacional (energético)
Planteamiento a nivel
global Enfoques
Ecs. integrales
Residuos ponderados
 Tema 1. Introducción

2.- Diferentes enfoques del planteamiento del MEF

La mejor forma de resolver cualquier problema físico gobernado por

una ecuación diferencial es obtener la solución analítica.

Problemas:
•Frontera
•Multiplicidad de materiales
•Problemas no lineales
•Materiales anisotrópicos

Solución:
Desarrollar un método numérico para hallar la solución a un
problema cuando no se puede hacer por medios analíticos.

La solución es aproximada, en puntos discretos del dominio. Pero es

mejor que no tener nada.
 Tema 1. Introducción

2.- Diferentes enfoques del planteamiento del MEF

2.1 - Planteamiento variacional

Formulación en la que una ecuación diferencial se recalcula en una integral

equivalente ponderando la ecuación diferencial de la variable dependiente y
una función de prueba (buscando un máximo o un mínimo).

Evaluación de un funcional en
ˆ (a1 ,..,am )
 j 1,..,m
función de los coeficientes incógnita
ˆ

Condición de estacionariedad del funcional 0 j 1,..,m
ai

Enfoque variacional en el problema elástico. Opciones:

 Por el Principio de Trabajos Virtuales

Obtención de  Minimizando la Energía Potencial Total
desplazamientos
 Otros
 Tema 1. Introducción

2.- Diferentes enfoques del planteamiento del MEF

2.2.- Residuos ponderados

u campo incógnita exacto
Ecuación diferencial L(u ) f 0 L operador diferencial
A Dominio
f función definida en A
Residuo R(u ) L(u ) f

Forma integral W (u )   R (u )dA Y función de ponderación

A

Residuo y forma R(uˆ ) L(uˆ ) f û campo incógnita aproximado

integral para û W (uˆ )  R(uˆ )dA
A

Forma integral para W j (uˆ )  j R (uˆ )dAW j (a1 ,..,am ) a1,…am coeficientes incógnita
selección de Y’s A
j 1,..,m Yj funciones de ponderación

Anulación de forma W j ( a1 ,..,am )0

j 1,..,m
integral Yj
 Tema 1. Introducción

2.- Diferentes enfoques del planteamiento del MEF

2.2.- Residuos ponderados

Colocación  ( x)  x  xi

Subdominio  ( x)  1,  i  

Galerkin  ( x )  h( x )

Mínimos cuadrados  ( x)  R( x)
H
 ( x)   [ R( x)]2 dx
0
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros

1940 Courant Prager y Hrenikoff
Schoenberg Synge McHenry
Newmark Génesis conceptos
1950 Polya Synge Langefors MEF e Introducción
Hersch McMahon Argyris ordenadores
Weisberger Turner,
Greenstadt Clough,
Martin, y
Topp
1960 Friedrichs Años Formativos
Clough
White Melosh;
Bessenling
Jones;
Fraeijs de
Veubeke
Zienkiewicz
y Cheung
1970 Años maduración
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Courant

1940 Courant Prager y Hrenikoff
Resolución del
Schoenberg Synge McHenry
Newmark Génesis conceptos problema de torsión
1950 Polya MEF e Introducción con polinomios
Synge Langefors
Hersch McMahon Argyris ordenadores lineales en regiones
Weisberger Turner, trianguladas
Greenstadt Clough,
Martin, y
Topp
1960 Friedrichs Años Formativos Schoenberg
Clough
White Melosh;
Bessenling Teoría de SPLINES.
Jones;
Fraeijs de Recomienda el uso
Veubeke de polinomios
Zienkiewicz definidos a tramos
y Cheung para aproximar e
1970 Años maduración interpolar
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Synge

1940 Courant Prager y Hrenikoff
Schoenberg Synge McHenry
Newmark Génesis conceptos
1950 Polya Synge Langefors MEF e Introducción
Hersch McMahon Argyris ordenadores
Weisberger Turner,
Greenstadt Clough,
Martin, y
Topp
1960 Friedrichs Clough
Años Formativos Synge y Prager
White Melosh;
Bessenling Desarrollan el método del hipercírculo:
Jones;
Fraeijs de Interpretación geométrica para los
Veubeke principios de mínimo de la teoría de
Zienkiewicz elasticidad
y Cheung
1970 Años maduración
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Hrenikoff

1940 Courant Prager y Hrenikoff
El comportamiento elástico de una placa
Schoenberg Synge McHenry
Newmark Génesis conceptos es equivalente al de un conjunto de
1950 Polya MEF e Introducción elementos viga conectados entre si
Synge Langefors
Hersch McMahon Argyris ordenadores  Resolución con métodos de
Weisberger Turner, estructuras de barras
Greenstadt Clough,
Martin, y
Topp
1960 Friedrichs Años Formativos McHenry y Newmark
Clough
White Melosh;
Bessenling Refinan esta idea
Jones;
Fraeijs de
Veubeke
Zienkiewicz
y Cheung
1970 Años maduración
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Polya

1940 Courant Prager y Hrenikoff
Schoenberg Synge McHenry
Newmark Génesis conceptos
1950 Polya Synge Langefors MEF e Introducción
Hersch McMahon Argyris ordenadores
Weisberger Turner,
Greenstadt Clough,
Martin, y Polya, Herchs y Weisberger
Topp
1960 Friedrichs Años Formativos Ideas similares a Courant (polinomios
Clough
White Melosh; lineales sobre regiones trianguladas)
Bessenling para evaluar límites de valores propios
Jones;
Fraeijs de Greenstradt
Veubeke
Zienkiewicz Divide un dominio en células,
y Cheung
Años maduración asigna una función diferente a cada una, y
1970
aplica un principio variacional
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Synge

1940 Courant Prager y Hrenikoff
Schoenberg Synge McHenry Usa funciones
Newmark Génesis conceptos lineales sobre una
1950 Polya Synge Langefors MEF e Introducción región triangulada y
Hersch McMahon Argyris ordenadores un procedimiento
Weisberger Turner, variacional de Ritz
Greenstadt Clough,
Martin, y
Topp
1960 Friedrichs Años Formativos McMahon
Clough
White Melosh;
Bessenling Resuelve problema electrostático 3D
Jones;
con tetraedros y funciones lineales
Fraeijs de
Veubeke
Zienkiewicz
y Cheung
1970 Años maduración
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Argyris

1940 Courant Prager y Hrenikoff
Schoenberg Synge McHenry
Newmark Génesis conceptos
1950 Polya Synge Langefors MEF e Introducción
Hersch McMahon Argyris ordenadores
Weisberger Turner,
Greenstadt Clough, Langefors y Argyris
Martin, y
Topp Reformulan el problema de análisis de
Años Formativos estructuras a una forma matricial
1960 Friedrichs Clough
adaptada al calculo mediante ordenador
White Melosh;
Bessenling
Jones; Turner, Clough, Martin y Top
Fraeijs de Modelan aviones (3D) mediante
Veubeke
ensamblado de paneles triangulares.
Zienkiewicz
y Cheung Considerable extension of the material
1970 Años maduración presented in this paper is possible
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Friedrichs

1940 Courant Prager y Hrenikoff
Schoenberg Synge McHenry
Newmark Génesis conceptos
1950 Polya Synge Langefors MEF e Introducción
Hersch McMahon Argyris ordenadores
Weisberger Turner,
Greenstadt Clough,
Martin, y
Topp
Años Formativos Friedrichs y White
1960 Friedrichs Clough
White Melosh; Usan elementos triangulares y
Bessenling principios variacionales
Jones;
Fraeijs de para desarrollar ecuaciones en
Veubeke diferencias finitas
Zienkiewicz
y Cheung
1970 Años maduración
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Clough

1940 Courant Prager y Hrenikoff
Aparece por primera vez el nombre de
Schoenberg Synge McHenry
Newmark Génesis conceptos ELEMENTOS FINITOS
1950 Polya Synge Langefors MEF e Introducción
Hersch McMahon Argyris ordenadores Fraejis de Veubeke
Weisberger Turner,
Greenstadt Clough,
Martin, y
Topp
1960 Friedrichs Años Formativos
Clough
White Melosh;
Bessenling
Jones; Melosh y Besseling
Fraeijs de
Veubeke Jones y Fraejis de Veubeke
Zienkiewicz
y Cheung Muestran el MEF como un método
1970 Años maduración variacional de Ritz usando
funciones definidas a tramos
 Tema 1. Introducción

3.- CRONOLOGÍA

Matemáticos Físicos Ingenieros Zienkiewicz

1940 Courant Prager y Hrenikoff
Schoenberg Synge McHenry
Newmark Génesis conceptos
1950 Polya Synge Langefors MEF e Introducción
Hersch McMahon Argyris ordenadores
Weisberger Turner,
Greenstadt Clough,
Martin, y
Topp
1960 Friedrichs Años Formativos
Clough
White Melosh;
Bessenling
Jones; Zienkiewicz y Cheung
Fraeijs de
Veubeke Muestran que el MEF es aplicable a
Zienkiewicz todos los problemas que se puedan
y Cheung
definir en forma variacional
1970 Años maduración