tAdded shear-test - 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 3c4f2087dab043e3e6bf921b8022150c8f78de51
 (DIR) parent 2a487de72d1a16ff8631293498605819772b8f15
 (HTM) Author: Anders Damsgaard <adc@geo.au.dk>
       Date:   Wed, 19 Dec 2012 11:52:07 +0100
       
       Added shear-test
       
       Diffstat:
         A python/shear-test.py                |     128 +++++++++++++++++++++++++++++++
       
       1 file changed, 128 insertions(+), 0 deletions(-)
       ---
 (DIR) diff --git a/python/shear-test.py b/python/shear-test.py
       t@@ -0,0 +1,128 @@
       +#!/usr/bin/env python
       +
       +# Import sphere functionality
       +from sphere import *
       +
       +### EXPERIMENT SETUP ###
       +initialization = False
       +consolidation  = False
       +shearing       = True
       +rendering      = True
       +plots               = True
       +
       +# Number of particles
       +np = 2e3
       +
       +# Common simulation id
       +sim_id = "shear-test"
       +
       +# Deviatoric stress [Pa]
       +devs = 10e3
       +#devs = 0
       +
       +### INITIALIZATION ###
       +
       +# New class
       +init = Spherebin(np = np, nd = 3, nw = 0, sid = sim_id + "-init")
       +
       +# Save radii
       +init.generateRadii(radius_mean = 0.05)
       +
       +# Use default params
       +init.defaultParams(gamma_n = 0.0, mu_s = 0.4, mu_d = 0.4)
       +
       +# Initialize positions in random grid (also sets world size)
       +init.initRandomGridPos(gridnum = numpy.array([12, 12, 1000]), periodic = 1, contactmodel = 2)
       +
       +# Set duration of simulation
       +init.initTemporal(total = 5.0)
       +
       +if (initialization == True):
       +  # Write input file for sphere
       +  init.writebin()
       +
       +  # Run sphere
       +  init.run()
       +
       +  if (plots == True):
       +    # Make a graph of energies
       +    visualize(init.sid, "energy", savefig=True, outformat='png')
       +
       +  if (rendering == True):
       +    # Render images with raytracer
       +    init.render(method = "angvel", max_val = 0.3, verbose = False)
       +
       +
       +### CONSOLIDATION ###
       +
       +# New class
       +cons = Spherebin(np = init.np, nw = 1, sid = sim_id + "-cons")
       +
       +# Read last output file of initialization step
       +lastf = status(sim_id + "-init")
       +cons.readbin("../output/" + sim_id + "-init.output{:0=5}.bin".format(lastf), verbose=False)
       +
       +# Setup consolidation experiment
       +cons.consolidate(deviatoric_stress = devs, periodic = init.periodic)
       +
       +
       +# Set duration of simulation
       +cons.initTemporal(total = 1.5)
       +#cons.initTemporal(total = 0.0019, file_dt = 0.00009)
       +#cons.initTemporal(total = 0.0019, file_dt = 1e-6)
       +#cons.initTemporal(total = 0.19, file_dt = 0.019)
       +
       +cons.w_m[0] *= 0.001
       +
       +
       +
       +if (consolidation == True):
       +  # Write input file for sphere
       +  cons.writebin()
       +
       +  # Run sphere
       +  cons.run(dry=True) # show values, don't run
       +  cons.run() # run
       +
       +  if (plots == True):
       +    # Make a graph of energies
       +    visualize(cons.sid, "energy", savefig=True, outformat='png')
       +    visualize(cons.sid, "walls", savefig=True, outformat='png')
       +
       +  if (rendering == True):
       +    # Render images with raytracer
       +    cons.render(method = "pres", max_val = 2.0*devs, verbose = False)
       +
       +
       +### SHEARING ###
       +
       +# New class
       +shear = Spherebin(np = cons.np, nw = cons.nw, sid = sim_id + "-shear")
       +
       +# Read last output file of initialization step
       +lastf = status(sim_id + "-cons")
       +shear.readbin("../output/" + sim_id + "-cons.output{:0=5}.bin".format(lastf), verbose = False)
       +
       +# Setup shear experiment
       +shear.shear(shear_strain_rate = 0.10, periodic = init.periodic)
       +
       +# Set duration of simulation
       +shear.initTemporal(total = 5.0)
       +
       +if (shearing == True):
       +  # Write input file for sphere
       +  shear.writebin()
       +
       +  # Run sphere
       +  shear.run(dry=True)
       +  shear.run()
       +
       +  if (plots == True):
       +    # Make a graph of energies
       +    visualize(shear.sid, "energy", savefig=True, outformat='png')
       +    visualize(shear.sid, "shear", savefig=True, outformat='png')
       +
       +  if (rendering == True):
       +    # Render images with raytracer
       +    shear.render(method = "pres", max_val = 2.0*devs, verbose = False)
       +