LIGGGHTS® - User Forum

So i'm still quite new to liggghts and am having trouble with getting values for contact area and overlap. I found the compute command and realised i need to dump the data somewhere. so far im using the commands:

compute oLap_cArea all pair/gran/local delta contactArea

dump dmp2 all local/gran/vtk ${dumpstep} calculations/calc.pvtu oLap_cArea

the error im receiving is: ERROR: dump local/gran/vtk no longer allow parallel writing by setting the % character. Instead use a filename with suffix .pvtX (X = {u, p}). (../dump_local_gran_vtk.cpp:109)

Hi everyone

I am trying to install LIGGGHTS in Centos but I am unable to compile it.
------------------------ERROR--------------------
Makefile:815: *** 'Could not obtain VTK_MAJOR_VERSION. As AUTOINSTALL_VTK is not set to "ON" we cannot continue'. Stop.
make[1]: Leaving directory `/home/software/LIGGGHTS/LIGGGHTS-PUBLIC/src/Obj_auto'
make: *** [auto] Error 2
--------------------------------------------------------

I tried to install VTK manually by following steps but it did NOT work

Hello LIGGGHTS forum,
Really thank you for your help.

I have one more problem in my simulation.
I am trying to increase the radius of my particle from 0.01 to 0.2 or 0.3 but I encountered following error.

"""
ERROR: Unable to calculate region volume. Possible sources of error:

I would like to read in initial positions, velocities, etc. for multispheres rather than using insert/pack to get the multisphere particles into my simulation. It is easy enough to read data for single spheres or a rigid body following the tips in the LIGGGHTS documentation. However, there is no mention of how to do this for multispheres. It seems as though it should be possible since multisphere restarts are possible from binary files.

Hi,
I have been trying to apply a stress with the help of if-statement i.e., when the wall reaches a certain amount of force it should stop applying force.
I am not getting a proper way to implement this logic. Please let me know how can I extract dumped force and put it into if-statement.

Thanks in advance.

Hello i'm using a servo wall for a ring shear cell. Now i want to calculate the forces acting on the servo wall during a shear run. I have not been able to output the forces data. Can somebody help me? How can i output the data over time into a single .txt file. Thanks a lot!

Hello all,

We (my tutor and I) would like to study faults with DEM. We are interested in applying time- and temperature dependent contact laws. We understand we need to create our own pair_style. This one will be as the Hertz law. But friction coefficient (and maybe some other parameters) depends on time and temperature.

1. How is it possible to put the temperature of a particle into the new contact law ?

Is it possible to have different atom types for individual spheres in a multisphere?

Hello,

till now i have been using a pre compiled windows version of liggghts(not made by me) and wanted to do the switch over to ubuntu to hopefully speed up my simulations, since that version just uses 25% of my Cpu.
Now when i use the make auto command i get the following message:

Makefile:456: *** 'Could not compile a simple MPI example. Test was done with MPI_INC="" and MPICXX="mpicxx"'

I tried setting the MPICXX_USR to my mpi folder and tried setting CXX_USR=icpc, nothing changed.

Dear Ladies and Gentlemen,

is it possible to assign in Liggghts time dependent temperature profiles to the particles?

My plan ist to calculate the thermal expansion of my particles in dependence ofhe thermal expansion coefficent and the timedependent temperature.
So far I assign the temperature in that way:

fix 1 all heat/g
ran initial_temperature 373.15 contact_area overlap

Is it possible to print the temeprature of each particle in a text file?
Maybe in that way?