UNIDAD III.

- SISTEMAS DE UN GRADO DE LIBERTAD

CON EXCITACIÓN ARMONICA

Sistema forzado.
Un sistema forzado es aquel que se encuentra sujeto a fuerzas o excitaciones externas. Estas
excitaciones pueden ser

1.- Armónicas
2.- Periódicas
3.- Constantes
4.- Aleatorias
5.- Choques

Las fuerzas armónicas son de las más comunes, y se representan por

F (t )  Fo sent .

La fuente más común de excitación armónica es el desbalance en las máquinas rotatorias. La

excitación armónica puede ocurrir en la forma de una fuerza o desplazamiento de algún punto
del sistema.

3.1.- Análisis de un sistema sujeto a fuerza armónica externa .

En general un sistema en vibración forzada se representa como sigue:

La ecuación diferencial para el sistema anterior es

mx  cx  kx  Fo sent -------------- (3.1)

La solución general de esta ecuación es

x  xH  xP -------------------------------- (3.2)

en donde xH = solución homogénea conocida

xP = solución particular a determinar
Para la solución particular supondremos lo siguiente:

xP  Xsen (t   ) ----------------------- (3.3)

en donde X y  son incógnitas a determinar, siendo X la amplitud de oscilación, mientras

que  es la fase del desplazamiento respecto a la fuerza excitatriz.

La amplitud y la fase de excitación se calculan sustituyendo la ecuación (3.3) en la ecuación

(3.1). Recordando que en el movimiento armónico la velocidad y la aceleración se encuentran
adelantadas con respecto al desplazamiento en 90 o y 180o respectivamente, los términos de la
ecuación diferencial se pueden representar gráficamente por:

Del diagrama anterior encontramos que

  
  2
2 2 2
Fo2  kX  m 2 X   c X   X 2  k  m 2   c   
 

Fo
X
----------------- (3.4)
 k m 2 
2 2
 c 

tan   c X  c 
( k  m 2 ) X k m 2

  tan 1  c
k m 2
 ----------------------- (3.5)
Las expresiones anteriores se pueden representar en forma adimensional considerando que
c Fo
n  k
m
, 2n  m y X est  k
.

Sustituyendo en (3.4) y (3.5) se obtiene lo siguiente:

X  1
X est 2 ---------------- (3.6)
1( / ) 2   (2 / )2
 n  n

2 /n
  tan 1 ---------------------------------- (3.7)
1( /n )2
Las ecuaciones (3.6) y (3.7) se representan gráficamente como se indica a continuación:

Figura (3.1)

Las curvas anteriores nos muestran que el factor de amortiguamiento  tiene gran influencia
sobre la amplitud y el ángulo de fase, en la región de la frecuencia próxima a la resonancia
( / n  1) .

De acuerdo con lo anterior se tienen tres casos límite:

a).-  / n  1 , (   0o )
En este caso las fuerzas de inercia y amortiguamiento son pequeñas, por lo que se traduce en
un pequeño ángulo de fase  , siendo la magnitud de la fuerza global casi igual a la fuerza del
resorte, por lo que F  kx .

b).-  / n  1 , (   90o )
En este caso la fuerza de inercia que ahora es mayor, es equilibrada por la fuerza de resorte;
mientras que la fuerza aplicada supera la fuerza de amortiguación. La amplitud de resonancia se
determina por la ecuación (3.6) quedando

F F
X res  co  2ok ------------------- (3.8)

c).-  / n  1 , (   180o )
En este caso la inercia se encarga de equilibrar la fuerza, por lo que F  mx .
En resumen la solución general de la ecuación diferencial (3.1) es

  F sen (t  )
x(t )  X1ent sen  1   2 nt  1   o
  k 1  /  2    2 /  2 ---------------- (3.9)
2
 n  n

xtran  X 1e nt sen  

1   2 nt  1 ------------------------------------------------------ (3.10)

Fo sen(t  )
xestac  
k
1  /   2    2 /   2 ---------------------------------------------------------- (3.11)
2
n n
 

3.2.- Desbalance rotatorio.

El desbalance rotatorio es una de las causas más comunes de vibración en las máquinas, y se
debe a que el centro de gravedad no coincide con el eje de rotación. Esto se puede observar en
la siguiente figura:

Consideremos el siguiente sistema resorte-masa restringido a moverse en la dirección vertical y

excitado por una masa rotatoria no balanceada tal y como se muestra en la siguiente figura:

De la figura tenemos que m = masa que gira

