This work presents a detailed three-dimensional finite element based model for wave propagation, combined with a postprocessing procedure to determine the fracture intensity caused by blasting. The data generated during this project includes output files of all simulations with detailed fields, geometries and meshes. The model incorporates the Johnson-Holmquist-2 constitutive model, which is designed for brittle materials undergoing high strain rates and high pressures and fracturing, and a tensile failure model. Material heterogeneity is introduced into the model through variation of the material properties at the element level, ensuring jumps in strain. The algorithm for the combined Johnson-Holmquist-2 and tensile failure model is presented and is demonstrated to be energy-conserving, with an open-source MATLABTM implementation of the model. A range of sub-scale numerical experiments are performed to validate the modelling and postprocessing procedures, and a range of materials, explosive waves and geometries are considered to demonstrate the model's predictive capability quantitatively and qualitatively for fracture intensity. Fracture intensities on 2D planes and 3D volumes are presented. The mesh dependence of the method is explored, demonstrating that mesh density changes maintain similar results and improve with increasing mesh quality. Damage patterns in simulations are self-organising, forming thin, planar, fracture-like structures that closely match the observed fractures in the experiments. The presented model is an advancement in realism for continuum modelling of blasts as it enables fully three-dimensional wave interaction, handles damage due to both compression and tension, and relies only on measurable material properties. The uploaded data are the specific simulation outputs for four explosion models occurring on two different rock types, and the specific fracture patterns generated.