Archivo de la etiqueta: movimiento libre amortiguado

Ecuaciones Diferenciales I: Oscilaciones mecánicas

Por Omar González Franco

Las matemáticas comienzan a parecerse demasiado a resolver un puzzle. La física también,
pero son puzzles creados por la naturaleza, no por la mente del hombre.
– Maria Goeppert Mayer

Introducción

En esta entrada estudiaremos algunos tipos de oscilaciones mecánicas con el propósito de poner en práctica los métodos desarrollados hasta este momento de la segunda unidad.

Consideraremos varios sistemas dinámicos lineales en los que cada modelo matemático será una ecuación diferencial de segundo orden con coeficientes constantes acompañada de condiciones iniciales especificadas en un tiempo que tomaremos como $t = 0$.

$$a \dfrac{d^{2}x}{dt^{2}} + b \dfrac{dx}{dt} + cx = g(t); \hspace{1cm} x(0) = x_{0}, \hspace{0.5cm} x^{\prime}(0) = x_{1} \label{1} \tag{1}$$

Con respecto a la notación, denotaremos con $x$ a la variable dependiente que físicamente representará la posición de un objeto, mientras que $t$ será la variable independiente y representara al tiempo, ya que nuestro propósito es describir el movimiento oscilatorio de un objeto a través del tiempo.

A la función $g(x)$ de (\ref{1}) la llamaremos entrada o función forzada del sistema. Una solución $x(t)$ de (\ref{1}) en un intervalo $\delta$ que contiene a $t = 0$ y satisface las condiciones iniciales se le llama salida o respuesta del sistema.

El sistema dinámico que estudiaremos será el de resorte – objeto y los tipos de movimiento que describiremos será el movimiento libre no amortiguado, el movimiento libre amortiguado y el movimiento forzado.

Movimiento libre no amortiguado

Consideremos un resorte de longitud $l$ suspendido verticalmente de un soporte rígido y en la parte inferior del resorte se encuentra un objeto de masa $m$, el peso del objeto hace que el resorte se elongue una distancia $s$. En la posición de equilibrio establecemos que $x = 0$, tal como se muestra en la siguiente figura.

Resorte sin objeto y resorte con el objeto de masa $m$ en la posición de equilibrio.

Es claro que la cantidad de alargamiento o elongación del resorte depende de la masa, además el resorte mismo ejerce una fuerza restauradora $F$ opuesta a la dirección de elongación y proporcional a la cantidad de elongación $s$, esta característica corresponde a la ley de Hooke y matemáticamente se expresa como

$$F = ks \label{2} \tag{2}$$

donde $k$ es una constante de proporcionalidad llamada constante de resorte.

Una vez colocado el objeto de masa $m$, el resorte se alarga una distancia $s$ y mantiene una posición de equilibrio en el que el peso $W$ del objeto se equilibra con la fuerza restauradora $F$ del resorte. Recordando que el peso de un objeto es

$$W = mg \label{3} \tag{3}$$

con $m$ la masa del objeto y $g$ la aceleración de la gravedad, podemos establecer que en el equilibrio ocurre que

$$W = F \label{4} \tag{4}$$

o bien,

$$mg -ks = 0 \label{5} \tag{5}$$

Si el objeto se desplaza una cantidad $x$ de su posición de equilibrio, la fuerza restauradora del resorte será

$$F_{x} = k(s + x) \label{6} \tag{6}$$

Objeto en reposo y objeto en movimiento desplazado una distancia $x$.

Como estamos analizando un movimiento no amortiguado, vamos a suponer que no hay fuerzas restauradoras que actúen sobre el sistema y que el objeto oscila libre de otras fuerzas externas. Entonces podemos igualar la segunda ley de Newton con la fuerza resultante de la fuerza restauradora y el peso.

$$m \dfrac{d^{2}x}{dt^{2}} = -k(s + x)+ mg = -kx + mg -ks$$

Considerando (\ref{5}) obtenemos que

$$m \dfrac{d^{2}x}{dt^{2}} = -kx \label{7} \tag{7}$$

El signo negativo indica que la fuerza restauradora del resorte actúa en dirección opuesta a la dirección del movimiento, además se toma la convención de que la dirección hacia abajo de la posición de equilibrio es positiva.

Si dividimos entre $m$ la ecuación (\ref{7}) y reordenamos obtenemos la ecuación diferencial

$$\dfrac{d^{2}x}{dt^{2}} + \dfrac{k}{m}x = 0 \label{8} \tag{8}$$

Veremos más adelante la razón por la que es conveniente definir la constante

$$\omega^{2} = \dfrac{k}{m} \label{9} \tag{9}$$

Usando esta definición podemos escribir la ecuación (\ref{8}) como

$$\dfrac{d^{2}x}{dt^{2}} + \omega^{2}x = 0 \label{10} \tag{10}$$

La ecuación diferencial (\ref{10}) se dice que describe el movimiento armónico simple o movimiento libre no amortiguado. Dos condiciones iniciales claras son el desplazamiento inicial

$$x(0) = x_{0}$$

y la velocidad inicial

$$x^{\prime}(0) = x_{1}$$

del objeto. Por ejemplo, si $x_{0} > 0$, entonces indica que el objeto parte de un punto por debajo de la posición de equilibrio lo que provocará una velocidad impartida hacia arriba, es decir, $x_{1} < 0$. Cuando $x^{\prime}(0) = 0$ el objeto se libera a partir del reposo. Y si por ejemplo $x_{0} < 0$ y $x_{1} = 0$, entonces indica que el objeto se libera desde el reposo pero desde una posición arriba de la posición de equilibrio.

La ecuación (\ref{10}) representa el modelo matemático que describe el fenómeno, pero ahora estamos interesados en conocer la ecuación de movimiento, así que es momento de aplicar lo aprendido y resolver la ecuación diferencial.

Se trata de una ecuación diferencial de segundo orden con coeficientes constantes, así que proponemos una solución de la forma

$$x(t) = e^{rt} \label{11} \tag{11}$$

Utilizamos la letra $r$ y no $k$ para no confundirnos con la constante de resorte. Al sustituir esta solución y su segunda derivada en (\ref{10}) se obtiene la ecuación auxiliar

$$r^{2} + \omega^{2} = 0 \label{12} \tag{12}$$

Las raíces son

$$r_{1} = i \omega \hspace{1cm} y \hspace{1cm} r_{2} = -i \omega$$

Identificamos que $\alpha = 0$ y $\beta = \omega$, Por lo tanto, la solución general es

$$x(t) = c_{1} \cos(\omega t) + c_{2} \sin(\omega t) \label{13} \tag{13}$$

Esta solución corresponde a la ecuación general de movimiento del objeto sujeto al resorte considerando que no hay amortiguación.

Si se aplican las condiciones iniciales y se determinan las constantes $c_{1}$ y $c_{2}$, entonces habremos encontrado la ecuación de movimiento del sistema en particular. Notemos que la solución efectivamente describe un movimiento oscilatorio ya que se encuentran presentes las funciones seno y coseno.

La constante

$$\omega = \sqrt{\dfrac{k}{m}} \label{14} \tag{14}$$

se llama frecuencia circular del sistema y nos permite definir algunas cantidades. $\omega$ se mide en radianes por segundo.

La cantidad

$$T = \dfrac{2 \pi}{\omega} \label{15} \tag{15}$$

determina el periodo del movimiento descrito por (\ref{13}), es decir, representa el tiempo que tarda el objeto en hacer un ciclo de movimiento. Un ciclo es una oscilación completa del objeto. Podemos decir que el periodo $T$ es el tamaño del intervalo de tiempo entre dos máximos sucesivos (o mínimos sucesivos) de $x(t)$. De acuerdo a nuestra convención, un máximo es el desplazamiento positivo del objeto en el que alcanza su distancia máxima debajo de la posición de equilibrio, mientras que un mínimo es el desplazamiento negativo en el que alcanza su altura máxima arriba de la posición de equilibrio. En cualquier caso decimos que hay un desplazamiento extremo del objeto.

La cantidad

$$f = \dfrac{1}{T} = \dfrac{\omega}{2 \pi} \label{16} \tag{16}$$

es la frecuencia de movimiento y representa el número de ciclos completados cada segundo.

Existe una forma alterna de la solución (\ref{13}) en la que se hace explícita la amplitud $A$ de las oscilaciones. Si en la solución (\ref{13}) $c_{1} \neq 0$ y $c_{2} \neq 0$, se define la amplitud como

$$A = \sqrt{c^{2}_{1} + c^{2}_{2}} \label{17} \tag{17}$$

y se define el ángulo de fase $\phi$, tal que

$$\sin(\phi) = \dfrac{c_{1}}{A}, \hspace{1cm} \cos(\phi) = \dfrac{c_{2}}{A} \hspace{1cm} \Rightarrow \hspace{1cm} \tan(\phi) = \dfrac{c_{1}}{c_{2}} \label{18} \tag{18}$$

Notemos lo siguiente

\begin{align*}
x(t) &= c_{1} \cos(\omega t) + c_{2} \sin(\omega t) \\
&= A \dfrac{c_{1}}{A} \cos(\omega t) + A \dfrac{c_{2}}{A} \sin(\omega t) \\
&= [A \sin(\phi)] \cos(\omega t) + [A \cos(\phi)] \sin(\omega t) \\
&= A \sin(\omega t) \cos(\phi) + A \cos(\omega t) \sin(\phi)
\end{align*}

Si en la última expresión aplicamos la identidad trigonométrica

$$\sin(a + b) = \sin(a) \cos(b) + \cos(a) \sin(b) \label{19} \tag{19}$$

entonces obtenemos la solución (\ref{13}) en una forma alterna más simple.

$$x(t) = A \sin(\omega t + \phi) \label{20} \tag{20}$$

En resumen, la ecuación que describe el movimiento armónico simple o movimiento libre no amortiguado es

$$\dfrac{d^{2}x}{dt^{2}} + \omega^{2}x = 0$$

Y las soluciones que representan el movimiento del objeto son

$$x(t) = c_{1} \cos(\omega t) + c_{2} \sin(\omega t)$$

o bien,

$$x(t) = A \sin(\omega t + \phi)$$

Movimiento libre amortiguado

Es claro que el movimiento libre no amortiguado es un movimiento ideal, pues el movimiento descrito por (\ref{13}) o (\ref{20}) supone que no hay fuerzas retardadoras actuando sobre el objeto y sabemos que, a menos que el objeto este suspendido en un vacío perfecto, siempre habrá por lo menos una fuerza de resistencia debido al medio circundante, por ejemplo la resistencia del aire.

El propósito ahora, al igual que antes, es determinar la ecuación diferencial o modelo matemático que describe al sistema cuando existen fuerzas de amortiguamiento, para posteriormente determinar la ecuación general de movimiento.

Consideremos nuevamente un objeto de masa $m$ suspendido sobre un resorte con constante $k$, pero en esta ocasión consideremos que existe una fuerza externa de amortiguamiento actuando sobre el objeto. En el estudio de la mecánica, las fuerzas de amortiguamiento que actúan sobre un cuerpo se consideran proporcionales a una potencia de la velocidad instantánea $\dfrac{dx}{dt}$. En nuestro caso supondremos que la fuerza de amortiguamiento esta dada por un múltiplo constante de la velocidad, esto es

$$F_{am} = \eta \dfrac{dx}{dt} \label{21} \tag{21}$$

donde $\eta$ es una constante de amortiguamiento positiva. De esta manera, cuando ninguna otra fuerza actúa sobre el sistema, de la segunda ley de Newton se tiene que

$$m \dfrac{d^{2}x}{dt^{2}} = -kx -\eta \dfrac{dx}{dt} \label{22} \tag{22}$$

El signo negativo en la fuerza de amortiguamiento indica que dicha fuerza actúa en dirección opuesta al movimiento.

Si dividimos la ecuación diferencial (\ref{22}) por $m$ y reordenamos, obtenemos

$$\dfrac{d^{2}x}{dt^{2}} + \dfrac{\eta}{m} \dfrac{dx}{dt} + \dfrac{k}{m}x = 0 \label{23} \tag{23}$$

Recordemos que

$$\omega^{2} = \dfrac{k}{m}$$

y por convención definimos

$$2\rho = \dfrac{\eta}{m} \label{24} \tag{24}$$

Así podemos reescribir la ecuación (\ref{23}) como

$$\dfrac{d^{2}x}{dt^{2}} + 2 \rho \dfrac{dx}{dt} + \omega^{2} x = 0 \label{25} \tag{25}$$

Esta ecuación corresponde al modelo matemático que describe al fenómeno. La utilidad de $2\rho$ se hace evidente al momento de intentar resolver la ecuación, pues si se considera la solución

$$x(x) = e^{rt}$$

y se sustituye en la ecuación (\ref{25}) junto con las derivadas correspondientes se obtiene la ecuación auxiliar

$$r^{2} + 2\rho r + \omega^{2} = 0 \label{26} \tag{26}$$

De donde se obtienen las siguientes dos raíces.

$$r_{1} = -\rho + \sqrt{\rho^{2} -\omega^{2}} \hspace{1cm} y \hspace{1cm} r_{2} = -\rho -\sqrt{\rho^{2} -\omega^{2}} \label{27} \tag{27}$$

Hay tres casos posibles dependiendo del valor del discriminante

$$\Delta = \rho^{2} -\omega^{2} \label{28} \tag{28}$$

Estudiemos cada caso.

Caso 1: $\rho^{2} -\omega^{2} > 0$

En este caso decimos que el sistema se encuentra sobreamortiguado porque el coeficiente de amortiguamiento $\rho$ es más grande que la constante del resorte $k$. Como las raíces son reales y distintas, la solución de la ecuación (\ref{25}), en este caso, es

$$x(t) = c_{1} e^{r_{1}t} + c_{2} e^{r_{2}t}$$

Si sustituimos los valores de (\ref{27}) podemos reescribir la solución como

$$x(t) = e^{-\rho t} \left( c_{1} e^{\sqrt{\rho^{2} -\omega^{2}}t} + c_{2} e^{ -\sqrt{\rho^{2} -\omega^{2}}t} \right) \label{29} \tag{29}$$

Esta ecuación representa un movimiento uniforme y no oscilatorio.

Caso 2: $\rho^{2} -\omega^{2} = 0$

En este caso cualquier ligera disminución en la fuerza de amortiguamiento daría como resultado un movimiento oscilatorio, decimos que el sistema está críticamente amortiguado. Como las raíces son reales e iguales, la solución de la ecuación (\ref{25}) es

$$x(t) = c_{1} e^{r_{1}t} + c_{2}t e^{r_{1}t}$$

Si sustituimos $r_{1} = -\rho$, la solución se puede reescribir como

$$x(t) = e^{-\rho t} \left( c_{1} + c_{2}t \right) \label{30} \tag{30}$$

Caso 3: $\rho^{2} -\omega^{2} < 0$

En este caso se dice que el sistema esta subamortiguado ya que el coeficiente de amortiguamiento es más pequeño que la constante del resorte. Las raíces son complejas y están dadas de la siguiente forma

$$r_{1} = -\rho + i\sqrt{\omega^{2} -\rho^{2}} \hspace{1cm} y \hspace{1cm} r_{2} = -\rho -i\sqrt{\omega^{2} -\rho^{2}} \label{31} \tag{31}$$

Identificamos que

$$\alpha = -\rho \hspace{1cm} y \hspace{1cm} \beta = \sqrt{\omega^{2} -\rho^{2}}$$

Entonces la solución está dada por

$$x(t) = e^{-\rho t} \left[ c_{1} \cos \left( \sqrt{\omega^{2} -\rho^{2}}t \right) + c_{2} \sin \left( \sqrt{\omega^{2} -\rho^{2}}t \right) \right] \label{32} \tag{32}$$

El movimiento descrito por (\ref{32}) es oscilatorio, pero debido al coeficiente $e^{-\rho t}$ las amplitudes de oscilación tienden a cero cuando $t \rightarrow \infty$.

En todos los casos la solución contiene el factor de amortiguamiento $e^{-\rho t}$, $\rho > 0$, lo que indica que los desplazamientos del objeto se vuelven despreciables conforme el tiempo $t$ aumenta.

De manera totalmente análoga que en el caso sin amortiguamiento, cualquier solución de la forma (\ref{32}) se puede escribir de forma alterna como

$$x(t) = Ae^{-\rho t} \sin \left( \sqrt{\omega^{2} -\rho^{2}} t + \phi \right) \label{33} \tag{33}$$

donde $A$ es la amplitud de las oscilaciones y el ángulo de fase $\phi$ se determina de las ecuaciones de (\ref{18}). El coeficiente $Ae^{-\rho t}$ se denomina amplitud amortiguada de oscilaciones y debido a que (\ref{33}) no es una función periódica, el número

$$T_{c} = \dfrac{2 \pi}{\sqrt{\omega^{2} -\rho^{2}}} \label{34} \tag{34}$$

se llama cuasi periodo y es el intervalo de tiempo entre dos máximos sucesivos de $x(t)$, así mismo, el número

$$f_{c} = \dfrac{\sqrt{\omega^{2} -\rho^{2}}}{2 \pi} \label{35} \tag{35}$$

se llama cuasi frecuencia.

En resumen, la ecuación que describe el movimiento libre amortiguado es

$$\dfrac{d^{2}x}{dt^{2}} + 2 \rho \dfrac{dx}{dt} + \omega^{2} x = 0$$

Y las soluciones que representan el movimiento del objeto, dependiendo del valor del discriminante de la ecuación auxiliar, son

  • Si $\rho^{2} -\omega^{2} > 0 \hspace{1cm} \rightarrow \hspace{1cm} x(t) = e^{-\rho t} \left( c_{1} e^{\sqrt{\rho^{2} -\omega^{2}}t} + c_{2} e^{ -\sqrt{\rho^{2} -\omega^{2}}t} \right)$
  • Si $\rho^{2} -\omega^{2} = 0 \hspace{1cm} \rightarrow \hspace{1cm} x(t) = e^{-\rho t} \left( c_{1} + c_{2}t \right)$
  • Si $\rho^{2} -\omega^{2} < 0 \hspace{1cm} \rightarrow \hspace{1cm} x(t) = Ae^{-\rho t} \sin \left( \sqrt{\omega^{2} -\rho^{2}} t + \phi \right)$

Movimiento forzado

Imaginemos que ahora, adicional a las situaciones anteriores, se ejerce una fuerza externa sobre el soporte del resorte. En los dos casos anteriores considerábamos al soporte fijo, pero en esta ocasión pensamos en una fuerza motriz que causa un movimiento vertical oscilatorio del soporte del resorte. Sea $F_{ext}(t)$ dicha fuerza externa, usando la segunda ley de Newton, la ecuación diferencial queda de la siguiente forma.

$$m \dfrac{d^{2}x}{dt^{2}} = -kx -\eta \dfrac{dx}{dt} + F_{ext}(t) \label{36} \tag{36}$$

Si dividimos la ecuación por $m$ y definimos

$$g(t) = \dfrac{F_{ext}(t)}{m} \label{37} \tag{37}$$

además de considerar las definiciones anteriores (\ref{9}) y (\ref{24}), podemos escribir la ecuación diferencial como

$$\dfrac{d^{2}x}{dt^{2}} + 2\rho \dfrac{dx}{dt} + \omega^{2}x = g(t) \label{38} \tag{38}$$

La ecuación (\ref{38}) representa el modelo matemático que describe al sistema con movimiento forzado. Esta ecuación es no homogénea, de manera que puede resolverse usando el método de coeficientes indeterminados o el de variación de parámetros.

Cuando $g$ es una función periódica, como

$$g(t) = g_{0} \sin(\lambda t) \hspace{1cm} o \hspace{1cm} g(t) = g_{0} \cos(\lambda t)$$

con $\lambda$ una constante, la solución general de (\ref{38}) para $\rho > 0$ es la suma de una función no periódica $x_{np}(t)$ (solución complementaria o solución de la ecuación homogénea asociada) y una función periódica $x_{p}(t)$ (solución particular de la ecuación no homogénea), en la que $x_{np}(t)$ se desvanece a medida que el tiempo incrementa, es decir,

$$\lim_{t \to \infty} x_{np}(t) = 0 \label{39} \tag{39}$$

Esta propiedad nos indica que para valores grandes de tiempo, los desplazamientos del objeto se aproximan mediante la solución particular $x_{p}(t)$.

La función complementaria $x_{np}(t)$ se denomina término transitorio o solución transitoria, mientras que la solución particular $x_{p}(t)$ se denomina término de estado estable o solución de estado estable.

Realicemos un ejemplo en el que apliquemos cada caso

Ejemplo: Considerar un resorte sujeto de manera vertical a un soporte. El resorte se estira $50 cm$ al aplicarle una fuerza de $4N$. En la parte inferior del resorte se coloca un objeto con peso de $19.6 N$. Al objeto se le aleja de su posición de equilibrio jalándolo $1 m$ hacia abajo, si se suelta sin aplicarle una velocidad inicial, estudiar el movimiento del objeto en los siguientes casos:

  • No hay resistencia del aire (movimiento libre no amortiguado).
  • Hay resistencia del aire y es de $F_{am} = 8\dfrac{dx}{dt}$ (movimiento libre amortiguado).
  • Además de la resistencia del aire, hay una fuerza aplicada al soporte de $F_{ext}(t) = 80\sin(2t)$ (movimiento forzado).

Solución: El peso del objeto es

$$W = 19.6 N$$

entonces su masa es

$$m = \dfrac{W}{g} = \dfrac{19.6}{9.8}$$

Es decir, $m = 2 kg$. Por otro lado, si el resorte se estira

$$s = 0.5 m$$

aplicando una fuerza de

$$F = 4N$$

por la ley de Hooke tenemos que la constante del resorte es

$$k = \dfrac{F}{s} = \dfrac{4}{0.5}$$

Es decir, $k = 8 N/m$. Las condiciones iniciales son $x(0) = 1$ (posición fuera de la posición de equilibrio) y $x^{\prime}(0) = 0$ (sin velocidad inicial).

Para la primera situación sabemos que

$$F_{am} = 0 \hspace{1cm} y \hspace{1cm} F_{ext} = 0$$

De manera que la ecuación que describe al sistema es (\ref{10}) con

$$\omega^{2} = \dfrac{k}{m} = \dfrac{8}{2} = 4$$

Así, la ecuación a resolver es

$$\dfrac{d^{2}x}{dt^{2}} + 4x = 0$$

La ecuación auxiliar es

$$r^{2} + 4 = 0$$

Las raíces son $r_{1} = i2$ y $r_{2} = -i2$. Identificamos que $\alpha = 0$ y $\beta = 2$. Por lo tanto, la solución general es

$$x(t) = c_{1} \cos(2t) + c_{2} \sin(2t)$$

Para aplicar las condiciones iniciales debemos conocer la expresión de la primera derivada de la solución, dicha expresión es

$$\dfrac{dx}{dt} = -2c_{1} \sin(2t) + 2c_{2} \cos(2t)$$

Aplicando las condiciones iniciales, tenemos

$$x(0) = c_{1} = 1 \hspace{1cm} y \hspace{1cm} x^{\prime}(0) = 2c_{2} = 0$$

De donde $c_{1} = 1$ y $c_{2} = 0$. Por lo tanto, la ecuación de movimiento es

$$x(t) = \cos(2t)$$

Esta solución representa un movimiento armónico de amplitud

$$A = 1 m$$

periodo

$$T = \dfrac{2 \pi}{2} = \pi seg$$

y frecuencia

$$f = \dfrac{1}{\pi } = 0.318 \dfrac{ciclos}{segundo}$$

A continuación se muestra una gráfica con el movimiento descrito por el objeto.

Función de movimiento del objeto.

De la gráfica observamos que el objeto siempre se mantendrá oscilando de la misma manera para $t \to \infty$ y tiene sentido ya que no existe ninguna fuerza exterior que lo amortigüe.

Consideremos ahora la resistencia del aire

$$F_{am} = 8 \dfrac{dx}{dt}$$

En este caso la ecuación a resolver es de la forma (\ref{22}) y es

$$m \dfrac{d^{2}x}{dt^{2}} + kx + 8\dfrac{dx}{dt} = 0$$

que adaptando a nuestros datos se tiene

$$\dfrac{d^{2}x}{dt^{2}} + 4 \dfrac{dx}{dt} + 4x = 0$$

La ecuación auxiliar es

$$r^{2} + 4r + 4 = 0$$

Las raíces son $r_{1} = r_{2}= -2$, como son iguales, entonces la solución es de la forma

$$x(t) = e^{-2t}(c_{1} + c_{2}t)$$

La derivada es

$$\dfrac{dx}{dt} = -2e^{-2t}(c_{1} + c_{2}t) + c_{2}e^{-2t}$$

Apliquemos las condiciones iniciales.

$$x(0) = c_{1} = 1 \hspace{1cm} y \hspace{1cm} x^{\prime}(0) = -2c_{1} + c_{2} = 0$$

de donde obtenemos que $c_{1} = 1$ y $c_{2} = 2$. Por lo tanto, la ecuación de movimiento es

$$x(t) = e^{-2t}(1 + 2t)$$

El factor de amortiguamiento es $e^{-2t}$.

A continuación se muestra una gráfica con el movimiento descrito por el objeto.

Función de movimiento del objeto.

De la gráfica observamos que no hay movimiento oscilatorio, sino que el objeto llega a la posición de equilibrio y se mantiene, esto se debe al factor de amortiguamiento.

Para la situación final tenemos un movimiento forzado con una fuerza externa

$$F_{ext}(t) = 80 \sin(2t)$$

La ecuación diferencial que tenemos en este caso es

$$m \dfrac{d^{2}x}{dt^{2}} = -kx -8 \dfrac{dx}{dt} + 80 \sin(2t)$$

o bien

$$\dfrac{d^{2}x}{dt^{2}} + 4 \dfrac{dx}{dt} + 4x = 40 \sin(2t)$$

La solución de la ecuación homogénea ya la conocemos, ya que corresponde a la solución transitoria obtenida anteriormente.

$$x_{np}(t) = e^{-2t}(c_{1} + c_{2}t)$$

Para el caso no homogéneo se puede aplicar variación de parámetros o coeficientes indeterminados, apliquemos el segundo método.

Estamos en condiciones del punto 3. Proponemos una solución de la forma

$$x(t) = A \cos(2t) + B \sin(2t)$$

La primera y segunda derivada están dadas de la siguiente forma.

$$\dfrac{dx}{dt} = -2A \sin(2t) + 2B \cos(2t) \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}x}{dt^{2}} = -4A \cos(2t) -4B \sin(2t)$$

Sustituimos en la ecuación diferencial.

\begin{align*}
\left[ -4A \cos(2t) -4B \sin(2B) \right] &+ 4 \left[ -2A \sin(2t) + 2B \cos(2t) \right] + 4 \left[ A \cos(2t) + B \sin(2t) \right] \\
&= 40 \sin(2t)
\end{align*}

\begin{align*}
-8A \sin(2t) + 8B \cos(2t) &= 40 \sin(2t) \\
-A \sin(2t) + B \cos(2t) &= 5 \sin(2t)
\end{align*}

Para que se cumpla la igualdad debe ocurrir que $A = -5$ y $B = 0$, entonces la solución de estado estable es

$$x_{p}(t) = -5 \cos(2t)$$

Por lo tanto, la solución general es

$$x(t) = x_{np}(t) + x_{p}(t) = e^{-2t}(c_{1} + c_{2}t) -5 \cos(2t)$$

Vemos que

$$\dfrac{dx}{dt} = -2e^{-2t}(c_{1} + c_{2}t) + c_{2}e^{-2t} + 10 \sin(2t)$$

Apliquemos las condiciones iniciales.

$$x(0) = c_{1} -5 = 1 \hspace{1cm} y \hspace{1cm} x^{\prime}(0) = -2c_{1} + c_{2} = 0$$

de donde $c_{1} = 6$ y $c_{2} = 12$. Entonces la ecuación de movimiento es

$$x(t) = e^{-2t}(6 + 12t) -5 \cos(2t)$$

A continuación se muestra una gráfica con el movimiento descrito por el objeto.

Función de movimiento del objeto.

Inicialmente el resorte sufre un estiramiento muy grande generando un movimiento transitorio y procede a amortiguarse hasta llegar al equilibrio entre la fuerza externa y la fuerza amortiguadora describiendo un movimiento estable.

La parte $e^{-2t}(6 + 12t)$ representa el movimiento transitorio, mientras que $-5 \cos(2t)$ representa el movimiento estable.

También se puede observar que las amplitudes (el estiramiento del resorte) son bastante grandes comparado con las dos situaciones anteriores.

$\square$

Resonancia

Resolvamos un problema de valores iniciales que nos permitirá definir el concepto de resonancia.

Resolver la ecuación

$$\dfrac{d^{2}x}{dt^{2}} + \omega^{2}x = g_{0} \sin(\lambda t) \label{40} \tag{40}$$

donde $g_{0}$ y $\lambda \neq \omega$ son constantes y los valores iniciales son $x(0) = 0$ y $x^{\prime}(0) = 0$.

Solución: Resolviendo la ecuación homogénea puedes verificar que la solución complementaria es

$$x_{c}(t) = c_{1} \cos(\omega t) + c_{2} \sin(\omega t) \label{41} \tag{41}$$

Para obtener una solución particular proponemos una solución de la forma

$$x_{p}(t) = A \cos(\lambda t) + B \sin(\lambda t)$$

y aplicamos el método de coeficientes indeterminados. Vemos que

$$\dfrac{dx_{p}}{dt} = -A \lambda \sin(\lambda t) + B \lambda \cos(\lambda t) \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}x_{p}}{dt^{2}} = -A \lambda^{2} \cos(\lambda t) -B \lambda^{2} \sin(\lambda t)$$

Sustituyamos en la ecuación diferencial.

\begin{align*}
\left[ -A \lambda^{2} \cos(\lambda t) -B \lambda^{2} \sin(\lambda t) \right] &+ \omega^{2} \left[A \cos(\lambda t) + B \sin(\lambda t) \right] \\
&= A (\omega^{2} -\lambda^{2}) \cos(\lambda t) + B (\omega^{2} -\lambda^{2}) \sin(\lambda t) \\
&= g_{0} \sin(\lambda t)
\end{align*}

Es claro que para que se cumpla la igualdad se debe satisfacer que

$$A = 0 \hspace{1cm} y \hspace{1cm} B = \dfrac{g_{0}}{\omega^{2} -\lambda^{2}}$$

Por tanto, la solución particular es

$$x_{p}(t) = \dfrac{g_{0}}{\omega^{2} -\lambda^{2}} \sin(\lambda t) \label{42} \tag{42}$$

Y la solución general es

$$x(t) = c_{1} \cos(\omega t) + c_{2} \sin(\omega t) + \dfrac{g_{0}}{\omega^{2} -\lambda^{2}} \sin(\lambda t) \label{43} \tag{43}$$

Verifica que aplicando las condiciones iniciales se obtiene que

$c_{1} = 0 \hspace{1cm} y \hspace{1cm} c_{2} = -\dfrac{\lambda g_{0}}{\omega (\omega^{2} -\lambda^{2})}$

Por lo tanto, para $\lambda \neq \omega$, la solución es

$$x(t) = \dfrac{g_{0}}{\omega (\omega^{2} -\lambda^{2})} \left[-\lambda \sin(\omega t) + \omega \sin(\lambda t) \right] \label{44} \tag{44}$$

Este resultado no esta definido para $\lambda = \omega$, sin embargo podemos obtener su valor límite conforme $\lambda \rightarrow \omega$, esto produciría en (\ref{44}) un incremento de forma sustancial de las amplitudes de oscilación.

Para $\lambda = \omega$ se define la solución como el límite $\lambda \to \omega$ de la ecuación (\ref{44}).

