Question #271394 on Yade changed:
https://answers.launchpad.net/yade/+question/271394
Status: Open => Answered
Luc Sibille proposed the following answer:
Hi,
1/ I think you're right.
2/ It should be
M_t = k_t * O_t - twisting
M_b = k_t * O_t - bending
where O_t is the relative rotation of particles at contact from the time of
contact creation.
However I have not checked in the code. So when you want to be sure what the
code is doing, the best and safest way to be sure is to look directly in the
code. For instance about the latter point you can read the code here:
http://bazaar.launchpad.net/~yade-pkg/yade/git-trunk/view/head:/pkg/dem/CohesiveFrictionalContactLaw.cpp#L173
3/ No there is angular velocities but M_t and M_b are computed according to 2/
see:
http://bazaar.launchpad.net/~yade-pkg/yade/git-trunk/view/head:/pkg/dem/CohesiveFrictionalContactLaw.cpp#L164
Best,
Luc
--
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
