A hybrid potential of mean force approach for simulation of fracture in heterogeneous media