Yade

Histogram of crack orientation

Asked by Luis Barbosa

Hi,

I'm recording all cracks from JCFpmPhys_JointedCohesiveFrictionalPM. What I want to plot is the histogram of crack orientation in diferent time steps.

From the saved file we have break iteration, the 3 coordinates of the contact point, the type, the mean radius of the 2 spheres, the 3 coordinates of the contact normal and nrg (I don't know what is it).

I think that the 3 coordinates of the contact normal is what I want, but I am not sure how it is calculated. For instance, It is from the moment of bond break?

Any suggestion to find out how is it calculated?
Thanks,
Luis

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Luis Barbosa
Solved:
Last query:
Last reply:
Revision history for this message
Luis Barbosa (luis-pires-b) said : 		#1

Actually, I have to calculate the angle of each crack in relation to a base.

Revision history for this message
Luc Scholtès (luc) said : 		#2

Hi Luis,

The 3 coordinates of the contact normal is what you are looking for. They give the coordinates of the unit vector normal to the crack plane. The crack plane corresponds to the surface perpendicular to the branch vector joining the particle centers (it is located right at the location where the breakage occurs, as indicated by the 3 coordinates of the contact point). Based on the unit vector, you should be able to orientate the cracks in a reference frame (most probably, the same reference frame in which the vector coordinates are defined).

If you come up with a clean postprocessing script doing what you want, would it be possible to share it with the community? I believe it would be of great benefit for JCFPM users.

"nrg" gives the energy released by the breakage of the bond. This energy corresponds to the elastic energy stored in the bond before it breaks (in both normal and shear springs).

Luc

Revision history for this message
Luis Barbosa (luis-pires-b) said : 		#3

Ok, thanks. I will work on it.