Triaxial test using polyhedral particles

Hello,

I am new user of YADE and I’m trying to do a Triaxial test using polyhedral shapes.
To prepare the sample I created a cloud of spheres and then replace them with polyhedral particles.
Therefore, my question is:
1- Is it possible to use polyhedral shapes using periodic boundary conditions?
If you have any suggestions about using polyhedral particles please tell me! (the following is my code)
sp0=pack.SpherePack();
NbParticles=sp0.makeCloud(mn,mx,psdSizes=[0.125,0.160,0.2,0.25,0.315,0.4,0.5,0.63], psdCumm=[0,0.0325,0.0725,0.1875,0.365,0.605,0.85,1.0],distributeMass=True, num=num_spheres, porosity=0.7,seed=1,periodic=True)
#### to replace the spheres with polyhedral particles ####
for center,radius in sp0:
 t=polyhedra_utils.polyhedra(matPoly,v=[],fixed=True,size=(radius,radius,radius) ,mask=1) ### to give the polyhedral particle the same "radius" as the sphere
 t.state.pos=center ### to give the polyhedral particle the same position as the sphere
 O.bodies.append(t)

2- I am also trying to prepare a Triaxial test with polyhedral particles and walls:
I created walls using aabbWalls
Use the TriaxialStressController engine to control the stress
However, when I run the simulation I get a warning saying that the verletDist is set to 0 because no spheres are found… and the walls move with no contacts with the particles. I used the following engine and laws:
O.engines=[
 ForceResetter(),
 InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb(),Bo1_Facet_Aabb()]),
 InteractionLoop(
  [Ig2_Facet_Polyhedra_PolyhedraGeom(),Ig2_Polyhedra_Polyhedra_PolyhedraGeom(),Ig2_Wall_Polyhedra_PolyhedraGeom()],
  [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys(),Ip2_FrictMat_PolyhedraMat_FrictPhys()],
  [Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()]
 ),
O.dt=0.25*polyhedra_utils.PWaveTimeStep()
...... I tried to change the time step but its the same.
When I run a previous example “examples_polyhedra_free-fall.py” I get the same warning. but its working.

Thank you for your help!