Ver el código de «Masa en superficie horizontal con masa colgando verticalmente»

= Enunciado =
[[Imagen:F1GIC_masaHorizontalMasaVerticalMuelle_enunciado.png|right|250px]]
Dos masas puntuales <math>m_1</math> y <math>m_2</math> están unidas por una cuerda sin masa y longitud <math>L</math>, que desliza sobre una polea también sin masa, como se indica en la figura. La masa <math>m_1</math> está conectada a un muelle de constante elástica <math>k</math> y longitud natural nula. El otro extremo del muelle está anclado en el origen del sistema de referencia <math>OXY</math>. 
#Suponiendo que no hay rozamiento, encuentra la posición de equilibrio del sistema.
#Si hay rozamiento entre <math>m_1</math> y la superfice horizontal, caracterizado por un coeficiente de rozamiento estático <math>\mu</math>, encuentra el rango de posibles posiciones de equilibrio. Suponiendo <math>m_1= kL/4g</math> y <math>m_2= kL/2g</math> ¿que condición debe cumplir <math>\mu</math> para que toda la superfice sea una posible posición de equilibrio?
#Encuentra la ecuación de movimiento de las dos masas suponiendo que no hay rozamiento.
#Calcula el vector de posición de las dos masas en función del tiempo si las condiciones iniciales son <math>x_1(0)=m_2g/k</math> y <math>\dot{x}_1(0)=v_0</math>.

= Solución =

== Posición de equilibrio sin rozamiento ==
[[Imagen:F1GIC_masaHorizontalMasaVerticalMuelle_fuerzasSinRozamiento.png|right|300px]]
La figura de la derecha muestra las fuerzas que actúan sobre cada una de las masas. Para la masa 1 son: su peso <math>\vec{G}_1</math>, la fuerza del muelle <math>\vec{F}_k</math>, la de la cuerda <math>\vec{T}_1</math> y la fuerza vincular ejercida por la superficie horizontal <math>\vec{N}_1</math>.
Para la masa 2 son su peso <math>\vec{G}_2</math> y la fuerza de la cuerda <math>\vec{T}_2</math>. Usando la base cartesiana asociada a los ejes de la figura estas fuerzas pueden expresarse así

Masa 1
<center>
<math>
\begin{array}{l}
\vec{G}_1 = -m_1 g \vec{\jmath}, \\
\vec{F}_k = -k\,\overrightarrow{OP}_1 = -kx_1\,\vec{\imath},\\
\vec{N}_1 = N\,\vec{\jmath},\\
\vec{T}_1 = T\,\vec{\imath}.
\end{array}
</math>
</center>
Masa 2
<center>
<math>
\begin{array}{l}
\vec{G}_2 = -m_2 g \vec{\jmath}, \\
\vec{T}_2 = T\,\vec{\jmath}.
\end{array}
</math>
</center>
Hemos usado que, al ser la polea y la cuerda sin masa, la tensión se transmite a lo largo de la cuerda. 

La condición de equilibrio es que la suma de las fuerzas sobre cada masa sea nula. Esto nos da 4 ecuaciones, pues el problema es bidimensional. Tenemos
<center>
<math>
\vec{G}_1 + \vec{F}_k + \vec{N}_1 + \vec{T}_1 = \vec{0}
\longrightarrow
\left\{
\begin{array}{lcll}
X) & \to & -kx_1 + T = 0, & (1)\\
Y) & \to & N-m_1g = 0. & (2)
\end{array}
\right.
</math>
</center>
<center>
<math>
\vec{G}_2 + \vec{T}_2 = \vec{0}
\longrightarrow
\left\{
\begin{array}{lcll}
X) & \to & 0 = 0, & (3)\\
Y) & \to & T-m_2g = 0. & (4)
\end{array}
\right.
</math>
</center>
En este caso la ecuación (3) es superflua porque la masa 2 sólo tiene movimiento vertical.

Las incógnitas son <math>\{x_1, N, T\}</math> y tenemos 3 ecuaciones. El problema tiene solución única y es
<center>
<math>
x_1^{eq} = m_2 g/k, \qquad N = m_1g, \qquad T = m_2g.
</math>
</center>

