Using cohesion instead of Clumping

Asked by Seti

Hi All,

As per my previous questions I am thinking to define cohesion between particles instead of clump them together. However I am wondering If I can introduce cohesion just between some percentage of whole grains. I mean just some of the particles stick to each other and the others just have frictional contacts. And Can I control this ration same as what we are doing in clump?
Is there any scripts or documentation to help me in this regards.

Regards
Seti

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
Chareyre (bruno-chareyre-9) said :
#1

Hi,
The user manual explains how you can create individual interactions
manualy. On this basis you can do whatever you like in terms of type of
contacts.
I don't think there is a realistic example of this unfortunately.
Bruno

Revision history for this message
Seti (seti) said :
#2

Hi Bruno,

Sorry for any inconvenience, I have tried to put more detailed in my question as per below.

Hi All,

I am trying to model a cemented sample for compression test. The sample is made by two materials (1- Loose sand, 2- bonder , bounding the sand grains), obviously each material has own cohesion.
Based on YADE documentation the cohesion between two material is calculated as per below
a_i=C_i min(R_1,R_2)^2, i=n,s.

However I need to choose C _max for C_i instead of C_min.
I was told "MatchMaker" in Ip2_CohFrictMat_CohFrictMat_CohFrictPhys should do what I am looking for. I have tried below engine, however it does not work. Do you have any idea why the code is not working the error is as per below:

"AttributeError: No such attribute: shearCohesion."

I Just put a snip of code :

O.engines=[
        ForceResetter(),
        InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()]),
        InteractionLoop(
        [Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()],
                [Ip2_FrictMat_FrictMat_FrictPhys (),
                [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(normalCohesion=MatchMaker(matches=((mat1,mat2,500),(mat2,mat2,1e6))),shearCohesion=MatchMaker(matches=((mat1,mat2,50e6),(mat2,mat2,80e6))),label="cohesiveIp")], [Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True,includeAdhesion=True)],
                [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment(
   useIncrementalForm=True,
   always_use_moment_law=False,
   label='cohesiveLaw')]

Thanks,

Seti

Revision history for this message
Bruno Chareyre (bruno-chareyre) said :
#3

The error is rather explicit and unrelated to your particular problem:
"AttributeError: No such attribute: shearCohesion."

That's why asking "it doesn't work, help me" in the original post was like asking blind people to look for your keys.
It should be enough to grep "shearCohesion" in your script and fix it.
Bruno

Revision history for this message
Bruno Chareyre (bruno-chareyre) said :
#4

by "original post" I actually meant https://answers.launchpad.net/yade/+question/295814, sorry.

Can you help with this problem?

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

To post a message you must log in.