Enunciado

Un sistema está formado por dos varillas homogéneas, ambas de masa y longitud , situadas sobre un plano horizontal (“sólido 1”). La varilla “2” está articulada por su extremo O a un punto fijo del plano, mientras que por su extremo A está articulada a la varilla “3”.

  1. Escriba la lagrangiana del sistema, empleando como coordenadas generalizadas los ángulos que ambas varillas forman con el eje , ϕ (para la varilla 2) y ψ (para la 3).
  2. Obtenga las ecuaciones de movimiento para estos dos ángulos.
  3. ¿Es cíclica alguna de estas coordenadas?

Si en lugar de esas coordenadas se usan el ángulo ϕ que la varilla OA forma con OX y el ángulo θ que AB forma con la prolongación de OA

  1. ¿Cómo queda la lagrangiana?
  2. ¿Y las ecuaciones de movimiento para estos ángulos?
  3. ¿Es cíclica alguna de estas coordenadas?
  4. Determine dos constantes de movimiento para este sistema.
  5. Con ayuda de estas constantes, reduzca el problema a una única ecuación de movimiento para el ángulo θ.

Suponga ahora que la varilla 2 es forzada a girar con velocidad angular constante Ω en torno a O.

  1. Escriba la lagrangiana para este sistema en función del ángulo θ.
  2. Obtenga la ecuación de movimiento para θ. ¿Es la misma que en el apartado (8)?
  3. ¿Se conserva la energía en este sistema? ¿Hay alguna otra constante de movimiento?

Lagrangiana del sistema

La lagrangiana del sistema se calcula como

siendo la energía cinética y U la potencial. En este caso, que no hay ninguna fuerza externa conservativa actuando sobre las varillas, esta energía potencial es nula.

La energía cinética es la suma de las de las dos varillas

La de la varilla 2 se calcula como la de una barra que gira en torno a un eje perpendicular por su extremo

Para la varilla 3 aplicamos el teorema de König

siendo G el centro de masas de la varilla 3.

No hace falta añadir más términos en la energía cinética de rotación porque el eje OZ es uno principal de inercia para ambas varillas y por tanto el momento cinético es paralelo a la velocidad angular

Este término de rotación vale

La velocidad del centro de la varilla 3 la calculamos mediante la expresión del campo de velocidades de un sólido

La velocidad de A la calculamos usando que es una articulación

Sustituimos la posición relativa y la velocidad angular

siendo una base ligada al sólido 2 (aunque por ser un sistema plano , por lo que el subíndice es superfluo). Por otro lado,

Obsérvese que la expresión contiene vectores de dos bases diferentes, por lo que hay que ser cuidadoso a la hora de hallar la energía cinética de traslación

El producto escalar de dos vectores unitarios es igual al coseno del ángulo que forman

Sumando todos los términos en la energía cinética obtenemos la total, que coincide con la lagrangiana

Ecuaciones de movimiento

Estas dos coordenadas son independientes, por lo que podemos emplear las ecuaciones de Lagrange.

Para la coordenada ϕ

La ecuación para esta coordenada es

Calculamos en primer lugar el momento conjugado

siendo su derivada temporal

Por otro lado

Si sustituimos en la ecuación de Lagrange queda

que, simplificada, se reduce a

Para la coordenada ψ

La ecuación para esta coordenada es

El momento conjugado de esta coordenada vale

siendo su derivada temporal

Por otro lado

Sustituimos en la ecuación de Lagrange correspondiente

que, simplificada, queda

Sistema de ecuaciones

Estas dos ecuaciones son acopladas en el sentido de que ambas combinan las dos aceleraciones angulares. Podemos despejar cada una de estas tratando las ecuaciones de Lagrange como un sistema de dos ecuaciones con dos incógnitas ( y ). El resultado es

siendo

Vemos que el sistema de ecuaciones no es en absoluto trivial. Su solución requiere del uso de software para computar la evolución en el tiempo de estas variables.

Cambio de coordenadas

El sistema anterior no permite identificar de forma sencilla la existencia de constantes de movimiento. Sí se puede apreciar que la energía mecánica es una constante, ya que no hay fuerzas que realicen trabajo sobre el sistema, pero no es evidente que haya constantes adicionales.

