← Back to team overview

yade-users team mailing list archive

Re: [Question #199704]: getStress for ScGeom

 

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

Chareyre proposed the following answer:
Something like this in getStress:

if (b1->clumpMember && b2->clumpMember){
l = O.bodies[b2->clumpId]->pos - O.bodies[b1->clumpId]->pos
} else {
//like now
}

It needs also to adapt this to clump vs. standalone.

For the periodic case, it needs to be shifted, but the current line should work without any change:
if (isPeriodic) branch-= scene->cell->hSize*I->cellDist.cast<Real>();

You can use triaxial test on clumps to validate (if boxes are used for
boundaries, it gives an alternate calculation of stress). There is an
example script from Janek IIRC.

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