Entrar Página Discusión Historial Go to the site toolbox

Oscilaciones amortiguadas (GIE)

De Laplace

(Diferencias entre revisiones)
(Amortiguamiento crítico)
(Amortiguamiento crítico)
Línea 260: Línea 260:
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.
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==
+
==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>
==Energía en un oscilador amortiguado==
==Energía en un oscilador amortiguado==
[[Categoría:Movimiento oscilatorio (GIE)]]
[[Categoría:Movimiento oscilatorio (GIE)]]

Revisión de 19:21 6 dic 2012

Contenido

1 El oscilador no amortiguado

En otras secciones se estudia la cinemática y la 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

F = -kx\,

siendo x 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

ma = -kx\qquad\rightarrow\qquad a= \ddot{x}=-\frac{k}{m}x

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

\ddot = -\omega_0^2x

siendo en este caso la frecuencia natural

\omega_0=\sqrt{\frac{k}{m}}

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

x(t) = A\cos(\omega_0 t+\varphi)\,
Archivo:Muelle.gif

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

x(t+T) = x(t)\qquad\qquad T = \frac{2\pi}{\omega}

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

x(t) = b_1\cos(\omega_0 t)+b_2\,\mathrm{sen}(\omega_0 t)

con

b_1=A\cos(\varphi) \qquad\qquad b_2 = -A\,\mathrm{sen}(\varphi)

Los valores de las constantes b1 y b2 pueden calcularse también a partir de las condiciones iniciales del movimiento

b_1= x_0 = x(t=0)\qquad\qquad b_2=\frac{v_0}{\omega_0}=\frac{v(t=0)}{\omega_0}

2 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í. 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.

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

Archivo:amortiguador-dentro.gif

Si un coche no tuviera suspensión (es decir, si el chasis estuviera unido rígidamente unido 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.

Archivo:bache.gif

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

\vec{F}=-\gamma\vec{v}

y para el caso particular del movimiento rectilíneo

F = -\gamma v = -\gamma\dot{x}

3 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

ma = -kx-\gamma v\,

Pasando todo al primer miembro

ma + \gamma v + k x = 0\,

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

m \ddot{x} + \gamma \dot{x}+k x = 0

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

\ddot{x}+2\beta \dot{x} +\omega_0^2 x = 0

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

\omega_0 = \sqrt{\frac{k}{m}}

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

\beta = \frac{\gamma}{2m}

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

Tanto la frecuencia propia ω0 como la constante de amortiguamiento β tienen dimensiones de inversa de un tiempo y se miden en s−1 en el SI.

4 Solución de la ecuación

4.1 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, debemos esperar que no llegue a oscile, sino que simplemente se vaya 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 β como ω0 tienen las mismas dimensiones y por tanto se pueden comparar. Establecemos entonces el criterio

  • Rozamiento débil: β < ω0
  • Rozamiento intenso: β > ω0

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.

4.2 Solución matemática

Debemos resolver la ecuación diferencial

\ddot{x} + 2\beta \dot{x}+\omega_0^2 x = 0

con ciertas condiciones iniciales

x(t=0) = x_0\qquad\qquad \dot{x}(t=0)=v_0

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

x = eλt

Derivando esta función

\dot{x} = \lambda \mathrm{e}^{\lambda t}\qquad\qquad \ddot{x}=\lambda^2 \mathrm{e}^{\lambda t}

y sustituyendo en la ecuación diferencial

\left(\lambda^2+2\beta\lambda + \omega_0^2\right)\mathrm{e}^{\lambda t} = 0

Puesto que la exponencial nunca puede anularse debe cumplirse que

\lambda^2+2\beta\lambda + \omega_0^2 = 0\,

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

x(t) = c_1 \mathrm{e}^{\lambda_1 t}+ c_2 \mathrm{e}^{\lambda_2 t}

donde c1 y c2 son dos constantes cuyos valores se calculan a partir de las condiciones iniciales.

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

\lambda_1 = -\beta + \sqrt{\beta^2-\omega_0^2}\qquad\qquad \lambda_2 = -\beta-\sqrt{\beta^2-\omega_0^2}

Vemos que, como se dijo antes, dependiendo del valor de β hay tres posibilidades, dependiendo del signo de lo que hay dentro de la raíz cuadrada

  • Si β > ω0 las dos soluciones son reales y diferentes
  • Si β = ω0 existe una solución real doble.
  • Si β < ω0 las dos soluciones son complejas conjugadas.

Cada uno de estos casos merece un estudio por separado.

5 Caso sobreamortiguado (β > ω0)

Consideraremos en primer lugar el caso de rozamiento intenso

\beta > \omega _0\,

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

\lambda_1 = -\beta + \sqrt{\beta^2-\omega_0^2}\qquad \qquad \lambda_2 = -\beta -\sqrt{\beta^2-\omega_0^2}

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

x(t) = c_1 \mathrm{e}^{-|\lambda_1| t}+c_2 \mathrm{e}^{-|\lambda_2| t}

Puesto que | λ2 | > | λ1 | 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 \omega_0=1\,\mathrm{s}^{-1} y que \beta=1.25\,\mathrm{s}^{-1}. En ese caso resultan

\lambda_1=-0.5\,\mathrm{s}^{-1}\qquad \qquad \lambda_1=-2.0\,\mathrm{s}^{-1}

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

Archivo:exponenciales-01.png

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

Archivo:exponenciales-02.png

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.

Archivo:exponenciales-03.png

6 Caso subamortiguado (β < ω0)

El caso opuesto al anterior lo obtenemos cuando el rozamiento es débil (incluyendo el caso en que no hay rozamiento).

Si llamamos

\omega = \sqrt{\omega_0^2-\beta^2}

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

\lambda_1 = -\beta+\mathrm{j}\omega \qquad\qquad \lambda_2 = -\beta -\mathrm{j}\omega\,

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

x(t) = c_1 \mathrm{e}^{(-\beta + \mathrm{j}\omega)t}+ c_2 \mathrm{e}^{(-\beta - \mathrm{j}\omega)t}\,

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

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

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

\mathrm{e}^{\mathrm{j}\omega t} = \cos(\omega t) + \mathrm{j}\,\mathrm{sen}(\omega t)

que transforma la solución en

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

con

b_1 = c_1+c_2 \qquad \qquad b_2 = \mathrm{j}(c_1-c_2)

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

x(t) = A_0\mathrm{e}^{-\beta t}\cos(\omega t+\varphi)

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

x(t) = A(t) \cos(\omega t+\varphi)

con una amplitud que decae exponencialmente

A(t) = A_0\mathrm{e}^{-\beta t}\,
Archivo:exponenciales-04.png

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

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

El tiempo que tarda en decaer la amplitud no los da el factor de decaimiento β. En un tiempo

\tau = \frac{1}{\beta}

la amplitud se reduce en un factor e (a un 36.8% de la que tuviera). Tenemos entonces dos escalas de tiempo: T0 nos mide el tiempo que tarda en oscilar, τ el tiempo que tarda en amortiguarse. El cociente adimensional

\frac{\tau}{T_0}=\frac{\omega_0}{2\pi \beta}

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

Archivo:exponenciales-05.png

Si el rozamiento es grande, de forma que \beta \simeq \omega_0 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.

Archivo:exponenciales-06.png

En el ejemplo de la figura (β = 0.9ω0) el primer mínimo está por debajo del eje solo -0.0006m y es inapreciable en la gráfica.

7 Amortiguamiento crítico, β = ω0

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

\beta = \omega_0 \qquad\Rightarrow\qquad \frac{\gamma}{2m}=\sqrt{\frac{k}{m}}\qquad\Rightarrow\qquad \gamma = 2\sqrt{km}

8 Energía en un oscilador amortiguado

Herramientas:

Herramientas personales
TOOLBOX
LANGUAGES
licencia de Creative Commons
Aviso legal - Acerca de Laplace