Proyecto 2 de Ecuaciones Diferenciales
Proyecto 2 de Ecuaciones Diferenciales
Proyecto 2 de Ecuaciones Diferenciales
ECUACIONES DIFERENCIALES
ALIAS: HOMBRES ROCA
ALUMNOS :
AREQUIPA - PERÚ
2018
𝑑2 𝑦 𝑑𝑦
𝑎 +𝑏 + 𝑐𝑦 = 𝑔(𝑡)
𝑑𝑡 2 𝑑𝑡
Las matemáticas de, digamos un circuito eléctrico en serie son idénticas a las de un sistema
vibratorio masa/resorte. Las formas de esta ED de segundo orden se presentan en el análisis de
problemas en diversas áreas de la ciencia e ingeniería.
En esta parte consideraremos varios sistemas dinámicos lineales en los que cada modelo
matemático es una ecuación diferencial de segundo orden con coeficientes constantes junto
con condiciones iniciales especificas en un tiempo que tomaremos como 𝑡 = 0
𝑑2 𝑦 𝑑𝑦
𝑎 2
+𝑏 + 𝑐𝑦 = 𝑔(𝑡), 𝑦(0) = 𝑦0 , 𝑦 ′ (0) = 𝑦1.
𝑑𝑡 𝑑𝑡
Recordemos que la función g es la entrada, función de conducción o función forzada del sistema.
Una solución 𝑦(𝑡) de la ecuación diferencial en un intervalo 𝐼 que contiene a 𝑡 = 0 que satisface
las condiciones iniciales se llama salida o respuesta del sistema.
LEY DE HOOKE
Ley de Hooke Supongamos que, como en la figura (b), una masa m1 está unida a un resorte
flexible colgado de un soporte rígido. Cuando se reemplaza rn1 con una masa distinta m2, el
estiramiento, elongación o alargamiento del resorte cambiará.
Según la ley de Hooke, el resorte mismo ejerce una fuerza de restitución, F, opuesta a la
dirección del alargamiento y proporcional a la cantidad de alargamiento s. En concreto, F = Rs,
donde k es una constante de proporcionalidad llamada constante del resorte. Aunque las masas
con distintos pesos estiran un resorte en cantidades distintas, 6ste está caracterizado
esencialmente por su numero k; por ejemplo, si una masa que pesa 10 libras estira i pie un
resorte, entonces 10 = k(i) implica que k = 20 lb/ft. Entonces, necesariamente, una masa cuyo
peso sea de 8 libras estirará el resorte f de pie.
Después de unir una masa 𝑀 a un resorte, ésta lo estira una longitud 𝑠 y llega a una posición de
equilibrio, en la que su peso, 𝑊, está equilibrado por la fuerza de restauración 𝑘𝑠 . Recuérdese
que el peso se define por W = mg 𝑊 = 𝑚𝑔, donde la masa se expresa en slugs, kilogramos o
𝑓𝑡 𝑚
gramos y 𝑔 = 32 𝑠2 , 9.8 𝑠2 , respectivamente. Como se aprecia en la figura (b), la condición de
equilibrio es 𝑚𝑔 = 𝑘𝑠 ó 𝑚𝑔 − 𝑘𝑠 = 0. Si la masa se desplaza una distancia x respecto de su
posición de equilibrio, la fuerza de restitución del resorte es 𝑘(𝑘 + 𝑠). Suponiendo que no hay
fuerzas de retardo que actúen sobre el sistema y que la masa se mueve libre de otras fuerzas
externas (movimiento libre), entonces podemos igualar la segunda ley de Newton con la fuerza
neta, o resultante, de la fuerza de restitución y el peso:
𝑑2 𝑥
𝑚 = −𝑘(𝑠 + 𝑥) + 𝑚𝑔 = −𝑘𝑥 + 𝑚𝑔 − 𝑘𝑠 = −𝑘𝑥 (1)
𝑑𝑡 2
El signo negativo de la ecuación (1) indica que la fuerza de restitución del resorte actúa en la
dirección opuesta del movimiento. Además, podemos adoptar la convención que los
desplazamientos medidos abajo de la posición de equilibrio son positivos.
1. MOVIMIENTO NO AMORTIGUADO
Donde:
K1 y K2 : Constantes de deformación de los resortes.
M : Masa del bloque
Cuando no hay movimiento el sistema está en equilibrio por ende el peso es igual a al
normal (N) cual es la fuerza que se genera en contra del peso. Al empezar el movimiento
a una distancia (x) a la derecha del punto de equilibrio, existen dos fuerzas ambas en
contra del movimiento F1 que es la fuerza que el resorte de la izquierda al ser elongado
tiende a regresar a su estado original y la F2 cuál es la fuerza que se genera también en
contra del movimiento es la fuerza que el resorte de la derecha genera al ser
comprimido tiende a regresar a su estado original.
Lo mismo sucede si se empieza el movimiento a la izquierda.
∑ 𝐹 = 𝑚a
𝐹1 + 𝐹2 = 𝑚. 𝑎
−𝑘1 . 𝑥 − 𝑘2 . 𝑥 = 𝑚. 𝑎
𝑑2 𝑥
−(𝑘1 + 𝑘2 ). 𝑥 = 𝑚. 2
𝑑𝑡
Dividiendo ambos lados entre m, la ecuación diferencial seria del sistema con
movimiento no amortiguado:
𝑑2 𝑥 (𝑘1 + 𝑘2 )
+ 𝑥=0
𝑑𝑡 2 𝑚
De donde el cuadrado de la frecuencia angular es: (k1+k2)/m
𝑑2 𝑥
+ 𝑤 2𝑥 = 0
𝑑𝑡 2
Donde
Por la 2° ley de Newton que establece que la fuerza es directamente proporcional a la a masa y
a la aceleración de un cuerpo. En este caso también se aplicará en un sistema como si fuera
movimiento no amortiguado, pero con la diferencia que solo existe una fuerza de
amortiguamiento generado por el líquido viscoso. Ley de Hooke que enuncia cuando debido a
una fuerza externa, un resorte sufre una deformación respecto a su posición de equilibrio la
fuerza recuperadora que ejerce este muelle o resorte es directamente proporcional al cambio
de longitud
∑𝐹 = 𝑚 ∗ 𝑎
𝐹1 + 𝐹2 = 𝐹𝑎𝑚𝑜𝑟𝑡𝑖𝑔𝑢𝑎𝑐𝑖𝑜𝑛 = 𝑚 ∗ 𝑎
𝑑2
−𝑘1 (𝑥 ) − 𝐾2 (𝑥 ) − 𝐵 = 𝑚 2
𝑑𝑡
𝑑2
−(𝑘1 + 𝐾2 )(𝑥 ) − 𝐵 𝑑𝑥/𝑑𝑡 = 𝑚 2 𝑑𝑖𝑣𝑖𝑑𝑖𝑒𝑛𝑑𝑜 𝑒𝑡𝑟𝑒 𝑚
𝑑𝑡
Por lo tanto, la ecuación de frecuencia del movimiento amortiguado es
𝑑2 𝑥 𝐵 𝑑𝑥 (𝑘1 + 𝑘2 )
+ + 𝑥=0 (1)
𝑑𝑡 2 𝑚 𝑑𝑡 𝑚
3
La ecuación de la ecuación (1) es 𝑟 2 + 𝑚 𝑟 + 𝑤 2 = 0 cuyas raíces pendiente de formula general
√𝑏2 −4𝑎𝑐 𝐵 𝐵 𝐵
𝑟 = −𝑏 ± 2𝑎
son: 𝑟𝑘 = − 2𝑚 ± √(2𝑚)2 − 𝑤 2 con 2𝑚 > 0
𝐵 𝐵 𝐵 𝐵
𝑟1 = − 2𝑚 + √(2𝑚)2 − 𝑤 2 y 𝑟1 = − 2𝑚 − √(2𝑚)2 − 𝑤 2 .
𝑥(𝑡) = 𝑐1 𝑒 𝑟1 𝑡 + 𝑐2 𝑒 𝑟2 𝑡
2 𝐵 2
−
𝐵
𝑡 √( 𝐵 ) −𝑤 2 𝑡 −√( ) −𝑤 2
𝑥(𝑡) = 𝑒 2𝑚 (𝑐1 𝑒 2𝑚 + 𝑐2 𝑒 2𝑚
ASPECTOS RELEVANTES DE LA GRAFICA DEL MOVIMENTO SOBRE AMORTIGUADO
Ejemplo
- Una masa de un slug estira a un resorte en 1.28 pies una fuerza de amortiguamiento
numéricamente igual a 𝐵 veces la velocidad instantánea que actua sobre el contrapeso.
Deducir la ecuación del movimiento si la masa se suelta en 1 pie por debajo de la
posición de equilibrio con una velocidad hacia abajo en 1pie por segundo cuando 𝐵 =
12 Trace la gráfica en este caso.
𝑐1 + 𝑐2 = 1
𝑥(𝑡) = 𝑐1 𝑒 𝑟𝑡 + 𝑐2 𝑡𝑒 𝑟𝑡
La masa puede pasar por la posición de equilibrio a lo más una sola vez.
CASO N°3 MOVIMIENTO SUBAMORTIGUADO
𝐵 2
Ocurre cuando el discriminante ( ) − 𝑤 2 < 0,
2𝑚
𝐵 𝐵
𝑟1 = + √𝑤 2 − ( )2 𝑖
2𝑚 2𝑚
𝐵 𝐵
𝑟1 = − √𝑤 2 − ( )2 𝑖
2𝑚 2𝑚
𝐵
𝐵 = √𝑤 2 − (2𝑚)2 es la parte imaginaria
𝐵 𝐵 2 𝐵 2
𝑥(𝑡) = 𝑒 −2𝑚𝑡 (𝑐1 𝑐𝑜𝑠√𝑤 2 − ( ) 𝑡 + 𝑐2 𝑠𝑒𝑛√𝑤 2 − ( ) 𝑡)
2𝑚 2𝑚
∆ꙍ 1
=
ꙍ𝑜 𝑄
Esto expresa que el valor inverso de Q es igual al cociente entre la anchura ∆ꙍ y la frecuencia
ꙍ𝑜
Cuando en un sistema masa-resorte hay una fuerza externa 𝑓𝑒(𝑡) (fuerza de excitación). Se le
conoce como sistema con movimiento forzado. Cuya ecuación diferencial es la siguiente:
𝐹 = 𝑚𝑎
𝑑2 𝑥 𝑑𝑥
𝑚 = −𝐵 𝑘 𝑇 𝑥 + 𝑓𝑒 (𝑡) ; 𝑘 𝑇 = 𝑘1 + 𝑘2 , 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒𝑠 𝑑𝑒 𝑙𝑜𝑠 𝑑𝑜𝑠 𝑟𝑒𝑠𝑜𝑟𝑡𝑒𝑠
𝑑𝑡 2 𝑑𝑡
𝑑2 𝑥 𝑑𝑥
𝑚 2 +𝐵 𝑘 𝑇 𝑥 = 𝑓𝑒 (𝑡) ;
𝑑𝑡 𝑑𝑡
𝑓𝑒 (𝑡) = 𝑓𝑢𝑒𝑟𝑧𝑎 𝑒𝑥𝑡𝑒𝑟𝑛𝑎 𝑑𝑒 𝑙𝑎 𝑟𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖 𝑑𝑒𝑙 𝑙𝑖𝑞𝑢𝑖𝑑𝑜
𝑚 𝑑2 𝑥 𝐵 𝑑𝑥 𝑘𝑇 𝑓𝑡
+ + (𝑥 ) = ;𝐵 = 𝑓𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑎𝑚𝑜𝑟𝑡𝑖𝑔𝑢𝑎𝑚𝑖𝑒𝑛𝑡𝑜
𝑚 𝑑𝑡 2 𝑚 𝑑𝑡 𝑚 𝑚
𝑘
= 𝑊2
𝑚
𝑑2 𝑥 𝐵𝑑𝑥
2
+ + 𝑊𝑇2 (𝑥 ) = 𝐹(𝑡)
𝑑𝑡 𝑚𝑑𝑡
𝑑2 𝑥 𝐵𝑑𝑡
2
+ + 𝑊12 (𝑥 ) + 𝑊22 (𝑥 ) = 𝐹(𝑡)
𝑑𝑡 𝑚𝑑𝑡
Diferencias entre el Movimiento Armónico Simple y Movimiento Forzado
Ejemplo
En el sistema mostrado la energía cinética del sistema se almacena en la masa rígida “m”, la
energía potencial elástica en el resorte sin masa de constante “k”, y la capacidad de disipación
de energía en el amortiguador viscoso que se mueve con velocidad proporcional a la fuerza, con
constante de amortiguamiento “c”. El sistema queda totalmente definido mediante la
coordenada “X” (Figura 1.7). Para que el sistema sea lineal los parámetros “k”, “m”, y “c” deben
ser constantes y no depender de la variable “X”. Si se aplica una fuerza F(t) sobre la masa “m”,
en la dirección positiva de X tendríamos el sistema de la Figura 1.9, además podemos obtener
la ecuación (1.1) que representa el movimiento del sistema discreto básico, común a todos los
sistemas lineales de 1 GDL, puede establecerse aplicando el principio D’Alembert como se
muestra en la Figura 1.10, introduciendo la fuerza de inercia y estableciendo el equilibrio de
fuerzas en la dirección X:
Todos los sistemas lineales con 1 GDL conducen a la ecuación diferencial ordinaria de segundo
orden antes mostrada. Normalizando la ecuación (1.1) tenemos la ecuación (1.2).
Donde:
𝑐
𝛽=
2𝑚𝜔𝑛
𝐾
𝜔𝑛 = √
𝑚
𝑥𝑦 1
𝐻(𝜔)𝑀𝐴𝑋 = ( )𝑀𝐴𝑋 =
𝐹𝑜 2𝛽
𝛽= Factor de amortiguamiento
c = Coeficiente de amortiguamiento.
ω = Frecuencia de excitación.
𝑋(𝑡) = 𝑥𝐻 + 𝑥𝑃 (1.6)
Donde:
𝑥𝐻 = 𝐴ℯ −δt cos(ϖ𝐷 𝑡 − 𝜑)
𝐹0 1 − 𝑟2
𝑥𝑝 = . 𝑠𝑒𝑛(ϖt)
𝐾 (1 − 𝑟 2 )2 + (2𝑟𝛽)2
δ
𝜔𝐷 = 𝜔𝑛 √1 − ( )2
𝜔𝑛
𝐶 𝐶
𝛽= =
𝐶𝐶 2𝑚𝜔𝑛
𝜔
𝑟=
𝜔𝑛
XH = Solución homogénea.
XP = Solución particular.
𝛿= Ángulo de fase.
δ = Desplazamiento estático
La Figura 1.12 representa el Factor de Amplificación Dinámica H(ω) en función de “r” = ω/ωn,
para varios valores del amortiguamiento relativo 𝛽.
Figura 1.12 Gráfica del Factor De Amplificación Vs. Relación de Frecuencias [I-17].
Para un valor del amortiguamiento relativo que puede considerarse normal de 𝛽 = 0.2, la Figura
1.12 muestra, para frecuencias de excitación próximas a la frecuencia natural (r =1), la amplitud
resultante del desplazamiento puede ser hasta 5 veces el que seobtendría aplicando
estáticamente una fuerza de la misma amplitud. Sin embargo, para frecuencias de excitación
que excedan en más de un 50% la frecuencia natural, el desplazamiento dinámico es mucho
menor que el estático. De ahí la importancia de hacer un diseño dinámico adecuado y escoger
los parámetros k y m de modo que las posibles frecuencias de excitación esté lejos de la
frecuencia natural del sistema. Cuando la frecuencia de excitación coincide con la frecuencia
natural (r =1), se dice que se está en la “Condición de Resonancia”. Por otro lado, la Figura 1.13
representa el Desfase 𝜑 de la respuesta del sistema (la vibración) respecto a la excitación y
permite aprecia en la resonancia, el desfase es siempre 90º, independientemente del valor del
amortiguamiento relativo 𝛽 .
2𝛽
𝜑 = 𝑇𝑎𝑛−1 (1−𝑟2 ) (1.13)
Los valores máximos del factor de amplificación dinámica H(ω) se obtienen derivando respecto
a r e igualando a cero, obteniéndose que el máximo factor se produce para 𝑟 2 = 1 − 2β2
ligeramente inferior a 1 y su valor esta dado en la expresión (1.14):
xk 1
H(ω)MAX = ( )MAX = (1.14)
Fo 2β√1−β2