Gravity Between Planets and Planetary Elliptical Orbit Simulation

1. Introduction

In 1609, Kepler[1] published the first and second laws of planetary motion. The first law points out that the orbit of a planet is an ellipse, and the sun is located at a focal point of the ellipse; the second law points out that the connection between the planet and the sun swept through equal areas within equal time intervals. The essence of Kepler's second law is the law of conservation of angular momentum for planetary operations [2]. In 1687, Isaac Newton expounded the law of universal gravitation in "Mathematical Principles of Natural Philosophy" [3]. The law of universal gravitation pointed out that there is gravity between any two objects in the universe, and its size is proportional to their mass and inversely proportional to the square of their distance. Newtonian gravity provides the centripetal force of planets, explaining the source of power for planets to operate in elliptical orbits.

This paper analyzes the gravity between planets and successfully simulates the elliptical orbit of planets.

2. The Equivalent Radius of the Planet's Graviton Emission

"Graining, gravitational field and graviton - inference about the frequency of gravitational energy waves" [8] believes that the most basic unit of matter is nucleons (collective name of protons and neutrons). All nucleons emit gravitons. The energy carried by gravitons is the Planck constant h, and the value is 6.626×10^ (-34) J·s. Gravitators propagate in space with gravitational energy waves. The gravitational energy wave resonates with other nucleons, transmitting energy to form gravity. For planets [9], gravitons emitted by nucleons inside the planet interact with other nucleons inside the planet to form the cohesion of the planet. The gravitons emitted by nucleons near the outside of the planet are partially emitted outside the sphere and propagate in space with gravitational energy waves. The gravitational energy waves encounter nucleons from other planets and resonate with them to form gravitational force between the planets.

Figure 1 is a schematic diagram of the equivalent spherical surface of a planet's graviton emitted. In thefigure, rs is the planet's radius, and rso is the thickness of the graviton shell sent to the outside of the sphere. There should be a spherical layer rse in the middle. It can be considered that all gravitons on the planet are emitted by this spherical layer. If this spherical layer is used as the equivalent spherical layer emitted by the planet's gravitons, for general circumstances, it can be considered that this spherical layer is in the middle of the planet's graviton emission shell."On the nuclear force is the manifestation of gravity at the microscopic distance" [10] The article calculates the number of gravitons emitted by a single nucleon based on the binding energy of hydrogen. On the relationship between atomic structure and basic force" [11] The article, based on the analysis that the resonance of gravitational energy waves and nucleons conforms to the normal distribution, the number of gravitons emitted per second is:

$$n_{ng}=2.227\times {10}^{22}$$

(1)

The ratio of gravitons passing through nucleons can be absorbed by nucleons is:

$$k_{ng}=0.378$$

(2)

"The attempt to correct the universal gravitational formula from the proportion of the planet emitted to the outside of the ball - the ratio of the outside of the ball graviton in the deflection gravitational theory" [12] article calculates the number of gravitons sent to the outside of the ball. The thickness of the outer layer where the planet can emit gravitons outside the ball is:

$$r_{so}={\displaystyle \frac{6m_0}{k_{ng}{r}_0^2{\rho }_s}}={\displaystyle \frac{k_{s\rho }}{{\rho }_s}}$$

(3)

In the formula, kng is the ratio in which gravitons passing through the nucleus can be absorbed by the nucleus, m0 is the mass of the nucleus, r0 is the radius of the nucleus, and ρs is the density of the shell matter of the planet. The value after the data is corrected is:

$$k_{s\rho }={\displaystyle \frac{6m_0}{k_{ng}{r}_0^2}}={\displaystyle \frac{6\times 1.6749\times {10}^{-27}}{0.378\times {\left(0.8\times {10}^{-15}\right)}^2}}=41540$$

(4)

The number of gravitons sent by the planet to the outside of the ball is:

$$n_{go}={\displaystyle \frac{32{\pi }^2{n}_{n\mathrm{g}}{k}_o}{k_{ng}{r}_0^2}}\left[{r_s}^2-{\displaystyle \frac{4m_0}{k_{ng}{r}_0^2}}{\displaystyle \frac{r_s}{{\rho }_s}}+6{\left({\displaystyle \frac{m_0}{k_{ng}{r}_0^2}}\right)}^2{\displaystyle \frac{1}{{\rho }_s^2}}\right]$$

