Técnicas de Aproximación

Ing. Nelson M. Lafontaine, Ph. D

FUNDEXTENSUS
Fundación para la Educación y el Desarrollo

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 1 / 40

 Índice

1 Técnicas de Aproximación
Forma Fuerte del Problema
Collocation Method
Least Square Method
Método de Galerkin
Mejoramiento de la solución
Formulación Débil de una ecuación diferencial

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 2 / 40

 Técnicas de Aproximación

Sección 1

Técnicas de Aproximación

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 3 / 40

 Técnicas de Aproximación

Técnicas de Aproximación
Métodos de los Residuos Ponderados

Objetivos
Resolución aproximada de ecuaciones diferenciales mediante Residuos
Ponderados empleando los métodos de Colocación, Mı́nimos
Cuadrados y Galerkin.

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 4 / 40

 Técnicas de Aproximación Forma Fuerte del Problema

Técnicas de Aproximación
Métodos de los Residuos Ponderados

Forma Fuerte del Problema

Considere la siguiente ecuación diferencial:

d 2u
− u = −x , en dominio (0, 1) (1.1)
dx 2
u(0) = 0 y u(1) = 0 Condición de frontera (1.2)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 5 / 40

 Técnicas de Aproximación Forma Fuerte del Problema

Técnicas de Aproximación
Métodos de los Residuos Ponderados

Asuma que funciones de prueba o de test contienen constantes

desconocidas αi , que aproxime la solución de la ecuación 1.1, como por
ejemplo:

ū = αx (1 − x ) (1.3)

donde ū es solución aproximada de la ecuación diferencial 1.1.

Requisitos de la función de prueba
La función de test se selecciona de tal manera que cumpla con las
condiciones de frontera de la ecuación 1.2.

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 6 / 40

 Técnicas de Aproximación Forma Fuerte del Problema

Técnicas de Aproximación
Métodos de los Residuos Ponderados

Se define el Residuo R como:

d 2 ū
R= − ū + x (1.4)
dx 2
= −2α − αx (1 − x ) + x (1.5)

Problema a resolver
Encontrar α que aproxime mejor la solución exacta y que el Residuo sea
nulo.

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 7 / 40

 Técnicas de Aproximación Forma Fuerte del Problema

Técnicas de Aproximación
Métodos de los Residuos Ponderados

Para determinar α se introduce una Función de Peso w

Z 1
I= wRdx = 0 (1.6)
0
!
d 2 ū
Z 1
= w − ū + x =0 (1.7)
0 dx 2
Z 1
= w (−2α − αx (1 − x ) + x )dx = 0 (1.8)
0

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 8 / 40

 Técnicas de Aproximación Collocation Method

Técnicas de Aproximación
Método de Colocación

El método de los Residuos Ponderados se basa en como se selecciona la

Función de Peso w . Se explican tres métodos a continuación:
Método de Colocación (Collocation Method CM)
w = δ (x − xo ) (1.9)

Propiedad de la función delta de Dirac

Z b  
f (xo ) si a < xo < b;
δ (x − xo ) f (x )dx = (1.10)
a
0 de otro modo

Sustituyendo en xo = 0.50
Z 1
I= δ (x − 0.50) (−2α − αx (1 − x ) + x )dx = 0 (1.11)
0
α = 2/9 = 0.2̄ (1.12)

u(x ) ' ū(x ) = 0.2̄x (1 − x ) Solución aproximada CM (1.13)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 9 / 40

 Técnicas de Aproximación Least Square Method

Técnicas de Aproximación
Métodos de los Mı́nimos Cuadrados

Método de Mı́nimos Cuadrados (Least Square Method LSM)

dR
w= = −2 − x (1 − x ) (1.14)
dα
Sustituyendo
Z 1
I= (−2 − x (1 − x )) (−2α − αx (1 − x ) + x )dx = 0 (1.15)
0
α = 0.2305 (1.16)
u(x ) ' ū(x ) = 0.2305x (1 − x ) Solución aproximada LSM (1.17)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 10 / 40

 Técnicas de Aproximación Método de Galerkin

Técnicas de Aproximación
Métodos de Galerkin

Método de Galerkin (Galerkin Method GM)

d ū
w= = x (1 − x ) (1.18)
dα
Sustituyendo
Z 1
I= (x (1 − x )) (−2α − αx (1 − x ) + x )dx = 0 (1.19)
0
α = 0.2272 (1.20)
u(x ) ' ū(x ) = 0.2272x (1 − x ) Solución aproximada GM (1.21)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 11 / 40

 Técnicas de Aproximación Método de Galerkin

