A simple question on interactionDetectionFactor

Asked by Huihuang Xia

Hi,

https://yade-dem.org/doc/yade.wrapper.html?highlight=interaction%20detection#yade.wrapper.Ig2_Sphere_Sphere_ScGeom.interactionDetectionFactor
Here gives a parameter that create interaction range for particle packing, but I have no idea on how to determine this parameter ?
Because in this paper (A DEM model for soft and hard rocks: Role of grain interlocking on strength), it was stated that if Rmax/Rmin=2, interactionDetectionFactor must not be larger than 1.5. In my simulations, Rmax/Rmin=9, thus, I have no idea to determine the interactionDetectionFactor. I have tried interactionDetectionFactor=1.5 to do some simulations, it got a relatively good result. But this selection (interactionDetectionFactor=1.5) seems go against the theory put in that paper.

Can anyone help me solving this confusion ?

Thanks in advance,

Huihuang

Question information

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

Hello,

We mentioned intR<1.5 for Rmax/Rmin=2 so the near neighbour interactions keep conceptually "representative" of a rock microstructure even though, I agree with you, it was not clearly explained. It is not easy to explain but I'll try anyway:

Lets consider a line of 3 particles 1, 2, 3 such that R1=2, R2=1 and R3=2. If you create bonds such as intR>1.5, then 1 and 3 will be bonded together. A bond will thus be created "over" 2.

In granular materials such as granite or sandstone (where only contacting grains interact one with each other), this is physically not acceptable. Then, if your idea is to discuss the influence of the microstructure of the medium on its overall behavior, this might be problematic. I must precise that this sort of long range interactions might nonetheless exist in some materials with exotic grain to grain interactions (e.g. Van Der Waals forces).

Now, even though this can be conceptually not exactly representative of the way grains interact in a geomaterial, it might nonetheless be an acceptable way to overcome some difficulties faced by classic bonded particle models to deal with certain types of cohesive materials (under the condition that no direct link is done between the numerical particles and the real grains). Another way to deal with such "non physically acceptable" interactions would be to use an enhanced contact law that would increased the degree of interlocking of grains (e.g. flat joint contact model in PFC, or moment contact law in YADE).

Luc

Revision history for this message
Best Robert Caulk (rcaulk) said :
#2

I think the idea of an interaction range is brilliant, it opens up a huge area of functionality that did not previously exist in yade.

I too found the sentence in that paper a little brief, so I took some time to derive the equation. Here was my thought process:

Worst case scenario is that two maximum-sized particles might be interacting. We want to avoid allowing a minimum sized particle between two interacting particles. So we know the maximum allowable Deq<2*R_min+2*R_max (distance between 2 max sized particle centers with a minimum sized particle in between).

And the maximum possible Deq for a given packing

Deq = gamma_int*(R_max+R_max).

Solving for gamma_int to satisfy both conditions:

gamma_int < (R_min+R_max)/(R_max)

So in your case (R_max/R_min = 9), your max gamma_int would be 1.11 to avoid overlapping small particles.

It is probably worth saying that in my experience you can go a bit above the maximum without seeing any negative consequences. Actually, I have noticed that I [i]need[/i] to go above the maximum to calibrate toughness...hmmm.

Cheers,

Robert

Revision history for this message
Huihuang Xia (huihuangxia) said :
#3

Thanks Luc and Robert, you solved my problem.

Revision history for this message
Huihuang Xia (huihuangxia) said :
#4

Thanks Robert Caulk, that solved my question.