PRODUCING A ROCK CONSIST OF DIFFERENT CLUSTERS (EXTENDING RICHTI83'S LIGGGHTS-WITH-BONDS)

Submitted by Tannaz on Tue, 03/19/2024 - 20:52

Hi everyone,

I am trying to simulate a rock under point load, similar to what is described in the LIGGGHTS-WITH-BONDS documentation in.bond. My model consists of different clusters, and I intend to start with two types of atoms (type 1 and type 2) that are randomly distributed within the model. Additionally, I want to define three types of bonds: bond1 between atom type 1 and atom type 1, bond2 between atom type 2 and atom type 2, and bond3 between atom type 1 and atom type 2. I have a couple of questions:

Is it possible to define two types of atoms that are randomly distributed in a model and then assign different properties to each of them using LIGGGHTS-WITH-BONDS?
My second question concerns the 'fix bondcr' command (fix bondcr all bond/create/gran 10 1 1 0.0051 1 6) used in bond.in. How does this command work?

Thanks
Tannaz