LIGGGHTS® - User Forum

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

Urgent "point_at_outlet" explanation

Submitted by nicolasoviedoc on Sat, 11/28/2015 - 14:52

Greetings,

I'm writing because i can't understand how to use the point_at_outlet keyword. Specifically, i don't know how it works. What i understand is that you define a point in the "outside" part of the mesh, but what simbolize that point?. For example, my cylinder-shaped mesh is in the coordinates 4.48,0,0. I would have to set the point in 4.9,0,0 in order to the keyword works?

Another thing, the inside_out keyword, how it works too?

I hope that someone can answer my doubts.

Nicolás.

Potential Energy of the system always zero while deforming the box

Noel682's picture
Submitted by Noel682 on Thu, 11/26/2015 - 18:28

Hi

## version 3.3.0

I am running this code, which inserts two layers made of two sets of random generated packings.
Since I need to compress the volume in order to reach a relatively high volume fraction (higher than 0.8), I apply an isotropic compression to the box, reducing its final dimension.

The problem I have is that the potential energy of the system during compression and relax, is always zero, while the kinetic energy has quite reasonable values.

region mesh/tet implementation

Submitted by partnow on Sat, 11/21/2015 - 13:56

The LIGGGHTS 3.3.0 release announcement states: "Also made an important improvement to region mesh/tet. The all_in yes feature was implemented, and some algorithms necessary for particle insertions have been improved considerably. You should see a speed-up when using tet meshes for insertion."

Is this feature fully implemented ? All_in seems to be only an option for "fix insert".

Inspecting the code reveals that a mesh/tet region always enforces the caluclation of the corresponding surface of the region

Total force on some mesh elements

Submitted by rob1 on Thu, 11/19/2015 - 15:27

Hello,

I would like to compute the total force acting on only some portion of a stl mesh.
I know that dump mesh/vtk output the stress components.

What I am looking for is to obtain an output like
fix ... ave/time ... f_cad[1] f_cad[2] f_cad[3] file.dat
but not on the whole cad, only in some parts of it.

Maybe should I split the cad file in different parts and import different stl files?

Thanks,
rob1

Inertia tensior of multisphere caculated correctly?

Submitted by zhou_yu on Tue, 11/17/2015 - 11:52

The particle clump is made up by some overlapping spheres ,I have confirmed that the volume and mass can be caculated correctly in liggghts. Now I wonder IF THE INERTIA TENSOR CAN BE CACULATED CORRECTLY ? since it is concerned with shape and density distribution. It says the caculated inertia tensor is based on the mass center of the clump, does that means the coordinate origin is moved to the masscenter automaticly ,or should I give the spheres locations based on mass center ?
PS: does anyone knows how to find the masscenter of a irregular shape clump using MATLAB?

Pages

Subscribe to RSS - LIGGGHTS® - User Forum