LIGGGHTS® - User Forum

I have adapted the compute conts function to make contact calculations for a model that depends on the surface contact area between contacting particles. For superquadrics however, the rsq <= radsumsq condition in compute_contact_atom.cpp overestimates the contact time when the superquadric contact point is not at the maximum radius for both particles. The correct calculation of deltan (total deformation) for superquadrics is needed to make this work properly. For particle-wall contacts this is much easier to find, see :

Hello friends, greetings.

I'm looking for the rotating drum tutorial in which the wall is made of particles. By any chance, any of you have a copy of in.* file?
Or
Any such files where we have a wall of particles moving / rotating? I want to understand the setup. Eagerly waiting for your reply.

Thank you,
Jagan Mohan.

Hi
I would like to count the collision frequency between different types of particles in each time step.
i have 16 types of particles :
i did manage to count for each type, but what i need is to count for example :
Between Type 1 and Type 1, between Type 1 and Type 2 ...., between Type 16 and Type 16.
here is my code :

Hi,

I am wondering if developers are on their way to implementing SJKR forces on superquadric particles.

This would be really useful for some of my simulations.

Best regards,

Oier.

Hi
In my simulation, I am using three types of atoms (say 1,2,3) I need to define the contact stiffness between each of them (E1,E2, E3)
For that I am using the following
fix m1 all property/global youngsModulus peratomtype ${E1} ${E2} ${E3}

As far as I can understand, this command has assigned a young modulus to each of the particles. How it takes care say when atom 1 and 3 come in contact. Is there a way to assign contact modulus like that? I tried peratomtypepair but it did not work.

Thanks
Shubham

Hi
I am using LIGGGHTS - flexible fibers version
I am defining a prism region using the following command:
variable dummy2 equal 1.01*${beam_length}

region bc_blade prism ${dummy2} ${dummy2}+${blade_particle_diameter} -${blade_particle_radius} ${blade_particle_radius} 0 ${blade_length} 0 0 ${blade_particle_radius} units box

All the variable are defined tin the code, however I am encountering the following error;
Illegal region prism command

How it can be solved??

Hi
I would like to count the collision frequency between different types of particles in each time step.
i have 16 types of particles :
i did manage to count for each type, but what i need is to count for example :
Between Type 1 and Type 1, between Type 1 and Type 2 ...., between Type 16 and Type 16.
here is my code :

Hello everyone,

I was wondering about one of the implementation details of liquid bridge models, namely the definition of the minimum separation distance ratio, which is required to have values > 1 and is recommended to have a value of around 1.01 according to the documentation

I am creating rigid bodies from groups of spheres using, for example,

fix 1 all rigid molecule

From doing some simple tests with two spheres setup as a diatomic molecule, it appears that LIGGGHTS will always output the total energy - translational + rotational - for ke. I want to separate this energy into rotational and translational pieces in the thermo output. LAMMPS seems to have had ways to do this with compute temp, but that does not appear to be a thing in LIGGGHTS. I tried using erotate, but that does not seem to correlate with what I need at all.

Hello
I am currently using LIGGGGHTS-flexible fibers version for granular bonds.
In my simulation I have 2 groups (web & blade), and I want to calculate the force that web exerts on the blade. For that I tried the following command:
compute force_inter blade group/group web

However, I am encountering the error saying that "Pair style does not support group/group".
In the contact style, I am using the following:
"pair_style gran model hertz tangential history "

Can someone please help me to resolve this

Thanks