← Back to team overview

yade-users team mailing list archive

Re: [Question #678710]: measure force in periodic boundary


Question #678710 on Yade changed:

    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..


[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.