Entrar Página Discusión Historial Go to the site toolbox

Estudio analítico de dos masas unidas por un muelle

De Laplace

Contenido

1 Enunciado

Como en el problema “Dos masas unidas por un muelle” tenemos dos masas m1 y m2 se mueven a lo largo del eje OX unidas por un resorte de constante k y longitud natura \ell_0. Inicialmente las dos masas se encuentran en reposo en x10 = 0 y x_{20}=\ell_0. Entonces se le comunica a la masa m1 una velocidad v0 en el sentido positivo del eje.

  1. Determine la lagrangiana del sistema en función de las posiciones de las dos partículas.
  2. Obtenga las ecuaciones de movimiento para x1 y x2

Realice el cambio de variables a las coordenadas generalizadas xG = (m1x1 + m2x2) / (m1 + m2), x=x_2-x_1-\ell_0.

  1. ¿Cómo queda la lagrangiana en función de estas coordenadas?
  2. Obtenga las ecuaciones de movimiento para xG y x.
  3. Determine dos constantes de movimiento para este sistema.

2 Lagrangiana en función de x1 y x2

Energía cinética

Las dos partículas describen un movimiento rectilíneo, por lo que la energía cinética es simplemente

T=\frac{1}{2}m_1\dot{x}_1^2+\frac{1}{2}m_2\dot{x}_2^2
Energía potencial

Esta es la elástica debida a la elongación del resorte

U=\frac{1}{2}k(x_2-x_1-\ell_0)^2
Lagrangiana

Restamos las dos cantidades anteriores

\mathcal{L}=T-U=\frac{1}{2}m_1\dot{x}_1^2+\frac{1}{2}m_2\dot{x}_2^2-\frac{1}{2}k(x_2-x_1-\ell_0)^2

3 Ecuaciones de movimiento para x1 y x2

Aplicamos la ecuación de Lagrange a cada una de las coordenadas

Para la coordenada x1
Tenemos que
p_1=\frac{\partial\mathcal{L}}{\partial \dot{x}_1}=m_1\dot{x}_1\qquad\qquad
\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial \dot{x}_1}\right)=m_1\ddot{x}_1
y
\frac{\partial\mathcal{L}}{\partial {x}_1}=k(x_2-x_1-\ell_0)
lo que da
m_1\ddot{x}_1=k(x_2-x_1-\ell_0)
Para la coordenada x2
De la misma manera
p_2=\frac{\partial\mathcal{L}}{\partial \dot{x}_2}=m_2\dot{x}_2\qquad\qquad
\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial \dot{x}_2}\right)=m_2\ddot{x}_2
y
\frac{\partial\mathcal{L}}{\partial {x}_2}=-k(x_2-x_1-\ell_0)
lo que nos da
m_2\ddot{x}_2=-k(x_2-x_1-\ell_0)

Estas son las mismas ecuaciones a las que se llega de forma simple partiendo de la segunda ley de Newton.

4 Lagrangiana en función de x y xG

Ninguna de las dos coordenadas es cíclica, ya que ambas aparecen en la energía potencial. Esto produce como resultado dos ecuaciones diferenciales acopladas.

Podemos separarlas mediante el cambio de variables sugerido. Si hacemos

x = x_2-x_1-\ell_0\qquad\qquad x_G=\frac{m_1x_1+m_2x_2}{m_1+m_2}

las coordenadas originales se escriben en función de las nuevas como

x_1=x_G-\frac{m_2(x+\ell_0)}{m_1+m2}\qquad x_2=x_G+\frac{m_1(x+\ell_0)}{m_2+m_2}

y las velocidades

\dot{x}_1=\dot{x}_G-\frac{m_2}{m_1+m2}\dot{x}\qquad \dot{x}_2=\dot{x}_G+\frac{m_1}{m_2+m_2}\dot{x}

Llevando esto a la energía cinética, tras operar nos queda

T=\frac{1}{2}(m_1+m_2)\dot{x}_G^2+\frac{m_1m_2}{2(m_1+m_2)}\dot{x}^2

que podemos leer como un caso particular del teorema de König para la energía cinética.

En cuanto a la potencial, su nueva expresión es

U=\frac{1}{2}kx^2

lo que nos da la lagrangiana

\mathcal{L}=\frac{1}{2}(m_1+m_2)\dot{x}_G^2+\frac{m_1m_2}{2(m_1+m_2)}\dot{x}^2-\frac{1}{2}kx^2

5 Ecuaciones de movimiento para x y xG

Para la coordenada xG
p_G=\frac{\partial\mathcal{L}}{\partial \dot{x}_G}=(m_1+m_2)\dot{x}_G\qquad\qquad
\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial \dot{x}_G}\right)=(m_1+m_2)\ddot{x}_G
y
\frac{\partial\mathcal{L}}{\partial {x}_G}=0
lo que da
(m_1+m_2)\ddot{x}_G=0
Para la coordenada x
De la misma manera
p=\frac{\partial\mathcal{L}}{\partial \dot{x}}=\frac{m_1m_2}{m_1+m_2}\dot{x}\qquad\qquad
\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial \dot{x}}\right)=\frac{m_1m_2}{m_1+m_2}\ddot{x}
y
\frac{\partial\mathcal{L}}{\partial {x}}=-kx
lo que nos da
\frac{m_1m_2}{m_1+m_2}\ddot{x}=-kx

Ahora sí resultan ecuaciones desacopladas, cada una de las cuales puede integrarse por separado. La primera es trivial

\ddot{x}_G=0\qquad\Rightarrow\qquad x_G=x_{G0}+v_{G0} t

y la segunda es la de un oscilador armónico

\ddot{x}=-\omega^2 x\qquad\qquad \omega=\sqrt{\frac{k(m_1+m_2)}{m_1m_2}}

con solución

x=c_1\cos(\omega t)+c_2\,\mathrm{sen}(\omega t)

Los valores de las constantes salen de las condiciones iniciales. Para la posición del CM

No se pudo entender (Falta el ejecutable de <strong>texvc</strong>. Por favor, lea <em>math/README</em> para configurarlo.): x_G(t=0)=\frac{m_1x_{10}+m_2x_{20}}{m_1+m_2}=\frac{m_2\ell_0}{m_1+m_2}\qquad\qquad \dot{x}_{G}(t=0)=\frac{m_1\dot{x}_{10}+m_2\dot{x}_{20}}{m_1+m_2}=\frac{m_2v_0}}{m_1+m_2}

lo que nos da, para todo instante

x_G=\frac{m_2}{m_1+m_2}(\ell_0+v_0t)

Para la elongación del muelle

x(t=0)=x_{20}-x_{10}-\ell_0=0\qquad\qquad \dot{x}(t=0)=\dot{x}_{20}-\dot{x}_{10}=v_0

siendo por tanto la elongación para todo instante

x=\frac{v_0}{\omega}\,\mathrm{sen}(\omega t)

Conocidas las dos coordenadas obtenemos las posiciones de las dos masas

x_1=\frac{m_2}{m_1+m_2}(\ell_0+v_0t)-\frac{m_2}{m_1+m_2}\left(\ell_0+\frac{v_0}{\omega}\,\mathrm{sen}(\omega t)\right)=\frac{m_2}{m_1+m_2}\left(v_0t-\frac{v_0}{\omega}\,\mathrm{sen}(\omega t)\right)

 

x_2=\frac{m_2}{m_1+m_2}(\ell_0+v_0t)+\frac{m_1}{m_1+m_2}\left(\ell_0+\frac{v_0}{\omega}\,\mathrm{sen}(\omega t)\right)=\ell_0+\frac{v_0}{\omega(m_1+m_2)}\left(m_2\omega t+m_1\,\mathrm{sen}(\omega t)\right)

6 Constantes de movimiento

Empleando las nuevas coordenadas, vemos que xG es cíclica, esto es, no aparece en la lagrangiana. Por tanto su momento conjugado es una constante de movimiento

p_G=(m_1+m_2)\dot{x}_G=m_1\dot{x}_1+m_2\dot{x}_2=\mathrm{cte.}

Esta magnitud no es otra que la cantidad de movimiento del sistema.

La otra constante procede de que la lagrangiana no depende del tiempo. Como además la energía cinética es una función cuadrática de las velocidades, la cantidad conservada es la energía mecánica

E=\frac{1}{2}(m_1+m_2)\dot{x}_G^2+\frac{m_1m_2}{2(m_1+m_2)}\dot{x}^2+\frac{1}{2}kx^2=\frac{1}{2}m_1\dot{x}_1^2+\frac{1}{2}m_2\dot{x}_2^2+\frac{1}{2}k(x_2-x_1-\ell_0)^2

Herramientas:

Herramientas personales
TOOLBOX
LANGUAGES
licencia de Creative Commons
Aviso legal - Acerca de Laplace