yade-users team mailing list archive

[Robert Caulk <question659177@xxxxxxxxxxxxxxxxxxxxx>]

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

Robert Caulk posted a new comment:
> the very cause of the issue was the degree of compaction of the
generated sphere pack.

It is worth it for me to correct the thread conclusion, maybe for your
sake, but more importantly for Yade documentation's sake. The cause of
the problem is not the degree of compaction of the sphere pack. It does
not require any special "calm()" functions to work properly, you just
need to use the correct functors and engines. Here is the fully
functional script:

from yade import pack

container=geom.facetCylinder((0,0,0), 1.0, 4.0, segmentsNumber=20, wallMask=6)
O.bodies.append(container)
pred=pack.inCylinder((0,0,0), (0,0,1), radius=.99)
sp=pack.randomDensePack(pred, spheresInCell=2000, radius=.05)
O.bodies.append(sp)

O.engines=[
	ForceResetter(), 
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
	[Ig2_Sphere_Sphere_ScGeom(),
		Ig2_Facet_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	NewtonIntegrator(gravity=(0,0,-9.81)),
]
O.dt=PWaveTimeStep()

from yade import qt
qt.Controller()
qt.View()
O.saveTmp()

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