Movimiento de m : x  e sent
M  m = masa que no gira
Movimiento de M  m : x
La ecuación diferencial del movimiento es:
d2
( M  m) x  m 2 ( x  e sent )  kx  cx 
dt
( M  m) x  mx  kx  cx  me 2 sent 

x  Mc x  Mk x  m
M
e 2 sent ----------------------- (3.12)

Esta ecuación tiene una excitación armónica de amplitud me 2 .

La solución estacionaria de la ecuación (3.12) se representa por

X me 2
--------------------- (3.13)
 k M  
2
2
 ( c ) 2

M X   /n  2

m e 2
1  /  2   2 /  2
-------- (3.13.1)
n  n
 

2 /n
tan   ----------------------------- (3.14)
1( /n )2

Las gráficas de las ecuaciones anteriores se representan como sigue:

Figura (3.2)
La solución de la ecuación diferencial (3.12) es

x(t )  X1e nt sen  1   2 nt  1  

me 2 sen (t  )
-------------- (3.15)
 
 k M  2  (c )2
2

3.3.- Cabeceo de flechas rotatorias.

El cabeceo (whirling) es la rotación del plano realizado por el eje flexionado con respecto a la
línea de centros de los cojinetes.

Esto puede representarse como sigue:

Analizando el disco de masa m de la figura anterior tenemos:

Posición de s : ( xs , ys )
Posición de G : ( xs  e cos t , ys  e sent )

En el cabeceo sincronizado O , s y G se mantienen fijos entre sí para  constante.

d2
Ecuación del movimiento en dirección de x : m 2 ( xs  e cos t )  kxs  cx s 
dt
2
mxs  cx s  kxs  me cos t -------------- (3.16)

k = constante elástica debida a la flexión en la flecha

 d2
Ecuación del movimiento en dirección de y : m 2 ( ys  esent )   kys  cy s 
dt
2
mys  cy s  kys  me sent --------------- (3.17)

La solución estacionaria para las ecuaciones anteriores es:

me 2 cos(t  )
Xs 
--------- (a)
 k m 2 
2
(c )2

me 2sen (t  )

Ys 
---------- (b)
 k m 2 
2
(c )2

Por otro lado r  X s2  Ys2 

r me 2
----------- (3.18)
 k m 2 
2
( c )2

r no depende del tiempo

3.4.- Excitación armónica en la base.

Cuando un sistema dinámico es excitado en la base, se tiene lo siguiente:

y  Ysent = movimiento prescrito de la base

Y y  son valores conocidos

El diagrama de cuerpo libre del sistema anterior es

La ecuación diferencial del movimiento es

mx  cx  kx  cy  ky ----------------- (3.19)

Usando el álgebra compleja se tiene que

y  Ye jt ----------- ( i )
x  Xe j (t  )  Xe jt e  j --------- ( ii )

X y  son parámetros a determinar

Sustituyendo ( i ) y ( ii ) en (3.19) se obtiene

( k  jc )e j ----------- ( iii )

(k  m 2  jc) Xe j (t  )  (k  jc )Ye jt 
X 
Y ( k  m 2 ) j c

X k 2  (c )2
 
Y
 
2
k  m 2  (c ) 2

2
X 1 2 /n 
Y
 2 ---------------- (3.20)
1( / ) 2  (2 / )2
 n  n

Esta expresión representa la magnitud de estado estacionario.

Sabiendo que e j  cos   jsen , y sustituyendo en la ecuación ( iii ), ordenando y

reduciendo, se obtiene lo siguiente

tan   mc 3 

k k m 2 (c )2 
2 ( /n )3
tan   ------------------ (3.21)
1( /n )2 (2 /n )2
Las ecuaciones (3.20) y (3.21) se representan gráficamente como sigue.

Figura (3.3)

3.5.- Aislamiento de vibración y cimentación de maquinaria .

Las fuerzas vibratorias generadas por máquinas y motores son a menudo inevitables, sin
embargo su efecto en un sistema dinámico puede reducirse considerablemente mediante
resortes diseñados apropiadamente, llamados aisladores

Consideremos la siguiente figura:

Figura (3.4).- Fuerza perturbadora transmitida por los resortes y el amortiguador.

La fuerza transmitida a través de los resortes y del amortiguador es

FT   kX  2   c X  2  kX 1   c / k  
2

FT FT
kX   ---------- (a)
2 2
1 c / k  1 2 /n 
Por otro lado
X 1 Fo / k
X est
 X  
1( / ) 2 2 2 1( / ) 2 2 2
 n   (2 /n )  n  (2 / n )

Fo
kX 
2
1  /  2   2 /  2 --------------- (b)
 n  n


