LIGGGHTS® - User Forum

How can I extract information on the contact model properties in my simulation (e.g., Young's modulus, coefficient of friction, etc.)? For debugging purposes, I want to be able to check that the right contact properties are being applied to different groups of particles and that the coefficient of friction changes when I put in a new 'fix property/global' partway through the simulation. I already know how to extract forces, overlaps, etc. using compute pair/gran/local but this doesn't give me any information on the actual contact model parameters.

Hello,
I am simulating the filling of a box with a large number of particles (10^5 to 10^6): first I insert the particles in a bigger volume and then I let them pack under the action of gravity.

Which is the best choice of processors to parallelize it? I would like that the Z axis is divided into non-constant grid sizes... e.g., so that more processors compute the particles at the bottom of the box and only 1-2 at the top (when the particles will stay for only a short time before packing).

Do you have any suggestion?
Thanks a lot,

rob1

Hallo,

I have a question concerning superquadric particles:

http://link.springer.com/article/10.1007/s40571-016-0131-6

Will this code be published in the public version?

Best wishes

Johannes

Hi everyone

I've got a problem by using the command massflow/mesh.

I want to simulate a micro-sandglass and would like to count the particelflow that passes the MassFlow.stl mesh (see Example_1).

I also made a simple simulation with only four particles (Example_2) trying to get more knowledge about this mysterious file called massflow.dat, whitch is written in this line:

Hello,
I am trying to increment a user variable every thermo output.
Some background, I am trying to get a soil bed to a certain void ratio. I am doing this by stepping 5000 and then checking what my void ratio is. If it is my tolerance, I continue the simulation. If not, I step another 5000 and check again.
What I would like to have though is a continuous simulation time in my thermo output. I was doing this using

variable sim_time equal dt*elapsed

but every time I would reuse the "run 5000" command, elapsed. I then tried the following

Hi,
I am trying to run a triaxial test on a cube shaped sample. I currently inset stl files (3 files for the 3 faces) using the fix_mesh_surface_stress_servo command and another 3 stl files using the fix_mesh_surface_stress command. Once I have finished the consolidation stage of the test, I want to change the top mesh (currently servo controlled) to a non-servo wall so that I can apply a constant velocity to it. Is there a way of doing this?
Many thanks,
Rector

Hi,

I am wondering how could I get the previous version (3.5.0)? I ran a simulation with 3.6.0 using a restart file generated from 3.5.0 but get different results. Anything I could after git pull the repository?

Many thanks

Does liggghts include Burger's model? If not, anyone has any suggestion how to incorporate the Burger's contact model? Thanks in advance!

Hi everyone,

Greetings! I'm a student and a new user to Liggghts. My current work is to build a model which is a plate moving horizontally in a pool full of granular particles. Now my plate can move very well, but I don't know how to extract the total reaction force on the plate, can anyone give me a hint please?
Any help is very much appreciated!

Warm regards,
Carol

Hallo,

I have created one layer of particle by using the following statements: