Ver el código de «Barra oscilando respecto a uno de sus extremos (GIC)»

== Enunciado ==
Una barra homogénea de longitud <math>L</math> y masa <math>M</math> cuelga por uno de sus extremos de modo que se encuentra en equilibio en posición vertical. Analiza el movimiento de la barra si se separa de la vertical un ángulo <math>\theta_0</math> pequeño.

== Solución ==
[[Imagen:F1_GIC_pendulo_fisico_momento_angular.png|right]]
La figura muestra la barra colgando del punto <math>O </math>. En un péndulo ideal se puede considerar que la masa está concentrada en un extremo, pero aquí eso no es posible. Hay que tratar toda la barra como un sólido rígido. Este tipo de sistema se llaman péndulos físicos o compuestos.

La barra realiza un movimiento de rotación respecto al eje <math>OZ </math>. Para describir este movimiento vamos a usar un vector rotación <math>\vec{\omega} = \omega\,vec{k} </math> y el momento angular de la barra respecto al punto <math>O </math>, <math>\vec{L}^O </math>. Como el eje de rotación es fijo, estas dos magnitudes están relacionadas
<center>
<math>
\vec{L}^O = I\vec{\omega}
</math>
</center>
donde <math>I </math> es el momento de inercia de la barra respecto al eje <math>OZ </math>.

La ecuación de movimiento se obtiene de aplicar la "Segunda Ley de Newton para la rotación"
<center>
<math>
\dfrac{\mathrm{d}\vec{L}^O}{\mathrm{d}t} = \vec{M}^O
</math>
</center>
siendo <math>\vec{M}^O </math> el momento neto de las fuerzas que actúan sobre la barra. En este caso esas fuerzas son el peso, aplicado en el centro de masas de la barra, <math>G </math>, y la fuerza de reacción vincular en <math>O </math>, que sujeta a la barra. El momento de esta fuerza de reacción vincular respecto a <math>O </math>. La única fuerza que ejerce momento es el peso
<center>
<math>
\vec{M}^O = \overrightarrow{OG}\times(m\vec{g})
=
\left(d\cos\theta\,\vec{\imath} + d\,\mathrm{sen}\,\theta\,\vec{\jmath}\right)\times (-mg\vec{\imath})
=
-mgd\,\mathrm{sen}\,\theta\,\vec{k}
</math>
</center>
El valor <math>d</math> es la distancia del punto <math>O</math> al centro de masas <math>G</math>. En este caso es la mitad de la longitud de la barra, <math>d=R/2</math>, por lo que
<center>
<math>
\vec{M}^O
=
-\dfrac{1}{2}mgR\,\mathrm{sen}\,\theta\,\vec{k}
</math>
</center>

Por otro lado, la derivada temporal del momento cinético es
<center>
<math>
\dfrac{\mathrm{d}\vec{L}^O}{\mathrm{d}t} 
=
I\dot{\vec{\omega}} = I\ddot{\theta}\,\vec{k}
</math>
</center>
Igualando las dos expresiones obtenemos
<center>
<math>
\ddot{\theta} = -\dfrac{mgR}{2I}\,\mathrm{sen}\,\theta
</math>
</center>
Esta es la ecuación diferencial para <math>\theta(t) </math>. No sabemos resolver esta ecuación, pero si el ángulo de oscilación es pequeño podemos hacer la aproximación <math>\mathrm{sen}\,\theta\approx\theta </math>, y la ecuación queda
<center>
<math>
\ddot{\theta} = -\dfrac{mgR}{2I}\,\theta
</math>
</center>
Esta es la ecuación de un oscilador armónico con frecuencia angular
<center>
<math>
\omega = \sqrt{\dfrac{mgR}{2I}}
</math>
</center>
Es decir, si el ángulo con la vertical es pequeño, la barra oscila con un período
<center>
<math>
T = \dfrac{2\pi}{\omega} = 2\pi\sqrt{\dfrac{2I}{mgR}}
</math>
</center>
El momento de inercia de una barra respecto de un eje perpendicular a ella que pasa por su extremo es
<center>
<math>
I = \dfrac{1}{3}mR^2
</math>
</center>
con lo que el período de oscilación es
<center>
<math>
T = 2\pi\sqrt{\dfrac{2R}{3g}}
</math>
</center>

[[Imagen:F1_GIC_pendulo_fisico_momento_angular_patata.png|right]]
Hay que destacar que, hasta que hemos sustituido el momento de inercia concreto para la barra, el tratamiento es válido para un cuerpo de forma general, por ejemplo el de la figura de la derecha. Lo único que hay que cambiar es dejar el resultado en función de la distancia entre el punto <math>O</math> y el centro de masas del cuerpo, <math>d</math>. Así pues, si el ángulo <math>\theta </math> es pequeño, este cuerpo realizará oscilaciones periódicas con período
<center>
<math>
T = 2\pi\sqrt{\dfrac{I}{mgd}}
</math>
</center>

En particular podemos recuperar el resultado para un péndulo ideal de longitud <math>d </math>. En este caso la masa está concentrada en el extremo de la cuerda, por lo que el momento de inercia es <math>I=md^2 </math>. Sustituyendo en la expresión del período recuperamos el período de un péndulo ideal.


[[Categoría:Problemas de dinámica del sólido rígido]]
[[Categoría:Problemas de movimiento oscilatorio]]
[[Categoría:Movimiento oscilatorio]]