Javier Valcarce's Homepage

VYO - Dinámica del MAS. Amortiguación. Resonancia.

Dinámica de los sistemas que se mueven con movimiento armónico simple (MAS). Amortiguación. Resonancia.

Definiciones

Tensión mecánica
Fuerza / Área, $\sigma = F / A$
Módulo de Young
Caracteriza el comportamiento de un material elástico, se define como $Y = \frac{\sigma}{\epsilon} = \frac{F/S}{\Delta L/L}$
siendo $\sigma$ la tensión mecánica aplicada y $\epsilon$ la deformación unitaria o relativa del material
Oscilador libre no amortiguado
$\ddot{x} + \omega_0^2 x = 0$
Oscilador libre amortiguado
$\ddot{x} + 2\gamma\dot{x} + \omega_0^2 x = 0$
Oscilador forzado no amortiguado
$\ddot{x} + \omega_0^2 x = \frac{F_0}{m}\cos(\omega t + \alpha_0)$
Oscilador forzado amortiguado
$\ddot{x} + 2\gamma\dot{x} + \omega_0^2 x = \frac{F_0}{m}\cos(\omega t + \alpha_0)$

No es físicamente posible un oscilador no amortiguado. Todos tienen en mayor o menor grado un coeficiente de amortiguamiento. Los modelos no amortiguados tienen un interés principalmente matemático.

Método general para EDOs lineales homogéneas de grado 2 con coeficientes constantes

$$ \lambda^2 + 2\gamma\lambda + \omega_0^2 = 0 \\ \lambda = -\gamma \pm i\omega_\gamma \longrightarrow \lambda_1, \lambda_2 $$ $$ x(t) = C_1e^{\lambda_1 t} + C_2e^{\lambda_2 t} $$ aplicando las condiciones iniciales obtenemos un sistema de ecuaciones algebráicas con el que obtenemos $C_1$ y $C_2$

Oscilador libre no amortiguado

Es una caso particular del oscilador amortiguado (con $\gamma = 0$)

Oscilador libre amortiguado

Muelle vs amortiguador

Un muelle ejerce una fuerza directamente proporcional ($k$) al desplazamiento relativo de los extremos mientras que un amortiguador ejerce una fuerza directamente proporcional ($b$) a la velocidad relativa de sus extremos. El vector fuerza es siempre opuesto al vector velocidad. $$ \begin{align} F_k &= -k (x_2 - x_1) \\ F_b &= -b \frac{d}{dt}(x_2 - x_1) = -b (v_2 - v_1) \\ \end{align} $$

Factor de calidad $Q$ y amortiguamiento $\zeta$

Está relacionado con la energía que se pierde en cada periodo del oscilador: $$ Q = 2\pi\frac{E_{total}}{E_{perdida-por-periodo}} = \frac{1}{2}\frac{\omega_0}{\gamma} = \frac{1}{2}\frac{1}{\zeta} $$ En un oscilador armónico la energía total es $E_0 = \frac{1}{2}kA^2$ y el decaimiento de esta energía total por disipación (solo en el caso de amortiguamiento subcrítico, si no no habría oscilación) es $E(t) = E_0e^{-2\gamma t}$

Niveles de amortiguamiento $\gamma$

  • Subcrítico: $\gamma^2 < \omega_0^2$
    En este caso la oscilación se desvanece paulatinamente, los cruces por cero no están exactamente en el punto medio entre el mínimo y el máximo. La solución tiene la forma $$ x(t) = C e^{-\gamma t} \sin(\omega_\gamma t + \delta) = C e^{-\gamma t} \cos\left(\omega_\gamma t + \delta - \frac{\pi}{2}\right) \\ \omega_\gamma = \sqrt{ \omega_0^2 - \gamma^2 } \\ C^2 = x_0^2 + \frac{(v_0 + \gamma x_0)^2}{\omega_\gamma^2} \\ \delta = arctan\left(\frac{\omega_\gamma x_0 }{ v_0 + \gamma x_0}\right) \\ $$
    • tiempo de relajación $\tau$
    • decremento logarítmico de la amplitud $\Delta$
    • coeficiente de amortiguación $\zeta = \gamma/\omega_0$
  • Crítico: $\gamma^2 = \omega_0^2$
    En este caso no hay oscilación. La solución tiene la forma $$ x(t) = (A_1 + A_2 t)e^{-\gamma t} $$
  • Supercrítico: $\gamma^2 > \omega_0^2$
    En este caso no hay oscilación. La solución tiene la forma $$ \begin{align} x(t) &= A_1 e^{-(\gamma + \beta)t} + A_2 e^{-(\gamma - \beta)t} \\ \beta &= \sqrt{\gamma^2 - \omega_0^2} \\ \end{align} $$
En todos los casos tenemos 2 constantes ajustables, que se determinan con las condiciones iniciales del sistema.

Oscilador forzado no amortiguado

Es una caso particular del oscilador forzado amortiguado (con $\gamma = 0$)

Forzamiento: régimen transitorio y régimen permanente

En el modelo del oscilador forzado amortiguado la solución es realmente la suma de dos términos

  • la solución de la ecuación homogénea $x_h(t)$
  • una solución particular $x_p(t)$
Una de ellas es un término evanescente (multiplicado por $e^-{\gamma t/2}$) que desaparece al cabo de un tiempo, dicho término es transitorio. El otro corresponde al régimen permanente y tiene una frecuencia de oscilación igual a la del generador externo, el oscilador "aprende" tras un tiempo a seguir la fuerza externa.

Oscilador forzado amortiguado

