Universidad Tecnológica de Panamá

Facultad de Ingeniería Eléctrica

Ingeniería Electromecánica
Dinámica Aplicada
Anthony Vigil

VIBRACIÓN LIBRE DE UN SISTEMA MASA-RESORTE.

Nombre: Héctor Murgas

Cédula: 4-805-135
e-mail: [email protected]
Nombre: Jesús Castillo
Cédula: 4-808-2303
e-mail: [email protected]
Nombre: Angel Justavino
Cédula: 8-954-2425
e-mail: [email protected]
Nombre: Angélica Calderón
Cédula: E-8-184053
e-mail: [email protected]

Resumen En esta experiencia de laboratorio se analiza el modelo físico y matemático de un oscilador armónico simple, para ello
se realizó el modelado matemático de forma manual y posteriormente se modeló el sistema mediante el software Scilab,
encontrando el gráfico de su comportamiento oscilatorio, este procedimiento se realizó tres veces, uno por cada constante de
elasticidad del resorte, proporcionada por el instructor de laboratorio; una vez completado el modelado matemático y encontrada
la ecuación de desplazamiento del sistema, se procede a encontrar la frecuencia natural circular, el periodo y la frecuencia
natural del sistema, debido a cada valor de K proporcionado, finalmente, se encuentran las constantes de la ecuación solución
de la ecuación diferencial del sistema, haciendo uso de las condiciones inciales proporcionadas por el instructor, esto es esencial
para encontrar los tiempos que ocasionan valores máximos en los gráficos de desplazamiento, velocidad y aceleración, una vez
determinados estos valores de tiempo, se encuentran los máximos de cada una de estas funciones.

Descriptores: frecuencia natural, grado de libertad, modelo matemático, periodo, vibración libre

1. Introducción.

El movimiento de un sistema, o la dinámica de un sistema, gracias a las leyes de la mecánica que pudo expresar Newton y las
matemáticas, se puede describir con un modelo y con este modelo, en algunos casos, podemos determinar los estados del sistema
en cualquier momento en el tiempo, esto se conoce como sistema determinista. En este caso el sistema masa resorte, que es también
un modelo de un oscilador armónico simple. El modelo matemático del oscilador armónico simple, lo describe una ecuación
diferencial, que es la segunda ley del movimiento de Newton, de la cual para su solución necesitaremos condiciones iniciales las
cuales asumiremos como 0 o nulas, luego de resolver la ecuación diferencial que describe nuestro sistema, obtendremos la posición
en función del tiempo en forma de onda senoidal, lo que responde a la naturaleza oscilatoria de las vibraciones o del sistema masa
resorte.

2. Materiales y métodos.

• SciLab: Se utilizó este software pare encontrar el gráfico oscilatorio del sistema, para ello se utilizaron distintas
constantes de elasticidad, proporcionadas por el instructor de laboratorio, así como los valores correspondientes a las
condiciones iniciales del sistema.
• Los métodos matemáticos utilizados en este documento, son la segunda ley de newton, para encontrar la ecuación
diferencial para la cual, en un oscilador armónico simple, la ley de conservación de la energía, conocimientos de
resolución de ecuaciones diferenciales de segundo grado, homogéneas, con coeficientes constantes, las ecuaciones de
frecuencia natural, periodo y frecuencia natural del sistema, para cada valor de constante de elasticidad proporcionada,
y conocimientos en resolución de sistemas de ecuaciones.
• Guia proporcionada por el instructor de laboratorio.

3. Resultados y discusión

Conociendo el sistema masa-resorte, estudiado en los laboratorios anteriores, se procede modelar matemáticamente el sistema:
Condiciones Iniciales:
𝑥(0) = 0.02 𝑚 ; 𝑥(0) = 0
Haciendo uso de la tercera ley de Newton, presentada en la ecuación (1) es posible modelar matemáticamente el
comportamiento del sistema, de la siguiente forma:

