yade-users team mailing list archive

Thread
Date

Re: [Question #679516]: How to get the Volume in the simulation with periodic boundary condition?

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Jan Stránský <question679516@xxxxxxxxxxxxxxxxxxxxx>
Date: Mon, 25 Mar 2019 07:47:32 -0000
Reply-to: question679516@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
> when I output the position of all the particle (by state.pos) and
observed it in paraview, it become strange.

export not b.state.pos, but O.cell.wrap(b.state.pos), see [1]. state.pos
may have arbitrary value in space, likely outside the periodic cell.
O.cell.wrap shifts the positions to be something like cuboid.

> But I want it to like a cuboid, because I want to calculate the stress
inside and the unevenness of the boundary cause lots of troubles.

you can compute stress as getStress() [2] or PeriTriaxController.stress
[3]. The "boundary" (once more, no boundary exist) is irrelevant here.

cheers
Jan

[1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Cell.wrap
[2] https://yade-dem.org/doc/yade.utils.html#yade._utils.getStress
[3] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.PeriTriaxController.stress

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