LIGGGHTS® - User Forum

I'm trying to simulate a triaxial test (compression) on a cubic sample surrounded by 6 mesh faces. Only the top face will move in (-ve) z direction to compress sample while the other 5 faces are stationary. I want to draw the stress-strain curve for this sample to get the yield and ultimate stresses. My problem is that during compression the stress is only increasing and never decrease. I think this is because the lateral deformation for the sample in the other two directions (x,y) is not allowed.

Dear Support Team, I am following the installation/compilation instructions as provided in the documents section.
At --2.2.2. Steps to build a LIGGGHTS(R)-PUBLIC executable-- when I do "make auto", I am getting the following error:

$ make auto
Creating list of contact models completed.
/usr/bin/sh: -c: line 1: unexpected EOF while looking for matching `"'
/usr/bin/sh: -c: line 2: syntax error: unexpected end of file
make: *** [Makefile:114: auto] Error 1

I am trying to change the group particles are assigned to based on their location in the simulation. I can do this as a one time thing, but I want to keep doing it throughout the simulation. Consider just a periodic box with particles randomly moving. If I want to separate particles into different groups based on y location, I can use

variable loc1 atom y
variable loc2 atom "y > 0.5"
variable loc3 atom "y < 0.5"
group top_group variable loc2
group bottom_group variable loc3

Hi everyone,

I am trying to perform a quasi-static uniaxial compression of an assembly made of a cubic box with inside 5000 particles randomly distributed with a radius of 0.00025 meters.

I am realizing this by using the "fix deform" command and a double loop. The external loop is for compression, the internal one for dissipation of kinetic energy realized with the "fix viscous" command.

I try to simulate an assembly with two group of particles. However, there is only one group of particles (ins1) added to the simulation while the other (ins2) is ignored. My code is given below. Could anyone tell where I did wrong? Thank you!

atom_style granular
atom_modify map array
boundary m m m
newton off

communicate single vel yes

units si

region reg block 0.0 2.0 0.0 0.1 0.0 1.0 units box
create_box 2 reg
group nve_group region reg

neighbor 0.001 bin
neigh_modify delay 0

Dear all,

Hello,

I am implementing heat transfer for granular media. I need to set two different temperatures: one for the boundaries(stl file) and another for the particles.
Is there any tutorial that can help or any other suggestion?

Thanks in advance,
Clara

hello everyone i am trying to apply a radial electrostatic force on particles in the documentation we can use only x y z axis but the force vectors are always pointing in a single direction i need them to be radial like applying constant stress to a cylinder and inclined if possible like using a conical electrode any tip is welcome thank you all.

Hello everyone, I need some help.

Currently, I am doing simulation in LIGGGHTS. My particles insertion area and surface have height approximately 0.02 in si units.

I would like to add some random height so that my particles drop from different height in the same simulations. Purpose of this is to bring noise on heights.

So I want to drop particles randomly at the height of 0.02 +- 10 % or 0.02 +- 20 %while running simulation.

Is it possible to add this idea in liggghts script?

What will be the command?

Thank you.

Hello
I am using "LIGGGHTS-flexible fibres"
For my study, I need to calculate energy stored in bonds such as normal, shear and bending. Is there a command to directly calculate that ?
If not, is there a way to find bond forces in normal, shear, moment, so that I can find energies by integration? As I can see the 'compute property/local ' can give bond forces in X, Y, Z not in normal and shear.

Thanks
Shubham