I'm posting what I have done to incorporate damping into a flexible fiber model.
I am using Yu Guo's paper "Granular Shear Flows of Flexible Rod-like Particles"
This release currently works for LIGGGHTS 3.7.0
The only thing that i would like to explore further is this bond model interaction with the fix_template_multiplespheres command.
As this would allow the user to insert directly into the simulation rather than read in the particle and bond information from a text file.
I will also be uploading a FORTRAN code that allows the insertion of flexible fibers into a cylindrical region.
Here is a link to richti83's original bond code
Here is my link