yade-users team mailing list archive

[Klaus Thoeni <question696028@xxxxxxxxxxxxxxxxxxxxx>]

Question #696028 on Yade changed:
https://answers.launchpad.net/yade/+question/696028

    Status: Open => Answered

Klaus Thoeni proposed the following answer:
Hi Fabio,

welcome!

You can look at the code to see how the threshold is calculated [1], but
it's basically the value you define in CohFrictMat * r^2, which in your
case is 10. Applying 5+5 makes 10 (note that you are doing a dynamic
calculation, the force will oscillate between 0 and 10, plot the force
and you will see it). Now, I would suggest you fix one node and apply a
velocity at the other node and plot displacement vs normalForce in the
interaction, i.e. replace eps with fN =
O.interactions[0,1].phys.normalForce[0].

Klaus 
 
[1] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/dem/CohesiveFrictionalContactLaw.cpp#L313

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.