Ver el código de «Estudio analítico de dos masas unidas por un muelle»

==Enunciado==
Como en el problema &ldquo;[[Dos masas unidas por un muelle]]&rdquo; tenemos dos masas <math>m_1</math> y <math>m_2</math> se mueven a lo largo del eje OX unidas por un resorte de constante <math>k</math> y longitud natura <math>\ell_0</math>. Inicialmente las dos masas se encuentran en reposo en <math>x_{10}=0</math> y <math>x_{20}=\ell_0</math>. Entonces se le comunica a la masa <math>m_1</math> una velocidad <math>v_0</math> en el sentido positivo del eje. 
<center>[[Archivo:Dos-masas-muelle-horizontal.png|400px]]</center>
# Determine la lagrangiana del sistema en función de las posiciones de las dos partículas.
# Obtenga las ecuaciones de movimiento para <math>x_1</math> y <math>x_2</math>
Realice el cambio de variables a las coordenadas generalizadas <math>x_G=(m_1 x_1+m_2 x_2 )/(m_1+m_2 )</math>, <math>x=x_2-x_1-\ell_0</math>.
<ol start="3">
<li>¿Cómo queda la lagrangiana en función de estas coordenadas?</li>
<li>Obtenga las ecuaciones de movimiento para <math>x_G</math> y <math>x</math>.</li>
<li>Determine dos constantes de movimiento para este sistema.</li>
</ol>

==Lagrangiana en función de x<sub>1</sub> y x<sub>2</sub>==
;Energía cinética:
Las dos partículas describen un movimiento rectilíneo, por lo que la energía cinética es simplemente

<center><math>T=\frac{1}{2}m_1\dot{x}_1^2+\frac{1}{2}m_2\dot{x}_2^2</math></center>

;Energía potencial:
Esta es la elástica debida a la elongación del resorte

<center><math>U=\frac{1}{2}k(x_2-x_1-\ell_0)^2</math></center>

;Lagrangiana:
Restamos las dos cantidades anteriores

<center><math>\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</math></center>
==Ecuaciones de movimiento para x<sub>1</sub> y x<sub>2</sub>==
Aplicamos la ecuación de Lagrange a cada una de las coordenadas

;Para la coordenada x<sub>1</sub>:Tenemos que

<center><math>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</math></center>

:y

<center><math>\frac{\partial\mathcal{L}}{\partial {x}_1}=k(x_2-x_1-\ell_0)</math></center>

:lo que da

<center><math>m_1\ddot{x}_1=k(x_2-x_1-\ell_0)</math></center>

;Para la coordenada x<sub>2</sub>:De la misma manera

<center><math>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</math></center>

:y

<center><math>\frac{\partial\mathcal{L}}{\partial {x}_2}=-k(x_2-x_1-\ell_0)</math></center>

:lo que nos da

<center><math>m_2\ddot{x}_2=-k(x_2-x_1-\ell_0)</math></center>

Estas son las mismas ecuaciones a las que se llega de forma simple partiendo de la segunda ley de Newton.
==Lagrangiana en función de x y x<sub>G</sub>==
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

<center><math>x = x_2-x_1-\ell_0\qquad\qquad x_G=\frac{m_1x_1+m_2x_2}{m_1+m_2}</math></center>

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

<center><math>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}</math></center>

y las velocidades

<center><math>\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}</math></center>

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

<center><math>T=\frac{1}{2}(m_1+m_2)\dot{x}_G^2+\frac{m_1m_2}{2(m_1+m_2)}\dot{x}^2</math></center>

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

<center><math>U=\frac{1}{2}kx^2</math></center>

lo que nos da la lagrangiana

<center><math>\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</math></center>
==Ecuaciones de movimiento para x y x<sub>G</sub>==
;Para la coordenada x<sub>G</sub>:
<center><math>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</math></center>

:y

<center><math>\frac{\partial\mathcal{L}}{\partial {x}_G}=0</math></center>

:lo que da

<center><math>(m_1+m_2)\ddot{x}_G=0</math></center>

;Para la coordenada x:De la misma manera

<center><math>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}</math></center>

:y

<center><math>\frac{\partial\mathcal{L}}{\partial {x}}=-kx</math></center>

:lo que nos da

<center><math>\frac{m_1m_2}{m_1+m_2}\ddot{x}=-kx</math></center>

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

<center><math>\ddot{x}_G=0\qquad\Rightarrow\qquad x_G=x_{G0}+v_{G0} t</math></center>

y la segunda es la de un oscilador armónico

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

con solución

<center><math>x=c_1\cos(\omega t)+c_2\,\mathrm{sen}(\omega t)</math></center>

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

<center><math>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}</math></center>

lo que nos da, para todo instante

<center><math>x_G=\frac{m_2}{m_1+m_2}(\ell_0+v_0t)</math></center>

Para la elongación del muelle

<center><math>x(t=0)=x_{20}-x_{10}-\ell_0=0\qquad\qquad \dot{x}(t=0)=\dot{x}_{20}-\dot{x}_{10}=v_0</math></center>

siendo por tanto la elongación para todo instante

<center><math>x=\frac{v_0}{\omega}\,\mathrm{sen}(\omega t)</math></center>

Conocidas las dos coordenadas obtenemos las posiciones de las dos masas

<center><math>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)</math></center>
&nbsp;
<center><math>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)</math></center>

==Constantes de movimiento==
Empleando las nuevas coordenadas, vemos que x<sub>G</sub> es cíclica, esto es, no aparece en la lagrangiana. Por tanto su momento conjugado es una constante de movimiento

<center><math>p_G=(m_1+m_2)\dot{x}_G=m_1\dot{x}_1+m_2\dot{x}_2=\mathrm{cte.}</math></center>

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

<center><math>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</math></center>


[[Categoría:Problemas de mecánica analítica (CMR)]]