tcontact area instead of contact radius - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
  HTML git clone git://src.adamsgaard.dk/sphere
   DIR Log
   DIR Files
   DIR Refs
   DIR LICENSE
       ---
   DIR commit ac4846ed72cbc22f07e21a53f2c6a61daca68c4e
   DIR parent d1fbd86f51a9281a5e16d282462a3c097212bcc0
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 19 Feb 2015 13:48:18 +0100
       
       contact area instead of contact radius
       
       Diffstat:
         M python/sphere.py                    |       3 ++-
       
       1 file changed, 2 insertions(+), 1 deletion(-)
       ---
   DIR diff --git a/python/sphere.py b/python/sphere.py
       t@@ -4315,10 +4315,11 @@ class sim:
                r_i = self.radius[i]
                r_j = self.radius[j]
                d = r_i + r_j - overlap
       -        return 1./(2.*d)*(
       +        contact_radius = 1./(2.*d)*(
                        (-d + r_i - r_j)*(-d - r_i + r_j)*
                        (-d + r_i + r_j)*( d + r_i + r_j)
                        )**0.5
       +        return numpy.pi*contact_radius**2.
        
            def findAllContactSurfaceAreas(self):
                '''