Dimensions of packing
Hi All,
I have a code that generate a cloud of spheres in a cylinder facet and then compact the spheres to a certain porosity. After compaction, I would like to know the height of the compacted spheres. Can you please give some suggestions on how to do that?
Below is my code.
Thank you,
Othman
from yade import plot,pack, export, ymport
targetp = .3 ##specify the targeted porosity
##specimen geometry
radiuscyl=.05
heightcyl=.203*4
SphereRadius = .0045
# material parameters
O.materials.
#######
sp=pack.
sp.makeCloud(
#### cylinder extraction
pred=pack.
spFilter=
spFilter.
print ('porosity = ', utils.porosity())
#######
facets=
cylinder=
yade.qt.View()
##creating disks
d1=geom.
d2=geom.
disk1IDs= O.bodies.append(d1)
disk2IDs= O.bodies.append(d2)
for i in disk1IDs:
body= O.bodies[i]
body.state.vel = (0,0,-10)
for n in disk2IDs:
body= O.bodies[n]
body.state.vel = (0,0,10)
#######
O.dt=.5*
enlargeFactor=1.5
O.engines=[
ForceResetter(),
InsertionSortC
Bo1_Sphere_
Bo1_Facet_Aabb()
]),
InteractionLoop(
[
Ig2_
Ig2_
],
[
Ip2_
Ip2_
],
[
Law2_
],
),
#PyRunner(
NewtonIntegrat
PyRunner(
PyRunner(
]
# reset interaction detection enlargement
bo1s.aabbEnlarg
def P():
print ('porosity = ', utils.porosity())
def stop():
if utils.porosity(
O.pause()
print 'Finished'
for i in disk1IDs: O.bodies.erase(i)
for i in disk2IDs: O.bodies.erase(i)
for i in cylinder: O.bodies.erase(i)
print ('unbalanced forces = ', utils.unbalance
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Othman Sh
- Solved:
- Last query:
- Last reply: