yade-users team mailing list archive
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:
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  and using periodic cell dimensions , 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..
You received this question notification because your team yade-users is
an answer contact for Yade.