UNIDAD PROFESIONAL INTERDISCIPLINARIA DE BIOTECNOLOGÍA
Aplicaciones Matemáticas (Taller)
Proyecto - Parcial 3
Equipo: 6
Cervantes Conde Estefanía
De la Cruz Zúñiga Fernanda Paola
González Álvarez Diana Karina
Guerrero Zacarías Fernanda
Gutiérrez Garay Braian Gael
Grupo: 3BM1
Fecha de entrega: 09/01/2023
Problema 1
Encuentre el desplazamiento u(x , t) en una cuerda elástica que está fija en sus
extremos y se pone en movimiento punteándola en su centro.
2 2
2 d u d u
a 2
= 2 , 0< x< L, t >0
dx dt
u ( 0 , t ) =u ( L, t )=0 , t ≥ 0 ,
u ( x , 0 )=f ( x ) , 0≤ x ≤ L ,
du
| =u ( x , 0 )=0 ,0 ≤ x ≤ L
dt t =0 t
Ax , 0 ≤ x ≤ L /2
f ( x )=
A ( L−x ) , L/2< x ≤ L
La solución está dada por la forma:
[ ]
∞
nπa nπa nπ
u ( x , t )=∑ A n cos t+ B n sen t sen x
n=1 L L L
Donde:
L
2 nπ
An = ∫ f ( x ) sen x dx
L0 L
L
2 nπ
B n= ∫
nπa 0
g ( x ) sen
L
xd x
Para obtener Bnidentificamos las condiciones iniciales respecto a “t” dadas por el
problema
du
| =u ( x , 0 )=0 ,0 ≤ x ≤ L
dt t =0 t
Donde:
ut ( x , 0 )=g ( x )=0
L
2 nπ
B n= ∫
nπa 0
0∙ sen
L
x dx=0
Bn=0
Obteniendo An
L
2 L
2 nπ 2 nπ
An = ∫ Ax sen x dx+ ¿ ∫ A( L−x)sen x dx ¿
L0 L LL L
2
Resolviendo la primera integral
u=Ax
nπ −L nπx
dv =sen x dx du= A dx v= cos
L nπ L
[ ]
L
2
2 −ALx
L nπ
cos
nπx − AL
L
−
nπ ( )∫ cos nπxL dx
0
L nπ [
2 −ALx
cos
nπx AL L
+
L nπ nπ
sen
L ( )
nπx L/2
0 ]
[ −2 Ax
nπ
cos
nπx 2 AL
+ 2 2 sen
L n π
nπx L/2
L 0 ]
[ −2 A(L /2)
nπ
cos
nπ (L/2) 2 AL
L
+ 2 2 sen
n π
nπ (L /2)
L
−[ 0+0 ]
]
−AL nπ 2 AL nπ
cos + 2 2 sen
nπ 2 n π 2
2 AL nπ AL nπ
sen − cos
2 2
n π 2 nπ 2
Resolviendo la segunda integral:
u=AL− Ax
du=− A dx
dv =sen ( nπxL ) dx
−L nπx
v= cos
nπ L
[ ]
L
2
L
AL−Ax
−L
nπ
cos ( )
nπx AL
L
−
nπ ( )∫ cos nπxL dx
L
2
L
[ −2 AL
nπ
cos
nπx 2 Ax
L
+
nπ
cos
nπx 2 AL
L
− 2 2 sen
n π
nπx
L
L
2
]
[ −2 AL
nπ
cos nπ +
2 AL
nπ
2 AL
cos nπ− 2 2 sen nπ −
n π
−2 AL
nπ ][ nπ AL
cos +
2 nπ
nπ 2 AL
cos − 2 2 sen
2 n π
nπ
2 ]
sen nπ=0
[
[ 0−0 ] − − AL cos nπ − 22AL2 sen nπ
nπ 2 n π 2 ]
AL nπ 2 AL nπ
cos + 2 2 sen
nπ 2 n π 2
Por lo tanto…
2 AL nπ AL nπ AL nπ 2 AL nπ
An = sen − cos + cos + 2 2 sen
2 2
n π 2 nπ 2 nπ 2 n π 2
4 AL nπ
An = sen
2 2
n π 2
Sustituyendo los valores obtenidos en la solución
[( ) ]
∞
4 AL nπ nπa ( ) nπa nπ
u ( x , t )=∑ 2 2
sen cos t+ 0 sen t sen x
n=1 n π 2 L L L
∑ [( ) ]
∞
4 AL nπ nπa nπ
u ( x , t )= 2 2
sen cos t sen x
n=1 n π 2 L L
Problema 2
Investigar una aplicación de ecuaciones diferenciales parciales en tu carrera,
resolverla y explicar que usos tiene y su relevancia.
Aplicación
Alumnos del IPN crearon bolsas a partir de un bioplástico fabricado con
pseudotallo de M. acuminata AAA y M. sapientum ABB
(plátano). Al ser bolsas, requieren ser resistentes, por lo
tanto, para probar su calidad, piezas cuadradas de este
bioplástico se sometieron a una serie de pruebas.
Suponiendo que la tela se encuentra situada como en la
figura 12.5, y que las aristas están fijas y son de longitud
unitaria, la ecuación diferencial parcial que describe el
movimiento de las vibraciones está dada por:
( )
2 2 2
∂ Z 2 ∂ Z ∂ Z
2
=a 2
+ 2 … (1)
∂t ∂x ∂ y
Donde:
Z ( x , y , t ) es el desplazamiento de cualquier punto ( x , y ) de la tela desde su
posición de equilibrio, es decir, el plano xy en cualquier tiempo
2
a es una constante que depende de la tensión y la densidad de la tela.
El hecho de que las aristas estén fijas implica tener cuatro condiciones de frontera
Z ( 0 , y , t )=0 Z ( 1 , y , t )=0 Z ( x , 0 ,t )=0 Z ( x ,1 , t )=0
El hecho de darle a la tela una forma inicial específica conduce a:
Z ( x , y , 0 ) =f ( x , y )
Finalmente, el hecho de soltar la tela después de haberle dado esta forma nos
dice que
Z t ( x , y ,0 )=0
Estas condiciones se cumplen por supuesto para 0< x <1 , 0< y <1 , t> 0.
Asumiendo que (1) tiene una solución separable de la forma Z=XYT ,donde
X , Y ,T son funciones de x , y , t ,respectivamente. Entonces sustituyendo en (1) se
obtiene:
'' '' ''
T❑ X❑ Y ❑
XY T '' =a2 ( X ' ' YT + X Y ' ' T ) o 2
= +
a T X Y
Después de dividir por a 2 XYT .
Puesto que un lado es una función solo de t, mientras que el otro lado es una
función de x y y, sigue que cada lado debe ser una constate, la cual denotamos
por −λ 2, esto es
'' '' ''
T❑ 2 X❑ Y ❑ 2
2
=− λ , + =−λ …... (2)
a T X Y
Si ahora escribimos la segunda ecuación de (2) en la forma
( )
'' ''
X❑ Y❑ 2
=− +λ
X Y
Vemos que cada lado debe ser una constante, la cual denotamos por −μ2 .
Entonces
'' ''
X❑ 2 Y❑ 2
=−μ , =−ν , …..(3)
X Y
Donde ν 2=λ2−μ 2, así que λ 2=μ2 +ν 2
De la ecuación (3) y la primera ecuación en (2) tenemos entonces
X + μ X=0 , Y + ν Y =0 , T +a ( μ + ν ) T =0
'' 2 '' 2 '' 2 2 2
Notamos que nos salen raíces complejas al realizar: (k toma los valores de µ, v y
a √ μ +v
2 2
2 2
r +k =0
❑ 2
r =−k
r =±i k
Y sabemos que si las raíces son complejas, entonces:
y=c1 y 1 +c 2 y 2
α+ β i
y=c1 x +c 2 y 2
α+ βi α −βi
y=c1 x + c2 x
iθ
e =cosθ +isenθ
α
y=x ¿
De modo que X =c 1 cosμx+ c 2 senμx ,
Y =c 3 cosνx+ c 4 senνx ,
T =c 5 cos a √ μ2 +ν 2 t+ c6 sena √ μ 2+ ν 2 t
De la primera y tercera condiciones de frontera encontramos c 1=0 , c 3=0 de modo
que
Z=c 2 c 4 senμxsenνy ( c 5 cos a √ μ 2+ ν 2 t+c 6 sena √ μ2 + ν 2 t )
De la segunda y cuarta condiciones de frontera encontramos
μ=mπ , ν=nπ , donde m , n=1 , 2 ,3 , 4 …
De modo que Z=c 2 c 4 senmπxsennπy ( c5 cos a √ m2 +n2 πt +c 6 sena √ m2 +n2 πt )
La condición Z t ( x , y ,0 )=0 conduce a c 6=0 de modo que
Z=Bsenmπxsennπy cos a √ m2 +n 2 πt …..(4)
Donde B=c2 c 4 c 5
Para satisfacer Z ( x , y , 0 ) =f ( x , y ) tendremos que usar el principio de superposición,
esto es, la sumatoria de soluciones de la forma (4) sobre todos los valores enteros
de m y n .Esto conduce a la solución de serie doble
∞ ∞
Z= ∑ ∑ Bmn senmπxsennπy cos a √ m2 +n2 π t
m=1 n=1
Donde hemos reemplazado Ben (4) por Bmn puesto que cada solución puede tener
un coeficiente diferente. Esto satisfará Z ( x , y , 0 ) =f ( x , y ) si
∞ ∞
f (x , y )= ∑ ∑ Bmn senmπxsennπy …..(5)
m=1 n=1
Podemos escribir la serie doble (5) en la forma
(∑ )
∞ ∞
f (x , y )= ∑ Bmn senmπx senmπx
m=1 n=1
esto es,
∞
f (x , y )= ∑ b m senmπx …..(6)
m=1
donde
∞
b m=∑ Bm sennπx …..(7)
n=1
Puesto que (6) representa una serie seno de Fourier en x tenemos por los
métodos acostumbrados
1
2
b m=
10
∫ f ( x , y ) senmπxdx…(8)
Puesto que (7) representa una serie de Fourier en y, tenemos similarmente
1
2
Bmn= ∫ b m sennπy dy …(9)
10
Al sustituir 8 en 9 obtenemos:
1 1
Bmn=4∫ ∫ ( f ( x , y ) senmπxdx ) sennπy dx dy …(10)
0 0
Entonces, la solución es:
∞ ∞
Z= ∑ ∑ Bmn senmπxsennπy cos a √ ( m2 +n2 ) πt
m=1 n=1
1 1
Donde Bmn=4∫ ∫ ( f ( x , y ) senmπxdx ) sennπy dx dy
0 0
Interpretación
Para examinar el posible significado físico de los términos en la serie asumimos
que f(x,y) es tal que todos los términos de la serie son cero, excepto aquel para el
cual m=3 y n=3. Así, este término será:
Z=B 33 sen 3 πx sen 3 πy cos π a √ ( 3 +3 ) t
2 2
Para todos los valores de tiempo t, el desplazamiento será 0 a lo largo de las
líneas x=1/3, 2/3; y=1/3, y=1/3, 2/3 (líneas nodales).
En la siguiente gráfica se indica el movimiento en una dirección con los cuadros
sombreados y los cuadros en blanco indican otra dirección.
Referencias:
Spiegel, M. R. (1996). Ecuaciones Diferenciales Aplicadas. Prentice Hall.
Jiménez Coello, C. (2017). Evaluación de polímeros en pseudotallos de Musa
acuminata.
AAA, Musa sapientum ABB y Musa paradisiaca AAB para elaboración de
bioplástico [Tesis de maestría]. El Colegio de la Frontera Sur.
