Triaxial Test on a cylinder

Hello,

> I am a newbie in yade

welcome :-)

please be more specific on "generate a cylinder".
If it means to extract particles within a cylindrical predicate, you can try:
#############
from yade import pack
pred = pack.inCylinder((0.5,0.5,0),(0.5,0.5,1),.4) #[1]
cyl = [b for b in O.bodies if isinstance(b.shape,Sphere) and pred(b.state.pos,b.shape.radius)]
# predicate(center,radius) returns True if it contains the sphere(center,radius), False otherwise

# test
sphs = [(b.state.pos,b.shape.radius) for b in cyl]
O.reset()
O.bodies.append([sphere(pos,r) for pos,r in sphs])
#############

if you need more help to understand the example (it is a bit compressed commands, could be rewritten to more understandable way), let us know
cheers
Jan

[1] https://yade-dem.org/doc/yade.pack.html#yade._packPredicates.inCylinder

Yes, I meant a cylindrical predicate. Now, where I have to put this snippet of code or modify the existing one in the trunk repository?

you place/execute the snippet just after your commands :-)

#############
TriaxialTest(numberOfGrains=500).load()
O.run(1000,True) # True argument is important, prevent further commands to be executed before 1000 iterations are finished

from yade import pack
pred = pack.inCylinder((0.5,0.5,0),(0.5,0.5,1),.4) #[1]
cyl = [b for b in O.bodies if isinstance(b.shape,Sphere) and pred(b.state.pos,b.shape.radius)]
# predicate(center,radius) returns True if it contains the sphere(center,radius), False otherwise

# test
sphs = [(b.state.pos,b.shape.radius) for b in cyl]
O.reset()
O.bodies.append([sphere(pos,r) for pos,r in sphs])
#############

cheers
Jan