EJERCICIOS 12.
En los problemas 1 a 10, resuelva la ecuación de Laplace (1) para una placa rectangular sujeta a
las condiciones de frontera dadas.
PROBLEMA 9
u ( 0 , y )=0 u ( 1 , y )=0
u ( x , 0 )=100 u ( x , 1 )=200
Este problema de valores de frontera tiene la forma del problema 1 que se desarrollo en esta
sección de ejercicios, teniendo en cuenta a=b=1 , f ( x )=100 y g ( x )=200. Por tanto
Tomando en cuenta la Ecuación de Laplace tendremos:
2
∇ u=0
2 2
∂ u ∂ u
2
+ 2 =0
∂x ∂y
0< x < a
0< y< b
-Resolvemos el caso más simple.
2 2
∂ u ∂ u
2
+ 2 =0
∂x ∂y
0< x < a
0< y< b
Sabiendo q tenemos las condiciones:
u ( 0 , y )=0 , 0< y< a
u ( 1 , y )=0, 0< y< b
u ( x , 0 )=100, 0< y< a
u ( x , 1 )=200, 0< y< b
APLICANDO SEPARACIÓN DE VARIABLESS
u ( x , y )= M ( x ) N ( y)
2
∂ u
2
=M ( x ) N(y
∂x
2
∂u
2
=M ( x ) N ( y )
∂y
SUSTITUYENDO EN LA E.D.P. OBTENDREMOS:
M ( x ) N(y) + M ( x ) N ( y )=0
M ( x ) N(y) = -M ( x ) N ( y )
M (x)} over {M(x)} =- {N( y )
=λ
N ( y)
Analizaremos los 3 posibles casos:
λ=0
λ> 0
λ< 0
M (x)} over {M(x)} =λ ^ {N( y )
=λ
N ( y)
*Transcribimos las condiciones iniciales y condiciones de frontera a M ( x ) y N ( y ).
u ( x , y )= M ( x ) N ( y)
*Para:
u ( 0 , y )=0
u ( x , y )= M ( x ) N ( y)
u ( 0 , y )=M ( 0 )∗N ( y ) =0
M ( 0 )=0
*Para:
u ( 1 , y )=0
u ( x , y )= M ( x ) N ( y)
u ( 1 , y )=M ( 1 )∗N ( y )=0
M (1 )=0
Para:
M (x)} over {M(x)} =λ → M left (x right ) N ( y) – M left ( x right ) N left (
→ M (x)} over {M(x)} = ¿
CASO 1:
Si λ=0, entonces:
- M (x) = 0
- M ' (x )=C1
- M ( x ) =C1 x+C 2
M ( x ) =C1 (0)+C 2
C 2=0
M (1 )=C1 ( 1 ) +C 2=0
C 1 ( 1 ) +0=0
C 1=0
CASO 2:
Si λ> 0, entonces λ=∝2
M (x)} over {M(x)} =- {N ( y )
=∝2
N ( y)
2
M (x )= {∝ } ^ {2 } ∗M (x)→ M ( x ) −∝ ∗M ( x )=0
Solucionamos por medio de la Ecuación Característica:
(D −∝ ¿ ¿ 2) x=0 ¿
2
2 2
r −∝ =0
( ∝−r )( ∝+r )=0
∝−r =0 ∨∝+r=0
∝=r ∨∝=−r
∝x −∝ x
M (x) = C 3∗e +C 4∗e
Por uno de los datos dados por el problema nos quedará la siguiente expresión:
(0 ) x −( 0) x
C 3∗e +C 4∗e =0
C 3+ C4 =0 … ( I )
Por otro de los datos dados por el problema nos quedará la siguiente expresión:
(1 ) x − ( 1) x
C 3∗e + C4∗e =0 … ( II )
Ahora despejamos C 3 de la Ecuación ( I ):
C 3=−C 4
Con el dato obtenido reemplazamos en la Ecuación ( II ):
( 1) x − ( 1) x
C 4∗e +C 4∗e =0
C 4 (−e x +e− x ) =0
Sabemos que lo resaltado nunca podrá ser CERO por lo que podemos llegar a la conclusión de:
C 4=0 que a su vez da por consecuencia que C 3=0
CASO 3:
Si λ< 0, entonces λ=−∝2
M (x)} over {M(x)} =- {N ( y )
=∝2
N ( y)
2
M (x )= − {∝ } ^ {2 } ∗M (x)→ M ( x )+∝ ∗M ( x )=0
Solucionamos por medio de la Ecuación Característica:
(D +∝ ¿ ¿ 2) x=0 ¿
2
2 2
r +∝ =0
r =± √ ∝2
r =±i ∝
M (x) = C 5∗cos (∝ x)+C 6∗sen(∝ x)
Por uno de los datos dados por el problema nos quedará la siguiente expresión:
M ( 0 ) = C 5∗cos ( ∝ x)+C 6∗sen (∝ x)
C 5∗cos ( 0 ) +C6∗sen ( 0 )=0
Sabemos que:cos ( 0 )=0 sen ( 0 ) = 1
Con la aclaración previa obtendremos que: C 5=0
Por uno de los datos dados por el problema nos quedará la siguiente expresión:
M ( 1 ) = C 5∗cos ((1) ∝)+C 6∗sen((1)∝)
Como sabemos que C 5=0 nos queda la siguiente expresión:
M ( 1 ) = C 6∗sen(∝)
Deducimos que C 6 no puede ser CERO, otro caso sucedo respecto a sen(∝) ya que para ser
CERO se debe cumplir que ∝ sea múltiplo entero de π .
Por lo que tendremos:
∝=n π , n=1 , 2 ,3 , … .
Ahora nos quedaría:
M n ( x ) =A n∗sen (n π x )
M n : Hace referencia que tendremos una solución para cada uno de los diferentes valores de n
Para:
−N (y)} over {N left (y right )} = {∝} ^ {2 ¿
N (y) = {∝} ^ {2} * N left (y right
N (y) - {∝} ^ {2} *N left (y right ) =
Solucionamos por medio de la Ecuación Característica:
( D −∝ ¿ ¿ 2) x=0 ¿
2
2 2
r −∝ =0
( r −∝ )( r −∝ )=0
∝−r =0 ∨∝+r=0
r =∝∨r =−∝
N ( y ) = C7∗cos h (∝ y )+C 8∗sen h(∝ y)
Por dato inicial dado por el problema tendremos la siguiente expresión:
N ( 0 ) = C7∗cosh ( ( 0 ) ∝) + C8∗senh ( ( 0 ) ∝) =100
Sabemos que: cosh ( 0 )=1 sen ( 0 )=0
Por la aclaración previa y el dato brindado tendremos como resultado la siguiente expresión:
C 7=100
Después de las aclaraciones dadas y el dato hallado tenemos la siguiente expresión:
N ( y ) = C8∗senh (∝ y ) , recordemos que ∝ cumple la condición :∝=n π
N n ( x )=Bn∗sen( nπx)
N n : Hace referenciaque tendremos una solución para cada uno de los diferentes valores de n
Por lo tanto, la respuesta final será:
∞
u ( x , y )=∑ ( A n∗cosh nπy+ Bn senh nπy )∗sen nπy
n=1
Que cumplirá para un:
( )
1 n
1−(−1 )
An =2∫ 100∗sen nπx dx=200
0
nπ
Tendremos también:
[ ]
1
1
B n= 2 ∫ 200∗sen nπx dx −A n cosh nπ
sen nπ 0
[ ( ) ( ) ]
n n
1 1−(−1 ) 1−(−1 )
¿ 400 −200 cosh nπ
sen nπ nπ nπ
[ ]
n
1−(−1 )
¿ 200 [ 2∗csc nπ −coth nπ ]
nπ
Reemplazando los valores hallados tendremos la expresión final:
( ( ) ( [ ] ) )
∞
1−(−1 )n 1−(−1 )n
u ( x , y )=∑ 200 ∗cosh nπy + 200 [ 2∗csc nπ −coth nπ ] senh nπy ∗sen nπy
n=1 nπ nπ