The discontinuous deformation analysis(DDA)method is widely used in the field of geotechnical engineering,among which the three-dimensional spherical particle DDA(3D-SDDA)theory can numerically simulate the motion of spherical object models.The author innovatively applies it to the simulation of the operation of celestial bodies in the solar system,treating each celestial body as a spherical object,only considering the gravitational force between celestial bodies and their own inertial force.By establishing the overall balance equation of the system,deriving the load matrix equation,solving the overall equation,and writing a Matlab calculation pro-gram,the motion parameters of each celestial body are obtained,and the numerical simulation of celestial body motion is carried out.The results indicate that the use of 3D-SDDA theory can simulate the operation of the sun-earth system and the sun-earth-moon system in space,and can obtain the spatial position and velocity information of celestial bodies at any time.The numerical calculation results of the sun-earth system were compared with the Kepler calculation method,and the results of the sun-earth-moon system were compared with the high-precision data of the virtual observatory software Stellarium,and the results were in good agreement.Simultaneously inputting the limiting conditions for the occurrence of a total lunar eclipse can accurately predict the timing of the eclipse.
discontinuous deformation analysis(DDA)motion of celestial bodynumerical simulationsun-earth systemsun-earth-moon system