LIGGGHTS® - User Forum

With LIGGGHTS 3.8.0 the script below returns the error

"ERROR on proc 0: Random number generation: It is required that the random seed value is > 10000 and a prime number.
The random seed used was 38914"

What's happening? The random seed provided is not 38914, that's the number of particles! Is it a bug?

Hi there,
I was trying to run one of the LIGGGHTS examples: LIGGGHTS/Tutorials_public/ParScale/heatTransferBed_paScale/.
I installed the ParScale following the instructions and no error appears .
When i trying to run LIGGGHTS/Tutorials_public/ParScale/heatTransferBed_paScale/Allrun_convective, the following error apppears:

Total # of neighbors = 73
Ave neighs/atom = 0.09125
Neighbor list builds = 1
Dangerous builds = 0
ERROR: Invalid fix style: "couple/pascal" (../modify.cpp:850)
data2Plot = heat.json

Hi,

Hi all,

Just wondering if there is a way to directly specify values for normal and tangential elastic constants and viscous damping constants in an input script. If there is could someone direct me to the documentation where this is explained or an example script showing how someone has done this would be awesome.

Regards,

Alex.

Hi all,

I am a little confused about one issue: In liggghts, we can output the stress of the an atom using command "compute stress/atom". If we introduce the seepage force into the simulation using CFDEM engine, e.g. drag force, buoyant force. Does the calculated stress of atom include these seepage force, or it only includes the particle-particle and particle-wall interaction force?

Thank you very much.
Zheng

Hello all,

I would like to ask whether LIGGGHTS can be used for modelling continuous materials ? For example, modelling tensile test of aluminium specimen, then compare the stress-strain diagram with experimental data. Has anybody had experience with it ?

Thank you.
Vinh

Hi all,

I want to delete all particles which have the stress larger than 10kPa, does anyone know which commands should I use?

Thanks very much.

Best,
Rachel

Dear all,

Good morning.

In my case, a mesh (plane, generated by gmesh) wall push the particle bed. However, I find the inconsistency of the total contact force on the wall by using two methods:

A) 1) fix wallleft all mesh/surface/stress; 2) variable fwlx equal f_wallleft[1]; 3) export it into the csv file
B) 1) compute cpw all wall/gran/local id pos force; 2) dump all components of cpw; 3) sum all forces in excel

The force by mesh/surface/stress is much different with the one by compute/wall/gran/local.

Dear colleagues,

I am facing a problems about deleting particles who have the interaction force with a stress wall (*.stl file) larger than 10kPa. My questions are:

1. How to monitor whether a particle have interaction with a stl wall and how to obtain the interaciton force or stress?

2. How to delete the particle when it has a particle/wall interaction stress larger than 10kPa?

Could anyone please give me any clues?

Best,
Zheng

Hi all,

I want to delete a certain fraction of particles in a region. I use the command "delete_atoms porosity" to realize this function.

From the documentation, I notice that 10% of particles inside a given region will be deleted if I use "delete_atoms porosity cube 0.1 123457". I am just not sure how to determine whether a particle belongs to the region? Based on the centre of the particle or based on othere criteria? Thanks.

Best,
Rachel