UNIVERSIDAD NACIONAL DE SAN AGUSTÍN
Escuela de ing. Civil
Asignatura: Calculo IV
Tema: Ecuaciones Diferenciales en Vigas
Integrantes:
-Ahuate Toribio Carlos Daniel
-Ayma Choquehuayta Alfredo
-Asto Maldonado Gustavo
-Cornejo Sulla Alexis
-Lacastro Ancco Carlos
-Lazo Quispe Oscar
-Quinta Estrada Fiorella
Arequipa,30 de octubre del 2019
Objetivos
Objetivo General
Aplicar ecuaciones de diferenciales en movimientos oscilatorios
Objetivo Especifico
-Aprender a identificar la ecuación diferencial para cada caso tipo de movimiento amortiguado,
MAS, Movimiento Amortiguado y Movimiento Vibratorio Forzado.
-Aplicar los métodos estudiados para la resolución de ecuaciones diferenciales planteadas en
este tipo de movimientos.
Movimiento Armónico Simple
1) Un cuerpo que pesa 64Lb sujeto al extremo de un resorte lo estira 0.32 pies. Cuando
el cuerpo está a 8 plg. Sobre la posición de equilibrio se le comunica una velocidad
dirigida hacia debajo de 5 pies/seg.
a) Halle la ecuación del movimiento.
Peso (w): 64
S: 0.32
G: 32pie/seg²
w 64
K= k= =200
s 0.32
Ecuación diferencial de un movimiento armónico
d ²x k
+ x =0
dt ² m
Reemplazando tenemos:
200
x + x=0
64
32
x +100 x=0
−2
Con un X(0)= X ( 0 )=5
3
K
W ²=
M
200
W=
√ 64
32
=10
Entonces la solución general es
x ( t )=C 1 Coswt +C 2 senwt
Reemplazando W=10 tenemos:
x ( t )=C ₁cos 10 t +C ₂ sen 10 t
−2
Como x(0) =
3
−2
=C ₁ cos ( 0 ) +C ₂ sen (0)
3
−2
=C ₁+0
3
-2/3=C₁
Como x`(0)=5
x ( t )=−10 C 1 sen 10 t+10 C 2 cos 10 t
5=10C ₂
C₂=1/2
Reemplazando en la ecuación tenemos:
x ( t )=−2/3 cos 10t +1/2 sen 10 t
4 1
x ( t )=
√ + sen(10 t+ φ)
9 4
−2/3
Pero φ¿ arctg φ=5.3558
1/2
5
sen(10 t +φ)
6
5
sen(10 t +5.3558)
6
b) Calcular la Amplitud, Periodo y Frecuencia.
2π
T=
W
2π
T=
10
π
T = segundos
5
6
A=√ 2 c ₁2+ c ₂2 A=
5
Frecuencia:
1
F=
2π
10
Periodo:
π
T = segundos
5
c) Cuantas oscilaciones completas ha dado:
Entonces en un ciclo hay un segundo por lo tanto
5
(3π)= 15 ciclos por segundo.
π
d) En que instante pasa el cuerpo por la posición de equilibrio en dirección hacia
abajo por segunda vez?
Si x(t)=0 , entonces
5
sen(10 t +5.3358)=0
6
Esto solo si:
(10 t+5.3358)=n π , n
Si
( 10 t+5.3358 )=0
Entonces
−5.3358
t= < 0 , tiempo negativo no existe
10
Si
( 10 t+5.3358 )=π
Entonces:
−π −5.3358
t= <0
10
Si.
( 10 t+5.3358 )=2 π
Entonces:
−2 π −5.3358
t= 0.09
10
Si:
( 10 t+5.3358 )=3 π
Entonces:
−3 π −5.3358
t= =0.407
10
Si:
( 10 t+5.3358 )=4 π
Entonces:
−4 π −5.3358
t= =0.721
10
Como el cuerpo pasa por la posición de equilibrio por segunda vez y en
dirección hacia abajo , el tiempo que estamos buscando es t 3.
e) En que instante alcanza el cuerpo sus desplazamientos extremos hacia uno u otro
lado de la posición de equilibrio.
Como x(t)=5/6sen (10t+5.3558) es la ecuación del movimiento, entonces la
velocidad instantánea es x´(t)= 50/6 cos(10t+5.35558).
Los desplazamientos extremos se dan cuando cos(10t+5.3558)=0 si y solo si
10t+5.3558=π/2
Los diferentes valores están dados por:
Si: n=0,
10t+5.3558=π/2
π
−5.3558
Entonces: 2
t 0= <0
10
Si n=1 , 10t+5.3558=π+π/2
Entonces:
3π
−5.3558
2
t 0= <0
10
Si n=2 , 10t+5.3558=2π
Entonces:
5π
−5.3558
2
t 0= =0.2499
10
Si n=3, 10t+5.3558=3π
Entonces:
7π
−5.3558
2
t 0= =0.564
10
Por lo tanto ,los tiempos en los que el cuerpo alcanza sus desplazamientos
( 2 k +1 ) π−10.7116
extremos son:t= para valores k=2,3….
20
f) t=3
5 5
x ( t )= sen ( 10 t+ 5.3558 )= sen ( 10 ( 3 ) +5.3558 )=−0.597 pies el signo indica que
6 6
está por encima de la posición de equilibrio
5 5 −5.816 pies
v ( t )=x ' ( t )= ∗10 cos ( 10 t+ 5.3558 )= ∗10 cos ( 10 ( 3 ) +5.3558 ) =
6 6 seg
−5 −5
a ( t )= ∗102 sen ( 10 t+ 5.3558 )= ∗102 sen ( 10 ( 3 )+5.3558 )=59.67 pies/se g 2
6 6
g) x ( t )=0→ t=0.09 seg v ( t ) =?
50 50
v ( t )=x ' ( t )= cos ( 10 t+ 5.3558 )= cos ( 10 ( 0.09 )+5.3558 )=8.33 pies/seg
6 6
5
h) x ( t )=5 plg= pies
12
5 5
x ( t )= sen ( 10 t+ 5.3558 )=
6 12
10 t+5.3558=arcsen ( 12 )
π 5π
10 t+5.3558= +2 nπ o 10 t+5.3558= + 2nπ
6 6
π 5π
t=
( 6
+2 nπ −5.3558)
o t=
6
+2 nπ −5.3558
10 10
i) Sube cuando el tiempo es
5π
+ 2nπ −5.3558
6
t=
10
Problema de movimiento vibratorio con amortiguación
Una masa de 1 slug estira a un resorte en 1,28 pie, una fuerza de
amortiguamiento numéricamente igual a β veces la velocidad instantánea que
actúa sobre el contrapeso. Deduzca la ecuación del movimiento si la masa se
suelta a 1 pie por debajo de la posición de equilibrio con una velocidad hacia
abajo de 1 pie por segundo cuando β = 8, 10, 12, trace la gráfica en cada caso.
Solución
∂2 x ∂x
m 2 =−kx−β
∂t ∂t
Para aplicar la ecuación debemos determinar k.
Como F=ks (ley de Hooke) y F=mg,
Por lo que ks=mg,
mg
Entonces k = ;
s
De donde k=32/1.28; entonces
k=25
Reemplazando en la ecuación
∂2 x ∂x
2
=−25 x−β
∂t ∂t
∂2 x ∂x
2
+β +25 x=0
∂t ∂t
Con x(0)=1, x'(0)=1 ; la ecuación auxiliar es:
r 2 + βr +25=0
Cuya solución:
−β ± √ β−100
r=
2
Caso I
Cuando β=8, la ecuación auxiliar
r 2 +8 r +25=0
Tiene raíces:
r =(−4 ± 3i)
Es un caso subamortiguado y la ecuación de movimiento tiene la forma
x ( t )=C 1 e−4 t cos 3 t+C 2 e−4 t sin 3 t
x ' ( t )=C1 (−4 e ¿ ¿−4 t cos 3 t−3 e−4 t sin 3 t)+C 2 ¿ ¿)
Como x(0)=1 y x'(0)=1;
Reemplazando C 1=1,
−4 C1 +3 C 2=1
5
Entonces C 2=
3
Ahora expresamos x(t) en su forma alternativa y para ello tenemos
34
A=√ C 21 +C22 = √
3
C1 3
3 tan θ= =
C2 5
−1 3
θ= tan = 0.5404
5
De donde θ es un ángulo del primer cuadrante.
Por lo tanto la solución es:
34 −4 t
x ( t )= √ e sin(3 t +0.5404)
3
Caso II: CRITICAMENTE AMORTIGUADO
De la forma x ( t )=e−⋋ t (C 1+C 2 t)
β=10
La ecuación r 2 +10 r +25=0
⋋ 2−w 2=0
−β 1 2
r 1 , r 2=−5 r= ± √ β −100
2m 2
Esto es por que las fuerzas de rozamiento son iguales a la fuerza ejercida por el
resorte
x ( t )=e−5 t (C1 +C 2 t)
Condiciones del problema
X(0)=1
X´(0)=1
Reemplazando
C1=1 Y C2=6
x ( t )=e−5 t (1+6 t)
Caso III
Cuando 12 , el sistema esta sobreamortiguado debido a que las raíces de la
ecuación auxiliar son r 6 11 ,luego la ecuación del movimiento es
x(t ) C1e ( 6 11) t
C 2 e ( 6 11)t
,con x(0) 1 , x '(0) 1
C1 C2 1
Con estas condiciones se tiene el sistema (6 11)C1 (6 11)C2 1
Resolviendo:
C1 C2 1..........................................(1)
(6 11)C1 (6 11)C2 1........(2)
De (1) se obtiene C2 1 C1
Reemplazando en (2)
(6 11)C1 (6 11)(1 C1 ) 1
( 11 6)(1 C1 ) ( 11 6)C1 1
( 11C1 6C1 6 11) ( 11C1 6C1 1
2 11C1 6 11 1
2 11C1 7 11
7 11
C1
2 11
11 7 11
C1
22
Reemplazando en C2
11 7 11
C2 1
22
11 7 11
C2
22
Los valores de C1 y C2 son
11 7 11 11 7 11
C1 C2
22 , 22
Reemplazando en la ecuación de movimiento:
(11 7 11) ( 6 (11 7 11) ( 6
x (t ) e 11) t
e 11) t
22 22
Movimiento Vibratorio Forzado
8
1. Un cuerpo que pesa 16lb estira un resorte de
pie. Inicialmente, el peso parte del
3
resorte desde un punto que está 2 pies debajo de la posición de equilibrio y el
movimiento posterior se realiza en
un medio que opone una fuerza de amortiguación numéricamente si el peso es
impulsado por una fuerza exterior igual a f(t)=10cos(3t)
Se trata de un movimiento vibratorio amortiguado en donde el peso es impulsado por una
d 2 x β dx k f ( t)
Fuerza exterior, cuya ecuación diferencial es 2
+ + x= . Para usar está
d t m dt m m
ecuación diferencial, hallemos los valores de k, m y β.
8 lb w 16 1
Por la ley de Hooke 16=k , entonces k=6 . Como m= entonces m= =
3 pie g 32 2
1
slug. Por dato del problema β= .
2
Sustituyendo estos resultados en la ecuación diferencial del movimiento tenemos
d 2 x dx '
2
+ +12 x=20 cos (3 t ) , con x ( 0 ) 02 , x ( 0 )=0
d t dt
I) Hallo la ecuación auxiliar de la parte homogénea:
d 2 x dx
+ +12 x=0
d t 2 dt
m2 +m+ 12=0
−1 ± √ 12−4∗1∗12
m=
2∗1
−1 ± √ 47 i
m=
2
Por lo tanto la solución complementaria es:
−t
x ( t )=e 2
( c 1∗cos ( √247 t )+c sen( 472 t ))
2
Por el método de los coeficientes indeterminados, suponemos una solución particular
del tipo
x p ( t )= Acos ( 3 t )+ Bsen ( 3 t ) , en tal caso , sustituyendo
x 'p ( t )=−3 Asen ( 3 t ) +3 Bcos ( 3t )
x ''p ( t ) =−9 Acos (3 t )−9 Bsen ( 3 t ) , luego se tiene
x ''p + x 'p +12 x p= (3 A+ 3 B ) cos ( t )+ (−3 A +3 B ) sen ( 3 t )=20 cos ( 3t )
3 A +3 B=20(ii)
Esto nos da {
−3 A +3 B=0(i)
20−3 B
(i) A=
3
20−3 B
Reemplazo (i) en (ii) -3*( )+3B=0
3
-20+3B+3B=0
10
B=
3
3∗10
20− 10
(i)A= 3 =
3
3
10 10
, cuya solución es A= , B= . Entonces la solución general es
3 3
−t
(
x ( t )=e 2 c 1∗cos ( √247 t )+c sen( 472 t ))+ 103 (cos ( 3t ) + sen(3 t))
2
−4
Como x(0)=2 , entonces c 1= , evaluando x ' ( 0 )=0 en
3
−t
'( t ) 2
x =e ¿
−64
Obtengo que c 2= . Por lo tanto laecuación de movimiento es
3 √ 47
−t
'( t )
x =e 2
( (
−4
3
) cos ( √247 t )+( 3−64
√ 47
47
2
10
)sen ( t ) ) + ( cos ( 3 t )+ sen ( 3 t ) )
3
2) Una masa de 1 slug se sujeta a un resorte cuya constante es 5lb/pie. Inicialmente, la
masa se sujeta desde un punto que está 1 pie abajo de la posición de equilibrio con
una velocidad dirigida hacia debajo de 5pies/seg y el movimiento posterior se realiza
en un medio que opone una fuerza de amortiguación numéricamente igual a dos veces
la velocidad instantánea.
a) Encuentre la ecuación del movimiento si la masa es impulsada por una fuerza
exterior igual a f ( t )=12 cos 2 t+3 sen 2 t
b) Trace la gráfica de la solución transitoria y de la solución estacionaria usando
el mismo sistema de ejes para ambas gráficas.
c) Trace la gráfica de la ecuación del movimiento.
Solución:
Para encontrar la función que describa el movimiento de la masa sujeta al resorte
debemos de resolver la siguiente ecuación diferencial:
d 2 x β dx k f (t)
2
+ + x= (1)
d t m dt m m
Dónde:
f ( t ) : fuerza externa
m: masa delbloque
k :constante de rigidez del resorte
β : constante de amortiguacion
a) Con la información del problema obtenemos :
5 lb
m=1 slug ; k = ; β=2 y f ( t )=12 cos 2t +3 sen 2t
pie
Reemplazamos estos valores en la ecuación (1)
d 2 x 2 dx 5 12 cos 2 t+3 sen 2t
2
+ + x=
d t 1 dt 1 1
Tenemos lista la ecuación para resolverla:
1.- Primero hallamos la solución complementaria, para ello tomamos la
ecuación homogénea:
d2 x dx
2
+ 2 +5 x=0
dt dt
2.- planteamos la ecuación:
p ( r )=r 2 +2 r +5
Cuyas raíces son:
r 1=−1+2 i y r 2=−1−2i
De aquí obtenemos la solución complementaria:
x ( t )=C 1e (−1+2 i) t +C 2 e(−1−2 i)t
Pero la solución la transformaremos aplicando la siguiente:
e iθ =cosθ +isenθ
Quedando como resultado:
x ( t )=e−t (C 1 cos 2 t+C 2 sen 2t )
Para hallar la solución particular aplicaremos el método de coeficientes
indeterminados, para ello suponemos que la solución particular es:
x ( t )= Acos 2 t+ Bsen 2t
Como solución de la ecuación no homogénea, la aplicamos en la ecuación
diferencial
d2 x dx
2
+ 2 +5 x=12 cos 2 t+3 sen 2t ¿
dt dt
Para hallar las constantes A y B derivamos x(t):
x ´ (t)=−2 Asen 2t +2 Bcos 2 t
x ´ ´ ( x )=−4 Acos 2t−4 Bsen 2t
Reemplazando en (**) obtenemos:
( A+ 4 B ) cos 2t +(−4 A+ B)sen 2 t=12cos 2t +3 sen 2 t
De aquí obtenemos el sistema:
A+ 4 B=12 y−4 A + B=3
A=0 y B=3
Entonces: x ( t )=e−t (C 1cos 2 t +C 2 sen 2 t )+3 sen 2 t
Para hallas las constantes aplicamos condiciones iniciales:
Para C1 x ( 0 )=1 ; C1=1
0
x ( 0 )=e (C 1cos 0+C 2 sen 0 )+ 3 sen 0=1
Para C2 x ' ( 0 )=5 ; C2=0
Entonces la ecuación del movimiento es:
x ( t )=e−t cos 2t +3 sen 2 t
Donde
e−t cos 2 t : solucion transitoria
3 sen 2t : solucion estacionaria
b) grafica
Solución transitoria Solución Estacionaria
c)
x ( t )=e−t cos 2t +3 sen 2 t
Conclusiones
- Se identificó cada una de las ecuaciones diferenciales que modelan el
movimiento de un bloque unido a un resorte.
- Se reforzó lo estudiado, esto resolviendo las ecuaciones diferenciales que
explican el movimiento de un bloque unido a un resorte.
