LIGGGHTS® - User Forum

LIGGGHTS® related topics can be discussed here: discussion about models, installation, feature requests and general discussion

Contact model - hertz vs hertz/stiffness -

Submitted by ykusaka on Mon, 08/06/2018 - 07:37

I'm new to LIGGGHTS and enjoying DEM simulation. I have the following questions.

1. Can I use a cohesion model such as sjkr2 with hertz/stiffness? namely,

fix boxwall all wall/gran model hertz/stiffness tangential history cohesion sjkr2 primitive type 1 xplane 1

is possible?

2. If I use stiffness model with kn, kt, gamman, gammat specified, the values of youngsModulus, poissonsRatio and damping coefficient defined by

fix d1 damp viscous ...

Different mesh rotations between LIGGGHTS versions

40044600's picture
Submitted by 40044600 on Fri, 08/03/2018 - 15:07

Hello,

I'm running simulation on both LIGGGHTS 3.3.1 and LIGGGHTS 3.8.
When I run a simulation which includes a mesh rotation about 2 axis I get a different rotations.

The command's I'm using are:

fix movemesh1 all move/mesh mesh cad1 rotate origin 0.0 0.0 0.0 axis 0. 0. 1. period 8
fix movemesh2 all move/mesh mesh cad1 rotate origin 0.0 0.0 0.0 axis 0. 1. 0. period 2

Using two "fix move" commands

Submitted by Sounik on Thu, 08/02/2018 - 00:17

Hi Liggghts Forum,

I have a compressed box of particles. The top layer of the particles needs to be moved horizontally in one direction(+y), while the bottom layer in the opposite direction(-y).

I created two groups (edge & edge1) for the particles in the top and the bottom layer respectively. My command for moving the particles are:

fix edge all move linear 0.0 0.08 0.0 units box
fix edge1 all move linear 0.0 -0.08 0.0 units box

Different VTK stress dump file formats

40044600's picture
Submitted by 40044600 on Tue, 07/31/2018 - 17:47

Hello,

I am currently running simulations on 2 computers and have found a clash in the output of the same the dump command.

The dump command i'm using is "dump dumpstress all mesh/gran/VTK 1000 dump/dump*.vtk stress cad1"
However, the vtk dump files are formatted differently for each computer I use (and this is affecting my post processing work).

Periodic boundary in LIGGGHTS and LAMMPS: differences in simulation setup?

Submitted by rangr1 on Thu, 07/26/2018 - 12:55

I am trying to setup a simple one particle simulation with periodic boundary conditions (PBC) in LIGGGHTS and LAMMPS. I am finding that the atoms are lost in LIGGGHTS-PUBLIC (v3.8.0 ), whereas the same simulation except for a minor change in the communicate command seems to work fine in LAMMPS (17 Nov 2016). Please let me know if I have to include additional commands when I try to set up PBC in LIGGGHTS. I have attached the minimal examples for your reference.

Question about density unit in cgs

Submitted by qq935074410 on Wed, 07/25/2018 - 18:36

Hi,

When I was using the si unit system, the simulation can run without error. However, when I want to change to cgs unit system and convert the density constant from 2500 to 2.5, there is an error saying "Mesh elements have been lost / left the domain". When I set the density constant to 2500 the error is fixed.

So, I am wondering what is the unit for density constant in cgs system?

Here is my input script.
# Moving mesh example

atom_style granular
atom_modify map array
boundary f f f
newton off

Particle Wall Stress Calculations - Clarification on Methods in LIGGGHTS

Submitted by dwhite on Tue, 07/24/2018 - 02:47

Hi,

I was wondering if someone can clarify for me how the stress between a particle and the wall is calculated? I am outputting the normal stress after setting the mesh to track the stress, and I ultimately want to use the stress to calculate the total force on the mesh. I specify the mesh as follows:

fix m1 all mesh/surface/stress file meshes/filename type 1 stress on

And then I dump the resulting force using this command:

dump dumpstress1 all mesh/gran/VTK 100 post/filename*.vtk output face area stress m1

Particles penetrate through the mesh

Submitted by qq935074410 on Mon, 07/23/2018 - 23:46

Hi,

I am new to the LIGGGHTS and I just want to simulate the powder spreading process. I have the STL file of a roller, and a surface to store the powder. But when I was doing the simulation, the particles always pass through the wall I imported from an STL file. Can anyone help me with this?

Thanks,
Mill

The input file:
# Moving mesh example

atom_style granular
atom_modify map array
boundary f f f
newton off

communicate single vel yes
units cgs

Pages

Subscribe to RSS - LIGGGHTS® - User Forum