In this paper, a thermomechanical coupled phase field method is developed to model cracks with frictional contact. Compared discrete methods, the can represent arbitrary crack geometry without an explicit representation of surface. The two distinguishable features proposed are: (1) for mechanical phase, no specific algorithm needed imposing contact constraints on fracture surfaces; (2) thermal ...