Tp.8 Ej.15.b

Calcule la masa de los siguientes cuerpos:

b) cuerpo definido por z \geq |y|, x^2+y^2+z^2 \leq 1 si la densidad en cada punto es proporcional a la distancia desde el punto al plano xy

Solución:

La función densidad es \delta(x,y,z) = k|z|, y como z \geq 0 nos queda:
\delta(x,y,z) = kz

Como la esfera esta cortada por dos planos, voy a tomar una variante de coordenadas esféricas (que sería equivalente a rotar x \to y, y \to z, z\to x):

\begin{cases} x = \rho\cos(\beta) \\ y = \rho\cos(\phi)\sin(\beta) \\ z = \rho\sin(\phi)\sin(\beta) \end{cases}

|J| = \rho^2 \sin(\beta)

Transformo la función densidad:
\delta(\phi, \beta, \rho) = k \rho \sin(\phi)\sin(\beta)

Por lo tanto la masa es:
k \int_{\pi/4}^{3\pi/4} \sin(\phi) d\phi \int_0^{\pi} \sin^2(\beta) d\beta \int_0^1 \rho^3 d\rho

\frac{k}{4} [-\cos(\phi)]_{\pi/4}^{3\pi/4} [\frac{\beta}{2} - \frac{\sin(2\beta)}{4}]_{0}^{\pi}

\frac{k}{4} [\frac{\sqrt{2}}{2} + \frac{\sqrt{2}}{2}] [\frac{\pi}{2}]

= \frac{k}{8} \sqrt{2} \pi

El siguiente es el gráfico del cuerpo:


draw3d(surface_hide=true,
xlabel = "x", ylabel = "y", zlabel = "z",
color=blue,
parametric_surface(cos(b), cos(p)*sin(b), sin(p)*sin(b), p,%pi/4,3*%pi/4, b,0,%pi),
color=light-blue,
parametric_surface(u*cos(b), u*sqrt(2)/2*sin(b), u*sqrt(2)/2*sin(b), u,0,1, b,0,%pi),
parametric_surface(u*cos(b), -u*sqrt(2)/2*sin(b), u*sqrt(2)/2*sin(b), u,0,1, b,0,%pi)
);

Esta entrada fue publicada en TP08 - Integrales Múltiples. Guarda el enlace permanente.

