New question #694837 on Yade:
https://answers.launchpad.net/yade/+question/694837
Hello,
I want to pack spheres in a box and then to apply gravity. For this I use my
own particle distribution. My question is how can I increase the fill factor?
My ideas were to increase the dimension of the box in order to generate more
spheres. But it did not help. Then I wanted to add the porosity parameter in
makeCloud function and again no significant changes in fill factor. I also
variied the timestep. It helps a little bit, but I am still far from the
desired result. The code that I used for packing and gravity applying and the
formula for fill factor you can find below:
################
from yade import export
# create box with free top, and ceate loose packing inside the box
from yade import pack, plot
import numpy as np
psdSizes = np.array(np.genfromtxt("r.txt", dtype=float, unpack=False))*2
psdCumm = np.genfromtxt("cdf.txt", dtype=float, unpack=False)
O.bodies.append(geom.facetBox((50,50,50),(50,50,50),wallMask=63))
sp=pack.SpherePack()
sp.makeCloud((0,0,0),(100,100,100),porosity=.65,
psdSizes=psdSizes,psdCumm=psdCumm, distributeMass=False)
sp.toSimulation()
O.engines=[
ForceResetter(),
# sphere, facet
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.1),
]
#O.dt=utils.PWaveTimeStep()
O.dt=0.02
O.run(1000)
###############
The fill factor is calculated as ratio of the sum of the volume of all
generated spheres to the volume of the box. The maximum achieved value was 46%.
The desired value is 65-75%
Best wishes
Milly
--
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
