yade-users team mailing list archive

[lingran <question234129@xxxxxxxxxxxxxxxxxxxxx>]

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.