yade-users team mailing list archive

[j0hnn1e <question672684@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

j0hnn1e is still having a problem:
Hello Robert,
>>In an effort to keep the forums organized for future Yade users, please take a moment to review these guidelines [1] (note the external links rules).
I am sorry for using external link. Please find the script below:
from yade import pack, export, plot
from yade import utils, ymport
import math,numpy
#from yade import plot
r_cyl=75
h_cyl=250
#psdSizes,psdCumm=[2.36,4.75,6.3,9.5,12.5,19,25],[.365,.512,.596,.77,.864,.974,1] #Category 4, N100
psdSizes,psdCumm=[2.36,4.75,6.3,9.5,12.5],[0.347,0.657,0.76,0.973,1] #Category 3, N70
 
 
cyl1=O.bodies.append(geom.facetCylinder((0,0,0),r_cyl,h_cyl,wallMask=6))
fill_size=100
pred=pack.inCylinder(centerBottom=(0,0,-(h_cyl/2)+1),centerTop=(0,0,fill_size),radius=r_cyl-3)
idAggregate=O.materials.append(FrictMat(density=2600,young=30e9,poisson=.4,frictionAngle=1.0472,label="aggregate")) #frictionangle in rad
sp=pack.SpherePack()
sp.makeCloud((-r_cyl,-r_cyl,-(h_cyl/2)),(r_cyl,r_cyl,h_cyl/2),psdSizes=psdSizes,psdCumm=psdCumm,distributeMass=False)
sp2 = pack.filterSpherePack(pred,sp,returnSpherePack=True)
sp2.toSimulation()
print (len(O.bodies)-1) #no of spheres generated
O.dt=0.005*PWaveTimeStep()
######################
I tried using 3/ method in [1] but I am still getting quite high %age of spheres of radius less than 2.36 (0.347 is input, 0.654 is the output). I have checked by toggling the bool variable 'distributeMass' (True and False) and I still get very high answers in both the cases. 


1. https://yade-dem.org/doc/yade.pack.html#yade._packSpheres.SpherePack.makeCloud

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