We present an assumed enhanced strain finite element framework for the simulation of tensile fracturing processes in transversely isotropic rocks. Fractures along the weak bedding planes and through the anisotropic rock matrix are treated with distinct enrichment, and a recently proposed dualmechanism tensile failure criterion for transversely isotropic rocks is adopted to determine crack initiation for the two failure modes. The cohesive crack model is adopted to characterize the response of embedded cracks. As for the numerical implementation of the proposed framework, both algorithms for the update of local history variables at Gauss points and of the global finite element system are derived. Four boundary-value problem simulations are carried out with the proposed framework, including uniaxial tension tests of Argillite, pre-notched square loaded in tension, three-point bending tests on Longmaxi shale, and simulations of tensile cracks induced by a strip load around a tunnel in transversely isotropic rocks. Simulation results reveal that the proposed framework can properly capture the tensile strength anisotropy and the anisotropic evolution of tensile cracks in transversely isotropic rocks. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

In recent decades, numerous geotechnical hazards, including landslides, foundation settlements, and tunnel collapses, have been linked to block-in-matrix rock (bimrock), resulting in substantial damage. The investigation and analysis of the engineering properties and mechanical behavior of both the rock and soil have become increasingly intriguing research areas. However, analyzing bimrock remains a formidable challenge due to its inherent heterogeneity. In this study, to investigate the tensile behavior of the bimrock, the numerical method that couples discontinuous deformation analysis (DDA) and smoothed particle hydrodynamics (SPH) is improved. First, a tension damage model is implemented in SPH for simulating the tensile behavior of the soil. The effectiveness of the presented model is verified through the direct tensile test model and the Brazil disc split model as an indirect tensile test. The contact algorithm of DDA-SPH is then modified by adding a tensile contact spring to introduce tensile strength at the interface between the matrix and the rock in the bimrock. Through a simple pulling numerical model, the accuracy of this modification has been verified, and the appropriate tensile contact stiffness is discussed. Furthermore, using the proposed method, the overall tensile strength of the bimrock with respect to the interface tensile strength is investigated. Finally, the proposed numerical method is applied to the simulation of geoengineering problems, demonstrating the capacity to analyze the stability and large deformation of bimrocks.