yade-users team mailing list archive

[Leonard <question683489@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,
When I add spheres to a simulation, sometimes I can not see all the particles (situation 1) in every direction, sometime I can only see all of them from some certain directions (situation 2), and sometimes a part of spheres are shown like a hole (situation 3), all these situations I have checked through "len(O.bodies) " to make sure that these spheres are existing. My yade version is 2018.02b with Ubuntu 18.04.
My script is as follow:
script1 which should be run firstly :
from yade import pack, plot
O.bodies.append(geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=31))
sp=pack.SpherePack()
sp.makeCloud((0,0,0),(1,1,1),rMean=.05,rRelFuzz=.5)
sp.toSimulation(color=[1,1,1])
Gl1_Sphere.quality=3
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),damping=0.4),
	PyRunner(command='checkUnbalanced()',realPeriod=2),

]
O.dt=.5*PWaveTimeStep()
O.trackEnergy=True
def checkUnbalanced():
	if unbalancedForce()<.05:
		O.pause()
        O.save('final.yade.bz2')
O.run()

script2 which should be run after script1, and you are expected to meet situation2 and situation3:
Gl1_Sphere.quality=3
O.load('final.yade.bz2')
for i in range(10):
	O.bodies.append(sphere((i,i,i), radius=0.5, color=[1, 0, 0]))

Thanks in advance!
Leonard

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