LIGGGHTS® - User Forum

LIGGGHTS® related topics can be discussed here: discussion about models, installation, feature requests and general discussion

particles rolling with high omega using rolling friction model epsd

Submitted by oranjesr on Sun, 04/19/2020 - 05:16

Dear LIGGGHTS users,

I recently use rolling frition model epsd for triaxial test simultation. When I create a bulk of particles, and let them settle under gravity, high omega is observed for alomst every particles, up to 500 rad/s. My particles have a radius from 1..10mm.

Can anybody help me with this.

neighbor 0.0005 bin
neigh_modify delay 0

hard_particles yes
# Material properties required for granular pair styles

Particles disappear during increase in gravity

Submitted by AzamatSalamat on Sat, 04/18/2020 - 07:47

Hello!

I am doing simulations of powder compaction using artificial gravity. I have been simulating them under 2000g. So once this gravity is applied, particles start to disappear as in the screenshot attached. You can see that number of particles reduced from 7884 to 7508 and then it stays constant. Reducing the time step from 1e-6 to 1e-7 seconds has increased the final number of particles, but only slightly, so some particles are always disappearing.

Read_data for bonded particle model

Submitted by PerB on Wed, 04/15/2020 - 11:46

Hi!

Im having trouble with read_data for a model with bonded particles (Schramm/Richti83) using adhesive particles (washino/capillary/viscous).

I cant seem to get the format right in the data-file. Is there a way to check what format my particles have (e.g. ID x y z ...) When I use write_data to create a file I get the following format:

Atoms

1 1 9.0780778811098348e-03 8.4565641784780238e-03 1.9864876941739838e-03 1.8000000000000000e-05 1.0000000000008616e+03 0 0 0

LIGGGHTS youngsModulus limitation for CGS units is not consistent with SI units

Min Zhang's picture
Submitted by Min Zhang on Tue, 04/14/2020 - 00:34

Hello all,

I am using LIGGGHTS version 3.6.0 and CGS units. I am getting this error:
ERROR: youngsModulus <= 1e8 required for CGS units

So I think the unit here is dyn.cm−2 (=0.1 Pa), which means the limitation here is 1e7Pa. Am I understanding it correctly?
If it is yes, then it is not consistent with the youngsModulus limitation in SI units (1e9Pa). So I am wondering why?

Thanks!

Tube filter is giving opposite scaling

Submitted by Sounik on Sun, 04/12/2020 - 08:23

Hi everybody,

I am showing the force chains of a granular assembly obtained from LIGGGHTS using Paraview.
Problem1: Scaling is occurring in the reverse order, i.e. the weaker force chains are thicker and the strong force chains are thinner.
Problem2: Horizontal lines can be seen joining the lateral boundaries.

PS: I am using periodic boundaries at the lateral boundaries

Any suggestions will be very helpful.

Thanks,
Sounik

Different forces from particles dump and computed from contacts dump

Submitted by iluvatar on Thu, 04/09/2020 - 23:41

Hello,
I need to compute some parameter that relates the forces on particles to the force on contacts. I assumed that if particle i have n contacts, by adding up the forces at each contact I will recover the total force on the particle (without external forces). Unfortunately, that does not seem to be the case.

How to insert initial positions for multispheres?

Submitted by Andrea Mellado on Mon, 04/06/2020 - 17:10

Hi,

I would like to know if it is possible to define the initial position of a multisphere in LIGGGHTS (eg. placing a cube in the bottom of a tank) instead of inserting it at random places in a domain?
Also, in the "multisphere_stone_restitution" example, I did not understand how the atoms' positions were defined: in the "stone1.multisphere" file, the x, y, z coordinates are all positives, but when I output the results in Paraview, the x and y coordinates are negative. How are the atoms' positions actually defined?

Thank you for your time,
Andrea M.

Pages

Subscribe to RSS - LIGGGHTS® - User Forum