11 respuestas a Tp.8 Ej.15.b

  1. Gabriel dijo:

    Damian, yo podria proyectar ese cuerpo en el plano XY?
    Intente calcular el volumen con cilindricas del lado de la Y positiva (angulo de -pi/2 a pi/2) y luego multiplicar por 2 pero no me da.

    En caso de no poder, no entiendo bien el criterio para elegir el plano al cual proyectar en casos como este

    gracias
    gabriel

    • damidami dijo:

      Hola Gabriel,
      En este caso te conviene usar coordenadas esféricas como en la resolución, en vez de cilíndricas, a pesar de ser proyectable en el plano xy
      Por otro lado, tené en cuenta que el ejercicio pide calcular la masa y no el volúmen, con lo cual no te conviene usar propiedades de simetría ya que además del cuerpo tendría que ser simétrica la función densidad.
      Saludos,
      Damián

  2. Aoshido dijo:

    Hola que tal?
    Yo lo proyecte pero en el plano ZY y me quedaria algo asi : (><). Estando en el medio el eje Z, entonces se podrian aplicar las propiedades de simetria para la masa, (se me ocurrió, ya que varia igual para ambos lados).
    El tema es que no me avive de usar esfericas y usé cilindricas. Los limites me quedaron asi: (siendo mis polares: z= ρ sen (φ) y=ρ cos (φ) ):

    D1: 7/4π<=φ<=1/4π
    0<=ρ<=1
    -(1-ρ^2)^1/2<= x <= (1-ρ^2)^1/2

    D2:3/4π<=φ<=5/4π
    0<=ρ<=1
    -(1-ρ^2)^1/2<= x <= (1-ρ^2)^1/2

    Y Calcule la primera integral triple (D1) y me termino dando (Dps de una integral de tablas bastante fea) (k3π/16) multiplicado por 2 (por la otra semi esfera) me da (K3π/8) MUY cerca del resultado. Tengo mal algun limite?
    Gracias, Aoshido.

    Pd: (en el grafico dibujaste 1 de las "semiesferas" pero son las 2 no?)

    • dami dijo:

      Hola Aoshido,
      No entiendo del todo como te quedó tu región de integración.
      Por un lado, a mi me queda una sola región, como si fuera un gajo de una mandarina, en el dibujo se ve bien, y no dos semiesferas.
      Por otro lado, en tu región D1 ponés \frac{7}{4}\pi \leq \phi \leq \frac{1}{4}\pi y ningún \phi real satisface eso.
      De todas formas proyectar sobre el plano yz si bien no es imposible no parece ser una forma sencilla de resolver este ejercicio. Me parece que en tu caso te debería quedar una región sola, y el ángulo ir de \frac{\pi}{4} \leq \phi \leq \frac{3\pi}{4}, habría que hacerlo y ver si da lo mismo.
      Suerte,
      Damián.

  3. Aoshido dijo:

    Hola Damian,gracias por tu pronta respuesta.
    Si, me exprese mal, no quise decir 2 semi-esferas, si no 2 gajos. Ya resolvi el problema, tristemente habia abierto mal el modulo, entonces me quedaban los dos «gajos» de los costados jajaja, no se como me dio tan cerca el resultado, muchas gracias por la ayuda!
    Ah por cierto, es verdad lo que decis de (7/4π<=φ<=1/4π), cuando lo defini no me di cuenta, pero entendiste que región quise definir al menos? jajaja.
    D1 (mal definido) era el gajo del costado derecho, pero creo que para quedara bien definido (creo) tuve que haberlo expresado como 7/4π<=φ<=9/4π.
    Devuelta muchas gracias por la ayuda!
    Saludos!

    • dami dijo:

      La verdad que no entendí como lo resolvistes. Supongo que si te da el resultado correcto las chances que esté bien son buenas, pero igual ojo cuando aplicás propiedades de simetría, yo nunca lo recomiendo cuando lo que se calcula es masa porque no alcanza con que el cuerpo sea simétrico, y no cuesta nada hacer la integral completa.
      No me lograron cerrar los ángulos que tomastes, \frac{9}{4}\pi > 2\pi lo cual es medio extraño.

      En cilíndricas sobre el eje x te debería quedar algo así como

      M = k \int_{\frac{\pi}{4}}^{\frac{3\pi}{4}} \sin(\phi) d\phi \int_0^1 \rho^2 d\rho \int_{-\sqrt{1-\rho^2}}^{+\sqrt{1-\rho^2}} dx

      donde la densidad quedó como
      \delta = k \rho \sin(\phi)
      (son todos positivos en la región del cuerpo)

      donde el cambio fué
      x = x
      y = \rho \cos(\phi)
      z = \rho \sin(\phi)

      siguiendo,
      M = 2k \left[ -\cos(\phi) \right]_{\pi/4}^{3\pi/4} \int_0^1 \rho^2 \sqrt{1-\rho^2} d\rho

      la integral que queda para \rho es horrible, pero se puede resolver por tabla.
      De todas formas podemos verificar que da el resultado correcto usando el wolfram alpha.

  4. Aoshido dijo:

    Tal cual!!!
    Asi me quedo la ecuacion de la masa, hice exactamente los mismos cambios, sin embargo tuve dos errores yo, abrí mal el modulo mi (region me habia quedado -y<=z<=y) por eso es que yo hablaba de "2 gajos" y ahi te queda definido el "gajo de la derecha" entre 7/4π<=φ<=9/4π.
    Es verdad que la integral que quedaba era horrible yo la resolvi por tabla y aun así fue bastante molesto.
    Respecto a las propiedades de simetría, no suelo aplicarlas pero la verdad que en este caso me significaba hacer 2 veces el gajo (por que yo tenia mal definida el area) por eso es que esta vez, lo utilice jajaja.
    Espero ahora me haya podido expresar mejor jajaja muchisimas gracias por la ayuda!

  5. agustin dijo:

    damian una pregunta como sacas el que angulo vaya de pi/4 a 3/4pi no me quedo claro eso

    • dami dijo:

      Hola agustín,
      Pensá en la proyección sobre el plano yz. El ángulo \phi parte del eje y, y se mueve hacia el eje z.
      Por la condición z \geq |y| sale que \pi/4 \leq \phi \leq 3/4 \pi, pensá en el gráfico del módulo. Se puede hacer más algebráicamente pero es más trabajo.
      Saludos,
      Damián.

  6. Juani dijo:

    Hola Damian, no entiendo porque haces esa rotacion en la coordenadas esfericas, me podrias explicar. Gracias

    • dami dijo:

      Hola Juani,
      Porque de esa forma es más sencillo. Si mirás las superficies coordenadas de esféricas te das cuenta que de esa forma quedan los intervalos de integración entre constantes, mientras que si no lo rotás es mucho más difícil (¿lo intentastes hacer sin rotar?)
      Saludos,
      Damián.