Laboratorio 3. Modelo Matemático de Un Sistema Masa-Resorte
Laboratorio 3. Modelo Matemático de Un Sistema Masa-Resorte
Laboratorio 3. Modelo Matemático de Un Sistema Masa-Resorte
Laboratorio 3
MODELO MATEMATICO DE UN SISTEMA MASA RESORTE
Docente:
Teófilo Gómez
Integrantes:
Ariel Beitia 4-786-2344
Eynar Cubilla 4-844-492
Francisco Delgado 4-778-1994
I Semestre, 2020
Objetivos Generales
1.1. Desarrollar y analizar el modelo matemático de un sistema masa-resorte, de
un grado de libertad, bajo vibración libre, sin amortiguamiento, con
movimiento de traslación puro.
Objetivos Específicos
1.2. Obtener la ecuación de equilibrio estático de un sistema masa-resorte.
1.3. Determinar la deformación del resorte en equilibrio estático.
1.4. Determinar la ecuación diferencial de movimiento de un sistema masa-
resorte en función de la variable 𝑥. 𝑥 es la posición de la masa con respecto
a la posición de equilibrio estático. Identifique el tipo de ecuación diferencial
que resulta. Defina las condiciones iniciales de movimiento.
1.5. Determinar la ecuación diferencial de movimiento de un sistema masa-
resorte en función de la variable 𝑦. 𝑦 es la posición de la masa con respecto
a la posición del resorte no deformado. Identifique el tipo de ecuación
diferencial que resulta. Defina las condiciones iniciales de movimiento.
1.6. Encontrar la solución de la ecuación diferencial de movimiento desarrollada
en el punto 2.3 para 𝑥(0) = 𝑋0 y 𝑥̇ (0) = 0. Obtenga expresiones para la
posición 𝑥(𝑡), la velocidad 𝑥̇ (𝑡) y la aceleración 𝑥̈ (𝑡). Grafique los resultados,
utilice EXCELL. Para dos ciclos de movimiento.
1.7. Calcule la frecuencia natural, el periodo del movimiento y la frecuencia
natural angular de oscilación.
1.8. Encontrar la solución de la ecuación diferencial de movimiento desarrollada
en el punto 2.3 para 𝑥(0) = −𝛿 y 𝑥̇ (0) = 0. Obtenga expresiones para la
posición 𝑥(𝑡), la velocidad 𝑥̇ (𝑡) y la aceleración 𝑥̈ (𝑡). Grafique los resultados,
utilice EXCELL. Para dos ciclos de movimiento.
1.9. Obtenga la solución de la ecuación diferencial, punto 2.5, mediante MatLab
ó SciLab. Grafique la posición, velocidad y aceleración. Para dos ciclos de
movimiento.
1.10. Desarrolle el diagrama de bloques correspondiente, punto 2.5, y
obtenga la solución mediante SIMULINK ó XSICO. Grafique la posición,
velocidad y aceleración. Para dos ciclos de movimiento.
1.11. Analice la posición, velocidad y aceleración de la masa m. ¿Qué puede
concluir respecto a la amplitud y ángulo de fase de cada movimiento?
Contenido
Un sistema masa-resorte vibrará libremente al desplazarse de su posición de equilibrio
estático y liberarse. El sistema es conservativo, no está sujeto a fuerzas no-conservativas
ni a excitaciones externas. La ecuación gobernante del movimiento oscilatorio es una
ecuación diferencial de segundo grado, homogénea con coeficientes constantes. La
solución de dicha ecuación corresponde a la solución complementaria en donde las
constantes dependen de las condiciones iniciales del sistema.
𝑑𝐸
𝑑𝑡
= 0, 𝑑𝑒𝑏𝑖𝑑𝑜 𝑎 𝑞𝑢𝑒 𝐸 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 (2.3)
𝑚𝑥̈ + 𝑘𝑥 = 0 (2.5)
∑ 𝐹 = 𝑚𝑥̈ (2.6)
−𝑘𝛿 + 𝑚𝑔 = 0 (2.9)
𝑚𝑥̈ + 𝑘𝑥 = 0 (2.10)
𝑥 = 𝐶𝑒 𝑠𝑡 (2.11)
𝑥̈ = 𝐶𝑒 𝑠𝑡 𝑠 2 (2.12)
(𝑚𝑠 2 + 𝑘) = 0 (2.14)
𝜔𝑛 = √𝑘⁄𝑚 (2.16)
Remplazando
𝐶1 = 𝑋 cos 𝜙 y 𝐶2 = 𝑋 sin 𝜙 (2.18)
Resulta
𝑥(𝑡) = 𝑋 sin 𝜔𝑛 𝑡 cos 𝜙 + 𝑋 cos 𝜔𝑛 𝑡 sin 𝜙 (2.19)
𝐶
𝜙 = tan−1 ( 2 ) (2.22)
𝐶1
Metodología y Procedimiento
1 Escoja el resorte con la menor constante y la masa correspondiente.
Obtener el modelo matemático del sistema masa-resorte de un grado de
libertad construido. Para punto 2.3 y punto 2.5.
Determinar esta constante
previamente, para la masa de 100
gr. Si ya la tiene del lab. anterior
puede usarla. Debe ser la constante
con respecto a la posición del
bloque color aqua mostrado en la
imagen.
Figura 2. Constante para el resorte 1.
2 Resolver la ecuación diferencial de movimiento para diferentes condiciones
iniciales.
3 Obtener analíticamente y experimentalmente la frecuencia angular natural,
la frecuencia natural y el periodo natural de movimiento.
Nota: los datos experimentales solo son para compararlos con los
analíticos.
4 Graficar la posición, la velocidad y la aceleración (datos analíticos paso 3
anterior). Para punto 2.3 y punto 2.5 (6 graficas).
5 Escoja un resorte con una constante mayor al anterior y la masa
correspondiente. Obtener el modelo matemático del sistema masa-resorte
de un grado de libertad construido.
Determinar esta constante
previamente, para la masa de 100
gr. Si ya la tiene del lab anterior
puede usarla. Debe ser la Cte con
respecto a la posición del bloque
color aqua mostrado en la imagen.
Figura 3. Constante para el resorte 2.
Procedimiento
a. Selecciones los parámetros (masa y constante de resorte).
b. Especifique las condiciones iniciales indicadas en los puntos 2.5 y 2.6 de los
objetivos específicos.
c. Obtenga la ecuación diferencial de movimiento en función de 𝑥. Obtener la
posición, velocidad y aceleración para:
• 𝑥(0) = 𝑋0 y 𝑥̇ (0) = 0. Graficar utilizando: EXCELL, MatLab ó SciLab y
SIMULINK ó XSICO.
• 𝑥(0) = −𝛿 y 𝑥̇ (0) = 0. Graficar utilizando EXCELL.
d. Obtenga la ecuación diferencial de movimiento en función de 𝑦. Obtener la
posición, velocidad y aceleración para:
• 𝑦(0) = 𝑋0 + 𝛿 y 𝑦̇ (0) = 0. Graficar utilizando: EXCELL, MatLab ó SciLab
y SIMULINK ó XSICO. 𝑦(0) = 0 y 𝑦̇ (0) = 0. Graficar utilizando EXCELL.
Recursos
Resorte de tensión
Discos
Marco para soporte
Base para los discos
Balanza
Cinta métrica
Resultados
• Determinación del modelo matemático para el sistema masa-resorte 1:
𝑚𝑔
𝐹 = 𝑘𝑥 = 𝑚𝑔 → 𝑘 =
𝑥
1 kg m
(100 g) ( )
1000 g (9.81 s2 ) N
𝑘= −2
= 6.13
16 × 10 m m
El modelo matemático de este sistema viene dado por la ecuación:
mẍ + kx = 0
0.1ẍ + 6.13x = 0
ẍ + 61.3x = 0
Esta es una ecuación diferencial de segundo grado, homogénea con coeficientes
constantes. Debido a la ausencia de un elemento amortiguador, la solución de esta
ecuación diferencial será armónica:
ẋ (0)
Dónde: A = x(0) y B = wn
x(t) = x0 ∙ cos wn ∙ t
Variable Amplitud
x(t) X (0)
ẋ (t) x(0) ∙ wn
ẍ (t) x(0) ∙ wn2
Medición (N°) 1 2 3
Tiempo (s) 2.38 2.44 2.45
τpromedio 2.42
τ= = = 0.8067 s
3 3
1 1
fn = = = 1.24 Hz
τn 0.8067
wn 7.83
fn = = = 1.24 Hz
2π 2π
1 1
τn = = = 0.8064s
fn 1.24
x(t) = x0 cos(7.83 ∙ 𝑡)
Así: x(0) = 10cm = 0.1m. Finalmente, las ecuaciones diferenciales para las
condiciones iniciales x(0) = X0 y ẋ (0) = 0:
x(t) = 0.1cos(7.83 ∙ 𝑡)
ẋ (t) = −0.783sencos(7.83 ∙ 𝑡)
ẍ (t) = −6.130 cos cos(7.83 ∙ 𝑡)
0.1ẍ + 10.9x = 0
ẍ + 109x = 0
Esta es una ecuación diferencial de segundo grado, homogénea con coeficientes
constantes. Debido a la ausencia de un elemento amortiguador, la solución de esta
ecuación diferencial será armónica.
ẋ (0)
Dónde: A = x(0) y B = , A y B se obtienen a partir de las condiciones iniciales
wn
del sistema: 𝑥(0) = 𝑋0 y 𝑥̇ (0) = 0. Así, reduciendo la solución y derivando:
x(t) = x0 ∙ cos wn ∙ t
Medición (N°) 1 2 3
Tiempo (s) 1.81 1.83 1.81
τpromedio 1.82
τ= = = 0.6067 s
3 3
1 1
fn = = = 1.65 Hz
τn 0.6067
wn 10.44
fn = = = 1.66 Hz
2π 2π
1 1
τn = = = 0.6025 s
fn 1.66
x(t) = x0 cos(10.44 ∙ 𝑡)
ẋ (t) = −10.44x0 sencos(10.44 ∙ 𝑡)
ẍ (t) = −109x0 cos cos(10.44 ∙ 𝑡)
x(t) = 0.1cos(10.44 ∙ 𝑡)
ẋ (t) = −1.044encos(10.44 ∙ 𝑡)
ẍ (t) = −10.90 cos cos(10.44 ∙ 𝑡)
Graficas
• ẍ + 61.3x = 0
y=x+α
m𝑦̈ + 𝑘𝑥 − 𝑚𝑔 = 0
k
ωn = √
m
rad
ωn = 2πf ( )
seg
1
τ= (seg)
f
Todos estos parámetros dependen de solo dos variables de un sistema, su masa y
su constante de rigidez equivalente.
Evaluación sugerida
35% Coherencia (tecnisismo)
35% Tecnica de desarrollo (destreza)
15% orden y estética
15% ortografía y puntualidad
100%
Referencias
Vibraciones Mecánicas. Singiresu S. Rao. Quinta edición. PEARSON EDUCATION,
México, 2012.