Ver el código de «Sistema de dos masas unidas por una varilla (CMR)»

==Enunciado==
Un sólido está formado por dos masas iguales,<math> m_1=m_2=m</math>, unidas por una varilla sin masa. En un instante dado, las dos partículas se hallan en <math>\overrightarrow{OA}=2b\vec{\imath}_2</math> y <math>\overrightarrow{OB}=b(3\vec{\imath}_2+2\vec{\jmath}_2 )</math>, respectivamente, siendo el sistema <math>OX_2Y_2Z_2</math> uno ligado al sólido (<math>OX_1Y_1Z_1</math> sería un sistema fijo). 
# ¿Cuál es la posición del CM del sistema?
# ¿Cuánto vale su tensor de inercia respecto al triedro <math>OX_2 Y_2 Z_2</math>?
# Determine los ejes principales de inercia y los momentos principales de inercia de este sólido.
# Si la varilla gira con velocidad angular constante <math>\Omega\vec{\imath}21</math> alrededor del eje <math>OX_1=OX_2</math>, ¿cuánto vale su momento cinético respecto a O? ¿Y su energía cinética?
# Para el movimiento anterior, ¿qué fuerza y que momento hay que aplicar en O para mantener el sistema en movimiento?
<center>[[Archivo:dos-masas-inercia.png|300px]]</center>

==Centro de masas==
Al ser las dos masas iguales, el centro de masas se halla en el punto medio

<center><math>\overrightarrow{OG}=\frac{\overrightarrow{OA}+\overrightarrow{OB}}{2}=\frac{5}{2}b\vec{\imath}_2+b\vec{\jmath}_2</math></center>

==Tensor de inercia==
Para un sistema de masas puntuales, basta con sumar las contribuciones de cada una a los diferentes momentos y productos de inercia.
===Respecto al eje OX===
La masa 1 se halla sobre el propio eje y la 2 a una distancia 2b, por tanto

<center><math>I_{xx}=0+m(2b)^2=4mb^2\,</math></center>
===Respecto al eje OY===
La masa 1 está a 2b y la 2 a 3b

<center><math>I_{yy}=m(4b^2+9b^2)=13mb^2\,</math></center>

===Respecto al eje OZ===
Se trata de una figura plana en el plano <math>OX_1Y_1</math> por lo que

<center><math>I_{zz}=I_{xx}+I_{yy}=17mb^2\,</math></center>

===Productos de inercia===
Por ser una figura plana en <math>OX_1Y_1</math>

<center><math>I_{xz}=I_{yz}=0\,</math></center>

El único elemento no diagonal no nulo es el xy

<center><math>I_{xy}=-\sum_i m_i x_iy_i=-m(2b\cdot 0)-m(3b)(2b)=-6mb^2\,</math></center>

===Forma matricial===
Reuniendo los resultados anteriores queda la matriz

<center><math>\bar{\bar{I}}=mb^2\begin{pmatrix}4 & -6 & 0 \\ -6 & 13 & 0 \\ 0 & 0 & 17\end{pmatrix}</math></center>

También se puede calcular esta matriz directamente en lugar de elemento a elemento, partiendo de la expresión

<center><math>\bar{\bar{I}}=\sum_i m_i\begin{pmatrix} y_i^2+z_1^2 & -x_iy_i & -x_iz_i \\ -x_iy_i & x_i^2+z_i^2 & -y_iz_i \\ -x_iz_i & -y_iz_i & x_i^2+y_i^2 \end{pmatrix}</math></center>

que en este caso nos da

<center><math>\bar{\bar{I}}=mb^2\begin{pmatrix} 0 & 0 & 0 \\ 0 & 4 & 0 \\ 0 & 0 & 4 \end{pmatrix}+mb^2\begin{pmatrix} 4 & -6 & 0 \\ -6 & 9 & 0 \\ 0 & 0 & 13 \end{pmatrix}=mb^2\begin{pmatrix}4 & -6 & 0 \\ -6 & 13 & 0 \\ 0 & 0 & 17\end{pmatrix}</math></center>
==Ejes y momentos principales==
Los momentos principales de inercia resultan de resolver el determinante

