Hi there,
I was just thinking about clumps in the periodic cell. When we compute
the stress tensor, we loop over all interactions, adding their force and
branch vector, then dividing by the cell volume. Now with clumps there
are invisible internal "forces", since the clump is rigid, but we don't
take intra-clump "interactions" in account. It follows that the stress
tensor will be artificially lower. How can I prescribe stress then?
Does someone have a suggestion how to make the computation "right"? Or
am I seeing an inconsistency where there is none?
Cheers, Václav
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp