post

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,num=14000,distributeMass=False)
#sp.makeCloud((-r_cyl,-r_cyl,-(h_cyl/2)),(r_cyl,r_cyl,h_cyl/2),rMean=3,rRelFuzz=.5,num=5000)#enforce positive mass
sp2 = pack.filterSpherePack(pred,sp,returnSpherePack=True)
sp2.toSimulation()
print (len(O.bodies)-1) #no of spheres generated
O.dt=0.005*PWaveTimeStep()