yade-users team mailing list archive

[Chareyre <question672684@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Chareyre proposed the following answer:
Hi,
I would suggest not using titles like «faulty something» since in most
cases problems come from faulty users. ;)

You distribute by mass then you check how it is distributed by number. Of
course it is not the same.

Next time please send the script as a whole.

Bruno


Le ven. 24 août. 2018 20:52, j0hnn1e <question672684@xxxxxxxxxxxxxxxxxxxxx>
a écrit :

> New question #672684 on Yade:
> https://answers.launchpad.net/yade/+question/672684
>
> Hello All,
> I am trying to generate cylindrical sphere packing according to a given
> PSD. First, I use the following script to generate sphere packing:
> >
> psdSizes,psdCumm=[2.36,4.75,6.3,9.5,12.5],[0.347,0.657,0.76,0.973,1]
> fill_size=100
>
> pred=pack.inCylinder(centerBottom=(0,0,-(h_cyl/2)+1),centerTop=(0,0,fill_size),radius=r_cyl-3)
> sp=pack.SpherePack()
>
> sp.makeCloud((-r_cyl,-r_cyl,-(h_cyl/2)),(r_cyl,r_cyl,h_cyl/2),psdSizes=psdSizes,psdCumm=psdCumm,num=14000,distributeMass=True)
> sp2 = pack.filterSpherePack(pred,sp,returnSpherePack=True)
> sp2.toSimulation()
> >
> Then after running the script, I check the % of spheres with radii less
> than 2.36 as follows:
> >
> yell=[]
> for b in O.bodies:
>   if isinstance(b.shape,Sphere):
>     yell.append(b.shape.radius)
> count=0
> for j in range(len(yell)):
>   if yell[j]<2.36:
>     count=count+1
> print (float(count)/len(yell))
> >
> Using latter method, I get my answer around 0.87 which should actually be
> around 0.347. Is there a mistake in my code or in the usage of makeCloud
> command? I am not able to figure out.
>
> Thanks in advance!
>
> --
> 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     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>
>

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.