Yade

Calculate the energy released when the bond between two particles breaks

Asked by Yor1

Hello

I want to calculate the energy released when a microcrack was created.
In fact, a microcrack was created mean that a bond between two particles was broken.
My question is how can we detect the particles which their bonds were broken in YADE?

Best regards
Jabrane

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Jérôme Duriez
Solved:
Last query:
Last reply:
Revision history for this message
Best Jérôme Duriez (jduriez) said : 		#1

Hello,

The answer to your question directly depends on which YADE contact model you're using. Assuming you're using the JCFpm model, the state of a given particle returns how many initial bonds on this particle have been broken, either in tension or in shear:

- https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.JCFpmState.tensBreak
- https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.JCFpmState.shearBreak
With their relative counterparts suffixed "Rel" (one line below in the doc)

Revision history for this message
Yor1 (jabrane-hamdi) said : 		#2

Thanks Jérôme Duriez, that solved my question.

Revision history for this message
loiseaurare (loiseaurare) said : 		#3

Hi Jerome,

Do you have any more informations about that state function ?
Like, what does mean "broken" from that point of view ? Does it consider the initial cohesive interactions only ?
I wasn't able to get that when looking at the yade code...

Revision history for this message
Jérôme Duriez (jduriez) said : 		#4

Hi,

A cohesive interaction in JCFpm is said to be broken in tension resp. shear when its tensile resp. shear force threshold is reached for the first time (afterwards, cohesion is lost: there is no more tensile strength or shear cohesion in these thresholds, see relevant papers quoted in the doc)

If a particle initially has 8 cohesive interactions, and than 3 of them break in tension during the time interval [0,t],
state.tensBreak = 3 at time t, and state.tensBreakRel = 3/8 (similar thing for shear variables)

The corresponding code is in JCFpm Law2 source code, you may search for "tensBreak" in gitHub Yade webpage

Jerome

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

Hi,

Have a look at the latest commit where all you need is available directly from the JCFPM law functor. Have a look at the example provided (examples/jointedCohesiveFrictionalPM/uniaxialCompressionTest_JCFPM.py).

Cheers

Luc