Introducción a la teoría matemática

del Método del elemento finito

Dandy Rueda Castillo

Noviembre 25, 2016

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Objetivo

Introducir la teoría matemática del elemento finito en 1D: Método de

Galerkin.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Espacios de funciones polinomiales por tramos
El espacio de polinomios lineales

Sean x0 y x1 puntos distintos tales que x0 < x1 y a los cuales nos

referiremos como nodos. Considerando el intervalo I = [x0 , x1 ], defini-
mos el espacio vectorial

P1 (I ) = {v : v (x ) = c0 + c1 x , x ∈ I , c0 , c1 ∈ R}

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...
El espacio de polinomios lineales

v queda determinado conociendo c0 y c1 , puesto que el conjunto

{1, x } es una base para P1 (I ) . Además, diremos que v tiene dos
grados de libertad.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...
El espacio de polinomios lineales

v queda determinado conociendo c0 y c1 , puesto que el conjunto

{1, x } es una base para P1 (I ) . Además, diremos que v tiene dos
grados de libertad.
Si v ∈ P1 (I ), dados los valores α0 = v (x0 ) y α1 = v (x1 ), que
llamaremos valores nodales, hay una única función v en P1 (I )
que cumple con las condiciones dadas.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...
El espacio de polinomios lineales

En efecto, el sistema
    
1 x0 c0 α0
= ,
1 x1 c1 α1

tiene solución única.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...
Base nodal

Introducimos una nueva base {λ0 , λ1 } para P1 (I ). Esta nueva base

llamada base nodal se define por

1, i =j
λj (xi ) = i , j = 0, 1
0, i 6= j

y puesto λ0 , λ1 ∈ P1 (I ) se tiene explícitamente que

x1 − x x − x0
λ0 (x ) = , λ1 (x ) = .
x1 − x0 x1 − x0

de esta manera cada función λj , j = 0, 1 es lineal y toma valores 0 y 1.

Luego, v se puede escribir como sigue

v (x ) = α0 λ0 (x ) + α1 λ1 (x ).

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
El espacio de funciones polinomiales lineales por tramos

Sea I = [a, b] y los puntos {xi }ni=0 que definen una partición de I tales
que
a = x0 < x1 < · · · < xn = b

(nos referiremos a esto como una malla o mesh.). Teniendo en cuenta

los subintervalos Ii = [xi −1 , xi ], i = 1, 2, . . . , n de longitud hi = xi − xi −1 ,
definimos
Vh = {v ∈ C 0 (I ) : v |Ii ∈ P1 (Ii )}

donde C 0 (I ) es el espacio de funciones continuas definidas sobre I .

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Toda función v en Vh está únicamente determinado por sus

valores nodales e inversamente para cualquier conjunto de
valores nodales existe una función en Vh con estos valores
nodales.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...
Funciones sombrero

Introducimos la base {φj }nj=0 para Vh tal que


1, i =j
ϕj (xi ) =
0, i 6= j

las funciones ϕi las llamamos funciones sombrero o “hat”.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Por construcción para una función v en Vh con valores nodales {αi }ni=0
se tiene
n
v (x ) = ∑ αi ϕi (x )
i =0

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Además, explícitamente tenemos

 x1 − x , x ∈ I

0
ϕ0 (x ) = x1 − x0
 0, otro caso
 x − xn−1 , x ∈ I

n
ϕn (x ) = xn − xn−1
 0, otro caso
 x −x
i −1
 , x ∈ Ii ,
 xi − xi −1


xi +1 − x
ϕi (x ) = , x ∈ Ii +1 , i = 1, . . . , n − 1.


 xi + 1 − xi
0, otro caso


Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Interpolación lineal

Teniendo en cuenta I = [x0 , x1 ] y dada una función continua f en

I, definimos:

π f (x ) = f (x0 )ϕ0 (x ) + f (x1 )ϕ1 (x )

es claro que π f ∈ P1 (I ).

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Interpolación lineal

Teniendo en cuenta I = [x0 , x1 ] y dada una función continua f en

I, definimos:

π f (x ) = f (x0 )ϕ0 (x ) + f (x1 )ϕ1 (x )

es claro que π f ∈ P1 (I ).
Error de interpolación: f − π f

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Interpolación lineal

Teniendo en cuenta I = [x0 , x1 ] y dada una función continua f en

I, definimos:

π f (x ) = f (x0 )ϕ0 (x ) + f (x1 )ϕ1 (x )

es claro que π f ∈ P1 (I ).
Error de interpolación: f − π f
Usamos la L2 (I )-norma definida para cualquier función cuadrado
integrable
Z 1/2
2
||v ||L2 (I ) = v dx
I

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Teorema 1.
π f satisface

