SOLUCIN DE LA ECUACIN DE LA PLACE EN UN RECTNGULO
(completo)
Considrese la ecuacin elptica de La Place:
u
xx +
u
yy =O
En un rectngulo O < x < I y O < x < I , con las condiciones de frontera
siguientes:
u(0 , y) = u(l , y) = O
u(x , 0) -/(je) y
u
(x , i] = g(x)
Esta ecuacin se resolver pur el procedimiento de separacin de variables, o sea que se propone
una solucin de a forma:
u (x , y) = M(x)W(y)
Derivando con respecto a x:
du ,. ,
=N(y)M\x)
ax
Volviendo a derivar:
Derivando respecto a y:
dU ,
Volviendo a derivar:
Sustituyendo en !a ecuacin diferencia!:
tf(y)M"Cr) + M(x)Ar"'(y) = O
" (y) _ n
Osea:
N"(y) =
Finalmente se obtiene la ecuacin diferencial para W(y):
N"(y) - AN( y) = O
Para M(x) se tiene:
"
M"(x) + AM(x) = O
Solucin de las ecuaciones:
Como la ecuacin para N (y) es una ecuacin lineal, se resolver de la manera clsica, calculando
i3s raices uci poinofmo caracterstico:
r
2
--1 = O cuyas races son r
x
=+ VI y r
2
= VA
As que la solucin general es:
/v'y) = aexpi vA yj + exp^MA y)
Donde a y b son constantes por determinar.
En cuanto a la ecuacin:
M"(x)+ M(x) = O
Se reconoce inmediatamente que tiene ia forma de ia ecuacin dei oscilador armnico cuya
solucin es:
MOO = Acos(JIx) + B sent^x}
Las condiciones de frontera para xson:
u(0 , y) -u( , y) =O
Entonces:
u(0 , y) = N(y)M(0) = O
Lo que implica que:
Af( 0) = ylcos(O) + Bsen(O) = O
2.
(y y) + b
k
cosh(y y)
Sin embargo, afn de evaluar las condiciones de frontera conviene expresar las funciones N
fe
(y)
como una suma de funciones seno. Empleando identidades para las funciones hiperblicas es
posibles expresar tales funciones como:
W
fe
(y) = a
h
senh (~y) + b
k
senh (y (y - )J
De esta forma, las soluciones para la ecuacin de Laplace, u
k
(x , y), se escriben como:
u
k
(x , y) = [a
k
senh [~yj + b
k
senh [y (y - )J] Jysen(yx)
La solucin formal general para la ecuacin de Laplace se obtiene como una superposicin lineal
Las constantes a
k
y b
k
se determinan a partir de las condiciones de frontera :
u(x, 0) -/(*) y u(x , 1) = g(x)
O sea:
--- wk el 2 fnk \
i(r,0) = V^
=
^h^senh I -- l /-sen \-r X ) =
u(x,T) = Z?= ia
k
senh [^ ] sen (^ x) =
5
W
A fin de determinar las constantes a
k
y b
k
se utiliza la ortonormalidad de las eigenfunciones de
M(x) o seaque (<p
k
(x) , <pj(x)) = O 517 ^ /cyesl si j = /c : as:
-T
f
]
sen
(T /o
se
" (T
f*
sm
(T *)
5en
(T
Entonces, usando la ortonormalidad se tiene:
b
k
senh [- y 1] = fj(x)sen (y x) dx
Osea:
3
Como cos(O) = 1 , esto implica A = 0. Como sen (0) = O, entonces B ^0. La constante
8 se asimila a las constantes de N(k~), as que la solucin tendr la forma:
u(x , y) = /V(y)sen (VA x]
Considrese ahora la condicin de frontera :
u(l , y) = JV(y)M(/) = JV(y) sen (VI ) = O
Lo que implica que sen (VA /) = O, condicin que se cumple si:
VA l = nk , k = 1, 2 , . . .
Esta ecuacin produce e conjunto de eigenvaiores {/t}:
= j- sen
v'
Las eigenfunciones para M(x), ya normalizadas sern:
I i
x
)
Y en consecuencia, las eigenfunciones para N(y) sern:
(nk \ (nk \
N
k
(y) = a
k
exp y + b
k
exp! y 1
\ I / \L /
Obsrvese que puesto que:
senft(x) = --^ y
2cosh(x) = e
x
+e~
x
Sumando miembro a miembro:
2 [senh (x} + cosh(x)] = 2 e
x
e
x
= coshO) +sen/i (z)
Y restando miembro a miembro:
2[cosh{ - senh(x)] = 2 e~
x
e
~
x
= cosh(x) + sen/i(x)
Las funciones N
k
(y) se pueden escribir corno:
De igual manera se obtiene que:
a
fc ~ u(nk \
sent \-j- J
Observaciones.
La funcin u(x , y) un desplazamiento de estado continuo o de equilibrio de una membrana
estirada, como la de un tambor, sujeta por medio de un marco rectangular que sigue un
movimiento de acuerdo a ias condiciones de frontera que se fijaron en este problema y que se
repitan aqu como recordatorio:
u(0,y) = u(/, .y) =O
u(x, 0) =/(x) y u(x,) = g(x)
Las condiciones impuestas en x = O y x = I fueron determinante para calcular el conjunto de
eigenvalores. Sin embargo parece que tales condiciones son un poco artificiales y particulares. Se
pueden estables condiciones de frontera diferentes de cero en esos puntos por medio de dos
conjuntos de condiciones de frontera. El primer conjunto seria como ei anterior:
u(0,y) -w(/, y) =0
i(*,0) = /(*) y u(x,t) = g(x)
Con el que se obtendra un conjunto de eigenvalores y una serie de Fourier igual a ios obtenidos
aqu.
Y otro conjunto de condiciones de frontera igual a! siguiente:
u(0, y) = /i(y) y u( /, y) =m(y)
u(; c,0) = u( x, F ) = O
Con el que obtendra un conjunto de eigenvalores y una serie de Fourier diferente a la anterior,
pero con ia unin de ambos conjuntos de eigenvalores y !a suma de las dos series se obtendra una
solucin a un problema que tuviera condiciones arbitrarias de frontera, aunque el problema
tendra como dominio de definicin de la solucin, un cuadrado de la forma : O < x < I y
O < y < .
5
