LIGGGHTS® - User Forum

Hello everyone,
I just want to know more about the possibility of Lees-Edwards boundary condition implementation. Actually, this BC is a mixture of periodic boundary condition along with the simple shear flow. So how it can be done in liggghts?
Thank you, in advance for your help.

Hello,

I am working on discharging particles from a flat bottom hopper with an orifice. But even I keep particle frictions low one layer of particles seems to stick to the wall surface of the hopper. Even if I reduce time step and increase Young's modulus the result is the same.

Can someone please help me about that?

Here is my input script:

atom_style granular
boundary f f f
newton off
atom_modify map array
communicate single vel yes
units si #(meters, seconds, kg)
timestep 0.00001
processors * * *

Hi,

I'd like to modelling the cement bonds by using parallel bonds which consist of a finite-sized beam between the two particles. I do not know the structure of LIGGGHTS system, so here is the light questions.

#1. Is there any method that I can use the codes provided by someone else in github or somewhere? and non-specialist can do that? I think it is a quite difficult work because I have to change from the source code.

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?