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
mschramm | Mon, 03/01/2021 - 16:20
superquadric particles
This can be done with superquadric particles.
Mateus Kostka | Thu, 04/08/2021 - 16:34
Superquadric particles
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
mschramm | Fri, 04/09/2021 - 17:58
doiuble check the example simulation
Hello,
Please double check your compile options. You should see a "define" set when compiling the files (-DSUPERQUADRIC_ACTIVE_FLAG -DNONSPHERICAL_ACTIVE_FLAG).
If everything was compiled correctly, try running the example in tutorials.