yade-users team mailing list archive
Mailing list archive
Re: Strange behavior of periodic cell
> I'm trying to do some simulations using periodic cell. In one case, I prescribed some extension of the periodic cell, but the particles in the simulation behave really strange - the dispacement of particles si continuously increasing.
I had a look again. It looks as if kinetic energy is increasing, even
I tried with Law2_Dem3DofGeom_FrictPhys_Basic and also
Law2_ScGeom_FrictPhys_Basic (instead of the Cpm law) and the result is
always the same: spheres bumping into each other going crazy after a
short while. What is happening??
Can someone have a look as well? Especially tracking stored elastic and
kinetic energy could be useful. (there is no plasticity involved) The
script is attached (run with -j1, it might be faster).
from yade import pack, log, qt
# create material
# create 3x3x3 particles
for i,j,k in itertools.product(range(3),repeat=3):
## the same with other consitutive laws:
# forced compression of periodic cell along x
O.cell.trsf+=Matrix3(-.01,0,0, 0,0,0, 0,0,0)