About compute wall/gran/local contact area

Submitted by thecomedian on Thu, 11/20/2014 - 01:35

Hi! i'm relatively new to liggghts, and i want to compute the contact area in a cone crusher. So far i'm using the "compute wall/gran/local contactArea" command, to store the information, unfortunately, i don't know how to use only a single mesh wall to compute the contact area.

The code that does the wall gran, computation and the dumping is:

fix granwalls all wall/gran model hertz tangential history mesh n_meshes 7 meshes cad1 cad2 cad3 cad4 cad5 spdr tapa
...
compute 2 all wall/gran/local contactArea
...
dump data all local 200 post/data*.vtk c_2[1]

as you can see, the compute command stores information of every interaction between mesh and atoms, and i only want to store information of the first mesh, i would aprecciate if someone help me with this.

Kind Regars

ckloss's picture

ckloss | Fri, 12/12/2014 - 13:17

Hi comedian,

such a feature is currently not implemented - but would be pretty easy to do!
See compute_pair_gran_local.h/cpp. All you'd have to do is check the mesh ID in the add_wall_* functions

Cheers
Christoph

anandmds's picture

anandmds | Wed, 04/13/2016 - 11:50

hey comedian,

Just have a look at the " dump mesh/vtk" command form the documentation. From what I understood, you can write the id of the mesh whose info you need.

Regards
Anand M

Anand.M
M.S. Mechanical Engineering,
I.I.T. Madras