tdilation in mean grain diameters as default - 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 d4606c13b8efe33dfd5d4dc972f9221374aaef5e
DIR parent 565130ebc1f63424aff2f72cf769fa7ed75aa7f3
HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Mon, 31 Mar 2014 13:29:51 +0200
dilation in mean grain diameters as default
Diffstat:
M python/sphere.py | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
---
DIR diff --git a/python/sphere.py b/python/sphere.py
t@@ -4182,8 +4182,16 @@ class sim:
xdisp[i] = xdisp[i-1] + sb.time_file_dt[0] * shearvel
sigma_eff[i] = sb.w_force[0] / A
sigma_def[i] = sb.w_devs[0]
- #dilation[i] = sb.w_x[0] - w_x0 # dilation in meters
- dilation[i] = (sb.w_x[0] - w_x0)/w_x0 * 100.0 # dilation in percent
+
+ # dilation in meters
+ #dilation[i] = sb.w_x[0] - w_x0
+
+ # dilation in percent
+ #dilation[i] = (sb.w_x[0] - w_x0)/w_x0 * 100.0 # dilation in percent
+
+ # dilation in number of mean particle diameters
+ d_bar = numpy.mean(self.radius)*2.0
+ dilation[i] = (sb.w_x[0] - w_x0)/d_var
# Test if this was the max. shear stress
if (tau[i] > tau_u):