𝑑
∑𝐹 = (𝑚𝑥 ) = 𝑚𝑥̈ (1)
𝑑𝑡
Para el caso específico de un sistema masa resorte, la ecuación resultante sería la mostrada en la ecuación (2):

𝑑
∑ 𝐹 = (𝑚𝑥 ) = 𝑚𝑥̈
𝑑𝑡
𝑚𝑔 − 𝑘(𝑥 + 𝛿 ) = 𝑚𝑥̈ (2)

Debido a que en el equilibrio estático 𝑚𝑥̈ = 0, obtenemos la ecuación (3), y al aplicar conservación de la energía,
obtenemos la ecuación (4), finalmente, al derivarla encontramos la ecuación (5) que rige el sistema, la cual se utilizará para realizar
el modelado mediante scilab.
𝑚𝑔 − 𝑘(𝑥 + 𝛿 ) = 0 (3)
Aplicando conservación de la energía:
𝑑
(𝐸 ) = 0
𝑑𝑡 𝑡
𝑑 1 1
( 𝑘𝑥 2 + 𝑚𝑥 2 ) = 0 (4)
𝑑𝑡 2 2
𝑚𝑥̈ + 𝑘𝑥 = 0 (5)
La solución de esta ecuación diferencial homogénea, de segundo grado, con coeficientes constantes, es la mostrada en la
ecuación (6):
𝑥(𝑡) = 𝐶1 sin(𝜔𝑛 𝑡) + 𝐶2 cos(𝜔𝑛 𝑡) (6)

A continuación, se presenta la solución del sistema masa-resorte, mediante el software Scilab, utilizando la ecuación (5),
y conociendo que las condiciones iniciales del sistema son:

𝑥(0) = 0.02 𝑚 ; 𝑥(0) = 0

Con las constantes del resorte proporcionadas pos el instructor y mostradas en la tabla 1, se realiza el diagrama de bloques
y se encuentra el gráfico de respuesta del sistema, a continuación se presentan en la figura 1 y 2, el diagrama de bloques y el gráfico
para k1, en la figura 3 y 4, se presenta el diagrama de bloques y el gráfico para k2, Finalmente en las figuras 5 y 6, se presenta el
diagrama de bloques y el gráfico para k3.

Tabla 1. Valores de constantes de elasticidad del resorte proporcionados por el instructor.

Constante Valor
K1 2.9744 N/m
K2 7.9764 N/m
K3 12.076 N/m

• Resorte #1 con K= 2.9744 N/m

Fig 1. Diagrama de resorte 1 en Xcos.

 Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Dinámica Aplicada
Anthony Vigil

Fig 2. Grafica de resorte 1 en Xcos.

• Resorte #2 K= 7.9764 N/m

Fig 3. Diagrama de resorte 2 en Xcos.

Fig 4. Grafica de resorte 2 en Xcos.

 • Resorte #3 con K=12.076N/m

Fig 5. Diagrama de resorte 3 en Xcos.

Fig 6. Grafica de resorte 3 en Xcos.

Del laboratorio anterior, conocemos las ecuaciones para encontrar la frecuencia natural circular del sistema, presentada
en la ecuación (7), el periodo del sistema que corresponde a le ecuación (8), y finalmente la frecuencia natural, mostrada en la
ecuación (9), con estas ecuaciones, podremos encontrar estas incógnitas para cada constante de elasticidad y llenar la tabla 2.
𝑘
𝝎𝒏 = √𝑚 (7)
𝑚
𝝉𝒏 = 2𝜋√ 𝑘 (8)
𝟏
𝑓𝑛 = 𝝉 (9)
𝒏

A continuación, se presentan los cálculos correspondientes para la frecuencia natural circular del sistema, debido
a diferentes valores de K:

Resorte 1

2.9744 𝑁/𝑚
𝝎𝒏 = √ = 4.6093 𝑟𝑎𝑑/𝑠
0.14 𝑘𝑔
 Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Dinámica Aplicada
