← Back to team overview

yade-users team mailing list archive

[Question #697597]: PeriIsoCompressor.sigma different from getStress()


New question #697597 on Yade:


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,


## MWE ##
O.periodic = True
O.cell.setBox((1, 1, 1))

sp.makeCloud(minCorner=(0,0,0), maxCorner=(1, 1, 1), rMean=5e-3, rRelFuzz=5e-4, num=5000, periodic=True, seed=42)

iso_comp = PeriIsoCompressor(charLen=5e-3, stresses=[-1e5], doneHook="O.pause()")
O.engines = O.engines + [iso_comp]


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.