tcontact area instead of contact radius - sphere - GPU-based 3D discrete elemen… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit ac4846ed72cbc22f07e21a53f2c6a61daca68c4e | |
parent d1fbd86f51a9281a5e16d282462a3c097212bcc0 | |
Author: Anders Damsgaard <[email protected]> | |
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(-) | |
--- | |
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): | |
''' |