yade-users team mailing list archive

[Jérôme Duriez <question700181@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jérôme Duriez proposed the following answer:
Hi,

Thinking of Law2_ScGeom_FrictPhys_CundallStrack.plasticDissipation, the
latter does not care of interacting bodies' shape so you can not really
use it.

You can nevertheless probably get inspiration from corresponding source
code [*] to recompute, all by yourself e.g. from Python, the "elastic
trial" shear force and get yourself the plastic dissipation for your own
subset of contacts. After checking e.g.

isinstance(O.bodies[cont.id1].shape,Sphere) and
isinstance(O.bodies[cont.id2].shape,Sphere) # cont = a YADE interaction


[*] https://gitlab.com/yade-dev/trunk/-/blob/master/pkg/dem/ElasticContactLaw.cpp#L92 and above

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