yade-users team mailing list archive

[Bilal Al Tfaily <question694335@xxxxxxxxxxxxxxxxxxxxx>]

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

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 you help!


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