yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #26934
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.