yade-users team mailing list archive

[Sanny <question267328@xxxxxxxxxxxxxxxxxxxxx>]

New question #267328 on Yade:
https://answers.launchpad.net/yade/+question/267328

Hello,


I just started with this package.  My goal is to generate random (non-overlapping) spheres in accordance with particle size distribution.  I want to achieve a volume fraction of spheres equal to approximately 70%.  I followed the online manual and came up with the following code:

# 25 mm cube, with vol frac of spheres=0.7
from yade import pack, export, qt
from math import *
Lx=25e-3 
sp=pack.SpherePack()
sp.makeCloud(minCorner=(0,0,0),maxCorner=(Lx,Lx,Lx),porosity=0.3,psdSizes=[2e-3,4e-3,6e-3,8e-3,10e-3,12e-3,14e-3,16e-3],psdCumm=[0.01,0.23,0.4,0.6,0.75,0.93,0.98,1])
sp.toSimulation()
qt.Controller() 
qt.View()
export.text('spheres')

I have no problems with the run, I get the output file called 'spheres'.  But when I calculate the volume fraction of spheres within the 25 mm cube, I find it to be incredibly small (less than 20%) and the number of particles is way too small.

I am sure I am missing something in the makeCloud that is giving me such low values.  I would have expected at least 400 plus particles.

Would any of you be so kind to help me out with the above code.

The screen output looks like this:

-------------
Welcome to Yade 1.07.0 
TCP python prompt on localhost:9000, auth cookie `secuda'
XMLRPC info provider on http://localhost:21000
Running script test1.py
/usr/lib/python2.7/dist-packages/IPython/frontend.py:30: UserWarning: The top-level `frontend` package has been deprecated. All its subpackages have been moved to the top `IPython` level.
  warn("The top-level `frontend` package has been deprecated. "
[[ ^L clears screen, ^U kills line. F12 controller, F11 3d view (use h-key for showing help), F10 both, F9 generator, F8 plot. ]]
--------------

The spheres file has the following data:

#format x_y_z_r
0.00340762	0.00955894	0.018026	0.00334986
0.00696167	0.0172939	0.0189378	0.00417147
0.0111412	0.00602221	0.00908746	0.00583889
0.019339	0.0165512	0.00921438	0.00474594
0.0205392	0.0178629	0.0191574	0.00424171
0.0147442	0.0202884	0.00583614	0.00177721
0.0199264	0.00781992	0.0112529	0.00128409
0.0170094	0.00858523	0.0191854	0.00560053
0.022655	0.00516024	0.0125429	0.00145099
0.012582	0.0153632	0.00598169	0.00256647


Thank you,
Sanny

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.