Ver el código de «Centro de masas de sistemas continuos»

== Enunciado ==
Calcula la posición del centro de masas de estos sistemas
#Una barra homogénea delgada de longitud <math>a</math> y masa <math>m</math>.
#Una barra de longitud <math>a</math> y densidad lineal de masa <math>\lambda = Cx</math>, siendo <math>x</math> la distancia a un extremo de la barra y <math>C</math> una constante.
#Una barra homogénea delgada en forma de semicírculo de radio <math>a</math> y masa <math>m</math>.
#Dos esferas macizas de masas <math>m_1</math> y <math>m_2</math> y radios <math>R</math> y <math>2R</math>, unidas por un cilindro de masa <math>m_3</math> y longitud <math>L</math>.
#Una esfera maciza de radio <math>2b</math>, densidad homogénea <math>\rho_0</math>, con una cavidad también esférica, de radio <math>b</math>, cuyo centro se encuentra a una distancia <math>b</math> del centro de la esfera original.

== Solución ==
Para un sistema discreto la posición del centro de masas (CM) viene dada por la expresión
<center>
<math>
\displaystyle \vec{r}_{CM} = \frac{\sum\limits_{i=1}^n m_i\vec{r}_i}{\sum\limits_{i=1}^n m_i}
</math>
</center>
donde <math>m_i </math> es la masa de cada partícula y <math>\vec{r}_i </math> su vector de posición.
Cuando tratamos con un sistema continua, la expresión se transforma según el cambio
<center>
<math>
m_i\to\mathrm{d}m \qquad \qquad \qquad \vec{r}_i \to \vec{r}
</math>
</center>
Así, en un sistema continuo la posición del centro de masas viene dada por la expresión
<center>
<math>
\displaystyle \vec{r}_{CM} = \frac{\int \vec{r}\,\mathrm{d}m}{\int \mathrm{d}m}
</math>
</center>
siendo <math>\vec{r} </math> un vector que recorre cada uno de los puntos del sistema y <math>\mathrm{d}m </math> la masa infinitesimal asociada a cada uno de esos puntos. 

=== Barra recta homogénea ===
[[Imagen:Centro_masas_barra_recta.png|right]]
Consideramos el caso de una barra homogénea delgada de masa <math>m </math> y longitud <math>a </math>. Lo primero que hay que hacer es escoger un sistema de ejes para describir la posición de cada punto de la barra. Elegimos el eje  <math>OX </math> de modo que coincida con la barra y situamos el origen en su extremo izquierdo. Con esta elección la posición de un punto genérico de la barra viene dada por el vector de posición
<center>
<math>
\vec{r} = x\,\vec{\imath} \qquad \mathrm{con} \qquad 0<x<a
</math>
</center>
La variable <math>x </math> es la etiqueta que identifica a cada punto de la barra.

Ahora consideramos que en cada punto de la barra hay un pequeño trocito de barra de longitud <math>\mathrm{d}x </math> y masa <math>\mathrm{d}m </math>. La d delante de la <math>x </math> y la <math>m </math> sólo significa que la longitud del elemento y su masa son muy pequeñas. ¿Cuanto vale está masa?. 

Como la barra es homogénea, podemos definir una '''densidad lineal de masa''' como el cociente de su masa total por su longitud
<center>
<math>
\lambda = \frac{m}{a}
</math>
</center>
Con esto, si el trocito de barra tiene una longitud <math>\mathrm{d}x </math>, su masa es
<center>
<math>
\mathrm{d}m = \lambda\mathrm{d}x = \frac{m}{a}\,\mathrm{d}x
</math>
</center>
Ahora podemos calcular las integrales en la expresión de <math>\vec{r}_{CM} </math>. La integral en el denominador es la suma de
las masas de todos los puntos que podemos considerar en la barra, esto es, su masa completa
<center>
<math>
\int\mathrm{d}m = \int\limits_0^a\lambda\,\mathrm{d}x = \lambda\int\limits_0^a\mathrm{d}x = \lambda\,a = m
</math>
</center>
La integral en el numerador es
<center>
<math>
\int\vec{r}\,\mathrm{d}m = \int\limits_0^a (x\,\vec{\imath})\,\lambda\,\mathrm{d}x = \vec{\imath}\lambda\int\limits_0^a x\,\mathrm{d}x
=\vec{\imath}\,\lambda \left[\frac{x^2}{2}\right]_0^a = \frac{1}{2}\lambda\,a^2 = \frac{1}{2}m\,a
</math>
</center>
El vector <math>\vec{\imath} </math> y la densidad de masa <math>\lambda </math> pueden salir de la integral pues no dependen de 
<math>x </math>, es decir, son iguales no importa en que punto de la barra estemos.

Con esto el vector de posición del CM de la barra respecto de su extremo izquierdo es
<center>
<math>
\vec{r}_{CM} = \frac{ma/2}{m}\vec{\imath} = \frac{1}{2}a\,\vec{\imath}
</math>
</center>
Es decir, para encontrar el CM nos situamos en el extremo izquierdo de la barra y nos desplazamos hacia la derecha una longitud 
igual a la mitad de su longitud. El CM se sitúa en el centro de la barra. Esto es lógico, pues los ejes de simetría de la barra 
pasan todos por su centro, por lo que el CM debe situarse en él.

