Question #695457 on Yade changed:
https://answers.launchpad.net/yade/+question/695457
Status: Open => Answered
Jérôme Duriez proposed the following answer:
Hi,
I did not look so closely to your code in your text but, for what concerns the
line at [*] the "(1 / phys->ks) * (trialForce - shearForce)" term gives the
plastic/irreversible component of the increment in shear displacement.
(Not the total one, in an elasto-plastic wording)
It is then just fair to multiply it with current force to get a purely
dissipated work.
(and, in case that was your question, there is no need to speak about dashpot
;-) )
[*]
https://gitlab.com/yade-dev/trunk/-/blob/7ab9457c6e69c6e381b3e58ed39665ece60baf81/pkg/dem/ElasticContactLaw.cpp#L91
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
