Ver el código de «Oscilaciones amortiguadas y forzadas (CMR)»

==Oscilaciones amortiguadas==
===El oscilador no amortiguado===
En otras secciones se estudia la [[Cinemática_del_movimiento_rectilíneo_(GIE)#Arm.C3.B3nico_simple|cinemática]] y la [[Aplicaciones_de_las_leyes_de_Newton_(GIE)#Oscilador_arm.C3.B3nico|dinámica]] del oscilador armónico. Éste es un sistema ideal gobernado por la ley de Hooke. Típicamente esta ley se aplica a resortes mecánicos, aunque puede generalizarse a muchas otras situaciones. En el caso de un resorte que oscila en una sola dimensión la ley de Hooke se escribe

<center><math>F = -kx\,</math></center>

siendo <math>x</math> la ''elongación'' del resorte (distancia respecto a la posición de equilibrio)

Una partícula sometida exclusivamente a la ley de Hooke en una dimensión cumple la ecuación de movimiento

<center><math>ma = -kx\qquad\rightarrow\qquad a= \ddot{x}=-\frac{k}{m}x</math></center>

Este es un caso particular de la ecuación para un movimiento armónico simple

<center><math>\ddot = -\omega_0^2x</math></center>

siendo en este caso la frecuencia natural 

<center><math>\omega_0=\sqrt{\frac{k}{m}}</math></center>

La solución general de esta ecuación diferencial es una oscilación sinusoidal

<center><math>x(t) = A\cos(\omega_0 t+\varphi)\,</math></center>

<center>[[Archivo:Muelle.gif]]</center>

con <math>A</math> la amplitud de las oscilaciones, <math>\varphi</math> la fase inicial o constante de fase. Este movimiento es periódico, de forma que

<center><math>x(t+T) = x(t)\qquad\qquad T = \frac{2\pi}{\omega}</math></center>

Esta solución también se puede escribir como una combinación lineal de un seno y un coseno

<center><math>x(t) = b_1\cos(\omega_0 t)+b_2\,\mathrm{sen}(\omega_0 t)</math></center>

con

<center><math>b_1=A\cos(\varphi) \qquad\qquad b_2 = -A\,\mathrm{sen}(\varphi)</math></center>

Los valores de las constantes <math>b_1</math> y <math>b_2</math> pueden calcularse también a partir de las condiciones iniciales del movimiento

<center><math>b_1= x_0 = x(t=0)\qquad\qquad b_2=\frac{v_0}{\omega_0}=\frac{v(t=0)}{\omega_0}</math></center>

===Amortiguamiento===
El modelo de un oscilador mecánico sometido exclusivamente a la ley de Hooke no es realista pues desprecia la presencia del rozamiento. La experiencia nos muestra que un oscilador se va frenando progresivamente hasta llegar a detenerse en la posición de equilibrio.

Esta disminución  progresiva en la amplitud de las oscilaciones es debida a la presencia de rozamiento. Éste puede deberse a un roce con una superficie (rozamiento seco) o la fricción del aire o líquido que rodea al oscilador (rozamiento viscoso).

El caso del oscilador con rozamiento seco tiene un interesante análisis físico-matenático, pero no lo consideraremos aquí, sino en un [[Resorte_con_rozamiento_seco|problema]]. En su lugar nos centraremos en el caso del rozamiento viscoso. La razón es que, aparte de ser un modelo de muchas aplicaciones, representa más adecuadamente lo que ocurre en un amortiguador mecánico.

Un amortiguador es un dispositivo como el que puede encontrarse en la suspensión de un automóvil o en una puerta con cierre automático.

<center>[[Archivo:amortiguador.jpg|400px]]</center>

Un amortiguador consta de un resorte mecánico, pero también, en el interior de éste, de un cilindro con un pistón

<center>[[Archivo:amortiguador-dentro.gif]]</center>

Si un coche no tuviera suspensión (es decir, si el chasis estuviera unido rígidamente al eje de las ruedas), cada bache o irregularidad en el suelo se notaría como un golpe en el interior del vehículo lo cual, además de incómodo, pone en peligro su integridad. Por otro lado, si la suspensión consistiera simplemente en un resorte casi sin rozamiento, cada bache produciría oscilaciones en el coche, incluso mucho después de haber superado el bache.

<center>[[Archivo:bache.gif]]</center>

Por ello, se introduce el amortiguador. El objetivo es que el coche oscile al pasar por el bache, pero lo menos posible, de forma que retorne a la posición de equilibro en el menor tiempo posible. Esto se consigue introduciendo una fricción viscosa que disipe la energía mecánica de la oscilación. En la práctica consiste en que un fluido es obligado a pasar por una serie de válvulas de un lado a otro del pistón, frenándolo en el proceso.

La fuerza de rozamiento que experimenta el resorte se opone siempre a la velocidad de éste (si la masa va hacia la derecha, la fuerza apunta hacia la izquierda y viceversa). En primera aproximación es proporcional a la velocidad (en reposo no hay fuerza de rozamiento), por lo que se puede escribir

<center><math>\vec{F}=-\gamma\vec{v}</math></center>

y para el caso particular del movimiento rectilíneo

<center><math>F = -\gamma v = -\gamma\dot{x}</math></center>

===Ecuación del oscilador amortiguado===
La segunda ley de Newton para un oscilador armónico con amortiguamiento viscoso (en una dimensión) se escribe entonces

<center><math>ma = -kx-\gamma v\,</math></center>

Pasando todo al primer miembro

<center><math>ma + \gamma v + k x = 0\,</math></center>

Aplicando que la velocidad y la aceleración son las primera y segunda derivadas respecto al tiempo de la elongación nos queda la ecuación diferencial

<center><math>m \ddot{x} + \gamma \dot{x}+k x = 0</math></center>

Dividiendo por la masa de la partícula podemos escribirla como

<center><math>\ddot{x}+2\beta \dot{x} +\omega_0^2 x = 0</math></center>

Esta es la ecuación diferencial del oscilador armónico amortiguado. La constante 

<center><math>\omega_0 = \sqrt{\frac{k}{m}}</math></center>

es la ''frecuencia propia'' del oscilador. Equivale a la frecuencia natural con la que oscilaría el resorte si no tuviera rozamiento. Como veremos, la presencia de rozamiento reduce la frecuencia de las oscilaciones.

La segunda constante

<center><math>\beta = \frac{\gamma}{2m}</math></center>

es la ''constante de amortiguamiento''. Mide la magnitud de la fricción, siendo mayor cuanto más intensa sea ésta.

Tanto la frecuencia propia <math>\omega_0</math> como la constante de amortiguamiento <math>\beta</math> tienen dimensiones de inversa de un tiempo y se miden en s<sup>&minus;1</sup> en el SI.

===Solución de la ecuación===
====Caracterización de las soluciones====
Antes de examinar la solución matemática de la ecuación diferencial, podemos describir como deberían ser las soluciones.

* Si el rozamiento es pequeño, debemos esperar que el resorte oscile pero con una amplitud decreciente, hasta que pasado un cierto tiempo se quede en reposo en la posición de equilibrio.

* Si el rozamiento es muy grande, en cambio, esperamos que no llegue a oscilar, sino que simplemente se va moviendo lentamente hasta la posición de equilibrio.

En física siempre que una magnitud se considera grande o pequeña hay que decir comparada con qué, cuál es el patrón en que nos basamos para decir si algo es grande o pequeño. En este caso aprovechamos que tanto <math>\beta</math> como <math>\omega_0</math> tienen las mismas dimensiones y por tanto se pueden comparar. Establecemos entonces el criterio

* Rozamiento débil: <math>\beta < \omega_0</math>
* Rozamiento intenso: <math>\beta > \omega_0</math>

La solución matemática debe reflejar por tanto un cambio de comportamiento dependiendo de como sea la constante de amortiguamiento comparada con la frecuencia propia.

====Solución matemática====
Debemos resolver la ecuación diferencial

<center><math>\ddot{x} + 2\beta \dot{x}+\omega_0^2 x = 0</math></center>

con ciertas condiciones iniciales

<center><math>x(t=0) = x_0\qquad\qquad \dot{x}(t=0)=v_0</math></center>

Esta ecuación diferencial es de las llamadas lineales (la elongación y sus derivadas no están elevadas a ninguna potencia). Para buscar una solución de una ecuación de este tipo proponemos una solución exponencial

<center><math>x = \mathrm{e}^{\lambda t}\,</math></center>

Derivando esta función

<center><math>\dot{x} = \lambda \mathrm{e}^{\lambda t}\qquad\qquad \ddot{x}=\lambda^2 \mathrm{e}^{\lambda t}</math></center>

y sustituyendo en la ecuación diferencial

<center><math>\left(\lambda^2+2\beta\lambda + \omega_0^2\right)\mathrm{e}^{\lambda t} = 0</math></center>

Puesto que la exponencial nunca puede anularse debe cumplirse que

<center><math>\lambda^2+2\beta\lambda + \omega_0^2 = 0\,</math></center>

Esta ya no es una ecuación diferencial. Es una ecuación de segundo grado cuyas soluciones nos dan los valores posibles de <math>\lambda</math>. Puesto que existen dos valores, la solución de la ecuación diferencial se escribe como la combinación

<center><math>x(t) = c_1 \mathrm{e}^{\lambda_1 t}+ c_2 \mathrm{e}^{\lambda_2 t}</math></center>

donde <math>c_1</math> y <math>c_2</math> son dos constantes cuyos valores se calculan a partir de las condiciones iniciales.

Resolviendo la ecuación de segundo grado nos quedan las soluciones

<center><math>\lambda_1 = -\beta + \sqrt{\beta^2-\omega_0^2}\qquad\qquad \lambda_2 = -\beta-\sqrt{\beta^2-\omega_0^2}</math></center>

Vemos que, como se dijo antes, dependiendo del valor de <math>\beta</math> hay tres posibilidades, dependiendo del signo de lo que hay dentro de la raíz cuadrada

* Si <math>\beta > \omega_0</math> las dos soluciones son reales y diferentes (caso ''sobreamortiguado'').
* Si <math>\beta = \omega_0</math> existe una solución real doble (''amortiguamiento crítico'').
* Si <math>\beta < \omega_0</math> las dos soluciones son complejas conjugadas (caso ''subamortiguado'').

Cada uno de estos casos merece un estudio por separado.

===Caso sobreamortiguado (&beta; > &omega;<sub>0</sub>)===
Consideraremos en primer lugar el caso de rozamiento intenso

<center><math>\beta > \omega _0\,</math></center>

En este caso las dos raíces de la ecuación son reales y además negativas 

<center><math>\lambda_1 = -\beta + \sqrt{\beta^2-\omega_0^2}\qquad \qquad \lambda_2 = -\beta -\sqrt{\beta^2-\omega_0^2}</math></center>

(para ver que la primera es negativa basta con observar que la raíz es menor que <math>\beta</math>). La solución de la ecuación diferencial es entonces una suma de dos exponenciales decrecientes

<center><math>x(t) = c_1 \mathrm{e}^{-|\lambda_1| t}+c_2 \mathrm{e}^{-|\lambda_2| t}</math></center>

Puesto que <math>|\lambda_2| > |\lambda_1|</math> la segunda exponencial decae más rápidamente, y es la primera de las dos la que determina el tiempo en decaer.

Por dar un ejemplo numérico, supongamos que <math>\omega_0=1\,\mathrm{s}^{-1}</math> y que <math>\beta=1.25\,\mathrm{s}^{-1}</math>. En ese caso resultan

<center> <math>\lambda_1=-0.5\,\mathrm{s}^{-1}\qquad \qquad \lambda_1=-2.0\,\mathrm{s}^{-1}</math></center>

Esto quiere decir que la primera exponencial decae en un tiempo típico de 2 segundos (la inversa de <math>\lambda_1</math>) mientras que la segunda lo hace en medio segundo, por tanto al cabo de un segundo prácticamente ya solo tenemos la primera exponencial.

<center>[[Archivo:exponenciales-01.png]]</center>

La solución completa es la combinación de las dos aunque rápidamente se asemeja mucho a la primera

<center>[[Archivo:exponenciales-02.png]]</center>

El caso de esta figura representa la situación en que se tira del oscilador y se libera desde el reposo (siendo nula su velocidad inicial, lo que corresponde a una tangente horizontal). El resorte tiende lentamente a la posición de equilibrio.

El que la solución sea una combinación de exponenciales decrecientes no quiere decir que la solución sea decreciente en todo instante. Por ejemplo, imaginemos el caso de un muelle que es golpeado en la posición de equilibrio. La masa se aleja originalmente de la posición de equilibrio, para luego retornar lentamente a ella.

<center>[[Archivo:exponenciales-03.png]]</center>

===Caso subamortiguado (&beta; < &omega;<sub>0</sub>)===
El caso opuesto al anterior lo obtenemos cuando el rozamiento es débil (incluyendo el caso en que no hay rozamiento). 

Si llamamos

<center><math>\omega = \sqrt{\omega_0^2-\beta^2}</math></center>

podemos escribir las dos soluciones de la ecuación de segundo grado como complejos conjugados

<center><math>\lambda_1 = -\beta+\mathrm{j}\omega \qquad\qquad \lambda_2 = -\beta -\mathrm{j}\omega\,</math></center>

siendo <math>\mathrm{j}=\sqrt{-1}</math> la unidad imaginaria. La solución general de la ecuación diferencial queda entonces

<center><math>x(t) = c_1 \mathrm{e}^{(-\beta + \mathrm{j}\omega)t}+ c_2 \mathrm{e}^{(-\beta - \mathrm{j}\omega)t}\,</math></center>

Aquí podemos extraer como factor común la parte real de la exponencial y escribir

<center><math>x(t)=\mathrm{e}^{-\beta t}\left(c_1 \mathrm{e}^{\mathrm{j}\omega t}+ c_2 \mathrm{e}^{-\mathrm{j}\omega t}\right)</math></center>

Para ver que esta solución representa oscilaciones amortiguadas aplicamos la fórmula de Euler

<center><math>\mathrm{e}^{\mathrm{j}\omega t} = \cos(\omega t) + \mathrm{j}\,\mathrm{sen}(\omega t)</math></center>

que transforma la solución en

<center><math>x(t) = \mathrm{e}^{-\beta t}\left(b_1\cos(\omega t) + b_2\,\mathrm{sen}(\omega t)\right)</math></center>

con

<center><math>b_1 = c_1+c_2 \qquad \qquad b_2 = \mathrm{j}(c_1-c_2)</math></center>

Esta combinación de senos y cosenos puede reducirse a uno solo, como se hace el caso del oscilador sin rozamiento, y escribir la solución en la forma

<center><math>x(t) = A_0\mathrm{e}^{-\beta t}\cos(\omega t+\varphi)</math></center>

Podemos leer esta solución como una oscilación sinusoidal 

<center><math>x(t) = A(t) \cos(\omega t+\varphi)</math></center>

con una amplitud que decae exponencialmente

<center><math>A(t) = A_0\mathrm{e}^{-\beta t}\,</math></center>

<center>[[Archivo:exponenciales-04.png]]</center>

Este comportamiento se dice cuasiperiódico, porque no llega a repetirse (al completar una oscilación no se encuentra en la misma posición que al iniciarla). El cuasiperiodo es mayor que el del oscilador sin rozamiento

<center><math>\omega < \omega_0\qquad\Rightarrow\qquad T = \frac{2\pi}{\omega} > \frac{2\pi}{\omega_0} = T_0</math></center>

El tiempo que tarda en decaer la amplitud no los da el factor de decaimiento <math>\beta</math>. En un tiempo 

<center><math>\tau = \frac{1}{\beta}</math></center>

la amplitud se reduce en un factor e (a un 36.8% de la que tuviera). Tenemos entonces dos escalas de tiempo: <math>T_0</math> nos mide el tiempo que tarda en oscilar, <math>\tau</math> el tiempo que tarda en amortiguarse. El cociente adimensional

<center><math>\frac{\tau}{T_0}=\frac{\omega_0}{2\pi \beta}</math></center>

nos mide la importancia del amortiguamiento pues nos da el número de oscilaciones en un tiempo típico de decaimiento. Si este número es grande quiere decir que el oscilador es muy poco amortiguado.

Comparando las oscilaciones con y sin rozamiento vemos que si éste es pequeño se nota un cambio apreciable en la amplitud, pero muy pequeño en el peridoo

<center>[[Archivo:exponenciales-05.png]]</center>

Si el rozamiento es grande, de forma que <math>\beta \simeq \omega_0</math> el decaimiento es muy importante y el periodo de oscilación es tan largo que prácticamente la partícula no llega a realizar ninguna oscilación.

<center>[[Archivo:exponenciales-06.png]]</center>

En el ejemplo de la figura (<math>\beta = 0.9\omega_0</math>) el primer mínimo está por debajo del eje solo -0.0006m y es inapreciable en la gráfica.

===Amortiguamiento crítico (&beta; = &omega;<sub>0</sub>)===
El tercer caso es uno particular que se da muy raramente, ya que requiere unos valores concretos de los parámetros. Para el caso del muelle con rozamiento debe cumplirse

<center><math>\beta = \omega_0 \qquad\Rightarrow\qquad \frac{\gamma}{2m}=\sqrt{\frac{k}{m}}\qquad\Rightarrow\qquad \gamma = 2\sqrt{km}</math></center>

La constante de rozamiento debe tener este valor exacto. Si es un poco mayor ya el movimiento es sobreamortiaguado; si es un poco menor, subamortiguado.

En el caso del amortiguamiento crítico, puede demostrarse que la solución es de la forma

<center><math>x(t) = (c_1+c_2t)\mathrm{e}^{-\beta t}\,</math></center>

Gráficamente esta función presenta un decaimiento exponencial, similar al caso sobreamortiguado.
<center>[[Archivo:exponenciales-07.png]]</center>
El amortiguamiento crítico posee una propiedad que lo hace interesante desde el punto de vista técnico: es el caso en que se retorna más rápidamente a la posición de equilibrio. Si la fricción es menor, el oscilador va y viene y tarda más en pararse. Si es mayor, la fuerte fricción ralentiza el movimiento y tarda también más en pararse. El tiempo típico de parada en el caso subamortiguado es

<center><math>\tau = \frac{1}{\beta}</math></center>

y en el sobreamortiguado

<center><math>\tau = \frac{1}{|\lambda_1|}=\frac{1}{\beta-\sqrt{\beta^2-\omega_0^2}}</math></center>

Representando este tiempo como función de <math>\beta</math> para una frecuencia propia fijada de 1&thinsp;s<sup>&minus;1</sup> vemos como efectivamente es mínimo cuando se da el amortiguamiento crítico.

<center>[[Archivo:exponenciales-08.png]]</center> 

Por ello, los amortiguadores de los automóviles y demás maquinaria procuran ajustarse al valor crítico, ya que de esta forma se consigue el objetivo de detener las vibraciones en el menor tiempo posible.

===Comparación de los tres casos===
Supongamos, como ejemplo, una frecuencia propia de 1\,rad/s y que las condiciones iniciales son que liberamos la partícula desde una cierta distancia de 10&thinsp;cm, es decir <math>x_0=0.1\,\mathrm{m}</math> y <math>v_0=0\,\mathrm{m}/\mathrm{s}</math>. Comparando los caso de <math>\beta = 0.5\mathrm{s}^{-1}</math>, <math>\beta = 1.0\mathrm{s}^{-1}</math>, y <math>\beta = 1.5\mathrm{s}^{-1}</math> obtenemos las curvas siguientes:

<center>[[Archivo:exponenciales-09.png]]</center>

Vemos como la crítica y la sobreamortiguada decaen siendo la crítica más rápida en hacerlo, mientras que la subamortiguada presenta oscilaciones rápidamente decrecientes.

Si consideramos el caso de una masa que parte del equilibrio, pero con una cierta velocidad de 1m/s (es decir <math>x_0=0\,\mathrm{m}</math>, <math>v_0 =1\,\mathrm{m}/\mathrm{s}</math> obtenemos para los mismos valores

<center>[[Archivo:exponenciales-10.png]]</center>

En este caso las curvas no se limitan a decaer, pues hay un alejamiento inicial. La partícula termina volviendo a la posición de equilibrio, siendo el camino más rápido el del amortiguamiento crítico. Vemos también que aunque en todos los casos parte con la misma velocidad el máximo alejamiento es menor cuanto mayor sea el rozamiento.

===Energía en un oscilador amortiguado===
Una de las consecuencias del amortiguamiento es la disipación de energía mecánica. De acuerdo con el teorema de las fuerzas vivas

<center><math>\frac{\mathrm{d}K}{\mathrm{d}t}=P=P_\mathrm{c}+P_\mathrm{nc}</math></center>

siendo <math>K</math> la energía cinética

<center><math>K = \frac{1}{2}mv^2</math></center>

y <math>P_c</math> la potencia de las fuerzas conservativas (la elástica, en este caso) y <math>P_{nc}</math> la de las no conservativas (que sería la de rozamiento). La potencia de las fuerzas conservativas verifica

<center><math>P_\mathrm{c} = -\frac{\mathrm{d}U}{\mathrm{d}t}</math></center>

con <math>U</math> la energía potencial elástica

<center><math>U=\frac{1}{2}kx^2</math></center>

Pasando la energía potencial al primer miembro obtenemos la energía mecánica

<center><math>\frac{\mathrm{d}E}{\mathrm{d}t}=\frac{\mathrm{d}\ }{\mathrm{d}t}(K+U) = P_\mathrm{nc}</math></center>

La potencia de las fuerzas no conservativas la calculamos multiplicando la fuerza por la velocidad

<center><math>P_\mathrm{nc}=F_rv = -\gamma v^2\,</math></center>

de forma que nos queda la relación

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{1}{2}mv^2+\frac{1}{2}kx^2\right)=-\gamma v^2</math></center>

Puesto que el segundo miembro es siempre negativo, esta ecuación nos dice que la energía mecánica se va disipando progresivamente, aunque no a ritmo uniforme: el consumo es mayor cuando lo es la rapidez del movimiento.

En el caso del amortiguamiento muy débil, si representamos la energía mecánica en una curva de potencial vemos como la energía mecánica va descendiendo a medida que la partícula va y viene

<center>[[Archivo:exponenciales-14.png]]</center>

Si realizamos gráficas equivalentes para los ejemplos anteriores vemos un decaimiento muy rápido de la energía. Puesto que en el caso crítico y el sobreamortiguado no llega a realizar una oscilación completa, lo que vemos es que la curva de energía "cae en picado", siendo la disipación tanto más rápida cuanto mayor es la constante de amortiguamiento.

<center>[[Archivo:exponenciales-12.png]]{{qquad}}{{qquad}}[[Archivo:exponenciales-11.png]]</center>

==Oscilaciones forzadas==
===Introducción===
Un oscilador armónico amortiguado es aquel que, en adición a la fuerza recuperadora dada por la ley de Hooke, experimenta una fuerza de [[rozamiento viscoso]] proporcional a la velocidad. 

<center><math>m\vec{v}=-k\vec{r}-\gamma\vec{v}\,</math></center>

Si este oscilador se mueve a lo largo de una recta, la segunda ley de Newton se reduce a 

<center><math>ma = -kx - \gamma v\,</math></center>

donde <math>x</math> es la elongación (distancia a la posición de equilibrio dada por la longitud natural).

Si este oscilador amortiguado se encuentra sometido a una fuerza externa adicional, en general dependiente del tiempo, se dice que el oscilador está ''forzado'', siendo su ecuación de movimiento

<center><math>ma = -kx - \gamma v + F(t)\,</math></center>

De entre los posibles tipos de fuerza que se pueden aplicar, la más importante desde el punto de vista físico, es aquella que en sí misma es oscilante, esto es

<center><math>ma = -kx - \gamma v + F_0\cos(\omega t)\,</math></center>

donde la frecuencia de oscilación de la fuerza no tiene por qué coincidir con la frecuencia propia del oscilador armónico

<center><math>\omega\neq \omega_0=\sqrt{\frac{k}{m}}</math></center>

===El caso de una fuerza constante===
El caso más sencillo de un oscilador forzado sería aquél sometido a una fuerza que no depende del tiempo <math>F=F_0</math> (<math>\omega = 0</math>). Este sería el caso de un muelle que cuelga verticalmente y del cual se cuelga una masa <math>m</math>, siendo su ecuación de movimiento

<center><math>m\ddot{z}+\gamma \dot{z} + k z = -mg</math></center>

====Solución estacionaria====
La solución en este caso es fácil de imaginar. Si de una balanza colgamos 1kg de plátanos, el muelle se estira y en el estado final se queda en reposo en una posición <math>z_s</math> inferior a la inicial. En este estado equilibrio, al ser constante, las derivadas respecto al tiempo se anulan y nos queda

<center><math>z=z_\mathrm{s} \qquad \dot{z}=0\qquad \ddot{z}=0</math></center>

y sustituyendo obtenemos la posición de equilibrio

<center><math>m\cdot 0 + \gamma\cdot 0 + k z_s = -mg\qquad\Rightarrow\qquad z_s = -\frac{mg}{k}</math></center>

Esta ''estado estacionario'' es independiente de las condiciones iniciales. Ni la posición inicial ni la velocidad inicial aparecen en su expresión.

====Solución transitoria====
En el ejemplo de los plátanos que se cuelgan de la balanza, la experiencia nos muestra que la posición final no se alcanza instantáneamente, sino que la balanza oscila unas cuantas veces antes de detenerse. Se dice entonces que tenemos un proceso ''transitorio'' que representa el paso desde el estado inicial hasta el final.

Para describir matemáticamente el estado transitorio introducimos como variable la diferencia entre la posición instantánea y la de equilibrio

<center><math>s = z - z_s\qquad\Rightarrow\qquad z = s + z_s</math></center>

Teniendo en cuenta que la solución estacionaria es una de equilibrio

<center><math>\dot{z}=\dot{s}+\dot{z}_s = \dot{s}\qquad\qquad \ddot{z}=\ddot{s}+\ddot{z}_s = \ddot{s}</math></center>

y sustituyendo en la ecuación del oscilador forzado queda

<center><math>m\ddot{s}+\gamma\dot{s}+k s +k z_s = -mg</math></center>

pero <math>kz_s = -mg</math> así que esto se reduce a

<center><math>m\ddot{s}+\gamma\dot{s}+k s =0</math></center>

que es la ecuación del [[Oscilaciones amortiguadas (GIE)|oscilador armónico amortiguado]]. Dependiendo del grado de amortiguamiento podemos tener un caso subamortiguado, crítico o sobreamortiguado. En cualquiera de los tres casos el comportamiento general es el mismo: la solución decae exponencialmente a cero. 
Este estado transitorio sí depende de las condiciones iniciales, que determinan las constantes de la solución.

====Superposición====
Según esto, la solución de un oscilador forzado se compone de dos partes:

<center><math>z(t) = z_s(t) + z_t(t)\,</math></center>

* Un estado transitorio inicial, <math>z_t(t)</math>, que depende de las condiciones iniciales y que decae exponencialmente.
* Un estado estacionario final, <math>z_s(t)</math>, que solo depende de la fuerza aplicada y no del estado inicial.

Este mismo principio se aplica al caso de que tengamos una fuerza oscilante en lugar de una constante.
<center>[[Archivo:osc-forzadas-01.png]]</center>

===Caso de una fuerza sinusoidal===
Consideraremos ahora el caso de que la fuerza aplicada oscile sinusoidalmente. 

<center><math>F(t) = F_0\cos(\omega t)\,</math></center>

Este caso contiene el anterior en el límite <math>\omega\to 0</math>.

Si escribimos la ecuación de movimiento expresando la velocidad y la aceleración como derivadas de la posición nos queda la ecuación

<center><math>m\ddot{x}+\gamma \dot{x}+kx = F_0\cos(\omega t)</math></center>

Dividiendo por la masa podemos escribirla empleando los parámetros definidos al estudiar las oscilaciones amortiguadas

<center><math>\ddot{x}+2\beta\dot{x}+\omega_0^2x = \frac{F_0}{m}\cos(\omega t)</math></center>

donde

<center><math>\beta = \frac{\gamma}{2m}\qquad\qquad \omega_0 = \sqrt{\frac{k}{m}}</math></center>

El problema general consiste en determinar la elongación como función del tiempo, <math>x(t)</math> para una posición y velocidad iniciales.

Como en el caso de una fuerza constante, la solución se puede escribir como suma de dos, una estacionaria y una transitoria, con la diferencia de que el estado estacionario ahora no es de reposo, sino de movimiento.

Cuando se aplica una fuerza oscilante a una masa atada a un muelle, el resultado es que la masa oscila. Lo hace además con la frecuencia de la señal aplicada.  La solución puede escribirse entonces en la forma

<center><math>x = A\,\cos(\omega t + \varphi)</math></center>

siendo la amplitud de las oscilaciones y <math>\varphi</math> el desfase con la fuerza aplicada (que no tiene por qué ser nulo; el que oscile con la misma frecuencia no quiere decir que esté sincronizado, puede ocurrir que cuando la fuerza sea máxima el desplazamiento sea nulo, por ejemplo). La amplitud y el desfase pueden cambiar mucho según sea la frecuencia de la fuerza aplicada, aunque esta sea siempre de la misma intensidad

<center>[[Archivo:osc-forzadas.gif]]</center>

Esto quiere decir que para conseguir una amplitud mayor no necesitamos aplicar una fuerza más intensa. Basta con elegir adecuadamente la frecuencia de su variación.
Nuestro objetivo se reduce entonces a determinar <math>A</math> y <math>\varphi</math> como funciones de la intensidad de la fuerza aplicada y de su frecuencia.

La forma más sencilla de determinar estas constantes es mediante el uso de [[Fasor (GIE)|fasores]]. La fuerza aplicada puede escribirse, con ayuda de la fórmula de Euler, como

<center><math>F(t) = F_0\cos(\omega t)=\mathrm{Re}\left(\hat{F}\mathrm{e}^{\mathrm{j}\omega t}\right)</math></center>

donde

<center><math>\hat{F}=F_0</math></center>

es el fasor, o amplitud compleja, de la fuerza. Es una cantidad compleja cuyo módulo nos da la amplitud de las oscilaciones (<math>F_0</math> en este caso) y cuyo argumento nos da el desfase (nulo, en este caso).

En la solución estacionaria sinusoidal, la elongación admite una expresión análoga

<center><math>x(t) = A\cos(\omega t+\varphi)=\mathrm{Re}\left(\hat{x}\mathrm{e}^{\mathrm{j}\omega t}\right)</math></center>

donde ahora

<center><math>\hat{x}=A\mathrm{e}^{\mathrm{j}\varphi}\,</math></center>

Este fasor es una constante compleja contiene tanto la amplitud como la constante de fase de la solución, por lo que si determinamos el fasor, ya hemos resuelto nuestro problema.

La ventaja de usar fasores es que transforma las derivadas en multiplicaciones. La [[Fasor (GIE)#Velocidad_en_un_MAS|velocidad]] y la [[Fasor (GIE)#Aceleraci.C3.B3n_en_un_MAS|aceleración]] pueden expresarse también en forma fasorial, siendo sus amplitudes complejas

<center><math>\hat{v}=\mathrm{j}\omega\hat{x}</math>{{qquad}}{{qquad}}<math>\hat{a}=(\mathrm{j}\omega)^2\hat{x}=-\omega^2\hat{x}</math></center>

Puesto que en la ecuación de movimiento

<center><math>a + 2\beta v + \omega_0^2x = \frac{F_0}{m}\cos(\omega t)\,</math></center>

todas las cantidades (<math>F</math>, <math>x</math>, <math>v</math> y <math>a</math>) oscilan con la misma frecuencia, podemos escribirla como una relación entre fasores

<center><math>\hat{a} + 2\beta \hat{v} + \omega_0^2\hat{x} = \frac{\hat{F}}{m}</math></center>

Sustituyendo los fasores de la velocidad y la aceleración nos queda

<center><math>-\omega^2\hat{x} +2\mathrm{j} \beta \omega\hat{x} + \omega_0^2\hat{x} = \frac{\hat{F}}{m}</math></center>

Esta es una ecuación algebraica de la cual es inmediato despejar el fasor de la elongación

<center><math>\hat{x}=\frac{\hat{F}/m}{(\omega_0^2-\omega^2)+2\mathrm{j}\omega \beta}</math></center>

;Amplitud:La amplitud de las oscilaciones la hallamos como el módulo del fasor. A su vez, el módulo de un cociente de números complejos es el cociente de los módulos, por lo que

<center><math>A = \left|\hat{x}\right| = \frac{F_0/m}{\sqrt{(\omega_0^2-\omega^2)^2+4\beta^2\omega^2)}}</math></center>

;Desfase: El desfase entre las oscilaciones de la elongación y la fuerza aplicada lo da el argumento del fasor

<center><math>\phi=\mathrm{arg}(\hat{x})=-\mathrm{arctg}\left(\frac{2\beta\omega b}{\omega_0^2-\omega^2}\right)</math></center>

===Dependencia con la frecuencia===
Los fasores de la posición y de la velocidad (y por tanto sus amplitudes y desfases) son funciones de la frecuencia, lo cual quiere decir que la respuesta de un oscilador a una fuerza que varía sinusoidalmente no depende solo de la magnitud de la fuerza. Una fuerza pequeña puede tener un gran efecto, si la frecuencia es adecuada, o, por el contrario, una gran fuerza aplicada puede no tener efecto apreciable. Dependiendo del sistema que se esté estudiando, puede interesar una cosa o la otra. Así, para un receptor de radio, interesa que se amplifiquen las frecuencias de las señales de la frecuencia deseada, y no las del resto, esto se consigue ajustando principalmente dos parámetros, la frecuencia propia del oscilador y su coeficiente de amortiguamiento, definidos por

<center><math>\omega_0=\sqrt{\frac{k}{m}}</math>{{qquad}}{{qquad}}<math>\beta = \frac{b}{2m}</math></center>

====Amplitud====
La amplitud de la elongación, en términos de la frecuencia propia y del factor de amortiguamiento queda en la forma

<center><math>A(\omega) = \frac{F_0/m}{\sqrt{(\omega_0^2-\omega^2)^2+4\beta^2\omega^2}}</math></center>

;Bajas frecuencias: Cuando <math>\omega\to 0</math>, la amplitud de las oscilaciones tiende a 

<center><math> A(\omega)\simeq \frac{F_0/m}{\omega_0^2}=\frac{F_0}{k} \qquad (\omega\ll\omega_0)</math></center>

:que es el comportamiento que uno obtiene para una fuerza no oscilante, como el peso, para el cual la frecuencia sería estrictamente nula.

;Altas frecuencias: Si <math>\omega\to\infty</math> Para frecuencias altas la amplitud tiende a cero

<center><math> A(\omega) \simeq 0 \qquad(\omega\gg\omega_0)</math></center>

:Esto quiere decir que si un oscilador lo excitamos mediante una fuerza cuya frecuencia sea mucho mayor que la propia del oscilador, éste prácticamente no se ve afectado.

Para valores intermedios de la frecuencia, de valor comparable a la frecuencia propia <math>\omega_0</math>, podemos tener un máximo de amplitud o no tenerlo dependiendo del grado de amortiguamiento. La amplitud es máxima cuando lo que hay dentro de la raíz del denominador es mínimo, lo cual ocurre para

<center><math>\omega_\mathrm{max}=\sqrt{\omega_0^2-2\beta^2}</math></center>

como se comprueba sin más que derivar el radicando e igualar a cero. Gráficamente obtenemos una curva con un máximo que es más agudo cuanto menor sea el coeficiente de amortiguamiento <math>\beta</math>. Para representarlo conviene emplear una escala logarítmica para las frecuencias ya que de esta forma se puede barrer un rango más amplio sin comprimir la gráfica.

<center>[[Archivo:osc-forzadas-02.png]]{{qquad}}{{qquad}}[[Archivo:osc-forzadas-03.png]]</center>

Cuando existe este máximo (que puede ser muy destacado) se dice que el sistema posee una ''resonancia''. Físicamente este resultado nos dice que ajustando la frecuencia de la fuerza aplicada podemos obtener una respuesta de gran amplitud (infinita en el caso sin rozamiento). Por ejemplo, para el caso <math>\beta = 0.1\omega_0</math> obtenemos que las amplitudes de las oscilaciones resultantes es 5 veces lo que se estira en el caso estático. La altura del máximo es

<center><math>A_\mathrm{max}=\frac{F_0/m}{2\beta\sqrt{\omega_0^2-\beta^2}}</math></center>

Para valores de <math>\beta</math> pequeños la posición del máximo ocurre aproximadamente a la frecuencia propia (por lo que también se la llama '''frecuencia de resonancia''') y la altura es inversamente proporcional al amortiguamiento

<center><math> \omega_\mathrm{max}\simeq \omega_0\qquad\qquad A_\mathrm{max}\simeq \frac{F_0/m}{2\beta\omega_0}\qquad\qquad(\beta\ll\omega_0)</math></center>

Este efecto puede ser deseable (es la base de los sintonizadores de radio, por ejemplo), pero también puede ser indeseado. A menudo es la causa de los ruidos molestos que producen ordenadores y otros aparatos, por vibraciones de los ventiladores y otras partes móviles. La razón de que muchos de estos ruidos desaparezcan cuando se golpea el aparato es porque al hacerlo se modifica la frecuencia propia y el sistema deja de estar en una resonancia. 

Una resonancia de gran amplitud puede tener efectos dramáticos, provocando la fractura de partes del sistema. El caso más famoso de este fenómeno es el colapso del puente de  Tacoma Narrows en Washington (USA).

<center>[[Archivo:tacoma-narrows-bridge.jpg]]</center>

([http://www.youtube.com/watch?v=3mclp9QmCGs vídeo de la catástrofe en Youtube]).

====Fase====
El movimiento del oscilador no tiene por qué ir en fase con la fuerza. Es posible que para una frecuencia dada cuando la fuerza vaya hacia abajo la partícula se mueva hacia arriba (ya que la fuerza nos da la aceleración, pero no directamente la posición).

La diferencia de fase lo hallamos a partir del argumento del fasor. Usando la fórmula para el argumento de un cociente de números complejos

<center><math>\varphi = \mathrm{arg}(\hat{x})=\mathrm{arg}(F_0/m)-\mathrm{arg}\left((\omega_0^2-\omega^2)+4\mathrm{j}\beta\omega\right)=-\mathrm{arctg}\left(\frac{2\beta\omega}{\omega_0^2-\omega^2}\right)</math></center>

Para frecuencias bajas (<math>\omega\to 0</math>) se obtiene desfase 0. Esto corresponde a que si la fuerza es constante en un sentido (por ejemplo, el peso, hacia abajo), el resultado es un desplzamiento constante en el mismo sentido. Justo en la resonancia el desfase es de <math>\pi/2</math>, es decir, cuando la fuerza es máxima el desplazamiento es nulo y viceversa. Para frecuencias grandes el desfase tiende a <math>\pi</math> lo que implica que el desplazamiento va en sentido contrario a la fuerza.

<center>[[Archivo:osc-forzadas-04.png]]</center>

===Velocidad===
A menudo es más interesante conocer la velocidad del oscilador forzado, más que su posición. La velocidad también se puede expresar en forma fasorial

<center><math>v(t) = V\cos(\omega t + \theta)=\mathrm{Re}\left(\hat{v}\mathrm{e}^{\mathrm{j}\omega t}\right)\qquad\qquad \hat{v}=V\mathrm{e}^{\mathrm{j}\theta}</math></center>

Obtenemos el fasor de la velocidad multiplicando el de la elongación por <math>\mathrm{j}\omega</math>

<center><math>\hat{v}=\mathrm{j}\omega\hat{x}=\frac{\mathrm{j}\omega \hat{F}/m}{(\omega_0-\omega^2)+2\mathrm{j}\beta\omega}</math></center>

Dividiendo en el numerador y el denominador por <math>\mathrm{j}\omega</math>

<center><math>\hat{v}=\frac{\hat{F}/m}{2\beta+\mathrm{j}(\omega-\omega_0^2/\omega)}</math></center>

====Amplitud====
La amplitud de las oscilaciones de la velocidad es

<center><math>V(\omega) = |\hat{v}|=\frac{F_0/m}{\sqrt{4\beta^2+(\omega - \omega_0^2/\omega)^2}}</math></center>

De nuevo tenemos una función dependiente de la frecuencia. 

;Bajas frecuencias: Cuando <math>\omega\to 0</math>, el término <math>\omega_0^2/\omega</math> tiende a infinito y la amplitud se va a cero. Esto quiere decir que para una fuerza que varía muy lentamente, el oscilador se mueve muy lentamente, de forma que en todo momento está casi en la posición de equilibrio calculada anteriormente.

;Altas frecuencias: Cuando <math>\omega\to \infty</math>, de nuevo el dernominador tiende a infinito y la amplitud se va de nuevo a cero. Esto quiere decir que para una fuerza que varía muy rápidamente, la inercia del oscilador le impide responder y este se limita a vibrar ligeramente alrededor de su posición de equilibrio <math>x=0</math>.

;Resonancia: En la raíz que aparece en la amplitud hay una suma de dos términos positivos. El primero no se anula nunca, pero el segundo se hace cero si

<center><math>\omega = \frac{\omega_0^2}{\omega}</math>{{tose}}<math> \omega = \omega_0 = \sqrt{\frac{k}{m}}</math></center>

:esto es, si la frecuencia de la fuerza aplicada coincide con la frecuencia propia del oscilador, la amplitud de las oscilaciones en la velocidad es máxima. Vemos que, a diferencia de lo que ocurre con la elongación, el máximo en la amplitud de la velocidad está siempre centrado en la misma frecuencia (conocida como ''frecuencia de resonancia''). El valor máximo de esta amplitud será

<center><math>V_\mathrm{max}=V(\omega_0)=\frac{F_0}{\gamma}</math></center>

:que tiende a infinito cuando el rozamiento se hace nulo.

<center>[[Archivo:osc-forzadas-06.png]]</center>

Si se emplea una escala logarítmica para las frecuencias el resultado es una curva simétrica alrededor de la frecuencia de resonancia.

====Desfase====
El desfase entre las oscilaciones de la velocidad y la fuerza aplicada es igual a

<center><math>\theta=\mathrm{arg}(\hat{v})=\mathrm{arctg}\left(\frac{\omega_0^2-\omega^2}{2\beta\omega}\right)</math></center>

La curva es igual al desfase de la elongación, pero desplazada en <math>\pi/2</math>. En la frecuencia de resonancia la velocidad y la fuerza van en fase

<center>[[Archivo:osc-forzadas-05.png]]</center>

===Energía y potencia===
<!--
====Energía almacenada====
Por tratarse de un oscilador armónico, el oscilador posee una energía mecánica, que será una función oscilante con el tiempo

<center><math>E=\frac{1}{2}mv^2+\frac{1}{2}kx^2</math></center> 

Sustituyendo la posición instantánea

<center><math>x = A\cos(\omega t+\varphi)\qquad\qquad v = \dot{x}=-\omega A\,\mathrm{sen}\,(\omega t+\varphi)</math></center>

queda

<center><math>E=\frac{A^2}{2}\left(m\omega^2\,\mathrm{sen}^2(\omega t+\phi)+k\cos^2(\omega t +\phi)\right)</math></center>

El valor promedio en un periodo se define como

<center><math>\left\langle E\right\rangle = \frac{1}{T}\int_0^T E\,\mathrm{d}t</math></center>

Aplicando que

<center><math>\frac{1}{T}\int_0^T \mathrm{sen}^2(\omega t+\phi)\,\mathrm{d}t=\frac{1}{T}\int_0^T \cos^2(\omega t+\phi)\,\mathrm{d}t = \frac{1}{2}</math></center>

obtenemos el valor medio de la energía

<center><math>\left\langle E\right\rangle(\omega) = \frac{A^2}{4}\left(m\omega^2+k\right)=\frac{F_0^2(m\omega^2+k)}{4((m\omega^2-k)^2+\gamma^2\omega^2)}</math></center>

En particular, en la frecuencia de resonancia, <math>\omega_0</math> esta energía vale

<center><math>\left\langle E\right\rangle(\omega_0)= \frac{F_0^2m}{2b^2}</math></center>

====Energía disipada====
-->
Al existir una fuerza de rozamiento, la energía se está transformando continuamente en calor. Puesto que la energía almacenada, en promedio, se conserva, esta energía que se pierde debe ser aportada regularmente por un agente externo, en este caso, la fuerza aplicada sobre la partícula.

La [[Teorema_de_conservación_de_la_energía_mecánica#Presencia_de_fuerzas_no_conservativas|potencia con la que se disipa la energía]] es

<center><math>P_d=P_\mathrm{nc}=\mathbf{F}_{nc}\cdot\mathbf{v}=-\gamma v^2</math></center>

Sustituyendo la velocidad como función del tiempo

<center><math>P_d=-\gamma V^2\cos^2(\omega t+\theta)</math></center>

La energía total disipada durante un periodo será, aplicando que

<center><math>\int_0^{2\pi}\cos^2(x)\mathrm{d}x=\pi</math></center>

igual a

<center><math>W_d=\int_0^P P_\mathrm{nc}\,\mathrm{d}t = -\frac{\pi \gamma V^2}{\omega}</math></center>

Sustituyendo el valor de la amplitud en la velocidad

<center><math>W_d(\omega)=-\frac{2\pi \beta \omega F_0^2/m}{(\omega_0^2-\omega^2)^2+4\beta^2\omega^2}</math></center>

siendo la potencia promedio
<center><math>
\langle P_d\rangle = \frac{W_d}{T}=-\frac{\beta \omega^2 F_0^2/m}{(\omega_0^2-\omega^2)^2+4\beta^2\omega^2}</math></center>

Representando gráficamente esta potencia promedio disipada como función de la frecuencia se obtiene de nuevo una curva con un máximo en la frecuencia de resonancia

<center>[[Archivo:osc-forzadas-07.png]]</center>

Esto quiere decir que la frecuencia de resonancia es aquella en la que el sistema absorbe más energía de la fuerza aplicada, es decir, es la forma más efectiva de transmitir energía a un sistema. Si se emplea una frecuencia diferente el proceso ya no es tan eficiente.

En un ejemplo sencillo, esto nos está diciendo que si queremos impulsar un columpio, por ejemplo, debemos aplicar una fuerza oscilante cuya frecuencia debe ser lo más parecida posible a la propia del oscilador. No sirve de nada que unas veces empujemos el columpio cuando se va a lejando de nosotros y otras cuando se está acercando. Hay que hacerlo siempre cuando se está alejando. De esta forma se le comunica el máximo de energía posible y la amplitud de las oscilaciones resultantes es la mayor posible.

===Factor de calidad===
[[Imagen:Qfactor.gif|300px|thumb|Ancho de banda]]

El ''factor de calidad'', <math>Q</math>, de un oscilador mide cómo de agudo es el pico de una resonancia. Se define, en términos de la energía, como

<center><math>Q=2\pi\frac{\mbox{Energia almacenada promedio}}{\mbox{Energia disipada en un periodo}}= 2\pi\frac{\left\langle E\right\rangle}{|W_d|}</math></center>

donde las cantidades se evalúan en la frecuencia de resonancia <math>\omega_0</math>. Calculando y sustituyendo los valores de la energía almacenada y de la disipada queda

<center><math>Q= \frac{\sqrt{km}}{\gamma}</math></center>

Gráficamente, para resonancias muy agudas, el factor de calidad es inversamente proporcional al ''ancho de banda'', que es la anchura del pico de la curva de la amplitud como función de la frecuencia, medida a media altura

<center><math>Q\simeq \frac{\omega_0}{\Delta\omega}</math></center>
[[Categoría:Mecánica de la partícula y de los sistemas (CMR)]]