Universidad de Santiago de Chile

Facultad de Ingeniería
Departamento de Ingeniería Mecánica y Climatización

LABORATORIO C236
“VIBRACIONES EN UN GRADO DE LIBERTAD”

Nombre: Oscar Ruz Núñez

Profesor: Marcela Cruchaga

Asignatura: Tópicos III, Sólidos

 Índice
1. Resumen ...................................................................................................................................... 3
2. Objetivos ...................................................................................................................................... 4
3. Metodología de trabajo ................................................................................................................ 5
4. Presentación de datos y los resultados ....................................................................................... 6
a. Ecuación de movimiento para cada situación ......................................................................... 6
5. Conclusiones .............................................................................................................................. 13
6. Apéndice .................................................................................................................................... 20
a. Planteamiento ecuación general de movimiento .................................................................. 20
b. Solución analítica .................................................................................................................. 20
c. Bibliografía............................................................................................................................. 24

2
1. Resumen

En el presente informe se estudia analítica y numéricamente la respuesta de un

sistema particular sometido a vibración. En particular, se analiza el efecto de este
sistema al variar diferentes parámetros característicos. Por último, se concluye respecto
a la influencia de cada parámetro en estudio sobre la respuesta del sistema.

3
2. Objetivos

Objetivo general:

- Estudiar la respuesta de un sistema propuesto al variar parámetros

característicos de este.

Objetivos específicos:

- Establecer ecuación diferencial de movimiento para cada situación, indicando la

frecuencia natural en cada una de ellas.
- Indicar la frecuencia analítica para cada caso y contrastarla con el método
numérico de diferencias finitas.
- Analizar el comportamiento del movimiento de acuerdo con diversas variaciones
de parámetros.

4
 3. Metodología de trabajo

A continuación, se presenta el problema de movimiento que se aborda en este

laboratorio.

Se desea estudiar el movimiento de una viga que posee un resorte y un amortiguador

acoplado en sus extremos, en conjunto de un motor con masas descentradas ubicado a la
mitad de su longitud (véase figura 1).

Ilustración 1, Problema en estudio.

A partir del modelo propuesto en la figura 1, se analizan los siguientes casos:

Tabla 1, Descripción casos en estudio

Caso N° 𝑁 𝑁𝑠 𝑊𝑛 [𝑅𝑃𝑀]
𝑘[ ] 𝑐[ ]
𝑚 𝑚
a) 2500
I b) 4500 - -
c) 7500
a) 4500 100
II b) 4500 1000 -
c) 4500 Crítico
a) 4500 100
III b) 4500 - 335
c) 4500 650
IV a) 4500 3000 335
b) 4500 Crítico 335

5
 4. Presentación de datos y los resultados

A continuación, se desarrollan los puntos expuestos anteriormente.

a. Ecuación de movimiento para cada caso

La ecuación de movimiento del sistema completo, considerando la acción del motor, la

amortiguación y el resorte, se presenta a continuación (ecuación 1) (véase apéndice).

4𝐿2 𝑐 4𝐿2 𝑘 2 ∙ 𝐿 ∙ 𝑚 ∙ 𝑒 ∙ 𝑤𝑓2 (1)

𝜃̈ + ∙ 𝜃̇ + ∙𝜃 = ∙ cos⁡(𝑤𝑓 𝑡)
𝐽 𝐽 𝐽
De donde,

L, es la mitad de longitud de la viga.

m, es la masa excéntrica de los motores.
c, es el valor de amortiguación del sistema.
k, es la rigidez del sistema.
J, es el momento polar del sistema.
e, es la distancia de la masa excéntrica respecto del centro del motor.
𝑤𝑓 , es la frecuencia del motor.
𝜃, es la posición angular de la viga.

A partir de la ecuación general de movimiento, es posible desprender las diversas

ecuaciones de movimiento para las situaciones específicas. La siguiente tabla presenta las
ecuaciones diferenciales para cada situación en estudio.

Tabla 2, Ecuaciones movimiento generales.

Sistema en estudio Ecuación diferencial

Viga únicamente con 4𝐿2 𝑘 (2)
resorte. Caso I 𝜃̈ + ∙𝜃 =0
𝐽
Viga con resorte y 4𝐿2 𝑐 4𝐿2 𝑘 (3)
amortiguador. Caso II 𝜃̈ + ∙ 𝜃̇ + ∙𝜃 =0
𝐽 𝐽
Viga con resorte y 4𝐿2 𝑘 2 ∙ 𝐿 ∙ 𝑚 ∙ 𝑒 ∙ 𝑤𝑓2 (4)
motor con masas 𝜃̈ + ∙𝜃 = ∙ cos⁡(𝑤𝑓 𝑡)
𝐽 𝐽
descentradas
funcionando. Caso III

