LIGGGHTS® - User Forum

Hi,

I want to simulate the material behaviour on a vibrating tray. I search for a feature which allows me to insert the particles which leave one side of the simulation on the opposite site. For insrance if the material on the tray moves from left to right and reaches the right end of the simulation the particles are moved to the left side (e.g. if xpos > 100 then xpos = xpos - 100).

Is there anything like this available or do I need to modify code to reach simulate this? I do not have experience in modifying or add commands.

Thanks in advance,

Urs

Hi,

I'd like to compress a box includes particles in all 6 directions. For doing that I have to build six (6) 2-D plate (stl file) and use "mesh/surface/stress/servo" for every plate. Is it true?
Or, I can build a 3-D box and compress in all 6 directions. If it is possible, would you please inform me how can I do that.

Best Regards,
Rasoul

Is there a LIGGGHTS capability to take a particle under given conditions and split it / replace it with several smaller particles? I'm specifically looking for a solution within the "granular" atom-style world (but it could be a hybrid style).

I could have sworn I saw an example of this a few months ago, but I can't find it now. Was this the multisphere capability? (node/20 illustrates multiple particles bonded together, but I want one particle that splits into multiple independent particles).

Thanks!

I have a question regarding to the usage of a realistic Young's modulus in my simulations.
Up to now I have used a value of 5e6 N/mm² for my particles because this value is used in most of the examples. Now I want to implement a realistic value of 2.1e11 N/mm² (Steel). I know that I have to reduce the timestep, but I don't know exactly which correlation I have to consider.
In the documentation of the command "fix check/timestep/gran" there are two equations and with my parameters I get the following values for my dt_max:

I am trying to play with fcc, bcc, hcp and radom style packing. All are working fine except hcp. For all other packing the unit cell is cubic and the unit vectors in all cases possess same direction. However, it is different in case of hcp.

Hi,

I implemented mesh walls (stl file), but I faced an error. This is about unit system.
"ERROR: Mesh error: dimensions too small - use different unit system (multi_node_mesh_parallel_I.h536"

If anyone knows its reason, I am looking forward to know that.

Regards,
Rasoul

Hello all,

I cannot get the servo wall to move in the -Z direction, even if I use a negative target_val. This is how I issue the command:

variable lid_force equal -1400
fix lid0 all mesh/surface/stress/servo file LIGGGHTS_files/STL_files/in/p_lid.stl type 2 move 0 0 ${lid_height} com 0. 0. ${lid_height} ctrlPV force dim z target_val ${lid_force} vel_max 0.01

When I print the z coordinates of the wall, f_lid0[9], they are strictly increasing over time no matter what values I input for target_val.

Hello all,

I am a beginner in LIGGGHTS, I am trying to implement something similar to or the same as in the http://www.youtube.com/watch?v=muZuukjb_6k video. Did anyone work on similar use case using LIGGGHTS? I am now actually not sure if this can be implemented with LIGGGHTS. Any ideas/suggestions?

Thanks,
Pradeep

Hi everybody,

In the new version of liggghts I can not create a cylindrical wall with fix wall/gran/hertz/history. I am using this code to specify the cylwall:

fix cylwalls all wall/gran/hertz/history primitive type 1 zcylinder 0.1 0. 0.7

Is there any error in this line?

Regards,

David

Hi, all!
Today, I just saved a point data of my project from "Paraview". And i found that in my "output.csv" had something strange.
Normally, In my understanding one ID(of particle) can be assigned to only one particle so that means when you save point data from Paraview inside your output.csv file for each line shouldn't have the same ID i mean no any ID is duplicated in output.csv file (if it's duplicated that mean that the particle can be located on multiple places at one time and that's impossible)
What are you thinking about this?