== Posiciones de equilibrio con rozamiento ==
[[Imagen:F1GIC_masaHorizontalMasaVerticalMuelle_fuerzasConRozamiento.png|right|300px]]
El diagrama de fuerzas de la derecha es similar al del apartado anterior. La diferencia es que hay que añadir una fuerza más sobre la masa 1, la de rozamiento que ejerce la superficie horizontal. Sólo conocemos a priori su dirección no el sentido. Tampoco conocemos su magnitude. Esta fuerza se expresa 
<center>
<math>
\vec{F}_F = f\,\vec{\imath}.
</math>
</center>
Las condiciones de equilibrio y las ecuaciones correspondiente son
<center>
<math>
\vec{G}_1 + \vec{F}_k + \vec{N}_1 + \vec{T}_1 + \vec{F}_R = \vec{0}
\longrightarrow
\left\{
\begin{array}{lcll}
X) & \to & -kx_1 + T + f = 0, & (5)\\
Y) & \to & N-m_1g = 0. & (6)
\end{array}
\right.
</math>
</center>
<center>
<math>
\vec{G}_2 + \vec{T}_2 = \vec{0}
\longrightarrow
\left\{
\begin{array}{lcll}
X) & \to & 0 = 0, & (7)\\
Y) & \to & T-m_2g = 0. & (8)
\end{array}
\right.
</math>
</center>
Ahora las magnitudes desconocidas son <math>\{x_1, T, f, N\}</math>. Son 4. Y sólo tenemos 3 ecuaciones.

Lo que ocurre es que ahora la solución no es única. El rozamiento hace que haya varios valores posibles de <math>x_1</math> para los que las masas están en equilibrio. ¿Cómo afrontamos esta situación? Vamos a suponer que el valor de <math>x_1</math> es un dato. Vamos a llamarlo <math>x_1 = d</math>. Si ahora <math>d</math> es un dato, las incógnitas son <math>\{T, f, N\}</math>. Podemos entonces calcular su valor con las ecuaciones (5), (6) y (8).
<center>
<math>
f = kd - m_2g, \qquad N=m_1g, \qquad T = m_2g.
</math>
</center>
Con estos resultados las fuerzas vinculares que actúan sobre las masas, (usando los valores de <math>m_1</math> y <math>m_2</math> que se dan para este apartado)
<center>
<math>
\vec{F}_R = ( kd - m_2g)\,\vec{\imath}, \qquad \vec{N}_1 = m_1g\,\vec{\jmath}, \qquad \vec{T}_1 = m_2g\,\vec{\imath}, \qquad \vec{T}_2=m_2g\,\vec{\jmath}.
</math>
</center>
Esos son los valores de las fuerzas necesarios para que haya equilibrio. Ahora bien, esas cuatro fuerzas no se comportan igual. Las fuerzas <math>\vec{N}_1</math>, <math>\vec{T}_1</math> y <math>\vec{T}_2</math> pueden ser tan grandes como haga falta, pues suponemos que las ligaduras a las que están asociadas se mantienen siempre. Pero la fuerza de rozamiento <math>\vec{F}_R</math> no se comporta igual. Esta fuerza tiene un módulo máximo que vale
<center>
<math>
|\vec{F}_R^{max}| = \mu |\vec{N}_1| = \mu m_1g
</math>
</center>
Es decir, el valor de la fuerza de rozamiento tiene que cumplir la condición
<center>
<math>
|kd - m_2g| \leq \mu m_1g.
</math>
</center>
Utilizando los valores de <math>m_1</math> y <math>m_2</math> que se dan para este apartado esta condición queda
<center>
<math>
|2d - L| \leq \mu L/2. \qquad\qquad(9)
</math>
</center>
Debemos considerar dos situaciones posibles:

<math>2d>L</math>: En este caso se tiene que <math>(2d-L)>0</math>, y por tanto <math>|2d-L| = 2d-L</math>. Usándolo en (9) obtenemos
<center>
<math>
2d-L \leq \mu \dfrac{L}{2} \longrightarrow d\leq (2+\mu)\,\dfrac{L}{4}=x_{max}.
</math>
</center>

