“…The method requires a fine mesh along the crack path and a suitable definition for regularization parameters (see a discussion in [24] and recent internal length-insensitive formulations in [25,26]). However, due to the above-mentioned advantages, the phase field method has been widely developed and applied to many problems, such as, among many others: brittle fracture [23,27,28], composite delamination [29], dynamic fracture [30][31][32], hydraulic fracture [33][34][35][36], topology optimization for resistance to cracking [3,4], anisotropic material fracture [37][38][39], ductile fracture [40][41][42][43][44], ductile/fragile transition [45,46], fracture in micro tomography image-based models of microstructures [5,6,47] and more recently adapted in machine learning strategies in [48].…”