Ver código fuente

==Enunciado==
Se tiene un rotor formado por dos masas iguales de valor <math>m</math> situadas en los extremos de una barra ideal (sin masa) de longitud <math>H</math>. Cuando este rotor está equilibrado gira en torno a un eje perpendicular a la barra y que pasa por su centro. Este eje está anclado en dos rodamientos situados a una distancia <math>h</math> del centro de la barra (uno por encima y otro por debajo de ella).

Calcule el momento cinético y la energía cinética (respecto a un sistema fijo y respecto al CM) si el rotor gira con velocidad angular constante <math>\omega</math> en torno al eje cuando:

# Es horizontal y se encuentra centrado en el eje vertical.
# Es horizontal pero se encuentra descentrado de forma que el eje no pasa por el centro de la barra, sino a una distancia <math></math>b de éste.
# Está centrado pero la barra está inclinada respecto a la horizontal un ángulo &beta;.
# Es horizontal y se encuentra centrado en el eje vertical, pero las masas no son exactamente iguales, sino que valen <math>m_1</math> y <math>m_2</math>.

<center>[[Archivo:rotor-desequilibrado-01.png]]{{qquad}}{{qquad}}[[Archivo:rotor-desequilibrado-02.png]]</center>

==Introducción==
Esta es la primera parte de un problema de dinámica que continúa en [[Rotores desequilibrados|otro problema]]. Se trata de analizar el comportamiento de un rotor formado por dos masas. 

Aunque los casos son diferentes, especialmente en el aspecto dinámico de las fuerzas que se hallan en la segunda parte de este problema, el cálculo es similar en los cuatro casos, por lo que podemos hacer la mayor parte de los cálculos de forma general.

En todos los casos, tenemos dos masas <math>m_1</math> y <math>m_2</math> describiendo órbitas circulares en torno al eje. 

<center><math>\vec{v}_1=\vec{\omega}\times\vec{r}_1 \qquad\qquad \vec{v}_2=\vec{\omega}\times\vec{r}_2</math></center>

El momento cinético del sistema será la suma de los de las dos partículas

<center><math>\vec{L}_O = m_1\vec{r}_1\times\vec{v}_1 + m_2\vec{r}_2\times\vec{v}_2</math></center>

y de manera análoga para la energía cinética

<center><math>K = \frac{1}{2}m_1\left|\vec{v}_1\right|^2+\frac{1}{2}m_2\left|\vec{v}_2\right|^2</math></center>

Para hallar los valores respecto al CM podemos aplicar las descomposiciones

<center><math>\vec{L}_O = M\vec{r}_C\times\vec{v}_C + \vec{L}_C\qquad\qquad K = \frac{1}{2}M\left|\vec{v}\right|^2 + K'</math></center>

===Caso de una sola partícula===
En lo que sigue, en todos los casos tanto el momento como la energía cinética son suma de las contribuciones de dos partículas que describen un movimiento circular alrededor de un eje. Por ello, en lugar de analizar cada problema independientemente, repitiendo los mismos cálculos varias veces, es preferible hallar previamente los resultados para una sola partícula y posteriormente sumar las veces que haga falta.

Cuando tenemos una partícula en un movimiento circular de radio <math>R</math>, su velocidad se puede escribir como

<center><math>\vec{v}=\vec{\omega}\times\vec{r}</math></center>

siendo <math>\vec{r}</math> el vector de posición respecto a un punto O del eje (no necesariamente el centro de la circunferencia). El momento cinético de esta partícula respecto a este punto valdrá

<center><math>\vec{L}_O=m\vec{r}\times\vec{v}=m\vec{r}\times\left(\vec{\omega}\times\vec{r}\right)</math></center>

Desarrollamos el doble producto vectorial

<center><math>\vec{L}_O = m\left|\vec{r}\right|^2\vec{\omega}-2m(\vec{r}\cdot\vec{\omega})\vec{r}</math></center>

En el caso particular de que el vector de posición sea ortogonal a la velocidad angular (es decir, si O es el centro de la circunferencia), esto se reduce a

<center><math>\vec{\omega}\cdot\vec{r}=0\qquad\Rightarrow\qquad \vec{L}_O = mR^2\vec{\omega}</math></center>

En este caso el momento cinético es proporcional a la velocidad angular, siendo la constante de proporcionalidad el llamado momento de inercia. 

En general, O puede ser un punto cualquiera del eje. Podemos escribir el vector de posición como

<center><math>\vec{r}=\rho\vec{u}_\rho+z\vec{k}=r\cos(\beta)\vec{u}_\rho+r\,\mathrm{sen}(\beta)\vec{k}</math></center>

donde <math>r</math> no es el radio de la circunferencia, sino la distancia al punto O. El radio se relaciona con ésta como

