How to show contact force on particles or polyhedre using changing colors

Asked by Kun Zeng

Hi,

I was wondering if there is a method to show contact forces with changing colors on paricles or polyhedre in paraview? Thanks.

Regards,

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Jan Stránský
Solved:
Last query:
Last reply:
Revision history for this message
Luc Scholtès (luc) said :
#1

Hi,

I am not sure to understand your request ("contact forces changing colors on particles"?) but it is possible to see the principal stresses applying on each particle by adding the 'bstresses' attribute in the VTKRecorder (along with the 'spheres' attribute of course). Then, you can select the stress you are interested in and color the glyph associated to each point (e.g. a sphere) according to the stress magnitude. If you want to see the interaction network, it is also possible by adding the 'intr' attribute in the VTKRecorder. Then, by selecting the component of the forces you want to see (normal or shear) and coloring it according to their intensity you can see how the forces are distributed in your assembly.

I don't know if these options work for polyhedron though.

Luc

Revision history for this message
Jan Stránský (honzik) said :
#2
Revision history for this message
Kun Zeng (zkbread) said :
#3

Hi,

Thanks for your suggestions.

I mean if it is possible to show the contact force on the surface of the polyhedra particles not using force chain.

Thanks

Revision history for this message
Best Jan Stránský (honzik) said :
#4

There is export.VTKExporter.exportContactPoints [1], which could suit your purpose. You can export force something like (not tested):
##
vtkExporter.exportContactPoints(what=[('force','i.phys.normalForce+i.phys.shearForce')])
##
and then use arrow glyph in paraview
Jan

[1] https://yade-dem.org/doc/yade.export.html#yade.export.VTKExporter.exportContactPoints

Revision history for this message
Kun Zeng (zkbread) said :
#5

Thanks Jan Stránský, that solved my question.