This is a modified version of LIGGGHTS that includes granular bond formulations. The original version of LIGGGHTS produced by CFDEM can be found here: https://github.com/CFDEMproject/LIGGGHTS-PUBLIC
Hello,
That would be great.
Any help is appreciated.
Some things that I am looking at are
Refactoring bond_gran.cpp to move the bond-breaking code into its own class
work on removing the entire fiber once one sphere leaves the domain (started this in fix_remove_molecule.{cpp, h})
Look into "smooth" fibers and alter atom mass such that there is little difference between a fiber with "kissing spheres" and a fiber with highly overlapped spheres.
Redoing the jkr contact model to remove the look-up-table and replace it with a newton solver (its faster...)
have the bond model work (correctly) for superquadratic particles (this would probably be a thesis...)
Add a new contact model to model hollow fibers (Yu Guo has a paper on this)
There's just a few things and I am open to anything else that you would be passionate about.
Hi @schrummy14
I would like to contribute to your code. There is something I can do to help?
Regards.