yade-users team mailing list archive

[Leonard <question700258@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I am simulating a triaxial compression test using Hertz-Mindlin model, and I'd like to ask that if it is possible to include adhesion in Hertz-Mindline contact model so that I can simulate bonded material (or cement material).

I read from MindlinPhys[1], there seems to be an option to set adhesion (named adhesionForce in [1]).

#I use FrictMat: 
MatSand = O.materials.append(FrictMat(young=young,poisson=0.3,frictionAngle=radians(30),density=2650,label='sand'))
#The engines I used are as follow:
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,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,
	#PyRunner(iterPeriod=2000,command='history()',label='recorder'),
	#PyRunner(command='stopIfDamaged()',iterPeriod=500),
]

#I tried to use the following command to active adhesion:
ContactModel.isAdhesive=True
ContactModel.adhesionForce=adhesion

But I didn't make the adhesion come into effect. Thereby, I'd like to ask that:
(1) How to include adhesion in Hertz-Mindline contact model if it is possible.
(2) In [1], it explains adhesionForce as "Force of adhesion as predicted by DMT", what is DMT?
(3) Where I can find the source code of Hertz-Mindline contact model in Yade trunk?

Thanks!
Leonard

[1]https://yade-dev.gitlab.io/trunk/yade.wrapper.html?highlight=dmt#yade.wrapper.MindlinPhys

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