capítulo 3

Oscilaciones forzadas
y resonancia
1
2
Hemos visto que, en las oscilaciones amortiguadas, la energía se disipa continuamente y la
amplitud de las oscilaciones disminuye. Para mantener en movimiento un sistema que presente
amortiguamiento, se debe introducir energía en el sistema para que el sistema continúe oscilando,
cuando se lleva a cabo esto, se dice que el oscilador es forzado. La mayoría de los sistemas que
utilizamos y vemos a diario son ejemplos de movimientos forzados, entre ellos se tiene el
péndulo de un reloj que recibe impulsos para mantenerse oscilando; en una escala mucho mayor,
la luna ejerce una fuerza gravitacional sobre la tierra que origina una fuerza impulsadora
periódica en los océanos de la Tierra causante de las mareas; a nivel microscópico, la radiación
en un horno microondas agita los electrones de las moléculas de agua contenida en los alimentos
elevando así su temperatura. Un ejemplo más simple es el movimiento que realizamos cuando
nos sentamos en un columpio y lo hacemos oscilar, el “bombeo” lo realizamos moviendo el
cuerpo y las piernas hacia delante y atrás, de forma que se convierte en un oscilador forzado. Si
se introduce energía en el sistema a un ritmo mayor del que se disipa, la energía aumenta con el
tiempo, lo cual se aprecia por un aumento de la amplitud del movimiento. Si la energía se
introduce al mismo ritmo que se disipa, la amplitud permanece constante con el tiempo.

Figura 3.1 Se puede ejercer una fuerza externa sobre un objeto sujeto a un resorte desplazando el punto del soporte
hacia arriba y hacia abajo.
Una manera de suministrar energía a un sistema formado por un objeto que cuelga de un resorte
vertical es mover el punto de soporte hacia arriba y hacia abajo como se muestra en la figura 3.1.
De igual modo, es posible suministrar energía a un péndulo simple moviendo el punto de soporte
hacia adelante y atrás. Si el punto de soporte del sistema masa-resorte (figura 3.1) o de un
péndulo simple se desplaza hacia adelante y atrás con movimiento armónico simple de amplitud
pequeña y de frecuencia angular  , el sistema empezará a oscilar; se obtiene el mismo resultado
si sobre un punto del soporte fijo de estos dos sistemas ( o sobre las masas) se aplica una fuerza
externa periódica llamada “fuerza impulsora” con frecuencia angular ω . Al principio el
movimiento de los sistemas es complicado, pero finalmente alcanzan un estado estacionario en el
que oscilan con la misma frecuencia  de la fuerza externa impulsora y con amplitud constante
y, por tanto, con energía constante también. En el estado estacionario, la energía introducida en el
sistema por la fuerza impulsora durante un ciclo es igual a la disipada en un ciclo debido al
amortiguamiento.
Analicemos ahora un fenómeno de profunda importancia en toda la física, que se presenta cuando
un sistema de este tipo, es decir, un oscilador físico, se somete a una fuerza impulsora periódica
externa, la palabra clave es "resonancia". Todo el mundo por lo menos está familiarizado
cualitativamente con este fenómeno, y probablemente la característica más curiosa de un
oscilador impulsado es el modo con que una fuerza periódica de amplitud constante produce
diferentes resultados dependientes de su frecuencia; esto es, en el oscilador forzado, la amplitud
3
de la oscilación y, por tanto, la energía del sistema en estado estacionario ( cuando la frecuencia
del oscilador se ha acoplado con la frecuencia de la fuerza externa ), depende no sólo de la
amplitud de la fuerza externa impulsora sino también de su frecuencia. Tomemos como ejemplo
k
un sistema masa-resorte cuya frecuencia de vibración natural es ω0 = , recuerde que la
m
frecuencia natural es la frecuencia que tendría el sistema si no estuviesen presentes ni el
amortiguamiento ni el sistema impulsor. Si la frecuencia impulsora es igual o aproximadamente
igual, a la frecuencia natural del sistema, este oscilará con una amplitud mucho mayor que la
propia amplitud de la fuerza impulsora, este fenómeno se denomina resonancia. Una fuerza del
mismo valor aproximadamente, con una frecuencia mucho mayor o mucho menor que la
frecuencia de resonancia es mucho menos eficaz: la amplitud producida sigue siendo pequeña.
Empecemos considerando el caso sencillo, aunque físicamente no real, de un oscilador en el que
sea totalmente despreciable el efecto del amortiguamiento.

3.1 Oscilador no amortiguado con impulsión armónica

Consideremos un sistema formado por una masa m situada en un resorte de constante k e
k
imaginemos la aplicación de una fuerza impulsora sinusoidal F=Fo cost. El valor , que
m
representa la frecuencia angular natural del sistema, se designará por o; la ecuación de
movimiento, es:
d 2x
m 2   k x  Fo cos  t
dt
o bien,
d 2x
m 2  k x  Fo cos  t (3.1)
dt
Antes de estudiar con detalle esta ecuación diferencial de movimiento consideremos
cualitativamente la situación. Si el oscilador se ve impulsado a partir de la posición de equilibrio
y se le deja oscilar por sí mismo, oscilará con su frecuencia natural, o. Una fuerza impulsora
periódica, sin embargo, intentará imponer su propia frecuencia  al oscilador. Por lo tanto, es de
esperar que el movimiento real, en este caso, se deba a la superposición de oscilaciones
correspondientes a las dos frecuencias  y o . La etapa inicial, en que los dos tipos de
movimiento son importantes, se denomina transitoria. La situación descrita anteriormente es
válida cuando sobre el sistema actúa una fuerza amortiguadora, se observa que, después de un
tiempo suficientemente largo las oscilaciones amortiguadas decrecen exponencialmente y el
único movimiento presente es el debido a la fuerza impulsora, que continuará sin disminuir con la
frecuencia , cuando se alcanza esta condición tenemos lo que se denomina movimiento
estacionario del oscilador impulsado. Cuando se resuelven ecuaciones diferenciales, a estas
soluciones transitoria y estacionaria se les conoce como solución homogénea y particular,
respectivamente. Cuando se desprecia la fuerza amortiguadora, la etapa inicial del movimiento no
es transitoria, esto es, no es despreciable, por lo que a esta solución simplemente la llamaremos
solución homogénea.
La solución matemática completa de la ecuación (3.1) es una superposición o combinación de dos
soluciones, es decir, es una suma simple de dos soluciones:

4
 x (t )  x h ( 0 )  xe ( )
dónde, xh(0) es la solución homogénea y xe() es la solución estacionaria (o particular). La
solución homogénea se encuentra cuando no hay presencia de una fuerza impulsora, o sea en la
ecuación (3.1) el término F0 = 0, es decir:

d 2x
m +k x= 0
dt 2
cuya solución corresponde a un M.A.S., esto es:
xh (t )  A cos 0t  B sen 0t
Que se puede escribir de la forma
xh (t )  C cos (0t   ) (3.2)
B
donde C  A2  B 2 y tan   . Analicemos ahora la forma que tiene la solución estacionaria,
A
en este caso F0 es diferente de cero. Se ha observado que cuando ha transcurrido un tiempo
suficientemente grande después de haber iniciado el movimiento, la fuerza impulsora F impondrá
su frecuencia al oscilador, así el sistema se moverá periódicamente con la misma frecuencia de la
fuerza externa ω. Además, la amplitud de la oscilación depende de la frecuencia de la fuerza
externa que se hace máxima para un valor de frecuencia igual a la frecuencia de resonancia.
Aumentando aún más la frecuencia, la amplitud disminuye y cuando ω es muy grande la
amplitud es muy pequeña. Este comportamiento sugiere que la solución estacionaria tiene la
forma:
xe (t )  B cos ( t   ) (3.3)
donde B es la amplitud, siempre positiva, que depende de la frecuencia de vibración y  es el
ángulo entre la fuerza impulsadora y el desplazamiento, tiene un signo menos debido a que el
desplazamiento está atrasado con respecto a la fuerza impulsadora. Se espera que su valor sea
igual a cero para bajas frecuencias y  para altas frecuencias (cuando las direcciones de la fuerza
y el desplazamiento son contrarias). Los valores de B y  se encuentran al sustituir x e y xe en la
ecuación de movimiento (3.1) que toma la forma:
m[   2 B cos ( t   )]  k B cos ( t   )  F0 cos  t

 m  2

B  k B cos ( t   )  F0 cos  t

 m  B  k B cos  t cos   sen  t sen   F cos  t

2
0

De la expresión anterior se obtiene:

 m  B  k B cos   F
2
0

 m  B  k B sen  0
2

Dividiendo las expresiones anteriores entre la masa se obtiene

  2

 02 B cos  
F
m
0
(3.4a)

  2
 02  B sen  0 (3.4b)
5
Dividiendo las expresiones (3.4) se encuentra que:
  0
tan   0   (3.5)
  
De la ecuación (3.4a) y cuando  = 0, se obtiene
F0
B ( )  (3.6)
m(0   2 )
2

válida para  < 0 ; cuando  = , se obtiene

F0
B ( )  (3.7)
m(  02 )
2

para  > 0. La solución de la ecuación de movimiento (3.1) tendrá la forma:

F0
x (t )  A cos (0t   )  cos ( t   ) (3.8)
m(   2 )
2
0

En la figura 3.2a se muestra el comportamiento de B() en función de ; cuando  = 0 el valor

F0
de la amplitud es B (0)  , esta es una amplitud impuesta por la aplicación de una fuerza
k
externa que tiene amplitud constante. Al aumentar el valor de , B() aumenta, este
comportamiento se ve de la ecuación (3.6), válida para frecuencias  < 0, ya que cuando la
frecuencia de la fuerza externa aumenta el denominador de la ecuación (3.6) se hace más
pequeño y por tanto el valor de la expresión se hace cada vez mayor. Lo contrario ocurre cuando
 > 0, de la ecuación (3.7) se observa que la amplitud disminuye con el aumento de , ya que al
aumentar la frecuencia de la fuerza impulsora el denominador se hace más grande y el valor de la
amplitud disminuye. Para valores de  muy grandes la amplitud tiende a cero. Cuando  = 0,
B(0) diverge, tiende a valores muy grandes, por lo que suele decirse que tiende a infinito, esto se
observa claramente en las ecuaciones (3.6) y (3.7) donde el denominador se hace cero para ese
valor de la frecuencia. Claramente este comportamiento no tiene un significado físico y aparece
porque se ha despreciado la fuerza de fricción o amortiguamiento.
B() (a) (b)


Fo/k

0   0  
Figura 3.2 (a) amplitud de absorción de las oscilaciones formadas en función de la frecuencia de la fuerza impulsora, en el caso
de amortiguamiento nulo. (b) Retraso de fase del desplazamiento respecto a la fuerza impulsora en función de la frecuencia.

6
En la figura 3.2b, se muestra el comportamiento de la constante de fase  en función también de
la frecuencia de la fuerza externa. Para  < 0 y considerando que el valor de  = 0 se tiene que
la fuerza impulsora está en fase con el desplazamiento (tienen la misma dirección ). Para  > 0
se tiene que  =  lo que nos dice que la fuerza impulsora está desfasada un valor igual a  con
respecto al desplazamiento (tienen direcciones opuestas). Se observa una discontinuidad
alrededor de   , un salto abrupto en el valor de delta, esto se debe a que no se consideró la
fuerza amortiguadora. En sistemas con amortiguamiento se espera variaciones más suaves de 
alrededor de  y un máximo (no una divergencia) en el gráfico de la amplitud del
desplazamiento.
Resolvamos ahora la ecuación (3.1) utilizando el método del exponente complejo; imaginemos
que la fuerza impulsora F0 cos t se obtiene como proyección sobre el eje x de un vector rotatorio
F0 e j  t , como lo muestra la figura 3.3a, o sea F0 cos  t  Re F0 e j  t  e imaginemos que x es la
proyección de un vector Z que gira con la misma frecuencia  como se ve en la figura 3.3b.

Figura 3.3 (a) representación compleja de la fuerza impulsora sinusoidal; (b) representación compleja del vector
desplazamiento en la oscilación forzada.
La ecuación diferencial (3.1) toma la forma:
d 2z
m 2
 kz  F0 e jt (3.8)
dt
Supongamos la solución z  A e j ( t  ) , con x (t )  Re z  , sustituyendo el valor de z en (3.8)
tenemos:
 m  2 A  k A  e j (t  )  F0 e j t   m  2  k  A e j t e j   F0 e j t
Dividiendo entre m:
  2
 2
0 A 
F0  j 
m
e

Que se puede escribir como

F0  j
(02  2 ) A  e
m
F0 F
(02  2 ) A  cos   j 0 sen (3.9)
m m
En esta expresión se deben cumplir dos condiciones, la parte imaginaria y real de ambos lados de
la ecuación deben ser iguales, así:
F
(02  2 ) A  0 cos  (3.10)
m

7
 F0
    sen (3.11)
m
De la expresión (3.4) se obtiene el valor de la amplitud:
F0
A cos  (3.12)
m (02  2 )
Para que la amplitud A sea positiva, vemos que de la expresión (3.12) que  debe tomar los
valores:  = 0 en el rango     0 y  =  cuando   0 ; estas soluciones están
representadas en los dos gráficos de la figura (3.2); los valores de la fase inicial y de la amplitud
son idénticas a las obtenidas por el método convencional.
Ahora vamos a analizar la solución de la ecuación de movimiento, en el caso de que la frecuencia
de la fuerza externa sea igual a la frecuencia natural del sistema, vimos gráficamente que se
presenta una divergencia en el valor de la amplitud, ya que el valor de ésta, dada por las
expresiones (3.6), (3.7) y (3.12) se indetermina para   0 , para este valor de  la ecuación de
movimiento (3.1) dividida entre la masa, toma la forma:
d 2x F
2
 02 x  o cos  t (3.13)
dt m
La solución de la ecuación (3.13) será la suma de la solución de la ecuación homogénea y de la
particular o estacionaria. La solución a la ecuación homogénea tendrá la forma
xh (t )  A cos 0t  B sen0t , que puede ser escrita de forma idéntica a la ecuación (3.2). La
ecuación particular o estacionaria, no podrá tener la misma forma que la solución homogénea:
xe (t )  A cos 0t  B sen0t , pues se tendrían dos soluciones repetidas; de acuerdo con el Método
de Coeficientes Indeterminados (anexo al final del capítulo), la solución estacionaria estará
representada por:
xe (t )  t ( D cos 0 t  E sen0 t ) (3.14)
Sustituyendo la ecuación (3.14) y su segunda derivada en la expresión (3.13) se obtiene
F0
 2 0 D sen 0 t  20 E cos 0 t  02 t ( D cos 0 t  E sen0 t )  02 t ( D cos 0 t  E sen0 t )  cos 0 t
m
donde:
F0
E y D0
2 m 0
Así la solución particular o estacionaria será:
F0
t sen0 t
xe (t )  (3.15)
2 m 0
y la ecuación de movimiento cuando la frecuencia de la fuerza externa es igual a la frecuencia
natural del sistema estará dada por:

x (t )  C cos 0 t    
F0
t sen0 t (3.16)
2 m 0
Las ecuaciones (3.15) y (3.16) nos dicen que, cuando la fuerza externa tiene una frecuencia igual
a la frecuencia natural del sistema, la amplitud aumenta lineal con el tiempo.

8
 x(t)

Figura 3.4 representación de la solución de la solución estacionaria en resonancia, dada por la ecuación (3.15)
Este es un resultado extraordinariamente importante cuando se quiere usar o prevenir la
resonancia: aunque la amplitud de nuestro sistema se está acercando al infinito, ahora sabemos
que su dependencia con el tiempo es lineal. Después de dos oscilaciones, la amplitud se
duplicará, después de diez aumentará en un factor de diez, y así sucesivamente. En el caso de que
usemos resonancia para aumentar la amplitud de oscilación de un sistema, esto nos da una
estimación de la frecuencia con la que tenemos que "alimentarlo".

3.2 Oscilaciones forzadas con amortiguamiento

