[cig-commits] r7633 - in short/3D/PyLith/benchmarks/trunk: .
quasistatic quasistatic/strikeslipnog
brad at geodynamics.org
brad at geodynamics.org
Tue Jul 10 10:51:23 PDT 2007
Author: brad
Date: 2007-07-10 10:51:20 -0700 (Tue, 10 Jul 2007)
New Revision: 7633
Added:
short/3D/PyLith/benchmarks/trunk/quasistatic/
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/createbc.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xneg.in.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xneg.spatialdb.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xpos.in.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xpos.spatialdb.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_ypos.in.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_ypos.spatialdb.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_zneg.in.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_zneg.spatialdb.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_zero.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fourmaterials.odb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/genbcpts.py
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/geometry.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0250m.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0500m.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_1000m.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_tapered.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_uniform.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliprate_uniform.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliptime_uniform.spatialdb
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0250m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0250m.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0500m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0500m.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_10000m.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_1000m.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_1000m.exo.gz
Log:
Added files for strike-slip benchmark with hex8 (CUBIT) mesh.
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/README 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,23 @@
+Quasi-static strike-slip (no gravity) benchmark
+
+FILES
+
+pylithapp.cfg - General parameters for problem
+strikeslip_cubit.cfg - Parameters common to CUBIT generated meshes
+strikeslip_hex8_1000m.cfg - Parameters specific to hex8 1000m resolution sim
+strikeslip_hex8_0500m.cfg - Parameters specific to hex8 500m resolution sim
+strikeslip_hex8_0250m.cfg - Parameters specific to hex8 250m resolution sim
+
+TO RUN SIMULATIONS
+
+hex8 @ 1000m resolution
+
+ pylith strikeslip_cubit.cfg strikeslip_hex8_1000m.cfg
+
+hex8 @ 500m resolution
+
+ pylith strikeslip_cubit.cfg strikeslip_hex8_0500m.cfg
+
+hex8 @ 250m resolution
+
+ pylith strikeslip_cubit.cfg strikeslip_hex8_0250m.cfg
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/createbc.jou 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/createbc.jou 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,66 @@
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume elastic_xpos
+block 1 name "elastic x_pos"
+block 2 volume elastic_xneg
+block 2 name "elastic x_neg"
+block 3 volume visco_xpos
+block 3 name "viscoelastic x_pos"
+block 4 volume visco_xneg
+block 4 name "viscoelastic x_neg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" remove node with z_coord < -16.001e+3 or with y_coord > 16.001e+3
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fixed BC
+# ----------------------------------------------------------------------
+group "analytic_disp_xpos" add node in surface 20
+group "analytic_disp_xpos" add node in surface 28
+nodeset 11 group analytic_disp_xpos
+nodeset 11 name "analytic_disp_xpos"
+
+group "analytic_disp_xneg" add node in surface 30
+group "analytic_disp_xneg" add node in surface 38
+nodeset 12 group analytic_disp_xneg
+nodeset 12 name "analytic_disp_xneg"
+
+group "analytic_disp_ypos" add node in surface 21
+group "analytic_disp_ypos" add node in surface 27
+group "analytic_disp_ypos" add node in surface 33
+group "analytic_disp_ypos" add node in surface 35
+group "analytic_disp_ypos" remove node in analytic_disp_xneg
+group "analytic_disp_ypos" remove node in analytic_disp_xpos
+nodeset 13 group analytic_disp_ypos
+nodeset 13 name "analytic_disp_ypos"
+
+group "analytic_disp_zneg" add node in surface 12
+group "analytic_disp_zneg" add node in surface 16
+group "analytic_disp_zneg" remove node in analytic_disp_xneg
+group "analytic_disp_zneg" remove node in analytic_disp_xpos
+group "analytic_disp_zneg" remove node in analytic_disp_ypos
+group "analytic_disp_zneg" remove node with z_coord < -23.999e+3 and with x_coord > 11.999e+3 and with x_coord < 12.001e+3
+nodeset 14 group analytic_disp_zneg
+nodeset 14 name "analytic_disp_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for symmetry BC
+# ----------------------------------------------------------------------
+group "symmetry_disp" add node in surface 23
+group "symmetry_disp" add node in surface 25
+group "symmetry_disp" add node in surface 31
+group "symmetry_disp" add node in surface 37
+group "symmetry_disp" remove node in analytic_disp_xpos
+group "symmetry_disp" remove node in analytic_disp_xneg
+group "symmetry_disp" remove node in analytic_disp_ypos
+group "symmetry_disp" remove node in analytic_disp_zneg
+group "symmetry_disp" remove node in fault
+nodeset 15 group symmetry_disp
+nodeset 15 name "symmetry_disp"
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xneg.in.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xneg.in.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xneg.spatialdb.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xneg.spatialdb.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xpos.in.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xpos.in.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xpos.spatialdb.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_xpos.spatialdb.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_ypos.in.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_ypos.in.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_ypos.spatialdb.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_ypos.spatialdb.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_zneg.in.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_zneg.in.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_zneg.spatialdb.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_analytic_zneg.spatialdb.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_zero.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_zero.spatialdb 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fixeddisp_zero.spatialdb 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = dof-0 dof-1 dof-2
+ value-units = m m m
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 0.0 0.0 0.0
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fourmaterials.odb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fourmaterials.odb 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/fourmaterials.odb 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,90 @@
+#!/usr/bin/env python
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard
+# U.S. Geological Survey
+#
+# <LicenseText>
+#
+# ----------------------------------------------------------------------
+#
+
+## Container of four materials.
+
+# Parent class
+from pylith.utils.ObjectBin import ObjectBin
+
+class MyMatBin(ObjectBin):
+ """
+ User-defined materials container.
+
+ Factory: object_bin
+ """
+
+ # INVENTORY //////////////////////////////////////////////////////////
+
+ class Inventory(ObjectBin.Inventory):
+ """
+ Python object for managing MyBC facilities and properties.
+ """
+
+ import pyre.inventory
+
+ from pylith.materials.ElasticIsotropic3D import ElasticIsotropic3D
+ from pylith.materials.MaxwellIsotropic3D import MaxwellIsotropic3D
+
+ elasticNeg = pyre.inventory.facility("elastic_neg", family="material",
+ factory=ElasticIsotropic3D)
+ elasticNeg.meta['tip'] = "Elastic material on -x side of the fault."
+
+ elasticPos = pyre.inventory.facility("elastic_pos", family="material",
+ factory=ElasticIsotropic3D)
+ elasticPos.meta['tip'] = "Elastic material on +x side of the fault."
+
+ viscoelasticNeg = pyre.inventory.facility("viscoelastic_neg",
+ family="material",
+ factory=MaxwellIsotropic3D)
+ viscoelasticNeg.meta['tip'] = "Viscoelastic material on -x side of the fault."
+
+ viscoelasticPos = pyre.inventory.facility("viscoelastic_pos",
+ family="material",
+ factory=MaxwellIsotropic3D)
+ viscoelasticPos.meta['tip'] = "Viscoelastic material on +x side of the fault."
+
+
+ # PUBLIC METHODS /////////////////////////////////////////////////////
+
+ def __init__(self, name="mymatbin"):
+ """
+ Constructor.
+ """
+ ObjectBin.__init__(self, name)
+ return
+
+
+ # PRIVATE METHODS ////////////////////////////////////////////////////
+
+ def _configure(self):
+ """
+ Set attributes from inventory.
+ """
+ ObjectBin._configure(self)
+
+ self.bin = [self.inventory.elasticNeg,
+ self.inventory.elasticPos,
+ self.inventory.viscoelasticNeg,
+ self.inventory.viscoelasticPos]
+ return
+
+
+# FACTORIES ////////////////////////////////////////////////////////////
+
+def object_bin():
+ """
+ Factory associated with MyMatBin.
+ """
+ return MyMatBin()
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/genbcpts.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/genbcpts.py 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/genbcpts.py 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,48 @@
+#!/usr/bin/env python
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard
+# U.S. Geological Survey
+#
+# <LicenseText>
+#
+# ----------------------------------------------------------------------
+#
+
+dx = 250.0
+faces = [
+ {'label': "ypos",
+ 'xrange': [0.0, 24.0e+3],
+ 'yrange': [24.0e+3, 24.0e+3],
+ 'zrange': [-24e+3, 0.0]},
+ {'label': "xpos",
+ 'xrange': [24.0e+3, 24.0e+3],
+ 'yrange': [0.0, 24.0e+3],
+ 'zrange': [-24.0e+3, 0.0]},
+ {'label': "xneg",
+ 'xrange': [0.0, 0.0],
+ 'yrange': [0.0, 24.0e+3],
+ 'zrange': [-24.0e+3, 0.0]},
+ {'label': "zneg",
+ 'xrange': [0.0, 24.0e+3],
+ 'yrange': [0.0, 24.0e+3],
+ 'zrange': [-24.0e+3, -24.0e+3]}
+ ]
+
+for face in faces:
+ fout = open("bcpts_%s.in" % face['label'], "w")
+ numX = int(1 + (face['xrange'][1] - face['xrange'][0])/dx)
+ numY = int(1 + (face['yrange'][1] - face['yrange'][0])/dx)
+ numZ = int(1 + (face['zrange'][1] - face['zrange'][0])/dx)
+ for iX in xrange(numX):
+ x = face['xrange'][0] + iX*dx
+ for iY in xrange(numY):
+ y = face['yrange'][0] + iY*dx
+ for iZ in xrange(numZ):
+ z = face['zrange'][0] + iZ*dx
+ fout.write("%12.4e %12.4e %12.4e\n" % (x, y, z))
+ fout.close()
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/geometry.jou 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/geometry.jou 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,31 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+reset
+brick x 24000 y 24000 z 24000
+volume 1 move x 12000 y 12000 z -12000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 12000
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -12000
+surface 8 name "material_interface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interfaces surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+webcut volume 1 with plane surface material_interface
+webcut volume 4 with plane surface material_interface
+volume 1 name "elastic_xpos"
+volume 4 name "elastic_xneg"
+volume 5 name "visco_xpos"
+volume 6 name "visco_xneg"
+
+# ----------------------------------------------------------------------
+# Imprint curves on volumes
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mat_elastic.spatialdb 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = density vs vp
+ value-units = kg/m^3 m/s m/s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2500.0 3000.0 5291.502622129181
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0250m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0250m.jou 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0250m.jou 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,24 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 250.0
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "strikeslip_hex8_0250m.exo" dimension 3 overwrite
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0500m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0500m.jou 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_0500m.jou 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,24 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 500.0
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "strikeslip_hex8_0500m.exo" dimension 3 overwrite
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_1000m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_1000m.jou 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/mesh_hex8_1000m.jou 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,24 @@
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 1000
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "strikeslip_hex8_1000m.exo" dimension 3 overwrite
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/pylithapp.cfg 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,145 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+debug = 0
+
+[pylithapp.mesh_generator.importer]
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.0*s
+default_dt = 1.0*s
+dimension = 3
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCSixSides
+interfaces = pylith.faults.SingleFault
+materials = fourmaterials
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+# TEMPORARY
+viscoelastic_neg = pylith.materials.ElasticIsotropic3D
+viscoelastic_pos = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.elastic_pos]
+label = Elastic material +x
+db.iohandler.filename = mat_elastic.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+[pylithapp.timedependent.materials.elastic_neg]
+label = Elastic material -x
+db.iohandler.filename = mat_elastic.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+[pylithapp.timedependent.materials.viscoelastic_pos]
+label = Elastic material +x
+db.iohandler.filename = mat_elastic.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+[pylithapp.timedependent.materials.viscoelastic_neg]
+label = Viscoelastic material -x
+db.iohandler.filename = mat_elastic.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.y_neg]
+fixed_dof = [0, 2]
+db.label = Symmetry Dirichlet BC
+db.iohandler.filename = fixeddisp_zero.spatialdb
+db.query_type = nearest
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1, 2]
+db.label = Analytic displacement Dirichlet BC on +x face
+db.iohandler.filename = fixeddisp_analytic_xpos.spatialdb
+db.query_type = nearest
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1, 2]
+db.label = Analytic displacement Dirichlet BC on -x face
+db.iohandler.filename = fixeddisp_analytic_xneg.spatialdb
+db.query_type = nearest
+
+[pylithapp.timedependent.bc.y_pos]
+fixed_dof = [0, 1, 2]
+db.label = Analytic displacement Dirichlet BC on +y face
+db.iohandler.filename = fixeddisp_analytic_ypos.spatialdb
+db.query_type = nearest
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [0, 1, 2]
+db.label = Analytic displacement Dirichlet BC on -z face
+db.iohandler.filename = fixeddisp_analytic_zneg.spatialdb
+db.query_type = nearest
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault]
+id = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+[pylithapp.timedependent.interfaces.fault.eq_src.slip_function]
+slip.label = Final slip
+slip.iohandler.filename = slip_tapered.spatialdb
+slip.query_type = linear
+slip_rate.label = Peak slip rate
+slip_rate.iohandler.filename = sliprate_uniform.spatialdb
+slip_time.label = Slip initiation time
+slip_time.iohandler.filename = sliptime_uniform.spatialdb
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+
+pc_type = asm
+ksp_rtol = 1.0e-8
+ksp_gmres_restart = 5000
+
+#ksp_type = preonly
+#pc_type = lu
+#mat_type = aijmumps
+
+ksp_max_it = 500
+ksp_monitor = true
+ksp_view = true
+log_summary = true
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_tapered.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_tapered.spatialdb 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_tapered.spatialdb 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,54 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+ value-units = m m m
+ num-locs = 24
+ data-dim = 2
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 3
+ }
+}
+// y = 0
+0.0 0.0 0.0 -1.0 0.0 0.0
+0.0 0.0 -12.0 -1.0 0.0 0.0
+0.0 0.0 -16.0 0.0 0.0 0.0
+
+// y = 12
+0.0 12.0 0.0 -1.0 0.0 0.0
+0.0 12.0 -12.0 -1.0 0.0 0.0
+0.0 12.0 -16.0 0.0 0.0 0.0
+
+// y = 16
+0.0 16.0 0.0 0.0 0.0 0.0
+0.0 16.0 -12.0 0.0 0.0 0.0
+0.0 16.0 -16.0 0.0 0.0 0.0
+
+// y = 13
+0.0 13.0 -13.0 -0.75 0.0 0.0
+0.0 13.0 -14.0 -0.50 0.0 0.0
+0.0 13.0 -15.0 -0.25 0.0 0.0
+0.0 13.0 -16.0 0.0 0.0 0.0
+
+// y = 14
+0.0 14.0 -14.0 -0.50 0.0 0.0
+0.0 14.0 -15.0 -0.25 0.0 0.0
+0.0 14.0 -16.0 0.0 0.0 0.0
+
+// y = 15
+0.0 15.0 -15.0 -0.25 0.0 0.0
+0.0 15.0 -16.0 0.0 0.0 0.0
+
+// z = -13
+0.0 14.0 -13.0 -0.50 0.0 0.0
+0.0 15.0 -13.0 -0.25 0.0 0.0
+0.0 16.0 -13.0 0.0 0.0 0.0
+
+// z = -14
+0.0 15.0 -14.0 -0.25 0.0 0.0
+0.0 16.0 -14.0 0.0 0.0 0.0
+
+// z = -15
+0.0 15.0 -15.0 -0.25 0.0 0.0
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_uniform.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_uniform.spatialdb 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/slip_uniform.spatialdb 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+ value-units = m m m
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 -1.0 0.0 0.0
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliprate_uniform.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliprate_uniform.spatialdb 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliprate_uniform.spatialdb 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 1
+ value-names = slip-rate
+ value-units = m/s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 1.0e+6
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliptime_uniform.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliptime_uniform.spatialdb 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/sliptime_uniform.spatialdb 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 1
+ value-names = slip-time
+ value-units = s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 -1.0
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_cubit.cfg 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,47 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+importer = pylith.meshio.MeshIOCubit
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic_pos]
+id = 1
+
+[pylithapp.timedependent.materials.elastic_neg]
+id = 2
+
+[pylithapp.timedependent.materials.viscoelastic_pos]
+id = 3
+
+[pylithapp.timedependent.materials.viscoelastic_neg]
+id = 4
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.y_neg]
+label = 15
+
+[pylithapp.timedependent.bc.x_pos]
+label = 11
+
+[pylithapp.timedependent.bc.x_neg]
+label = 12
+
+[pylithapp.timedependent.bc.y_pos]
+label = 13
+
+[pylithapp.timedependent.bc.z_neg]
+label = 14
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+label = 10
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0250m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0250m.cfg 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0250m.cfg 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.importer]
+filename = strikeslip_hex8_0250m.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+filename = strikeslip_hex8_0250m.vtk
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0250m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0250m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0500m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0500m.cfg 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0500m.cfg 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.importer]
+filename = strikeslip_hex8_0500m.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+filename = strikeslip_hex8_0500m.vtk
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0500m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_0500m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_10000m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_10000m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_1000m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_1000m.cfg 2007-07-10 17:49:28 UTC (rev 7632)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_1000m.cfg 2007-07-10 17:51:20 UTC (rev 7633)
@@ -0,0 +1,14 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.importer]
+filename = strikeslip_hex8_1000m.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+filename = strikeslip_hex8_1000m.vtk
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_1000m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/strikeslip_hex8_1000m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
More information about the cig-commits
mailing list