yade-users team mailing list archive

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

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi James

> But I found there is no periodic option in facetCylinder.

just make the cylinder open and make it "long enough".
Also consider using boxes instead of facets.

A MWE. Notice, that the sphere fly away after a few bounces (read [1])
###
from yade import geom
r = 3 # cyl radius
h = 2 # cyl "height"
e = 3 # cyl height multiplicator to capture "overflown" spheres
#
cyl = geom.facetCylinder((r,r,.5*h),r,e*h,segmentsNumber=24,wallMask=4,wire=False)
O.bodies.append(cyl)
O.periodic = True
O.cell.setBox(10,10,h)
#
rs = 0.5
s = sphere((1.1*rs,r,0),rs)
s.state.vel = 1e-2*Vector3(1,1,.5)
O.bodies.append(s)
#
O.engines[1].allowBiggerThanPeriod = True # or set it to collider in O.engines definition
# (!) Read [1] carefully
# or divide the cylinder also along axis not to have large bodies
###

cheers
Jan

[1] #https://yade-
dem.org/doc/yade.wrapper.html#yade.wrapper.InsertionSortCollider.allowBiggerThanPeriod

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