Ver el código de «Masa conectada a un muelle horizontal (GIC)»

= Enunciado =
Una masa <math>m</math> puede deslizar sin rozamiento sobre el eje horizontal <math>OX</math>.  La masa está conectada a un muelle de constante elástica <math>k</math> y longitud  natural <math>l_0</math>. El otro extremo del muelle está conectado al origen de  coordenadas. Determina el movimiento <math>x(t)</math> de la masa en estos tres casos:
#Se suelta la masa desde el origen con velocidad nula.
#Se suelta la masa desde la posición de equilibrio con velocidad <math>v_0</math>.
#Se suelta la masa desde la posición <math>x_0</math> con velocidad <math>v_0</math>.

= Solución =

== Ecuación de movimiento ==
[[File:F1GIC_muelle_horizontal_fuerzas.png|right|250px]]
Las fuerzas que actúan sobre la masa son su peso <math>\vec{P}</math>, la fuerza ejercida por el muelle, <math>\vec{F}_k</math> y la fuerza vincular de la superficie horizontal de la superficie sobre la masa, <math>\vec{N}</math>. La ecuación diferencial del movimiento es proporcionada por la Segunda Ley de Newton
<center>
<math>
m\vec{a} = \vec{F}_k + \vec{P} + \vec{N}
</math>
</center>
Escogemos el eje <math>X</math> paralelo a la superficie y el eje <math>Y</math> perpendicular a ella. Entonces la expresión del término de inercia y las fuerzas es
<center>
<math>
\begin{array}{l}
m\vec{a} = m\ddot{x}\,\vec{\imath},\\
\\
\vec{F}_k = -k\overrightarrow{OA} = -k(x-l_0)\,\vec{\imath},\\
\\
\vec{P} = -mg\,\vec{\jmath},\\
\\
\vec{N} = N\,\vec{\jmath}.
\end{array}
</math>
</center>
Al aplicar la Segunda Ley obtenemos dos ecuaciones escalares, una por cada componente de la base vectorial
<center>
<math>
\begin{array}{lclr}
X) & \to & m\ddot{x} = -kx + kl_0, & (1)\\
&&&\\
Y) & \to & N-mg = 0. & (2)
\end{array}
</math>
</center>
La ecuación (2) nos dice que la fuerza vincular equilibra al peso para que la partícula no atraviese la superficie. Podemos entonces ignorar el peso y la fuerza normal. Si hubiera rozamiento no podríamos hacerlo, pues en régimen dinámico el módulo de la fuerza de rozamiento es proporcional a la normal y no podemos ignorar el valor de esta.

La ecuación (1) es una ecuación diferencial donde la incógnita es <math>x(t)</math>, la coordenada de la partícula sobre el eje <math>X</math>. Vamos a hacer un cambio de variable. Vamos la determinar la posición de la partícula respecto a la posición de equilibrio del muelle. Para ello definimos <math>s(t)</math> de modo que
[[File:F1GIC_muelle_horizontal_xs.png|right|250px]]
<center>
<math>
x(t) = l_0 + s(t). \qquad \qquad (3)
</math>
</center>
Cuando <math>s=0</math> el muelle está relajado y no ejerce fuerza sobre la partícula. Podemos derivar respecto al tiempo a ambos lados de la expresión (3). Teniendo en cuenta que <math>l_0</math> no depende del tiempo obtenemos
<center>
<math>
\dot{x} = \dot{s}, \qquad\qquad \ddot{x}=\ddot{s}. \qquad\qquad(4)
</math>
</center>
Sustituyendo las expresiones (3) y (4) en la ecuación (1) obtenemos una ecuación diferencial para <math>s(t)</math>
<center>
<math>
m\ddot{s} = -k\,s.
</math>
</center>
Vamos a pasar <math>m</math> al otro lado  
<center>
<math>
\ddot{s} = -\dfrac{k}{m}\,s.
</math>
</center>
Definimos <math>\omega_0^2=k/m</math>. La ecuación queda
<center>
<math>
\ddot{s} = -\omega_0^2s, \qquad\qquad \omega_0=\sqrt{k/m}.
</math>
</center>
Obtenemos una ecuación diferencial de segundo orden, porque aparece la derivada segunda de la función incógnita. Esta ecuación es especial y muy importante, por lo que recibe un nombre especial. Es la ecuación de un oscilador armónico, o ecuación del Movimiento Armónico Simple (MAS).

Para que una función verifique la ecuación debe ocurrir que al derivarla dos veces respecto al tiempo el resultado sea la misma función multiplicada por <math>\omega_0</math>. Hay dos funciones que cumplen esto: <math>\cos(\omega_0t)</math> y <math>\,\mathrm{sen}\,(\omega_0t)</math>. La solución general es una combinación lineal de las dos
<center>
<math>
s(t) = A\cos(\omega_0t) + B\,\mathrm{sen}\,(\omega_0t).
</math>
</center>
Vamos a comprobar que esta función verifica la ecuación diferencial. Derivando respecto al tiempo obtenemos
<center>
<math>
\dot{s}(t) = -A\omega_0\,\mathrm{sen}\,(\omega_0t) + B\omega_0\cos(\omega_0t).
</math>
</center>
Derivando otra vez
<center>
<math>
\ddot{s} = -A\omega_0^2\cos(\omega_0t) - B\omega_0^2\,\mathrm{sen}\,(\omega_0t) = -\omega_0^2\,(A\cos(\omega_0t) + B\,\mathrm{sen}\,(\omega_0t))
=
-\omega_0^2s(t).
</math>
</center>
Efectivamente esta función verifica la ecuación diferencial.

La solución general tiene dos constantes de integración <math>A</math> y <math>B</math>. Desde el punto de vista matemático esto es así porque la ecuación es de segundo orden. Desde el punto de vista físico esto significa que para determinar el movimiento de la partícula necesitamos dos condiciones iniciales: su posición y su velocidad en el instante inicial.

== Caso 1 ==
En este caso la posición y velocidad iniciales de la partícula son
<center>
<math>
x(0) = 0, \qquad v(0) = \dot{x}(0) = 0.
</math>
</center>
Usando las expresiones (3) y (4) estas condiciones se traducen en las siguientes para <math>s</math>
<center>
<math>
s(0) = -l_0, \qquad \dot{s}(0) = 0.
</math>
</center>
Con el fin de imponer estas condiciones sustituimos <math>t=0</math> en la solución general <math>s(t)</math> y su derivada <math>\dot{s}</math>
<center>
<math>
\begin{array}{lcl}
\left.\begin{array}{l}
s(0) = A\\
s(0) = -l_0
\end{array}
\right\}
& \to A = -l_0,\\
&&\\
\left.\begin{array}{l}
\dot{s}(0) = B\omega_0\\
\dot{s}(0) = 0
\end{array}
\right\}
& \to B = 0.\\
\end{array}
</math>
</center>
[[File:F1GIC_muelle_horizontal_xz.png|right|400px]]
Así pues, el movimiento de la partícula en el caso 1 está descrito por la función
<center>
<math>
s(t) = -l_0\cos(\omega_0t).
</math>
</center>
La coordenada <math>x</math> es
<center>
<math>
x(t) = l_0 + s(t) = l_0\,(1-\cos(\omega_0t)).
</math>
</center>
La figura de la derecha muestra las gráficas de <math>s(t)</math> y <math>x(t)</math>. La masa empieza su movimiento en el punto <math>O</math>. Se desplaza hacia la derecha, sobrepasa la posición de equilibrio <math>s=0</math>, hasta que llega hasta una distancia <math>l_0</math> a la derecha del punto de equilibrio <math>(s=l_0, \, x=2l_0)</math>. Ahí da la vuelta hasta que llega a <math>O</math> otra vez y se repite el ciclo. El movimiento es periódico. El período es el tiempo que la partícula tarda en ir y volver al punto <math>O</math>, es decir, el tiempo <math>T</math> para el cual el argumento del coseno es <math>2\pi</math>. Entonces
<center>
<math>
\omega_0 T = 2\pi 
\Longrightarrow
T = \dfrac{2\pi}{\omega_0} = 2\pi\sqrt{\dfrac{m}{k}}.
</math>
</center>
Ese es el período de oscilación de un muelle ideal.


[[File:F1GIC_muelle_horizontal_va.png|right|400px]]
La figura de la derecha muestra la evolución en el tiempo de la velocidad y de la aceleración de la partícula
<center>
<math>
\begin{array}{l}
v(t) = \dot{x} = \dot{s} = l_0\omega_0\,\mathrm{sen}\,(\omega_0 t),\\
\\
a(t) = \ddot{x} = \ddot{s} = l_0\omega_0^2\cos(\omega_0 t).
\end{array}
</math>
</center>
La rapidez es cero en los puntos extremos del movimiento y es máxima cuando la partícula pasa por el punto de equilibrio del muelle. Por el contrario, la aceleración es máxima en los puntos extremos del movimiento y es cero cuando la posición de la partícula coincide con la de punto de relajación del muelle.

== Caso 2 ==
Tenemos que calcular los valores de las constantes <math>A</math> y <math>B</math>  que describen esta situación.  Ahora las condiciones iniciales son 
<center>
<math>
x(0) = l_0, \qquad \qquad v(0) =  \dot{x}(0) = v_0.
</math>
</center>
Esto se traduce en estas condiciones iniciales para <math>s</math>
<center>
<math>
s(0) = 0, \qquad \qquad  \dot{s}(0) = v_0.
</math>
</center>
Sustituimos <math>t=0</math> en la solución general <math>s(t)</math> y su derivada <math>\dot{s}</math> e imponemos estas condiciones
<center>
<math>
\begin{array}{lcl}
\left.\begin{array}{l}
s(0) = A\\
s(0) = 0
\end{array}
\right\}
& \to A = 0,\\
&&\\
\left.\begin{array}{l}
\dot{s}(0) = B\omega_0\\
\dot{s}(0) = v_0
\end{array}
\right\}
& \to B = v_0/\omega_0.\\
\end{array}
</math>
</center>
Entonces la posición de la partícula es descrita por la función
<center>
<math>
s(t) = \dfrac{v_0}{\omega_0}\,\mathrm{sen}\,(\omega_0 t),
\qquad\qquad
x(t) = l_0 + s(t) = l_0 + \dfrac{v_0}{\omega_0}\,\mathrm{sen}\,(\omega_0 t).
</math>
</center>

== Caso 3 ==
Ahora las condiciones iniciales son 
<center>
<math>
x(0) = (x_0+l_0), \qquad \qquad v(0) =  \dot{x}(0) = v_0.
</math>
</center>
Esto se traduce en estas condiciones iniciales para <math>s</math>
<center>
<math>
s(0) = x_0, \qquad \qquad  \dot{s}(0) = v_0.
</math>
</center>
Sustituimos <math>t=0</math> en la solución general <math>s(t)</math> y su derivada <math>\dot{s}</math> e imponemos estas condiciones
<center>
<math>
\begin{array}{lcl}
\left.\begin{array}{l}
s(0) = A\\
s(0) = x_0
\end{array}
\right\}
& \to A = x_0,\\
&&\\
\left.\begin{array}{l}
\dot{s}(0) = B\omega_0\\
\dot{s}(0) = v_0
\end{array}
\right\}
& \to B = v_0/\omega_0.\\
\end{array}
</math>
</center>
Entonces la posición de la partícula es descrita por la función
<center>
<math>
s(t) = x_0\cos(\omega_0 t) + \dfrac{v_0}{\omega_0}\,\mathrm{sen}\,(\omega_0 t),
\qquad\qquad
x(t) = l_0 + s(t) = l_0 + x_0\cos(\omega_0 t) +  \dfrac{v_0}{\omega_0}\,\mathrm{sen}\,(\omega_0 t).
</math>
</center>

En los tres casos la partícula realiza oscilaciones periódicas con período
<center>
<math>
T = \dfrac{2\pi}{\omega_0} = 2\pi\sqrt{\dfrac{m}{k}}.
</math>
</center>
[[Categoría:Dinámica del punto material|1]]