yade-users team mailing list archive

[zhou qian <question683187@xxxxxxxxxxxxxxxxxxxxx>]

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

zhou qian posted a new comment:
Hi Robert Jan:
Here are my codes:

from libyade import yade
from yade import *
from yade import utils,polyhedra_utils
from minieigen import Vector3, Matrix3, Matrix6

# basic materials
matS = PolyhedraMat(young=1.e14,poisson=0.1,frictionAngle=0.2,density=1e7)
O.materials.append(matS)

POLY1=[
Vector3(-1.1252229043696602e-9,0.04044867969995926,-0.02500000014298641),
Vector3(-2.404624131217237e-10,0.0499995,-0.015723519897073306),
Vector3(0.003975620285304186,0.04082183933765579,-0.018135369824275445),
Vector3(0.005226423189199525,0.04972559,-0.02500000012427156)]

POLY2=[
Vector3(1.222315982854378e-10,0.0499995,-0.015723520027142578),
Vector3(-0.008409782917328579,0.03956479014801928,-0.02500000011219652),
Vector3(-0.006184036292741935,0.03689338966056221,-0.018402069868769105),
Vector3(-0.010395579966696341,0.0489073797989524,-0.015723519982404767)]

poly1 = utils.polyhedron(material=matS, vertices=POLY1,fixed=True)
poly2 = utils.polyhedron(material=matS, vertices=POLY2,fixed=True)
O.bodies.append(poly1)
O.bodies.append(poly2)

# basic YADE engines
O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Facet_Aabb(),Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb()],verletDist=0.01),
   InteractionLoop(
      [Ig2_Wall_Polyhedra_PolyhedraGeom(), Ig2_Polyhedra_Polyhedra_PolyhedraGeom(), Ig2_Facet_Polyhedra_PolyhedraGeom()], 
      [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], # collision "physics"
      [Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()]   # contact law -- apply forces
   ),

   NewtonIntegrator(damping=0.1,gravity=(0,0,0.0)),
]

O.dt = 0.00001

O.step()


POLY1 and POLY2 contain the coordinates of the 2 polys when it got
error. It got error in line 271-289 in function Simplify in
/pkg/dem/Polyhedra_support.cpp


I use yade-2017.01a on Ubuntu16.04, but I can't get my CGAL version.
CGAL is installed in 2015, so I think its version>4.5. I didn't enable
lbmflow linsolv pfvflow sph which don't have any influence on polyhedron
mudule.

Thanks for your patience!
Yours,
Joe

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