||f − π f ||L2 (I ) ≤ ch2 ||f 00 ||L2 (I )

||(f − π f )0 ||L2 (I ) ≤ ch||f 00 ||L2 (I )

c es una constante y h = x1 − x0 .

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Interpolación lineal continua por tramos

Dada una función continua f definida en I = [a, b] y una partición {xi }ni=0
de I tal que
a = x0 < x1 < · · · < xn = b

definimos
n
π f (x ) = ∑ f (xi )ϕi (x )
i =0

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Teorema 2.
π f satisface
n 2
||f − π f ||2L2 (I ) ≤ C ∑ hi4 f 00 L2 (Ii )
i =1
n
(f − π f )0 22
2
∑ hi2 f 00 L2 (Ii )

L (I )
≤ C
i =1

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
L2 -proyección

Dada f ∈ L2 (I ), la L2 -proyección Ph f ∈ Vh de f es definida por

Z
(f − Ph f )vdx = 0, ∀v ∈ Vh .
I

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Cálculo de Ph f

Para calcular Ph f observar que la definición es equivalente a

Z
(f − Ph f )ϕi dx = 0, i = 0, . . . , n
I

y puesto que Ph f ∈ Vh se tiene

n
Ph f = ∑ ξj ϕj , j = 0, 1, 2, . . . , n
j =0

donde ξj deben determinarse. Notemos que para i = 0, 1, . . . , n se

tiene Z Z Z
(f − Ph f )ϕi dx = 0 ⇔ f ϕi dx = Ph f ϕi dx
I I

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...
Matriz de masa y vector de carga

Trabajamos el lado derecho

!
Z Z n
Ph f ϕi dx = ∑ ξj ϕj ϕi dx
I I j =0
n Z
= ∑ ξj ϕi ϕj dx
j =0 I

Introducimos
Z
Mij = ϕi ϕj dx , i , j = 0, 1, . . . , n
ZI
bi = f ϕi dx , i = 0, 1, . . . , n
I

y tenemos el sistema
Mξ = b
donde nos referiremos a M como matriz de MASA y b vector de CARGA.
Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...
Ensamblaje de la matriz de masa

h1 hn
M0,0 = , Mn,n = .
3 3
hi hi + 1
Mi ,i = + , i = 1, 2, . . . , n − 1,
3 3
hi +1
Mi ,i +1 = , i = 1, 2, . . . , n − 1,
6
hi +1
Mi +1,i = , i = 1, 2, . . . , n − 1,
6
h1 h1
···
 
3 6
0 0 0
h1 h1

6 3
+ h32 h2
6
··· 0 0 
h2 h2 h3 h3
 

 0 6 3
+ 3 6
0 0 

M =
 0 0 

 .. .. .. .. .. .. 
 . . . . . . 
 hn−1 hn−1

 0 0 0 6 3
+ h3n 6 hn 
hn hn
0 0 0 ··· 6 3 del elemento finito
Dandy Rueda Castillo Introducción a la teoría matemática del Método
...

h1 h1
     
3 6
0 0 0 0
h1 h1 h2 h2

 6 3
 
  3 6



 0 0 

h2 h2
 0 0  
6 3
  0 0 
M = + +···+
     
.. .. .. .. .. .. 

 . .  
  . . 


 . . 

hn hn
 0 0   0 0  
3 6

hn hn
0 0 0 0 6 3
| {z } | {z } | {z }
M I1 M I2 M In

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
....

h1 h1
     
3 6
0 0 0 0
h1 h1 h2 h2

 6 3
 
  3 6



 0 0 

h2 h2
 0 0  
6 3
  0 0 
M = + +···+
     
.. .. .. .. .. .. 

 . .  
  . . 


 . . 

hn hn
 0 0   0 0  
3 6

hn hn
0 0 0 0 6 3
| {z } | {z } | {z }
M I1 M I2 M In

Cada matriz es obtenida al restringir la integración al subintervalo Ii .

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Notemos que sobre cada elemento Ii se tienen los bloques

 
hi 2 1
M Ii =
6 1 2

donde hi es la longitud del intervalo Ii . También M Ii es llamado ele-

mento de masa local y el ensamblaje es la suma de estos bloques y
forman M que es la matriz de masa del global.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Algoritmo

1 Crear malla
2 Calcular M y b
3 Resolver M ξ = b
4 Establecer Ph f = ∑nj=0 ξj ϕj .

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Estimador de error apriori

Teorema 3.
La L2 -proyección Ph f satisface el mejor resultado de aproximación

||f − Ph f ||L2 (I ) ≤ ||f − v ||L2 (I ) , ∀v ∈ Vh

