tdimensionless dilation - sphere - GPU-based 3D discrete element method algorit… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit 379ce413f9d61b1688b740e177cfda5d605f0285 | |
parent edf3cc179f0b80edbf759132bd207792f5708b8f | |
Author: Anders Damsgaard <[email protected]> | |
Date: Tue, 9 Sep 2014 10:22:53 +0200 | |
dimensionless dilation | |
Diffstat: | |
M python/shear-results.py | 6 ++++-- | |
1 file changed, 4 insertions(+), 2 deletions(-) | |
--- | |
diff --git a/python/shear-results.py b/python/shear-results.py | |
t@@ -69,16 +69,18 @@ ax2 = plt.subplot(212, sharex=ax1) | |
ax1.plot(shear_strain[0], friction[0], label='dry') | |
ax2.plot(shear_strain[0], dilation[0], label='dry') | |
+mean_diameter = numpy.mean(sim.radius)*2.0 | |
for c in numpy.arange(1,len(cvals)+1): | |
ax1.plot(shear_strain[c], friction[c], label='$c$ = %.2f' % (cvals[c-1])) | |
- ax2.plot(shear_strain[c], dilation[c], label='$c$ = %.2f' % (cvals[c-1])) | |
+ ax2.plot(shear_strain[c], dilation[c]/mean_diameter, \ | |
+ label='$c$ = %.2f' % (cvals[c-1])) | |
#plt.plot(dpdz[c], K[c], 'o-', label='$c$ = %.2f' % (cvals[c])) | |
#plt.semilogx(dpdz[c], K[c], 'o-', label='$c$ = %.2f' % (cvals[c])) | |
#plt.semilogy(dpdz[c], K[c], 'o-', label='$c$ = %.2f' % (cvals[c])) | |
#plt.loglog(dpdz[c], K[c], 'o-', label='$c$ = %.2f' % (cvals[c])) | |
ax2.set_xlabel('Shear strain [-]') | |
ax1.set_ylabel('Shear friction $\\tau/\\sigma\'$ [-]') | |
-ax2.set_ylabel('Dilation $\\Delta h$ [m]') | |
+ax2.set_ylabel('Dilation $\\Delta h/(2r)$ [-]') | |
plt.setp(ax1.get_xticklabels(), visible=False) | |
ax1.grid() | |
ax2.grid() |