yade-users team mailing list archive

[Chiara Modenese <question206415@xxxxxxxxxxxxxxxxxxxxx>]

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

Chiara Modenese proposed the following answer:
Hi Jan,

Thanks for reply. The definition of average stress tensor in Bagi's
paper, for instance, is equal to the one in Yade [1]. I was not sure
about its validity for arbitrary grain shape (Bagi specifies that the
method is valid for grains of convex shape but I cannot see why his
derivation would not apply to concave shape too). I have adapted the
definition of branch vector in my local branch [2] but it only handles
the case of clump-clump interaction whereas it would be nice to have in
trunk all the possible combinations (sphere-clump, etc.).

[1] https://yade-
dem.org/doc/yade.utils.html?highlight=getstress#yade._utils.getStress

[2]  // compute branch vector
      if ( Body::byId(I->getId2(),scene)->isClumpMember() && Body::byId(I->getId1(),scene)->isClumpMember()) { // if both bodies are clump mumbers       
         int id_clump_2 = Body::byId(I->getId2(),scene)->clumpId;
         int id_clump_1 = Body::byId(I->getId1(),scene)->clumpId;
         branch = Body::byId(id_clump_2,scene)->state->pos + scene->cell->hSize*I->cellDist.cast<Real>() - Body::byId(id_clump_1,scene)->state->pos;
      } else { // only spheres
         branch= Body::byId(I->getId2(),scene)->state->pos + scene->cell->hSize*I->cellDist.cast<Real>() - Body::byId(I->getId1(),scene)->state->pos;
      }

Chiara

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.