yade-users team mailing list archive

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

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

Jan Stránský proposed the following answer:
> the hammer does not appear to move at all

but it does :-) you can easily check by
O.engines += [PyRunner(iterPeriod=100,command="print hammer.state.displ()")]

> There is no motion of the hammer

see above

> Also in this case there has not been any force applied to the spheres
in the cylinder that would cause them to move.

not true.No external force was applied, but If two spheres are
overlapping, the default contact law computes repulsive force between
them. You can visually see the overlaps before the simulations starts.
That's why it "explodes". Use a different pacing strategy, like

###
...
sp=pack.randomDensePack(..., spheresInCell=200) # decreases overlaps *almost* to zero
...

O.step()
calm()
maxPD = max(i.geom.penetrationDepth for i in O.interactions) # get max overlap
factor = 1 - maxPD / 0.01
print maxPD, factor
for b in O.bodies:
   growParticle(b.id,factor) # and decrease size of each particle to be sure they do not overlap
###

cheers
Jan

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