Ph f es el minimizador de

min ||f − v ||L2 (I )

v ∈Vh

y se dice que Ph f aproxima a f en el sentido de los mínimos cuadrados.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Estimador de error apriori

Teorema 3.
La L2 -proyección Ph f satisface el mejor resultado de aproximación

||f − Ph f ||L2 (I ) ≤ ||f − v ||L2 (I ) , ∀v ∈ Vh

Ph f es el minimizador de

min ||f − v ||L2 (I )

v ∈Vh

y se dice que Ph f aproxima a f en el sentido de los mínimos cuadrados.

Teorema 4.
La L2 -proyección Ph f satisface el estimado
n
||f − Ph f ||2L2 (I ) ≤ c ∑ hi4 ||fL0022(Ii )
i =0

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Método del elemento finito

Consideremos el problema de valor de frontera

 −u 00 = f , x ∈ I = [0, L]


u (0) = u (L) = 0


Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Formulación variacional

Multiplicamos la ecuación por una función de prueba v , la cual se anula

en los extremos del intervalo, e integramos por partes
Z L Z L
fvdx = − u 00 vdx
0 0
Z L
= u 0 v 0 dx − u 0 (L) v (L) + u 0 (0) v (0)
0
Z L
= u 0 v 0 dx
0

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

La más grande colección de funciones v tal que v (0) = v (L) = 0 está

dada por el espacio
n o
V0 = v : ||v ||2L2 (I ) < ∞, ||v 0 ||2L2 (I ) < ∞, v (0) = v (L) = 0

Notemos que u pertenece a esta clase de funciones.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Problema variacional

Encontrar u ∈ V0 tal que

Z L Z L
u 0 v 0 dx = fvdx , ∀v ∈ V0
0 0

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Aproximación por elementos finitos

Introducimos el espacio

Vh,0 = {v ∈ Vh : v (0) = v (L) = 0}

n
donde una base es el conjunto de funciones sombrero {ϕj }j =1 . Como
Vh,0 ⊂ V0 , en la formulación variacional remplazamos V0 por Vh,0 y
el problema se convierte en el método de elemento finito: Encontrar
uh ∈ Vh,0 tal que
Z L Z L
uh0 v 0 dx = fvdx , ∀v ∈ Vh,0
0 0

y lo llamaremos el método de Galerkin.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Para hallar uh notemos que

Z Z
0 0
uh ϕi dx = f ϕi dx , i = 1, 2, . . . , n − 1
I I

y puesto que uh ∈ Vh,0 se tiene

n−1
uh = ∑ ξj ϕj
j =1

donde los coeficientes se deben determinarse.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

En efecto, para i = 1, 2, . . . , n − 1 se tiene

!
Z Z Z n−1 n−1 Z
0 0 0
f ϕi dx = uh ϕi dx = ∑ ξj ϕj ϕi0 dx = ∑ ξj ϕj0 ϕi0 dx
I I I j =1 j =1 I

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Introduciendo la notación
Z
Aij = ϕj0 ϕi0 dx
I
Z
bi = f ϕi dx
I

tenemos el sistema

Aξ = b

Nos referiremos a A como la matriz de rigidez y a b como el vector de

carga.

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Teorema 5.
La aproximación uh satisface la condición de ortogonalidad
Z
(u − uh )0 v 0 dx = 0, ∀v ∈ Vh,0
I

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Teorema 5.
La aproximación uh satisface la condición de ortogonalidad
Z
(u − uh )0 v 0 dx = 0, ∀v ∈ Vh,0
I

Teorema 6.
La solución uh satisface el mejor resultado de aproximación

(u − uh )0 2 ≤ (u − v )0 2 , ∀v ∈ Vh,0

L (I ) L (I )

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
...

Teorema 5.
La aproximación uh satisface la condición de ortogonalidad
Z
(u − uh )0 v 0 dx = 0, ∀v ∈ Vh,0
I

Teorema 6.
La solución uh satisface el mejor resultado de aproximación

(u − uh )0 2 ≤ (u − v )0 2 , ∀v ∈ Vh,0

L (I ) L (I )

Teorema 7.
La solución uh satisface el estimado
n
(u − uh )0 2 ≤ C ∑ hi2 u 00 2

L (I ) L (I )
i =1

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito
Referencias I

Malte Braack. Elemente Finite. Vorlesungsskript, 25.5.2012

Mats. G. Larson - Fredrik Bengzon. The Finite Element Method
Theory, Implementation and Applications. 2013
Susanne C. Brenner, Ridgway Scott. The Mathematical Theory of
Finite Element Methods. 2008

Dandy Rueda Castillo Introducción a la teoría matemática del Método del elemento finito