yade-users team mailing list archive

[Jan Stránský <question680334@xxxxxxxxxxxxxxxxxxxxx>]

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

Jan Stránský proposed the following answer:
Hi, sorry for late answer. A MWE:

###
O.materials.append(JCFpmMat(young=10e9,cohesion=1e20,tensileStrength=1e20))
O.engines = [
   ForceResetter(),
   InsertionSortCollider([Bo1_Polyhedra_Aabb()]),
   InteractionLoop(
      [Ig2_Polyhedra_Polyhedra_ScGeom()],
      [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
      [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM()],
   ),
   NewtonIntegrator(),
]
def poly(x0,**kw):
   pts = [(x0+x,y,z) for y in (0,1) for z in (0,1) for x in (0,1)]
   return polyhedron(pts,**kw)
ps = p1,p2,p3 = [poly(x,color=color,wire=False) for x,color in zip((0,1-1e-9,3),((1,0,0),(1,1,0),(0,0,1)))]
O.bodies.append(ps)
p3.state.vel = (-500,0,0)
O.dt = 1e-8
###

red and yellow particles are conented with a cohesive link, the blue one
is just an "impactor"

ScGeom6D is not implemented, but Ig2_Plyhedra_Polyhedra_ScGeom should be
a good inspiration if needed

HTH
Jan

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