The use of glass as structural material has highlighted the need for more reliable numerical approaches to analyze its mechanical behavior, especially in the accidental eventuality of fracture. Modelling the behavior of fractured laminated glass, in fact, is fundamental to assess the Post-Fracture load-bearing capacity. However, this is a highly challenging task because of the many interplaying factors, such as the viscoelastic and thermal-dependent behavior of the interlayer, the presence of a highly complex and variable crack pattern and the interaction among fragments. The objective of the present work is the development and testing of a robust numerical model that can naturally introduce the generated crack pattern into virtual specimens and manage the interaction among many fragments. The phase field fracture model is herein explored, by assigning the damage variable to fit the pre-existing crack pattern. Then, the specimen is loaded letting the phase field managing the fragments interaction. The dependence of the stress tensor with the damage variable is herein defined through the Cleavage-Deviatoric model, since it prevents fully damaged regions from transmitting tensile and shear stresses yet keeping their ability to bear compressive forces. Indeed, this model can asymptotically reproduce unilateral and frictionless contact conditions between the existing crack lips. Preliminary case studies are discussed to check the potentiality of the proposed approach.