yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #06152
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.