**4. Conclusions**

In this study, the time partial derivative of the elastic wave equation was discretized using the central difference scheme, and the spatial partial derivative was discretized using a meshless GFDM. When realizing the elastic wave numerical simulation via the GFDM in the time domain, the Poisson disk node generation algorithm and the centroid Voronoi node adjustment algorithm were combined to obtain an even and random distribution of nodes, and the Cerjan damping boundary condition was introduced. Through numerical calculation, the following understandings and conclusions were obtained:


**Author Contributions:** Conceptualization, S.L. and Z.Z.; methodology, S.L.; software, S.L.; validation, S.L.; formal analysis, W.Z.; review and editing, W.Z. and Z.Z.; supervision, Z.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the Fundamental Research Funds for the Central Universities of Central South University 2019ZZTS300.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data that support the findings of this study are available from the corresponding author upon reasonable request.

**Acknowledgments:** We are grateful for technical support from the High-Performance Computing Center of Central South University. We gratefully acknowledge the editors and reviewers for providing thoughtful and useful suggestions.

**Conflicts of Interest:** The authors declare no competing financial interests.
