Modeling of cubic or cylindrical particles

Submitted by Mateus Kostka on Mon, 03/01/2021 - 08:28

Dear Ladies and Gentlemen,

is it possible to create cubic or cylindrical particles and to assign to them the hertzian contact model with friction?

I would like to model the movement of a block or cylinder along a defined way and to analyze the velocity during this movement.

Best regards and thank you

Mateus Kostka

Mateus Kostka | Thu, 04/08/2021 - 16:34

Hello everybody,

i tried to implement superquadric particles in that way:

Contact model:

pair_style gran model hertz tangential history rolling_friction cdt surface superquadric gaussianCurvature yes

pair_coeff * *

Particle insertion:

fix pts1 all particletemplate/superquadric 123457 atom_type 1 density constant 8742.0 shape constant 0.0007 0.0007 0.0007 #Bronze
fix pts2 all particletemplate/superquadric 123457 atom_type 2 density constant 8742.0 shape constant 0.0007 0.0007 0.0007 #Bronze

fix pdd1 all particledistribution/discrete 22751 1 pts1 1.0
fix pdd2 all particledistribution/discrete 17903 1 pts2 1.0

region factory1 cylinder x 0. 0. 0.0125 -0.00497 0.00345 units box #Should be 4122
region factory2 cylinder x 0. 0. 0.0125 0.0044 0.006885 units box #Should be 2100
fix ins1 all insert/rate/region seed 27647 distributiontemplate pdd1 &
nparticles 27647 particlerate 50000 insert_every 100000000 overlapcheck yes &
vel constant 0 0 0 region factory1 ntry_mc 10000
fix ins2 all insert/rate/region seed 18307 distributiontemplate pdd2 &
nparticles 18307 particlerate 50000 insert_every 100000000 overlapcheck yes &
vel constant 0 0 0 region factory2 ntry_mc 10000

when I use these commands nothing happend. I assume that I need the boost library and adapt the make file.

Can anyone explain, how to do this and in which folder I have to insert the library?

Maybe there are other reasons for this problem?

Apart from that I would like to know if its possible to model complex geometries close to the shape of an octahedron. How can I realize this?

Best regrads and thanks in advance

Mateus Kostka