← Back to team overview

yade-users team mailing list archive

Re: harmonic average of stiffness

 

Damned, you are right!
I reflected that in the documentation of the class. No change in the code itself.
Thank you for this remark.

Bruno

p.s. Commiting energy in minutes.


chiara modenese a écrit :


On 21 April 2010 14:53, Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx <mailto:bruno.chareyre@xxxxxxxxxxx>> wrote:



    My simplest explanation is actually in the documentation of the
    class, let me know if it is not clear enough :

    /"The compliance of one sphere under symetric point loads is
    defined here as 1/(E.D), with E the stiffness of the sphere and D
    its diameter, and corresponds to a compliance 1/(2.E.D) from each
    contact point. The compliance of the contact itself will be the
    sum of compliances from each sphere, i.e. 1/(2.E.D1)+1/(2.E.D1) in
    the general case, or 1/(E.D) in the special case of equal sizes.
    Note that summing compliances corresponds to an harmonic average
    of stiffnesss, which is how kn is actually computed in the
    :yref:`Ip2_FrictMat_FrictMat_FrictPhys` functor."
    /

Hi Bruno, in your explanation you reference the particles diameters, in fact in Yade we use the radii so that if particles have same geometry k_eq=E.R. I guess it does not really change the substance of the matter.
cheers, Chiara

------------------------------------------------------------------------

_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp


--
_______________
Bruno Chareyre
Associate Professor
Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : 33 4 56 52 86 21
Fax : 33 4 76 82 70 43
________________




References