 |
tdilation in percent as default - sphere - GPU-based 3D discrete element method… |
 |
git clone git://src.adamsgaard.dk/sphere |
 |
Log |
|
Files |
|
Refs |
|
LICENSE |
|
--- |
|
commit 565130ebc1f63424aff2f72cf769fa7ed75aa7f3 |
|
parent cc45b6bed81acf7b2877065c2229171a9614edb8 |
|
Author: Anders Damsgaard <[email protected]> |
|
Date: Mon, 31 Mar 2014 13:27:38 +0200 |
|
|
|
dilation in percent as default |
|
|
|
Diffstat: |
|
M python/sphere.py | 7 ++++--- |
|
|
|
1 file changed, 4 insertions(+), 3 deletions(-) |
|
--- |
|
diff --git a/python/sphere.py b/python/sphere.py |
|
t@@ -4182,8 +4182,8 @@ 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 p… |
|
+ #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 pe… |
|
|
|
# Test if this was the max. shear stress |
|
if (tau[i] > tau_u): |
|
t@@ -4209,7 +4209,8 @@ class sim: |
|
# Plot dilation |
|
ax2 = plt.subplot2grid((2,1),(1,0)) |
|
ax2.set_xlabel('Shear strain [-]') |
|
- ax2.set_ylabel('Dilation [m]') |
|
+ #ax2.set_ylabel('Dilation [m]') |
|
+ ax2.set_ylabel('Dilation [%]') |
|
ax2.plot(xdisp/w_x0, dilation, '+-') |
|
ax2.grid() |
|