Anthony Vigil

Resorte 2

7.9764 𝑁/𝑚
𝝎𝒏 = √ = 7.5481 𝑟𝑎𝑑/𝑠
0.14𝑘𝑔

Resorte 3

12.076 𝑁/𝑚
𝝎𝒏, = √ = 9.2875 𝑟𝑎𝑑/𝑠
0.14𝑘𝑔

A continuación, se presentan los cálculos correspondientes para el periodo del sistema, debido a diferentes valores de k:

Resorte 1

0.14 𝑘𝑔
𝝉𝒏 = 2𝜋√ = 1.3632𝑠
2.9744 𝑁/𝑚

Resorte 2

0.14 𝑘𝑔
𝝉𝒏 = 2𝜋√7.9764 𝑁/𝑚 = 0.83242𝑠

Resorte 3

0.14 𝑘𝑔
𝝉𝒏 = 2𝜋√12.076 𝑁/𝑚 = 0.67652𝑠

continuación, se presentan los cálculos correspondientes para la frecuencia natural del sistema, debido a diferentes valores
de K:

Resorte 1
1
𝒇𝒏 = 1.3632𝑠 = 0.733595ℎ𝑧

Resorte 2
1
𝒇𝒏 = 0.83242𝑠 = 1.20132ℎ𝑧

Resorte 3
1
𝒇𝒏, = = 1.47815 ℎ𝑧
0.67652𝑠

Para encontrar las constantes, y tener los valores máximos y mínimos, derivamos la ecuación (6) y encontramos la
ecuación (7), evaluamos las ecuaciones (6) y (10) en sus condiciones iniciales y encontramos los valores de las constantes:
𝑥(𝑡) = 𝐶1 sin(𝜔𝑛 𝑡) + 𝐶2 cos(𝜔𝑛 𝑡) (6)
𝑥 (𝑡) = 𝐶1 𝜔𝑛 cos(𝜔𝑛 𝑡) − 𝐶2 𝜔𝑛 sin(𝜔𝑛 𝑡) (10)

Evaluando en 𝑥(0) = 0.02 𝑚 ; 𝑥(0) = 0, y solucionando el sistema de ecuaciones de dos incógnitas encontramos los
valores de C1 y C2
C1 = 0
C2 = 0.02

Por lo tanto, la función de desplazamiento de nuestro sistema, es la mostrada en la ecuación (11):

𝑥(𝑡) = 0.02 cos(𝜔𝑛 𝑡) (11)
 Esta función derivada, se muestra en la ecuación (12):
𝑥 (𝑡) = −0.02 𝜔𝑛 sin(𝜔𝑛 𝑡) (12)

La función de aceleración del sistema, se muestra en la ecuación (13):

𝑥̈ (𝑡) = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡) (13)

Finalmente, para encontrar los valores de 𝑥𝑚𝑎𝑥 , 𝑥𝑚𝑎𝑥 ,𝑥̈ 𝑚𝑎𝑥 , es necesario derivar la ecuación (13), la ecuación resultante
de esta derivación se muestra en la ecuación (14); posteriormente, se deben evaluar las ecuaciones (12), (13) y (14) en cero, para
encontrar los tiempos máximos para cada función, y posteriormente, evaluar la función original en estos tiempos, para cada valor
de frecuencia natural circular correspondiente para los distintos valores de “k”.

𝑥⃛(𝑡) = 0.02 𝜔𝑛 3 sin(𝜔𝑛 𝑡) (14)

Evaluando las ecuaciones (12), (13) y (14) en cero, obtenemos las ecuaciones (15), (16), (17), las cuáles nos permiten
obtener los valores de tiempo para los cuales las funciones son máximas, para cada frecuencia natural circular del sistema:
0 = −0.02 𝜔𝑛 sin(𝜔𝑛 𝑡1 ) (15)
0 = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡2 ) (16)
0 = 0.02 𝜔𝑛 3 sin(𝜔𝑛 𝑡3 ) (17)