Igualando (a) y (b) y ordenando se obtiene la siguiente expresión:

2
FT 1 2 /n 
Fo
 -------------- (3.22)
2
1  /  2   2 /  2
 n  n


Esta relación se conoce como “transmisibilidad” de fuerza ( TR ).

La comparación de las ecuaciones (3.22) y (3.20) nos indica que FT / Fo es idéntica a X / Y ,

por lo que el problema de aislar una masa del movimiento del punto de soporte es idéntico al de
aislar fuerzas perturbadoras. Las curvas mostradas en la figura (3.3) muestran que la
transmisibilidad es menor que la unidad solo para  / n  2 , por lo que el aislamiento
vibratorio solo es posible bajo ésta condición.

Es posible reducir la amplitud de la vibración apoyando la máquina de masa m , sobre una gran
masa M como se muestra en la siguiente figura:

Cuando el amortiguamiento es despreciable, la ecuación de a transmisibilidad se reduce a

TR  1
 n  2 1
 / ------------ (3.23)

en donde se entiende que el  / n  2 .

3.6.- Instrumentos de medición de vibraciones.
Las mediciones que se van a medir pueden clasificarse como

a).- Periódicas
b).- De choque
c).- Casual o aleatoria

De estos movimientos el periódico es el más conocido, y los instrumentos utilizados para medir
la frecuencia, amplitud, velocidad, aceleración, o pendiente de onda, están bien desarrollados.

En la medición de choques, solamente son de interés los valores pico. En el caso de los
movimientos casuales, es deseable el espectro de frecuencia del valor cuadrático medio, siendo
los instrumentos utilizados para estas mediciones de gran complejidad y de reciente desarrollo.

El elemento básico de muchos instrumentos de medición es la unidad sísmica que se describe a

continuación:

La ecuación diferencial del movimiento es

mx  c( x  y )  k ( x  y ) 

mx  c( x  y )  k ( x  y )  0 ------------- (3.24)

Si z es el movimiento relativo (lectura directa), y hacemos z  x  y , entonces x  z  y , por lo

que

m( y  z )  cz  kz  0 

mz  cz  kz   my -------------------------- (3.25)

Suponiendo un movimiento sinusoidal y  Ysent del cuerpo vibrante, se obtiene la siguiente

ecuación:

mz  cz  kz  mY  2 sent ---------------- (3.26)

Esta ecuación es idéntica al caso de desbalance rotatorio, si se cambia a z por x y a mY  2

por me 2 .
La solución estacionaria es z  Zsen(t   ) , en donde la amplitud
Z mY  2 
( k m 2 )2 ( c )2

Y ( /n ) 2
Z
2
1  /  2   2 /  2 --------------- (3.27)
n  n
 

El ángulo de fase se determina por

tan   c 2 
k  m

2 / n
tan   2 --------------------- (3.28)
1  /n 

A partir de las ecuaciones anteriores se obtienen las siguientes gráficas:

Sismómetro.
Es un instrumento de n pequeña con respecto a la frecuencia de vibración  que se va a
medir; esto es  / n  1 , por lo que Z / Y  1.0 , sin importar el valor del amortiguamiento.
La masa m permanece entonces estacionaria mientras que la caja portante se mueve con el
cuerpo vibrante.

Una desventaja del sismómetro es su gran tamaño. Como Z  Y , el movimiento relativo de la

masa sísmica debe ser del mismo orden de magnitud que el de la vibración que se va a medir.

El movimiento relativo z es generalmente convertido en un voltaje haciendo que la masa

sísmica sea un magneto que se mueve relativamente a bobinas fijas a la caja. Como el voltaje
generado es proporcional a la razón de corte del campo magnético, la salida del instrumento es
proporcional a la velocidad del cuerpo vibrante, razón por la cual se les conoce como velómetros.
El voltaje de salida se determina por

Es  Sv  S ( z ) ----------- (3.29)

en donde S = sensibilidad
v   z = velocidad del cuerpo vibrante

Un sismómetro tiene las siguientes características:

Frecuencia natural: n  (1 a 5 Hz)
Rango útil de frecuencias: 10 a 2000 Hz
Sensibilidad: 20 mV/cm/s a 350 mV/cm/s

Acelerómetro.
Es un instrumento de alta frecuencia natural comparada con la de la vibración que se va a medir,
por lo que indica aceleración.

2
Un examen de la ecuación (3.27) muestra que el factor 1    /   2    2 /   2  1
 n 
 n

para  / n  0 , de modo que

