LIGGGHTS® - User Forum

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

Error :Could not locate a fix/property storing value(s) for characteristicVelocity as requested by model hooke

Submitted by Rohith tvn on Tue, 04/09/2019 - 13:40

Im doing tutorial and got an error and not able to solve. here is my code can anyone help me with it.
units si
atom_style sphere
atom_modify map array
boundary f f f
newton off
communicate single vel yes
#processors 2 2 3

# Declare domain
region domain block -0.138 0.138 -0.138 0.138 -0.0045 0.43 units box
create_box 2 domain

# Neighbor listing
neighbor 0.003 bin
neigh_modify every 1 check no

multisphere partilce - limit motion

Submitted by Neiht on Fri, 04/05/2019 - 15:15

Hi all,
Im currently using LIGGGHTS to simulate the multisphere particle system. However, there are some difficulties in sample generation because of high Young modulus.
My question is: are there any functions similar to nve/limit to restrict the motion of particles within an expected value?
Thank you!
Cheers,
Neiht

On computing the non-contacting force between two interacting particles

Submitted by XIA Huihuang on Thu, 04/04/2019 - 16:39

Dear all,

I just implemented the Van Der Waals force model on LIGGGHTS. In order to validate my numerical model, I need to know the non-contacting van der Waals force between two interacting particles in LIGGGHTS.

Actually, I tried to compute the aforementioned force using 'compute 1 all pair/gran/local id force' in LIGGGHTS, however, the results only contain some contacting forces.

I would like to know whether LIGGGHTS can output or compute this type of non-contacting force now? If not, I need to write a 'fix' function to write this kind of forces?

WARNING: Fix mesh: Mesh contains highly skewed element, moving mesh (if used) will not parallelize well

Submitted by lcraul on Tue, 04/02/2019 - 04:21

Hi, all,
I'm running a simulation with some problems on the mesh. I use a STL file as the wall of my model. I got this STL from the CFD ICEM software. However, the STL exported from this software seems to be bad. The LIGGGHTS gives neither warning or error on this STL file.

The warning is like the following line. It just warns, but the script can run without a hitch. Does it matter? Does this warning have impact on the result?

WARNING: Fix mesh: Mesh contains highly skewed element, moving mesh (if used) will not parallelize well (../surface_mesh_I.h:657)

Question about fix massflow/mesh command

Submitted by lcraul on Tue, 04/02/2019 - 03:59

Hello, all,
I'm trying to use the fix massflow/mesh command to count the particle flowing trough the outlet of my model. I use the following code in my script.
# Mass flow measurement
fix massflow all mesh/surface file Outlet.stl type 3
fix mass all massflow/mesh mesh massflow vec_side 1.0 0.0 0.0 count once file ../DEM/post/Particlemassflow*.txt

Calculating Contact Force

Submitted by fabregas04 on Sun, 03/31/2019 - 23:04

Hi all,

I already read the documentation about the contact model but I still can't get it.

(https://www.cfdem.com/media/DEM/docu/gran_model_hertz.html)

I understand the equation about the contact force but how can we calculate the normal overlap, tangential overlap,

normal relative vel, tangential relative vel in the LIGGGHTS?

I searched the code 'contact_model_hertz.h' for detail but can't understand how those are calculated.

Can anyone explain this for me?

How to manage Stress/Force using Fix deform

furkhat's picture
Submitted by furkhat on Thu, 03/28/2019 - 13:21

Dear community!

I'm trying to conduct first part the Triaxial test -confining. So , For this I used Fix deform command to compress sample. I used Deform because walls should be periodic. And now I don't know how to manage compression to achieve constant Normal force or stress.

Does anybody know what kind command should I use ?

Count Wall/Sphere contacts

Submitted by johnr on Thu, 03/28/2019 - 08:58

I would like to calculate the number of particles that hit a wall. I have implemented the 'compute conts all contact/atom' with the fix that allows the collisions to be summed (https://www.cfdem.com/forums/counting-particle-particle-collisions), but this is limited to sphere/sphere contacts. Is there a simple way to count wall-sphere contacts?

Pages

Subscribe to RSS - LIGGGHTS® - User Forum