LIGGGHTS® - User Forum

Hi all,

I want to delete a certain fraction of particles in a region. I use the command "delete_atoms porosity" to realize this function.

From the documentation, I notice that 10% of particles inside a given region will be deleted if I use "delete_atoms porosity cube 0.1 123457". I am just not sure how to determine whether a particle belongs to the region? Based on the centre of the particle or based on othere criteria? Thanks.

Best,
Rachel

Hello,
I am running some tests using the python wrapper and would like to hide ALL
screen output from liggghts.

Is there any way to do this?

I am using the flags pre no and post no on my runs to hide some of the output
and echo none to hide line reads but I still have some output.
I also do not want to redirect the screen output to a file as the file would grow
rather large.

Thank you.

Hi,

To my understanding, Meff (effective mass) used in contact models can be calculated from "density" defined in

fix pts1 all particletemplate/sphere ...

but I have questions as follows:

1. I need to know how to define the mass for walls including mesh to be used in Meff.

2. please let me know If there is a way to directly define Meff used in hertz/stiffness or hooke/stiffness (I want to set it unity for simplicity).

Thank you

I am running simulations of a ring shear tester and I use the fix mesh/surface/stress/servo command.

I am not completely sure about the way to calculate the tangential shear stress from the values of forces and torques I obtain from LIGGGHTS. As result, I have a global vector with 9 components for access by various output commands. The first 3 components are equal to the total force on the mesh, the next 3 components store the total torque on the mesh.

Hello
while doing post processing, I wrote the commond : lpp dump*.chute
but it showed the following message
Command 'lpp' not found

Hey guys,

I use the fix bond/create/gran:
fix bondcr11 all bond/create/gran 10 1 1 0.5 1 100 prob 1 123457

The simulation works for a probability=1 but when I reduce the probability I get the segmentation fault.
*** Process received signal ***
Signal: Segmentation fault (11)
Signal code: Address not mapped (1)
Failing at address: 0xb8

I looked for this error but I don't know how to fix it. Does anybody have an idea?

Thanks for considering my request.

Hi all,

I am doing a triaxial test using LIGGGHTS. After I consolidated my sample in a cube to a constant confining pressure 50kPa, I want to set one boundary which can be penetrated by the small particles while keeping the big particles insides. During this process, the small particles will penetrate through the boundary suddenly and the cube sample will suddenly deform (the volumetric strain could be even 7%).

Hello everybody, I'm writing to know if someone has an example of an eccentric funnel flow simulation. My idea is to run a simulation and use the output to create an animation. I could fulfill a cylindrical flat-bottom silo and discharge it through an eccentric outlet, but the result resembles much more a mass flow.

So, I guess the physical properties are not corrected.

Any help are more the appreciated.

Greetings

Hi, all,
I'm trying to analyze powder spreading process of additive manufacturing.
In my case, a cylindrical roller rolls on a bed of deposited powder and spread it to make a new layer on the fusing bed like shown simplified in the attached figure.
My question is how to rotate and translate the roller at the same time.
I have not come up with an idea to move the center of rotation.
Thank you in advance.

Akio

WARNING: time-step is 28.194306 % of rayleigh time (../fix_check_timestep_gran.cpp:200)
WARNING: time-step is 14.155515 % of hertz time (../fix_check_timestep_gran.cpp:208)

I can't help me solving this problem with help of other discussions on forum related to this concern. Can anyone kindly help me with this problem in simpler way..?