tupdate shear test example - sphere - GPU-based 3D discrete element method algo… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit f77c6e06f125c71b64eceae7039e285429a7b9cc | |
parent 81d7b44047babc963125026004031b612e6de13d | |
Author: Anders Damsgaard <[email protected]> | |
Date: Fri, 8 Sep 2017 10:29:24 -0400 | |
update shear test example | |
Diffstat: | |
M python/shear-test.py | 27 ++++++++++++++++----------- | |
1 file changed, 16 insertions(+), 11 deletions(-) | |
--- | |
diff --git a/python/shear-test.py b/python/shear-test.py | |
t@@ -14,11 +14,10 @@ plots = True | |
np = 1e4 | |
# Common simulation id | |
-sim_id = "shear-test-devs3" | |
+sim_id = "shear-test" | |
# Deviatoric stress [Pa] | |
-devslist = [80e3, 10e3, 20e3, 40e3, 60e3, 120e3] | |
-#devs = 0 | |
+Nlist = [80e3] | |
### INITIALIZATION ### | |
t@@ -63,12 +62,13 @@ if (initialization == True): | |
# For each normal stress, consolidate and subsequently shear the material | |
-for devs in devslist: | |
+for N in Nlist: | |
### CONSOLIDATION ### | |
# New class | |
- cons = sphere.sim(np = init.np, nw = 1, sid = sim_id + "-cons-devs{}".form… | |
+ cons = sphere.sim(np = init.np, nw = 1, sid = sim_id + | |
+ "-cons-N{}".format(N)) | |
# Read last output file of initialization step | |
lastf = status(sim_id + "-init") | |
t@@ -78,7 +78,8 @@ for devs in devslist: | |
cons.periodicBoundariesXY() | |
# Setup consolidation experiment | |
- cons.consolidate(normal_stress = devs, periodic = init.periodic) | |
+ cons.consolidate(normal_stress = N, periodic = init.periodic) | |
+ cons.adaptiveGrid() | |
# Set duration of simulation | |
t@@ -108,23 +109,27 @@ for devs in devslist: | |
if (rendering == True): | |
# Render images with raytracer | |
- cons.render(method = "pres", max_val = 2.0*devs, verbose = False) | |
+ cons.render(method = "pres", max_val = 2.0*N, verbose = False) | |
### SHEARING ### | |
# New class | |
- shear = sphere.sim(np = cons.np, nw = cons.nw, sid = sim_id + "-shear-devs… | |
+ shear = sphere.sim(np = cons.np, nw = cons.nw, sid = sim_id + | |
+ "-shear-N{}".format(N)) | |
# Read last output file of initialization step | |
- lastf = status(sim_id + "-cons-devs{}".format(devs)) | |
- shear.readbin("../output/" + sim_id + "-cons-devs{}.output{:0=5}.bin".form… | |
+ lastf = status(sim_id + "-cons-N{}".format(N)) | |
+ shear.readbin("../output/" + sim_id + | |
+ "-cons-N{}.output{:0=5}.bin".format(N, lastf), | |
+ verbose = False) | |
# Periodic x and y boundaries | |
shear.periodicBoundariesXY() | |
# Setup shear experiment | |
shear.shear(shear_strain_rate = 0.05, periodic = init.periodic) | |
+ shear.adaptiveGrid() | |
# Set duration of simulation | |
shear.initTemporal(total = 20.0) | |
t@@ -144,4 +149,4 @@ for devs in devslist: | |
if (rendering == True): | |
# Render images with raytracer | |
- shear.render(method = "pres", max_val = 2.0*devs, verbose = False) | |
+ shear.render(method = "pres", max_val = 2.0*N, verbose = False) |