(5)

In the formula, ngo is the number of gravitons sent to the outside of the sphere by the planet, nng is the number of gravitons emitted by a single nucleus in one second, ko is the proportion of the outer part of the nucleus occupying the entire surface area, kng is the proportion of the graviton passing through the nucleus to be absorbed by the nucleus, r0 is the radius of the nucleus, rs is the radius of the planet, m0 is the mass of a single nucleus, ρs is the density of the shell matter of the planet, substituting the constant and calculating:

$$n_{go}=1.599\times {10}^{55}\left[{r_s}^2-2.679\times {10}^4{\displaystyle \frac{r_s}{{\rho }_s}}+2.876\times {10}^8{\displaystyle \frac{1}{{\rho }_s^2}}\right]$$

(6)

When the planet is relatively large, the above formula is approximate:

$$n_{go}\approx {\mathrm{k}}_{\mathrm{gr}}{r_s}^2=1.599\times {10}^{55}{r_s}^2$$

(7)

Table 1 is a statistical table of the number of extraspheric gravitons emitted from the solar system planet. The first column in the table is the planet name, and the second column is the planet mass.

The third column is the planet's radius, the fourth column is the planet's shell density, the sun is the photosphere density, the earth is the crust density, the fifth column is the thickness of the shell that the planet can send out of the ball, the sixth column is the ratio of the shell thickness sent out of the ball to the planet's radius, it can be seen that this ratio is on the order of 10^-6, the seventh column is the number of gravitons sent out of the ball calculated based on the above formula (6), and the eighth column is the number of gravitons sent out of the ball calculated based on the above formula (7), and it can be seen that the simplified formula has no effect on the calculation of the number of gravitons sent out of the ball, and the last column is several commonly used constants. From the above analysis, it can be seen that the thickness of the shell of the planet emitting gravitons outward is too small relative to the planet's radius and can be ignored. Therefore, the equivalent shell radius of the planet's emitting gravitons is approximately equal to the planet's radius.

$$r_{se}=r_s-{\displaystyle \frac{1}{2}}{r}_{so}\approx r_s$$

(8)

3. The Center of Mass and Gravity of the Planet

Figure 2 is an analysis diagram of the gravitational action between two planets. In the picture, the planet E rotates around the planet S. The center of mass of the central planet S is S0 and the radius is rs.The mass is ms, the center of mass of the planet E is E0, the radius is re, the mass is me, and the center of mass distance between the planets is R0. Now, the effect of planet S on planet E is analyzed. In the figure, only the gravitons emitted by planet S facing the nucleons on the sphere of planet E can form gravity on planet E's shell nucleons facing planet S.This phenomenon can be proved by the gravity double valley phenomenon during the solar eclipse [13]. Before the solar eclipse, objects on the ground are subjected to the dual gravity of the sun and the moon, and the gravity of the earth that the object receives will decrease; during the solar eclipse, the gravity of objects on the ground is measured, and the results show that the gravity of objects on the ground is the same as the gravity when there is only the sun, which means that the gravitons received by the sun on the ground are blocked by the moon; when the solar eclipse ends, the sun is no longer blocked by the moon, and the objects on the ground are superimposed by the gravity of the sun and the moon, and the gravity decreases again, which forms the gravity double valley phenomenon during the solar eclipse.

Accurately calculate the effect of gravitons emitted by planet S on planet E, and calculate the surface element on planet S on planet E.

When the distance between the planets is much greater than the planet's radius, it can be approximately considered that the energy transmitted between the two planets' spheres is equal, so that the energy transmitted from the planet S sphere to the planet E sphere is proportional to the area of the sphere.

In the picture:

$$\sin {\alpha }_e={\displaystyle \frac{r_s-r_e}{R_0}}$$

(9)

$${\alpha }_s={\displaystyle \frac{\pi }{2}}-{\alpha }_e$$

(10)

$$\cos {\alpha }_s=\cos \left({\displaystyle \frac{\pi }{2}}-{\alpha }_e\right)=\sin {\alpha }_e={\displaystyle \frac{r_s-r_e}{R_0}}$$

(11)

The surface area calculation formula of the spherical crown is S=2πRh. Obviously, the equivalent center of the nuclear nucleus that the planet S ball works against the planet E is on the connection line between the center of mass of the two planets. As shown in the figure, let the distance between the center of gravity action of the center of mass of the planet S S0 is hs, Ss1 is the area of the spherical crown above the gravitational action point, Ss2 is the area of the spherical crown that emits the graviton.

$$S_{S1}=2\pi r_s\left(r_s-h_s\right)$$

(12)

$$S_{S2}=2\pi r_s\left(r_s-r_s\cos {\alpha }_s\right)$$

(13)

$$S_{S2}=2S_{S1}$$

(14)

$$2\pi r_s\left(r_s-r_s\cos {\alpha }_s\right)=2\times 2\pi r_s\left(r_s-h_s\right)$$

(15)

$$r_s-r_s\cos {\alpha }_s=2r_s-2h_s$$

(16)

$$h_s={\displaystyle \frac{1+\cos {\alpha }_s}{2}}{r}_s={\displaystyle \frac{1+\frac{r_s-r_e}{R_0}}{2}}{r}_s={\displaystyle \frac{R_0+r_s-r_e}{2R_0}}{r}_s\approx {\displaystyle \frac{1}{2}}{r}_s$$

(17)

For planet E, the area affected by planet S is greater than half the sphere. Therefore, the equivalent center of gravity is closer to the center of mass of planet E. As shown in the figure, let the distance between the center of gravity and the center of mass is he, the area of ​​the spherical crown below the gravitational action point of planet E is Se1, and the entire spherical crown area that can receive gravitons is Se2:

$$S_{e1}=2\pi r_e\left(r_e-h_e\right)$$

(18)

$$S_{e2}=2\pi r_e\left[r_e+r_e\cos \left({\displaystyle \frac{\pi }{2}}-{\alpha }_e\right)\right]=2\pi r_e\left(r_e+r_e\sin {\alpha }_e\right)$$

(19)

$$S_{e2}=2S_{e1}$$

(20)

$$2\pi r_e\left(r_e+r_e\sin {\alpha }_e\right)=2\times 2\pi r_e\left(r_e-h_e\right)$$

(21)

$$r_e+r_e\sin {\alpha }_e=2r_e-2h_e$$

(22)

$$h_e={\displaystyle \frac{1-\sin {\alpha }_e}{2}}{r}_e={\displaystyle \frac{R_0-r_s+r_e}{2R_0}}{r}_e\approx {\displaystyle \frac{1}{2}}{r}_e$$

(23)

Suppose the distance between the equivalent gravity center of the two planets is Rg, and the distance between the gravity action is

$$R_g=R_0-h_e-h_s\approx R_0$$

(24)

The energy transmitted by planet S to planet E is the ratio of the area of the sphere occupied by planet E from planet S distance R0 to the total number of gravitons emitted by planet S:

$$F_{se}={\displaystyle \frac{\pi r_e^2}{4\pi R_0^2}}{n}_{sg}h={\displaystyle \frac{r_e^2}{4R_0^2}}{k}_{gr}{r}_s^{2}h={\displaystyle \frac{k_{gr}h{r}_e^2{r}_s^2}{4R_0^2}}=G_i{\displaystyle \frac{r_e^2{r}_s^2}{R_0^2}}$$

(25)

$$G_i={\displaystyle \frac{1}{4}}{k}_{gr}h$$

(26)

Here Gi is the gravitational coefficient, kgr is the proportional coefficient of the graviton sent to the outside of the ball by the planet, and h is the Planck constant. From the above analysis, we can find that the energy (gravity) transmitted by nucleons between planets is proportional to the planet's surface area, rather than to the planet's mass.

4. Analysis of the Effect of Gravity on Planetary Movement

Figure 3 is an analysis diagram of the planet's revolution and rotation caused by planet gravity. In the figure, S0 is the center of mass of the central planet S.E is a planet orbiting the central planet S. When the planet E is at E0, the gravitational action point of planet S on planet E is Eg0. Since the gravitational action point is separated from the planet's center of mass, gravity is divided into Fs0 related to the center planet S and Fe0 related to the planet E. When the same planet E is at position E1, gravity is divided into Fs1 related to the center planet S and Fe1 related to the planet E. When the planet E is at position E2, gravity is divided into Fs2 related to the center planet S and Fe2 related to the planet E.It can be seen that Fs0, Fs1, and Fs2 related to the center planet S act on the gravitational action point E, forming the driving force for the planet E to orbit the planet S. From the effect, Fs forms the centripetal force of planet E orbiting the planet S.

Since the gravitational action point Eg is separated from the planet's center of mass, another effect of gravity received by the gravitational action point Eg is to drag the planet E to run. Generally, external forces rotate around the center of the sphere. Without considering the planet E orbiting the central planet S, it can be considered that the planet E has always been affected by a gravitational Fe, and this point of action is fixed at the gravitational action point Eg. Therefore, this is not an external force rotates around the planet, but the center of mass of the planet rotates around the gravitational action point Eg. At the initial position E0, the planet's center of mass tends to run towards A0. With the increase of time, the angle of rotation of the planet's center of mass increases. Under the combined action of the two components of gravity, when the planet is at position E1, the center of mass moves to point A1, and when the planet is at position E2, the center of mass moves to point A2.

According to the above analysis, the total gravity Fse of the planet is:

$$F_{se}=F_s+F_e$$

(27)

The components related to the planet E revolution are:

$$F_s=m_e{\displaystyle \frac{v_s^2}{R_0}}=m_e{\mathrm{\Omega }}_s^2{R}_0$$

(28)

In the formula, gravity causes the linear velocity of the planet E to orbit the central planet S to be vs, unit m/s, the angular velocity of the revolution to be Ωs, unit radian/s, and me is the mass of the object (unit kg).

For a rotating body, the moment of inertia is a physical quantity that describes the magnitude of inertia when an object rotates about a certain axis. For a uniform sphere, its moment of inertia is:

$$I_c={\displaystyle \frac{2}{5}}m{r}^2$$

(29)

Where m is the mass of the sphere and r is the radius of the sphere.

For rotating bodies whose rotation axis does not coincide with the centroid, the parallel axis theorem considers the moment of inertia:

$$I=I_c+m{d}^2$$

(30)

Here m is the mass of the rigid body, Ic is the moment of inertia around the rotation axis of the centroid, and d is the distance between the axis of rotation and the centroid.

From this we can see that the moment of inertia of the planet E rotating around the gravitational action point Eg is:

$$I_e={\displaystyle \frac{2}{5}}{m}_e{r}_e^2+m_e{h}_e^2$$

(31)

Where me is the mass of the planet E, re is the radius of the planet E, and he is the distance between the gravitational action point Eg and the center of mass of the planet E.

According to the rotation law of rigid body fixed axis:

$$M_Z=I\alpha$$

(32)

Where Mz represents the external torque for a certain fixed axis, I represents the moment of inertia of the rigid body about a given axis, and α represents the angular acceleration. Here the torque is the component force of gravity and the rotation of the planet. The force arm is the distance he between the gravitational action point Eg and the center of mass of the planet. According to the definition of angular acceleration, there are:

$$\alpha ={\displaystyle \frac{d^2{\theta }_e}{d^{2}t}}={\displaystyle \frac{{\mathrm{\Delta }}^{2}t{\omega }_e^2}{{\mathrm{\Delta }}^{2}t}}={\omega }_e^2$$

(33)

At this time, the law of rotation of the rigid body fixed axis can be written as:

$$F_e{h}_e=I{\omega }_e^2$$

(34)

$$F_e={\displaystyle \frac{I{\omega }_e^2}{h_e}}=\left({\displaystyle \frac{2}{5}}{m}_e{r}_e^2+m_e{h}_e^2\right){\displaystyle \frac{{\omega }_e^2}{h_e}}$$

(35)

For the rotation of a planet, it refers to the component of rotation in the planet's rotation plane, which is generally not on the same plane as the actual rotation of the planet. Therefore, the angular velocity of the planet here is not the planet's rotation angular velocity that is usually observed.

For objects that move in a circular motion, the relationship between linear velocity and angular velocity is as follows:

$$v_0=R_0{\mathrm{\Omega }}_s$$

(36)

$$v_0\mathrm{\Delta }t\sin {\beta }_0={\theta }_s$$

(37)

According to the law of conservation of momentum, momentum cannot be generated and disappeared out of thin air. For planet E, the planet's revolution generates a positive momentum, and planet E rotates around the center of mass Eg to produce a reverse momentum, and these two should be equal:

$$p=m_e{v}_s=m_e{\mathrm{\Omega }}_s{R}_0=m_e{v}_e=m_e{\omega }_e{h}_e$$

(38)

$${\mathrm{\Omega }}_s{R}_0={\omega }_e{h}_e$$

(39)

$${\mathrm{\Omega }}_s={\displaystyle \frac{h_e}{R_0}}{\omega }_e$$

(40)

Bring the above result into formula 27:

$$F_{se}=m_e{\mathrm{\Omega }}_s^2{R}_0+\left({\displaystyle \frac{2}{5}}{m}_e{r}_e^2+m_e{h}_e^2\right){\displaystyle \frac{{\omega }_e^2}{h_e}}$$

(41)

$$F_{se}=m_e{R}_0{\left({\displaystyle \frac{{\omega }_e{h}_e}{R_0}}\right)}^2+\left({\displaystyle \frac{2}{5}}{m}_e{r}_e^2+m_e{h}_e^2\right){\displaystyle \frac{{\omega }_e^2}{h_e}}$$

(42)

$${\omega }_e=\sqrt{{\displaystyle \frac{F_{se}}{m_e\left[\frac{h_e^2{R}_0}{R_0^2}+\frac{1}{h_e}\left(\frac{2}{5}{r}_e^2+h_e^2\right)\right]}}}=\sqrt{{\displaystyle \frac{F_{se}}{m_e\left(\frac{h_e^2}{R_0}+\frac{2r_e^2}{5h_e}+h_e\right)}}}$$

(43)

The ratio of component force used for planet rotation to the entire gravity is:

$$\begin{array}{l}{k}_F={\displaystyle \frac{F_e}{F_{se}}}={\displaystyle \frac{\left(\frac{2}{5}{m}_e{r}_e^2+m_e{h}_e^2\right)\frac{{\omega }_e^2}{h_e}}{m_e{R}_0{\left(\frac{{\omega }_e{h}_e}{R_0}\right)}^2+\left(\frac{2}{5}{m}_e{r}_e^2+m_e{h}_e^2\right)\frac{{\omega }_e^2}{h_e}}}\\ ={\displaystyle \frac{\left(\frac{2}{5}{r}_e^2+h_e^2\right)}{\frac{h_e^3}{R_0}+\left(\frac{2}{5}{r}_e^2+h_e^2\right)}}\approx {\displaystyle \frac{\left(\frac{2}{5}{r}_e^2+\frac{1}{4}{r}_e^2\right)}{\frac{r_e^3}{8R_0}+\left(\frac{2}{5}{r}_e^2+\frac{1}{4}{r}_e^2\right)}}={\displaystyle \frac{26R_0}{5r_e+26R_0}}\approx 1\end{array}$$

(44)

From this we can see that almost all the gravity of the planet is used for the rotation of the planet, and only a little bit is left for the rotation of the planet. Here, it can also be considered that at the point of gravity, the torque of the planet's revolution is equal to the moment of the planet's rotation. Since the force arm of the planet's revolution is much larger than the force arm of the planet's rotation, the force used for the planet's revolution is much smaller than the force of the planet's rotation.

The movement of a planet is not a simple unforced system, it is subject to gravity at any time. Therefore, the angular momentum of a planet at any point includes: the angular momentum of the initial velocity of the planet relative to the angular momentum of the central planet, gravity forms the angular momentum of the planet's revolution, and gravity forms the angular momentum of the planet's rotation reflected to the sphere of the revolution:

$$L_0=m_e{R}_0{v}_0\sin {\beta }_0+m_e{R}_0{v}_s-m_e{R}_0{v}_e$$

(45)

Since vs is equal to ve, after the above formula is included

$$L_0=m_e{R}_0{v}_0\sin {\beta }_0$$

(46)

The angular momentum of the planet's revolution is only the angular momentum caused by the initial velocity. It can be seen that the offset of the angular momentum of the planet's revolution generated by gravity is the fundamental reason for the conservation of angular momentum in the planet's revolution system.

For different positions, according to the law of conservation of angular momentum, there are:

$$L=R_0{m}_e{v}_0\sin {\beta }_0=R_1{m}_e{v}_1\sin {\beta }_1$$

(47)

In the formula, L is the angular momentum, R0, v0, and β0 are the distance between the planet and the central planet at its initial position, the linear velocity of the planet, the angle between the gravity line and the linear velocity, R1, v1, and β1 are the distance between the central planet after the change in the position of the planet, the linear velocity of the planet, the angle between the gravity line and the linear velocity. From the above formula:

$$v_1={\displaystyle \frac{R_0{v}_0\sin {\beta }_0}{R_1\sin {\beta }_1}}$$

(48)

At present, the minimum orbital speed and maximum orbital speed of the moon have not been found. Kepler's area law is used to estimate the minimum speed of the moon's apogee. Assuming that the semi-major axis of the orbit of the planet E is a and the semi-major axis is b, then the elliptical area of the orbit of the planet E is:

$$S=\pi ab$$

(49)

According to Kepler's law of area, when the planet E moves in an elliptical orbit, the area it swepts through with the sun's line within an equal time. If the planet E runs into n parts, the ellipse area is also divided into n parts. Assuming the unit area is Sn, assuming the linear velocity of planet E running is v0, the angle between planet E and the gravitational line is β0, the distance between planet E and the central planet S is R0, and the time interval is Δt:

$$S_n={\displaystyle \frac{S}{n}}={\displaystyle \frac{1}{2}}{R}_0{v}_0{\displaystyle \frac{T}{n}}\sin {\beta }_0$$

(50)

When the moon is at an apogee, the direction of the moon's velocity is at an angle of 90° with the direction of gravity line. The above formula is simplified to:

$$S={\displaystyle \frac{1}{2}}{R}_0{v}_{0}T$$

(51)

The initial linear velocity of the moon's apogee is:

$$v_0={\displaystyle \frac{2S}{R_{0}T}}$$

(52)

Figure 4 is a planet operation analysis diagram, where the planet E operation cycle is divided into n parts by time, and the unit time is Δt.In the figure, the center S0 of the central planet S is the coordinate origin, R0 is the distance between the initial position of the planet and S0, v0 is the initial running speed of the planet E, β0 is the angle between the initial running direction of the planet and the gravitational line, θe is the angle at which the center of mass of the planet E rotates around the gravitational action point Eg in unit time Δt, and θs is the angle at which the planet E rotates around the central planet S under the action of the gravitational component, and moves v0 to point A0, and its end point is A1, and then rotates A1 to the angle of θs to position A2, A2 is the end point of the center of mass of the planet E through Δt time. here:

$$B_0{A}_0=h_e\mathit{sin}{\theta }_e$$

(53)

$$B_0{E}_g=h_e\mathit{cos}{\theta }_e$$

(54)

The coordinates of point A1 are:

$$x_{a1}=v_0\mathrm{\Delta }t\mathit{sin}{\beta }_0-h_e\mathit{sin}{\theta }_e$$

(55)

$$y_{a1}=R_0-v_0\mathrm{\Delta }t\mathit{cos}{\beta }_0-h_e+h_e\mathit{cos}{\theta }_e$$

(56)

The distance from point A1 to S0 is:

$$R_1=\sqrt{x_{a1}^2+y_{a1}^2}=\sqrt{{\left(v_0\mathrm{\Delta }t\mathit{sin}{\beta }_0-h_e\mathit{sin}{\theta }_e\right)}^2+{\left(R_0-v_0\mathrm{\Delta }t\mathit{cos}{\beta }_0-h_e+h_e\mathit{cos}{\theta }_e\right)}^2}$$

(57)

The angle α between point A1 and the initial position of planet E satisfies:

$$\mathit{tan}\alpha ={\displaystyle \frac{v_0\mathrm{\Delta }t\mathit{sin}{\beta }_0-h_e\mathit{sin}{\theta }_e}{R_0-v_0\mathrm{\Delta }t\mathit{cos}{\beta }_0-h_e+h_e\mathit{cos}{\theta }_e}}$$

(58)

When the planet E rotates through Δt time to the A2 position under the action of the gravitational component Fs, the distance between planet E and S0 remains unchanged to R1, and the angle between the position A2 of planet E and the initial position increases to:

$${\theta }_c={\theta }_s+\alpha ={\mathrm{\Omega }}_s\mathrm{\Delta }t+\alpha$$

(59)

It can be seen that the planet's rotation is inertia when it is running on the equilibrium planet.

5. Simulation of the Elliptical Orbit of the Planet

