yade-dev team mailing list archive

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

>> If the constitutive law is simple as
>>
>> // loop over interactions etc
>> Fn=contactGeom->strainN()*crossSection*E;
>> Ft=contactGeom->strainT()*crossSection*G;
>> // apply forces here
>>
>> the code for computing strain being elsewhere, then that code is still
>> in one place (and can be fixed at one place) and the constitutive law is
>> not a mess of flags.
>>     

Hi there, I implemented ElasticContactLaw2 (for ilustration) that
behaves like ElasticContactLaw (except that all contact are cohesive
(don't break) and there is no moment rotation law) and uses new
SpheresContactGeometry code; the code is about 15 lines.

This is the way I imagine all constitutive laws should eventually look
like: relationship from relative displacements (computed in the contact
geometry class) to forces, nothing else.

Vaclav