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 <[email protected]> wrote: > > > On 3 December 2010 14:57, Bruno Chareyre <[email protected]>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 : [email protected] >> Unsubscribe : >> https://launchpad.net/~yade-dev<https://launchpad.net/%7Eyade-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
