yade-users team mailing list archive

Thread
Date

[Question #270503]: Are the formulae for adhesions in CohesiveFrictionalContactLaw correct?

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Christopher Stanbridge <question270503@xxxxxxxxxxxxxxxxxxxxx>
Date: Tue, 18 Aug 2015 13:46:44 -0000
Reply-to: question270503@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

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

Hi all,

I have noticed that line 279 of CohesiveFrictionalContactLaw.cpp (here: https://github.com/yade/trunk/blob/48bb1a66dd0713de1a51a0c3562399133c87fd5e/pkg/dem/CohesiveFrictionalContactLaw.cpp#L279) gives the formula for normal adhesion calculation as follows:
contactPhysics->normalAdhesion = std::min(sdec1->normalCohesion,sdec2->normalCohesion)*pow(std::min(Db, Da),2);

If adhesion is being calculated between two spheres, this means that the effective area of contact between the spheres is taken to be equal to the square of the radius of the smallest sphere. Should it really be pi*r^2 instead of r^2?

Cheers,
Christopher Stanbridge

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.