Técnicas de Aproximación
Métodos de los Residuos Ponderados

Solución Exacta CM LSM GM

0.0566 0.0556 0.0576 0.0568
Cuadro 1: Comparación de los resultados en x = 0.50

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 12 / 40

 Técnicas de Aproximación Mejoramiento de la solución

Técnicas de Aproximación
Métodos de los Residuos Ponderados

Mejoramiento de la Solución
Función de prueba con (2) variables αi desconocidas a encontrar

ū = α1 x (1 − x ) + α2 x 2 (1 − x ) (1.22)
= α1 Φ1 (x ) + α2 Φ2 (x ) (1.23)

Residuo R
d 2 ū
R= − ū + x (1.24)
dx 2
= α1 (−2 − x + x 2 ) + α2 (2 − 6x − x 2 + x 3 ) + x (1.25)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 13 / 40

 Técnicas de Aproximación Mejoramiento de la solución

Técnicas de Aproximación
Métodos de los residuos ponderados

Collocation method

w1 = δ(x − x1 ) ; w2 = δ(x − x2 ) (1.26)

Least Square Method

dR dR
w1 = = −2 − x + x 2 ; w2 = = 2 − 6x − x 2 + x 3 (1.27)
α1 α2
Galerkin Method
d ū d ū
w1 = = Φ1 (x ) ; w2 = = Φ2 (x ) (1.28)
α1 α2

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 14 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Formulación Débil
Z 1
I= wRdx = 0 (1.29)
0
!
d 2 ū
Z 1
= w − ū + x =0 (1.30)
0 dx 2
Z 1
dw d ū d ū 1

= − − w ū + wx dx + w | =0 (1.31)
0 dx dx dx 0
Integración por partes

d 2 ū
Z 1 Z 1
d ū dw d ū
w = w |10 − (1.32)
0 dx 2 dx 0 dx dx

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 15 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulacón Débil

Observaciones
Las funciones de test tienen que tener una primera derivada y no
anularse
La solución aproximada dependerá enormemente de la función de
prueba seleccionada

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 16 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulacón Débil

Considérese las siguientes funciones lineales definidas a trozos

( (x −x ) )
i−1
hi para xi−1 < x < xi
Φ(x ) = (xi+1 −x ) (1.33)
hi+1 para xi < x < xi+1

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 17 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulacón Débil

Z 1
dw d ū d ū 1

− − w ū + wx dx + w | =0 (1.34)
0 dx dx dx 0
(1.35)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 18 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulacón Débil

Considérese las siguientes funciones definidas a trozos

 
 3x
 para 0 < x < 1/3 

Φ1 (x ) = 2 − 3x para 1/3 < x < 2/3 (1.36)

 0 
para 2/3 < x < 1 

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 19 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Considérese las siguientes funciones definidas a trozos

 
 0
para 0 < x < 1/3 

Φ2 (x ) = 3x − 1 para 1/3 < x < 2/3 (1.37)
 3 − 3x para 2/3 < x < 1
 


Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 20 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Función de prueba
 
 α1 (3x )
 para 0 < x < 1/3 

ū = α1 (2 − 3x ) + α2 (3x − 1) para 1/3 < x < 2/3 (1.38)
 α (3 − 3x )

para 2/3 < x < 1


2

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 21 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Funciónes de test. Método de Galerkin

 
d ū  3x
 para 0 < x < 1/3 

w1 = = Φ1 (x ) = 2 − 3x para 1/3 < x < 2/3 (1.39)
dα1 
 0 
para 2/3 < x < 1 

 
d ū  0
para 0 < x < 1/3 
 
w2 = = Φ2 (x ) = 3x − 1 para 1/3 < x < 2/3 (1.40)
dα2
3 − 3x para 2/3 < x < 1 

 

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 22 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Residuos ponderados
Z 1
dw1 d ū d ū 1

I1 = − − w1 ū + w1 x dx + w1
| =0 (1.41)
0 dx dx dx 0
Z 1
dw2 d ū d ū

I2 = − − w2 ū + w2 x dx + w2 |10 = 0 (1.42)
0 dx dx dx
Obsérvese que:

w1 (0) = w1 (1) = w2 (0) = w2 (1) = 0 (1.43)

d ū 1
w1 | =0 (1.44)
dx 0
d ū
w2 |10 = 0 (1.45)
dx

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 23 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Residuos ponderados
Z 1/3
I1 = (−3(3α1 ) − 3x (3α1 x ) + x (3x )) dx
0
Z 2/3
+ (3(−3α1 + 3α2 ) − (2 − 3x )(2α1 − 3α1 x + 3α2 x − α2 ) + x (2 − 3x )) dx
1/3
Z 1
+ 0dx
2/3
I1 = −6.222α1 + 2.9444α2 + 0.11111 = 0 (1.46)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 24 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Residuos ponderados
Z 1/3
I2 = 0dx
0
Z 2/3
+ (3(−3α1 + 3α2 ) − (3x − 1)(2α1 − 3α1 x + 3α2 x − α2 ) + x (3x − 1)) dx
1/3
Z 1
+ (3(−3α2 ) − (3 − 3x )(3α2 − 3α2 x ) + x (3 − 3x )) dx
2/3
I1 = 2.9444α1 − 6.222α2 + 0.22222 = 0 (1.47)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 25 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Formulación Débil

Resolviendo

I1 = −6.222α1 + 2.9444α2 + 0.11111 = 0 (1.48)

I2 = 2.9444α1 − 6.222α2 + 0.22222 = 0 (1.49)

Solución

α1 = 0.0448 y α2 = 0.0569 (1.50)

ū = 0.0448Φ1 (x ) + 0.0569Φ2 (x ) (1.51)

Observación
Si se hubiera utilizado la forma fuerte del problema, las funciones de peso
no fueran aplicables ya que su segunda derivada se anula.

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 26 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Formulación de elementos finitos basado en Galerkin

Considere
Solución Exacta

Solución EF

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 27 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Formulación de elementos finitos basado en Galerkin

Aproximación de la solución empleando variables nodales discretas

u ' uh = N1 (x )u1 + N2 (x )u2 + N3 (x )u3 + N4 u4 (1.52)

xi+1 − x
N1e (x ) = (1.53)
hi
x − xi
N2e (x ) = (1.54)
hi
hi = xi+1 − xi (1.55)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 28 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Formulación de elementos finitos basado en Galerkin

Correspondencia entre función de test global y local

N1 = N11 para 0 < x < 1/3 (1.56)

(
N21 para 0 < x < 1/3
N2 = (1.57)
N12 para 1/3 < x < 2/3
(
N22 para 1/3 < x < 2/3
N3 = (1.58)
N13 para 2/3 < x < 1
N4 = N23 para 2/3 < x < 1 (1.59)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 29 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Formulación de elementos finitos basado en Galerkin

Correspondencia entre función de test global y local

N11 = 1 − 3x para 0 < x < 1/3 (1.60)

N21 = 3x para 0 < x < 1/3 (1.61)
N12 = 2 − 3x para 1/3 < x < 2/3 (1.62)
N22 = 3x − 1 para 1/3 < x < 2/3 (1.63)
N13 = 3 − 3x para 2/3 < x < 1 (1.64)
N23 = 3x − 2 para 2/3 < x < 1 (1.65)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 30 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Formulación de elementos finitos basado en Galerkin

Empleando EF de interpolación lineal, el residuo ponderado puede
escribirse

n Z x i+1 
dw du du xi+1
X 
I= − − wu + wx dx + w | (1.66)
e=1 xi
dx dx dx xi

Particularizando para elemento e = 1

h i
w1 = [N1 , N2 ] = N11 , N21

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 31 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

dw dN1 dN2 dN3 dN4

 
= , , , (1.67)
dx dx dx dx dx
du duh dN1 dN2 dN3 dN4
' = u1 + u2 + u3 + u4 (1.68)
dx dx dx dx dx dx

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 32 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Particularizando para el elemento e = 1

Z 1/3
h iT h i " u1
#
dN1 dN2 dN1 dN2 T u2
− , , 0, 0 , , 0, 0 + [N1 , N2 , 0, 0] [N1 , N2 , 0, 0] dx
dx dx dx dx u3
0 u4
Z 1/3
T


+ [N1 , N2 , 0, 0] x dx
0

T du T du
+ [N1 , N2 , 0, 0] |x =1/3 |x =1/3 − [N1 , N2 , 0, 0] |x =0 |x =0 = 0 (1.69)
dx dx

 dN1 dN1 dN1 dN2

( ) ( )
Z 1/3  dx dx
+ N1 N1 dx dx
+ N1 N2 0 0  u1 Z 1/3 N1
dN2 dN1 dN2 dN2 u2 N2
− dx dx
+ N2 N1 dx dx
+ N2 N2 0 0 dx + x dx (1.70)
0 0 0 0 u3 0
0

0 0 0 0
 u4 0 0
( ) ( )
0 1
1 0 0 0
+ u (1/3) − u (0) (1.71)
0 0
0 0

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 33 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Particularizando para el elemento e = 1

−u 0 (0)
( )( ) ( ) ( )
−3.111 2.9444 0 0 u1 0.0185
2.944 −3.111 0 0 u2 0.0370 u 0 (1/3)
+ + (1.72)
0 0 0 0 u3 0 0
0 0 0 0 u4 0 0

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 34 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Particularizando para el elemento e = 2

Z 2/3
h iT h i " u1
#
dN2 dN3 dN2 dN3 T u2
− 0, , ,0 0, , , 0 + [0, N2 , N3 , 0] [0, N2 , N3 , 0] dx
dx dx dx dx u3
1/3 u4
Z 2/3
T


+ [0, N2 , N3 , 0] x dx
1/3

T du T du
+ [0, N2 , N3 , 0, ] |x =2/3 |x =2/3 − [0, N2 , N3 , 0] |x =1/3 |x =1/3 = 0 (1.73)
dx dx

 (
0 0 0 0 u1
) (
0
)
Z 2/3  dN2 dN2 dN2 dN3
 Z 2/3
0 dx dx
+ N2 N2 dx dx
+ N2 N3 0 u2 N2
− dN3 dN2 dN3 dN3 dx + x dx (1.74)
0 + N3 N2 + N3 N3 0 u3 N3
1/3
 dx dx dx dx  u4 1/3 0
0 0 0 0
( ) ( )
0 0
0 0 1 0
+ u (2/3) − u (1/3) (1.75)
1 0
0 0

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 35 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Particularizando para el elemento e = 2

( )( ) ( ) ( )
0 0 0 0 u1 0 0
0 −3.111 2.9444 0 u2 0.0741 −u 0 (1/3)
+ + (1.76)
0 2.944 −3.111 0 u3 0.0926 u 0 (2/3)
0 0 0 0 u4 0 0

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 36 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Particularizando para el elemento e = 3

Z 1
h iT h i " u1
#
dN3 dN4 dN3 dN4 T u2
− 0, 0, , 0, 0, , + [0, 0, N3 , N4 ] [0, 0, N3 , N4 ] dx
dx dx dx dx u3
2/3 u4
Z 1
T


+ [0, 0, N3 , N4 ] x dx
2/3

T du T du
+ [0, 0, N3 , N4 ] |x =1 |x =1 − [0, 0, N3 , N4 ] |x =2/3 |x =2/3 = 0 (1.77)
dx dx

 (
0 0 0 0 u1
) (
0
)
Z 2/3   Z 2/3
0 0 0 0 u2 N2
− 0 0
dN3 dN3
+ N3 N3
dN3 dN4
+ N3 N4 dx + x dx (1.78)
dx dx dx dx u3 N3
1
 dN4 dN3 dN4 dN4  u4 1/3 0
0 0 dx dx
+ N4 N3 dx dx
+ N4 N4
( ) ( )
0 0
0 0 0 0
+ u (1) − u (2/3) (1.79)
0 1
1 0

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 37 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Particularizando para el elemento e = 3

( )( ) ( ) ( )
0 0 0 0 u1 0 0
0 0 0 0 u2 0 0
+ + (1.80)
0 0 −3.111 2.9444 u3 0.1296 −u 0 (2/3)
0 0 2.944 −3.111 u4 0.1481 u 0 (1)

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 38 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Sistema de Resolución
0.0185 − u 0 (0)
( )( ) ( ) ( )
−3.111 2.9444 0 0 u1 = 0 0
2.9444 −6.2222 2.9444 0 u2 0.1111 0
+ = = 0 (1.81)
0 2.9444 −6.2222 2.9444 u3 0.2222 0
0 0 2.9444 −3.111 u4 = 0 0.1481 + u 0 (1) 0

Solución

( )
u1 =0
u2 = 0.0448
(1.82)
u3 = 0.0569
u4 =0

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 39 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

Técnicas de Aproximación
Método de Galerkin

Gracias por su atención

Preguntas, dudas o aclaraciones?

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 40 / 40

 Técnicas de Aproximación Formulación Débil de una ecuación diferencial

F. P. Beer, E. R. Johnston Jr, D. F. Mazurek, P. J. Cornwell, E. R. Eisenberg, and S. Sanghi.

Vector mechanics for engineers, volume 1.
Tata McGraw-Hill Education, 1972.

Ing. Nelson Lafontaine, Ph. D Método de los Elementos Finitos 40 / 40