LIGGGHTS® - User Forum

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

How to read data from a lookup table

Submitted by khalifa on Sat, 05/04/2019 - 00:52

Hello everyone.

I am trying to implement the standard JKR cohesion model into LIGGGHTS. I modified the equations as required in normal_force_hertz.h. For speed and efficiency reasons I created an external file (a lookup table), which provides the value of the contact radius based on the overlap distance (deltan). So with the available "deltan" value, I need to access this lookup table and interpolate the corresponding value of the contact radius. As mentioned before, I need this interpolated value of the contact radius in order to ease the calculations in normal_force_hertz.h.

confusion about the definition of key word "seed"

Submitted by 1026727265@qq.com on Mon, 04/29/2019 - 12:28

I've used a DEM software, i try to learn the "liggghts", but I 'm still confused about the key word "seed" in fix command , which could not be found a corresponding value or definition in the DEM software. So is there a clear definition or meaning about this key word? Hope helps from you!

Add constant force

Submitted by Rahul Kumar on Mon, 04/29/2019 - 11:15

Hello
I am doing wall friction test simulation. The procedure is as follow
1. Fill a cylinder with particle
2. Add a constant force of 100 N and move the cylinder against wall with constant velocity
3. While cylinder is moving against the wall the force is constant
4 I want a graph between normal force and shear force on the wall

Can any body help me in this?

Atom creation with different id

Submitted by sourav0921 on Fri, 04/26/2019 - 08:16

Hi all,
I want to create two particles with two different id using the following command
create_atoms 1 single 0.95 0.4 0.5
set atom 1 type 1 diameter 0.02 density 2500 vx 2 vy 0 vz 0
create_atoms 2 single 0.99 0.4 0.5
set atom 1 type 2 diameter 0.02 density 2500 vx 0 vy 0 vz 0

but when I trying to compile the script file it says
ERROR: Invalid atom type in create_atoms command (../create_atoms.cpp:102)

Can anyone please enlighten me that waht I am doing wrong?

Thanks and regards,
Sourav

insert/stream low particlerate

Submitted by johnr on Fri, 04/19/2019 - 15:41

I have been using fix insert/stream to add particles into a simulation through an insertion face. However, the specified particle rate is not obtained if I set a high value of particle rate and a low value for the inflow velocity. I presume the absolute limit of this would be when the spheres are tightly packed, but when I visualise the inlet box using the correct particle size in paraview (glyph), the particles are well spaced out. Can anyone explain what is happening and what is the limiting criteria for correct particle rates to be achieved?

Many thanks!

faults in adapt 'scale yes'

Submitted by tjleps on Thu, 04/18/2019 - 06:13

I'm trying to initialize a simulation with a particle size distribution, created using particledistribution/discrete -> insert/pack, then grown using adapt to increase the volume fraction to 40%. If I try to run the simulation multithreaded, I get MPI faults, and when I run it single threaded I get a segmentation fault. Previously, when using a single particle size, I used the 'scale no' option and had no problems. Has anyone else had trouble with fix adapt 'scale yes'?

Pages

Subscribe to RSS - LIGGGHTS® - User Forum