yade-users team mailing list archive

[Bettina Suhr <question290933@xxxxxxxxxxxxxxxxxxxxx>]

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

Dear all,

I would like to generate a packing of clumps inside a box, where at front and back of the box periodic boundary conditions are applied. At first, I defined the periodic cell and then added a facet box (without the front and back side) inside the cell. The problem is that when the clumps fall down due to gravity, some clumps go through the bottom of my facet box. I think that this problem is not related to the young modulus of the material or the time step size: When I define my box with facets at front and back, the problem does not occur. Can anybody explain me, what is the problem? Below you find my minimal working example.

Thanks a lot,
Bettina


from yade import pack

O.materials.append(FrictMat(young=1e10, poisson=0.28, frictionAngle=0.2, density=2600.))
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(), Bo1_Wall_Aabb(), Bo1_Facet_Aabb()],allowBiggerThanPeriod=True),
	InteractionLoop(
	  [Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom(), Ig2_Wall_Sphere_ScGeom()],
	  [Ip2_FrictMat_FrictMat_FrictPhys()],
	  [Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	NewtonIntegrator( gravity=(0,0,-9.81),label='newtonEng')# NO DAMPING 
]

Len=700e-3
Wid=120e-3
H=450e-3+100e-3
O.periodic=True
O.cell.setBox(Wid,Len+1e+0,H+1e+0)

#add box
#wall mask: -x (1), +x (2), -y (4), +y (8), -z (16), +z (32).
O.bodies.append(geom.facetBox(center=(Wid/2.,Len/2.+0.5e+0,H/2.+0.5e+0),extents = ( Wid/2.,Len/2.,H/2.), wallMask=3+4+8+16+32, color=(1,0,0)))

#create pack of clumps in box 
sp=pack.SpherePack()
Num=100
c1=pack.SpherePack([ ((0.,0.,0.),37.6e-3/2.), ((37.6e-3/2.,0.,0.),25e-3/2.) ])# overlap (!) between both spheres
sp.makeClumpCloud( (0,0+0.5e+0,0+0.5e+0), ( Wid,Len+0.5,H+0.5), [c1],num=Num,seed=1)
ballastID=sp.toSimulation()
O.bodies.updateClumpProperties(discretization=10)# correct mass, volume, inertia!!


O.dt=0.5*PWaveTimeStep() 



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