The relevant parameters of the moon and the earth [14] are as follows: the average radius of the moon is about 1737.10km, the mass is 7.342×10^22kg, the average radius of the earth is 6371.393km, the perigee distance of the moon: 363300km; the apogee distance is 405696km; the average revolution period is 27.32 days; the average revolution speed is 1.023 kilometers/sec; the rotation period is: 27 days, 7 hours, 43 minutes, 11.559 seconds (27.32 days, synchronous rotation); the inclination angle of the rotation axis varies between 3.60° and 6.69°, the semi-major axis of the moon orbit orbit is 384403km, and the criterion rate is 0.0549.

Based on the above derivation, the moon's orbit can be simulated. Table 2 is a partial screenshot of the simulation data table for the moon's orbit around the earth.The most column in the table is the correlation constant, re is the radius of the moon, me is the mass of the moon, rs is the radius of the earth, and Gr is the gravitational coefficient for the application of the planet's semi-compassing calculation. Unlike the commonly used gravitational coefficient G that uses the mass of the planet to calculate gravity, T is the orbital period, which refers to the time when the moon orbits the earth, unit seconds, n is the number of equal parts of the period. Here is 10,000, and Table 1 above is just a few of the data. Δt is unit time, S is the orbital area calculated based on the moon's semi-major axis and eccentricity, and is used to calculate the initial velocity of the moon's apogee. The first column R in Table 2 shows the distance between the moon and the earth when the moon is at different positions. The initial value is the apogee. The second column he is the distance between the gravitational action point Eg, which acts on the moon and the center of mass. The third column is the angle between the direction of the moon and the gravitational line. The fourth column is the velocity of the moon. The calculation of the initial velocity of the apogee uses Kepler's area law. The fifth column is the gravity of the earth to the moon. Here the earth and the radius of the moon are used to calculate the gravity of the earth to the moon. The sixth column ωe is the angular velocity generated by gravity to cause the moon to be driven by the autobiography of the moon. The seventh column Ωs is the angular velocity generated by the force used for the moon's revolution. The eighth column α is the angle in which the moon's initial inertia v0 and gravity causes the moon to be rotated. The tenth column θc is the angular displacement actually generated by the moon unit time. It is the sum of θs and α. Column 11 θ is the accumulation of angular displacement per unit time, and columns 12-13 are the rectangular coordinates used when plotting simulated data.

Figure 5 is a simulation diagram of the moon orbit directly generated in Table 2. The apogee is 4.056960E+08 m, given by the initial simulation value, and the perigee is 3.5547541E+8m.

Returns the ariel point 4.051655E+08m, where E is the commonly used power base in the table, 4.056960E+08m is the same as 4.056960×10^8m, and it can be seen from the parameters that it is an ellipse. The angle between the moon's running direction and the gravitational line is at a given value of 1.57080 radians at the aurora point. As the moon moves from avera to perigee, this angle gradually decreases. After decreasing to the minimum value of 1.5056601 radians, it begins to gradually increase. After passing through 1.57080, it continues to increase. When it increases to the maximum value of 1.6359258, it begins to gradually decrease, and finally returns to the aurora, with the return value of 1.5711103 radians. The velocity value of the moon's apogee is 968.061m/s, the perigee velocity is 1107.073m/s, and the return value of the apogee is 969.329m/s. Under the action of the earth's gravity, the minimum rotation angular velocity of the moon on the white path is 7.355E-5 radian/s, the maximum angular displacement is 8.409E-5 radian/s, the average revolution angular velocity is 2.662E-6 radian/s, and the simulated average is 2.725E-6 radian/s. The angular velocity of the moon generated by the earth's gravity is 29.49 times greater than the revolution angular velocity. More than 99.9% of the earth's gravity is used for the rotation of the moon.

The Earth's orbit around the sun can also be simulated. The average radius of the earth [15] is 6.371×10^6m, the earth's mass is 5.972×10^24kg, the sun's radius is 6.955×10^8m, the earth's orbital period is 365.256363 days (stellar day), 23:56:4.100 seconds (stellar hour) (stellar day), the aphrodisiac distance is 1.52097597×10^11m, the perihelion distance is 1.4709845×10^11m, the average revolution speed is 29783 m/s, the maximum revolution speed is 30287 m/s, the minimum revolution speed is 29291 m/s, the semi-major axis of the orbit 1.49598×10^11m, and the semi-short axis of the orbit 1.4958×10^11m.