<center><math>R = r\cos(\beta)\,</math></center>

La velocidad para esta partícula, en este caso general, apunta en la dirección acimutal,

<center><math>\vec{v}=\vec{\omega}\times\vec{r}=\omega r\cos(\beta)\vec{u}_\varphi=\omega R\vec{u}_\varphi</math></center>

por lo que la rapidez es

<center><math>\left|\vec{v}\right|=\omega R</math></center>

El momento cinético para este caso general se puede escribir

<center><math>\vec{L}_O=m\vec{r}\times\vec{v}=m\left(r\cos(\beta)\vec{u}_\rho+r\,\mathrm{sen}(\beta)\vec{k}\right)\times\left(\omega r\cos(\beta)\vec{u}_\varphi\right) = mr^2\cos(\beta)\left(-\mathrm{sen}(\beta)\vec{u}_\rho+\cos(\beta)\vec{k}\right)</math></center>

Vemos que en el caso general el momento cinético no es paralelo a la velocidad angular (que va en la dirección de <math>\vec{k}</math>).

La energía cinética de esta partícula es

<center><math>K = \frac{1}{2}m|\vec{v}|^2 = \frac{1}{2}m\omega^2R^2 = \frac{1}{2}m\omega r^2\cos^2(\beta)</math></center>

Se cumple la relación general

<center><math>K = \frac{1}{2}\omega\cdot\vec{L}_O</math></center>

==Rotor equilibrado==
En el primer caso, las dos masas son iguales

<center><math>m_1=m_2=m\,\qquad\qquad M = m_1+m_2=2m</math></center>

las posiciones son simétricas respecto al eje
<!--\vec{r}_1=\frac{H}{2}\cos(\omega t)\vec{\imath}+\frac{H}{2}\mathrm{sen}(\omega t)\vec{\jmath}\qquad\qquad -->
<center><math>\vec{r}_2=-\vec{r}_1</math></center>

y derivando, resultan también velocidades opuestas

<center><math>\vec{v}_2=-\vec{v}_1</math></center>

Puesto que la barra es horizontal, el ángulo de inclinación es nulo, <math>\beta = 0</math> y el centro de reducción coincide con el centro de la circunfrencia, cuyo radio es

<center><math>R = \frac{H}{2}</math></center>

Esto implica para el momento cinético que los dos sumandos son iguales

<center><math>\vec{L}_O = 2m\vec{r}_1\times\vec{v}_1</math></center>

y aplicando la expresión para una sola partícula en movimiento circular queda

<center><math>\vec{L}_O = \frac{MH^2}{4}\vec{\omega}</math></center>

Puesto que el centro de masas coincide con el centro de la varilla también se cumple

<center><math>\vec{L}_C=\vec{L}_O = \frac{MH^2}{4}\vec{\omega}</math></center>

Para la energía cinética el cálculo es similar. También son iguales las energías de las dos partículas.

<center><math>K = 2\left(\frac{1}{2}m|\vec{v}_1|^2\right)</math></center>

Por tratarse de un movimiento circular uniforme la rapidez es igual a la velocidad angular multiplicada por el radio de giro

<center><math>K = \frac{1}{2}\left(\frac{MH^2}{4}\right)\omega^2</math></center>

y, por coincidir el origen con el CM

<center><math>K' = K = \frac{1}{2}\left(\frac{MH^2}{4}\right)\omega^2</math></center>

==Rotor descentrado==
En el segundo caso, las dos masas son de nuevo iguales

<center><math>m_1=m_2=m\,\qquad\qquad M = m_1+m_2=2m</math></center>

Las dos partículas describen de nuevo circunferencias horizontales, pero de distinto radio:

<center><math>R_1 = \frac{H}{2}+b\qquad\qquad R_2 = \frac{H}{2}-b</math></center>

Esto provoca que los dos sumandos en el momento cinético ya no sean iguales, pero cada uno se calcula como antes, empleando el radio correspondiente

<center><math>\vec{L}_O = m_1R_1^2\vec{\omega}+m_2R_2^2\vec{\omega}=\frac{M}{2}\left(R_1^2+R_2^2\right)\vec{\omega}</math></center>

siendo

<center><math>R_1^2+R_2^2 =\left(\frac{H}{2}+b\right)^2+\left(\frac{H}{2}-b\right)^2=\frac{H^2}{2}+2b^2</math></center>

Nos queda finalmente

<center><math>\vec{L}_O=M\left(\frac{H^2}{4}+b^2\right)\vec{\omega}</math></center>

El centro de masas sigue estando en el centro de la varilla (que ya no coincide con el punto del eje). Este punto a una distancia <math>b</math> del eje, describiendo una circunferencia con velocidad angular <math>\vec{\omega}</math>. Por tanto el momento cinético por moverse con el CM valen

