yade-dev team mailing list archive

[Chiara Modenese <c.modenese@xxxxxxxxx>]

Bruno,
I understand now the formulation of Plassiard. In the paper he is used to
define the contact stiffnesses in a different way than we do. He is using
the mean radius to get them whereas we employ the radius of each particle in
the calculation. Hence, in our case, it makes sense to use Da*Db as you were
doing. I will delete Modenese's formulation (not consistent at this point).
One last thing I need to know: how would you define the plastic moment?
Plassiard uses
Mmax=eta*Rmean*|Fn|
but in place of Rmean we could perhaps simply take the min(Da,Db)? It does
not make a big difference anyway, just we would avoid to use Rmean. Let me
know.
Chiara


On 3 December 2010 15:16, Chiara Modenese <c.modenese@xxxxxxxxx> wrote:

>
>
> On 3 December 2010 14:57, Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>wrote:
>
>>
>> >     - I restored Plassiard definition of Kr. Your definition is still
>> >     optionaly available
>> >     (flag useMeanRad).
>> >
>> > Actually, I got to change it because I was following Plassiard's
>> > definition. In his paper (2009), he is taking the average of the two
>> > radii. But it is fine, since I agree with what you say below.
>>
>> Huh?!! Thanks for info, then we have to move the comment with Plassiard's
>> paper reference
>> to the variant, and call default formula Modenese variant!
>>
> Ahah! That is funny. I will update you before to commit it (sorry for the
> regression test, I completely forgot about them, but it is nice to have
> them). Also, if you do not use the radius mean, there would be the problem
> on how to define the maximum plastic value (again, Plassiard was using the
> radius mean x Fn). We need to make it consistent either we use the radius
> mean or not. I let you shortly know.
>
>>
>> >     - There used to be some physical parameters (cohesion and others) in
>> >     Ip2_2xCohFrictMat_CohFrictPhys, but I recently moved them all to
>> >     FrictMat. If you have a
>> >     moment, could you please do the same for alphas and etaRoll (and
>> >     probably harmonic-average
>> >     them in the functor)? Only flags should be left in the functor.
>> >
>> > No, problem. I am doing it now.
>>
>> Cool.
>> I'm running back to teaching.
>>
>> Cheers.
>>
>> B.
>>
>> _______________________________________________
>> Mailing list: https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
>> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
>> Unsubscribe : https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-dev>
>> More help   : https://help.launchpad.net/ListHelp
>>
>
>