yade-users team mailing list archive

[Jérôme Duriez <question692146@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jérôme Duriez proposed the following answer:
Hello,

The answer in C++ language is at [1].

In a maybe more human language, the interactions that carry the forces
entering the stress expression are assumed to have a geom [2] component
of GenericSpheresContact [3] (or child) type (if not, YADE would
probably crash)

As such, there necessarily is a reference radius assigned to each
interacting body, see refR* [4]. For wall or facets, the radius seems to
be chosen as once or twice the sphere's one [5], hence the volume.

How much sense this makes might be a next question ;-)

[*] https://gitlab.com/yade-dev/trunk/-/blob/b87224b638ed162dd8865d097332286b4b2c9263/pkg/dem/Shop_02.cpp#L885, most recent version as of now
[2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Interaction.geom
[3] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.GenericSpheresContact
[4] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.GenericSpheresContact.refR1 and refR2
[5] https://gitlab.com/yade-dev/trunk/-/blob/b87224b638ed162dd8865d097332286b4b2c9263/pkg/dem/Ig2_Facet_Sphere_ScGeom.cpp#L122 for Facet-Sphere and https://gitlab.com/yade-dev/trunk/-/blob/b87224b638ed162dd8865d097332286b4b2c9263/pkg/dem/Ig2_Facet_Sphere_ScGeom.cpp#L220 for Wall-Sphere

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