Una forma de buscar nuevas constantes es eligiendo un sistema distinto de coordenadas generalizadas, de manera que las integrales primeras aparezcan de forma natural.

Tal como se indica en el enunciado, vamos a emplear como coordenadas el mismo ángulo φ que forma el eje con el y el ángulo θ que forma el con el , es decir

Para obtener la nueva lagrangiana no hace falta calcular la energía cinética desde el principio. Basta sustituir en la expresión final. El resultado es

Nuevas ecuaciones de movimiento

Operamos de nuevo para obtener las ecuaciones de movimiento.

Para la coordenada ϕ

Calculamos el momento conjugado

Obsérvese que no tiene el mismo valor que con la elección inicial de coordenadas. Si agrupamos términos y abreviamos el coseno como C

La derivada temporal del momento conjugado vale ahora

Por otro lado

con lo que la ecuación de movimiento, simplificada, es

Para la coordenada θ

La ecuación para esta coordenada es

El momento conjugado de esta coordenada vale

Agrupamos y abreviamos

Derivamos respecto al tiempo

Por otro lado

Sustituimos en la ecuación de Lagrange correspondiente

que, simplificada, queda

Constantes de movimiento

Asociada a ϕ

En la segunda forma de la lagrangiana que hemos escrito, la coordenada ϕ es ciclica, esto es, no aparece explícitamente en la lagrangiana

lo cual implica, por la ecuación de Lagrange, que

y por tanto el momento conjugado

es una constante de movimiento.

Esta cantidad está asociada a la invariancia del sistema ante una rotación, es decir, que si giramos todo el sistema un cierto ángulo Δϕ el sistema no se ve modificado.

En términos de las magnitudes vectoriales, esta constante representa la componente z del momento cinético del sistema respecto al origen O.

Independencia del tiempo

Por otro lado, la lagrangiana tampoco depende explícitamente del tiempo

lo cual implica que la función hamiltoniana

es otra constante de movimiento.

Dado que la energía cinética es una función cuadrática de las velocidades, esta cantidad coincide con la energía mecánica del sistema

Como además la energía potencial es nula, se deduce que la energía cinética

es otra constante en este problema.

Reducción del sistema

La constancia asociada a ϕ permite expresar las derivadas de esta variable en función de cantidades constantes. La velocidad generalizada vale

y la aceleración generalizada, derivando esta ecuación,

Podemos llevar estas dos expresiones a la ecuación de movimiento para θ

Resulta, tras un laborioso cálculo realizado con ayuda de Mathematica

Este resultado ilustra que no siempre la obtención de una constante de movimiento produce un resultado más sencillo.

Podemos aprovechar también la conservación de la energía sustituyendo directamente el valor de en esta constante. En ese caso sí llegamos a una ecuación más sencilla

Lagrangiana con varilla rotatoria

En el caso de que la varilla 2 tenga una rotación fijada nos vincula la coordenada ϕ

Obtenemos la nueva lagrangiana por simple sustitución

Obsérvese que ahora tenemos un solo grado de libertad, el del ángulo θ

Ecuación con varilla rotatoria

Obtenemos la ecuación de movimiento de la forma habitual.

Primero hallamos el momento conjugado

Derivamos respecto al tiempo

Por otro lado tenemos

Lo cual nos da la ecuación

que simplificada nos queda

Vemos que, aunque en la lagrangiana nos hemos limitado a sustituir, resulta una ecuación completamente diferente al caso anterior.

En este caso hemos obtenido que verifica la misma ecuación que un péndulo, es decir, que la barra 3 realiza un movimiento oscilatorio (o de rotación) alrededor de la articulación A.

Magnitud conservada

La lagrangiana sigue siendo independiente del tiempo, por lo que la función hamiltoniana es una constante de movimiento. No obstante, su expresión no es la misma que cuando no estaba fijada.

y resulta

Esta cantidad no coincide con la energía mecánica del sistema (que en este caso sería igual a la propia lagrangiana). La razón es que ahora la energía cinética no es una función cuadrática de la velocidad, sino que es de la forma

Cuando la lagrangiana es de esta forma, la cantidad coinservada es