LIGGGHTS® - User Forum

Hello,
i am trying to simulate a shear test. In the simulation, i need to apply a constant force on the coverwall. does anyone know how to do it?

Thanks
Ada

Hi Christoph,

When I check the package status, the Granular package was not installed and the messages appeared as follows. Please give me some advice to handle it.

Thanks.
Jianwen

username@ubuntu:~/liggghts_1p3p1/src$ make package-status
Installed YES: package ASPHERE
Installed YES: package CLASS2
Installed YES: package COLLOID
Installed NO: package DIPOLE
Installed NO: package DSMC
Installed NO: package GPU
test: 18: ../angle.cpp: unexpected operator
test: 18: ../add_fix_property.h: unexpected operator
Installed NO: package GRANULAR
Installed YES: package KSPACE

Hi Christoph,

I have noticed that the GPULAMMPS (LAMMPS Cuda) has already allowed the gran/hooke/cuda since Feb 3, 2011. And some guys still work on the Hertz contact model.

Is is possible to get this gpu-acceleration for hooke in LIGGGHTS?

Best,
Tong

Hi Christoph,

I want to ask about what rolling friction model is used in LIGGGHTS. And what is the typical value of rolling friction coefficient for sand?

Thanks in advance!

Best,
Tong

Hi,
Could anybody give us a list of published papers (or thesis) using the LIGGGHTS code ?
Thanks a lot .
Alex

Hi all,

I'm new to Ligggghts and this might seem a little 'out-of-liggghts' question but I thought you might help.

I'm running simulations on a supercomputer provided by my university.

Hey

easy question, I've forgotten how to do, so I thought posting on here might be quicker.

How do I define two different atom types?

Essentially I want to pour atoms in of one type, and define meshes to be the atom style of another.
With one atom type I have never had to previously define types as the following did the job]

region incyl cylinder z 0 0 0.035 0.04 0.13 units box
group nve_group region incyl
fix ins nve_group pour 25000 1 2 region incyl diam uniform 0.001 0.002 dens uniform 1500 1700

with fixes before applied to all for atom properties

Hi,

actually I tried to figure out in the LAMMPS manual, how I could access the kinetic energy of my granular system.
I understood, that LAMMPS is calculating the kinetic energy in a MD sence (ke = 1/2 Kb T).

I wonder if the ke in my thermo-output is calculated like 1/2mv^2 for all particles if I use the LIGGGHTS package?

Thanks in advance,

Sebastian

Hi,
I'm trying to use Kn and Kt instead Young Modulus. I wrote:

pair_style gran/hooke/history/stiffness 3 0
pair_coeff * *
...
fix m9 all property/global kn peratomtypepair 1 5000000
fix m10 all property/global kt peratomtypepair 1 5000000
fix m11 all property/global gamman peratomtypepair 1 0.7
fix m12 all property/global gammat peratomtypepair 1 0.7

But LIGGGHTS crashes without error message.
What I'm missing?

Thanks for any help

Andrea

Hello liggghts-team,

I want to simulate a role(cylinder) that moves linear along the systems X-Axis and it should rotate around its roll-axis.
Is there a way to do this?

1) I tried using 2 fix commands, but as expected only the late fix-command wins.
2) Rotate the role and moving the assembling linear around the role, is not a good alternative, producing wrong results.

Is it possible to script something like this:
even step = linear move
odd step = rotate role