yade-users team mailing list archive

[Bruno Chareyre <question694531@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Bruno Chareyre proposed the following answer:
Hi,
In the context of CPM the espT can be seen as a proxy to current shear force Ft, since the model is linear. You could replace epsT by 
Ft (and "rotate(Ft)") everywhere and the result would be the same.

In the incremental algorithms (Law2_ScGeom6D_CohFrictPhys_CohesionMoment
and actually nearly all functors other than CPM, including non-linear
ones such as Hertz-Mindlin or Lubrication), epsT is not defined. But  Ft
is, and so we reflect rigid body motion on current Ft directly. It does
not imply hard limit on the magnitude of shear or rotation.

What you need is to formulate your hyperelastic material in terms of
dFt/dut as a function of current _Ft_ (not current ut), which should be
possible (hopefully). If you have this you can implement directly into
an incremental algorithm.

I hope it helps.

Bruno

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