Resolucion mediante el MEF de la ecuacion de

Laplace
Motivacion:
Problema estacionario de conduccion del calor.
Balance de la energa termica:
_
_
Tasa de cambio
de energa termica
en un cuerpo
_
_
=
_
_
Calor generado
por fuentes
en un cuerpo
_
_
+
_
_
Flujo de calor
entrante en
el cuerpo
_
_
Energa termica:
_

c(x)(x)(x, t)dV
c(x) : capacidad calorca
(x) : densidad
(x, t) : temperatura
Calor generado por fuentes:
_

f(x, t)dV
f(x, t) : densidad de calor generado
Flujo de calor entrante:
Figura 1: Dominio de estudio y normal saliente
1

q ndS
q : calor
: contorno de
n : normal saliente de
Balance:
d
dt
_

c(x)(x)(x, t)dV =
_

f(x, t)dV
_

q ndS
Por el teorema de la divergencia:
_

q ndS =
_
V
qdV
Para el caso estacionario:
d
dt
_

c(x)(x)(x, t)dV = 0
Por tanto,
_

f(x, t)dV
_

q(x, t)dV = 0
_

[f(x, t) q(x, t)] dV = 0

Esto se debe cumplirse para cualquier por peque no que sea. Por tanto
q = f en
Ley de Fourier (ecuaci on constitutiva)
q = k
k : coeciente de conductividad termica
Por tanto,
(k) = f en
k = f en
donde, para R
2
,
= =
_

2
x
2
1
,

2
x
2
2
_
(Laplaciano)
Entonces, se obtiene la siguiente expresion:
2
k = f en
Si k = 1 y f = 0,
= 0 en (ecuacion de Laplace)
Problema de valores de contorno
= 0 en (ecuacion de gobierno)
=

en
D
(condiciones de Dirichlet)

n
= () n = q

en
N
(condiciones de Neumann)
Forma debil
Multiplicando ambos miembros de la ecuaci on de Laplace e integrando sobre
el dominio , se obtiene
_

wdV = 0 w
Por otro lado, se tiene que
(w) = w + w
= w + w
Entonces,
w = (w) w
Reemplazando lo anterior, se obtiene
_

wdV =
_

(w)dV
_

w dV = 0 w
_

w dV
_

n (w)dS = 0 w
Se consideran funciones de peso que cumplen w = 0 en
D
y, por tanto,
_

w dV
_
N
w(n )dS = 0 w
_

w dV
_
N
wq

dS = 0 w
donde en la ultima expresion se han considerado las condiciones de contorno de
Neumann. A esta expresion se la conoce como forma debil de la ecuaci on de
Laplace.
3
De la forma debil a la forma fuerte
Se tiene que
_

w dV =
_

(w)dV
_

wdV w
=
_

n (w)dS
_

wdV w
=
_
N
(n )wdS
_

wdV w
Por tanto, teniendo en cuenta que

n
= n ,

wdV +
_
N
_

n
_
wdS
_
N
wq

dS = 0 w
y

wdV +
_
N
_

n
q

_
wdS = 0 w
Considerese el subespacio W
0
= {w : w|

= 0} del espacio de funciones ad-

misibles. Entonces,

wdV = 0 w W
0
Considerese una funci on > 0 en y = 0 en . Adoptese w = W
0
.
Entonces _

()
2
dV = 0
= 0 en
Con lo que se obtiene
_
N
_

n
q

_
wdS = 0 w

n
q

= 0 en
N
Aproximaci on de la forma debil
A partir de
_

w dV
_
N
wq

dS = 0 w
Se construye el problema aproximado
_

w
h

h
dV
_
N
w
h
q

dS = 0 w
h
4
con

h
(x) =
1
N
1
(x) + . . . +
n
N
n
(x)
w
h
(x) = w
1
N
1
(x) + . . . + w
n
N
n
(x)
Se puede ver facilmente que

h
(x) =
1
N
1
(x) + . . . +
n
N
n
(x) =

i
N
i
(x)
w
h
(x) = w
1
N
1
(x) + . . . + w
n
N
n
(x) =

j
w
j
N
j
(x)
Entonces,
_