Reemplazando los valores de 𝜔𝑛 , para cada valor de K, podemos obtener los tiempos para los cuales las funciones son
máximas para cada función, donde:
𝑡1 : representa el tiempo para el cual 𝑥(𝑡) es máximo.
𝑡2 : representa el tiempo para el cual 𝑥 (𝑡) es máximo.
𝑡3 : representa el tiempo para el cual 𝑥̈ (𝑡) es máximo.

Para la frecuencia natural circular de 4.60931 rad/s, correspondiente a k=2.9744N/m:

Para t1:
0 = −0.02 𝜔𝑛 sin(𝜔𝑛 𝑡1 )
0 = −0.02 (4.60931) sin(4.60931𝑡1 )
𝑡1 = 0s

Para t2:
0 = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡2 )
0 = −0.02 (4.60931)2 cos(4.60931𝑡2 )
𝑡2 = 0.34078 𝑠

Para t3:
0 = 0.02 𝜔𝑛 3 sin(𝜔𝑛 𝑡3 )
0 = 0.02 (4.60931)3 sin(4.60931𝑡1 )
𝑡3 = 0𝑠

Para la frecuencia natural circular de 7.5481 rad/s, correspondiente a k=7.9764N/m:

Para t1:
0 = −0.02 𝜔𝑛 sin(𝜔𝑛 𝑡1 )
0 = −0.02 (7.5481) sin(7.5481 𝑡1 )
𝑡1 = 0 s

Para t2:
0 = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡2 )
0 = −0.02 (7.5481)2 cos(7.5481𝑡2 )
𝑡2 = 0.208105 s

Para t3:
0 = 0.02 𝜔𝑛 3 sin(𝜔𝑛 𝑡3 )
0 = 0.02 (7.5481)3 sin(7.5481𝑡1 )
𝑡3 = 0s
 Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Dinámica Aplicada
Anthony Vigil

Para la frecuencia natural circular de 9.28747 rad/s, correspondiente a k=12.076 N/m:

Para t1:
0 = −0.02 𝜔𝑛 sin(𝜔𝑛 𝑡1 )
0 = −0.02 (9.28747 ) sin(9.28747 𝑡1 )
𝑡1 = 0 s

Para t2:
0 = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡2 )
0 = −0.02 (9.28747 )2 cos(9.28747 𝑡2 )
𝑡2 = 0.169131 s

Para t3:
0 = 0.02 𝜔𝑛 3 sin(𝜔𝑛 𝑡3 )
0 = 0.02 (9.28747 )3 sin(9.28747 𝑡1 )
𝑡3 = 0 𝑠

Una vez encontrados los valores de tiempo, para los cuales el desplazamiento, velocidad y aceleración son máximos; se
procede a evaluar las ecuaciones (11), (12) y (13), para los valores de tiempo t1, t2, t3 , respectivamente, y se reemplazan los valores
para cada frecuencia natural circular correspondiente a cada valor de “k”.

Para la frecuencia natural circular de 4.60931 rad/s, correspondiente a k=2.9744N/m:

Para t1:
𝑥(𝑡) = 0.02 cos(𝜔𝑛 𝑡1 )
𝑥(0.651561) = 0.02 cos[(4.60931)(0)]
𝑥(0.651561) = 0.02𝑚

Para t2:
𝑥 (𝑡) = −0.02 𝜔𝑛 sin(𝜔𝑛 𝑡1 )
𝑥 (0.34078) = −0.02 (4.60931) sin[(4.60931)(0.34078)]
𝑥 (0.34078) = −0.0921862 m/s

Para t3:
𝑥̈ (𝑡) = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡2 )
𝑥̈ (0.681561) = −0.02 (4.60931)2 cos[(4.60931)(0)]
𝑥̈ (0.681561) = −0.4249 𝑚⁄𝑠 2

Para la frecuencia natural circular de 7.5481 rad/s, correspondiente a k=7.9764N/m:

Para t1:
𝑥(𝑡) = 0.02 cos(𝜔𝑛 𝑡1 )
𝑥(0.416208) = 0.02 cos[(7.5481)(0)]
𝑥(0.416208) = 0.02 𝑚

Para t2:
𝑥 (𝑡) = −0.02 𝜔𝑛 sin(𝜔𝑛 𝑡1 )
𝑥 (0.416208) = − 0.02 (7.5481) sin[(7.5481)(0.2081)]
𝑥 (0.416208) = −0.150963m/s

Para t3
𝑥̈ (𝑡) = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡2 )
𝑥̈ (0.416208) = − 0.02 (7.5481)3 sin[(7.5481)(0)]
𝑥̈ (0.416208) = −1.13948 𝑚⁄𝑠 2
 Para la frecuencia natural circular de 9.28747 rad/s, correspondiente a k=12.076 N/m:
Para t1:
𝑥(𝑡) = 0.02 cos(𝜔𝑛 𝑡1 )
𝑥(0.33824) = 0.02 cos[(9.28747)(0)]
𝑥(0.33824) = 0.02 𝑚

Para t2:
𝑥 (𝑡) = −0.02𝜔𝑛 sin(𝜔𝑛 𝑡2 )
𝑥 (0.33826) = −0.02 (9.28747) sin[(9.28747)(0.1691)]
𝑥 (0.33826) = −0.185749𝑚/𝑠

Para t3:
𝑥̈ (𝑡) = −0.02 𝜔𝑛 2 cos(𝜔𝑛 𝑡2 )
𝑥̈ (0.33826) = −0.02(9.28747)2 cos[(9.28747)(0)]
𝑥̈ (0.33826) = −1.72514198 𝑚⁄𝑠 2

Una vez determinados los valores máximos para el desplazamiento, velocidad y aceperación, se completa la tabla 1.

Tabla 2. Respuesta de un sistema masa-resorte

Respuesta de un sistema masa-resorte
𝒌𝒊 𝒎𝒊 𝝎𝒏 𝒇𝒏 𝝉𝒏 𝒙𝒎𝒂𝒙 𝒙𝒎𝒂𝒙 𝒙̈ 𝒎𝒂𝒙
2.9744 N/m 0.14kg 4.60931 rad/s 0.73359 hz 1.36315s 0.02 m -0.0921862 m/s -0.4249m/s2
7.9764 N/m 0.14kg 7.54813 rad/s 1.20132 hz 0.83242𝑠 0.02m -0.150963m/s -1.13948 m/s2
12.076 N/m 0.14kg 9.28747 rad/s 1.47815 hz 0.67652𝑠 0.02m -0.185749m/s -1.72514198m/s2

Preguntas.
1. Que concluye respecto a las frecuencias angulares, naturales y periodos naturales de oscilacion, para los sistemas Masa-
Resorte estudiados?
De manera general podemos concluir que según mayor sea la constante del resorte mayores valores tendremos para la
frecuencia naturales y frecuencias angulares y menores valores del periodo.
Para distintas constantes del resorte “k” y masas iguales en cada sistema las frecuencias angulares se elevean dando
resultados de mayor velocidad para la masa del sistema, los periodos de manera logica y experimental se logro observar
que transcurria menos tiempo para el recorrido de una oscilacion (para la misma dstancia y mayor velocidad, menor
tiempo) y finalmente las frecuencias naturalres del sistema el aumento de la constante del resorte provoca una mayor
ratio de frecuencia donde la onda se repetia mas veces.

2. Como compara las amplitudes de posicion, velocidad y aceleracion para los tres sistemas estudiados?

