yade-dev team mailing list archive

[Bruno Chareyre <1301443@xxxxxxxxxxxxxxxxxx>]

@Jerome
Unfortunately, your hyperlink now shows the fixed version... (and mine too in #6).
[1] points to the hold one.
And no, the fact that axis!=Vector3r::Zero() does not imply that axis.norm()>0, since MIN_REAL² < MIN_REAL.

[1]
https://github.com/yade/trunk/blob/6319de511e8118da4ca120ee1725f34060186a2e/pkg/dem/NewtonIntegrator.cpp#L235


** Changed in: yade
       Status: New => Fix Committed

-- 
You received this bug notification because you are a member of Yade
developers, which is subscribed to Yade.
https://bugs.launchpad.net/bugs/1301443

Title:
  Spheres/Box contact : spheres ori suddenly set to NaN.

Status in Yet Another Dynamic Engine:
  Fix Committed

Bug description:
  Hi,

  Please see and launch the attached script.
  Some spheres are falling into a box because of gravity. The simulation starts well, but when all particles are almost steady, wait for a few seconds and you will see that some of them who are in contact with the bottom wall will disapear. I spent some time to investigate and that's my preliminary conclusions :
  - it is fully reproducible
  - particles don't really disappear : their orientation is suddenly turned to Quaternion(Nan,Nan,Nan,2*pi). It's not linked with a high torque on them. All other attributes remains okay, and the other particles still can interact with it. If you manually set the orientation to a correct value after it disappears, it will appear again at the right place.
  - you can reproduce this with one sphere (uncomment the last line of the script and comment the previous paragraph) : let it fall then apply to it a little angVel : "O.bodies[-1].state.angVel=(0.,0.,0.001)". It will disappear after a few timesteps.
  - it's independent on the timestep value.
  - it's dependant on the interpenetration of the contact : if you increase the gravity or decrease the young modulus, you will be able to avoid this behaviour.

  Does somebody have an idea ? I think it's a rotation calculation
  issue, in a very specific geometric case.

To manage notifications about this bug go to:
https://bugs.launchpad.net/yade/+bug/1301443/+subscriptions