yade-users team mailing list archive

[Jan Stránský <question678710@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi James,

using periodic boundary conditions (I prefer to call it periodic contact detection), the sample is periodic and there is no surface..
You can compute average stress [1] and using periodic cell dimensions [2], you can extract something like "surface force":
###
stressTensor = getStress()
s_xx = stressTensor[0,0]
s_yy = stressTensor[1,1]
s_zz = stressTensor[2,2]
dx,dy,dz = O.cell.refSize
a_xx = dy*dz
a_yy = dx*dz
a_zz = dx*dy
f_xx = s_xx / a_xx
f_yy = s_yy / a_yy
f_zz = s_zz / a_zz
###

but **personally** I would use just the stress values..

cheers
Jan

[1] https://yade-dev.gitlab.io/trunk/yade.utils.html#yade._utils.getStress
[2] https://yade-dev.gitlab.io/trunk/yade.wrapper.html#cell

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