yade-users team mailing list archive

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

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
> I found that when there are only three polyhedrons, it can fall freely
normally

then the initial overlaps are most suspicious

> During the initial generation of particles, enough space is left so
that the particles do not touch each other.

In your original code, if I add this at the end
###
O.step()
for i in O.interactions:
    pd = i.geom.penetrationDepth
    b1,b2 = [O.bodies[id] for id in (i.id1,i.id2)]
    c1,c2 = [b.shape.__class__.__name__ for b in (b1,b2)]
    print(pd,c1,c2)
###
I got something like (depending on the run and makeCloud mood, sometimes there is no interaction, sometimes there are more interactions)
###
3.365690143044668e-05 PotentialBlock PotentialBlock
0.00021375937087687118 PotentialBlock PotentialBlock
###
meaning that you DO have some non-negligible initial interactions (and therefore the space left is NOT enough), disturbing the perfect free fall.
Under the presence of interactions, time step of course does have some influence.

cheers
Jan

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