Numerical Simulation of Crack Propagation and Damage Behavior of Glass Plates under Impact Loading

Crack initiation and propagation is a long-standing difficult problem in solid mechanics, especially for elastic-brittle material. To explore the damage and crack propagation behavior of glass plates under impact loading, the element deletion, discontinuous Galerkin peridynamic (DG-PD), and meshless peridynamic (M-PD) methods are used to conduct numerical simulations, respectively. The JH-2 material model, and the maximum principal stress and maximum principal strain failure criteria are adopted in the element deletion method. The node separation operation and the critical energy release rate criterion are used in the DG-PD method. In the M-PD method, a self-programmed particle discretization method is utilized along with an appropriate computational domain, and a critical elongation criterion is imposed. The simulation results show that: (1) the element deletion method can roughly simulate the damage morphology of glass under impact loading, but it is insufficient in capturing crack bifurcation and penetration. (2) In the DG-PD method, circumferential cracks and radial cracks are observed, and the cracks are of high symmetry. In addition, there are a lot of glass fragments splashing at the impact point and the frame. (3) Radial cracks and circumferential cracks can be captured in the M-PD method, and the symmetry of the cracks is good. The size of horizon and the impact velocity show great influence on the dynamic responses of the glass plates. As far as the damage form is concerned, the M-PD method and the DG-PD method yield consistent results.