2 ( aceleración )
Z   2Y  ------------- (3.30)
n n2
Acelerómetros de cristal piezoeléctrico .
Este tipo de acelerómetros utilizan cristales de cuarzo o titanato de bario, y se utilizan para
medidas de alta frecuencia natural. Los cristales están montados de manera que bajo
aceleraciones, se comprimen o se deflectan para generar carga eléctrica.

La frecuencia natural de tales acelerómetros puede hacerse muy alta ( f n  50, 000 Hz ), lo
que permite medidas hasta de 3000 g . (g=gravedad).
El tamaño del acelerómetro es muy pequeño, cerca de 1 cm de diámetro y de altura, siendo muy
resistente ya que puede soportar choques hasta de 10,000 g.

La sensibilidad S del acelerómetro está dada en términos de carga (picocoulombs/g), o en

términos del voltaje (milivolts/g). Si S  Q (carga en coulombs), y CT es la capacitancia del
cristal, incluyendo la capacitancia “ Shunt” del cable conector, entonces la salida de voltaje en el
circuito abierto se reduce a

Es  S  S
CT C1 C2 ----------------- (3.31)

en donde S = sensibilidad del acelerómetro

C1 = capacitancia del cristal
C2 = capacitancia del cable
Es = voltaje de salida en mV/g (milivolts/gravedad)

Ejemplo 1.- Un medidor de vibraciones tiene una sensibilidad S  30 mV/cm/s . Suponiendo

que la precisión límite del aparato es de 3 mV (rms), determinar el límite superior de frecuencia
del instrumento para una excitación de una gravedad ( 1 g). ¿Qué voltaje debe generarse a 200
Hz?

Solución:
Aceleración: a  g  981 cm/s 2
Sensibilidad: S  30 mV/cm/s

Es 3   z  0.1 cm/s = v
En el límite de precisión Es  S ( z )  z  S
 30 lím
g g
vlím    2 f  f  2981
(0.1)
 1561.31 cps 

f  1561.31 cps

981
Para f  200 cps , v  2 (200)  0.78065 cm/s
Es  Sv  30  0.78065  23.42 mV 

Es  23.42 mV
Ejemplo 2.- Un motor de 100 kg está montado sobre cuatro resortes, cada uno de constante
k  90 kN/m , y está conectado al suelo por medio de un amortiguador con c  6500 N.s/m .
El motor está restringido a moverse verticalmente, y se observa que la amplitud de su
movimiento es de 2.1 mm a una velocidad de 1200 rpm. Si la masa del rotor es de 15 kg,
determinar la distancia entre el centro de masa del rotor y el eje de la flecha.

  40 rad/s
Solución: M  100 kg , m  15 kg , x  2.1 mm ,   1200
30
Constante elástica equivalente: ke  4  90  360 kN/m
ke 360,000
n  M
 100
 60 rad/s

   
2
  40  2.094    4.38648
n 60 n
c 6,500
2n  M
  210060
 0.54167

M X   / n  2 2 2
MX 1  /n    2 /n 
2
m e
   
2
1  /  2   2 /  2 e 2

 n 
 n m   / n 
2 2
1000.0021  1 4.38648   20.541672.09439  e  13.01 mm
e 154.38648
 0.01301 m 

Ejemplo 3.- Un vibrómetro cuyo amortiguamiento es despreciable, tiene una frecuencia natural
de 31.4 rad/s y se emplea para medir la amplitud de vibración de una máquina. Si el vibrómetro
da una lectura del desplazamiento relativo de 0.08 pul, ¿Cuál es la amplitud de vibración de la
parte de la máquina considerada, si se encuentra girando a 200 rpm?

 
2
Solución: n  31.4 rad/s ,   2 60
200  20.9438 rad/s
,   0.44489 ,   0
n
Y ( /n )2
Z  10.44489Y 
0.55510.08
1  /  2 2 2
0.44489 Y 0.44489
 0.09981 pul 
n   2 /n 


Y  0.09981 pul

Ejemplo 4.- Una parte de una máquina de masa m  4 kg , vibra en un medio viscoso.
Determine el coeficiente de amortiguamiento c , si una fuerza excitatriz armónica de 30 N genera
una amplitud de resonancia de 1.5 cm, con periodo de 0.3 seg. Encontrar el valor de k .
 
Solución: Para n  1 , X res  1.5 cm = 0.015 m ,   n  2  0.3
2  20.944 rad/s

k  mn2  4(20.944) 2  1754.59 N/m  k  1754.59 N/m

F F 30
X res  co  c   Xo  20.9440.015
 95.5  c  95.5 N.s/m
res

Ejemplo 5.- Obtenga la constante del resorte equivalente para la base de una rectificadora de
corriente directa que gira a 300 rpm y pesa 50 lb, la cual reduciría la fuerza transmitida a la base
a una quinta parte de su valor. Suponga   0 .

50  1.55279 lb.s 2 /pie (slugs)   2 300  31.4159 rad/s

Solución: m  32.2 , 60 ,
TR  15  0.2
2
Para   0 , TR   / 2 1    / n   1  TR   / n  1  5 
1 1
 n
n  31.4159  12.825 rad/s
6

k  mn2  1.55279  (12.825) 2  255.4 lb/pie  k  21.28 lb/pul

Ejemplo 6.- Un sistema resorte masa es excitado por una fuerza F  Fo sent . En resonancia
la amplitud medida es de 0.6 pul. A 0.5 de la frecuencia de resonancia, la amplitud es 0.4 pul.
Hallar el factor de amortiguamiento  del sistema.

Solución:

X
A resonancia X res  2o  X o  2 X res  2 (0.6)  X o  1.2 ------- (a)

 
2 2
X
Xo
 1  X o  0.4 1  0.52   2  0.5  
A 0.5 de resonancia 2 2

1

 
n  (2 /n )

2

X o  0.4 0.5625   2 ----------- (b)

Igualando (a) y (b) se obtiene lo siguiente:

1.2  0.4 0.5625   2  9 2  0.5625   2   2  0.0703125    0.26516

Ejemplo 7.- Una máquina que pesa 200 lb es soportada por cuatro resortes de rigidez k lb/pul
cada uno. Si la unidad opera a 1200 rpm, ¿cual debe ser el valor de la constante k si solo el
20% de la fuerza excitatriz de la unidad debe ser transmitida a la estructura que lo soporta?
Solución:   0 , m  200
386
 0.5181347 lb.s 2 / pul ,   2 60
1200  125.6636 rad/s

TR  1    / n  2  1  0.2
1  6   /   2.44949 
n n  51.3 rad/s
  /n  2 1
kT  mn2  0.5181347  (51.3)2  1363.57 lb/pul

La rigidez para cada resorte es k  1363.57

4
 340.89 lb/pul  k  340.89 lb/pul

Ejemplo 8.- Un medidor de vibraciones sin amortiguamiento, con una frecuencia natural de 1 cps
es utilizado para medir vibraciones armónicas de 4 cps. Si la amplitud indicada por el medidor
(amplitud relativa entre la masa del medidor y el marco) es 0.052 cm, ¿Cuál es la amplitud
correcta?

Solución: Z  0.052 cm ,   0
Z 1  /n  
2
Y ( /n )2 ( /n )2 Y  
Z  2
Y 2

2
1  /  2   2 /  2 1  /n  ( /n )
 n  n

 1  2 
 
  
2
Y  Z  n  1  0.052    1  0.04875 cm  Y  0.04875 cm
   4  

Ejemplo 9- Un instrumento sensible de 113 kg debe instalarse en un sitio donde la aceleración es

de 15.24 cm/s2 a una frecuencia de 20 cps. El instrumento se coloca en una cama de caucho con
las siguientes propiedades: k  2802 N/cm y   0.1 , ¿Qué aceleración se transmite al
instrumento?

Solución: n  k  280200  49.796 rad/s ,   2  20  125.6636 rad/s

m 113

 
125.6636 2
   2.52356 , 
n n
 6.3684
49.796
Aceleración   2Y  15.24 cm/s 2  Y  15.24  0.000965 cm
(125.6636)2
2 2
X 1 2 /n  1 20.12.52356 
Y
 2
  0.207 
1( / )2  (2 / )2
n  n
 16.3864  2 (20.12.52356)2


X  0.207Y  0.207  0.000965  0.000199772 cm

Aceleración transmitida al instrumento:

atrans   2 X  (125.6636)2  0.000199772  3.15468 cm/s 2  atrans  3.15468 cm/s 2

Ejemplo 10- Para el sistema mostrado, determinar la amplitud de la vibración, y el ángulo de

fase, de acuerdo con los siguientes datos: k  800 N/m , m  8 kg , c  20 N.s/m
Solución:

Ecuación diferencial del movimiento: mx  cx  kx  Fo sent

n  k  800  10 rad/s
m 8
c 20  0.125
2n  m   2108

Fo / k 10/800
X   0.00411 m 
2 2 2
X  4.11 mm

1

 
n  (2 /n )

2
 
1 20 2  
10 

 20.12520/10  2

tan  
2    20.1252  0.16666 

n
  9.46o
1 
214
n