>> 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 _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
