Question #234129 on Yade changed: https://answers.launchpad.net/yade/+question/234129
lingran posted a new comment: Hi all, As promised, I try to add bendingElastEnergy and twistElastEnergy to Law2_ScGeom6D_CohFrictPhys_CohesionMoment. As a result, the corresponding cpp and hpp files are slightly modified. In addition, the new source codes are recompiled to update yade. Several simple 2D tests are conducted to check this extended function. The results show that when Law2_ScGeom6D_CohFrictPhys_CohesionMoment is only used, bendingElastEnergy and twistElastEnergy values are verified and consistent with calculated values, while in the other hand, when Law2_ScGeom6D_CohFrictPhys_CohesionMoment is combined with Law2_ScGeom_FrictPhys_CundallStrack, bendingElastEnergy is given as ‘nan’ or ‘inf’ in yade. Here are the diff files: cpp_diff_for_cohesiveLaw: http://pastebin.com/eWzpLdGC hpp_diff_for_cohesiveLaw: http://pastebin.com/XD18Wktn Here are the modified source codes: cpp_add bending_to_cohesiveLaw: http://pastebin.com/XQADy5Vt hpp_add bending_to_cohesiveLaw: http://pastebin.com/GuUiL2g0 Hope it helps. Best regards. Lingran -- You received this question notification because you are a member of yade-users, which is an answer contact for Yade. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp

