yade-users team mailing list archive

[José Mercado <question268936@xxxxxxxxxxxxxxxxxxxxx>]

New question #268936 on Yade:
https://answers.launchpad.net/yade/+question/268936

Hi everybody!

I have problems trying to eliminate a facet. The geometry that i need is the intersection of a facetbox with a facetcylinder. With the wallMask command inside the facetCylinder i removed the top "circle" of the cylinder but i have problem with the bottom wall of the box. The balls cant enter in the cylinder because i need a "circular hole" in the bottom wall of the box. There is a possibility of generating a geometry through the intersection of the cylinder and the box? (like predicates).

I tried the following code:

from __future__ import division
from yade import pack
from yade import plot,qt
import math,numpy

### SPHERES 

sp=pack.SpherePack()
sp.makeCloud((-5,-5,-5),(5,5,5),rMean=0.05,rRelFuzz=.5,num=10000)
sp.toSimulation()

### MATERIAL OF SPHERES AND WALLS
 
O.materials[0].frictionAngle=atan(0.5)
O.materials[0].poisson=.5
O.materials[0].young=1e6
O.materials[0].density=1000

id_wall = O.materials.append(FrictMat(density=10000,young=1e11,poisson=.5,frictionAngle=atan(0.5),label="wallmat"))
O.materials[id_wall]

### WALLS

O.bodies.append(yade.geom.facetCylinder((0,0,-6),radius=1,height=2,segmentsNumber=20,wallMask=6))

O.bodies.append(yade.geom.facetBox((0,0,0),(5,5,5)))

# ENGINE

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()],
   ),
   NewtonIntegrator(damping=0.4,gravity=(0,0,-9.81)),
   # call the checkUnbalanced function (defined below) every 2 seconds
]
 
O.dt=utils.PWaveTimeStep()




-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.