Ver el código de «Anilla ensartada en un aro rodante»

==Enunciado==
Se tiene un sistema formado por un aro &ldquo;2&rdquo; de masa <math>m_2</math> que puede rodar sin deslizar sobre una superficie horizontal &ldquo;1&rdquo;. Ensartado en este aro se encuentra una pequeña anilla &ldquo;3&rdquo; de masa <math>m_3</math> que puede deslizarse sin fricción a lo largo del aro &ldquo;2&rdquo;.

<center>[[Archivo:particula-aro.png]]</center>

Empleando como coordenadas la posición x del centro C del aro &ldquo;2&rdquo; a lo largo del eje horizontal y el ángulo &theta; que la posición de la anilla &ldquo;3&rdquo; forma con la vertical, determine:
# Las ecuaciones de movimiento para estas dos coordenadas.
# Dos constantes de movimiento no triviales.
# Suponiendo que el aro y la anilla están en reposo, se sitúa la anilla formando un pequeño ángulo <math>\theta_0</math> con la vertical. Determine la frecuencia de las oscilaciones que realiza el sistema.

==Ecuaciones de movimiento==
Cada uno de los dos cuerpos &ldquo;2&rdquo; y &ldquo;3&rdquo; obedece las ecuaciones dadas por el teorema de la cantidad de movimiento y el teorema del momento cinético. En cada caso debemos incluir las posibles fuerzas de reacción sobre cada uno de los cuerpos.
===Para el aro===
====Teorema de la cantidad de movimiento====
El aro experimenta las siguientes fuerzas externas:
* Su peso <math>-m_2g\vec{\jmath}_1</math>
* La reacción del suelo en el punto de contacto A, la cual tiene dos componentes
** Una vertical <math>\vec{F}_{Ay}=F_{Ay}\vec{\jmath}_1</math> que impide que el aro atraviese el suelo, y mantiene el centro a una altura constante R.
** Una horizontal, <math>\vec{F}_{Ax}=F_{Ax}\vec{\imath}_1</math> que impide que el aro deslice, de manera que solo ruede.
* Una fuerza de reacción debida a la anilla. El vínculo de que la anilla permanezca en el aro se consigue realizando una fuerza normal sobre ésta <math>\vec{F}_P</math> y, por la tercera ley de Newton, una fuerza igual y de sentido contrario <math>-\vec{F}_P</math> se aplica en el aro. 

Esto nos da la sigueinte relación

<center><math>m_2\vec{a}^C_{21}=m_2\vec{g}+\vec{F}_{Ax}+\vec{F}_{Ay}-\vec{F}_P</math></center>

La aceleración del centro del aro es puramente horizontal

<center><math>\vec{a}^C_{21}=\ddot{x}\vec{\imath}_1</math></center>

Para expresar la fuerza <math>\vec{F}_P</math> usamos un sistema &ldquo;3&rdquo; de ejes que está girado un ángulo &theta; respecto al sistema &ldquo;1&rdquo;, de manera que <math>\vec{\imath}_3</math> es tangente al aro y <math>\vec{\jmath}_3</math> es radial y hacia adentro. 

<center><math>\vec{\imath}_3 = C\vec{\imath}_1+S\vec{\jmath}_1\qquad\qquad \vec{\jmath}_3=-S\vec{\imath}_1+C\vec{\jmath}_1</math></center>

(<math>S=\mathrm{sen}(\theta),\ C=\cos(\theta)</math>). Estos vectores se relacionan directamente con los correspondientes a las coordenadas polares.

<center>[[Archivo:aro-anilla-ejes.png]]</center>
Tenemos entonces

<center><math>m_2 \ddot{x}\vec{\imath}_1=-m_2g\vec{\jmath}_1+F_{Ax}\vec{\imath}_1+F_{Ay}\vec{\jmath}_1-F_P\vec{\jmath}_3</math></center>

Separando por componentes

