yade-users team mailing list archive

[Jan Stránský <question696515@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hello,

Ig2 is only one component in the interaction creation chain. Ig2_ Sphere_ Sphere_ ScGeom.interactionDetectionFactor "just" considers also "distant interactions".
The [1] check would return false (= no interaction) for the case of no overlap, but it does the computation and returns true (= interaction for further consideration) for interactionDetectionFactor>1.
It computes penetrationDepth [2] as 
Real norm = normal.norm(); // distance of sphere centers
Real penetrationDepth = s1->radius + s2->radius - norm;
so for "distant interactions" i.geom.penetrationDepth is negative.

Then, also Ip2 and Law2 are evaluated.
E.g. for "classical" Law2_ScGeom_FrictPhys_CundallStrack, the interaction is deleted for the case of negative penetrationDepth [3] (i.e. "the Law2 functor in use permitting" in [4]).

cheers
Jan

[2] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/dem/Ig2_Sphere_Sphere_ScGeom.cpp#L51
[3] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/dem/ElasticContactLaw.cpp#L62
[4] https://yade-dem.org/doc/user.html#creating-interactions

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