Particle packing - How to avoid overlap
I am generating a pack of spheres (lets call them aggregates), using makeCloud. Each of those spheres consists of a packing of smaller particles glued together using cohesive material model. I am using randomDensePack to generate this packing of particles within each sphere.
Now, when I visualize the packing, I notice that sometimes there is an overlap between the particles of two different aggregates. I want to avoid this, since it results in internal stresses in the initial condition. Is there a better way to generate this kind of packing ?
Here is the minimal working example script:
from yade import pack,qt,
from math import *
#######
# Set up run
#######
run_name=
data_root_dir="."
#######
# Materials
#######
sample_
young=4e9
,poisson=0.25
,density=1400
,frictionAng
,normalCohes
,shearCohesi
,momentRotat
,label=
O.materials.
#######
# Component dimensions and operating condition
#######
# Sample dimension
sample_radius=5e-4
# Sub-particle dimension
particle_
#######
# Geometry
#######
xExt=4*
yExt=4*
zExt=xExt
xLim=xExt
yLim=yExt/4
zLim=zExt
#######
# Particle Packing
#######
min_corner= (0,yLim,0)
max_corner= (xLim, yExt-yLim, zLim)
sp=pack.
sp.makeCloud( min_corner,
print "Generated ",len(sp)," aggregates"
#######
# Sample
#######
for s in sp:
sphere=
sp1=pack.
sphere
,spheresIn
,radius=
,memoizeDb
,returnSph
)
sp1.toSimulat
print 'Generated ',len(sp1),' particles'
Gl1_Sphere(
Thanks
Varun
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask VG for more information if necessary.