yade-users team mailing list archive

[Jan Stránský <question630281@xxxxxxxxxxxxxxxxxxxxx>]

Question #630281 on Yade changed:
https://answers.launchpad.net/yade/+question/630281

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi Tina,

see [1], Jerome's answer #4. The problem is that after InteractionLoop,
the particles are moved by NewtonIntegrator, so i.geom.penetrationDepth
is the old value and setting it as i.phys.unp is nto enough.

One option is setting O.dt=0 for the first step (don't know how it works
with GlobalStiffnessTimeStepper, probably should be deactivated for the
first step?)

O.dt = 0
for i in O.interactions:
  i.phys.unp = i.geom.penetrationDepth
O.dt = ...

then I had no explosion

cheers
Jan

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.