yade-users team mailing list archive

Thread
Date

Re: [Question #697432]: Why doesn't the strain equal to zero

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Jan Stránský <question697432@xxxxxxxxxxxxxxxxxxxxx>
Date: Mon, 07 Jun 2021 07:50:50 -0000
Reply-to: question697432@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
a) "modifying the position directly is likely to break Yade’s algorithms" [2]. O.cell.trsf is similar as setting position directly.
b) Some values are "propagated" to correct values during the iteration.

PeriTriaxController.strain is only updated inside the engine's action [3,4]. So you can:
a) try to run one more step to let triax be updated
b) set cell.velGrad to get trsf=I in the next step

Cheers
Jan

[2] https://yade-dem.org/doc/user.html#imposed-velocity
[3] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/dem/PeriIsoCompressor.cpp#L126
[4] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/dem/PeriIsoCompressor.cpp#L195

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.