Phase-Field Fracture Modelling of Thin Monolithic and Laminated Glass Plates under Quasi-Static Bending

A phase-field description of brittle fracture is employed in the reported four-point bending analyses of monolithic and laminated glass plates. Our aims are: (i) to compare different phase-field fracture formulations applied to thin glass plates, (ii) to assess the consequences of the dimensional reduction of the problem and mesh density and refinement, and (iii) to validate for quasi-static loading the time-/temperature-dependent material properties we derived recently for two commonly used polymer foils made of polyvinyl butyral or ethylene-vinyl acetate. As the nonlinear response prior to fracture, typical of the widely used Bourdin–Francfort–Marigo model, can lead to a significant overestimation of the response of thin plates under bending, the numerical study investigates two additional phase-field fracture models providing the linear elastic phase of the stress-strain diagram. The typical values of the critical fracture energy and tensile strength of glass lead to a phase-field length-scale parameter that is challenging to resolve in the numerical simulations. Therefore, we show how to determine the fracture energy concerning the applied dimensional reduction and the value of the length-scale parameter relative to the thickness of the plate. The comparison shows that the phase-field models provide very good agreement with the measured stresses and resistance of laminated glass, despite the fact that only one/two cracks are localised using the quasi-static analysis, whereas multiple cracks evolve during the experiment. It was also observed that the stiffness and resistance of the partially fractured laminated glass can be well approximated using a 2D plane-stress model with initially predefined cracks, which provides a better estimation than the one-glass-layer limit.