tcall adjustUpperWall in shear function - sphere - GPU-based 3D discrete elemen… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit 7342f87b19e26af65e0f6855c4d4d7161ce0f2bd | |
parent a1f59ee7e5e8915d203ddc5d7203700eb8884dd2 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Mon, 31 Mar 2014 11:05:46 +0200 | |
call adjustUpperWall in shear function | |
Diffstat: | |
M python/sphere.py | 4 ++++ | |
1 file changed, 4 insertions(+), 0 deletions(-) | |
--- | |
diff --git a/python/sphere.py b/python/sphere.py | |
t@@ -2064,6 +2064,7 @@ class sim: | |
#self.w_m[idx] = numpy.array([self.rho[0]*self.np*math.pi \ | |
# *(cellsize/2.0)**3]) | |
self.w_m[idx] = numpy.array([self.rho*self.L[0]*self.L[1]*d_max]) | |
+ print(self.w_m[idx]) | |
def consolidate(self, normal_stress = 10e3): | |
''' | |
t@@ -2165,6 +2166,9 @@ class sim: | |
# zero kinematics | |
self.zeroKinematics() | |
+ # Adjust grid and placement of upper wall | |
+ self.adjustUpperWall() | |
+ | |
# Fix horizontal velocity to 0.0 of lowermost particles | |
d_max_below = numpy.max(self.radius[numpy.nonzero(self.x[:,2] < | |
(z_max-z_min)*0.3)])*2.0 |