A su vez, es posible apreciar de las ecuaciones (1), (2), (3) y (4), que independientemente
del sistema, la frecuencia natural permanece constante, y está dada por la siguiente relación:

(2)
4𝐿2 𝑘
𝑤𝑛 = √
𝐽

6
 Dependiendo del valor de los coeficientes, así como de la naturaleza misma de la ecuación,
existen diversos resultados posibles para las diversas ecuaciones diferenciales [1]. En la siguiente
tabla se presentan las soluciones analíticas para cada caso en estudio.

Tabla 3, Resumen ecuaciones generales sistemas de vibración.

Caso Naturaleza Ecuación General

Viga con resorte, Ecuación diferencial (1) no
amortiguación y carga homogénea, con c = 3000, 𝜃(𝑡) = 𝐴 ∙ 𝑒 −𝜉𝑤𝑛 𝑡 sen(𝑊𝑑 ∙ 𝑡 + 𝛼)
desbalanceada raíces reales. + 𝑋𝑐𝑜𝑠(𝑤𝑛 𝑡 − 𝜃)
(Vibración forzada Ecuación diferencial (1) no
com amortiguación) homogénea, con c = critico
raíces reales iguales.
Viga con únicamente Ecuación diferencial (2) 𝜃(𝑡) = 𝐴 ∙ cos⁡(𝑊𝑛 ∙ 𝑡 − 𝛼)
resorte homogénea, con raíces
(Vibración libre) imaginarias.
Ecuación diferencial (3)
homogénea, con c = 100
raíces imaginarias.
Viga con resorte y Ecuación diferencial (3)
amortiguación homogénea, con c = 1000 𝜃(𝑡) = 𝐴 ∙ 𝑒 −𝑛𝑡 ∙ cos⁡(𝑊𝑑 ∙ 𝑡 − 𝛼𝑑 )
(Vibración libre raíces reales.
amortiguada) Ecuación diferencial (3)
homogénea, con c = critico
raíces reales iguales.
Viga con resorte y carga Ecuación diferencial (4) no 𝜃(𝑡)
desbalanceada homogénea, con raíces −1
𝑄 𝑊𝑓 ⁡ 2
(Vibración forzada sin imaginarias. = ∙ ⁡[1 − ( ) ] ∙ 𝑠𝑒𝑛(𝑊𝑓 ∙ 𝑡)
𝑘 𝑊𝑛 ⁡
amortiguación)
𝑥̇ 𝑜
+ 𝑥𝑜 ∙ cos(𝑊𝑛 ∙ 𝑡) + sin(𝑊𝑛 ∙ 𝑡)
𝑊𝑛
𝑄
𝑚𝑜 𝑊𝑓
+ 2 2 [sin(𝑊𝑓 ∙ 𝑡) − 𝑊
(𝑊𝑛 − 𝑊𝑓 ) 𝑛

∙ sin(𝑊𝑛 ∙ 𝑡)]

b. Método numérico

En contraste con la solución analítica, se realiza una solución numérica para cada caso
estudiado, a través del método de diferencias finitas.

Discretizando la ecuación general (1), a través del método de diferencias finitas, queda la
siguiente expresión:

𝜃𝑖+1 − 2𝜃𝑖 + 𝜃𝑖−1 4𝐿2 𝑐 𝜃𝑖+1 − 𝜃𝑖−1 4𝐿2 𝑘 2 ∙ 𝐿 ∙ 𝑚 ∙ 𝑒 ∙ 𝑤𝑓2 (7)

2
+ ∙ + ∙ 𝜃𝑖 = ∙ cos⁡(𝑤𝑓 𝑡)
∆𝑡 𝐽 2∆𝑡 𝐽 𝐽

7
 Es posible reordenar la ecuación (7), de tal manera de dejar como incógnita el valor futuro de
𝜃, expresando así una ecuación de tipo explícita, la cual calcula el valor de 𝜃 para cada paso de
tiempo directamente. En el apéndice se describe la programación del método numérico utilizado.

c. Caso I

A partir del cálculo de la amplitud, la frecuencia y el desfase según las fórmulas indicadas en el
apéndice, y comprendiendo que se trata de un caso de vibración libre (tabla 1), en el siguiente
recuadro se expresan las correspondientes ecuaciones de movimiento:

Tabla 4, Ecuaciones de movimiento caso I

Caso I: Vibración libre

N° 𝐴⁡[𝑚] 1 𝐹𝑢𝑛𝑐𝑖ó𝑛
𝑊𝑛 [ ]
𝑠
a) 0.041 10 𝑥(𝑡) = 0,041 ∙ 𝑐𝑜𝑠(10𝑡 − 1,326)
b) 0.031 13.416 𝑥(𝑡) = 0,031 ∙ 𝑐𝑜𝑠(13,416𝑡 − 1,247)
c) 0.025 17.321 𝑥(𝑡) = 0,025 ∙ 𝑐𝑜𝑠(17,321𝑡 − 1,162)

En la siguiente figura se representan las soluciones analíticas, según el valor de la constante de

rigidez del sistema, expresadas en la tabla 4.

Ilustración 2, comparación solución analítica diferentes k.

En tanto, es posible contrastar los resultados numéricos con los resultados analíticos para cada
subcaso. En particular, la elección del tamaño de paso es de vital importancia para que ambos
resultados sean idénticos. La importancia de la adecuada elección de un tamaño de paso se ilustra

8
en la figura 3, donde se compara la solución analítica del caso k = 2500 con la solución numérica
para un tamaño de paso grande.

Ilustración 3, Comparación solución numérica versus analítica paso de tiempo grande.

Refinando el tamaño del paso, es posible encontrar soluciones numéricas lo suficientemente

cercanas a la solución analítica. En la siguiente figura, se compara la solución analítica con la
numérica de cada subcaso para un paso de tiempo pequeño.

(a)

9
 S

(b)
S

(c)
Ilustración 4, comparación solución analítica versus numérica. a) k= 2500 b) k=4500 c) k=7500

10
d. Caso II

El segundo caso consta de vibración libre amortiguada. De tal manera, dependiendo de la

naturaleza del amortiguamiento, es que existirán diferentes ecuaciones que describen el
movimiento.

Tabla 5, Ecuaciones de movimiento caso II

Caso II: Vibración libre amortiguada

N° 𝜉 𝐸𝑠𝑡𝑎𝑑𝑜 𝐹𝑢𝑛𝑐𝑖ó𝑛
a) 0.149 Amortiguamiento 𝑥(𝑡) = 0,033 ∙ 𝑒 −2𝑡 ∙ cos⁡(13,266 ∙ 𝑡 − 1,265)
subcrítico
𝑥(𝑡)
= 𝑒 −20∙𝑡 ∙ [0,01
∙ 𝑐𝑜𝑠 (13,416 ∗ √1 − 1,4912
b) 0.031 Sobre 0,6
amortiguamiento +( )
13,416 ∗ √1 − 1,4912
∙ 𝑠𝑒𝑛 (13,416 ∗ √1 − 1,4912 𝑡)]

c) 0.025 Amortiguamiento 𝑥(𝑡) = 𝑒 −13,416∙𝑡⁡ (0,01 + (0,53) ∙ 𝑡)

crítico

La figura siguiente muestra el movimiento del sistema versus el tiempo, según las ecuaciones de
movimiento expresadas en la tabla 5.

Ilustración 5, comparación analítica para diferentes casos de amortiguación

11
A su vez, considerando un paso de tiempo pequeño, el contraste entre solución numérica y
analítica se presenta en la siguiente figura.

(a)
S

(b)

12
 S

(c)
Ilustración 6, comparación solución analítica versus numérica. a) c= 100 b) c=1000 c) c=crítica

e. Caso III

El caso tres es sobre vibración forzada sin amortiguación. La siguiente tabla muestra los
parámetros más influyentes del problema actual. (Debido a su extensión, en la tabla 9 del apéndice
se presentan los valores de los coeficientes de la ecuación diferencial completa).

Tabla 6, Ecuación de movimiento caso III

Caso III: Vibración forzada sin amortiguación

N° 𝑊𝑓 [𝑅𝑃𝑀] 𝑁 𝐹𝑢𝑒𝑟𝑧𝑎
𝑘[ ]
𝑚
a) 100 4500 0,0005 ∙ (10,47)2 ∙ 𝑐𝑜𝑠(10,47 ∙ 𝑡)
b) 335 4500 0,0005 ∙ (35,08)2 ∙ 𝑐𝑜𝑠(35,08 ∙ 𝑡)
c) 650 4500 0,0005 ∙ (35,08)2 ∙ 𝑐𝑜𝑠(35,08 ∙ 𝑡)