Consideremos un sistema masa-resorte como en la sección anterior, sobre el cual además de
actuar una fuerza impulsora, actúa una fuerza resistiva proporcional a la velocidad y una fuerza
elástica restauradora que tiende a llevar el bloque a su posición de equilibrio. El sistema se mueve
sobre una superficie horizontal y el desplazamiento del bloque se mide respecto de su posición de
equilibro, donde x = 0, sobre el bloque actúa una fuerza impulsora de la forma Fo cos  t , donde
 es la frecuencia angular de la fuerza externa; aplicando la segunda ley de Newton, obtenemos:
d 2x dx
m   kx  b  Fo cos  t
dt 2 dt
o bien
d 2x b d x k F
2
  x  o cos  t
dt m dt m m
llamemos
b k
  o2 
m m
2
d x dx F
2
 2  02 x  o cos  t (3.17)
dt dt m
La solución de la ecuación (3.17) se puede encontrar utilizando el método del exponente
complejo, donde:

z  A e j ( t   ) (3.18)
con x=Re (z), A representa la amplitud del movimiento que es siempre positiva y  es el ángulo
entre la fuerza impulsora y el desplazamiento, el signo menos es debido a que el desplazamiento
9
está atrasado con respecto a la fuerza impulsora; al igual que para un sistema forzado no
amortiguado, se espera que el valor de  sea cero para bajas frecuencias (donde la fuerza
impulsora tiene la misma dirección del desplazamiento y tiene un valor de  para altas
frecuencias (cuando las direcciones de la fuerza impulsora y el desplazamiento son contrarias);
reemplazando la expresión (3.18) en la ecuación (3.17), ésta toma la forma:
d 2z dz F
2
 2   02 z  0 e jt (3.19)
dt dt m
La solución de la ecuación (3.19) es la suma de dos soluciones, una solución homogénea xh(t) o
transitoria y una solución particular xp(t) o estacionaria, la solución homogénea se encuentra
cuando la ecuación (3.19) toma la forma:
d 2z dz
2
 2   02 z  0
dt dt
y como vimos en el capítulo 2, la ecuación diferencial anterior representa un sistema
amortiguado, cuyas soluciones pueden tener diferentes comportamientos dependiendo del signo
del radical que aparece al plantear la solución de la ecuación auxiliar, o sea el signo de  2   02
en m     2  02 ; así se puede tener que si:

 2  02  0 xh (t )  C e   t cos( t   )

 2  02 xh (t )  C1  C 2t  e   t

 2  02  0 x h (t )  C1 e m1 t  C 2 e m2 t
Estas soluciones corresponden a los movimientos llamados subamortiguado, críticamente
amortiguado y sobreamortiguado, respectivamente, es de notar que cuando el tiempo transcurrido
es muy grande xh(t) tiende a cero, dada su dependencia con la función exponencial decreciente,
de ahí que a esta solución homogénea de le conozca como solución transitoria.
Analicemos ahora la solución de la ecuación particular o estacionaria y calculemos los valores de
A y , para esto sustituimos la expresión dada en (3.18) en la ecuación (3.19), obteniéndose:

  2
A  2 j   A  02 A e j  t    
F0 j  t
e    2

 2 j     02 A e j  t e  j  
F0 j  t
m
e
m

  2

 2 j     02 A 
F0 j 
m
e

Por lo tanto,
F0 j
( 02   2 ) A  2 j   A  e (3.20)
m
Podemos considerar la ecuación (3.20) como una condición geométrica. El primer miembro nos
dice que dibujemos un vector de longitud (o2 - 2) A, y que luego perpendicularmente a él
dibujemos otro vector de longitud 2 A. El segundo miembro nos dice que dibujemos un vector
F
de longitud o formando un ángulo  con el eje real. La ecuación exige que estas dos
m
operaciones nos lleven al mismo punto, de modo que los vectores formen un triángulo cerrado,
como se ve en la figura 3.4(a).
10
 Figura 3.4 Representación geométrica de la ecuación (3.20)
Igualando los términos reales e imaginarios de la ecuación (3.20), se obtiene

   A  Fm cos 
2
o
2 o
(3.21)
Fo
   Asen | (3.22)
m
Elevando al cuadrado las ecuaciones (3.21) y (3.22) y sumando, se puede obtener la amplitud:
Fo
A()  m
     2  1/ 2
(3.23)
2 2 2 2
o

Dividiendo la ecuación (3.22) entre la ecuación ((3.21) se puede calcular el valor del ángulo .
2 
tg  ()  (3.24)
o2  2
En la figura 3.5 se puede observar la relación que existe entre la amplitud A y el ángulo de fase 
respecto a la frecuencia  en el caso de un valor constante de Fo. Estas curvas tienen un aspecto
general claro muy parecido al de la figura 3.2 para el caso del oscilador no amortiguado. En la
figura 3.5(a) vemos el comportamiento de la amplitud en función de la frecuencia de la fuerza
F
impulsora, cuando =0 el valor de la amplitud será o ; al aumentar , pero manteniéndose en el
k
rango  < r, dónde r es la frecuencia en la que se observa un máximo en la amplitud, el
denominador se hace cada vez más pequeño y el valor de la amplitud aumenta; cuando  sigue
aumentando pero considerando ahora el rango donde  > r el valor del denominador de la
ecuación (3.13) aumenta haciendo que la amplitud disminuya; al seguir aumentando el valor de la
frecuencia de la fuerza externa la amplitud tiende a cero. Cuando  = r (r < 0) el sistema está
en resonancia, donde se presenta la amplitud máxima; sin embargo, a diferencia de un sistema
forzado sin amortiguamiento ya no se observa una divergencia en la amplitud para esta
frecuencia de resonancia debido a que ahora tuvimos en cuenta la fuerza amortiguadora. A
continuación, se realiza un resumen del comportamiento de la amplitud para los diferentes
valores de la frecuencia de la fuerza impulsora:
F
=0 A()  o
k
 = m A()  Amax

11
 Fo
 = 0 A() 
2 m 0
→ A() 0
En la figura 3.5(b) vemos el comportamiento de la contante de fase  que varía suavemente entre
0 y , no aparece la discontinuidad observada en el caso del oscilador forzado sin
amortiguamiento, debido a la presencia de una fuerza amortiguadora; para entender el
comportamiento de la constante de fase analicemos la ecuación (3.24) en ella observamos que
cuando =0, el valor tan  = 0 y por tanto la constante de fase vale cero, cuando la frecuencia de
la fuerza impulsora aumenta el valor de la tan  también aumenta ya que el valor del
denominador se hace pequeño. Cuando  se aproxima a 0, el valor de tan  →  y en

consecuencia   , cuando  se hace mayor que 0 el valor de la tangente de la constante de
2
fase se hace negativa y cuando el valor de la frecuencia de la fuerza impulsora se hace muy
1
grande decimos que tan    y    ; el análisis anterior lo podemos resumir como:

=0 tan   0  0

 = 0 tan    
2
1
= tan     

Analizando los dos comportamientos, observamos que en la región donde  < 0,  es pequeño y
podemos decir que la fuerza impulsora está en fase con el desplazamiento (tienen la misma
 
dirección). Para  = 0 la constante de fase   , o sea que el desplazamiento está atrasado
2 2
con respecto a la fuerza impulsora y para  >> 0 se tiene que    , es decir, la fuerza
impulsora está desfasada un ángulo  con el desplazamiento (tienen direcciones opuestas).

A()




F0/k

0 m   0 
0

Figura 3.5. (a) Relación entre la amplitud y la frecuencia impulsora en el caso de oscilaciones forzadas con
amortiguamiento. (b) Fase de1 desplazamiento respecto a 1a fuerza impulsora en función de 1a frecuencia
impulsora.
Teniendo en cuenta la ecuación (3.23) donde vemos claramente que el valor de la amplitud
depende de , podemos calcular para que valor de la frecuencia de la fuerza impulsora el valor de
12
A() será máximo, esto se da, cuando su denominador:

y( )   2
o  2   4  
2
(3.25)
sea mínimo, lo que requiere que
dy
0
d
dy 1 2
d 2

 o  2   4  
2
 2
1/ 2
2
o  2  2   8  

dy


 4   o2   2  8  
0

d
 
1/ 2
 2 2
2 o    4 
2  
 
Para que la fracción anterior se anule, se requiere que su numerador sea igual a cero.
 4 o2   2   8    0   4 o2   2   2     0

 2
o   2   2   0

 2   o2  2   (3.26)
De esta última ecuación vemos que el valor de la frecuencia angular de la fuerza aplicada
ω 2 = ωr2 hace máxima la amplitud del desplazamiento y se le denomina frecuencia de
resonancia; de forma similar la frecuencia de resonancia en Hz, es:

r o2  2 
fr   (3.27)
 
Se obtiene el valor máximo de la amplitud A sustituyendo la expresión para  en la expresión
(3.23)
Fo Fo
A(r )  m  A(r )  m
 2
o  2   2
2 2
1/ 2
 2
o  o2  2 2 ) 2   2 o2  2 2 ) 
2

2 1/ 2

Fo
Amax  m (3.28)

2  o2   2  1/ 2

De las ecuaciones anteriores observamos que

 La frecuencia de resonancia difiere de la frecuencia natural o y de la frecuencia del oscilador
amortiguado.
 La frecuencia de resonancia disminuye cuando el factor de amortiguamiento crece.
 Para la luz láser ( pequeño) la curva de resonancia en aguda y alta, en cambio para la luz
natural ( grande) la curva es ancha y baja.
 Si b=0 la amplitud de resonancia es .

13
3.3. Potencia absorbida por un oscilador impulsado
Con frecuencia tendrá importancia e interés conocer a qué ritmo debe suministrarse energía a un
oscilador impulsado y amortiguado para mantener sus oscilaciones a una amplitud fija. Como en
cualquier otro fenómeno dinámico, podemos calcular la potencia instantánea P en función del
producto de la fuerza impulsora por la velocidad:
dW dx
P F Fv
dt dt
Tenemos que para el oscilador forzado con amortiguamiento, x(t) está dada por x  Re (Z ) , con
Z  A e j ( t  ) y por tanto x  A cos (t   )

 vt    A sen   t   
dx
dt
que se puede escribir en la forma
vt    vo sen  t   

en donde v0  A , es el valor máximo de v para valores cualesquiera de Fo y . Tomando el

valor de A de la ecuación (3.23) se tiene
Fo
v o ()  m 

 2  2
 o
2
 4   2  1 / 2


La expresión anterior se puede modificar, pasando la frecuencia como denominador del

denominador, así:
Fo
v o ()  m (3.29)
 
1/ 2
 2  2 2

 o 2  4  
  

El valor de vo pasa por un máximo para =o exactamente, fenómeno que podemos denominar
resonancia de velocidad, para este valor de frecuencia la expresión anterior se convierte en:
Fo
v o máx  (3.30)
2 m

Consideremos ahora el trabajo y la potencia necesaria para mantener las oscilaciones forzadas. Se
tiene que, la fuerza impulsora F= Fo cos t cuando hace que el sistema realice un desplazamiento
dx, realiza un trabajo
dW = Fo cos  t dx
En un ciclo completo el trabajo realizado será
2
dx
W  Fo cos t dx    Fo cos t
1 ciclo
0 dt
dt

14
 dx
Hallamos derivando la ecuación x  A cos (t   ) , con A dada por la ecuación (3.23), y
dt
reemplazamos en la ecuación anterior, obteniendo
2 2
W   Fo cos  t v(t ) dt  W   Fo cos  t   A sen ( t   ) dt
0 0

2
 F 
W    Fo cos t   o  sen t cos   sen cos t  dt
0
 my
2
 F 

W    Fo   o  cos  sen t cos t  sen cos 2 t dt (3.31)
0
 my
donde y está dada por (3.25); teniendo en cuenta el gráfico (3.6) que es una representación
geométrica para y, tenemos:
o2  2  
cos   y sen  
y y

Figura 3.6. Ayuda geométrica para definir las funciones trigonómetricas correspondientes al ángulo  que forma la
fuerza impulsora con el eje real
La ecuación (3.31) toma la forma

Fo2      2
2 
2

W 
my 

y 0

sen t cos t dt 
y 
0

cos 2 t dt 


Fo2        2
Fo2   2  
2
W 
m y 

y



0
 sen t cos t dt  
m y  y 
 0
 cos 2 t dt (3.32)

1
Para realizar la primera integral podemos usar la identidad sen  t cos  t  sen 2  t , así:
2

 sen 2 t dt  4   sen u du   4  cos 4    41 1 1  0

2 2
1 1 4 1 4
0

sen t cos t dt  0

0 0
2

Para realizar la segunda integral podemos usar la identidad cos 2  t 

1
1  cos 2  t  , así:
2

15
 2 2 2 2
 
1 cos 2 t  1  
cos 2 t 
0
 cos 2 t dt  
0 

2
 dt 


0

2
dt   0 

2
 dt

2
 
 
2
  1  2  1  2  
cos t dt   t   1 4cos   0  0 
1 1 4
 

u du    0  sen u 0  
2
0 2  4 0 2  4 2   
 0

Reemplazando, el valor anterior en la ecuación (3.32), se obtiene:
Fo2   2      2   Fo2 
W    
m y  y     m y2

Reemplazando el valor de y en la expresión anterior

2   Fo2  (3.33)
W

m ( 02   2 ) 2  4  2  2 
La expresión anterior representa la energía neta transmitida por la fuente en un ciclo, representa
además el trabajo hecho sobre la fuerza amortiguadora en un ciclo. La potencia está dada por
W W
P  ; la potencia media por ciclo será entonces P  y está dada por
t T
2   Fo2 

m (02   2 ) 2  4  2  2   P ( ) 
2   Fo2 
P ( ) 
2 
2  m (02   2 ) 2  4  2  2 

Fo2 

P ()  m (3.34)

o2  2  2

 4 
2
Esta potencia tiene su valor máximo cuando el primer término del denominador de (3.34) se hace
cero, o sea cuando =o que es exactamente la condición para la velocidad de resonancia. Así la
potencia de resonancia (máxima energía transferida) y la velocidad de resonancia ocurren a la
misma frecuencia; el valor de Pmáx está dado por

Fo2
Pmáx  (3.35)
4m
La potencia cae a la mitad de su valor cuando los términos del denominador de (3.34) son
iguales, en este caso, la Pmáx se puede escribir en la forma
Fo2 Q
Pmáx  (3.36)
2 m 

donde Q  ; volviendo a la expresión (3.34), ésta puede escribirse como
2
Fo2  1
P () 
m  2
o  2 
2

 4
2
16
 
