LIGGGHTS® - Developer Forum

He everybody,

I am writing my own rigid bond contact model and therefore, I am using the cohesion model.
It all works nicely and I want to write in a file all the bonds that have been set and the bonds that are broken at the respective timestep.

In general it is quite easy to do so, but I have hard coded the filename into the contact model. I was wondering if there is a way to give a filename in the input script. I know how to provide booleans like the tangential_reduce in "gran cohesion easo/capillary/viscous".

I am grateful for any idea.

Hello everyone,

when trying to compile the newly released 3.6.0 on my 64 bit Ubuntu, I get the following error message:

In file included from ../dump_mesh.cpp:43:0:
../dump:mesh.h:52:29: fatal error: vtkSmartPointer.h: cannot find file or folder
compilation terminated.

I tried to compile using:

make clean-all
make fedora

Furthermore, I unsuccessfully tried to compile using:
make clean-all
make auto

Does anyone know what the problem could be?

Best regards,
MiRa

Hello all,

I am using C++ interfacing with LIGGGHTS by following an example case given on github.com by the name of QUEST - (lmpqst) mentioned below:

https://github.com/CFDEMproject/LAMMPS/tree/master/examples/COUPLE/lammp...

I am not able to access the following variables from the file neigh_list.h:

int *ilist;
int *numneigh;
int **firstneigh;

Hello,

I am using C++ interfacing with LIGGGHTS by following an example case given on github.com by the name of QUEST - (lmpqst) mentioned below:

https://github.com/CFDEMproject/LAMMPS/tree/master/examples/COUPLE/lammp...

Hello everyone,

when looking at the documentation of the EPSD2 rolling friction model [1] it seems there is a missing reference. In the documentation it says "see Iwashita and Oda" but no such reference can be found on this page.

Best regards,
MiRa

[1] http://www.cfdem.com/media/DEM/docu/gran_rolling_friction_epsd2.html

The current multi-sphere model assumes a rigid body (non-deformable) for each multi-sphere particle. Is it possible to extend this multi-sphere model as to include bond between sphere particles within multi-sphere body, so that the multi-sphere particle can deform upon impact?

Hi CFDEMers/Liggghtsers,

Recently, I had scanned some non particles (like corn) by 3D laser scanner, and the output format of particles is .STL or others. But how to use the scanned particles in liggghts? This quesition puzzled me for a long time, even I searched and went through the related topics in the forum.

Is it possible to use the scanned particles rather than the multi spherical particles by clumping spheres?

Any comment or clue about this is appreciated,

Cheers!
Wei-Han

Will I mess up anything in the code if I change line 91 in fix_template_multispheres.cpp to read "if(nspheres < 1) error->fix_error(FLERR,this,"illegal number of spheres");" instead of "if(nspheres < 2 error->fix_error(FLERR,this,"illegal number of spheres");"?

I just want to create an assembly of spheres but get data that is only available to multisphere or aspherical particles. That way I can use multisphere computes to get quaternions and calculate particle rotations. When I use basic spheres and no multisphere fix I can't get particle rotation angles.

Hello everyone,

I'm trying to write a new compute script to throw in the src folder that sums up the magnitude of the contact forces in the system. I tried to basically copy compute_pair_gran_local.cpp since that file directly works with the contact forces, and in my new file I would put a compute_scalar command to output a single global value that I want. The problem I run into is that compute_pair_gran_local.cpp only allows a single simulation to have one PairGran object registered using the register_compute_pair_gran_local function.