Experimental and numerical failure mechanism evaluation of anisotropic rocks using extended finite element method

Rocks often show an anisotropic behavior due to microcracks, mineral-based structure, and unequal in-situ stresses. The crack tip's stress intensity factors (SIFs) that control rock failure heavily depend on anisotropy. Given the significance of extensive research to identify the effects of anisotropy on the mechanical behavior and strength of rocks, numerical modeling is considered an essential method in designing structures and analyzing their stability. This study experimentally and numerically evaluated the effects of the bedding angle, initial crack angle with the loading axis, crack length, and anisotropy ratio on the stress intensity factor of transversely isotropic rocks via the extended finite element method (XFEM). Cracked chevron-notched Brazilian disks (CCNBDs) were fabricated from phyllite under the International Society for Rock Mechanics (ISRM) standards and tested at different initial crack lengths and angles for validation. The XFEM approach was then used to predict toughness and crack propagation path in the tested and new specimens with varying ratios of anisotropy and bedding angles (ψ = 0, 30, 45, 60, and 90°). The results demonstrated that with the increase of the anisotropy ratio, in equal crack length and initial angle, the absolute value of the SIFs have larger values. In addition, with the increase of crack length, the first and second modes of SIF increased, which indicates the failure of the specimen at lower loads. Also, bedding angle is effective on the stress intensity factors and this effect becomes more intense as the crack length increases. Finally, the anisotropy properties not only affect the values of SIFs, but also can change the angle at which the specimen experiences pure mode I or mode II.