yade-users team mailing list archive

[Henry <question254640@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Henry is still having a problem:
    Many thanks to Dr. Bruno.
    I have read some references on the contact model of the cohesive material in DEM, for example the algorithm in PFC and other DEM code. And I have found that there are some different on the calculation of the normal and shear adhesion forces:
    In PFC the normal and shear adhesion fores are calculated according to the "area" of the contact ( which is understanded easily):
       normalAdhesion =min(normalCohesion_A,normalCohesion_B)*Pi*pow(min(D_B, D_A),2);
      shearAdhesion =min(shearCohesion_A,shearCohesion_B)*Pi*pow(min(D_B, D_A),2);
   However in YADE the adhension forces are calculated according the "square" of the contact diameter[1]:
    
      contactPhysics->normalAdhesion = std::min(sdec1->normalCohesion,sdec2->normalCohesion)*pow(std::min(Db, Da),2);
      contactPhysics->shearAdhesion = std::min(sdec1->shearCohesion,sdec2->shearCohesion)*pow(std::min(Db, Da),2);
    
    Is there any introductions on the reasons of this different?

Thanks a lot!

Henry



[1] http://bazaar.launchpad.net/~yade-pkg/yade/git-trunk/view/head:/pkg/dem/CohesiveFrictionalContactLaw.cpp

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