The estimation of accumulated local plastic strain combining with the Hertzian stress analysis can be further incorporated into professional multiscale failure modelling frameworks such as coarse-grained (CG) model [26] and extended finite element methods (XFEM) [27] to enhance their accuracy.