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. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
