yade-users team mailing list archive

[Milly Ro <question694837@xxxxxxxxxxxxxxxxxxxxx>]

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.