<center><math>\left\{\begin{array}{lrcl}x:\ & m_2\ddot{x}&=&F_{Ax}+F_PS \\ y:\ &0&=&-mg+F_{Ay}-F_PC\end{array}\right.</math></center>

Estas dos ecuaciones no son suficientes para determinar el movimiento del aro.

====Teorema del momento cinético====
La tercera ecuación para el aro la da la evolución del momento cinético. Respecto al centro del aro

<center><math>\vec{L}_C=I\vec{\omega}_{21} \qquad\Rightarrow\qquad L_C=I\omega_{21}</math></center>

Podemos usar escalares ya que ambos vectores apuntan en la dirección del eje Z, perpendicular al plano de movimiento. Podemos relacionar la velocidad angular con la velocidad de avance del centro por la condición de rodadura

<center><math>\vec{0}=\vec{v}^A_{21}+\omega_{21}\vec{k}\times\overrightarrow{GA}\qquad\Rightarrow\qquad \omega_{21}=-\frac{\dot{x}}{R}</math></center>

y de aquí

<center><math>\alpha_{21}=-\ddot{x}{R}</math></center>

La única fuerza que provoca un momento respecto al centro es la tangencial en el punto de contacto. Por ello

<center><math>I\alpha=F_{Ax}R\qquad\Rightarrow\qquad m_2R^2\left(-\frac{\ddot{x}}{R}\right)=F_{Ax}R\qquad\Rightarrow\qquad F_{Ax}=-m_2\ddot{x}</math></center>

Llevando esto a la componente horizontal del TCM podemos eliminar esta fuerza de reacción y queda

<center><math>2m_2\ddot{x}=F_PS</math></center>

Para completar el sistema necesitamos además las ecuaciones para la anilla.

===Para la anilla===
Considerando la anilla como una partícula, no es necesario considerar la rotación. En ese caso nos basta con el TCM, es decir, con la 2ª ley de Newton.

<center><math>m_3\vec{a}^P_{31}=m_3\vec{g}+\vec{F}_P</math></center>

La aceleración de la partícula se relaciona con la del centro del aro como

<center><math>\vec{a}^P_{31}=\vec{a}^P_{31}+\alpha_{31}\vec{k}\times\overrightarrow{CP}-\omega^2_{31}\overrightarrow{CP}</math></center>

siendo

<center><math>\vec{a}^C_{31}=\vec{a}^C_{21}=\ddot{x}\vec{\imath}_1</math></center>

(aquí estamos imaginando un sólido o sistema de referencia "3" cuyo eje <math>OY_3</math> sería el radio por C y P).

Además

<center><math>\alpha_{31}=\ddot{\theta}\qquad\qquad \omega_{31}=\dot{\theta}\qquad\qquad \overrightarrow{CP}=-R\vec{\jmath}_3</math></center>

Esto nos da

<center><math>m_3\left(\ddot{x}\vec{\imath}_1+R\ddot{\theta}\vec{\imath}_3+R\dot{\theta}^2\vec{\jmath}_3\right)=-m_3g\vec{\jmath}_1+F_P\vec{\jmath}_3</math></center>

Esta expresión combina dos bases vectoriales. para obtener una ecuación escalar debemos proyectar sobre una dirección dada. 

===Ecuaciones de movimiento===
Si proyectamos la ecuación vectorial de la anilla hacemos sobre <math>\vec{\imath}_3</math> eliminamos la fuerza de reacción <math>F_P</math> y queda

<center><math>m_3\left(\ddot{x}C+\ddot{\theta}R\right)=-m_3gS\qquad\Rightarrow\qquad \ddot{x}C+\ddot{\theta}R=-gS</math></center>

Si proyectamos sobre la dirección horizontal <math>\vec{\imath}_1</math>

<center><math>m_3\left(\ddot{x}+R\ddot{\theta}C-R\dot{\theta}^2S\right)=-F_PS</math></center>

Sumamos esta ecuación con la correspondiente al aro

<center><math>2m_2\ddot{x}=F_PS</math></center>

y nos queda

<center><math>(2m_2+m_3)\ddot{x}+m_3R\ddot{\theta}C-m_3R\dot{\theta}^2S=0</math></center>

Tenemos entonces el sistema
<center><math>
\left\{\begin{array}{rcccl}
\ddot{x}C&+&\ddot{\theta}R&=&-gS \\
(2m_2+m_3)\ddot{x}&+&m_3R\ddot{\theta}C&=&m_3R\dot{\theta}^2S
\end{array}\right.</math></center>

De aquí podemos despejar las segundas derivadas

<center><math>\left\{\begin{array}{rcl}
\ddot{x}& = & \dfrac{m_3S(gC+\dot{\theta}^2R)}{2 m_2 + m_3S^2} \\
&&\\
\ddot{\theta} & = & -\dfrac{(2 m_2 + m_3)gS  + m_3RCS\dot{\theta}^2 }{(2 m_2 + m_3S^2) R}
\end{array}\right.</math></center>

La segunda de las ecuaciones contiene solo a la variable &theta; y sus derivadas, con lo que puede ser resuelta sin tener en cuenta la posición x del anillo, que puede determinarse posteriormente.

Como límite de interés podemos ver qué pasa si <math>m_2\to\infty</math>, es decir, el aro es infinitamente masivo. En ese caso

<center><math>\ddot{x}=0\qquad\qquad \ddot{\theta}=-\frac{g}{R}\mathrm{sen}(\theta)\qquad\qquad (m_2\to\infty)</math></center>

es decir, el aro no se mueve, mientras que la anilla cumple la ecuación del péndulo.

==Constantes de movimiento==
===Cantidad de movimiento generalizada===
En este sistema no se conserva la cantidad de movimiento, 

* En su componente vertical, porque está el peso, que acelera a la anilla en esta dirección.
* En la componente horizontal, porque hay una fuerza externa <math>F_{Ax}</math> que acelera el CM del sistema.

No obstante, si examinamos la ecuación de movimiento

<center><math>(2m_2+m_3)\ddot{x}+m_3R\ddot{\theta}C-m_3R\dot{\theta}^2S=0</math></center>

vemos que esta ecuación equivale a

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t}((2m_2+m_3)\dot{x}+m_3 R\dot{\theta}C)=0</math></center>

y por tanto, una constante no trivial es

<center><math>(2m_2+m_3)\dot{x}+m_3 R\dot{\theta}C = B \mathrm{cte.}</math></center>

Esta magnitud generaliza la cantidad de movimiento, incorporando un término de rotación.

De hecho, puede integrarse una segunda vez y dar

<center><math>(2m_2+m_3)x+m_3 RS = A+Bt</math></center>

con lo cual, si quisiéramos, podríamos obtener una segunda constante no trivial como

<center><math>(2m_2+m_3)x+m_3 RS - Bt = (2m_2+m_3)x+m_3 RS-((2m_2+m_3)\dot{x}+m_3 R\dot{\theta}C) t</math></center>

No obstante, hay una segunda constante con mayor significado físico.
===Energía mecánica===
En este sistema se conserva la energía mecánica. La justificación para ello es que:
* El peso es una fuerza conservativa que deriva de una energía potencial.
* La fuerza de reacción <math>\vec{F}_A</math> se aplica sobre un punto de velocidad nula (el de contacto con el suelo) y por tanto no desarrolla potencia.
* La fuerza interna <math>\vec{F}_P</math> podría, en principio, afectar a la energía, sin embargo, su potencia es suma de la que hace sobre cada sólido

<center><math>P=\vec{F}_P\cdot\vec{v}^P_{31}+(-\vec{F}_P)\cdot\vec{v}^P_{21}=\vec{F}_P\cdot\left(\vec{v}^P_{31}-\vec{v}^P_{21}\right)=\vec{F}_P\cdot\vec{v}^P_{32}=0</math></center>

:ya que esta fuerza es radial y la velocidad de la anilla respecto al aro es tangencial.

Por tanto la energía mecánica se conserva. Ésta es suma de

* La energía cinética del aro, suma de la de traslación y la de rotación.

<center><math>T_2=\frac{1}{2}m_2\left|\vec{v}^C_{21}\right|^2+\frac{1}{2}I\left|\vec{\omega}_21\right|^2 = \frac{1}{2}m_2\dot{x}^2+\frac{1}{2}\left(m_2R^2\right)\left(-\frac{\dot{x}}{R}\right)^2=m_2\dot{x}^2</math></center>

* La energía cinética de la partícula

<center><math>T_2=\frac{1}{2}m_3\left|\vec{v}^P_{31}\right|^2=\frac{1}{2}m_3\left|\dot{x}\vec{\imath}_1+R\dot{\theta}\vec{\imath}_3\right|^2=\frac{1}{2}m_3\left(\dot{x}^2+2R\dot{x}\dot{\theta}C+R^2\dot{\theta}^2\right)</math></center>

* La energía potencial de la partícula, que podemos medir desde la altura del centro del aro

<center><math>U_3=-m_3gC\,</math></center>

Sumando estas tres cantidades resulta la energía

<center><math>E=m_2\dot{x}^2+\frac{1}{2}m_3\left(\dot{x}^2+2R\dot{x}\dot{\theta}C+R^2\dot{\theta}^2\right)-m_3gC</math></center>

Podemos comprobar que esta cantidad es una constante derivando respecto al tiempo y sustituyendo las ecuaciones de movimiento.

==Pequeñas oscilaciones==
La ecuación de movimiento para el ángulo &theta; es

<center><math>\ddot{\theta}  =  -\dfrac{(2 m_2 + m_3)gS  + m_3RCS\dot{\theta}^2 }{(2 m_2 + m_3S^2) R}</math></center>

Si consideramos pequeñas oscilaciones alrededor del punto más bajo podemos hacer las aproximaciones

<center><math>S=\mathrm{sen}(\theta)\simeq \theta\qquad\qquad C=\cos(\theta)\simeq 1\qquad\qquad \dot{\theta}^2\ll 1</math></center>

con esto queda la ecuación aproximada

<center><math>\ddot{\theta}  =  -\dfrac{(2 m_2 + m_3)g  }{2 m_2 R}\theta</math></center>

que es la de un oscilador armónico con frecuencia

<center><math>\omega=\sqrt{\dfrac{(2 m_2 + m_3)g  }{2 m_2 R}}</math></center>

Para <math>m_2\to\infty</math> se recupera de nuevo la ecuación del péndulo simple.