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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
