yade-users team mailing list archive

[Luc Sibille <question271394@xxxxxxxxxxxxxxxxxxxxx>]

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.