3D Discrete Element models on grain comminution in a granular material

This dataset contains the models and scripts used to create Figures 2 and 3 in the article Ioannidi et al., submitted. Impact of microscopic frictional parameters on grain comminution in a granular material: 3D Discrete Element Models. /repo/data contains six directories (cohesion_0 to cohesion_6), which are the models for each of the six cohesion values tested in the article. Each cohesion directory has six friction coefficient folders (frictionCoefficient_0 to frictionCoefficient_6), which are the models for each of the six friction coefficient values tested for each cohesion value (a total of 49 models). Each model has the following data: /repo/data/cohesion_*/friction_*/ Nbonds.csv: the total number of elastic bonds between grain aggregates per timestep out_force.csv: horizontal (x), vertical (y), horizontal (z) force measured on the top plate of the model, and the same x, y, z force components measured at the bottom plate /repo/data/cohesion_* contains the converted and combined numpy arrays of horizontal force (Fh.npy), vertical force (Fn.pny), total number of bonds (Nbonds.npy) and model names (Fname.txt) for each cohesion value. These are the data read by the script to produce the figures. /repo/Ioannidi_etal_2026_Figs2-3.py contains a python script that reads the numpy arrays, calculates the percentage of broken bonds during compaction and shearing, the macroscopic friction and saves them in the same directorty (file metrics.dat). It also produces figures 2 and 3 of the article (Fig2.pg and .eps, and Fig3.png and .eps, respectively).