yade-users team mailing list archive

[chiara modenese <c.modenese@xxxxxxxxx>]

On 28 July 2010 18:34, chiara modenese <c.modenese@xxxxxxxxx> wrote:

>
>
> On 28 July 2010 10:56, nejib hadda <nejib.hadda@xxxxxxxxxxx> wrote:
>
>> Hi,
>>
>>
>> I'm using the "ElasticContactLaw" for the interactions, I want to
>> extract the relative normal and shear displacements and relative
>> rotations in the local (contact) coordinate system.
>>
>
Hi Nejib,
Using contact law functors (Law2) based on ScGeom, you can access the total
normal and relative shear displacements in this way:

i=O.interactions[id1,id2] # get the interaction
un=i.geom.penetrationDepth # computed in Ig2s functors
us=i.geom.shearIncrement # computed inside function precompute() of ScGeom
class.

Note that above un is total, whereas us is the increment. You could also
easily calculate the total shear displacement into the contact law using
rotate() function of ScGeom, in case you need it.

Chiara