If the Flat Earth is infinitely large and has a thickness b and a homogeneous density ρ, all horizontal gravitational force components would cancel each other at every place on earth and only leave the sum of all vertical forces components. There would be no center of mass like on the globe.
But is the sum of the vertical force components of all volume elements of the infinitely large plane not infinite, assuming that Newton's law of universal gravitation is applied like on the Globe Earth?
Newton's law of universal gravitation
It turns out: NO
If the earth were an infinite plane of a certain thickness b and homogeneous density ρ, you would get a homogeneous gravity field that acts perpendicular to the surface everywhere.
Such a Flat Earth would not collapse into a sphere, theoretically. But this plane would not be stable. The slightest irregularities in the density would disturb the equilibrium and crunch the earth, except if it is infinitely rigid. An infinite plane would also mean an infinite mass for the earth.
The formula for calculating g for an infinite plane turns out to be really simple. The calculations are shown below:
To achieve g = 9.806 m/s2 gravitational acceleration, you can calculate how thick the Flat Earth has to be:
Note: In the equation (5) for the gravitational acceleration g for an infinite plane there is no term for the altitude h like on the Globe Earth! This implies that the gravitaional acceleration is the same at every altitude.
This is due to the homogeneous gravity field that such an infinte plane would create.
So what do we observe? Is g constant at each altitude h from the surface of the earth or is it changing with 1/(R+h)2?
We know the gravitational acceleration g is changing with 1/(R+h)2 from measuring it with a scale and a weight at different altitudes or by using a gravimeter. And it also depends on latitude because the earth rotates and is not a perfect sphere.
So Gravity measurements falsify the hypothesis of an infinite Flat Earth plane. Observations of gravity on earth confirm a rotating ellipsoid.
For the following calculations we divide the infinite Flat Earth plane into an infinite number of infinitely small volume elements of density ρ and calculate the gravitational acceleration from the sum of all volume elements for any location at an altitude of h from the surface.
I choose a cylindrical coordinate system with the origin at the observer, the horizontal radial distance x from a volume element pointing away from the observer at an azimuthal angle φ and the vertical down is the Z axis.
The distance r of any volume element from the observer is then:
The mass of any volume element in the cylindrical coordinate system is given by:
The magnitude of the gravitational acceleration from a certain volume element i in any direction φi is according to Newton:
Because of the symmetry around the Z axis all horizontal components gx of the gravitational acceleration cancel. We have only to calculate the contributions of all vertical accelerations gz of each volume element i:
The blue term is derived from the geometry. If we have a right angle triangle with height z and hypothenuse r, the z component of the acceleration is gz = g · z/r.
From here on I only refer to the vertical z component of the accelerations, so I will not write the index z explicitly anymore.
Lets insert the equation for
To calculate the sum of all vertical accelerations, we have to integrate the accelerations of all volume elements. This gives 3 nested integrals: one for the angle from φ = 0 to 2π, one for the horizontal distance from x = 0 to infinity and one for the vertical distance from the observer from z = h to h + b:
Lets solve the integrals one by one, starting with the inner most for the angle φ from 0 to 2π. Because the terms x and z do not depend on φ we can take them out of the integral, which simply results in 2π:
Inserting this result into (11) gives:
Now lets solve the next inner integral for x from 0 to infinity:
This is an amazing simple result. Inserting it into (13) leaves us with the following simple integral:
Now lets solve the last integral for z = h to h + b:
Inserting this result into (15) gives our final equation for the gravitational acceleration on an infinite Flat Earth plane: