LIGGGHTS® - User Forum

I am trying to convert a 3D simulation of granular flow that includes STL objects into some form of periodic or true 2D simulation in order to reduce run times. I'm hoping someone can point me to an example or to documentation specific to true 2D simulations or to simulations that include a single periodic boundary condition.

I had no luck with my attempts to set "dimension 2". LIGGGHTS complains that my grains are 3D objects.

Can one use LIGGGHTS for pharmaceutical coating application?

hi i cannot get the atom charge or electric charge in my dump files even after adding q to it all i get is zeros if someone could help i'd greately appreciate it thanks

dump dmp all custom/vtk 100 post/fall*.vtk id type type x y z ix iy iz vx vy vz fx fy fz q omegax omegay omegaz radius

Hi, i want to rotate a region i defined for post-processing purposes. I use:

region reg1 cylinder y 1.0 -1.0 0.35 4.75 5.25 units box rotate 1.0 0.0 0.0 0.0 0.0 1.0 0.0

However on execution i get the error "illegal region command". What is wrong, i thought i could rotate a region as described in "region command" documentation.

Thanks!

Hello all,

I am trying to model a confined cylidrical box which consists of uniform particles. The particles have to be compressed to a stress of 100N and then under the stressed state, I have to apply sinusoidal load. I am using a mesh to represent the plate at the top. I use mesh/surface/stress/servo to apply the force 100N on to the particle. I cross checked the force by looking at the force values of the mesh plate.

Hello Friends, I have found from the documentation that the following property definitions should be used in order to use hooke/stiffness model.

fix m1 all property/global kn peratomtypepair 1 1e03
fix m2 all property/global kt peratomtypepair 1 0.0
fix m3 all property/global gamman peratomtypepair 1 0.0
fix m4 all property/global gammat peratomtypepair 1 0.0
fix m5 all property/global coefficientFriction peratomtypepair 1 0.0
fix m6 all property/global coefficientRestitution peratomtypepair 1 1.0

So I'm trying to simulate a vibratory feeder using the "wiggle command" to move particles up a 3 degree incline. I used a cylindrical path to drop the particles into and its supposed to shake and move the particles upward the 3 degree incline. For some reason I was not able to model an inclined ramp so I modeled it flat and modified the gravity vector to simulate an incline. The particles climb up the ramp at a 5 to 6 degree incline but not at a 3 degree incline. Can someone look at my code and see if possible my gravity vector is incorrect or some other settings needs to be adjusted?

Hi!

Hello Friends, Greetings. We all know that particle simulations are high dissipation systems. Lets say we have a cuboid container defined by 6 primitive planes oriented along X-Y-Z axes, filled with spherical particles to a certain height, gravity acting along Z. When this system is fluidized (that is, particles are set into motion),

1. How can we calculate energy that is being lost due to (a) collisions, (b) friction?

I'm tracking instantaneous net energy of the system using compute totRotKE all erotate/sphere and compute totLnrKE all ke.

Hello Friends, Greetings.

Lets say we have a cuboid container defined by 6 primitive planes oriented along X-Y-Z axes, filled with spherical particles to a certain height, gravity acting along Z. When this system is at rest,

1. How can we calculate net downward force acting on bottom wall (as a single value)? This must be equal to weight of all particles combined.
2. Is it possible to get force field on this bottom wall? This must be uniform, vectors pointing downwards (assuming flat or constant hFill) and sum of all field must be equal to weight of all particles.