yade-users team mailing list archive

[Swapnil <question657126@xxxxxxxxxxxxxxxxxxxxx>]

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

Swapnil posted a new comment:
Hi Jan, now here is what I am trying.

from yade import pack,plot
pred=pack.inCylinder((0,0,0),(0,0,0.25),0.10)
sp=pack.randomDensePack(pred,radius=0.01,spheresInCell=300)
O.bodies.append(sp)
yade.qt.Controller()
idSteel=O.materials.append(FrictMat(young=210e9,poisson=.25,frictionAngle=.8,label="steel"))
hammer=sphere((0,0,0.30),0.04,material=idSteel)
hammerID=O.bodies.append(hammer)
hammer.state.blockedDOFs='z'
hammer.state.vel=Vector3(0,0,-0.005)
O.engines=[ForceResetter(),InsertionSortCollider([Bo1_Sphere_Aabb()]),InteractionLoop([Ig2_Sphere_Sphere_ScGeom()],[Ip2_FrictMat_FrictMat_FrictPhys()],[Law2_ScGeom_FrictPhys_CundallStrack()]),NewtonIntegrator(),PyRunner(iterPeriod=100,command="print hammer.state.displ()")]
O.dt=0.5*PWaveTimeStep()
O.saveTmp()

I understood what you explained regarding the interaction of the spheres and it is the reason behind the exploding phenomenon (repulsive forces).
When the body is alloted a velocity it should be seen moving, isn't it?

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