$$x(t) = \lim_{\lambda \to \omega} g_{0} \dfrac{-\lambda \sin(\omega t) + \omega \sin(\lambda t)}{\omega (\omega^{2} -\lambda^{2})} \label{45} \tag{45}$$

Para resolver el límite apliquemos la regla de L´Hôpital.

La derivada del numerador con respecto a $\lambda$ es

$$\dfrac{d}{d \lambda} \left[ -\lambda \sin(\omega t) + \omega \sin(\lambda t) \right] = -\sin(\omega t) + t \omega \cos(\lambda t)$$

Y la derivada del denominador con respecto a $\lambda$ es

$$\dfrac{d}{d \lambda} \left[ \omega^{3} -\omega \lambda^{2} \right] = -2 \omega \lambda$$

Sustituyendo en el límite (\ref{45}) obtenemos

\begin{align*}
x(t) &= g_{0} \lim_{\lambda \to \omega} \dfrac{-\sin(\omega t) + t \omega \cos(\lambda t)}{-2 \omega \lambda} \\
&= g_{0} \dfrac{-\sin(\omega t) + t \omega \cos(\omega t)}{-2 \omega^{2}} \\
&= \dfrac{g_{0}}{2 \omega^{2}} \sin(\omega t) -\dfrac{g_{0}}{2 \omega}t \cos(\omega t)
\end{align*}

Por lo tanto, para $\lambda = \omega$ la solución es

$$x(t) = \dfrac{g_{0}}{2 \omega^{2}} \left[ \sin(\omega t) -t \omega \cos(\omega t) \right] \label{46} \tag{46}$$


Conforme $t \to \infty$ los desplazamientos del objeto se vuelven más largos, de hecho, $|x(t_{n})| \to \infty$ cuando $t_{n} = \dfrac{n \pi}{\omega}$ para $n = 1, 2, 3, \cdots$.

Este fenómeno se conoce como resonancia pura.

Una gráfica que muestra el comportamiento de (\ref{46}) es la siguiente.

Resonancia pura.

No profundizaremos más en el concepto de resonancia, pero cabe mencionar que la resonancia pura es una situación ideal, pues físicamente las oscilaciones grandes del objeto forzarían en algún momento al resorte más allá de su límite elástico, además en el desarrollo realizado no se han toman en cuenta efectos retardadores de las fuerzas de amortiguamiento que siempre están presentes.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

Resolver los siguientes problemas:

  1. Un resorte cuelga verticalmente; su extremo superior está fijo y del inferior pende una caja que pesa $196 N$. Una vez en equilibrio se tira de la caja hacia abajo haciéndola desplazar $0.25 m$ y se suelta. Sabiendo que $k = 80 N/m$ y que la resistencia del aire es despreciable, hallar:
  • La ecuación de movimiento de la caja.
  • El tiempo necesario para que la caja se mueva desde la posición inicial hasta $0.0625 m$ por debajo de la posición de equilibrio.
  1. Una masa de $98 N$ de peso se cuelga de un resorte con lo que éste interrumpe su estado de reposo. Sabiendo que $k = 4.9 N/m$, hallar el movimiento de la masa si al soporte del resorte se le imprime una fuerza $F_{ext}(t) = \sin(\sqrt{2g}t)$.
  1. Se suspende una masa de $10 kg$ de un resorte, el cual se alarga $0.6533 m$. La masa se pone en movimiento desde la posición de equilibrio con una velocidad inicial $1 m/s$, dirigida hacia arriba. Hallar el movimiento resultante si la fuerza debida al aire es $F_{am} = 80 \dfrac{dx}{dt}$.
  1. De un resorte que tiene una constante $k = 50$ se suspende un peso de $49 N$. El peso se pone en movimiento desde el reposo estirándolo $0.98 m$ hacia abajo de la posición de equilibrio y aplicando una fuerza externa $F_{ext}(t) = 10\sin(2t)$. Si no hay resistencia del aire, hallar el movimiento del peso.
  1. Se cuelga de un resorte una masa de $2 kg$, de tal manera que el resorte se alarga $0.6125 m$. A esta masa se le aleja de su posición de equilibrio jalándola $1 m$ hacia arriba y se suelta. Hallar el movimiento resultante de la masa sabiendo que hay una resistencia del aire de $F_{am} = 16 \dfrac{dx}{dt}$.

Más adelante…

Es momento de estudiar el caso en el que los coeficientes no son constantes, es decir, estudiaremos ecuaciones diferenciales de la forma

$$a(x) \dfrac{d^{2}y}{dx^{2}} + b(x) \dfrac{dy}{dx} + c(x) = g(x)$$

Donde $a(x)$, $b(x)$ y $c(x)$ son funciones de la variable independiente $x$.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»