tadded shearStrainRate and inertiaParameterPlanarShear functions - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
 (HTM) git clone git://src.adamsgaard.dk/sphere
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
 (DIR) commit 460bcf1cd69bf981cffc443468bd21e168fbfadf
 (DIR) parent c76e7b620aab56cb71874ccea4c6a0f2ba6880b2
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon,  1 Sep 2014 14:44:02 +0200
       
       added shearStrainRate and inertiaParameterPlanarShear functions
       
       Diffstat:
         M python/sphere.py                    |      28 ++++++++++++++++++++++++++++
       
       1 file changed, 28 insertions(+), 0 deletions(-)
       ---
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -3468,6 +3468,8 @@ class sim:
        
                :returns: The shear velocity applied by the upper, fixed particles [m/s]
                :return type: float
       +
       +        :see also:`func:shearStrainRate()`,`func:shearStrain()`
                '''
        
                # Find the fixed particles
       t@@ -3484,6 +3486,8 @@ class sim:
        
                :returns: The total shear strain [-]
                :return type: float
       +
       +        :see also:`func:shearStrainRate()`,`func:shearVel()`
                '''
        
                # Current height
       t@@ -3495,6 +3499,30 @@ class sim:
                # Return shear strain
                return xdisp/w_x0
        
       +    def shearStrainRate(self):
       +        '''
       +        Calculates the shear strain rate (dot(gamma)) value of the experiment.
       +
       +        :returns: The value of $I$
       +        :return type: float
       +
       +        :see also:`func:shearStrain()`,`func:shearVel()`
       +        '''
       +        return self.shearStrain()/self.time_current[0]
       +
       +    def interiaParameterPlanarShear(self):
       +        '''
       +        Returns the value of the inertia parameter $I$ during planar shear
       +        proposed by GDR-MiDi 2004.
       +
       +        :returns: The value of $I$
       +        :return type: float
       +
       +        :see also:`func:shearStrainRate()`,`func:shearVel()`
       +        '''
       +        return self.shearStrainRate() * numpy.mean(self.radius) \
       +                * numpy.sqrt(self.rho[0]/self.w_devs[0])
       +
            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