Figure 6 is the Earth orbit simulation diagram directly generated by the Earth orbit data simulation table. The aphrodisiac point is 1.520976E+11m, given by the simulation initial value. Point 1.473349E+11m returns to the aphrodisiac point 1.516298E+11m. It can be seen from the parameters that it is An ellipse. The angle between the earth's direction of movement and the gravitational line: At the arising point, the given value is 1.57080 radians. As the earth moves from the arising point to the perihelion, this angle gradually decreases. After decreasing to the minimum value of 1.5551586 radians, it begins to gradually increase. When it reaches the perihelion, it continues to increase after passing through 1.57080 radians. When it increases to the maximum value of 1.5855145 radians, it begins to gradually decrease, and finally returns to the arising point, with the return value of 1.5706601 radians. The velocity value of the Earth's ahelion point is 29371.944m/s, the perihelion velocity is 30321.412m/s, and the return value of the Ahelion point is 29462.553m/s. The minimum rotation angular velocity of the Earth on the ecliptic surface produced by gravity is 6.022E-5 radian/s, the maximum value is 6.216E-5 radian/s, the return value is 6.040E-5 radian/s, and the average value is 6.122E-5 radian/s. 99.9992% of the sun's gravity on the earth is used for the earth's rotation.

Although the above simulation process is generally consistent with the actual situation, there are still many inconsistencies in the details, and the simulation data of the moon and the earth need to be further carefully adjusted.

6. Verification Method

The rotation of a planet formed by gravity can be verified by an artificial satellite that launches a solid rigid body. After the solid rigid body satellite enters orbit, no power or adjustment is added. At this time, the artificial satellite is only left with the effect of the earth's gravity, and then test whether the artificial satellite has rotation, whether the rotation direction is opposite to the direction of the revolution, and also test the rotation speed. Here, the artificial satellite will be affected by the sun and other planets.

7. In Conclusion

Gravity is the process in which nucleons emit gravitons and gravitons propagate in space with gravitational energy waves, and gravitational energy waves resonate with other nucleons and form energy transfer. For planets, gravitons emitted by nucleons inside the planet interact with other nucleons inside the planet, forming the cohesion of the planet. Gravitationaltons emitted by nucleons on the surface of the planet are partially emitted outside the ball. These gravitons emitted outside the ball meet the nucleons of other planets and resonate with them to form energy transfer. These transferred energy will cause the resonant nucleons to produce a vertical gravity line displacement, forming a vertical and gravitational line action force. The planets orbiting the central planet S, only spherical nucleons facing the central planet S can be received The gravitons emitted by the central planet S, the nucleus of the central planet S, will not receive the gravitons of the central planet S. In this way, for the entire planet, the equivalent gravitational action point Eg is not in the center of mass of the planet, but on the spherical surface of about 0.5 radius near the center of mass S. In this way, the central planet S will have two effects on the planet E that rotates around it. One is the centripetal force that rotates around the central planet S, and the other is the rotational force that rotates around the center of mass of the planet around the gravitational action point Eg. Within a certain time Δt, the initial velocity of planet E will cause the planet to move a uniform linear displacement. Planet E is subjected to the component force Fs of the gravity of the center planet S, which will cause the planet E to move for a distance in the arc. Planet E is subjected to another component force Fe of the gravity of the center planet S, which will cause the center of mass of planet E to move backwards on the arc for a distance around the gravitational point. Under the combined action of these three, planet E will form a standard elliptical orbit. According to the law of conservation of momentum, the linear velocity of a planet formed by gravity around the central planet is equal to the linear velocity of the planet rotation, and the angular velocity of the planet formed by gravity around the central planet is equal to the angular velocity of the planet’s rotation reflected to the central planet. At this time, for the planet’s revolution, the angular momentum generated by gravity cancels each other, leaving only the angular momentum formed by the initial velocity of the planet. This is the fundamental reason for the conservation of the angular momentum of the planet under the action of gravity. After derivation calculation and data simulation, more than 99.9% of the gravity of the central planet is used for the rotation of the planet. It can also be said that at the balanced inertia of the planet's movement, it can also be considered that at the point of gravity, the moment of the planet's revolution is equal to the moment of the planet's rotation. Since the force arm of the planet's revolution is much larger than the force arm of the planet's rotation, the force used for the rotation of the planet is much smaller than the force of the planet's rotation.