=== Barra recta con densidad variable ===
[[Imagen:Centro_masas_barra_recta.png|right]]
El procedimiento es similar al del apartado anterior. La posición de un elemento de longitud en la barra está dada por el vector
<center>
<math>
\vec{r} = x\,\vec{\imath} \qquad \mathrm{con} \qquad 0<x<a
</math>
</center>
La diferencia es que ahora la masa de cada elemento de longitud es
<center>
<math>
\mathrm{d}m = \lambda\mathrm{d}x = Cx\mathrm{d}x
</math>
</center>
La masa total de la barra es la suma de las masas de todos los elementos de longitud
<center>
<math>
m = \int\mathrm{d}m = \int\limits_0^aCx\,\mathrm{d}x = \dfrac{1}{2}Ca^2.
</math>
</center>

La integral en el numerador de la expresión que da el vector de posición del centro de masas es
<center>
<math>
\int\vec{r}\,\mathrm{d}m = \int\limits_0^a (x\,\vec{\imath})\,Cx\,\mathrm{d}x = \vec{\imath}\,C\int\limits_0^a x^2\,\mathrm{d}x
=\vec{\imath}\,C\, \left[\dfrac{x^3}{3}\right]_0^a = \frac{1}{3}C\,a^3\,\vec{\imath}
</math>
</center>
El vector <math>\vec{\imath} </math> y la constante <math>C </math> pueden salir de la integral pues no dependen de <math>x </math>, es decir, son iguales no importa en que punto de la barra estemos.

Con esto el vector de posición del CM de la barra respecto de su extremo izquierdo es
<center>
<math>
\vec{r}_{CM} = \frac{Ca^3/3}{Ca^2/2}\vec{\imath} = \frac{2}{3}a\,\vec{\imath}
</math>
</center>
Es decir, para encontrar el CM nos situamos en el extremo izquierdo de la barra y nos desplazamos hacia la derecha una longitud igual a dos tercios de su longitud. Ahora el CM se desplaza hacia la derecha porque la densidad de masa es mayor para valores mayores de <math>x</math>.

=== Barra semicircular ===

[[Imagen:Centro_masas_barra_circular.png|right]]
Consideramos ahora el caso en que la barra tiene forma de semicírculo. De nuevo, consideramos pequeños elementos de línea a lo largo de
la barra. Escogemos el origen del sistema de coordenadas en el centro del semicírculo, de modo que el eje <math>OX </math> pase por los dos
extremos de la semicircunferencia. Con estos ejes, la posición de un punto de la barra queda definida por un valor del ángulo <math>\theta </math>
<center>
<math>
\vec{r}(\theta) = a\cos\theta\,\vec{\imath} + a\,\mathrm{sen}\theta\,\vec{\jmath}
\qquad \qquad \mathrm{con} \qquad 0\leq\theta\leq\pi
</math>
</center>
Cada elemento de línea tiene una longitud
<center>
<math>
\mathrm{d}l = a\,\mathrm{d}\theta
</math>
</center>
Como la barra es homogénea su densidad de masa es uniforme e igual a su masa dividida por su longitud
<center>
<math>
\lambda = \frac{m}{\pi\, a}
</math>
</center>
Con esto, la masa de cada elemento de línea es
<center>
<math>
\mathrm{d}m = \lambda\,\mathrm{d}l = \lambda\,a\,\mathrm{d}\theta
</math>
</center>
Podemos calcular la posición del centro de masas de la barra usando la expresión del apartado anterior. El numerador es
<center>
<math>
\int\vec{r}\,\mathrm{d}m = \int\limits_0^{\pi}(a\cos\theta\,\vec{\imath} + a\,\mathrm{sen}\theta\,\vec{\jmath}) \lambda\,a\mathrm{d}\theta=
\lambda\,a^2\left[ \mathrm{sen}\theta\,\vec{\imath} - \cos\theta\,\vec{\jmath}\right]_0^{\pi} =
2\lambda\,a^2\vec{\jmath}
</math>
</center>
Sustituyendo el valor de <math>\lambda </math> obtenemos
<center>
<math>
\int\vec{r}\,\mathrm{d}m = \frac{2}{\pi}a\,m\,\vec{\jmath}
</math>
</center>
El vector de posición del centro de masas es
<center>
<math>
\displaystyle \vec{r}_{CM} = \frac{\int\vec{r}\,\mathrm{d}m}{\int\mathrm{d}m} = \frac{\frac{2}{\pi}a\,m}{m}\vec{\jmath}
=
\frac{2}{\pi}a\,\vec{\jmath}
</math>
</center>
Debido a la simetría, el CM está en el diámetro vertical de la semicircunferencia. Como <math>(2/\pi)=0.637 </math>, el CM está por debajo de 
la semicircunferencia, como se indica en la figura

