Nonlinear fracture analysis of rubber-like materials is computationally challenging due to a number of complicated numerical problems. The aim of this paper is to study finite strain fracture problems based on appropriate enrichment functions within the extended finite element method. Two-dimensional static and quasi-static crack propagation problems are solved to demonstrate the efficiency of the proposed method. Complex mixed-mode problems under extreme large deformation regimes are solved to evaluate the performance of the proposed extended finite element analysis based on different tip enrichment functions. Finally, it is demonstrated that the logarithmic set of enrichment functions provides the most accurate and efficient solution for finite strain fracture analysis. FINITE STRAIN FRACTURE ANALYSIS BY NEW XFEM ENRICHMENT FUNCTIONS 1317 Dolbow [9,10]; nevertheless, there is still lack of reliable commercial codes on the finite strain XFEM method, and the few available software [1,2] are extremely limited, and can only provide reliable results for very simple linear elastic large strain benchmarks [1]. Moreover, most of the existing commercial codes have serious inaccuracies in predicting the singularity of the stress field at crack tips, as studied by Gigliotti and Kroon [1], because the crack tip singularity enrichment functions are not included in the XFEM formulation. Consequently, the crack tip is forced to lay on the element edges for propagation purposes [1].Many researchers have attempted to use XFEM within the large deformation regime; among them Dolbow and Devan [11], Legrain et al. [12,13], Shen and Lew [14], Mergheim et al. [15], Fagerstrom and Larsson [16], Anahid and Khoei [17], and Khoei et al. [18] have presented various aspects of geometrically nonlinear XFEM analyses. Legrain et al. [12] proposed the first extended finite element stress analysis of the crack tips for plane stress conditions and investigated the importance of the right choice for the singularity enrichment functions. Moreover, they pointed out that due to the blunting process, neglecting the tip enrichments would be more acceptable than using inappropriate functions [12]. In this study, therefore, the goal is set to evaluate different available analytical solutions, to derive appropriate crack tip enrichment functions, and to assess their performance for analysis of a number of complex problems.Developing dedicated enrichment functions for various problems have been the subject of XFEM research in the past decade to improve the accuracy of XFEM for new evolving problems. Sukumar et al. [19] developed the tip enrichment functions from the analytical solution near an interface crack tip in isotropic bimaterials. Asadpour and Mohammadi [20] and Asadpour et al. [21] presented two different sets of orthotropic crack-tip enrichment functions, while Hattori et al. [22] proposed an alternative formulation to derive the orthotropic enrichment functions. Later, Esna Ashari and Mohammadi [23] derived the necessary enrichment functions...