LIGGGHTS® - User Forum

Hi, all

The input files for liggghts of an undrained triaxial test have been attached. Someone has sent me e-mail for this request. Actually, there is no inconvenience for me and I am pleased to share this.

I am a beginner and this example is not that good(you could see its q-p result in the ods file attached). This may be due to the velocity of compression is too fast or I did not use general equation to calculate q and p. If there are any problems you find or any improvement you give, please tell me cause I also need your help. Thanks for that.

Hi all,
Is there any way to set a parameter in LIGGGHTS to account for rolling resistance?
Thank you in advance,
Regards
Carlos

Hi, everyone, i'm new here. I used LIGGGHTS to simulate pouring problem. The balls are poured into a cylinder under gravity. After the calculation, i got the pressure on the bottom of wall. I would like to know this relationship between this pressure and the weight of balls. Based on soil mechanics, i know the pressure on the bottom wall should be equal the the weight of balls divided by the bottom area (Here i need to mention i am not familar with particle mechanics). Anyone can help me?

Hi
There is error when new beta version starts to write restart file in the end. I tired with both beta versions and also changed system. Error is:

Total # of neighbors = 216699
Ave neighs/atom = 5.28534
Neighbor list builds = 5595
Dangerous builds = 0
System init for write_restart ...
*** glibc detected *** lmp_fedora: free(): invalid pointer: 0xb0de6008 ***
*** glibc detected *** lmp_fedora: free(): invalid pointer: 0xb0e58008 ***
*** glibc detected *** lmp_fedora: double free or corruption (out): 0xb0dc8008 ***
======= Backtrace: =========
======= Backtrace: =========

Hi Christoph,
I am dumping force on base (mesh CAD) due to granular packing in cylinder, Base is moving downward. Practically the weight ( hence force along z-axis) felt by base should be constant during course of time but is it changing linearly. I am attaching the file herewith. Commands in the script are given as under:

...
fix cad2 all mesh/gran/stressanalysis base.stl 1 0.01 5. 5. 10 0. 0. 0.
...
fix movecad all move/mesh/gran linear 0. 0.0 -0.00082 cad2
....
dump dumpstl2 all mesh/gran/VTK 50 post/cad2/dump*.vtk stress cad2

Hi Christoph,
I have problem with new versions, i poured particles with version 1.1.9 and now I can run continue my simulation with new versions. Today i deleted the old version so I need LIGGGHTS 1.1.9. The following is error message when i run the restart file with 1.1.10 and 1.2 Beta:

......
WARNING: Restart file version does not match LAMMPS version
restart file = LIGGGHTS 1.1.9 based on lammps-10Mar10, LAMMPS = LIGGGHTS 1.1.10 based on lammps-10Mar10
orthogonal box = (-0.01 -0.01 -0.01) to (0.41 0.41 1.31)
1 by 1 by 4 processor grid
425350 atoms

Hi, all

How can I make that there's no friction between walls and particles while there's friction among particles?

Could someone tell me this in detail?

Thanks in advance.

Thanks!

Hi all,

I need to change a parameter for my walls. E.g. I want to change the damping flag (or the velocity amplitude if I am doing wiggle with a LAMMPS-style plane wall).

My impression is that I need to use the command fix_modify in my input script. But then I need to add the function modify_param in my wall class, which will take as inputs the arguments of fix_modify and reset the values I require.

Does this sound fine, or am I missing an easier way of doing this?

Thanks,

George

Hi all,
How can I use LJ units for density and yougsModulus?
Thanks in advance,
Ram

Hi Christoph,

I would want to write a new dump class (dump_gran) which saves the connectivity and the scalars as: normal force, damage parameter, radius. Thus, it will be possible to display with paraview the lattice beams or the chain force for granular simulations. But I need help to understand the class Dump_mesh and know where the functions are called: