New question #672531 on Yade: https://answers.launchpad.net/yade/+question/672531
Hello everyone, I am trying to generate a pack of spheres for a slope made up of granite rock using this script : from yade import pack, export, ymport import gts, os.path, locale def mat(): return JCFpmMat(type=1,young=1e8,poisson=0.3,frictionAngle=radians(30),density=3000,tensileStrength=1e6,cohesion=1e6,jointNormalStiffness=1e7,jointShearStiffness=1e7,jointCohesion=1e6,jointFrictionAngle=radians(20),jointDilationAngle=0.0) #### control parameters mesh='iit' #name f gts mesh sizeRatio=10. # defines discretisation (sizeRatio=meshLength/particleDiameter) #### import mesh locale.setlocale(locale.LC_ALL,'en_US.UTF-8') #gts is locale-dependend. If, for example, german locale is used, gts.read()-function does not import floats normally surface=gts.read(open(mesh+'.gts')) print 'closed? ', surface.is_closed() #### generate packing if surface.is_closed(): pred=pack.inGtsSurface(surface) # get characteristic dimensions aabb=pred.aabb() dim=pred.dim() center=pred.center() minDim=min(dim[0],dim[1],dim[2]) # define discretisation radius=minDim/(2*sizeRatio) print center, dim, ' | minDim=', minDim, ' | diameter=', 2*radius O.bodies.append(pack.randomDensePack(pred,radius=radius,rRelFuzz=0.3,useOBB=True,memoizeDb='/tmp/gts-triax-packings.sqlite',returnSpherePack=False,color=(0.9,0.8,0.6),material=mat)) O.bodies.append(pack.gtsSurface2Facets(surface,color=(0.8,0.8,0.8),wire=True)) #### export packing# in a text file export.text(mesh+'_'+str(int(sizeRatio))+'.spheres') But it is showing me this warning and any results are not generating Best-fit oriented-bounding-box computed for GTS surface, orientation is Quaternion((1,-6.131644428895787e-125,-2.9761932039277854e-113),4.641198870949644e-197) WARN /data/trunk/pkg/dem/SpherePack.cpp:104 makeCloud: porosity must be >0, changing it for you. It will be ineffective if rMean>0. WARN /data/trunk/pkg/dem/Shop_01.cpp:445 PWaveTimeStep: PWaveTimeStep has not found any suitable spherical body to calculate dt. dt is set to 1.0 -- 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 : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
