yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #03085
Re: harmonic average of stiffness
On 21 April 2010 14:53, Bruno Chareyre <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
Follow ups
References