Esta expresión puede escribirse en función de Q, donde 4  2 
Q2
Fo2    1
P (   
m 
 2 Q     2 
2



2 Q2
Fo2     1  1
P (      
 
m  2 Q   


    2 2


 

 
    Q 2 
Fo2 1
P (  (3.37)
2 m  Q     2  2  

    2
    Q 
La ecuación (3.37) se puede escribir en la forma
Fo2 1
P ( 
2 m  Q      
2

 
 
  2
    Q 
Es posible escribir una forma aproximada de la ecuación para P () basada en el siguiente
cálculo elemental:
     
 
     
    
       
De aquí que, si   o. podemos poner
  2     2   
  
   
Sustituyendo este valor en el denominador de la ecuación (3.26) se tiene
Fo2 1
P ( 
2 m  o Q  4  2 1 

  
   Q 2 
Multiplicando y dividiendo por  0
Fo2  1
P ( 
2 m  0 Q  4  2
2

  2
 02
Q 
F2  1
P (  o 

2 m Q 4  2  4    (3.38)

Las frecuencias o  , para las cuales P () disminuye a la mitad del valor máximo P (o),
dada por (3.38) se definen por
4  4  2  2   2
2

17
Es decir,
o
2   (3.39)
Q
Así pues, la anchura de la curva de resonancia para el oscilador impulsado, medida por la
potencia de entrada (figura 3.7), es igual al recíproco del tiempo necesario para que las
oscilaciones libres disminuyan a 1/e de su energía inicial.
Pmax
Potencia media de entrada

Anchura de la curva de resonancia

de potencia a la mitad de su altura
o muy proxima a 0 /Q
1/2 Pmax

0-1/2 0 0+1/2

Figura 3.7 “Agudeza” de la curva de resonancia determinada en función de la curva de potencia.

En la figura 3.8 se muestra la relación existente entre P y  para diversos valores de Q; vemos
que la potencia de entrada tiende a cero para frecuencias muy bajas y muy altas, y que excepto
para valores bajos de Q las curvas son casi simétricas alrededor del máximo. Es conveniente
definir una anchura para estas curvas de resonancia de potencia, tomando la diferencia entre
aquellos valores de  para los cuales la potencia de entrada es la mitad del valor máximo. Esto
puede hacerse de un modo particularmente claro y útil si Q es grande (como sucede en la mayoría
de los casos de interés). Esto significa que  de resonancia está efectivamente contenida dentro
de una banda estrecha de frecuencias próximas a o.
Potencia absorbida por un oscilador impulsado
P Q=50

Q=30

Q=20
Q=10
Q=3
Q=1
 

18
 Figura 3.8 Potencia media absorbida por un oscilador forzado en función de la frecuencia para diversos valore de Q

Podemos predecir que si un sistema tiene una respuesta de resonancia muy estrecha (medida bien
por la amplitud o por la absorción de potencia), la disminución de sus oscilaciones libres será
muy lenta e inversamente, como es natural, si las oscilaciones libres disminuyen rápida o
lentamente la respuesta del oscilador amortiguado es respectivamente ancha o estrecha. De las
ecuaciones (3.38) y (3.39) podemos decir que la resonancia es estrecha si la anchura es solamente
una pequeña fracción de la frecuencia de resonancia, es decir,
2 
 1 (3.40)


3.4 Analogías eléctricas

Los sistemas mecánicos pueden representarse y estudiarse por medio de sus circuitos eléctricos
equivalentes, los cuales se construyen más fácilmente que los modelos de los correspondientes
sistemas mecánicos; de aquí que sea más conveniente tomar los resultados experimentales de los
circuitos eléctricos equivalentes, que de los propios modelos mecánicos. Los circuitos eléctricos
equivalentes se obtienen comparando las ecuaciones de movimiento de los dos sistemas. Un
sistema mecánico y uno eléctrico son análogos (tabla 3.2) si sus ecuaciones diferenciales de
movimiento son matemáticamente iguales. Cuando esto ocurre, los términos correspondientes de
las dos ecuaciones diferenciales de movimiento también son análogos. Los circuitos eléctricos
equivalentes se pueden construir utilizando las leyes de Kirchhoff.

3.4.1 Leyes de Kirchhoff

 Ley de Kirchhoff de la tensión. En cualquier red, la suma algebraica de todas las tensiones
alrededor de un circuito cerrado cualquiera es igual a cero.
 Ley de Kirchhoff de la corriente. En cualquier circuito, la suma algebraica de las
intensidades de las corrientes alrededor de un nodo cualquiera es igual a cero.
Las aplicaciones técnicas más simples y difundidas de la resonancia están en la electricidad. En el
mundo eléctrico hay un gran número de dispositivos que se pueden conectar para construir
circuitos eléctricos. Entre el gran número de dispositivos eléctricos existentes estudiaremos tres
elementos que pueden hacer parte de un circuito. El primero se llama capacitor o condensador
(figura 3.9(a)); un ejemplo lo conforman dos placas planas metálicas separadas por un material
aislante una distancia muy pequeña. Cuando las placas se cargan se produce entre ellas una
diferencia de potencial; la misma diferencia de potencial aparece entre los terminales. Por lo
tanto, existe un cierto voltaje V entre las placas, si hay una cierta carga +q y -q en ellas. La
diferencia de potencial a través de un condensador es proporcional a la carga y está dada por:
q
VC  (3.41)
C
1
la constante de proporcionalidad es , donde C es la capacitancia del condensador, y su unidad
C
es el faradio. El segundo tipo de elemento de circuito se llama resistencia; y como su nombre lo
indica ofrece resistencia al paso de la corriente. Si hay una diferencia de potencial aplicado,

19
 dq
existe una corriente eléctrica i  que es proporcional a la diferencia del potencial eléctrico:
dt
dq
VR  R i  R (3.42)
dt
El coeficiente de proporcionalidad se llama resistencia R, y su unidad es el ohmio. Esta relación
es la ley de Ohm.
Si pensamos que la carga q en un capacitor es análoga al desplazamiento x de un sistema
dq 1
mecánico, vemos que la corriente i  es análoga a la velocidad, es análoga a la constante
dt C
k del resorte y R es análoga al coeficiente de la fuerza resistiva b. Es muy interesante que exista
además otro elemento de circuito que es análogo a la masa, es una bobina que forma un campo
magnético dentro de sí misma cuando hay una corriente circulando en ella. Un campo magnético
di
variable desarrolla en la bobina un voltaje que es proporcional a (de hecho, así funciona un
dt
transformador). El campo magnético es proporcional a la corriente y el voltaje inducido (así
llamado) en tal bobina es proporcional a la rapidez de variación de la corriente
di d 2q
VL  L L 2 (3.43)
dt dt
El coeficiente L es la inductancia y es análoga a la masa en un circuito mecánico oscilante. La
unidad de L es el henrio.

(a) capacitor (b) resistencia (c) inductancia

Figura 3.9 Los tres elementos pasivos de circuito.

Tabla 3.2 Analogía entre parámetros mecánicos y eléctricos

Propiedad mecánica Propiedad eléctrica variable
Tiempo (t) Tiempo (t)
Posición (x) Carga (q)
Masa (m) Inductancia (L)
Constante de fuerza viscosa (b) Resistencia (R)
Constante del resorte (k) Capacitancia (1/C)
Fuerza impulsora (F) Diferencia de Potencial Impulsor (V)
1
Frecuencia de resonancia ( k m ) Frecuencia de resonancia ( )
LC
Energía potencial (½ k x2 ) Energía de carga estática (½ q2/C)
Energía cinética (½ mv2) Energía electromagnética de la

20
 carga móvil (½ L C2 )
2
Potencia absorbida en la resonancia (F0 /2b) Potencia absorbida en la resonancia (V02/2R)

Uno de los sistemas resonantes más familiares e importantes es el sistema eléctrico compuesto
por un condensador y una bobina. El análisis de este sistema tiene gran semejanza con los
sistemas mecánicos vistos hasta ahora. Consideremos primero las oscilaciones libres ignorando
cualquier proceso disipativo.

3.4.2 Circuito LC
El circuito oscilante más simple está compuesto por un capacitor de capacitancia C conectado en
serie a un inductor de inductancia L, como se muestra en la figura 3.10. Supondremos que el
capacitor está inicialmente cargado, cuando se cierra el circuito una corriente i circula por este.
La diferencia de potencial en el capacitor y en el inductor están dadas por las expresiones (3.41) y
dq d i d 2q
(3.43). Recordemos que i = , por tanto, = 2 . Podemos analizar este circuito usando la
dt dt dt
primera ley de Kirchhoff, en este caso, la suma de todas las diferencias de potencial alrededor del
circuito cerrado debe ser cero, por lo que puede escribirse:

0  VL  VC
di q
0L  (3.44)
dt C
Se puede obtener una ecuación que sólo contenga la carga, resulta entonces:
d 2q 1
0L  q (3.45)
dt 2 C
Esta ecuación representa el cambio de la carga eléctrica en una de las placas del capacitor, note
que esta expresión es totalmente semejante a la ecuación diferencial básica del movimiento
1
armónico simple para un sistema masa-resorte, donde q desempeña el papel de x, L el de m y
C
el de k; de lo anterior podemos admitir la existencia de oscilaciones eléctricas libres tales que

1
0 
LC

Figura 3.10. Circuito LC

21
La solución a la ecuación (3.44) es:
q (t ) = q0 cosωt +  
La energía almacenada en el capacitor es
1 1 q2
EC  C VC2 
2 2 C
y la energía en el inductor es:
2
1 1  dq 
EL = L i 2 = L  
2 2  dt 
Sumando estas dos energías y derivando la energía total con respecto al tiempo, llegamos a la
ecuación diferencial de movimiento (3.44).

3.4.3 Circuito RLC

Ahora se considera un circuito como el representado en la figura 3.11, en el que se ha agregado
una resistencia R (figura 3.9(b)). Una resistencia es un elemento del circuito que siempre disipa
energía eléctrica como calor, por lo que es muy distinto de un capacitor o una bobina, que son
elementos que pueden almacenar energía.

Figura 3.11 Circuito RLC

La presencia de la resistencia origina un término adicional en la ecuación (3.44) que expresa una
diferencia de potencial VR=iR entre sus terminales. Por lo tanto, aplicando la ley de Kirchoff:
di q
0L   iR (3.46)
dt C
ó
d 2q dq 1
L 2
R  q0 (3.47)
dt dt C
Como se ve esta ecuación tiene la forma de un movimiento oscilatorio amortiguado que resulta
cuando se somete un oscilador armónico simple a una fuerza disipativa que es directamente
proporcional a la velocidad.
d 2 q R dq 1
2
  q0 (3.48)
dt L dt L C
Comparando con
22
 d2 x dx
2
 2  o2 x  0
dt dt
Comparando estas últimas dos ecuaciones diferenciales de movimiento, se llega a las relaciones
presentadas en la tabla 3.2.
Para escribir la solución de la ecuación (3.48), es cuestión de remplazar todos los símbolos en la
ecuación general del movimiento del movimiento amortiguado, obteniéndose así:

q (t )  e  t  A e 
 2 o2 t   2 o2 t
 Be (3.49)
 
en el caso en que  2   o2 < 0, o sea  2 < o2 , las raíces de la ecuación diferencial son diferentes
complejos conjugados y la solución general de la ecuación (3.48) corresponde a un sistema
subamortiguado:

q(t )  e  t (c1 cos o2   2 t  c2 sen o2   2 t )

donde c1 y c2 son constantes arbitrarias. Esta ecuación se puede escribir de la siguiente forma


q(t )  c e  t cos o2   2 t   

R  R2 
cos  o2  2 t   
t
q(t )  c e 2L
(3.50)
 4L 
 
Estas soluciones solo son válidas para resistencias tan pequeñas que:
L
R2
C
En otras palabras
R2 1
2

4L LC
En tanto que satisfaga esta desigualdad, la corriente (derivada de la ecuación (3.50) tendrá la
forma oscilatoria amortiguada. En el caso en que  2   o2 = 0, o sea  2 = o2 , la resistencia es
R2 1
tan grande que 2
 , entonces, la solución toma la forma:
4L LC
R

q t   c1  c2 t  e
t
2L

Ahora la corriente pierde totalmente su carácter oscilatorio y decrece en forma exponencial, en

estas condiciones se dice que el circuito está críticamente amortiguado. Si se aumenta todavía
más la resistencia, la solución debe escribirse como una suma de dos términos que decaen
exponencialmente.
   2  02 t    2  02 t
q(t )  c1 e c 2 e

Se dice el ahora que el circuito está sobreamortiguado. Desde el punto de vista de la

conservación de la energía, la disminución en la amplitud de la corriente asociada a los circuitos
oscilatorios amortiguados se debe a la conversión continua en calor de la energía almacenada en
los campos eléctricos y magnéticos del capacitor y del inductor. En cualquier intervalo de tiempo
23
dt, la cantidad de energía eléctrica disipada como calor seria I2R dt, y de acuerdo con la ley de
conservación de la energía, la suma de las energías de los campos eléctricos y magnéticos debe
disminuir en el mismo intervalo exactamente en esa cantidad.

3.4.4 Circuito RLC en serie con una FEM alterna

En la sección anterior analizamos las soluciones transitorias de un circuito RLC, la energía
almacenada en la inductancia y en el capacitor se disipa finalmente en esos casos, y por último la
corriente cae a cero. Para mantener una corriente oscilatoria permanente en un circuito RLC, debe
tenerse una fuente que también sea oscilatoria, por ejemplo, un generador de corriente alterna
puede proporcionar una FEM oscilatoria de la forma
V(t)= V0 cos  t (3.51)
De acuerdo con las leyes de Kirchhoff, este circuito RLC (figura 3.12) satisface la ecuación:
d 2q dq 1
L 2
R  q  V0 cos  t (3.52)
dt dt C
d 2 q R dq 1 V
2
  q  0 cos t (3.53)
dt L dt LC L

Figura 3.12. Circuito RLC en serie con una FEM

La ecuación (3.53) puede ser comparada con la ecuación de movimiento correspondiente a un

sistema masa-resorte donde existe una fuerza impulsora y una fuerza de amortiguamiento:
d 2x b d x k F
2
  x  o cos  t
dt m dt m m
Del mismo modo que Fdx es el trabajo realizado por la fuerza impulsora F en un desplazamiento
dx, Vdq es el trabajo realizado por la tensión impulsora V cuando pasa a través del circuito una
cantidad de carga dq.

3.4.5 Impedancia en el oscilador forzado y amortiguado

Un oscilador amortiguado se caracteriza por tres cantidades: su masa m, la constante elástica k, y
la constante de amortiguamiento b. Ahora recordemos la ecuación (3.18)
24
 Fo
m Fo
vo ()   vo () 
   
1/ 2 1/ 2
 2  2 2
  2 2  2 2
 b 
 o 2  4   m o
 m 2  2 
    2 m 

Fo Fo
vo ()  1/ 2
 vo ()  1/ 2
m2  k 2
2

   m  k  m  2  2 
 2    b           b  
   m       m     

Fo
vo ()  1/ 2
(3.54)
 k
2


 m     b 
  
La cantidad que aparece en el denominador de la ecuación (3.54), se denomina impedancia del
oscilador, y se designa por la letra Z:
1/ 2
 k
2

Z m   m     b   (3.55)
  
Se definen la reactancia X y la resistencia R, por:
k
X  m   R b (3.56)

De esta forma la ecuación (3.55) se puede escribir como:


Zm  X 2  R 
1/ 2
(3.57)
Haciendo la analogía entre los sistemas mecánicos y eléctricos tenemos que de la ecuación
(3.55), se puede obtener la impedancia asociada al sistema eléctrico:
1/ 2
 1 
2

Z e    L    R   (3.58)
 C  

Problemas
N s N
1. Considere un oscilador amortiguado con m = 0,2 kg, b = 4 y k= 80 ; suponga que este
m m
oscilador se ve impulsado por una fuerza F = F0 cos t, siendo F0=2 N y  = 30 s-1; (i) ¿cuáles
son los valores A y  de la respuesta de estado estacionario descrita por x = A cos (t - )?; (ii)
¿cuánta energía se ve disipada contra la fuerza resistente en un ciclo?; (iii) ¿cuál es la potencia
media de entrada? Respuesta: (i) 1,3 cm, 130º; (ii) 0,063 J, (iii) 0,30 W
2. Un objeto de masa 2 kg cuelga de un resorte de masa despreciable, el resorte se alarga en 2,5
cm cuando se le sujeta dicho objeto; el extremo superior del resorte se hace oscilar hacia arriba y
hacia abajo con un movimiento armónico simple de amplitud 1 mm; si el valor de Q del sistema
es 15; (i) ¿cuál es o para este sistema?; (ii) ¿cuál es la amplitud de la oscilación forzada para  =

25
o; (iii) ¿cuál es la potencia media de entrada para mantener la oscilación forzada a una
frecuencia 2% mayor que o ? Respuesta: (i) 19,8 s-1; (ii) 1,5 cm; (iii) 0,086 W
3. El gráfico P3-1 muestra la curva de resonancia de potencia de un determinado sistema
mecánico cuando se ve accionado por una fuerza F=F0 sen t, donde F0 = constante y  es
variable; (i) hallar los valores numéricos de o y Q para este sistema, (ii) si se suprime la fuerza
1
impulsora, después de cuántos ciclos de oscilación libre ha descendido la energía a 5 de su
e
valor inicial? (con buena aproximación, el período de la oscilación libre puede hacerse igual a
2/o).

Figura P3-1
4. Se observa que las oscilaciones libres de un sistema mecánico tienen una cierta frecuencia
angular 1. El mismo sistema, cuando se ve accionado por una fuerza F=F0 cos t (donde F0 =
constante y  es variable), tienen una curva de resonancia de potencia cuya anchura de la
1
frecuencia angular a un valor de potencia igual a la mitad del valor máximo es ; (i) ¿para qué
5
frecuencia angular se produce la máxima potencia de entrada?, (ii) ¿cuál es el valor Q del
sistema?; (iii) si el sistema se compone de una masa m unido a un resorte de constante k; en
función de m y k, ¿cuál es el valor de la constante b?
5. La figura P3-2 muestra la potencia media de entrada P en función de la fuerza impulsora en el
caso de una masa unida a un resorte, donde se presenta amortiguamiento; la fuerza impulsora
tiene la forma F=F0 sent, donde F0=constante y  es variable. El valor de Q es suficientemente
elevado para que la potencia media de entrada, que es máxima para o, disminuya hasta la mitad
del máximo para las frecuencias 0,98 0 y 1,02 0 ; (i) ¿cuál es el valor numérico de Q?; (ii) Si se
suprime la fuerza impulsora, la energía disminuye de acuerdo con la expresión E=E 0e-t; ¿cuál es
el valor de ?; (iii) si se suprime la fuerza impulsora, ¿qué fracción de la energía se pierde por
ciclo?

26
 Figura P3-2

6. Un peso de 6 libras está unido al extremo inferior de un resorte que está suspendido del techo,
libra
la constante del resorte es 27 . El peso se encuentra en reposo en su posición de equilibrio
pie
y, al empezar en t = 0 se aplica al sistema una fuerza externa dada por F(t) = 12 cos 20t.
Determine el desplazamiento resultante como una función del tiempo, suponiendo que el
cos 12 t  cos 20 t
amortiguamiento es despreciable. Respuesta: x 
4
7. Un peso de 16 libras está unido al extremo inferior de un resorte que está suspendido del techo;
el peso se encuentra en reposo en su posición de equilibrio y el resorte está alargado 0,4 pies;
después, al empezar en t = 0, se aplica al sistema una fuerza externa dada por F(t) = 40 cos 16t.
dx dx
El medio ofrece una resistencia en libras numéricamente igual a 4 , donde es la
dt dt
pies
velocidad instantánea en . Determine el desplazamiento del peso como una función del
s
tiempo.
8. Un peso de 10 libras cuelga del extremo inferior de un resorte que está suspendido del techo, la
libras
constante del resorte es 20 ; el peso se encuentra en reposo en su posición de equilibrio y,
pie
al empezar en t = 0 se aplica al sistema una fuerza externa dada por F(t) =10 cos 8t. El medio
dx dx
ofrece una resistencia en libras numéricamente igual a 5 donde es la velocidad
dt dt
pies
instantánea en . Determine el desplazamiento del peso como una función del tiempo.
s
sen 8 t
Respuesta: x (t )   2 t e 8t 
4
9. Un peso de 4 libras cuelga del extremo inferior de un resorte que está suspendido de una viga,
el peso se encuentra en reposo en su posición de equilibrio y el resorte está alargado 3 pulgadas.
Después, el peso se desplaza 6 pulgadas hacia abajo de esta posición y se suelta en t = 0, en este
instante, se aplica al sistema una fuerza externa dada por F(t) =13 sen 4t. La resistencia del

27
 pies
medio en libras es numéricamente igual a dos veces la velocidad instantánea medida en; (i)
s
determine el desplazamiento del peso como una función del tiempo, (ii) observe que el
desplazamiento es la suma de un término transitorio y un término de estado permanente, y
determine la amplitud del término de estado permanente.
10. Un peso de 6 libras cuelga del extremo inferior de un resorte que está suspendido del techo, el
peso se encuentra en reposo en su posición de equilibrio y el resorte está alargado 4 pulgadas.
Después, al empezar en t = 0 se aplica al sistema una fuerza externa dada por F(t) =27 sen 4t -3
cos 4t. Si el medio ofrece una resistencia en libras numéricamente igual a tres veces la velocidad
pies
instantánea, medida en , determine el desplazamiento en función del tiempo. Respuesta:
s
 2 
x (t )  e  8t  sen 4 2 t  cos 4 2 t   sen 4 t  cos 4 t
 2 
libras
11. Un resorte que tiene una constante de 20 está suspendido del techo; un peso de 32
pie
libras está unido al extremo inferior del resorte que se encuentra en reposo en su posición de
equilibrio. Al empezar en t = 0 se aplica al sistema una fuerza externa dada por F(t) = 40 cos 2t;
después, el efecto de esta fuerza se mantiene hasta t = , y en ese instante cesa de aplicarse,
cuando t > , no actúa ninguna fuerza externa. El medio ofrece una resistencia en libras
dx dx pies
numéricamente igual a 4 , donde es la velocidad instantánea en . Determine el
dt dt s
desplazamiento del peso como una función del tiempo para todo t  0. Respuesta:
 3 sen 4t 
x  e 2 t    2 cos 4 t   sen 2 t  2 cos 2 t , 0  t   ;
 2 
   3 sen 4t
x  e 2  1 e 2 t 

 2 cos 4 t  , t  
 2 
12. Un peso de 12 libras está unido al extremo inferior de un resorte que está suspendido del
techo, el peso se encuentra en reposo en su posición de equilibrio y el resorte está alargado 6
pulgadas. Al empezar en t = 0 se aplica al sistema una fuerza externa dada por F(t) = 2 cos t; (i)
dx dx
si la fuerza de amortiguamiento en libras es numéricamente igual a 3 , donde es la
dt dt
pies
velocidad instantánea en , (i) determine la frecuencia de resonancia del movimiento
s
resultante y determine el desplazamiento como una función del tiempo, cuando la fuerza externa
está en resonancia con el sistema; (ii) suponiendo que no existe amortiguamiento, determine el
valor de  que da lugar a una resonancia no amortiguada, y encuentre en este caso la posición
como una función del tiempo.
13. Un peso de 20 libras está unido al extremo inferior de un resorte que está suspendido del
techo, el peso se encuentra en reposo en su posición de equilibrio y el resorte está estirado 6
pulgadas. Se aplican al sistema varias fuerzas externas de la forma F(t) = cos t y se encuentra
ciclos
que la frecuencia de resonancia es 0.5 . Suponiendo que la resistencia del medio en libras
s
28
 dx dx
es numéricamente igual a a , donde es la velocidad instantánea en pies por segundo,
dt dt
determine el coeficiente de amortiguamiento a.
14. La ecuación diferencial para el movimiento de una masa unitaria en un cierto resorte bajo la
acción de una fuerza externa de la forma F(t)=30 cos t es
d2 x dx
2
a  24 x  30 cos  t
dt dt

donde a  0 es el coeficiente de amortiguamiento; (i) Si a = 4, determine la frecuencia de

resonancia y determine también la amplitud de la vibración de estado permanente, cuando la
función de fuerza está en resonancia con el sistema; (ii) proceda como en la parte (i) si a = 2.
15. Un circuito tiene en serie una fuerza electromotriz dada por E(t) =100 sen200t V, un
resistor de 40 , un inductor de 0.25 H y un capacitor de 410-4 faradios. Si la corriente inicial es
cero y la carga inicial en el capacitor es 0.01 coulombs; calcule la corriente en cualquier tiempo t
> 0.
16. Un circuito tiene en serie una fuerza electromotriz dada por E(t) = 5 sen 100 t V, una
resistencia de 40 , un inductor de 0.05 H y un capacitor de 210-4 faradios. Si la corriente
inicial y la carga inicial en el capacitor son ambas cero, calcule la carga en el capacitor en
cualquier tiempo t > 0.
17. Un circuito tiene en serie una fuerza electromotriz dada por E(t) = 200e-100t V, un resistor
de 80 , un inductor de 0.2 H y un capacitor de 510-6 farads. Si la corriente inicial y la carga en
el capacitor son cero, calcule la corriente en cualquier tiempo t > 0.
18. La potencia promedio absorbida por un oscilador forzado es igual a la potencia disipada por
la fuerza amortiguadora. Obtenga la expresión (3.34) considerando que la potencia instantánea
está dada por P = b v2

Anexo 1: Método de los coeficientes indeterminados

Consideremos el caso especial de una ecuación diferencial lineal no homogénea de n-ésimo orden
en la que todos los coeficientes son constantes reales, es decir, analizaremos la ecuación:
dn y dn1 y dy
an n
+ a n 1 n1
+ ... + a 1 + a0 y = g (x) (3A-1)
dx dx dx
La solución de esta ecuación diferencial está dada por la suma de dos términos, uno
correspondiente a la solución de la ecuación homogénea yh(x) y otro llamado solución particular
yp(x); luego la solución general de (3A-1) es:
𝑦(x) = yh(x) + yp(x)
En el capítulo 2, vimos cuál es el procedimiento para calcular yh(x) ahora veremos un método que
nos ayuda a obtener yp(x) y es el método de los coeficientes indeterminados.
La idea básica de éste método es una propuesta coherente acerca de la forma de yp(x) originada
por los tipos de funciones que forman la función de entrada g(x). El método es básicamente
directo, pero está limitado a ecuaciones lineales no homogéneas como la ecuación (3A-1), donde:
29
 los coeficientes 𝑎𝑖, para 𝑖 = 0, 1, 2…𝑛; son constantes
 g(𝑥) puede ser una constante 𝑘, una función polinomial de cualquier grado, una función
exponencial 𝑒𝑎𝑥, una función seno o coseno, como seno(𝑏𝑥) o coseno(𝑏𝑥) o sumas y productos
finitos de estas funciones.
En la tabla 3A-1 se muestran algunos ejemplos específicos de g(x) con la forma correspondiente
de su solución particular yp(x); naturalmente suponemos que ninguna función en la solución
particular propuesta yp(x) esta duplicada por una función en la solución homogénea yh(x).
Veamos dos casos que se pueden presentar al aplicar este método propuesto.
Caso I: ninguna función en la solución particular propuesta yp(x) es solución de la ecuación
diferencial homogénea asociada yh(x).
Caso II: una función en la solución particular propuesta yp(x) también es solución de la ecuación
diferencial homogénea asociada yh(x).
Tabla 3A-1 listado de funciones correspondientes a la solución particular, yp(x)
g(x) yp(x)
Constante A
3x7 AxB
4 x2  3 A x2  B x  C
x3  x  2 A x3  B x 2  C x  D
3 sen 4 x A sen 4 x  B cos 4 x
cos 4 x A sen 4 x  B cos 4 x
e5 x A e5 x
(9 x  2) e5 x ( A x  B) e 5 x
x 2 e5 x ( A x 2  B x  C ) e5 x
e3 x sen 4 x A e3 x sen 4 x  B e3 x cos 4 x
5 x 2 cos 4 x ( A x 2  B x  C ) sen 4 x  ( D x 2  E x  F ) cos 4 x
x e3 x sen 4 x ( A x  B) e3 x sen 4 x  (C x  D) e3 x cos 4 x

Ejemplo 1: Resuelva la ecuación diferencial no homogénea:

d2y dy
2
+4  2 y = 2 x2  3 x  6 (3A-2)
dx dx
Solución: primero establecemos la solución de la ecuación homogénea yh(x), es decir, aquella que
satisface la ecuación homogénea
d2y dy
2
+4 2y = 0
dx dx
para la cual establecemos su ecuación auxiliar: 𝑚2+4𝑚−2=0, y la resolvemos utilizando la
fórmula general:

30
 m
1
2

 4  16  4(2)  m  2  2 6

las raíces son reales y diferentes por lo que la solución de la ecuación homogénea tendrá la
forma:
yh ( x)  Ae ( 2  2 6)x
 Be ( 2  2 6) x

yh ( x)  e2 x ( Ae 2 6x
 Be  2 6x
)

Ahora como g ( x) = 2 x 2  3 x  6 es un polinomio de orden 2, entonces de acuerdo a la tabla

3A-1 proponemos que su solución particular también tenga la forma de una ecuación cuadrática:
y p ( x)  A x 2  B x  C
donde el objetivo es encontrar los valores de las constantes 𝐴, 𝐵 𝑦 𝐶 que satisfacen (3A-2). Aquí
observemos que ningún término de la solución particular propuesta 𝑦𝑝(x) está presente en 𝑦h(x),
por lo que tenemos un ejemplo del caso I. Para lograr el objetivo mencionado, calculamos la
primera y segunda derivada de 𝑦𝑝(x), es decir:
yp ( x)  2 A x  B yp ( x)  2 A

Reemplazando estos valores en (3A-2), tenemos

2 A  4 (2 A x  B)  2( A x 2  B x  C )  2 x 2  3x  6
Agrupando se tiene:
 2 A x 2  (8 A  2 B) x  (2 A  4 B  2C )  2 x 2  3x  6
Para que esta última igualdad se cumpla, los coeficientes de las diferentes potencias de x a ambos
lados de la ecuación deben ser iguales, es decir:
 2A 2  A  1
5
8 A  2B   3  B 
2
2 A  4 B  2C  6  C  9
Por lo que la solución particular es:
5
y p ( x)   x 2  x 9
2
y por lo tanto, la solución general es:
5
y ( x)  e 2 x ( Ae 2 6x
 Be  2 6x
)  x2  x 9
2

Ejemplo 2. Resuelva la ecuación diferencial no homogénea:

d2y dy
2
5  4 y = 8 ex (3A-3)
dx dx

31
Solución: primero establecemos la solución de la ecuación homogénea yh(x), es decir, aquella que
satisface la ecuación homogénea
d2y dy
2
5  4y = 0
dx dx
para la cual establecemos su ecuación auxiliar: 𝑚2-5𝑚+4=0, y la resolvemos utilizando la
fórmula general:

m
1
2

5  25  4(4)  1
m  (5  3)
2
las raíces son reales y diferentes por lo que la solución de la ecuación homogénea tendrá la
forma:
yh ( x)  Ae 4 x  Be x

Para hallar la solución particular, tenemos que el término no homogéneo g(x) = 8 ex es una
constante multiplicada por una función exponencial, entonces de acuerdo a la tabla 3A-1
proponemos que su solución particular también tenga la forma de una constante multiplicada por
una exponencial
y p ( x)  A e x

donde en este caso, el objetivo es encontrar el valor de la constante 𝐴 que satisface la ecuación
x
(3A-3). Aquí observemos que Ae de la solución de la ecuación particular propuesta yp(x), es
idéntica al segundo término de la ecuación homogénea yh(x), por lo que tenemos un ejemplo del
caso II. Ninguna función en la solución particular propuesta yp(x) debe estar contenida en la
función de la solución homogénea yh(x), por lo que nuestra solución particular propuesta yp(x) es
incorrecta. Para corregirla aplicaremos la siguiente regla:
Regla de multiplicación para el caso II. Sí algún término de la solución particular propuesta
yp(x) está contenido en la ecuación homogénea yh(x), entonces dicho término se debe multiplicar
por 𝑥𝑛, donde 𝑛 es el entero positivo más pequeño que elimina esa duplicación.
x
Entonces, de acuerdo a la regla anterior, es suficiente multiplicar Ae por 𝑥 para que sea diferente
de Be x . Por lo que la solución particular propuesta correcta es:
y p ( x)  A x e x

yp ( x)  A e x  A x e x yp ( x)  2 A e x  A x e x
Sustituyendo en la ecuación (3A-3), tenemos:
2 Aex  A x ex  5 ( Aex  A x ex )  4 A x ex  8ex
8
3 A e x  8 e x A
3
Por lo tanto, la ecuación particular será:
8
y p ( x)   x e x
3

32
y la solución general:
8
y ( x)  A e 4 x  B e x  x e x
3
Referencias
 French, A. P. Vibraciones y ondas, Editorial Reverté. S.A., 1974
 Hecht, E. Física en perspectiva, Editorial Addison-Wesley Iberoamericana, 1987
 King, G. C. Vibrations and waves, Wiley, 2009
 Ross, S. L. Introducción a las Ecuaciones Diferenciales, Editorial McGraw-Hill, 1989
 Seto, W. Vibraciones Mecánicas, Editorial McGraw-Hill, 1970
 Thomson, W. T. Teoría de vibraciones. Aplicaciones, Editorial Prentice-Hall
Hispanoamericana, S.A., 1982

33