LIGGGHTS® - User Forum

hello i have question regarding to the shape of the simulation box

as far as i know the simulation box can be defined with

"create_box N reg" command with previously defined geometric domain

is it possible to define a simulation box with mesh/tet region domain or combination of block and cylinder shapes?

thanks.

Dear all,

I am doing a simulation which needs a quasi-static state during my simulation, i.e a triaxial compression test.

Hi everybody,

Hi,

Hi,

Has anyone experienced problems with simulations stopping in what appears to be random points, without returning any errors? To make things more confusing, the point it stops at is not repeatable. All I am using is a modified version of the conveyor tutorial fine included with LIGGGHTS. I am currently using 3.7.0, but had a similar problem with 3.3.0. The issue with 3.3.0 seemed to appear just after upgrading to Ubuntu 16.04, but the current 3.7.0 has been recompiled after upgrading.

Hi,
I have some questions regarding the implementation of the “cohesion_model_easo_capillary_viscous.h” (or for that matter washino) which I would really appreciate your reply to help me understand the algorithm and to properly apply the code. My main questions are with regards to the criteria to evaluate the three possible cases of no-bridge, active, breakage of the liquid bridge force.
Thank you for your time and consideration.
Ali

Dear all,

Dear all,

I am a little confused about the usage of the variable type of f_a and ${a}. To give an example, when I use "fix a all mesh/serface/stress/servo ...", I want to obtain the displacement of the servo wall, thus I use "variable loct euqal f_a[9]-0.1" where 0.1 is the original of the servo stl wall. When I want to applied this variable in another command (e.g. fix b all mesh/surface/stress ... move 0. 0. ${loct}), the ERROR message shows: invalid fix ID in variable formula.

When trying to create bonds in "LIGGGHTS-WITH-BONDS" , I get the following error:

creating bond btw atoms 726 and 1741 (i has now 12 bonds) at step 1001
creating bond btw atoms 726 and 1743 (i has now 13 bonds) at step 1001
ERROR on proc 0: New bond exceeded bonds per atom in fix bond/create (../fix_bond_create_gran.cpp:485)

Dear all,

I am using command "read_restart" to restart my simulation and used the same servo wall command as the original code which is applied to create the restart file (shown below).

"fix topwall all mesh/surface/stress/servo file meshes/topwall.stl type 2 com 0. 0. 0.1 ctrlPV force axis 0. 0. -1 target_val 300 vel_max 0.03 kp 1"