yade-users team mailing list archive

[behzad <question259096@xxxxxxxxxxxxxxxxxxxxx>]

New question #259096 on Yade:
https://answers.launchpad.net/yade/+question/259096

Hi,

I'm running the following script. It makes an orthogonal packing in a cylindrical shape. I'm wondering why the system has such overlaps which gives the particles an explosive velocity. In this case I need to run the calm() function to avoid sample collapsing. 
Even with calm() function running at every iteration, if we choose O.dt=1.0e-4, particles move. 
Why this is happening? 


O.reset()
from yade import utils, plot
from yade import pack, qt

O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.0),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]),
NewtonIntegrator(damping=0.7,gravity=[0,0,0])
]

Co=O.materials.append(FrictMat(young=1e7,poisson=0.3,density=1377,frictionAngle=1))

pred=pack.inCylinder((0.02,0,0),(0.02,0,0.06),0.02)

O.bodies.append(pack.regularOrtho(pred,radius=0.001,gap=0.0, material=Co))



-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.