yade-users team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

Hi all (Chiara, check this out :)),

I reorganized the code for moment/rotation that was used in 
CohesiveFrictionalLaw,  Law2_ScGeom_CohFrictPhys_CohesionMoment, and 
other moment laws using more or less the same code, in order to 
precompute bending and twisting in the intercation geometry.

There is now a ScGeom6D that you can use to avoid duplicating all this 
rotational code in contact laws. You can see example usage in  
CohesiveFrictionalContactLaw.cpp:204 :

/currentContactPhysics->moment_twist = 
(currentContactGeometry->getTwist()*currentContactPhysics->kr)*currentContactGeometry->normal;
currentContactPhysics->moment_bending = 
currentContactGeometry->getBending() * currentContactPhysics->kr;/

Using this code implies usage of Ig2_Sphere_Sphere_ScGeom_6D_, instead 
of Ig2_Sphere_Sphere_ScGeom.

Some renaming will happen in other commits (e.g. 
Law2_ScGeom_CohFrictPhys_CohesionMoment => Law2_ScGeom6D_...).

Cheers.

Bruno



--

_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : +33 4 56 52 86 21
Fax : +33 4 76 82 70 43
________________