yade-users team mailing list archive

[Luis Barbosa <question693834@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Luis Barbosa is still having a problem:
Thank you Jan,

I was able to create the pack using:

sp = pack.randomDensePack(...., returnSpherePack=True)
sp.toSimulation()

However, the particles inside the pack are not even touching each other.

Before an actualiyzation of YADE (Change from Python2 to 3), I was using
the same script to creat dense pack with overlaping particles. Much
different from what I am getting now using the same parameters though.


Here E you have the link to the GTS surface:
https://drive.google.com/file/d/1w_r_q438jZAqzgL5DaQEprZr9Mroazq6/view?usp=sharing


Here you have the MWE:

from yade import plot
import gts,yade.pack
######################################Material#################################################
O.materials.append(JCFpmMat(type=1,young=70e4,poisson=0.3,frictionAngle=radians(30),density=0.0025,tensileStrength=194,cohesion=194,jointNormalStiffness=0,jointShearStiffness=0,jointCohesion=194,jointFrictionAngle=radians(30),jointDilationAngle=0.0,label='spheres'))
######################################Selector#################################################

surf=gts.read(open('agg18.gts'))
pred=pack.inGtsSurface(surf)
sp=pack.randomDensePack(pred,radius=0.025,returnSpherePack=True)#memoizeDb='/tmp/gts-triax-packings.sqlite'spheresInCell=200,
#O.bodies.append(sp)
sp.toSimulation()

O.bodies.append(pack.gtsSurface2Facets(surf,wire=True))

######################################Engine#################################################
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, Key="Wei", recordCracks=True)]
	),

#  TranslationEngine(ids=p2,translationAxis=[0,0,-1],velocity=0.03),
  GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=500,timestepSafetyCoefficient=0.5),
  NewtonIntegrator(damping=0.1,gravity=[0,0,-981]),#cm/s2
]

######################################Exporting#################################################

from yade import qt
qt.View()
qt.Controller()
#export.text("agg18.txt")

#exit()

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