LIGGGHTS® - User Forum

Hello,

Hello everyone,
This post is to follow up with a problem of "isolated particles" (particles having no contacts with other particles) having such a high and unreasonable omega, especially around the base of the pile. This problem was mentioned a few years ago in this thread: https://www.cfdem.com/forums/constant-high-omega-few-particles-close-bas....
I just wonder:
1) Does anyone know how to fix or work around this problem? Any ideas are highly appreciated.

I am currently doing my research on metal 3-d printing and I was hoping if anyone can help me set up a powder bed using liggghts.
Thanks in advance.

I am currently doing my research on metal 3-d printing and I was hoping if anyone can help me set up a powder bed using liggghts.
Thanks in advance.

I'm trying to simulate particles with radius 0.0002m (200 micrometers). If I used "SI" units the required dt is about 1e-6 second.

Does using "micro" units in liggghts instead of "SI" units accelerate the run time since the time unit in "micro" unit system is microseconds not second, so dt will be equal to 1 microsecond?

Is there any differences in the output results if simulating the same problem with different unit system?

I'm generating particles in script (1) then save a restart file. Then apply some forces on the particles in script (2) after reading the restart file.

My question is can I define different inter-particle interactions in script (2) without changing it in (1) or they must be the same since script (2) reads the restart file of (1).

I want to do that to avoid running script (1) each time I need to change these parameters to save time.

Hi,

Could you please explain how the hard particles (Event-driven) approach is implemented in LIGGGHTS? I couldn't find the details in the documentary and source code.

In the event-driven approach, usually the overlap of the particles doesn't happen, and the cohesion models are dependent on the overlap. So how can we use the cohesion models in hard particles?

I request to please provide the relationship between cohesion energy and surface energy used in literature in JKR models?

Thanks and regards
David

I'm trying to make a triaxial compression test in two scripts. First, generating particles using insert/stream in a cubic region then unfix all walls and save a restart file. Second, fix new servo walls and apply equal pressure from all walls except the base wall and after reaching the target confining pressure, I only increase the pressure of the upper wall to apply vertical compression on the sample. I use kp = 1, ki = 0, kd = 1e-6.

My problem is that some steps takes too much time especially the compression steps. Here is the output timings:

I am still trying to run a simple shear flow. Since nvt/sllod was not carried over from LAMMPS, I am wondering if anyone knows any alternatives. Following LAMMPS guide for a simple shear, I can run a simulation with a deforming triclinic box with remap set to v. Normally, the nvt/sllod integration would update the velocities, positions, and control the temperature.

Hi all,

i have to questions regarding friction coeffcient. First, why is the upper limit for the friction coeffcient = 1?
There are plenty of materials that have a friction coeffcient greater 1 in literature.

Second question is a bit more special. How can i implement a model that allows for friction coeffcient's greater than 1? And also would it be possible to implement a friction model that allows to set the friction coeffcient in relation to the normal force that act's on individual particles?

Thanks a lot for hints!