Yade

Why the interaction(force chain) between particles and facet all focusing on one point on the facet

Asked by Kun Zeng

I set up a model by using polyhedra partcles and use export.VTKExporter(name).exportInteractions() to export the interaction to VTK file. However, I saw the result in paraview and also in yade, and find the interaction between particles and on certain facet all concentrate on point of one facet.

Does anyone know what's the problem with my simulaton?

Here is my input code and .VTK file

https://drive.google.com/open?id=14UARn-UNa4yPsFmGZeJFyaNtX0I19l1P

Question information

Language:
English Edit question
Status:
Answered
For:
Yade Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Jan Stránský (honzik) said : 		#1

Hello,
the easiest approach (what you are experienceing) is to connect state.pos of two bodies. In the case of facet, its pos is in the center of inscribed circle.
If you want something more detailed, you can implement your own postprocessing (projection of the sphere's center to the facet) or if the overlaps are small, contact points can also be useful for your problem.
Cheers
Jan

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

Hi Jan,

I used polyhedra particles.

[1] As you suggested, how can I do by using what "command" in Yade for dealing with the facet case ?

[2] You mean contact points? Is this one? vtkExporter.exportContactPoints(what=[('force','i.phys.normalForce+i.phys.shearForce')]?

How to show the contact tube in paraview? is that same operation with vtkExporter.exportInteraction() ?

Thanks.

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

Hi Kun,
In yade, currently there is no complete tool to do what you want.
But you have access to all necessary information (contact points of interactions, vertices of facets, vertices of polyhedrons), so you can create your own postprocessing tool (save it to vtk file for instance).
cheers
Jan

Revision history for this message
De zhang (dzlearnyade) said : 		#4

Hi Kun,
I met the same problem of yours, I think Jan supplied a good method, and I had another way that I dispersed the facet into small ones so that you can get more detail interations of polyhedra and facet, but the O.dt will be vary small as the massive number of facets.
cheers
De Zhang

Can you help with this problem?

Provide an answer of your own, or ask Kun Zeng for more information if necessary.

To post a message you must log in.