Question #694335 on Yade changed:
https://answers.launchpad.net/yade/+question/694335
Description changed to:
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!
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
