How can I define different cutoff distances for assigning bonds in the fix create/bond/gran command?

Submitted by Tannaz on Tue, 09/10/2024 - 01:59

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