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.

Intense surface eruptions are observed along the curved surface of a confined cylindrical film of hydrogel subject to laser-induced converging-diverging shock loading. Detailed numerical simulations are used to identify the dominant mechanisms causing mechanical instability. The mechanisms that produce surface instability are found to be fundamentally different from both acoustic parametric instability and shock-driven Richtmyer-Meshkov instability. The time scale of observed and simulated eruption formation is much larger than that of a single shock reflection, in stark contrast to previously studied shock-driven instabilities. Moreover, surface undulations are only found along external, as opposed to internal, soft solid boundaries. Specifically, classic bubble surface instability mechanisms do not occur in our experiments and here we comment only on the new surface undulations found along the outer boundary of solid hydrogel cylinders. Our findings indicate a new class of impulsively excited surface instability that is driven by cycles of internal shock reflections.