=== Dos esferas con cilindro ===
[[Archivo:F1-CM-EsferasCilindro.png|center]]
Usaremos la técnica de asociación (o descomposición) para resolver este caso. Para calcular la posición del CM de un sistema continuo hay que hacer una integral definida. Recordemos que, en el fondo, una integral definida es una suma. Entonces se puede aplicar la propiedad distributiva de la suma. Esto significa que podemos dividir el sistema completo en subsistemas, y substituir cada subsistema por una masa puntual situada en el CM del subsistema y cuya masa sea igual a la masa total del subsistema.

En el caso del sistema que nos ocupa, cada una de las esferas se sustituye por una masa puntual situada en su centro. El cilindro se sustituye por una masa puntual situada en su centro. Sabemos que los CM de los subsistemas tienen que estar en esas posiciones por la simetría de la esfera y el cilindro.

Como se observa en la figura, escogemos el eje <math>X</math> de modo que coincida con el eje de simetría del sistema y el origen <math>O</math> coincidente con el centro de la esfera de la izquierda. Con esta elección el vector de posición de cada masa puntual es
<center>
<math>
\begin{array}{l}
\vec{r}_1 = (3R + L)\,\vec{\imath},\\
\vec{r}_2 = \vec{0},\\
\vec{r}_3 = (R+L/2)\,\vec{\imath}.
\end{array}
</math>
</center>
Por tanto, el vector de posición del CM del sistema en el sistema de referencia elegido es
<center>
<math>
\vec{r}_{CM} = 
\dfrac{m_1\vec{r}_1 + m_2\vec{r}_2 + m_3\vec{r}_3}{m_1+m_2+m_3}
=
\left(
(3R+L)m_1 + \left(R + \dfrac{1}{2}L\right)m_3
\right)\,\vec{\imath}.
</math>
</center>
La posición aproximada del CM se muestra también en la figura.

=== Esfera con hueco esférico ===
[[Archivo:F1-CM-EsferaHueco.png|right]]
Aplicamos la misma técnica para resolver esta configuración, aunque usando un pequeño truco. El sistema original tiene un eje de simetría que coincide con el eje <math>X</math> mostrado en la figura. Por tanto, el CM tiene que estar en esa línea. Y probablemente estará a la izquierda de <math>O</math>, pues hay mas masa hacia la izquierda. Pero no podemos saber nada mas. La idea es sustituir el sistema por dos esferas homogéneas completas (sin huecos), de radios <math>2b</math> y <math>b</math>, y con densidades volumétricas de masa <math>\rho_1</math> y <math>\rho_2</math>. El truco está en el valor asignado a estas densidades de masa, a saber
<center>
<math>
\rho_1=\rho_0, \qquad \rho_2=-\rho_0.
</math>
</center>
Puede resultar sorprendente considerar una esfera de masa negativa. La clave es darse de cuenta que este sistema con dos esferas da el mismo valor para la masa de un pequeño elemento de volumen que el sistema original. Para el elemento de volumen marrón de la figura tendríamos
<center>
<math>
\mathrm{d}m = \rho_1\,\mathrm{d}V = \rho_0\,\mathrm{d}V,
</math>
</center>
que es el valor que debe tener. Si es el elemento de volumen verde (en el hueco del sistema original) tendríamos
<center>
<math>
\mathrm{d}m = \rho_1\,\mathrm{d}V + \rho_2\,\mathrm{d}V = (\rho_0-\rho_0)\,\mathrm{d}V = 0,
</math>
</center>
que es el valor correcto para el hueco.

A continuación, cada esfera completa homogénea se sustituye por una masa puntual en su centro, por simetría. El valor de cada masa se obtiene multiplicando su densidad volumétrica de masa por su volumen
<center>
<math>
\begin{array}{l}
m_1 = \rho_1\,V_1 = \rho_0 \dfrac{4\pi}{3}(2b)^3 = \dfrac{32\pi}{3}\rho_0b^3,\\
\\
m_2 = \rho_2\,V_2 = -\rho_0 \dfrac{4\pi}{3}b^3 = -\dfrac{4\pi}{3}\rho_0b^3.
\end{array}
</math>
</center>
Con el eje <math>X</math> elegido en la figura el vector de posición de cada masa puntual es
<center>
<math>
\vec{r}_1 = \vec{0}, \qquad \vec{r}_2 = b\,\vec{\imath}.
</math>
</center>
Entonces, el vector de posición del CM del sistema original es
<center>
<math>
\vec{r}_{CM} = \dfrac{m_1\,\vec{r}_1+m_2\,\vec{r}_2}{m_1+m_2}
=
-\dfrac{4\pi\rho_0b^4/3}{28\pi\rho_0b^3/3}\,\vec{\imath}
=
-\dfrac{1}{7}\,b\,\vec{\imath}.
</math>
</center>
Como esperábamos, el CM del sistema está a la izquierda del punto <math>O</math>.

[[Categoría:Problemas de dinámica de un sistema de partículas|0]]
[[Categoría:Dinámica de un sistema de partículas|1]]