Ir al contenido principal

¿Cómo hallar el centro de masa de una lámina con integrales dobles? 2

 A petición de un suscriptor y seguidor voy a mostrar dos ejercicios para el cálculo del centro de masa con integrales dobles, comencemos:

El segundo ejercicio es:

2. Calcular el centro de masa de una lámina representada por la región $R$ que se encuentra por encima del eje $x$ y entre las líneas

$y=x$; $y=-x$, $x^{2}+y^{2}=4y$; $x^{2}+y^{2}=6y$; $y>0$; donde la densidad es $\sqrt{x^{2}+y^{2}}$

Primero vamos a representar el área, al cuál vamos a hallar el centro de masa (voy a utilizar Geogebra).

Este problema resulta más simple si lo resolvemos mediante coordenadas polares, así que la densidad $\rho(x,y)$ en coordenadas polares de acuerdo a las reglas de transformación $x=rcos\theta$ y $y=rsen\theta$ es:
\[\rho(x,y)=\rho(rcos\theta,rsen\theta)=\sqrt{x^{2}+y^{2}}=\sqrt{(rcos\theta)^{2}+(rsen\theta)^{2}}=\sqrt{r^{2}}=r\]
Como vimos en el primer ejercicio, las integrales que tiene cada punto del centro de masa se pueden representar en coordenadas polares, escogemos este tipo de coordenadas porque las rectas tienen relación con los ángulos $\frac{\pi}{4}$ y $\frac{3\pi}{4}$, esto por la coordenada $\theta$, respecto a la coordenada $r$, vamos a ver los límites de integración respecto a r, por comparación de las funciones:
\[x^{2}+y^{2}=4y \quad \wedge \quad x^{2}+y^{2}=6y\]
Igualamos a cero y pasamos a coordenadas polares:
\[x^{2}+y^{2}-4y=0 \quad \wedge \quad x^{2}+y^{2}-6y=0\]
\[(rcos\theta)^{2}+(rsen\theta)^{2}-4rsen\theta=0\quad \wedge \quad(rcos\theta)^{2}+(rsen\theta)^{2}-6rsen\theta=0 \]
\[r^{2}-4rsen\theta=0\quad \wedge \quad r^{2}-6rsen\theta=0\]
Por lo tanto:
\[r^{2}-4rsen\theta<r^{2}-6rsen\theta\]
\[4rsen\theta<6rsen\theta\]
\[4<6\]
\[2<3\]
Es a partir de está ultima desigualdad que integramos desde $r^{2}-4rsen\theta$ hasta $r^{2}-6rsen\theta$, que está de acuerdo con la condición $y>0$, y corresponde a nuestra área en verde, así la integral respecto a la densidad en coordenadas polares queda de la forma:
\[\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}\int_{r^{2}-4rsen\theta}^{r^{2}-6rsen\theta}r^{2}drd\theta\]
Integrando respecto a $r$:
\[\frac{1}{3}\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}[(r^{2}-4rsen\theta)^{3}-(r^{2}-4rsen\theta)^{3}]d\theta\]
Para integrar respecto a theta, vamos a encontrar integrales del tipo $\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}sen\theta d\theta$, $\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}sen^{2}\theta d\theta$, y $\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}sen^{3}\theta d\theta$, por lo tanto obtenemos el siguiente resultado para la integral doble respecto a la densidad:
\[\int\int_R \rho(x,y)dxdy=r^{3}\left(-2\sqrt{2}r^{2}-5(\pi+2)r-\frac{76\sqrt{2}}{9}\right)\]
Ahora calculamos:
\[\int\int_R x\rho(x,y)dxdy=\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}\int_{r^{2}-4rsen\theta}^{r^{2}-6rsen\theta}r^{3}cos\theta drd\theta\]
Integrando nuevamente respecto a r:
\[\frac{1}{4}\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}[(r^{2}-4rsen\theta)^{4}cos\theta-(r^{2}-4rsen\theta)^{3}cos\theta]d\theta\]
Sí hacemos un cambio de variable para $u=r^{2}-6rsen\theta$, o para $u=r^{2}-4rsen\theta$, realizando los cambios necesarios, nos damos cuenta que el valor de la integral respecto a $\theta$ es cero:
\[\frac{1}{4}\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}[(r^{2}-4rsen\theta)^{4}cos\theta-(r^{2}-4rsen\theta)^{3}cos\theta]d\theta=0\]
Ahora realizamos la siguiente y última doble integral:
\[\int\int_R y\rho(x,y)dxdy=\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}\int_{r^{2}-4rsen\theta}^{r^{2}-6rsen\theta}r^{3}sen\theta drd\theta\]
Integramos respecto a $r$:
\[\frac{1}{4}\int_{\frac{\pi}{4}}^{\frac{3\pi}{4}}[(r^{2}-4rsen\theta)^{4}sen\theta-(r^{2}-4rsen\theta)^{3}sen\theta]d\theta\]
\[\int\int_R y\rho(x,y)dxdy=r^{4}\left(-\left(\frac{\pi+2}{2}\right)r^{3}+25\sqrt{2}r^{2}-\frac{19}{8}(3\pi+8)r+\frac{1456\sqrt{2}}{3}\right)\]
Finalmente podemos hallar los puntos donde se encuentran el centro de masa de nuestra lámina:
\[\bar{x}=0\quad \wedge \quad \bar{y}=\frac{r\left(-\left(\frac{\pi+2}{2}\right)r^{3}+25\sqrt{2}r^{2}-\frac{19}{8}(3\pi+8)r+\frac{1456\sqrt{2}}{3}\right)}{\left(-2\sqrt{2}r^{2}-5(\pi+2)r-\frac{76\sqrt{2}}{9}\right)}\]


Comentarios

Entradas populares de este blog

¿Cómo resolver la ecuación diferencial de un oscilador armónico simple con gravedad?

En este caso vamos a solucionar la ecuación diferencial para un oscilador armónico unidimensional que solo tiene el movimiento en el eje y, y además está en presencia de la gravedad, su ecuación diferencial es la siguiente: \[\ddot{y}+\omega^{2}y=g\] Esta es una ecuación diferencial no homogénea y la vamos a resolver por el método de variación de parámetros. Para este método necesitamos resolver la ecuación diferencial homogenea: \[\ddot{y}+\omega^{2}y=0\] Que ya la hemos resuelto en este enlace  aunque resuelta con la coordenada $x$, pero eso no va a importar, porque va a ser la misma solución, solo que cambiamos de $x$ a $y$: \[y(t)=c_{1}cos(\omega t)+c_{2}sen(\omega t)\] Para resolver la ecuación diferencial debemos hallar el wronskiano: \[W(f(t),g(t))=\begin{vmatrix}f(t) & g(t)\\f'(t) & g'(t)\end{vmatrix}\] Identificamos las funciones $f(t)=c_{1}cos(\omega t)$ y $g(t)=c_{2}sen(\omega t)$, hallamos las derivadas $f'(t)=-c_{1}\omega sen(\omega t)$ y $g'(t)=c_{...

¿Cuál es la integral de sen(2x)?

Nuestra integral a resolver está vez es la siguiente: \[\int sen(2x)dx\] Podemos realizar la siguiente sustitución $u=2x$, $du=2dx$, entonces $\frac{du}{2}=dx$, así la integral nos queda de la forma: \[\frac{1}{2}\int sen(u)du\] Que corresponde a una integral fundamental: \[\frac{1}{2}\int sen(u)du=-\frac{1}{2}cos(u)+C\] Deshacemos el cambio de variable y obtenemos la respuesta a nuestra integral: \[-\frac{1}{2}cos(2x)+C\]

¿Cómo resolver la ecuación diferencial de un oscilador armónico simple?

En este caso vamos a solucionar la ecuación diferencial para un oscilador armónico unidimensional que solo tiene el movimiento en el eje x, y su ecuación diferencial es la siguiente: \[\ddot{x}+\omega^{2}x=0\] Podemos ver que la ecuación corresponde a una ecuación diferencial homogénea, así que proponemos una solución de la forma $x=e^{mt}$ y hallamos sus derivadas $\dot{x}=me^{mt}$, $\ddot{x}=m^{2}e^{mt}$, reemplazamos en nuestra ecuación diferencial: \[m^{2}e^{mt}+\omega^{2}e^{mt}=0\] Sacando factor común $e^{mt}$: \[e^{mt}(m^{2}+\omega^{2})=0\] Donde $e^{mt}$ no puede ser cero, entonces lo será el termino entre parentesis y hallamos sus respectivas raíces: \[m^{2}+\omega^{2}=0\] Donde las soluciones son imaginarias y corresponden a: \[m_1=i\omega \quad m-2=-i\omega\] La solución queda expresada como: \[x(t)=e^{i\omega t}+e^{-i\omega t}\] Aplicamos la identidad de Euler para simplificar mas elegante esta solución: \[e^{i\theta}=cos\theta+isen\theta\] Nos queda ahora ...