yade-users team mailing list archive

[Nait <question212480@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I have encountered something strange when conducting gravity deposition: some particles will go outside the box.
Can you give some suggestions for solving this problem?

Thanks.

Nait


My script is as following:

####SCRIPT####

from yade import pack


rm = 0.8
rr = 0.6
compFricDegree = 30
finalFricDegree = 30 

stabilityThreshold=0.01

young=5e6
mn,mx=Vector3(0,0,0),Vector3(20,80,200)
mxp = Vector3(20,80,150)
mnp = Vector3(20,80,90)

thick = 0.01
targetporosity = 0.5

O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls'))


walls=utils.aabbWalls([mn,mx],thickness=thick,material='walls')
wallIds=O.bodies.append(walls)


sp=pack.SpherePack()
sp.makeCloud(mn,mxp,rMean=rm,rRelFuzz=rr,porosity=targetporosity)


O.bodies.append([utils.sphere(center,rad,material='spheres') for center,rad in sp])


O.dt=.5*utils.PWaveTimeStep() # initial timestep, to not explode right away
O.usesTimeStepper=True


O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
   InteractionLoop(
      # handle sphere+sphere and facet+sphere collisions
      [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()]
   ),
   GravityEngine(gravity=(0,0,-9.81)),
   NewtonIntegrator(damping=0.6),
]

O.run(2000, True)

while 1:
  O.run(500, True)
  unb=unbalancedForce()
  print 'unbalanced force:',unb
  if unb<0.01:
    break

yade.qt.Controller(), yade.qt.View()

####END OF SCRIPT####



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