<p>This study investigates the equidistant rupture phenomenon in layered rock masses, which significantly influences the mechanical integrity and safety of geotechnical structures. An advanced peridynamic approach is employed to address the limitations of classical method in accurately modeling the stress distribution and the effects of damage at crack sites. The proposed model incorporates non-local differential operator and strain energy density criterion for strain softening, establishing a link between mechanical parameters at crack tips and residual strain energy. The GPU parallel computing is incorporated to significantly enhance computational efficiency, enabling large-scale processing of peridynamic simulations. The results reveal that the ratio of crack spacing to layer thickness is identified as a critical determinant of the overall mechanical state of the layered rock mass. Additionally, increasing surface layer thickness reduces the number of saturated cracks while increasing the average networked crack area. Conversely, thinner layers exhibit a denser crack distribution, underscoring the significant influence of layer thickness on fracture dynamics. This study emphasizes the importance of understanding fracture mechanisms in layered rock masses to enhance the design and the safety of engineering structures. In general, this research contributes to the field of rock mechanics by providing a robust computational framework for analyzing fracture behavior in layered materials, thereby facilitating better predictions of structural performance under various loading conditions. The findings underscore the significance of considering material interfaces and evolving mechanical properties in the assessment of fracture processes in geotechnical applications.</p>