-------- Original Message -------- Subject: Re: [Yade-dev] Is shear incremental calculation verified Date: Sun, 17 Jan 2010 08:44:47 +0100 From: eudoxos <[email protected]> To: Bruno Chareyre <[email protected]>
I was actually thinking about implementing contact geometry classes from scratch, both incremental and non-incremental versions, in such way that deformations (normal, shear, twist, bending) would be precomputed in the Ig2_* functor, and then just accessed directly from the constitutive law. The InteractionGeometry classes themselves would be the same, only the functor would make the distinction of absolute/incremental. There already was 6dof contect implemented (Dem6Dof_Sphere_Sphere), but it is commented out currently; and guess what, it is much cleaner code than the horrible mess (sorry) in CohesiveFrictionalContactLaw. I might do just some part of it, let's see. Cheers, Vaclav > >> From my simulations, updateShearForce is correct (I am able to >> reproduce the results obtained by Plassiard et al's paper). > Always good to hear this sort of verification! It's also interesting to > have rotational friction implemented in Yade. > >> When I use the codes in updateShear, the deviator stress I get is >> much lower. >> > Sure it was not due to unconsistent copy/pastes? updateShear is > connected to other classes where the computed values are used. The two > versions are based on different logics, you can't just switch from the > one to the other and expect similar results I think. > > Bruno > > > _______________________________________________ > Mailing list: https://launchpad.net/~yade-dev > Post to : [email protected] > Unsubscribe : https://launchpad.net/~yade-dev > More help : https://help.launchpad.net/ListHelp _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
