 |
tadd function to find loaded contacts - sphere - GPU-based 3D discrete element … |
 |
git clone git://src.adamsgaard.dk/sphere |
 |
Log |
|
Files |
|
Refs |
|
LICENSE |
|
--- |
|
commit c71ae0205444e3fe0a860127912ea69c98d4a924 |
|
parent c5cca9ab322bea5b1b0e09da0416e677f93a2e3a |
|
Author: Anders Damsgaard <[email protected]> |
|
Date: Fri, 20 Feb 2015 09:08:13 +0100 |
|
|
|
add function to find loaded contacts |
|
|
|
Diffstat: |
|
M python/sphere.py | 14 ++++++++++++++ |
|
|
|
1 file changed, 14 insertions(+), 0 deletions(-) |
|
--- |
|
diff --git a/python/sphere.py b/python/sphere.py |
|
t@@ -4347,6 +4347,20 @@ class sim: |
|
self.findAllContactSurfaceAreas() |
|
self.sigma_contacts = self.f_n_magn/self.contact_area |
|
|
|
+ def findLoadedContacts(self, threshold): |
|
+ ''' |
|
+ Finds the indices of contact pairs where the contact stress magnitude |
|
+ exceeds or is equal to a specified threshold value. This function calls |
|
+ :func:`findContactStresses()`. |
|
+ |
|
+ :param threshold: Threshold contact stress [Pa] |
|
+ :type threshold: float |
|
+ :returns: Array of contact indices |
|
+ :return type: array of ints |
|
+ ''' |
|
+ self.findContactStresses() |
|
+ return numpy.nonzero(self.sigma_contacts >= threshold) |
|
+ |
|
def forcechains(self, lc=200.0, uc=650.0, outformat='png', disp='2d'): |
|
''' |
|
Visualizes the force chains in the system from the magnitude of the |
