Hello,
I am using the LIGGGHTS Flexible Fibers package in my work. In the Potyondy & Cundall 2004 paper where they introduce the parallel bond model implemented in LIGGGHTS, they describe a methodology (Appendix 1) to compute stress in the body arising from the bonds, which should be the same as the compute stress function in LIGGGHTS. From my experiments, it doesn't seem like the bond stresses are factored into the LIGGGHTS compute stress. This brings me to two questions:
1. Does anyone have any recommendations for computing stress arising from the bonds efficiently in the postprocessing? Would it be better to just update the compute to include bond forces?
2. In PFC, the forces and moments arising from the bonds are added to the contact forces which are used in combination with the external forces to update the particle positions. The absence of the bond stress in the compute makes me believe LIGGGHTS uses a different methodology. At what stage in the update are bond forces and moments used to update particle positions and angular velocities?
Thank you,
Lexi