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. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
