Going back to the dynamic case (energy.py) :
It can be verified that there is an unbalanced energy evolution
independant of friction and damping.
To make plastic dissipation 0, it is enough to set frictionAngle close
to pi/2.
For damping, make it 0.
The results are attached. ScGeom and L3 give
Hi,
For info : there is a small fix of stress definition in periodic BCs in
revision 2589 that
slightly modify energy results.
- Before fix, the (normalized) diff between total energy and boundary work was
around 5e-3
(0.5%) for a compression εyy=0.1 (~60k iterations), as in attachment of
Hi Vaclav,
This energy test contains a wall. I realize it can link to a
long-lasting (I showed you in Prague) but unreported bug. I just
reported (https://bugs.launchpad.net/yade/+bug/675955). Conserving
energy in this situation would have been surprising.
Considering this bug, and also this
Bruno is right, if you exclude friction between spheres and wall, energy
will be conserved. Moreover, the same happens if you include boxes rather
than walls and you set friction equal to zero for boxes only (see
http://www.mail-archive.com/yade-us...@lists.launchpad.net/msg02395.html).
In this
This energy test contains a wall. I realize it can link to a
long-lasting (I showed you in Prague) but unreported bug. I just
reported (https://bugs.launchpad.net/yade/+bug/675955). Conserving
energy in this situation would have been surprising.
Look more carefully, it uses
Hi Vaclav,
why not trying to calculate the velocities at the right time? I remember I
was doing it and energy was conserved (without damping, either local or
global). We just need to shift them by adding acc*dt/2 (btw, we need to
_subtract_ this term in order to shift in the right way, as
6 matches
Mail list logo