<math>2d<L</math>: En este caso se tiene que <math>(2d-L)<0</math>, y por tanto <math>|2d-L| = -(2d-L)=L-2d</math>. Usándolo en (9) obtenemos
<center>
<math>
L-2d \leq \mu \dfrac{L}{2} \longrightarrow d\geq (2-\mu)\,\dfrac{L}{4}=x_{min}.
</math>
</center>
Es decir, para que la fuerza de rozamiento sea capaz de mantener el equilibrio deben cumplirse estas dos condiciones a la vez, lo que nos dice que el valor de <math>d</math> (que da el valor posible de equilibrio de <math>x</math>) debe cumplir
<center>
<math>
d \in [x_{min}, x_{max}] 
</math>
</center>
¿Que pasa si <math>2d=L</math>?. Esta posición corresponde al valor de <math>x_1^{eq}=L/2</math> con los valores de las masas usados en este apartado. En este caso, la fuerza de rozamiento es cero pues no hace falta para mantener el equilibrio. Podemos entender ahora que significan las dos situaciones que hemos considerado. Si <math>d<L/2</math> la masa 1 estaría a la izquierda de la posición de equilibrio sin rozamiento. Entonces la fuerza del muelle sería menor que la fuerza que ejerce la cuerda y la fuerza de rozamiento apuntaría hacia la izquierda para compensar. Si <math>d>L/2</math> la masa 1 estaría a la derecha de la posición de equilibrio sin rozamiento y la fuerza ejercida por el muelle sería mayor que la de la cuerda y la fuerza de rozamiento apuntaría hacia la derecha para compensar.

=== Condición para que toda la superficie horizontal sea de equilibrio ===
Para que esto ocurra debe cumplirse
<center>
<math>
\begin{array}{l}
x_{min} \leq 0 \longrightarrow 2-\mu \leq 0 \longrightarrow \mu\geq 2, \\
x_{max} \geq L \longrightarrow 2+\mu \geq 4 \longrightarrow \mu\geq 2.
\end{array}
</math>
</center>
Por tanto, la condición sobre el coeficiente de rozamiento es
<center>
<math>
\mu \geq 2.
</math>
</center>

== Movimiento sin rozamiento ==
[[Imagen:F1GIC_masaHorizontalMasaVerticalMuelle_distancias.png|right|300px]]
Si no hay rozamiento el diagrama de fuerzas del primer apartado vuelve a ser válido. Lo que cambia ahora es que la ley física que tenemos que aplicar es la Segunda Ley de Newton. Hay que hacerlo para las dos masas por separado
<center>
<math>
\begin{array}{l}
m_1\vec{a}_1 = \vec{F}_{g1} + \vec{N}_1 + \vec{F}_k  + \vec{T}_1,\\
m_2\vec{a}_2 = \vec{F}_{g2} + \vec{T}_2.
\end{array}
</math>
</center>
Tenemos que encontrar expresiones para las aceleraciones de las masas. Del dibujo de la derecha vemos que 
<center>
<math>
\begin{array}{l}
\overrightarrow{OP}_1 = x_1\,\vec{\imath},\\
\overrightarrow{OP}_2 = L\,\vec{\imath} - h\,\vec{\jmath}.
\end{array}
</math>
</center>
El enunciado dice que la longitud de la cuerda es <math>L</math>. Entonces se cumple
<center>
<math>
(L-x_1) + h = L \longrightarrow h = x_1.
</math>
</center>
Entonces tenemos
<center>
<math>
\begin{array}{l}
\overrightarrow{OP}_1 = x_1\,\vec{\imath},\\
\overrightarrow{OP}_2 = L\,\vec{\imath} - x_1\,\vec{\jmath}.
\end{array}
</math>
</center>
Derivamos una vez respecto al tiempo para obtener las velocidades
<center>
<math>
\begin{array}{l}
\vec{v}_1 = \dot{\overrightarrow{OP}}_1 = \dot{x}_1\,\vec{\imath},\\
\vec{v}_2 = \dot{\overrightarrow{OP}}_2 = - \dot{x}_1\,\vec{\jmath}.
\end{array}
</math>
</center>
Derivamos otra vez respecto del tiempo para obtener las aceleraciones
<center>
<math>
\begin{array}{l}
\vec{a}_1 = \dot{\vec{v}}_1 = \ddot{x}_1\,\vec{\imath},\\
\vec{a}_2 = \dot{\vec{v}}_2 = - \ddot{x}_1\,\vec{\jmath}.
\end{array}
</math>
</center>
A partir de la Segunda Ley de Newton obtenemos
<center>
<math>
\begin{array}{l}
m_1\vec{a}_1 = \vec{G}_1 + \vec{N}_1 + \vec{T}_1 + \vec{F}_k 
\left\{
\begin{array}{lclr}
X) & \to & m_1\ddot{x}_1 = T - kx_1, & (10)\\
Y) & \to & 0 = -m_1g + N_1. & (11)
\end{array}
\right.\\
m_2\vec{a}_2 = \vec{G}_2 + \vec{T}_2 \to
\left\{
\begin{array}{lclr}
X) & \to & 0 = 0 & (12)\\
Y) & \to & -m_2\ddot{x}_1 = -m_2g + T. & (13)
\end{array}
\right.
\end{array}
</math>
</center>
Sumando las ecuaciones (10) y (13) tenemos
<center>
<math>
(m_1+m_2) \ddot{x}_1 = -kx_1 + m_2g
\longrightarrow
\ddot{x}_1 = -\dfrac{k}{m_1+m_2}x_1 + \dfrac{m_2}{m_1+m_2}g. \qquad (14)
</math>
</center>
Esta es la ecuación de movimiento. Las ecuaciones (11) y (13) nos dan los valores de <math>N_1</math> y <math>T</math>:
<center>
<math>
N_1 = m_1g, \qquad T = m_2g - m\ddot{x}_1.
</math>
</center>
Es interesante destacar que <math>T</math> no es igual al peso de la masa 2. Esto sólo ocurre en situación de equilibrio estático, es decir, cuando <math>\ddot{x}_1=0</math>.

