Rock fracture in geo-materials is a complex phenomenon due to its intrinsic characteristics and the potential external loading conditions. As a result, these materials can experience intricate fracture patterns endowing various cracking phenomena such as: Branching, coalescence, shielding, and amplification, among many others. In this article, we present a numerical investigation concerning the applicability of an original bulk-interface fracture simulation technique to trigger such phenomena within the context of the phase field approach for fracture. In particular, the prediction of failure patterns in heterogenous rock masses with brittle response is accomplished through the current methodology by combining the phase field approach for intact rock failure and the cohesive interface-like modeling approach for its application in joint fracture. Predictions from the present technique are first validated against Brazilian test results, which were developed using alternative phase field methods, and with respect to specimens subjected to different loading case and whose corresponding definitions are characterized by the presence of single and multiple flaws. Subsequently, the numerical study is extended to the analysis of heterogeneous rock masses including joints that separate different potential lithologies, leading to tortuous crack paths, which are observed in many practical situations.Ministerio de Economía y Competitividad MAT2015-71036-
