Una introducción a los elementos finitos
Ricardo Prato Torres
Universidad del Norte
Departamento de matemáticas y estadı́stica
Grupo de investigación en Matemáticas
Aplicadas
28.02.15
Página 1 FEM 28.02.15 R. Prato
Método de los elementos finitos (FEM)
Una aproximación numérica para la solución de ecuaciones
diferenciales parciales, problemas de valor iniciales (PVI) o problemas
de valores en la frontera (BVP) -en el caso de dominios acotados-
Página 2 FEM 28.02.15 R. Prato
Método de los elementos finitos (FEM)
Una aproximación numérica para la solución de ecuaciones
diferenciales parciales, problemas de valor iniciales (PVI) o problemas
de valores en la frontera (BVP) -en el caso de dominios acotados-
◮ Elementos finitos.
◮ Método de diferencias divididas.
◮ Método de volumenes finitos.
◮ Método de elementos de frontera.
◮ Métodos meshless.
Página 2 FEM 28.02.15 R. Prato
Método de los elementos finitos (FEM)
Una aproximación numérica para la solución de ecuaciones
diferenciales parciales, problemas de valor iniciales (PVI) o problemas
de valores en la frontera (BVP) -en el caso de dominios acotados-
◮ Elementos finitos.
◮ Método de diferencias divididas.
◮ Método de volumenes finitos.
◮ Método de elementos de frontera.
◮ Métodos meshless.
Ventajas del método
Página 2 FEM 28.02.15 R. Prato
Método de los elementos finitos (FEM)
Una aproximación numérica para la solución de ecuaciones
diferenciales parciales, problemas de valor iniciales (PVI) o problemas
de valores en la frontera (BVP) -en el caso de dominios acotados-
◮ Elementos finitos.
◮ Método de diferencias divididas.
◮ Método de volumenes finitos.
◮ Método de elementos de frontera.
◮ Métodos meshless.
Ventajas del método
◮ Soluciones aproximadas para problemas complejos.
Página 2 FEM 28.02.15 R. Prato
Método de los elementos finitos (FEM)
Una aproximación numérica para la solución de ecuaciones
diferenciales parciales, problemas de valor iniciales (PVI) o problemas
de valores en la frontera (BVP) -en el caso de dominios acotados-
◮ Elementos finitos.
◮ Método de diferencias divididas.
◮ Método de volumenes finitos.
◮ Método de elementos de frontera.
◮ Métodos meshless.
Ventajas del método
◮ Soluciones aproximadas para problemas complejos.
◮ Sistemas lineales de fácil solución.
Página 2 FEM 28.02.15 R. Prato
FEM 1d
Dominio I := (0, 1)
−(τ (x)u′ (x))′ + q(x)u(x) − f (x) = 0
u(0) = u(1) = 0
donde τ, q y f son funciones dadas, definidas en I. Asumimos además que
τ y q son positivas y acotadas, es decir
0 < τmin ≤ τ (x) ≤ τmax , 0 < qmin ≤ q(x) ≤ qmax .
Página 3 FEM 28.02.15 R. Prato
− (τ (x)u′ (x))′ + q(x)u(x) = f (x)
Página 4 FEM 28.02.15 R. Prato
− (τ (x)u′ (x))′ + q(x)u(x) · v(x) = f (x)·v(x)
Página 4 FEM 28.02.15 R. Prato
Z 1 Z 1
′ ′
− (τ (x)u (x)) + q(x)u(x) · v(x) dx = f (x)·v(x) dx
0 0
Página 4 FEM 28.02.15 R. Prato
Z 1h i Z 1
′ ′
− (τ (x)u (x)) v(x) + q(x)u(x)v(x) dx = f (x)v(x) dx
0 0
Página 4 FEM 28.02.15 R. Prato
Z 1h i Z 1
′ ′
− (τ (x)u (x)) v(x) + q(x)u(x)v(x) dx = f (x)v(x) dx
0 0
Z 1h i Z 1
τ (x)u′ (x)v ′ (x) + q(x)u(x)v(x) dx − [τ (x)u′ (x)v(x)]|10 = f (x)v(x) dx
0 0
Página 4 FEM 28.02.15 R. Prato
Z 1h i Z 1
− (τ (x)u′ (x))′ v(x) + q(x)u(x)v(x) dx = f (x)v(x) dx
0 0
Z 1h i Z 1
τ (x)u′ (x)v ′ (x) + q(x)u(x)v(x) dx−[τ (x)u′ (x)v(x)]|10 = f (x)v(x) dx
0 | {z } 0
=0
Página 4 FEM 28.02.15 R. Prato
Z 1h i Z 1
− (τ (x)u′ (x))′ v(x) + q(x)u(x)v(x) dx = f (x)v(x) dx
0 0
Z 1h i Z 1
τ (x)u′ (x)v ′ (x) + q(x)u(x)v(x) dx−[τ (x)u′ (x)v(x)]|10 = f (x)v(x) dx
0 | {z } 0
=0
Z 1h i Z 1
τ (x)u′ (x)v ′ (x) + q(x)u(x)v(x) dx = f (x)v(x) dx
0 0
Página 4 FEM 28.02.15 R. Prato
Z 1h i Z 1
− (τ (x)u′ (x))′ v(x) + q(x)u(x)v(x) dx = f (x)v(x) dx
0 0
Z 1h i Z 1
τ (x)u′ (x)v ′ (x) + q(x)u(x)v(x) dx−[τ (x)u′ (x)v(x)]|10 = f (x)v(x) dx
0 | {z } 0
=0
Z 1h i Z 1
τ (x)u′ (x)v ′ (x) + q(x)u(x)v(x) dx = f (x)v(x) dx
0 0
B(u, v) = (f, v), for all v ∈ H01 (I)
Página 4 FEM 28.02.15 R. Prato
Propiedades de la forma bilineal
1 Simetrı́a:
B(u, v) = B(v, u),
Página 5 FEM 28.02.15 R. Prato
Propiedades de la forma bilineal
1 Simetrı́a:
B(u, v) = B(v, u),
2 bilinealidad:
B(αu + βv, w) = αB(u, w) + βB(v, w)
B(u, αv + βw) = αB(u, v) + βB(u, w)
Página 5 FEM 28.02.15 R. Prato
Propiedades de la forma bilineal
1 Simetrı́a:
B(u, v) = B(v, u),
2 bilinealidad:
B(αu + βv, w) = αB(u, w) + βB(v, w)
B(u, αv + βw) = αB(u, v) + βB(u, w)
3 Continuidad:
|B(u, v)| ≤ C1 kuk1 kvk1 ,
Página 5 FEM 28.02.15 R. Prato
Propiedades de la forma bilineal
1 Simetrı́a:
B(u, v) = B(v, u),
2 bilinealidad:
B(αu + βv, w) = αB(u, w) + βB(v, w)
B(u, αv + βw) = αB(u, v) + βB(u, w)
3 Continuidad:
|B(u, v)| ≤ C1 kuk1 kvk1 ,
4 Coercividad:
B(u, u) ≥ C2 kuk21 ,
Página 5 FEM 28.02.15 R. Prato
Theorem (Lax-Milgram)
Sea H un espacio de Hilbert, B(·, ·) : H × H → IR una forma bilineal en
un continua y H-coerciva, y F : H → IR un funcional lineal continuo en
H. Entonces, existe un único u ∈ H que satisface
B(u, v) = F (v), for all v ∈ H.
Más aun, kuk ≤ C2−1 kF k.
Página 6 FEM 28.02.15 R. Prato
Sea S un subespacio finito dimensional de H01 (I). La aproximación de
Galerkin U ∈ S a la solución u se define por
B(U, v) = (f, v), for all v ∈ S ⊂ H01 (I)
esto da lugar a la condición de ortogonalidad
B(u − U, v) = 0, for all v ∈ S ⊂ H01 (I)
Página 7 FEM 28.02.15 R. Prato
Sea S un subespacio finito dimensional de H01 (I). La aproximación de
Galerkin U ∈ S a la solución u se define por
B(U, v) = (f, v), for all v ∈ S ⊂ H01 (I)
esto da lugar a la condición de ortogonalidad
B(u − U, v) = 0, for all v ∈ S ⊂ H01 (I)
Sea {φj }N j=1 una base para el espacio S; Por lo tanto, existen constantes
c1 , . . . , cN tales que
XN
U (x) = cj φj (x).
j=1
Página 7 FEM 28.02.15 R. Prato
Sea S un subespacio finito dimensional de H01 (I). La aproximación de
Galerkin U ∈ S a la solución u se define por
B(U, v) = (f, v), for all v ∈ S ⊂ H01 (I)
esto da lugar a la condición de ortogonalidad
B(u − U, v) = 0, for all v ∈ S ⊂ H01 (I)
Sea {φj }N j=1 una base para el espacio S; Por lo tanto, existen constantes
c1 , . . . , cN tales que
XN
U (x) = cj φj (x).
j=1
Reemplazando esto en la forma debil y fijando v = φk :
N
X
B(φj , φk )cj = (f, φk ), k = 1, . . . , N
j=1
Página 7 FEM 28.02.15 R. Prato