Hello,
I have problems with understanding a few commands:
The simulation issue is a box filled with particles where I want to analyse stress on different heights of the walls.
I use these commands
fix cad all mesh/gran/stressanalysis cube 1 1 0 0 0 0 0 0 finnie yes
fix GlobalforceOnMesh all ave/time 50 2 100 f_cad[1] f_cad[2] f_cad[3] file filename.txt
dumpstress all mash/gran/VTK 200 post/dumpCAD*.vtk stress wear cad
The manual says, that the fix mesh/gran/stressanalysis stores a global vector with 6 components, where the first 3 components are equal to the total force on the mesh in x-, y- and z-direction.
I have a primitive stl file for my box and I can visualize the pressure on the (in total 12) triangles in paraview.
1) What I don't understand is, if f_cad would calculate just a global force vector (3 scalars each timestep for the 3 directions of a 3d simulation), it shouldn't be possible for paraview to visualize the varying pressure on the two triangles of one of the y-walls of the box...
2) How can I access the forces on the mesh depending on the height of the box?
Thanks in advance,
kind regards,
Silias
P.S.: I assume, that it's not possible to calculate these force by looking at the per-atom forces, because they are equal zero, while the bulk material is more or less settled.
raguelmoon | Tue, 05/10/2011 - 04:07
Hi, <>fix GlobalforceOnMesh
Hi,
<>fix GlobalforceOnMesh all ave/time 50 2 100 f_cad[1] f_cad[2] f_cad[3] file filename.txt
Because this statement calculate force on whole CAD.
>>2) How can I access the forces on the mesh depending on the height of the box?
If you want to calculate force on each segment of height then you have to break CAD in as many pieces as forces are required and then import all pieces hight-wise, and use this statement to calculate force on each pieces:
fix GlobalforceOnMesh all ave/time 50 2 100 f_cad[1] f_cad[2] f_cad[3] file filename.txt
I think Christoph can answer in better way.
Cheers,
Ram
Silias | Tue, 05/10/2011 - 11:11
How are the stresses calculated?
Thanks Christoph and Ram for your replys!
Ok, I'll try the two possibilies you suggested:
1. Access the forces of each segment by cutting the CAD in several pieces and use the fix ave/time to output the forces on each CAD-piece
2. Access the forces in paraview via vtk-files (which sounds more straight-forward)
Trying to understand, what kind of output the commands
fix mesh/gran/stressanalysis and
dump mesh/gran/vtk
produce exactly, I looked at the vtk-files.
Using the simple 1 particle model, I realized that stresses (unequal zero) are written to the vtk-files just in these timesteps, while the particles lies on the floor (with negligible motion/bouncing) and there is one value for each triangle of the mesh.
What I don't understand is, in what way is the value of the stress calculated?
It is not the gravitational force divided by the surface of the triangle...
Thanks in advance,
Sebastian
ckloss | Tue, 05/10/2011 - 11:15
The stress is simply the
The stress is simply the force on the triangle divided by its area
Christoph
Silias | Tue, 05/10/2011 - 11:40
sorry, calculation error, i
sorry,
calculation error, i took the wrong radius...
thanks! :)
ckloss | Tue, 05/10/2011 - 08:44
Hi Sebastian, 1) as Ram said,
Hi Sebastian,
1) as Ram said, it is the total force on the mesh that is calculated by tthe fix command, the dump mesh/gran/VTK command dumps the per-triangle stress
2) I am pretty sure there is some way in paraview to extract that, or you can write your own postprocessing routine
Christoph
xpt | Mon, 06/20/2016 - 07:45
Paraview force
Hi Sebastian!
Do you find a way to get the force in Paraview?