Las gráficas de movimiento, según las soluciones de la tabla 6, se presentan en la figura 7.

13
 S

Ilustración 7, Comparación solución analítica diferentes frecuencias motor.

Además, considerando un paso de tiempo pequeño, en la siguiente figura se comparan los

resultados numéricos respecto a los resultados analíticos.

(a)

14
 S

(b)
S

(c)
Ilustración 8, Comparación analítica versus numérica. a) wf = 100 b) wf = 335 c) wf = 650

f. Caso IV

El caso cuatro trata sobre vibración forzada con amortiguamiento. La tabla siguiente
representa los parámetros más relevantes para el cálculo de las ecuaciones (debido a su
extensión, en la tabla 9 del apéndice se presentan los valores de los coeficientes de la ecuación
diferencial completa).

15
Tabla 7, Ecuación de movimiento caso IV

Caso IV: Vibración forzada con amortiguación

N° 𝑊𝑓 [𝑅𝑃𝑀] 𝑁𝑠 𝐹𝑢𝑛𝑐𝑖ó𝑛
𝑐[ ]
𝑚
a) 3000 0,0005 ∙ (35,08)2 ∙ 𝑐𝑜𝑠(35,08 ∙ 𝑡)
b) 335 Crítico 0,0005 ∙ (35,08)2 ∙ 𝑐𝑜𝑠(35,08 ∙ 𝑡)

A partir de los resultados previamente descritos, la gráfica de posición versus tiempo para el
sistema completo se representa en la siguiente figura:

Ilustración 9, Comparación solución analítica diferentes coeficientes amortiguación motor.

En tanto, la respuesta estacionaria del problema, a saber, en un tiempo lejano, se representa en

la figura siguiente.

16
 S

Ilustración 10, Comparación solución analítica diferentes coeficientes amortiguación motor; Solución
estacionaria.

El presente problema posee complicaciones adicionales en comparación a los problemas previos

numéricamente hablando. Visto desde otra manera, se necesita un paso de tiempo aún más
preciso para lograr describir el problema completamente. A continuación, se presenta la
comparación de la solución analítica contra la solución numérica para el problema con c = 3000
con un paso de tiempo moderado.

Ilustración 11, Comparación solución numérica versus analítica paso de tiempo mediano.

17
Por último, la comparación entre los resultados numéricos (para un paso de tiempo más
pequeño) contra el resultado analítico se representa en la figura 12.

(a)
S

(b)
Ilustración 12, Comparación analítica versus numérica. a) c = 3000 b) c = crítico

18
 5. Conclusiones

En lo concerniente a la viabilidad de
aplicación del método de diferencias finitas, los resultados demuestran que dicho método
produce soluciones idénticas a la solución analítica para un paso de tiempo pequeño. Sin
embargo, para pasos de tiempo mayores, dicho método se aleja de los resultados
esperables; Es debido a lo previamente descrito que se desprende la importancia de una
correcta elección de paso de tiempo, con el objeto de obtener resultados apegados a la
realidad. Dicho de otro modo, posterior a cualquier simulación, es menester cerciorarse de
haber elegido el paso de tiempo adecuado.

Respecto a los resultados analíticos para el caso I, es posible observar que, mientras menor
sea el valor de la constante de rigidez, mayor será la amplitud del movimiento del sistema.
En contraparte, a mayor rigidez, menor será el periodo del mismo, o mayor su frecuencia
natural. Por último, se aprecia un comportamiento sinusoidal marcado, lo cual concuerda
con los resultados expresados en la tabla 4.

En tanto, a partir de la figura 5, es posible desprender las siguientes aseveraciones para el

caso II. En una comparación entre los coeficientes de amortiguación c = 100 y c = 1000,
es posible observar que, a mayor valor, la respuesta será completamente exponencial. En
contra parte, para el caso de c = 100, se aprecia que se poseen oscilaciones mayores antes
de alcanzar el amortiguamiento completo. Para el caso crítico, cuyo valor de c es menor
que para c = 1000, se aprecia un decaimiento brusco de la respuesta, inclusive mayor que
para el caso con mayor valor (c critico = 670.8). Visto desde un ámbito ingenieril, lo anterior
refleja que, para aplicaciones donde se desee que la respuesta se amortigüe
inmediatamente, es aconsejable poseer un dispositivo con amortiguamiento cercano o igual
al crítico.

En cuanto a los resultados expresados en la figura 7, se observa que la respuesta se

comporta de manera similar para las diferentes rotaciones del motor con cargas
desbalanceadas, contrario a lo que se tendería a pensar, suponiendo que, a mayor rotación
del motor, mayor será la amplitud (o frecuencia natural) del sistema. Dicho de otra manera,
si se desea perturbar un sistema, la rotación del motor acoplado no es el parámetro más
influyente, sino la masa desbalanceada del mismo.

Por último, a partir de la figura 12, relacionada con el caso IV, es posible observar que existe
una respuesta exponencial al comienzo para ambos casos, para luego en una situación
amortiguada (o “estacionaria”), el sistema se comporte de forma sinusoidal. Respecto a los
resultados relacionados con los valores de los coeficientes, es posible concluir que, aún
con un valor de amortiguamiento elevado, el sistema decae (o amortigua) mucho más
rápido para el caso de amortiguamiento crítico. Visto desde un punto de vista económico y
practico, lo anterior implica que, si se desea una aplicación con amortiguamiento elevado,
no es necesario buscar dispositivos con un alto coeficiente de amortiguación, sino uno que
posea un coeficiente cercano al valor crítico del sistema.

19
 6. Apéndice

a. Planteamiento ecuación general de movimiento

Aplicando un momento desde el soporte de la viga, es posible plantear la siguiente ecuación:

𝐽𝜃̈ = −𝐹𝑐 ∙ 2𝐿 ∙ cos(𝜃) − 𝐹𝑘 ∙ 2𝐿 ∙ cos(𝜃) + 𝐹𝑒 ∙ 𝐿

Considerando las siguientes relaciones geométricas, hipótesis y expresiones:

𝜃 → 0; cos⁡(𝜃) → 1

𝐹𝑐 = 𝑥̇ ∙ 𝑐; 𝐹𝑘 = 𝑥 ∙ 𝑘

𝐹𝑒 = 2 ∙ 𝑚 ∙ 𝑒 ∙ 𝑤𝑓2

𝑥
𝜃 → 0; sen(𝜃) → 𝜃; 𝑠𝑒𝑛(𝜃) = → 𝑥 = 𝜃 ∙ 2𝐿⁡;⁡𝑥̇ = 𝜃̇ ∙ 2𝐿
2𝐿

Reemplazando los anteriores términos en la expresión general, y reordenando, se obtiene:

4𝐿2 𝑐 4𝐿2 𝑘 2 ∙ 𝐿 ∙ 𝑚 ∙ 𝑒 ∙ 𝑤𝑓2

𝜃̈ + ∙ 𝜃̇ + ∙𝜃 = ∙ cos⁡(𝑤𝑓 𝑡)
𝐽 𝐽 𝐽

b. Solución analítica

Al ser ecuaciones diferenciales de segundo orden, la solución de cada ecuación dependerá

de si su naturaleza es o no homogénea y del valor de los coeficientes en la ecuación característica.

De tal manera, en principio, se emplea el siguiente programa de Matlab para determinar si

los valores de la ecuación característica dan valores reales o imaginarios, a partir del cual se
determinan las soluciones características para cada caso en estudio:

20
 Ilustración 13, Código determinador de naturaleza de raíces ecuación característica

A partir de los valores obtenidos previamente con el programa, es posible representar la siguiente
tabla con la solución de las ecuaciones diferenciales para cada caso.

Tabla 8, Tabla resumen soluciones generales ecuaciones diferenciales.

Caso Naturaleza Ecuación General

Viga con resorte, Ecuación diferencial (1) no 𝜃(𝑡) = 𝑐1 𝑒 𝜆1𝑡 + 𝑐2 𝑒 𝜆2𝑡 + 𝑣1 (𝑥)𝑒 𝜆1𝑡
amortiguación y carga homogénea, con c = 3000, + 𝑣2 (𝑥)𝑒 𝜆2𝑡
desbalanceada raíces reales.
Ecuación diferencial (1) no 𝜃(𝑡) = 𝑒 𝜆𝑡 (𝑐1 + 𝑐2 𝑡 + 𝑣1 (𝑥)
homogénea, con c = critico + 𝑣2 (𝑥)𝑡)
raíces reales iguales.
Viga con únicamente Ecuación diferencial (2) 𝜃(𝑡) = 𝑒 𝑎𝑡 (𝑐1 cos(𝑏𝑡) + 𝑐2 sen(𝑏𝑡))
resorte homogénea, con raíces
imaginarias.
Ecuación diferencial (3) 𝜃(𝑡) = 𝑒 𝑎𝑡 (𝑐1 cos(𝑏𝑡) + 𝑐2 sen(𝑏𝑡))
homogénea, con c = 100 raíces
imaginarias.
Viga con resorte y Ecuación diferencial (3) 𝜃(𝑡) = 𝑐1 𝑒 𝜆1𝑡 + 𝑐2 𝑒 𝜆2𝑡
amortiguación homogénea, con c = 1000
raíces reales.

