← Back to team overview

yade-users team mailing list archive

Re: [Question #700258]: How to involve adhesionForce in Hertz-Mindlin contact model

 

Question #700258 on Yade changed:
https://answers.launchpad.net/yade/+question/700258

    Status: Answered => Open

Leonard is still having a problem:
Hi Karol,

Thanks for your reply which really helps me.

I added gamma in Ip2_FrictMat_FrictMat_MindlinPhys as you suggested:
O.engines=[
 ForceResetter(),
 InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
 InteractionLoop(
  [Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom()],
  [Ip2_FrictMat_FrictMat_MindlinPhys(betan=betan,betas=betas,gamma=gamma,label='ContactModel')],
  [Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=False,includeAdhesion=True,label='Mindlin')]
 ),
 GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
 triax,
 TriaxialStateRecorder(iterPeriod=100,file='WallStresses'),
 newton,
]

If I understood correctly, the value of gamma controls the adhesion force between particles. And move one step forward, I'd like to ask that if it is possible to assign different gamma to particles with different material?
I remember MatchMaker may be able to do this, so that I tried:
ContactModel.isAdhesive=True
ContactModel.gamma=MatchMaker(matches=((Mat1,Mat1,0),(Mat1,Mat2,0),(Mat2,Mat2,gammaValue)))

But I got:
ArgumentError: Python argument types in
    None.None(Ip2_FrictMat_FrictMat_MindlinPhys, MatchMaker)
did not match C++ signature:
    None(Ip2_FrictMat_FrictMat_MindlinPhys {lvalue}, double)


I guess it is because that gamma is not a MatchMaker item, if so, I'd like to ask that how to assign different gamma to different material interactions.

Thanks.
Leonard

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.