trun consolidation with heavy wall instead of a normal stress value - 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 123496aa3901d83e7f8084934e3b23867383eff3
(DIR) parent 92973723eddc4dc1b88d2e40ae0338cd0ba96ed9
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 16 Sep 2014 09:19:18 +0200
run consolidation with heavy wall instead of a normal stress value
Diffstat:
M python/sphere.py | 7 ++++++-
1 file changed, 6 insertions(+), 1 deletion(-)
---
(DIR) diff --git a/python/sphere.py b/python/sphere.py
t@@ -2391,7 +2391,12 @@ class sim:
# Set the top wall BC to a value of normal stress
self.wmode = numpy.array([1])
- self.w_devs = numpy.ones(1) * normal_stress
+ #self.w_devs = numpy.ones(1) * normal_stress
+
+ # Set top wall to a certain mass corresponding to the selected normal
+ # stress
+ self.w_devs = numpy.zeros(1)
+ self.w_m[0] = numpy.abs(normal_stress*sim.L[0]*sim.L[1]/sim.g[2])
def uniaxialStrainRate(self, wvel = -0.001):
'''