yade-dev team mailing list archive

[eudoxos <eudoxos@xxxxxxxx>]

-------- Original Message --------
Subject: Re: [Yade-dev] Is shear incremental calculation verified
Date: Sun, 17 Jan 2010 08:44:47 +0100
From: eudoxos <eudoxos@xxxxxxxx>
To: Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>

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     : yade-dev@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-dev
> More help   : https://help.launchpad.net/ListHelp