LIGGGHTS® - User Forum

Hi LIGGGHTS Users!

I have run into this problem where LIGGGHTS fails to restart in some cases (I guess when the particle-tri contacts are high) and exits with the following error:


Reading restart file ...
orthogonal box = (-0.00065 -0.00072 0) to (0.00065 0.00072 0.002)
1 by 1 by 1 processor grid
2198 atoms
Resetting global state of Fix zwallbottom Style wall/gran/hertz/history from restart file info
Resetting per-atom state of Fix zwallbottom Style wall/gran/hertz/history from restart file info

Importing STL file 'WireAssembly.stl'

Hi all,

I use the version 1.2.5 and I have this warning when spheres are in contact with triangular facets:

WARNING: Dangerous build in triangle neighbor list.

Please, can you give more information about this error message ?

Moreover, please can you indicate me the function where the neighbor list between spheres and facets is built. I use the function fix_wall_gran_hertz_history

Regards
JF

Hi!

I would like to know, if "Liggghts" works on the Linux version "Open Suse"? Or does it only run on "ubunto" and "fedora"?

Christian

Hi All

Hi!

I was wondering if there is a way to create composite objects in LIGGGHTS made out of other objects. For example, building a tetrahedron out of 4 spherical atoms.

If this is not possible, are there other ways to create objects with more complex geometry that can move and interact like atoms? I would like to create these objects but still have the "smoothed" spring-like model for interaction between each other.

Stefan

Hi All

I am trying to run the LIGGGHTS granular pair style with an "old" LAMMPS pair style. To do this I am using the "hybrid/overlay" command to include the effects of both pair styles. My input script runs fine with the plain granular system, but as soon as I use the hybrid/overlay command I get a segmentation error as given below (the first line shows the last line that was reached in my input script):

fix wall all wall/gran/hertz/history 1 0 mesh/gran 4 bottom top inner outer
[dwain-dunn:03288] *** Process received signal ***
[dwain-dunn:03288] Signal: Segmentation fault (11)

Hi,
I'm creating atoms region by lattice command inside a .stl container.
This works ok but I have a "problem" when I create the .vtk: the particles are as scaled respect to my container.
In detail:
- I created a cubic box region of b = 0.3 m edge length.
- My .stl container is in the middle of box with l = b/3 m edge length.
- I defined the lattice region inside this container (you can see the file in the .txt attached file).

Hi,

We're using LIGGGHTS to study a problem in materials handling for the minerals processing industry. Particles enter into a rectangular containment vessel from several defined feedport regions near its roof - this is to simulate material flowing into the vessel through feedports from storage hoppers above (for various reasons we don't want to model the particles inside the hoppers themselves). The idea is that the vessel slowly fills up with particles until they physically block the feedports, at which point the feed rate into the vessel slows down or stops ("choke-fed" condition).

Hello Christoph, I would like to include a flag that changes its state when a particle touchs for first time the wall. This is like the:
touch = firsttouch[i];

in the pair_gran_hooke_history.cpp compute function. I want to see when the contact is lost in order to make some variable assignments. So in: void FixWallGran::post_force(int vflag), I can modify the code and make something like:

if (wallstyle != MESHGRAN_FWG){
rsq = dx*dx + dy*dy + dz*dz;
if (rsq > radius[i]*radius[i])
{
if (touch[i]){
touch[i]=0;

I encounter difficulties in running a 2D simulation. A rigid rod consisting of 200 atoms is attached to a spring. The rod should only move in the xy-plane. Therefor I put 'dimension 2d' in the input file, then setup a simulation box with a small z-dimension with periodic BCs in z and fix enforce 2d. The atoms are treated as a rigid rod via the fix rigid command with force components in x and y direction and torque around the x-axis. NVE integration should be done with fix rigid.

Although I apply a strong spring force, the rod doesn't move at all?