yade-users team mailing list archive

[Jan Stránský <question267328@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
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
>

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