8.2: Oscilador armónico amortiguado
- Page ID
- 129614
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Hasta ahora no hemos tenido en cuenta la amortiguación de nuestros osciladores armónicos, lo que por supuesto no es muy realista. La principal fuente de amortiguación para una masa en un resorte se debe al arrastre de la masa cuando se mueve a través del aire (o cualquier fluido, ya sea gas o líquido). Para velocidades relativamente bajas, las fuerzas de arrastre en un objeto escalan linealmente con la velocidad del objeto, como lo ilustra la ley de Stokes (Ecuación 2.2.5). Para un objeto de forma arbitraria que se mueve a través de un fluido arbitrario escribiremos\(F_{\mathrm{drag}}=-\gamma \dot{x}\), con\(\gamma\) el coeficiente de arrastre, y por supuesto opuesto a la dirección del movimiento. Añadiendo esto a la fuerza del resorte da para la ecuación de movimiento del oscilador armónico amortiguado:
\[m \ddot{x}=-\gamma \dot{x}-k x \label{dampedharmosc}\]
Ahora tenemos dos números que determinan el movimiento: la frecuencia no amortiguada\(\omega_{0}=\sqrt{k / m}\) y la relación de amortiguación\(\zeta=\gamma / 2 \sqrt{m k}\). En términos de estos parámetros, podemos reescribir la ecuación\ ref {dampedharmosc} como:
\[\ddot{x}+2 \zeta \omega_{0} \dot{x}+\omega_{0}^{2} x=0 \label{rewritten}\]
La solución de la Ecuación\ ref {reescrito} depende fuertemente del valor de\(\zeta\), ver Figura 8.2.1. Podemos encontrarlo 2 sustituyendo el Ansatz\(x(t)=e^{\lambda t}\), que da una ecuación característica para\(\lambda\):
\[\lambda^{2}+2 \zeta \omega_{0} \lambda+\omega_{0}^{2}=0\]
por lo
\[\lambda=-\zeta \omega_{0} \pm \omega_{0} \sqrt{\zeta^{2}-1} \label{lambda}\]
Para\(\zeta < 1\), existen dos soluciones complejas para\(\lambda\), y encontramos que\(x(t)\) sufre una oscilación con una amplitud exponencialmente decreciente:
\[x(t)=e^{-\zeta \omega_{0} t}\left[A \cos \left(\omega_{d} t\right)+B \sin \left(\omega_{d} t\right)\right]\]
donde\(\omega_{d}=\omega_{0} \sqrt{1-\zeta^{2}}\) y A y B se derivan de las condiciones iniciales. Debido a que todavía hay una oscilación, este tipo de movimiento se llama subamortiguado. En contraste, si\(\zeta > 1\), las raíces\(\lambda_{ \pm}\) en la Ecuación\ ref {lambda} son reales, y obtenemos un comportamiento cualitativamente diferente, sobreamortiguado, en el que x vuelve a 0 con una decadencia exponencial sin oscilaciones:
\[x(t)=A e^{\lambda_{+} t}+B e^{\lambda_{-} t}=e^{-\zeta \omega_{0} t}\left[A e^{\Omega t}+B e^{-\Omega t}\right]\]
donde\(\Omega=\omega_{0} \sqrt{\zeta^{2}-1}\). Naturalmente, el caso límite es cuando\(\zeta = 1\), que es un oscilador amortiguado críticamente, el retorno más rápido a 0 sin oscilaciones. Porque en este caso la Ecuación\ ref {lambda} solo tiene una raíz, nuevamente obtenemos una solución cualitativamente diferente:
\[x(t)=(A+B t) e^{-\omega_{0} t}\]
Los tres casos diferentes y la oscilación no amortiguada se muestran en la Figura 8.2.1.
2 Véase el apéndice A.3.2 para conocer los detalles matemáticos sobre cómo resolver ecuaciones generales de este tipo.