LIGGGHTS® - User Forum

I have been wondering if it is logical to perform fix nve or fix nve/sphere on a pack of frictional grains? since in the case of friction energy can not be conserved and is lost in each collision. but all the sample input files that I have seen in the example directory of liggghts are integrated in with NVE.

Hi everybody

Im currently modifying the movingMesh example to learn how to use Liggghts
Till now I have being able to modify the region geometry, size of particles, particles insertion region, etc

But now I got stuck when I try to apply the wall/gran fix to the stl file that I defined. The program simply crashes.

Does anyone detect where is my mistake?

Here is the scrypt:

atom_style granular
boundary m m m
newton off

communicate single vel yes
units si

#region reg block -.1 0.6 -.1 0.2 -0.3 0.25 units box
region reg cylinder z 0 0 .06 -0.1 .3 units box

Hi,
Does anyone is now using the gran/hooke/history/stiffness command?
I am using it but system warned:
[localhost:13166] *** Process received signal ***
[localhost:13166] Signal: Segmentation fault (11)
[localhost:13166] Signal code: Invalid permissions (2)
[localhost:13166] Failing at address: 0x7fc06aa59000
[localhost:13166] [ 0] /lib64/libpthread.so.0(+0xf500) [0x7fc06fe13500]
[localhost:13166] [ 1] liggghts() [0x51681c]
[localhost:13166] [ 2] liggghts() [0x59761b]
[localhost:13166] [ 3] liggghts() [0x4fe93b]
[localhost:13166] [ 4] liggghts() [0x426a34]

Hi

Im able to run the simulation and in paraview observe the mesh created through stl file.
Also I changed the scalefactor from 1 and 0.1 and was able to appreciate the difference in size of the mesh.

My problem is the offset. I tried to move the mesh by changing the xoff, yoff and zoff, but nothing happens. The mesh always stays in the same place.

Have someone had the same kind of problem?

Hi everybody,

I'm looking for an effective way to generate a dense bed of spherical particles in a 2d simulation.
I hoped that fix_pour_dev_packing could be used for that - also the doc says that - but since one has to use fix_particledistribution_discrete together with fix_particletemplate_sphere which cannot be used in 2d simulation it seems to be not possible.

Is there any other way to generate dense packing in 2d simulations effectively instead of simulating all the pouring process?

Thanks for the help,
kind regards,

Sebastian

Hi guys,

I am doing the simulation of shear box test, I split the simulation into packing (based on the packing example)and shearing.

As a trial run, about 200 particles are inserted and it works just fine for the packing part. But in the shearing part,
I used the:
1. fix rigid command to create a rigid body to push the particles contained in the "upperbox and lowerbox"
2. fix move/mesh/gran to move the lower box to do the shearing.
However, just when I added force on the rigid body to push the particles, system stopped and printed:

Dear community,

weeks ago I posted some questions regarding a set-up of a mill simulation.

Now, after much help from richti83, the ball mill simulation progressed.

But one problem still remains: the balls are "glued" to the lower(!) part of the lifting bar (see attachment).

I couldn't fix this problem. I tried to mesh it via GMSH, I drew it with ProE with zero wall thickness (bars had still thickness) then exported it to STEP and finally to STL, I tried to increase the wall thickness but it did not help either - the balls are still under the lifting bar when rotating (nonsense).

I am a beginner of liggghts and need to make a simulation of oedometer test which is a one-dimensional test and the clay sample is immersed in water, the top wall and bottom wall allow the passage of water during consolidation. The followings are my questions :

1. Should there be two types of particles , one is the clay particle and the other one is water particle?
2. How can I create the wall which allow the passage of water particles but not for clay particles?
3. Could anyone send me example of the script of simulation of oedometer test?

Hi guys,

I am now still testing the packing of particles and I used the packing example.
Differently, I used the pair_style of gran/hooke/history/stiffness 1 0 and tried to used slightly smaller particle size (I tried 3mm but failed), I also import mesh from stl files generated by Gmsh.
My problem is:

1. the Pour/dev function actually insert much less particles into the region, (less than 100), but the region is big enough to contain a lot more particles.

Hello all,