treconfigured shear() to not touch upper wall parameters - 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 366c91b40daa1975adce5251c8a51bd3cf3a6ee6
(DIR) parent 64e76c5f86212627f079b05da0ded770660c1dff
(HTM) Author: Anders Damsgaard <adc@geo.au.dk>
Date: Thu, 29 Nov 2012 21:17:39 +0100
reconfigured shear() to not touch upper wall parameters
Diffstat:
M python/sphere.py | 13 ++++++++-----
1 file changed, 8 insertions(+), 5 deletions(-)
---
(DIR) diff --git a/python/sphere.py b/python/sphere.py
t@@ -696,7 +696,7 @@ class Spherebin:
self.w_vel = numpy.array([wvel])
- def shear(self, deviatoric_stress = 10e3,
+ def shear(self,
shear_strain_rate = 1,
periodic = 1):
""" Setup shear experiment. Specify the upper wall
t@@ -705,18 +705,21 @@ class Spherebin:
initial height per second.
"""
- # Compute new grid, scaled to fit max. and min. particle positions
+ # Find lowest and heighest point
z_min = numpy.min(self.x[:,2] - self.radius)
z_max = numpy.max(self.x[:,2] + self.radius)
+
+ # the grid cell size is equal to the max. particle diameter
cellsize = self.L[0] / self.num[0]
- self.adjustUpperWall(z_adjust = 1.3)
- # Initialize upper wall
- self.w_devs[0] = deviatoric_stress
+ # make grid one cell heigher to allow dilation
+ self.num[2] += 1
+ self.L[2] = self.num[2] * cellsize
# zero kinematics
self.zeroKinematics()
+ # set the thickness of the horizons of fixed particles
#fixheight = 2*cellsize
fixheight = cellsize