LIGGGHTS® - User Forum

Hi,
I want to create 2D layer of beads. fix pour only works in 3D while lattice is not working with 2D. Any help will be appreciated.
Thanks,
Ram

i download the new version and followed the same steps to install it as i followed for the last versions. It was installed successfully and did not give any error in the end.The message in the end is as

"size ../lmp_fedora
text data bss dec hex filename
2677691 2572 10416 2690679 290e77 ../lmp_fedora
make[1]: Leaving directory `/home/sajjad/Desktop/liggghts_1p2p3-beta /src/Obj_fedora'
[sajjad@localhost src]$ "

But when i write the command " lmp_fedora" it gives

bash: lmp_fedora: command not found

Or when i execute any input file it gives

Hi,
I have imported rectangular mesh wall vertically in simulation box and now want to find total pressure experted by the particles. To do this, i have made mesh with no triangles, just one triangle of whole CAD. I don't want to measure pressure as function of height but total pressure at a time. Is this right way of doing?

Thanks,
Ram

Hello everyone,

I want to know if it is possible to use two different neighbor lists in my simulation. I use one neighbor list to load once my interactions via a txt file to define two bodies and now I want to detect the new interactions between these two bodies created during my simulation. To do this, I need to use a second neighbor list linked to a second pair_style.

Please, can you help me to do that ?

Hi, all

Could someone please tell me how to use the compute pair/gran/local and the dump local command to output pair forces in detail? I do not know how to use them. Actually, after dump the pair forces, it will be possible for me to extend the pizza.py to realize them in Paraview.

Could I add these two lines in the input file?
.....
compute 2 all pair/gran/local force
.....
dump dmpf all local 1000 dump.force c_2 #or c_2[1] c_2[2] c_2[3]
................

Hi, Christoph. I have a question about Rayleigh time:

Question 1:
in code:

rayleigh_time_i=M_PI*r[i]*sqrt(density[i]/shear_mod)/(0.1631*nu->values[type[i]-1]+0.8766)

in doc
dt_r =PI*r*sqrt(rho/Y)/(0.1631*n+0.8766),

but in doc Young Modulus instead of shear_mod, I want to known, which was right?

Question 2:

what are the Rayleigh time and Hertz time relationship with timestep of run, and how timestep should be determined.

Question 3:
Are there some references about this topic?

Thank you
xz

Hi every one,
I am a new user of LIGGGHTS. I have installed it and working on it for a couple of weeks. Now when I am running my input script, it is giving me the following error.
"error while loading shared libraries: libimf.so: cannot open shared object file: No such file or directory"
What this libimf.so shared object is? Why it is missing now while every thing was fine a few days before?
I have tried to recompile it but is giving errors
Thanks in advance.
Naved

Hi,
Why I can not dump correctly the stresses on moving mesh? How frequent should I dump? Can I dump like this frequency :
dump dump all mesh/gran/VTK 1 dump*.vtk stress cad2

Thanks in advance,
Ram

Hello everyone,

I have some troubles to use the command fix move linear.

The function fix move wiggle works fine: fix 3 top move wiggle 0.0 0.01.0 1.0 units box

However when I try with : fix 3 top move linear 0.0 0.0 1.0 or fix 3 top move linear 0.0 0.0 1.0 1.0 units box

The simulation does not work. Please, can you correct my error so as to I use this command ?

Regards
JF

Hi,
I want to change the density of already poured particles in the simulation box. Can I change it with fix adapt command or set command? Or, I should re-pour the particles (with new density) if it is not valid way of doing the simulation.
Thanks in advance,
Ram