Thanks Axel for the reply.
Do you know how to dump the interaction forces for granular pair_style.
It calculates the interactions during the simulation and also in compute stress/atom.
What is the problem in dumping these interaction forces?
it seems strange.
Regards,
Samantaray
Someone could write and add a single() function
to the granular pair styles. No one has done so
yet. The problem is
more conceptual about how to interpret what
they mean. Granular forces have history
effects, so they aren't unique to the current
configuration. And they are non-conservative.
And if you were the one who asked about
normal vs tangential forces, the interface to
the single() routine only returns 3 forces, not
6, so that would have to be worked around
somehow.