La aceleración de la gravedad en el interior y en el exterior de una distribución esférica y uniforme de masa

La aceleración de la gravedad en un punto P situado a una distancia r de la masa puntual M se define como la fuerza sobre la unidad de masa, situada en dicho punto. Es un vector de módulo

g=G M r 2

dirección radial y sentido hacia la masa puntual M, tal como se muestra en al figura.

El punto P está en el exterior de la esfera de radio R.

Supongamos que el punto P está situado a una distancia r>R a lo largo del eje Z del centro de un planeta de masa M. Dividimos la esfera en discos (de color azul claro) de radio variable y de espesor dz, tal como se muestra en la figura. Para calcular la fuerza que ejerce uno de estos discos sobre la unidad de masa situada en P, dividimos cada disco en anillos (en color amarillo) de radio x, de anchura dx y espesor dz.

Por simetría, como vemos en la figura, las componentes horizontales (a lo largo del eje X e Y) de la fuerza que ejercen los elementos de masa del anillo (considerados como masas puntuales, en color rojo) se anulan de dos en dos, quedando solamente la componente Z de la fuerza resultante que ejerce la masa contenida en el anillo.

Si ρ es la densidad constante de la distribución esférica y uniforme de masa. La masa contenida en el anillo es ρ·(2πx·dx)·dz.  La fuerza resultante que ejerce la masa contenida en el anillo sobre la sobre la unidad de masa situada en P es

G ρ·2πx·dx·dz x 2 + (rz) 2 ·cosα=G ρ·2πx·dx·dz x 2 + (rz) 2 rz x 2 + (rz) 2

Integramos respecto de la variable x, entre los límites 0 e y para calcular resultante de las fuerzas en P debidas a la distribución de masa contenida por el disco de radio y y de espesor dz.

0 y Gρ·2π(rz)dz xdx ( x 2 + (rz) 2 ) 3 = Gρ·2π(rz)dz( 1 rz 1 y 2 + (rz) 2 )

Relacionamos la variable y y z para integrar respecto de la variable z entre los límites –R y +R

z2+y2=R2

La aceleración de la gravedad en el punto P, se obtendrá integrando

g= R +R Gρ·2π( 1 r R 2 + r 2 2rz + z R 2 + r 2 2rz )dz

Integrando por partes el tercer término entre paréntesis

g=Gρ·2π { z+( 1 z r ) R 2 + r 2 2rz 1 3 r 2 ( R 2 + r 2 2rz) 3 } R +R = Gρ·2π{ ( R+ (rR) 2 r (rR) 3 3 r 2 )( R+ (r+R) 2 r (r+R) 3 3 r 2 ) }=Gρ·2π 2 R 3 3 r 2 = G M 4 3 π R 3 ·2π 2 R 3 3 r 2 =G M r 2

En vez de dividir la esfera en discos de radio variable y y espesor dz, dividimos la esfera en capas esféricas concéntricas de radio x y de espesor dx.

Dividimos dicha capa en anillos. Por simetría, como vemos en la figura, las componentes horizontales (a lo largo del eje X e Y) de la fuerza que ejercen los elementos de masa del anillo (considerados como masas puntuales, en color rojo) se anulan de dos en dos, quedando solamente la componente Z de la fuerza resultante que ejerce la masa contenida en el anillo.

Si ρ es la densidad constante de la distribución esférica y uniforme de masa. La masa contenida en el anillo es ρ·(2πxsinθ)·(x·dθ)·dx.  La fuerza resultante que ejerce la masa contenida en el anillo sobre la sobre la unidad de masa situada en P

G ρ·2πx·sinθ·x·dθ·dx (x·sinθ) 2 + (rxcosθ) 2 ·cosα= G ρ·2π x 2 ·sinθ·dx·dθ (x·sinθ) 2 + (rxcosθ) 2 rxcosθ (x·sinθ) 2 + (rxcosθ) 2

Integramos respecto de la variable θ, entre los límites 0 y π para calcular la fuerza sobre la unidad de masa situada en P ejercida por la distribución de masa contenida en la capa esférica de radio x y d espesor dx.

Gρ2π x 2 dx 0 π sinθ(rxcosθ)dθ ( r 2 + x 2 2rxcosθ) 3 = Gρ2π x 2 dx( r 0 π sinθ·dθ ( r 2 + x 2 2rxcosθ) 3 x 0 π sinθcosθ·dθ ( r 2 + x 2 2rxcosθ) 3 )

La primera integral es inmediata, integramos por partes la segunda. El resultado es

Gρ2π x 2 dx ( 1 x r 2 + x 2 2rxcosθ + cosθ r r 2 + x 2 2rxcosθ + r 2 + x 2 2rxcosθ r 2 x ) 0 π = Gρ2π x 2 dx{ ( 1 x(r+x) 1 r(r+x) + (r+x) r 2 x )( 1 x(rx) + 1 r(rx) + (rx) r 2 x ) }= Gρ4π x 2 dx r 2

Gρ4π x 2 r 2 dx

Esta expresión nos da la aceleración de la gravedad producida por una capa esférica de radio x y de espesor dx en un punto P situado a una distancia r>x del centro de la capa esférica.

Calculamos la fuerza ejercida por la esfera de masa M sobre la unidad de masa situada en P, sumando la fuerza que ejerce cada una de las capas esféricas en la que hemos dividido la esfera. El módulo de la aceleración de la gravedad vale.

g= 0 R Gρ4π x 2 r 2 dx = Gρ4π r 2 R 3 3 = G4π r 2 M 4 3 π R 3 R 3 3 =G M r 2

Si tenemos en cuenta que el vector g apunta hacia el centro de la Tierra

g=G M r 2 r>R

El punto P está en el interior de la esfera de radio R.

El punto P está a una distancia r<R del centro de la esfera. Dividimos la esfera en dos, una esfera hueca de radio interior r y radio exterior R. y una esfera maciza de radio r.

Fuerza que ejerce la esfera interior

La fuerza que ejerce la distribución de masa contenida en la esfera maciza de radio r sobre la unidad de masa situada en el punto P a una distancia r del centro, la hemos calculado en el apartado anterior, solamente hemos de cambiar el límite superior de la integral R por r.

g= 0 r Gρ4π x 2 r 2 dx = Gρ4π r 2 r 3 3 =G4π M 4 3 π R 3 r 3 =GM r R 3

Fuerza que ejerce la esfera exterior

Volvemos a calcular la fuerza que ejerce una capa esférica de radio r<x<R y de espesor dx, sobre la unidad de masa situada en P. Ahora bien, la fuerza en P debida a la masa contenida en los anillos que está por encima de P tiene la dirección del eje Z y sentido positivo, mientras que la fuerza en P debida a los anillos que están por debajo de P tiene la misma dirección pero sentido contrario.

Porción de la capa esférica por encima de P

Si ρ es la densidad constante de la distribución esférica y uniforme de masa. La masa contenida en el anillo es ρ·(2πxsinθ)·(x·dθ)·dx.  La fuerza resultante que ejerce la masa contenida en el anillo sobre la sobre la unidad de masa situada en P

G ρ·2πx·sinθ·x·dθ·dx (x·sinθ) 2 + (xcosθr) 2 ·cosα =G ρ·2π x 2 ·sinθ·dx·dθ (x·sinθ) 2 + (xcosθr) 2 xcosθr (x·sinθ) 2 + (xcosθr) 2

Integramos respecto de la variable θ, entre los límites 0 y θp=arccos(r/x) para calcular la fuerza sobre la unidad de masa situada en P ejercida por la porción de capa esférica de radio x y de espesor dx que está por encima de P.

Gρ2π x 2 dx 0 θ p sinθ(xcosθr)dθ ( r 2 + x 2 2rxcosθ) 3 = Gρ2π x 2 dx( r 0 θ p sinθ·dθ ( r 2 + x 2 2rxcosθ) 3 +x 0 θ p cosθ·sinθ·dθ ( r 2 + x 2 2rxcosθ) 3 )

La segunda integral es inmediata, integramos la primera por partes. El resultado es

Gρ2π x 2 dx ( 1 x r 2 + x 2 2rxcosθ cosθ r r 2 + x 2 2rxcosθ r 2 + x 2 2rxcosθ r 2 x ) 0 θ p = Gρ2π x 2 dx{ ( 1 x x 2 r 2 1 x x 2 r 2 x 2 r 2 r 2 x )( 1 x(xr) 1 r(xr) (xr) r 2 x ) }= Gρ2π x 2 dx( x 2 r 2 r 2 x + 1 r 2 )

Porción de la capa esférica por debajo de P

Si ρ es la densidad constante de la distribución esférica y uniforme de masa. La masa contenida en el anillo es ρ·(2πxsinθ)·(x·dθ)·dx.  La fuerza resultante que ejerce la masa contenida en el anillo sobre la sobre la unidad de masa situada en P

G ρ·2πx·sinθ·x·dθ·dx (x·sinθ) 2 + (rxcosθ) 2 ·cosα= G ρ·2π x 2 ·sinθ·dx·dθ (x·sinθ) 2 + (rxcosθ) 2 rxcosθ (x·sinθ) 2 + (rxcosθ) 2

Integramos respecto de la variable θ, entre los límites θp=arccos(r/x) y π para calcular la fuerza sobre la unidad de masa situada en P ejercida por la distribución de masa contenida en la capa esférica de radio x y d espesor dx.

Gρ2π x 2 dx θ p π sinθ(rxcosθ)dθ ( r 2 + x 2 2rxcosθ) 3 = Gρ2π x 2 dx( r θ p π sinθ·dθ ( r 2 + x 2 2rxcosθ) 3 x θ p π sinθcosθ·dθ ( r 2 + x 2 2rxcosθ) 3 )

La primera integral es inmediata, integramos la segunda por partes. El resultado es

Gρ2π x 2 dx ( 1 x r 2 + x 2 2rxcosθ + cosθ r r 2 + x 2 2rxcosθ + r 2 + x 2 2rxcosθ r 2 x ) θ p π = Gρ2π x 2 dx{ ( 1 x(r+x) 1 r(r+x) + (r+x) r 2 x )( 1 x x 2 r 2 + 1 x x 2 r 2 + x 2 r 2 r 2 x ) }= Gρ2π x 2 dx( 1 r 2 x 2 r 2 r 2 x )

Como la fuerza ejercida en P por la distribución de masa contenida en la capa semiesférica que está por encima de P y la fuerza ejercida en P por la distribución de masa contenida en por la parte de la capa semiesférica que está por debajo de P tienen el mismo valor pero signos contrarios. La fuerza neta en P debida a la capa semiesférica completa de radio x y de espesor dx es cero.

Fuerza total

Por tanto, la fuerza ejercida sobre la unidad de masa situada en P por la distribución de masa contenida en la esfera hueca de radio interior r y exterior R es cero.

La fuerza ejercida sobre la unidad de masa situada en P a una distancia r del centro de la esfera de radio R, solamente es debida a la distribución de masa contenida en la esfera maciza de radio r<R

Si tenemos en cuenta que el vector g apunta hacia el centro de la Tierra

g=GM r R 3                r<R

La aceleración de la gravedad g aumenta linealmente de 0 a GM/R2, en el interior de una distribución esférica y uniforme de masa M y radio R. En el exterior de dicha distribución, la aceleración de la gravedad disminuye en proporción inversa al cuadrado de la distancia al centro de dicha distribución GM/r2

En la figura, se muestra la gráfica del módulo de la aceleración de la gravedad g en función del cociente r/R para el planeta Tierra, G=6.67·10-11 Nm2/kg2, R=6.37·106 m, M=5.98·1024 kg. La aceleración de la gravedad en la superficie de la Tierra vale g=9.83 m/s2

Energía potencial

Una partícula de masa m experimenta una fuerza F=mg. La energía potencial Ep(r) correspondiente a esta fuerza conservativa se calcula del siguiente modo

F= d E p dr r d E p = r F(r)·dr E p (r)= r F(r)·dr

Tomamos como nivel cero de energía potencial Ep(∞)=0

En la figura, se representa la energía potencial Ep(r) en función del cociente (r/R)

El principio de conservación de la energía

1 2 m v 2 + E p (r)=cte

Conocida la posición r0 inicial y la velocidad inicial v0 de la partícula, calculamos la velocidad v de la partícula para una determinada posición r.

Ejemplo: Se excava un túnel a través de un diámetro de la Tierra, y se suelta una partícula de masa m en uno de sus extremos. Determinamos la velocidad de la partícula cuando se encuentra a una distancia r del centro de la Tierra

Se iguala la energía en la posición inicial r=R, v=0, con la energía en la posición considerada r.

1 2 m v 2 3 2 GMm R + GMm R 3 r 2 2 = GMm R v 2 = GM R GM R 3 r 2

Cuando pasa por el centro de la Tierra r=0, v=7913.0 m/s.

Los datos son: masa de la Tierra M=5.98·1024 kg, radio de la Tierra R=6.37·106 m, G=6.67·10-11 Nm2/kg2

Presión en el centro de la Tierra

Consideremos un elemento de volumen (en color rojo) de área dA y de espesor dr. Las fuerzas sobre el elemento de masa dm=ρ·dA·dr son

En el equilibrio

p(r)·dA-p(r+dr)·dA-g(r)·dm=0

La ecuación fundamental de la hidrostática es

dp=-ρg(r)dr

La presión a una distancia r del centro de la Tierra es

p 0 p= r R M 4 3 π R 3 GM r R 3 dr = 3G M 2 8π R 6 ( R 2 r 2 )

donde p0=1.013·105 Pa es la presión para r=R es decir, la presión atmosférica

En el centro de la Tierra r=0.

p= p 0` + 3G M 2 8π R 4 =1.73· 10 11 Pa

Con G=6.67·10-11 Nm2/kg2, R=6.37·106 m, M=5.98·1024 kg

La presión atmosférica es despreciable frente a la presión en el centro de una distribución esférica y uniforme de masa.