Question #291326 on Yade changed: https://answers.launchpad.net/yade/+question/291326
Status: Answered => Open Luis Barbosa is still having a problem: Hi Jan, Sure I can. #!/usr/bin/python # -*- coding: utf-8 -*- #====================================== from yade import pack import gts from yade import utils surf=gts.read(open('tetrahedron.gts')) pred=pack.inGtsSurface(surf) vol = surf.volume() #====================================== # Spheres O.materials.append(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,label='sphereso')) #====================================== O.bodies.append(pack.gtsSurface2Facets(surf,color=(0,0,1))) sp = pack.randomDensePack(pred,radius=0.08,rRelFuzz=0.01, material = 'sphereo',useOBB=False,returnSpherePack=True) sp.toSimulation(color=(0.9,0.8,0.6)) #====================================== O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()], [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)], [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=False)] ), GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.2), NewtonIntegrator(damping=0.9,gravity=[0,0,-9.81]), ] ve=utils.getSpheresVolume() print vol,ve #======================================= from yade import qt qt.View() qt.Controller() ######################################################### And here follows the gts surface... https://drive.google.com/open?id=0B9LhrwAf_vasanZ4bVc4b0FHVUE Thanks, Luis -- 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@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp