The Discontinuous Deformation Analysis (DDA) method is parallel to the finite element method and is a numerical integration method for solving the deformation and motion analysis of discrete block systems. In the process of applying the DDA method to celestial body simulation, using a constant acceleration time integration scheme for long-term simulation will accumulate errors. This paper applies the DDA method of linear acceleration method and average acceleration method in the Newmark integration scheme for celestial body simulation. The calculated results are compared with the data of the virtual planetarium software Stellarium. The improved DDA method achieves error reduction in the celestial body simulation process, verifying the effectiveness and practicality of the improved DDA method in celestial body simulation.