<center><math>\left|\begin{matrix}4-I & -6 & 0 \\ -6 & 13-I & 0 \\ 0 & 0 & 17-I \end{matrix}\right|=0</math></center>

(donde en el resultado final habrá que multiplicar el número que salga por <math>mb^2</math>). Desarrollamos el determinante

<center><math>\Delta = (17-I)(I^2-17I+16)=0\,</math></center>

Uno de los autovalores es inmediato, <math>I_{ZZ}=17mb^2</math>. Esto era de esperar pues ya sabíamos que el eje OZ era principal. Los otros dos salen de una ecuación de segundo grado

<center><math>I_{XX}=mb^2 \qquad\qquad I_{YY}=16mb^2</math></center>

Para hallar cada autovector resolvemos la ecuación vectorial. Para <math>I_{ZZ}</math> es inmediato

<center><math>mb^2\begin{pmatrix}4 & -6 & 0 \\ -6 & 13 & 0 \\ 0 & 0 & 17 \end{pmatrix}\cdot\begin{pmatrix}u_x \\u_y\\ u_z\end{pmatrix}=17mb^2\begin{pmatrix}u_x \\u_y\\ u_z\end{pmatrix}</math></center>

que da el sistema

<center><math>4u_x-6u_y=17u_x\qquad\qquad -6u_x+13u_y=17u_y \qquad\qquad 17u_z=17u_z</math></center>

con solución 

<center><math>u_x=0,\quad u_y=0\,\quad u_z=1\qquad\qquad\vec{k}_3=\vec{k}_2</math></center>

Para el autovalor <math>I_{XX}</math> operamos de la misma forma

<center><math>mb^2\begin{pmatrix}4 & -6 & 0 \\ -6 & 13 & 0 \\ 0 & 0 & 17 \end{pmatrix}\cdot\begin{pmatrix}u_x \\u_y\\ u_z\end{pmatrix}=mb^2\begin{pmatrix}u_x \\u_y\\ u_z\end{pmatrix}</math></center>

que da el sistema

<center><math>4u_x-6u_y=u_x\qquad\qquad -6u_x+13u_y=u_y \qquad\qquad 17u_z=u_z</math></center>

con solución 

<center><math>u_x=2,\quad u_y=1\quad u_z=0\qquad\qquad\vec{\imath}_3=\frac{2}{\sqrt{5}}\vec{\imath}_2+\frac{1}{\sqrt{5}}\vec{\jmath}_2</math></center>

donde en el última paso se ha normalizado el vector. 

Por último, para el autovalor <math>I_{YY}=16mb^2</math> resulta el tercer vector del triedro
<center><math>
\vec{\jmath}_3=-\frac{1}{\sqrt{5}}\vec{\imath}_2+\frac{2}{\sqrt{5}}\vec{\jmath}_2</math></center>

En esta base, el tensor de inercia se escribe

<center><math>\bar{\bar{I}}=mb^2\begin{pmatrix}1 & 0 & 0 \\ 0 & 16 & 0 \\ 0 & 0 & 17 \end{pmatrix}</math></center>
==Momento y energía cinética==
===Momento cinético===
Para el momento cinético no necesitamos los momentos o ejes principales (aunque en problemas de dinámica donde haya que resolver las ecuaciones de Euler si son muy útiles). Nos basta con multiplicar el tensor de inercia por la velocidad angular

<center><math>\vec{L}_O=\bar{\bar{I}}\cdot\vec{\omega}=mb^2\begin{pmatrix}4 & -6 & 0 \\ -6 & 13 & 0 \\ 0 & 0 & 17\end{pmatrix}\cdot\begin{pmatrix}\Omega \\ 0 \\ 0 \end{pmatrix}=mb^2\Omega\begin{pmatrix}4 \\ -6 \\ 0 \end{pmatrix}</math></center>

