Blind prediction of curved fracture surfaces in gypsum samples under three-point bending using the Discontinuous Galerkin Cohesive Zone method

A computational framework based on a Discontinuous Galerkin (DG)/Cohesive Zone formulation is utilized to simulate the experiments of the Purdue Damage Mechanics Modeling Challenge. The inelastic response of the additively-manufactured gypsum material used in the experimental tests is modeled via a dilatational plasticity model. The constitutive and fracture model parameters are calibrated using the load–displacement curves corresponding to three-point bending tests initially provided by the Challenge organizers. The test samples contained initial notches especially designed to force specific types of mixed fracture modes. The calibrated computational modeling framework is used to blindly simulate the more complex configuration of the Challenge experiments. The numerical predictions of the load–displacement curve and the shape of the curved fracture surface are compared to the experimental data provided a posteriori. It is found that the computational method is able to quantitatively describe the fracture response of the material including crack propagation, plastic wake, and the curved geometry of the fracture surface that results from the evolving fracture mode mixity with significant fidelity.