Ver el código de «Estudio analítico de una barra apoyada»

==Enunciado==
Supongamos que tenemos una barra de masa m y longitud b apoyada en el suelo y en una pared vertical, sometida a la acción del peso (vertical y hacia abajo) y a las fuerzas de reacción en los puntos de contacto. No hay rozamiento con las superficies

<center>[[Archivo:Barra-apoyada.png|250px]]</center>
# Determine la lagrangiana del sistema.
# Halle la ecuación de movimiento para el ángulo θ.
# Determine una constante de movimiento no trivial.
# Añadiendo una coordenada x que representaría la separación de la barra respecto de la pared vertical, calcule la fuerza de reacción ejercida por la pared.
# Existe un valor de θ para el cual la barra se separa de la pared. Determine este valor.
# Halle la ecuación de movimiento para la barra una vez que se ha separado de la pared.

==Lagrangiana del sistema==
===Energía cinética===
Podemos hallar la energía cinética a partir de la expresión general para un sólido

<center><math>T=\frac{1}{2}m|\vec{v}_G|^2+\frac{1}{2}\vec{\omega}\cdot\bar{\bar{I}}\cdot\vec{\omega}</math></center>

En el caso de un movimiento plano esta expresión se reduce a

<center><math>T=\frac{1}{2}n(\dot{x}_G^2+\dot{y}_G^2)+\frac{1}{2}I_{zz}\omega^2</math></center>

La posición del CM es, en función de &theta;

<center><math>x_G=\frac{b}{2}\mathrm{sen}(\theta)=\frac{b}{2}S\qquad\qquad y_G=\frac{b}{2}\mathrm{cos}(\theta)=\frac{b}{2}C</math></center>

y su velocidad

<center><math>\dot{x}_G=\frac{b}{2}\dot{\theta}C\qquad\qquad \dot{y}_G=-\frac{b}{2}\dot{\theta}S</math></center>

Esto nos da la energía cinética

<center><math>T=\frac{1}{2}\left(\frac{b^2}{4}\dot{\theta}^2\right)+\frac{1}{2}\left(\frac{1}{12}mb^2\right)\dot{\theta}^2=\frac{mb^2}{6}\dot{\theta}^2</math></center>

También podemos llegar a este resultado empleando la expresión para una barra conocidas las velocidades de sus extremos A y B

<center><math>T=\frac{m}{6}(|\vec{v}_A|^2 +|\vec{v}_B|^2+\vec{v}_A\cdot\vec{v}_B)</math></center>

===Energía potencial===
La energía potencial gravitatoria la da la altura del CM

<center><math>U=mgy_G=\frac{mgb}{2}C</math></center>

===Lagrangiana===
Restamos las dos cantidades

<center><math>\mathcal{L}=T-U=\frac{mb^2}{6}\dot{\theta}^2-\frac{mgb}{2}C</math></center>

==Ecuación de movimiento==
Aplicamos la ecuación de Lagrange a la coordenada &theta;

<center><math>\frac{\partial\mathcal{L}}{\partial\dot{\theta}}=\frac{mb^2}{3}\dot{\theta}\qquad\qquad \frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial\dot{\theta}}\right)=\frac{mb^2}{3}\ddot{\theta}\qquad\qquad \frac{\partial\mathcal{L}}{\partial{\theta}}=\frac{mgb}{2}S</math></center>

lo que nos da la ecuación de movimiento

<center><math>\frac{mb^2}{3}\ddot{\theta}-\frac{mgb}{2}S=0\qquad\Rightarrow\qquad \ddot{\theta}=\frac{3g}{2b}\mathrm{sen}(\theta)</math></center>

Vemos que empleando la mecánica analítica llegamos al resultado de forma mucho más corta que mediante mecánica vectorial, ya que no precisamos considerar la acción de las fuerzas de reacción.

==Constante de movimiento==
La lagrangiana de este problema no depende del tiempo y la energía cinética es una función cuadrática de la velocidad generalizada. Por tanto se conserva la energía mecánica

<center><math>E=T+U=\frac{mb^2}{6}\dot{\theta}^2+\frac{mgb}{2}C=\mathrm{cte.}</math></center>
==Fuerza de reacción==
La fuerza de reacción está asociada al hecho de que el extremo superior B de la barra está adosado a la pared. Si x es la distancia de B a la pared,

<center><math>x=0\,</math></center>

Para obtener la fuerza de reacción que garantiza este vínculo debemos suponer que la barra podría atravesar la pared, pero si no lo hace es porque hay una fuerza que se lo impide. Esta fuerza se calcula mediante la ecuación de Lagrange para la coordenada x

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial\dot{x}}\right)-\frac{\partial \mathcal{L}}{\partial x}=\lambda = F_x</math></center>

Sin embargo, para llegar a este resultado debemos calcular de nuevo la lagrangiana, ya que antes la hallamos presuponiendo el vínculo <math>x=0</math>. Si deshacemos este vínculo y suponemos que la barra puede estar a una distancia x de la pared, la posición del CM es

<center><math>x_G=x+\frac{b}{2}S\qquad\qquad y_G=\frac{b}{2}C</math></center>

y su velocidad

<center><math>\dot{x}_G=\dot{x}+\frac{b}{2}\dot{\theta}C\qquad\qquad \dot{y}_G=-\frac{b}{2}\dot{\theta}S</math></center>

lo que llevado a la lagrangiana nos da

<center><math>\mathcal{L}=T-U=\frac{1}{2}m\dot{x}^2+\frac{mb}{2}\dot{x}\dot{\theta}C+\frac{mb^2}{6}\dot{\theta}^2-\frac{mgb}{2}C</math></center>

Calculamos ahora el momento conjugado de la coordenada x y su derivada temporal

<center><math>p_x=\frac{\partial\mathcal{L}}{\partial\dot{x}}=m\dot{x}+\frac{mb}{2}\dot{\theta}C\qquad\qquad \frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial\dot{x}}\right)=m\ddot{x}+\frac{mb}{2}\ddot{\theta}C-\frac{mb}{2}\dot{\theta}^2S</math></center>

x es una coordenada cíclica, ya que no aparece en la lagrangiana, por lo que la fuerza en x es igual a

<center><math>F_x=m\ddot{x}+\frac{mb}{2}\ddot{\theta}C-\frac{mb}{2}\dot{\theta}^2S</math></center>

Ahora bien, esta fuerza debe ser tal que se cumpla el vínculo de que x se anule, y lo mismo con sus derivadas respecto al tiempo, lo que nos da finalmente

<center><math>F_x=\frac{mb}{2}\ddot{\theta}C-\frac{mb}{2}\dot{\theta}^2S</math></center>
==Separación de la pared==
Al ser el vínculo unilateral, la barra se separa de la pared en el momento en que la fuerza <math>F_x</math> se anule, ya que no puede ser negativa. Por tanto, la condición de separación es

<center><math>\frac{mb}{2}\ddot{\theta}C-\frac{mb}{2}\dot{\theta}^2S =0</math></center>

Para ver en qué ángulo se produce la separación debemos sustituir la velocidad y la aceleración angular en función del ángulo. La aceleración la obtenemos de la ecuación de movimiento

<center><math>\ddot{\theta}=\frac{3g}{2b}S</math></center>

y la velocidad angular al cuadrado de la conservación de la energía mecánica

<center><math>E=\frac{mb^2}{6}\dot{\theta}^2+\frac{mgb}{2}C\qquad\Rightarrow\qquad \dot{\theta}^2=\frac{6E}{mb^2}-\frac{3g}{b}C</math></center>

Llegamos entonces a la ecuación

<center><math>\frac{3g}{2b}SC-\left(\frac{6E}{mb^2}-\frac{3g}{b}C\right)S=0</math></center>

Simplificamos aquí y llegamos al ánguloa de separación

<center><math>\theta_s=\arccos\left(\frac{4E}{3mgb}\right)</math></center>

En el caso particular de que la barra se halle inicialmente vertical, apoyada en la pared

<center><math>(\theta_0=0)\qquad\qquad E=\frac{mgb}{2}\qquad\Rightarrow\qquad\theta_s=\arccos\left(\frac{2}{3}\right)</math></center>

==Movimiento tras la separación==
Una vez que separa, desaparace la fuerza <math>F_x</math>. A partir de ese momento la ecuación de movimiento que se cumple es

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial\dot{x}}\right)-\frac{\partial \mathcal{L}}{\partial x}=m\ddot{x}+\frac{mb}{2}\ddot{\theta}C-\frac{mb}{2}\dot{\theta}^2S</math></center>

La ecuación para el ángulo &theta; debe calcularse de nuevo empleando la lagrangiana con dos grados de libertad

<center><math>\mathcal{L}=\frac{1}{2}m\dot{x}^2+\frac{mb}{2}\dot{x}\dot{\theta}C+\frac{mb^2}{6}\dot{\theta}^2-\frac{mgb}{2}C</math></center>

y queda

<center><math>\frac{\partial\mathcal{L}}{\partial\dot{\theta}}=\frac{mb}{2}\dot{x}C+\frac{mb^2}{3}\dot{\theta}\qquad\qquad \frac{\mathrm{d}\ }{\mathrm{d}t}\left(\frac{\partial\mathcal{L}}{\partial\dot{\theta}}\right)=\frac{mb}{2}\ddot{x}C-\frac{mb}{2}\dot{x}\dot{\theta}S+\frac{mb^2}{3}\ddot{\theta}\qquad\qquad \frac{\partial\mathcal{L}}{\partial{\theta}}=-\frac{mb}{2}\dot{x}\dot{\theta}S+\frac{mgb}{2}S</math></center>

lo que nos da la ecuación

<center><math>\frac{3}{2b}\ddot{x}C+\ddot{\theta}-\frac{3g}{2b}S=0</math></center>

Combinando esta ecuación con la de movimiento en x, obtenemos una ecuación exclusivamente para &theta;

<center><math>\ddot{\theta}=\frac{3 S (2 g - b C \dot{\theta}^2)}{b (4 - 3 C^2)}</math></center>
[[Categoría:Problemas de mecánica analítica (CMR)]]