21
 Ecuación diferencial (3) 𝜃(𝑡) = 𝑒 𝜆𝑡 (𝑐1 + 𝑐2 𝑡)
homogénea, con c = critico
raíces reales iguales.
Viga con resorte y carga Ecuación diferencial (4) no 𝜃(𝑡) = 𝑒 𝑎𝑡 (𝑐1 cos(𝑏𝑡) +
desbalanceada homogénea, con raíces 𝑐2 sen(𝑏𝑡) + 𝑣1 (𝑥) cos(𝑏𝑡) +
imaginarias. 𝑣2 (𝑥)sen(𝑏𝑡))

Cabe señalar que, a partir de manipulaciones algebraicas, se desprenden las ecuaciones

expresadas en la tabla 1. De este modo, dichas soluciones son generales.

A su vez, a través de un programa en Matlab, es posible encontrar las constantes para cada
ecuación, la cual se resume en la siguiente tabla.

Tabla 9, Tabla con valor coeficientes constantes ecuaciones diferenciales tabla 8.

Caso 𝜆1 /𝑎 𝜆2 /𝑏 𝑐1 𝑐2
a) 0 10 0.025 0.1
I b) 0 13.41 0.025 0.074
c) 0 17.32 0.025 0.057
a) -2 13.26 0.025 0.079
II b) -5.16 -34.83 0.063 -0.0381
c) -13.41 0.025 1.33
a) 0 13.41 0.025 0.062
III b) 0 13.41 0.025 0.105
c) 0 13.41 0.025 0.127
IV a) -2.55 -117.44 0.025 -0.008
b) -13.41 -

A continuación, se resumirá las fórmulas empleadas para calcular los diferentes parámetros
respecto a las ecuaciones de movimientos planteadas en la tabla 1.

Tabla 10, Tabla valores coeficientes ecuaciones genéricas de movimiento presentadas en tabla I.

Caso A 𝜙
Vibración libre 2 −1
𝑥̇ 𝑜
𝑥̇ 𝑜 tan ( )
√𝑥𝑜 2 +( ) 𝑥𝑜 ∙ 𝑊𝑛
𝑊𝑛

Vibración libre con 𝑥̇ 𝑜 + 𝑥𝑜 ∙ 𝑛

(𝑥̇ 𝑜 + 𝑥𝑜 ∙ 𝑛)2 tan−1 ( )
amortiguación √𝑥𝑜 +2 𝑥𝑜 ∙ 𝑊𝑑
𝑊𝑑 2
Vibración forzada 𝑞 2∙𝜉∙𝜓
𝑊𝑛 2 tan−1 (⁡ ⁡)
1 − 𝜓2
√4 ∙ 𝜉 2 ∙ 𝜓 2 + (1 − 𝜓 2 )2
Vibración forzada con 𝑥0 − 𝑋cosθ 𝑤𝑑 (𝑥0 − 𝑥 cos 𝜃)
amortiguación arctg ( )
sin 𝜙 𝑥̇ 0 + (𝑥0 − 𝑋 cos 𝜃)𝜉𝑤𝑛 − 𝑤𝑋𝑠𝑒𝑛𝜃

22
Para vibración forzada:

𝑓0
𝑋=
√(𝑤𝑛2 − 𝑤 2 )2 + (2𝜉𝑤𝑛 𝑤)2

2𝜉𝑤𝑛 𝑤
𝜃 = tan−1
𝑤𝑛2 − 𝑤 2

c. Programación método numérico

23
Ilustración 14, Código numérico de diferencias finitas para la resolución de las diferentes ecuaciones
de movimiento.

d. Bibliografía

[1] Singiresu Rao. (2012). Vibraciones mecánicas. México: Pearson.

[2] Viviana Meruane. (2016). Vibraciones mecánicas. Chile: Universidad de Chile.

24