“…Therefore, to overcome analytical limitations, numerical methods have been used to simulate the stress distribution and/or crack initiation and propagation of (transversely) isotropic rocks for the Brazilian test (see Li and Wong [2013] for a brief review). For example, these numerical efforts include the displacement continuity for investigating the influence of the loading contact on the initiation and propagation of fractures [e.g., Malan et al, 1994], the boundary element method for evaluating the effect of the initiated cracks on stress distribution within isotropic rock samples [e.g., Lanaro et al, 2009], the finite element method for 3-D stress analysis to show size/shape effects under the Brazilian test [e.g., Yu et al, 2006], and for analyzing failure process of heterogeneous isotropic rocks under static/dynamic loading [e.g., Zhu and Tang, 2006]. Recently, the numerical methods have been extended to utilize the extended finite element method for simulating crack propagation in the cracked Brazilian disc [e.g., Eftekhari et al, 2015], the discrete element method for the mechanical behavior of transversely isotropic rocks using embedded smooth joints (e.g., Park and Min [2015] in 2-D and Park et al [2016] in 3-D), and the hybrid finite element method-discrete element method for analyzing the influence of microscale heterogeneity and microcracks on the brittle fractures [e.g., Mahabadi, 2012] and for modeling the transition from continuum to discontinuum during fracturing process [e.g., An et al, 2016].…”