Las posiciones en los tres sistemas estudiados son iguales, obvio ya que restringimos nuestro sistema para una
elongacion en concreto.
Las velocidades son similares entre ellas para la primera constante del resorte obtuvimos la menor velocidad y es la que
es menos similar en comparacion para el segundo y tercer caso de la constante del resorte que tienen velocidades muy
parecidas.
En cuanto a la aceleracion contamos con dos casos particulares para el primer caso y el segundo contamos con el valor
de la aceleracion negativa, con esto podemos decir que la velocidad va decayendo en el tiempo, y para el tercer caso
tenemos lo contrario el cual la aceleracion aumenta.

3. Cual de los tres metodos de solucion considera mas conveniente? Explique.

El metodo mas eficiente y acertado, con menos margen de error, es simulink ya que con el manejo de estas herramientas
los problemas de sistema de masa – resorte se pueden simular de manera facil, rapida, y sencilla sin la necesidad de caer
caluclos manuales que podrian tener un nivel de error ya que siempre que conamos con la mano del hombre existe la
posibilidad ded un margen de error y cuando es un software hay mucho menos posibilidades.
 Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Dinámica Aplicada
Anthony Vigil

4. Determinar la masa, la constante de elasticidad, la frecuencia natural y las condiciones iniciales del problema indicado
por el instructor en el punto e.

Las condiciones iniciales son,

𝑥(0) =0.2 𝑥 (0) = −0.03 𝑥̈ (0) = −0.35

Para un sistema masa resorte tenemos la siguente ecuacion diferencial,

𝑚𝑥̈ + 𝑘𝑥 = 0

Obtenemos la siguente funcion solucion para esta ecuacion diferencial,

𝑥(𝑡) = 𝐶1 sin(𝜔𝑛 𝑡) + 𝐶2 cos(𝜔𝑛 𝑡)

Derivando dicha funcion,

𝑥 (𝑡) = 𝐶1 𝜔𝑛 cos(𝜔𝑛 𝑡) − 𝐶2 𝜔𝑛 sin(𝜔𝑛 𝑡)

Ahora evualuando condiciones inicales para encontrar 𝐶1 y 𝐶2

0.2 = 𝐶1 sin(0) + 𝐶2 cos(0)

−0.03 = 𝐶1 𝜔𝑛 cos(0) − 𝐶2 𝜔𝑛 sin(0)

Obtenemos que
−0.03
𝐶1 = 𝜔𝑛
𝐶2 =0.2
−0.03
𝑥 (𝑡 ) = sin(𝜔𝑛 𝑡) + 0.2 cos(𝜔𝑛 𝑡)
𝜔𝑛

De la gráfica, podemos observar que el periodo para el cuál se cumple una oscilación del sistema, es:

𝝉𝒏 = 𝟓𝒔

Por lo tanto, la frecuencia natural del sistema, a partir de la ecuación (9), es:
𝟏 𝟏
𝑓𝑛 = = = 𝟐𝟎𝑯𝒛
𝝉𝒏 𝟓𝒔
A partir de la ecuación (8), podemos encontrar el valor de k/m:

𝝉𝒏 𝑚
= √
2𝜋 𝑘

𝟐
5 𝑚
( ) =
2𝜋 𝑘
𝑚
0.63326 =
𝑘
Por lo tanto, la masa es:
𝑚 = 0.63326𝑘

Y la constante es:
𝑘 = 1.5791𝑚
5. Encontrar la solucion de la ecuacion diferencial de movimiento para 𝑥(0) = −𝛿 y 𝑥 (0) = 0. Obtenga las expresiones
para la posicion x(t), la velocidad 𝑥 (𝑡) y la aceleracion 𝑥̈ (𝑡). Grafique los resultados, utilice MatLab/Scilab y
Simulink/Xcos para dos ciclos de movimiento.
Para un sistema masa resorte donde la ecuacion diferencial es,

𝑚𝑥̈ + 𝑘𝑥 = 0

Condiciones inciales,

𝑥(0) = −𝛿 y 𝑥 (0) = 0

Obtenemos la siguente funcion solucion para esta ecuacion diferencial,

𝑥(𝑡) = 𝐶1 sin(𝜔𝑛 𝑡) + 𝐶2 cos(𝜔𝑛 𝑡)

Derivando dicha funcion,

𝑥 (𝑡) = 𝐶1 𝜔𝑛 cos(𝜔𝑛 𝑡) − 𝐶2 𝜔𝑛 sin(𝜔𝑛 𝑡)

Ahora evualuando condiciones inicales para encontrar 𝐶1 y 𝐶2

−𝛿 = 𝐶1 sin(0) + 𝐶2 cos(0)
0 = 𝐶1 𝜔𝑛 cos(0) − 𝐶2 𝜔𝑛 sin(0)

Obtenemos que

𝐶1 = 0
𝐶2 = −𝛿

Remplazando 𝐶2 en la función,

𝑥(𝑡) = −𝛿 cos(𝜔𝑛 𝑡)

Resorte 1

𝑥(𝑡) = −𝛿 cos(4.60931 𝑡)

Grafico 1. Posicion vs Tiempo para Resorte 1

 Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Dinámica Aplicada
Anthony Vigil

𝑥 (𝑡) = 4.60931𝛿 sin(4.60931 𝑡)

Grafico 2. Velocidad vs Tiempo para Resorte 1

𝑥̈ (𝑡) = 21.2485𝛿 cos(4.60931 𝑡)

Grafico 3. Aceleración vs Tiempo para Resorte 1

Resorte 2

𝑥(𝑡) = −𝛿 cos(7.54813 𝑡)

Grafico 4. Posicion vs Tiempo para Resorte 2

𝑥 (𝑡) = 7.54813𝛿 sin(7.54813 𝑡)

Grafico 5. Velocidad vs Tiempo para Resorte 2

𝑥̈ (𝑡) = 56.9743𝛿 cos(4.60931 𝑡)

Grafico 6. Aceleración vs Tiempo para Resorte 2

Resorte 3

𝑥(𝑡) = −𝛿 cos(9.28747 𝑡)

Grafico 7. Posicion vs Tiempo para Resorte 3

 Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Dinámica Aplicada
Anthony Vigil

𝑥 (𝑡) = 9.28747𝛿 sin(9.28747 𝑡)

Grafico 8. Velocidad vs Tiempo para Resorte 3

𝑥̈ (𝑡) = 86.2571𝛿 cos(4.60931 𝑡)

Grafico 9. Aceleración vs Tiempo para Resorte 3

Conclusion
.
Para concluir resaltamos como el movimiento del sistema masa resorte, que describe al oscilador armonico simple, responde a su naturaleza
vibratoria en forma de ondas en senos y cosenos, esto lo podemos describir gracias al modelado matematico que logramos el cual graficamos
con la ayuda de software como scilab, al modelar matematicamente un sistema podemos conocer varios aspectos del estado de este, como por
ejemplo al saber que responde a onda senoidal, podremos determinar la frecuencia, la amplitud, el periodo, que nos sirven para definir conceptos
como por ejemplo la frecuncia de resonancia, que es un fenomeno que se da cuando un objeto es sometido a vibraciones que son cercanas a su
frecuencia natural lo que puedo ocacionar que oscile de manera máxima, lo que puede ser peligroso ya que si es un objeto que cumpla una
función estructural puede colapsar o romperse. Para finalizar esta conclusión señalar que cuando determinamos la frecuencia natural del sistema
nos referimos a la vibración que tiene el sistema cuando se perturba desde el reposo.
5. Referencias Bibliográficas

1. William W. Seto. Vibraciones Mecánicas – Teoría y problemas resueltos. McGraw-Hill de México, 1970.
2. Yunus A. Cengel y John M. Cimbala-Mecánica de los fluidos, fundamentos y aplicaciones segunda edición,
Mc Graw Hill, México, 2010.
3. RAO, S. Vibraciones Mecánicas (Quinta ed.). Pearson, México, 2012.
4. Software PhET colorado, “Masses and springs”