Phase field modeling of mixed-mode crack in rocks incorporating heterogeneity and frictional damage