w
h

h
dV =
_

_
_

j
w
j
N
j
(x)
_
_
_

i
N
i
(x)
_
dV
Defnanse
B
T
=
_

_
N
1
x
1
N
1
x
2
.
.
.
.
.
.
N
n
x
1
N
n
x
2
_

_
; w =
_

_
w
1
.
.
.
w
n
_

_
; =
_

1
.
.
.

n
_

_
Entonces,
w
h
= Bw

h
= B
y
(w
h
) (
h
) = w
T
B
T
B
Deniendo
K =
_

B
T
BdV
se obtiene que
_

w
h

h
dV = w
T
K
5
Espacios de elementos nitos
La malla de elementos nitos
Dado un dominio poligonal se crea una partici on de elementos {
1
,
2
, ...,
e
, ...,
E
},
donde
E : n umero de elementos Esto es,
e

f
= 0 para e = f
E
_
e=1

e
=

,

=
Se considera
e
poligonal.
Se pide que cada lado de
e
sea parte del contorno o lado de alg un otro
elemento. Esto deja fuera a conguraciones de elementos como la mostrada en
la gura (2).
Figura 2: Conguraci on de elementos no permitida
Puntos nodales
Est an ubicadas por lo menos en los vertices del elemento, como se indica en
la gura (3).
Figura 3: Puntos nodales de un elemento triangular de tres nodos
Pero pueden ubicarse en otros puntos.
Existen siempre n nodos: 1, 2, ..., n cuyas posiciones son x
1
, x
2
, ..., x
n
El conjunto de nodos y elementos forman la malla de elementos nitos.
6
Funciones base N
i
1. Son continuas y acotadas
N
i
C(

)
2. Existen n funciones N
i
, una por cada nodo
N
1
, N
2
, ..., N
n
y
N
i
(x) = 0 x
e
si
x
i

e
Esto es, N
i
es cero en los elementos que no est an conectados al nodo i
3. N
i
(x
j
) =
ij
4. N
i
(x) es un polinomio en
e
si x
i

e
Observaci on:
N
(e)
i
(x
e
j
) =
ij
en
e
con N
(e)
i
(x) = N
TL(e,i)
(x)
Aproximaci on de la forma debil con elementos nitos
_

w
h

h
dV =
_
1
w
h

h
dV +. . .+
_
E
w
h

h
dV =
E

e=1
_
e
w
h

h
dV
donde _
e
w
h

h
dV = w
T
K
e

con
K
e
=
_
e
= B
T
BdV
Considerando una aproximacion por elementos nitos, se tiene que, para
e
,
B =
_

_
0 . . .
N
(e)
1
x
1
. . .
N
(e)
n
x
1
. . . 0
0 . . .
N
(e)
1
x
2
. . .
N
(e)
n
x
2
. . . 0
_

_
donde n
(e)
: denota el n umero de nodos de
e
Defnase
K
e
=
_

(e)
= B
T
e
B
e
dV
7
con
B
e
=
_

_
N
(e)
1
x
1
. . .
N
(e)
n
x
1
N
(e)
1
x
2
. . .
N
(e)
n
x
2
_

_
Entonces _

w
h

h
dV = w
T
e
K
e

e
con
w
e
=
_

_
w
(e)
1
.
.
.
w
(e)
n
_

_
;
e
=
_

(e)
1
.
.
.

(e)
n
_

_
Por otro lado,
_

w
h
q

dS =
n
(e)

e
_
e
w
h
q

dS
Donde
e
=
e

_
e
w
h
q

dS = w
T
_
e
Nq

dS
con
N =
_

_
0
.
.
.
N
(e)
1
(x)
.
.
.
N
(e)
n
(x)
.
.
.
0
_

_
Defnanse,
f
(e)
s
=
_
e
N
(e)
q

dS
y
N
(e)
=
_

_
N
(e)
1
.
.
.
N
(e)
n
_

_
Entonces,
_
e
w
h
q

dS = w
T
e
f
(e)
s
en
e
8
Globalmente,
_

w
h

h
dV
_

w
h
q

dS = 0 w
Implica
w
T
(Ku f
s
) = 0 w
Observaci on. Para encontrar K y f
s
hay que ensamblar K
(e)
y f
(e)
s
Elementos triangulares de tres nodos
Figura 4: Elemento triangular de tres nodos
Coordenadas triangulares
Figura 5: Coordenadas triangulares
Observese la gura (5). A partir de ella se puede denir la coordenada tri-
angular asociada al vertice i como

i
=
A

A

_

i
= 0 para los nodos j = i

i
= 1 para el nodo i
9
Figura 6: Graca de la coordenada triangular
Se tiene, entonces, un sistema de coordenadas (
1
,
2
,
3
) En este sistema,
las coordenadas de los nodos son
Nodo 1: (1, 0, 0)
Nodo 2: (0, 1, 0)
Nodo 3: (0, 0, 1)
Por otro lado, las doordenadas del centroide son
_
1
3
,
1
3
,
1
3
_
Las coordenadas
1
,
2
,
3
forman una partici on de la unidad:

1
+
2
+
3
= 1
Interpolacion lineal en R
2
Si f(x, y) vara linealmente, en
e
, entonces f(x, y) = a
0
+a
1
x +a
2
y en
e
.
Por otro lado, existen tres funciones lineales tales que

1
= N
1
(x, y)

2
= N
2
(x, y)

3
= N
3
(x, y)
Defnase, ademas,
f
1
= f(x
1
, y
1
)
f
2
= f(x
2
, y
2
)
f
3
= f(x
3
, y
3
)
De lo anterior se sigue que
f(x, y) = f
1
N
1
(x, y) + f
2
N
2
(x, y) + f
3
N
3
(x, y)
o, de manera equivalente,
f(x, y) =
_
f
1
f
2
f
3
_
_
_
_
N
1
(x, y)
N
2
(x, y)
N
3
(x, y)
_
_
_
10
Transformaci on de coordenadas
Se tiene que

1
+
2
+
3
= 1
Ademas,

1
x
1
+
2
x
2
+
3
x
3
= x

1
y
1
+
2
y
2
+
3
y
3
= y
Por tanto,
_
_
_
1
x
y
_
_
_
=
_
_
1 1 1
x
1
x
2
x
3
y
1
y
2
y
3
_
_
_
_
_

3
_
_
_
Se puede demostrar que
2A = det
_
_
_
_
1 1 1
x
1
x
2
x
3
y
1
y
2
y
3
_
_
_
_
El cambio de coordenadas inverso es
_
_
_

3
_
_
_
=
1
2A
_
_
x
2
y
3
x
3
y
2
y
2
y
3
x
3
x
2
x
3
y
1
x
1
y
3
y
3
y
1
x
1
x
3
x
1
y
2
x
2
y
1
y
1
y
2
x
2
x
1
_
_
_
_
_
1
x
y
_
_
_
De aqu,

i
x
=
N
i
(x, y)
x
= y
j
y
k
j = k, i = k, i = k

i
y
=
N
i
(x, y)
y
= x
k
x
j
j = k, i = k, i = k
Se puede entonces plantear las siguiente aproximacion lineal para la inc ognita
dentro de un elemento:

h
(x) = u
(e)
1
N
(e)
1
(x) + u
(e)
2
N
(e)
2
(x) + u
(e)
3
N
(e)
3
(x) en
e
= u
(e)
1

1
+ u
2

2
+ u
(e)
3

1
en
e
En este caso, la matriz B
(e)
tiene la siguiente forma:
B
(e)
=
_

_
N
(e)
1
x
N
(e)
2
x
N
(e)
3
x
N
(e)
1
y
N
(e)
2
y
N
(e)
3
y
_

_
=
1
2A
_
y
2
y
3
y
3
y
1
y
1
y
2
x
3
x
2
x
1
x
3
x
2
x
1
_
Por tanto,
K
(e)
=
_
e
_
B
(e)
_
T
B
(e)
dV =
_
B
(e)
_
T
B
(e)
A
e
donde A
e
es el area del elemento. Se obtiene nalmente que
_
e
w
h

h
dV =
_
w
(e)
_
T
K
11