yade-users team mailing list archive

[veerawat <question688643@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Open

veerawat gave more information on the question:
Dear jan,

this's link

https://www.youtube.com/watch?v=b_yLp0onOzg

i want to put the compressin like this video for my script :

from yade import pack, plot

# create rectangular box from facets
O.bodies.append(geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=63))

# create empty sphere packing
# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry

O.materials.append(FrictMat(density=1000,young=1e4,poisson=0.3,frictionAngle=radians(30),label='sphereMat'))

# sphere packing is not equivalent to particles in simulation, it contains only the pure geometry
sp=pack.SpherePack()
# generate randomly spheres with uniform radius distribution
O.materials.append(FrictMat(density=1000,young=1e4,poisson=0.3,frictionAngle=radians(30),label='sphereMat'))

sp.makeCloud((0,0,0),(1,1,1),rMean=.03,rRelFuzz=.01,num=10000,seed=10000)
# add the sphere pack to the simulation
sp.toSimulation()

O.engines=[
 ForceResetter(),
 InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
 InteractionLoop(
  # handle sphere+sphere and facet+sphere collisions
  [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),
 # call the checkUnbalanced function (defined below) every 2 seconds

]
O.dt=.5*PWaveTimeStep()

#
O.saveTmp()

thank for helping me

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