o, en forma vectorial,

<center><math>\vec{L}_O=mb^2\Omega\left(4\vec{\imath}_2-6\vec{\jmath}_2\right)</math></center>

Por no ser <math>OX_2</math> un eje principal, el momento cinético no es paralelo a la velocidad angular.

Nótese que el momento cinético es constante en el sistema de referencia 2, ligado al sólido, pero no en el 1, fijo. A medida que el sólido gira, el vector momento cinético lo acompaña, describiendo un cono. 
===Energía cinética===
Conocido el momento cinético, la energía cinética se calcula como

<center><math>T=\frac{1}{2}\vec{\omega}\cdot\vec{L}_O=2mb^2\Omega</math></center>

Aquí es importante el hecho de que O es un punto fijo respecto a este movimiento. Si no fuera así la expresión sería más complicada.
==Fuerza y par==
===Fuerza===
Por la ley de Newton

<center><math>\vec{F}=m_T\vec{a}_G</math></center>

aquí <math>m_T</math> es la masa total

<center><math>m_T=2m\,</math></center>

y <math>\vec{a}_G</math> es la aceleración del CM respecto a un sistema fijo.

El CM está describiendo un movimiento circular uniforme de radio <math>b</math>, por lo que su aceleración es puramente normal, dirigida hacia el eje de rotación (el OX)

<center><math>\vec{a}_G=-\Omega^2 b\vec{\jmath}_2</math></center>

y

<center><math>\vec{F}_O=-2mb\Omega^2\vec{\jmath}_2</math></center>

De manera más formal y general, podemos introducir la velocidad angular relativa entre los dos sólidos (que no será otra que la de rotación)

<center><math>\vec{\omega}_{21}=\Omega\vec{\imath}_1=\Omega\vec{\imath}_2</math></center>

El centro de masas está estacionario respecto al sistema de referencia 2. Por ello, su aceleración es 

<center><math>\vec{a}^G_{21}=\vec{a}^O_{21}+\vec{\alpha}_{21}\times \overrightarrow{OG}+\vec{\omega}_{21}\times\left(\vec{\omega}_{21}\times\overrightarrow{OG}\right)</math></center>

El punto O es fijo, por lo que su aceleración es nula. También se anula la aceleración angular, por lo que nos queda

<center><math>\vec{a}^G_{21}=(\Omega\vec{\imath}_2)\times \left((\Omega\vec{\imath}_2) \times \left(\frac{5}{2}b\vec{\imath}_2+b\vec{\jmath}_2\right)\right)=(\Omega\vec{\imath}_2)\times \left(\Omega b\vec{k}_2\right)=-\Omega^2b\vec{\jmath}_2</math></center>

y de aquí llegamos a la fuerza

<center><math>\vec{F}_O=-2mb\Omega^2\vec{\jmath}_2</math></center>

Nótese que, de nuevo, esta fuerza es constante en el sistema 2, pero en un sistema fijo la dirección de la fuerza va girando.

===Momento de las fuerzas===
Por ser O un punto fijo

<center><math>\vec{M}_O=\left.\frac{\mathrm{d}\vec{L}_O}{\mathrm{d}t}\right|_1</math></center>

En este caso, el momento cinético acompaña al sólido en su rotación, de forma que en el sistema 2 es constante

<center><math>\left.\frac{\mathrm{d}\vec{L}_O}{\mathrm{d}t}\right|_2=\vec{0}</math></center>

por lo que, según la fórmula de Poisson

<center><math>\vec{M}_O=\overbrace{\left.\frac{\mathrm{d}\vec{L}_O}{\mathrm{d}t}\right|_2}^{=\vec{0}}+\vec{\omega}_{21}\times\vec{L}_O=(\Omega\vec{\imath}_1)\times \left(mb^2\Omega\left(4\vec{\imath}_1-6\vec{\jmath}_1\right)\right)=-6mb^2\Omega\vec{k}_1</math></center>
[[Categoría:Problemas de dinámica del sólido rígido (CMR)]]