Hi everyone!
I'm new in LIGGGHTS and to begin I see the example chute_wear, in this , I have a question, how i can get the contact force between the chute and the particles?
With pressure and area I can calculate force per cell, but it's a scalar value, and I would like to sum all of this to get a resultant force on the chute per time step.
I've done this ( like the graph), but are in different directions, thath is wrong :(
in summary, how i can get vectorial forces from LIGGGHTS (not pressure)? or, how i can transform the pressure in the vectorial force with Paraview?
Thanks!
| Attachment | Size |
|---|---|
 paraview.png | 132.33 KB |
aaigner | Fri, 10/23/2015 - 10:33
mesh/surface/stress
Hi xpt,
since the chute is a mesh/surface/stress you can get the total force and torque acting on the chute directly.
For more details check the point Restart, fix_modify, output, run start/stop, minimize info: in the documentation.
Bests
Andreas
xpt | Mon, 06/20/2016 - 05:46
Hi!
Hi!
I saw the documentation, but i don't understand how to do it :(
It says:
This fix stores a global vector with 9 components for access by various output commands. The first 3
components are equal to the total force on the mesh, the next 3 components store the total torque on the body
exerted by the particles. Finally, the last 3 components are the coordinates (moved, scaled, rotated) of the
reference point. Other info see fix mesh
So, I added this line after the other dumps:
dump dumpmove all local 200 post/dumpmove*.move index f_cadAnd this doesn't work :(
xpt | Mon, 06/20/2016 - 07:28
I answer to myself
Just i must to use fix ave/time:
fix baseforce all ave/time 50 2 100 f_cad1[1] f_cad1[2] f_cad1[3] file file.dat