Al aplicar una fuerza oscilante $F_0\cos(\omega t + \alpha_0)$ en el término independiente de la ecuación del sistema la solución tiene dos partes: $$ x(t) = \color{blue}{x_h(t)} + \color{green}{x_p(t)} = \color{blue}{C e^{-\gamma t}\sin(\omega_\gamma t + \delta)} + \color{green}{ A \sin(\omega t + \alpha_0 + \beta) } $$

  • un término $x_h$ evanescente multiplicado por $e^{-\gamma t}$ que corresponde a la ecuación homogénea y por lo tanto es la misma que la del oscilador libre (amortiguado o no). Es importante durante el régimen transitorio. La amplitud $C$ y constante de fase $\delta$ depende de las condiciones iniciales ($t=0$), ojo, a dichas condiciones iniciales contribuye también la componente estacionaria $x_p(t)$
  • un término $x_p$ permanente que tiene la misma frecuencia $\omega$ que la frecuencia forzante. Es el término que corresponde al régimen sinusoidal permanente. Esta componente de la solución no depende de las condiciones iniciales ($t=0$)

Vamos a ocuparnos del término permanente:
La amplitud de oscilación $A$ depende de la relación entre la frecuencia natural $\omega_0$ y la frecuencia forzante $\omega$ $$ A = \frac{F_0/m}{\sqrt{ (\omega^2-\omega_0^2)^2 + 4\gamma^2\omega^2 }} \qquad \beta = -arctan\left( \omega^2-\omega_0^2 \over 2\gamma\omega \right) $$ Si la solución es en forma de coseno entonces hay que invertir la fracción de la expresión de $\beta$. En la literatura se pueden encontrar numerosas variantes de esta expresión dependiendo de si $\beta$ aparece sumando o restando, de si es seno o coseno, etc.

Concepto de ganancia o factor de amplificación de un oscilador forzado: $$ F_{amp} = \frac{A}{\delta_{est}} = \frac{A}{\frac{F_0/m}{\omega_0^2}} = \frac{1}{\sqrt{ \left(1-(\frac{\omega}{\omega_0})^2\right)^2 + \left(2\zeta\frac{\omega}{\omega_0}\right)^2 }} $$

curvas de resonancia

Energía y potencia en el oscilador forzado. Resonancia.

Potencia suministrada por la fuerza externa $$ P(t) = F(t) \dot{x}(t) $$ La potencia absorbida tiene ciclos en que es positiva y luego negativa, eso significa que la energía absorbida por el sistema luego se devuelve a la fuente. La potencia disipada es justamente la diferencia, por tanto, potencia absorbidad y potencia disipada no son lo mismo.

Tabla resumen de soluciones

Tipo Ecuación Solución
Oscilador libre no amortiguado $ \ddot{q} + \omega^2 q = 0 $ $$ q(t) = A\cos\left(\omega t + \phi\right) = A\sin\left(\omega t + \phi + \frac{\pi}{2}\right) $$ siendo $$ \begin{align} A &= \sqrt{q_0^2 + \dot{q}_0^2 / \omega^2 } \\ \phi &= \arccos\left( +\frac{q_0}{A} \right) = \arcsin\left( -\frac{\dot{q_0}}{\omega A} \right) = \arctan\left(-\frac{\dot{q}_0 }{q_0\omega}\right) \\ \end{align} $$
Oscilador libre amortiguado con amortiguación de tipo $\vec{F}=-b\vec{v}$ $\ddot{x} + 2\gamma\dot{x} + \omega_0^2 x = 0$

Dependiendo del nivel de amortiguamiento puede haber tres soluciones distintas:
$\gamma = \frac{b}{2m} \quad \omega_0^2=\frac{k}{m} \quad \omega_\gamma = \sqrt{\omega_0^2 - \gamma^2}$

  • subcrítico: $(\omega_0^2 > \gamma^2;\ \omega_\gamma \in \mathbb{R})$
    $x(t) = A_1 e^{-\gamma t} \sin(\omega_\gamma t + \delta)$
    $A_1^2 = x_0^2 + \frac{(v_0 + \gamma x_0)^2}{\omega_\gamma}$
  • crítico: $(\omega_0^2 = \gamma^2;\ \omega_\gamma = 0)$
    $x(t) = (A_1 + A_2 t)e^{-\gamma t}$
  • supercrítico: $(\omega_0^2 < \gamma^2;\ \omega_\gamma \in \mathbb{C})$
    $x(t) = A_1 e^{-(\gamma + \beta)t} + A_2 e^{-(\gamma - \beta)t}$

Oscilador forzado no amortiguado $\ddot{x} + \omega_0^2 x = \frac{F_0}{m}\cos(\omega t)$ $ x(t) = x_h(t) + x_p(t) $
$x(t) = B\cos(\omega_0 t + \beta) + C\cos(\omega t)$
$C=\frac{F_0/m}{\omega_0^2 - \omega^2}$
Oscilador forzado amortiguado $\ddot{x} + 2\gamma\dot{x} + \omega_0^2 x = \frac{F_0}{m}\cos(\omega t + \alpha_0)$ La solución es la suma de un término evanescente y de otro permanente

$ x(t) = \color{blue}{x_h(t)} + \color{green}{x_p(t)} $
$ x(t) = \color{blue}{C e^{-\gamma t}\sin(\omega_\gamma t + \delta)} +\color{green}{ A \sin(\omega t + \alpha_0 + \beta) } $
$\omega_\gamma = (\omega_0^2 - \gamma^2)^{1/2} $
$ A = \frac{F_0/m}{\sqrt{(\omega^2-\omega_0^2)^2 + 4\gamma^2\omega^2}} $
$\beta = -\arctan(\frac{\omega^2 - \omega_0^2}{2\gamma\omega})$
$C$ y $\delta$ son funciones de la posición y la velocidad inicial.

This page forms part of website https://javiervalcarce.eu