Hi,
I want to simulate the bonded particle model of a uniaxial compression test using a uniform distribution of particle radii. The radius range is from 0.8 to 2.5 millimeters. I would like to assign bonding properties to these particles using the fix create/bond/gran command, but I don’t want to use a fixed cutoff distance for assigning the bonding properties. Instead, I want the cutoff distance to be equal to the sum of the radii of each pair of neighboring or contacting particles. Does anyone have recommendations on how to do this? I was considering using the Python interface or modifying the source code, but since changing the source code would be quite difficult, I wanted to know if anyone has an easier suggestion.
Regards,
Tannaz