yade-dev team mailing list archive

[Vincent Richefeu <richefeu@xxxxxxxxx>]

> Note that the integral is computed on deformed contour.
> In <tblatex-1.png>, <tblatex-2.png> is the contact point, closer to center than radius. Hence, dividing by a volume of the Voronoi sort would not be consistent with the integration domain. I don't see what other volume could be used than the one of the sphere actually (exception for Hertz-Mindlin experts who know what is the exact very small volume change, as a function of Poisson ratio, in the vicinity of a contact).

The volume must obviously be greater or equal to the particle-volume.
A function that provide the moment tensor (\sum_k x_i^k f_j^k without division by a volume) for one particle 
can be used with any shape (avoid the volume computation by means of radii). 
It can also be summed and then divided by the volume you want...