LIGGGHTS® - User Forum

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

ERROR: Invalid compute style (../modify.cpp:982)

Submitted by amingheibi on Mon, 10/17/2016 - 18:22

Dear all.
I am new to LIGGGHTS and I want to run a code for a set of granular materials under compaction.
in compute section I have used the ''compute 7 all voronoi/atom'' and when I run the code I receive the error: invalid compute style (../modify.cpp:982).
I know that the error is for the "compute 7 all voonoi/atom". I think I have used the correct syntax for VORONOI compute, but still I get the error. Does any body know how can I solve this problem?
Compute section of my code is as follows:

Problem with mesh/tet region - custom vtk file

Submitted by zealakis on Sat, 10/15/2016 - 23:00

Hi!,

I am inserting atoms in a specific region using the region mesh / tet command. Now i insert the particles in two different places defined in the same VTK file. When I define two or more separate locations from each other to insert atoms, I must reduce the quality of the mesh to avoid problems like 'you have 6 neighs, but 5 expected'. If I want to generate several mini-regions, I have to make a very poor meshing and the figures lose their original form. How i could solve this problem ?, I can't use "heal auto_remove_duplicates" because it only works in 'fix mesh/surface'

Gravity for settling does not work

Claudio Wolfer's picture
Submitted by Claudio Wolfer on Wed, 10/12/2016 - 09:09

Dear particle friends

Based on a LIGGGHTS-Tutorial I have set up a planar sandglass. As first step I will fill in some particles and let them settle. A simple case I guess. But the gravity is totally inactive. I checked the input script several times tried some changes, checked the help, the forum an so on. It's no use.

Is there anybody who can quick check my script and check if it runs properly?

Thanks a lot for help.
Cilo

Unbalanced Forces

Submitted by estefan31 on Tue, 10/11/2016 - 06:01

Does LIGGGHTS have a compute style to calculate the total unbalanced force in the system? Or better yet, a ratio of the total unbalanced force to the sum of all the contact forces? I use such a ratio in other programs to determine whether my simulation has reached a stable state, and I'd like to move onto using LIGGGHTS now.

Bonded particle model with break and no bond breakage

Submitted by zealakis on Fri, 09/30/2016 - 02:58

hello,

i compiled the bondpackage model by Potyondy and Cundall and i could make the fracture simulation of a bonded particle model and also a simulation of a no bond breakage by specifying a very high bond-break radius.
But i haven't been able to simulate both behaviors at the same time for two bodies separately because i cant use a hybrid bond style with gran.

Can someone helpme with this?. i have little knowledge about c++ so its hard for me to hack de source code.

Any suggestion or question will be greatly appreciated.

Have a nice day,

Efield

Submitted by JohLoh on Thu, 09/29/2016 - 16:28

Hallo everyone,

I am interested in charged particles and their interaction with an external field. (comparable with the video "Electrostatic" on this website). Can I use the efield command in Liggghts? And how can I define a charged particle (attribute q)?

Thanks for your help. Best wishes

Johannes

finite plane wall

Submitted by llib11ik on Wed, 09/28/2016 - 14:22

I understand that using command
fix xwalls1 all wall/gran model hertz tangential history primitive type 2 xplane 0.0
can create a infinite X-plane wall at x=0.

Could anyone tell me how to create a finite plane wall? For example, a wall from (0, 0) to (0, 2) in 2D or in 3D formed by 4 points.
Or it can only be done by importing mesh file?

Thanks!

Jason

Pages

Subscribe to RSS - LIGGGHTS® - User Forum