yade-dev team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

> could you tell me which is in brief the difference between ScGeom and 
> Dem3DofGeom?
ScGeom uses incremental computation of shear, whereas Dem3DofGeom uses 
total formulation. ScGeom replaces facet with sphere in facet-sphere 
contact, Dem3DofGeom has proper representation of facet with 0 thickness.

I was suggesting to refactor those 2 (i.e. write it properly and have 
the same interface for both incremental and non-incremental version (so 
that the law would be the same for both), but there was not much reaction...

> BTW, about the naming for a new contact law, should I follow the convention 
> Law2_>>_>>_>>, shouldn't I?
Yes, please
.
> If I have new variables, could I include them in the FrictPhys class? 
Derive your own class from FrictPhys, just like e.g. CpmPhys does.

> Should I use the virtual function go? I see that in the ElasticContactLaw used 
> in the Triaxial we use the function action..
ElasticContactLaw is old engine that puts the interaction loop inside 
itself and internally calls the elastic LawFunctor anyway. It makes the 
code slower.

Instead, derive your class from LawFunctor and overrider the ::go 
method, which will always act on 1 interaction only. (Read 
https://www.yade-dem.org/sphinx/prog.html#multiple-dispatch, it can be 
useful although it gets in lots of details on the c++ side)

v