yade-users team mailing list archive

[Jérôme Duriez <question670487@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jérôme Duriez proposed the following answer:
Hi,

Resulting force and displacement of YADE bodies can be respectively accessed with :
- O.forces.f(yourBodyId), see [1,2]
- O.bodies [yourBodyId].state.displ(), assuming O.bodies[yourBodyId].state.refPos was correctly initialized (this may be up to you). See [3]


The "problem" is : in your uniaxial configuration using UniaxialStrainer engine you do not have a single boundary element whose such variables would be representative of the mechanics of the sample (there is no loading rigid platen)

Maybe you could use sum/average the above variables on all
UniaxialStrainer.posIds bodies but you have to check if this makes
sense.


Another possibility is to surround your numerical samples with e.g. YADE boxes that would serve as rigid loading platens (2 surrounding boxes may be enough for uniaxial conditions) as done in [Duriez2016] (see Fig. 3) from https://yade-dem.org/doc/publications.html, and elsewhere. Then you could directly use O.forces.f() and O.bodies[].state.displ()


A last possibility would be to re-derive with continuum mechanics the force and displacement from the stress and strain quantities you already know ;-)


[1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Omega.forces
[2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.ForceContainer.f
[3] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.State.displ

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