<center><math>M\vec{r}_C\times\vec{v}_C = M|\vec{r}_C^2|\vec{\omega}=Mb^2\vec{\omega}</math></center>

Restamos este término de la expresión completa y nos queda el momento cinético respecto al CM

<center><math>\vec{L}_C = \vec{L}_O-M\vec{r}_C\times\vec{v}_C = \frac{MH^2}{4}\vec{\omega}</math></center>

Podemos ver que el momento respecto al CM es el mismo que en el caso simétrico.

==Rotor inclinado==
En el caso del rotor inclinado, las posiciones vuelven a ser simétricas

<center><math>\vec{r}_2=-\vec{r}_1</math></center>

y también lo son las velocidades

<center><math>\vec{v}_2=-\vec{v}_1</math></center>

Por ello, de nuevo el vomento cinético es el doble del de una de las partículas

<center><math>\vec{L}_O = 2m\vec{r}_1\times\vec{v}_1</math></center>

pero ahora debemos tener en cuenta la inclinación de la varilla y usar la expresión general para un cierto ángulo <math>\beta</math>

<center><math>\vec{L}_O=\frac{MH^2}{4}\cos(\beta)\left(-\mathrm{sen}(\beta)\vec{u}_\rho+\cos(\beta)\vec{k}\right)</math></center>

En este contexto, <math>\vec{u}_\rho</math> es el unitario radial correspondiente a la masa 1 (el de la masa 2 sería el opuesto). Para evitar la posible confusión podemos usar la base cartesiana

<center><math>\vec{u}_\rho = \cos(\varphi)\vec{\imath}+\mathrm{sen}(\varphi)\vec{\jmath}= \cos(\omega t)\vec{\imath}+\mathrm{sen}(\omega t)\vec{\jmath}</math></center>

y queda finalmente

<center><math>\vec{L}_O=\frac{MH^2}{4}\cos(\beta)\left(-\mathrm{sen}(\beta)\cos(\omega t)\vec{\imath}-\mathrm{sen}(\beta)\mathrm{sen}(\omega t)\vec{\jmath}+\cos(\beta)\vec{k}\right)</math></center>

Puesto que de nuevo el centro de masas coincide con el punto del eje donde está unida la varilla

<center><math>\vec{L}_C=\vec{L}_O</math></center>

La energía cinética es más simple:

<center><math>K = K' = 2\left(\frac{m}{2}\omega^2\left(\frac{H}{2}\cos(\beta)\right)^2\right) = \frac{MH^2}{8}\omega^2\cos^2(\beta)</math></center>
==Masas desiguales==
Por último, cuando las masas son desiguales y la barra horizontal y centrada, se verifica

<center><math>M = m_1+m_2\qquad\qquad \vec{r}_2=-\vec{r}_1\qquad\qquad \vec{v}_2=-\vec{v}_1</math></center>

El momento cinético respecto al centro de rotación es de nuevo

<center><math>\vec{L}_O = (m_1+m_2)\vec{r}_1\times\vec{v}_1=\frac{MH^2}{4}\vec{\omega}</math></center>

En este caso, el momento cinético respecto al CM es diferente, ya que al tratarse de masas diferentes, su posición está descentrada. Se halla en

<center><math>\vec{r}_C = \frac{m_1(H/2)+m_2(-H/2)}{m_1+m_2}\vec{u}_\rho = \frac{m_1-m_2}{m_1+m_2}\,\frac{H}{2}\vec{u}_\rho</math></center>

El centro de masas gira con la misma velocidad angular que las dos masas, por lo que el momento cinético debido al movimiento del CM es

<center><math>M\vec{r}_C\times\vec{v}_C = M|\vec{r}|_C^2\vec{\omega} = M\left(\frac{m_1-m_2}{m_1+m_2}\right)^2\frac{H}{2}\vec{\omega}</math></center>

Restamos esta cantidad del momento completo y nos queda

<center><math>\vec{L}_C=\left(1-\frac{(m_1-m_2)^2}{(m_1+m_2)^2}\right)\frac{MH^2}{4}\vec{\omega} =\frac{m_1m_2 H^2}{M}\vec{\omega}</math></center>

De forma análoga tenemos la energía cinética de traslación

<center><math>\frac{M}{2}\left|\vec{v}_C\right|^2 =  \frac{M}{2}\left(\frac{m_1-m_2}{m_1+m_2}\right)^2\frac{H}{2}\omega^2</math></center>

y queda, para la energía cinética respecto al CM,

<center><math>K'=\frac{m_1m_2 H^2}{2M}\omega^2</math></center>




[[Categoría:Problemas de dinámica de los sistemas de partículas (GIE)]]