← 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

Chiara Modenese proposed the following answer:

On 23 Jun 2012, at 20:01, Chareyre wrote:

> 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>();
Bruno, why is it -=? In PeriIso*.cpp it used to be +=. Am I missing something?
Thank you, Chiara

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

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