yade-users team mailing list archive

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

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

    Status: Open => Answered

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

yes, for this case Ig2_ Polyhedra_ Polyhedra_ ScGeom is not working correctly.
penetrationDepth is computed as radius1 + radius2 - norm.
So it is correct only if the interaction is "concentric" - true contact point lies between centers of the particles.
The more excentric the interaction is (the more center-branches are "tirangular"), the bigger is the error.
true contact point ... i.geom.contactPoint is set between the centers, so using the value as it is now would not help much
Probably it is worth to create an issue on Yade gitlab pages.

You can:
- use something like Ig2_ Wall Polyhedra_ ScGeom (which you would have to implement)
- adjust Ig2_ Polyhedra_ Polyhedra_ ScGeom
both cases needs change of the source code :-(

Cheers
Jan

[1] https://gitlab.com/yade-
dev/trunk/-/blob/master/pkg/dem/Polyhedra_Ig2.cpp#L565

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