yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #11252
Re: [Question #267328]: Volume fraction of spheres appear to be far too less using Makecloud
Hi Sanny,
see documentation of makeCloud [1]: "Create random loose packing ..."
to create denser packing, have a look at randomDensePack [2] or
randomPeriPack [3].
cheers
Jan
[1]
https://yade-dem.org/doc/yade.pack.html#yade._packSpheres.SpherePack.makeCloud
[2] https://yade-dem.org/doc/yade.pack.html#yade.pack.randomDensePack
[3] https://yade-dem.org/doc/yade.pack.html#yade.pack.randomPeriPack
2015-05-24 20:36 GMT+02:00 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.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>
Follow ups
References