== Solución para las condiciones iniciales dadas ==
[[Imagen:F1GIC_masaHorizontalMasaVerticalMuelle_cambio.png|right|300px]]
Para poner la ecuación (14) en una forma mas reconocible, usamos la distancia <math>s</math> para describir la posición de la masa 1 (ver figura de la derecha)
<center>
<math>
x_1(t) = x_1^{eq} + s(t) = \dfrac{m_2g}{k} \longrightarrow
\dot{x}_1(t) = \dot{s}(t) \longrightarrow
\ddot{x}_1(t) = \ddot{s}(t).
</math>
</center>
Sustituyendo estas expresiones de <math>x_1</math> y <math>\ddot{x}_1</math> en la ecuación (14) obtenemos la ecuación diferencial
<center>
<math>
\ddot{s} = -\dfrac{k}{m_1+m_2}s
</math>
</center>
Esta ecuación tiene la forma
<center>
<math>
\ddot{s} = -\omega_0^2s, 
\qquad 
\omega_0 = \sqrt{\dfrac{k}{m_1+m_2}}.
</math>
</center>
Corresponde a un movimiento armónico simple con período
<center>
<math>
T = \dfrac{2\pi}{w_0} = 2\pi\sqrt{\dfrac{m_1+m_2}{k}}.
</math>
</center>
La solución general es de la forma
<center>
<math>
s(t) = A\cos(\omega_0t) + B\,\mathrm{sen}(\omega_0t).
</math>
</center>
Las constantes <math>A</math> y <math>B</math> hay que determinarlas a partir de las condiciones iniciales.
Las condiciones iniciales para <math>s(t)</math> son
<center>
<math>
x_1(0) = \dfrac{m_2g}{k} \longrightarrow s(0) = 0,
\qquad
\dot{x}_1(0) = v_0 \longrightarrow \dot{s}(0) = v_0.
</math>
</center>
Aplicando estas condiciones obtenemos
<center>
<math>
\begin{array}{l}
\left.
\begin{array}{l}
s(0) = 0\\
s(0) = A
\end{array}
\right\} \longrightarrow A = 0.\\
\left.
\begin{array}{l}
\dot{s}(0) = v_0\\
\dot{s}(0) = B\omega_0
\end{array}
\right\} \longrightarrow B = v_0/\omega_0.
\end{array}
</math>
</center>
Por tanto, el movimiento de la masa 1 queda descrito por 
<center>
<math>
s(t) = \dfrac{v_0}{\omega_0}\,\mathrm{sen}\,(\omega_0t)
\longrightarrow
x_1(t) = \dfrac{m_2g}{k} + \dfrac{v_0}{\omega_0}\,\mathrm{sen}\,(\omega_0t)
</math>
</center>
El movimiento de la masa 2 también queda determinado al conocer <math>x_1(t)</math>.

[[Categoría:Dinámica del punto material|1]]