Summary
The aim of this work is to efficiently implement the Park‐Paulino‐Roesler cohesive zone model with the objective of creating realistic high‐resolution simulations of material deformation, fracture, and postfracture behavior. Intrinsically, unstructured meshes can create more realistic fracture patterns in bulk material than structured meshes. Implicit methods, stable for much larger time steps, have greater potential to model both fracture and postfracture behavior without sacrificing speed of execution. Several technical contributions are presented, including (i) GPU‐accelerated implementation of the Park‐Paulino‐Roesler cohesive zone model, (ii) efficient creation of sparse matrix structure, and (iii) comparison of different unloading/reloading relations when using an implicit scheme. A potential‐based collision response scheme was implemented that allows one to model the interaction of fragmented material. Several test simulations are carried out to demonstrate the flexibility of the model and its ability to reproduce different materials under various loading conditions. Benchmarking results show that most of the computational time is spent by the third‐party solver library, meaning that other operations do not require optimization. The library is made available as open source.