Ver el código de «Oscilaciones amortiguadas (GIE)»

==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{x} = -\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 paro 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 debe seleccionarse modificando las propiedades del sistema, 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 x_0 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>

[[Categoría:Movimiento oscilatorio (GIE)]]