yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #25780
[Question #697597]: PeriIsoCompressor.sigma different from getStress()
New question #697597 on Yade:
https://answers.launchpad.net/yade/+question/697597
Hello,
I am trying to compact a cloud in a periodic cell using PeriIsoCompressor.
At first it seems to work well : the cloud gets denser until PeriIsoCompressor.sigma.mean() reaches the target mean stress. But when I compare PeriIsoCompressor.sigma to getStress(), the mean stress is completely different (~26 times lower).
Please see below a mwe and its output when using yade at git revision 4852d5601.
Am I correctly using PeriIsoCompressor ?
Thanks in advance,
Sacha
## MWE ##
O.periodic = True
O.cell.setBox((1, 1, 1))
sp=pack.SpherePack()
sp.makeCloud(minCorner=(0,0,0), maxCorner=(1, 1, 1), rMean=5e-3, rRelFuzz=5e-4, num=5000, periodic=True, seed=42)
sp.toSimulation()
iso_comp = PeriIsoCompressor(charLen=5e-3, stresses=[-1e5], doneHook="O.pause()")
O.engines = O.engines + [iso_comp]
O.run(wait=True)
print("\nPIC mean stress: ", iso_comp.sigma.mean())
print("getStress mean stress: ", getStress().trace()/3)
## Output with yade 4852d5601 ##
Welcome to Yade 2021-06-08.git-4852d56
Using python version: 3.6.9 (default, Jan 26 2021, 15:33:00)
[GCC 8.4.0]
TCP python prompt on localhost:9000, auth cookie `yecssk'
XMLRPC info provider on http://localhost:21000
Running script mwe.py
PIC mean stress: -99503.28987122637
getStress mean stress: -3824.356053577036
--
You received this question notification because your team yade-users is
an answer contact for Yade.