[cig-commits] r17129 - in short/3D/PyLith/branches/v1.5-stable: . tests/3d tests/3d/plasticity tests/3d/plasticity/threehex27 tests/3d/plasticity/threehex27/config tests/3d/plasticity/threehex27/mesh tests/3d/plasticity/threehex27/results tests/3d/plasticity/threehex27/spatialdb tests/3d/plasticity/threehex8 tests/3d/plasticity/threehex8/config tests/3d/plasticity/threehex8/mesh tests/3d/plasticity/threehex8/results tests/3d/plasticity/threehex8/spatialdb
willic3 at geodynamics.org
willic3 at geodynamics.org
Thu Aug 26 21:35:48 PDT 2010
Author: willic3
Date: 2010-08-26 21:35:47 -0700 (Thu, 26 Aug 2010)
New Revision: 17129
Added:
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/geometry.jou
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/threehex27.exo
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/results/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/results/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_bc_3.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_0015.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_012.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_03.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/initial_stress.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_elastic.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp1.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp2.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp3.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_vm1.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e01.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e02.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e03.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e04.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e05.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e06.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e07.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e08.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e09.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e10.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e11.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e12.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/elas_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp1_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp2_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp3_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_vm1_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/pylithapp.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/visc_e10_ng_ni_v03.cfg
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/geometry.jou
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/mesh_threehex8.jou
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/threehex8.exo
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/results/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/results/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/Makefile.am
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_bc_3.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_0015.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_012.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_03.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/initial_stress.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_elastic.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp1.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp2.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp3.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_vm1.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e01.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e02.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e03.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e04.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e05.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e06.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e07.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e08.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e09.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e10.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e11.spatialdb
short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e12.spatialdb
Modified:
short/3D/PyLith/branches/v1.5-stable/configure.ac
short/3D/PyLith/branches/v1.5-stable/tests/3d/Makefile.am
Log:
Added simple 3-cell tests for plasticity.
Both linear and quadratic hexes are being used.
Tests are incomplete. I still need to run and debug them all, and include
gravity and initial stresses.
Modified: short/3D/PyLith/branches/v1.5-stable/configure.ac
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/configure.ac 2010-08-27 01:00:49 UTC (rev 17128)
+++ short/3D/PyLith/branches/v1.5-stable/configure.ac 2010-08-27 04:35:47 UTC (rev 17129)
@@ -304,6 +304,17 @@
tests/3d/Makefile
tests/3d/matprops/Makefile
tests/3d/slipdir/Makefile
+ tests/3d/plasticity/Makefile
+ tests/3d/plasticity/threehex27/Makefile
+ tests/3d/plasticity/threehex27/config/Makefile
+ tests/3d/plasticity/threehex27/mesh/Makefile
+ tests/3d/plasticity/threehex27/results/Makefile
+ tests/3d/plasticity/threehex27/spatialdb/Makefile
+ tests/3d/plasticity/threehex8/Makefile
+ tests/3d/plasticity/threehex8/config/Makefile
+ tests/3d/plasticity/threehex8/mesh/Makefile
+ tests/3d/plasticity/threehex8/results/Makefile
+ tests/3d/plasticity/threehex8/spatialdb/Makefile
doc/Makefile
doc/developer/Makefile
doc/install/Makefile
Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/Makefile.am 2010-08-27 01:00:49 UTC (rev 17128)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -18,6 +18,7 @@
SUBDIRS = \
matprops \
+ plasticity \
slipdir
# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,23 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+SUBDIRS = \
+ threehex8 \
+ threehex27
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,25 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+SUBDIRS = \
+ config \
+ mesh \
+ results \
+ spatialdb
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA = \
+ elas_ng_ni_v03.cfg \
+ plas_dp1_ng_ni_v03.cfg \
+ plas_dp2_ng_ni_v03.cfg \
+ plas_dp3_ng_ni_v03.cfg \
+ plas_vm1_ng_ni_v03.cfg \
+ pylithapp.cfg \
+ visc_e10_ng_ni_v03.cfg
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,121 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is elastic.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The linear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith elas_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output managers.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic_material]
+
+[pylithapp.timedependent.materials.elastic_material]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_elastic.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda]
+output.cell_data_fields = [total_strain,stress]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/elas_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/elas_ng_ni_v03-elastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,126 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Drucker-Prager elastoplastic with a friction angle
+# of 1 deg, a cohesion of 50 Pa, and a dilatation angle of 0 deg.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_dp1_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp1.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_dp1_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_dp1_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,126 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Drucker-Prager elastoplastic with a friction angle
+# of 5 deg, a cohesion of 50 Pa, and a dilatation angle of 0 deg.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_dp2_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp2.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_dp2_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_dp2_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,126 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Drucker-Prager elastoplastic with a friction angle
+# of 10 deg, a cohesion of 50 Pa, and a dilatation angle of 0 deg.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_dp3_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp3.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_dp3_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_dp3_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,125 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is von Mises elastoplastic with a cohesion of 50 Pa.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_vm1_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_vm1.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_vm1_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_vm1_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,105 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies the general parameters common to the simulations
+# in this directory.
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+faultcohesivekin = 1
+fiatlagrange = 1
+pylithapp = 1
+materials = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1 ; uncomment to get very verbose mesh information
+
+# Change the default mesh reader to the CUBIT reader.
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+# Set filename of mesh to import.
+filename = ../mesh/threehex8.exo
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Nondimensionalization
+normalizer.length_scale = 1.0*m
+normalizer.shear_modulus = 4.0e5*Pa
+normalizer.relaxation_time = 1.0*s
+
+[pylithapp.timedependent.formulation.time_step]
+# Define the total time for the simulation and the default time step size.
+total_time = 0.0*s ; total time of simulation
+
+# ----------------------------------------------------------------------
+# PETSc and solution settings.
+# ----------------------------------------------------------------------
+# Set the solver options.
+
+[pylithapp.petsc]
+
+# Preconditioner settings.
+# ksp_type = gmres
+pc_type = asm
+# comment this out for fieldSplit
+sub_pc_factor_shift_type = nonzero
+
+# Convergence parameters.
+ksp_rtol = 1.0e-10
+ksp_atol = 1.0e-22
+ksp_gmres_restart = 150
+ksp_max_it = 1500
+
+# Linear solver monitoring options.
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+ksp_monitor_singular_value = true
+log_summary = true
+
+# Field split
+# [pylithapp.timedependent.formulation]
+# split_fields = True
+# use_custom_constraint_pc = True
+# matrix_type = aij
+
+[pylithapp.petsc]
+# fs_pc_type = fieldsplit
+# fs_pc_fieldsplit_real_diagonal =
+# fs_pc_fieldsplit_type = multiplicative
+# fs_fieldsplit_0_pc_type = ml
+# fs_fieldsplit_1_pc_type = ml
+# fs_fieldsplit_2_pc_type = ml
+# fs_fieldsplit_3_pc_type = ml
+# fs_fieldsplit_0_ksp_type = preonly
+# fs_fieldsplit_1_ksp_type = preonly
+# fs_fieldsplit_2_ksp_type = preonly
+# fs_fieldsplit_3_ksp_type = preonly
+
+# Nonlinear solver monitoring options.
+snes_rtol = 1.0e-12
+snes_atol = 1.0e-22
+snes_max_it = 25
+snes_monitor = true
+snes_view = true
+snes_converged_reason = true
+
+# Uncomment to launch gdb when starting PyLith.
+# start_in_debugger = true
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,128 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Maxwell viscoelastic with a viscosity of
+# 1.3e10 Pa-s.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The linear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith visc_e10_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Use adaptive time stepping.
+[pylithapp.timedependent.formulation]
+time_step = pylith.problems.TimeStepAdapt
+
+# Change the total simulation time to 10 years, and set the maximum time
+# step size to 0.1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+max_dt = 0.1*year
+stability_factor = 0.1
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [viscous_material]
+materials.viscous_material = pylith.materials.MaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.viscous_material]
+label = Viscous material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_viscous_e10.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,maxwell_time]
+output.cell_data_fields = [total_strain,stress,viscous_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/visc_e10_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.viscous_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/visc_e10_ng_ni_v03-viscous_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,24 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA = \
+ geometry.jou \
+ mesh_threehex27.jou \
+ threehex27.exo
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/geometry.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/geometry.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/geometry.jou 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,26 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 3 m x 1 m x 1 m
+# -1.5 m <= x <= 1.5 m
+# -0.5 m <= y <= 0.5 m
+# -0.5 m <= z <= 0.5 m
+reset
+brick x 3 y 1 z 1
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,82 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 1
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Create block for material and define cells as HEX27
+# ----------------------------------------------------------------------
+block 1 volume 1
+block 1 name "plastic"
+block 1 element type HEX27
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 6
+# We remove the nodes in face_ypos so there is no conflict in BC.
+group "face_xpos" remove node in surface 3
+nodeset 11 group face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 4
+# We remove the nodes in face_yneg so there is no conflict in BC.
+group "face_xneg" remove node in surface 3
+nodeset 12 group face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 5
+nodeset 13 group face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 3
+nodeset 14 group face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 2
+nodeset 15 group face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 1
+nodeset 16 group face_zpos
+nodeset 16 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "threehex27.exo" dimension 3 overwrite
+
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/threehex27.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/threehex27.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/results/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/results/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/results/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,22 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA =
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,43 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA = \
+ basal_bc_3.spatialdb \
+ basal_rate_bc_0015.spatialdb \
+ basal_rate_bc_012.spatialdb \
+ basal_rate_bc_03.spatialdb \
+ initial_stress.spatialdb \
+ mat_elastic.spatialdb \
+ mat_plastic_dp1.spatialdb \
+ mat_plastic_dp2.spatialdb \
+ mat_plastic_dp3.spatialdb \
+ mat_plastic_vm1.spatialdb \
+ mat_viscous_e01.spatialdb \
+ mat_viscous_e02.spatialdb \
+ mat_viscous_e03.spatialdb \
+ mat_viscous_e04.spatialdb \
+ mat_viscous_e05.spatialdb \
+ mat_viscous_e06.spatialdb \
+ mat_viscous_e07.spatialdb \
+ mat_viscous_e08.spatialdb \
+ mat_viscous_e09.spatialdb \
+ mat_viscous_e10.spatialdb \
+ mat_viscous_e11.spatialdb \
+ mat_viscous_e12.spatialdb
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_bc_3.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_bc_3.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_bc_3.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = displacement-x displacement-y displacement-z
+ value-units = mm mm mm
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Ux (mm)
+// (5) Uy (mm)
+// (6) Uz (mm)
+-1.5000 0.0 0.0 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.3 0.0 0.0
+ 1.5000 0.0 0.0 -0.3 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_0015.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_0015.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_0015.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = displacement-rate-x displacement-rate-y displacement-rate-z rate-start-time
+ value-units = mm/year mm/year mm/year year
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Vx (mm/year)
+// (5) Vy (mm/year)
+// (6) Vz (mm/year)
+// (7) t0 (year)
+-1.5000 0.0 0.0 0.0000 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.0000 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.0015 0.0 0.0 0.0
+ 1.5000 0.0 0.0 -0.0015 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_012.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_012.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_012.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = displacement-rate-x displacement-rate-y displacement-rate-z rate-start-time
+ value-units = mm/year mm/year mm/year year
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Vx (mm/year)
+// (5) Vy (mm/year)
+// (6) Vz (mm/year)
+// (7) t0 (year)
+-1.5000 0.0 0.0 0.000 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.000 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.012 0.0 0.0 0.0
+ 1.5000 0.0 0.0 -0.012 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_03.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_03.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/basal_rate_bc_03.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = displacement-rate-x displacement-rate-y displacement-rate-z rate-start-time
+ value-units = mm/year mm/year mm/year year
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Vx (mm/year)
+// (5) Vy (mm/year)
+// (6) Vz (mm/year)
+// (7) t0 (year)
+-1.5000 0.0 0.0 0.000 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.000 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.030 0.0 0.0 0.0
+ 1.5000 0.0 0.0 -0.030 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/initial_stress.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/initial_stress.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/initial_stress.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,36 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the initial stresses for the model
+// The stresses increase linearly with depth.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of stress components
+
+ // names of the stress components
+ value-names = stress-xx stress-yy stress-zz stress-xy stress-yz stress-xz
+
+ value-units = Pa Pa Pa Pa Pa Pa // units
+
+ num-locs = 2 // number of locations: 1 value each at the upper and lower
+ // surfaces. Linear interpolation should be used.
+
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) stress-xx (Pa)
+// (5) stress-yy (Pa)
+// (6) stress-zz (Pa)
+// (7) stress-xy (Pa)
+// (8) stress-yz (Pa)
+// (9) stress-xz (Pa)
+0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+0.0 -1.0 0.0 -100.0 -100.0 -100.0 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_elastic.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_elastic.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,27 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for an elastic material. In this case, the material
+// properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3 // number of material property values
+ value-names = density vs vp // names of the material property values
+ value-units = kg/m**3 m/s m/s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp1.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp1.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp1.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp2.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp2.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp2.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 5.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp3.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp3.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_dp3.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 10.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_vm1.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_vm1.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_plastic_vm1.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 0.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e01.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e01.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e01.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e1
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e02.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e02.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e02.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e2
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e03.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e03.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e03.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e3
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e04.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e04.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e04.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e4
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e05.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e05.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e05.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e5
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e06.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e06.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e06.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e6
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e07.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e07.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e07.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e7
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e08.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e08.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e08.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e8
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e09.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e09.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e09.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e9
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e10.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e10.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e10.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e10
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e11.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e11.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e11.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e11
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e12.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e12.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/spatialdb/mat_viscous_e12.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e12
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,25 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+SUBDIRS = \
+ config \
+ mesh \
+ results \
+ spatialdb
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA = \
+ elas_ng_ni_v03.cfg \
+ plas_dp1_ng_ni_v03.cfg \
+ plas_dp2_ng_ni_v03.cfg \
+ plas_dp3_ng_ni_v03.cfg \
+ plas_vm1_ng_ni_v03.cfg \
+ pylithapp.cfg \
+ visc_e10_ng_ni_v03.cfg
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/elas_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/elas_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/elas_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,121 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is elastic.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The linear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith elas_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output managers.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic_material]
+
+[pylithapp.timedependent.materials.elastic_material]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_elastic.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda]
+output.cell_data_fields = [total_strain,stress]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/elas_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/elas_ng_ni_v03-elastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp1_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp1_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp1_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,126 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Drucker-Prager elastoplastic with a friction angle
+# of 1 deg, a cohesion of 50 Pa, and a dilatation angle of 0 deg.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_dp1_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp1.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_dp1_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_dp1_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp2_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp2_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp2_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,126 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Drucker-Prager elastoplastic with a friction angle
+# of 5 deg, a cohesion of 50 Pa, and a dilatation angle of 0 deg.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_dp2_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp2.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_dp2_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_dp2_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp3_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp3_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_dp3_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,126 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Drucker-Prager elastoplastic with a friction angle
+# of 10 deg, a cohesion of 50 Pa, and a dilatation angle of 0 deg.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_dp3_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp3.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_dp3_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_dp3_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_vm1_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_vm1_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/plas_vm1_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,125 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is von Mises elastoplastic with a cohesion of 50 Pa.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The nonlinear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith plas_vm1_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# For this problem we must switch to a nonlinear solver.
+formulation.solver = pylith.problems.SolverNonlinear
+
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Change the total simulation time to 10 years, and set the time
+# step size to 1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [plastic_material]
+materials.plastic_material = pylith.materials.DruckerPrager3D
+
+[pylithapp.timedependent.materials.plastic_material]
+label = Plastic material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_plastic_vm1.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+output.cell_data_fields = [total_strain,stress,plastic_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/plas_vm1_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.plastic_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/plas_vm1_ng_ni_v03-plastic_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/pylithapp.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,105 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies the general parameters common to the simulations
+# in this directory.
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshiocubit = 1
+implicitelasticity = 1
+faultcohesivekin = 1
+fiatlagrange = 1
+pylithapp = 1
+materials = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1 ; uncomment to get very verbose mesh information
+
+# Change the default mesh reader to the CUBIT reader.
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+# Set filename of mesh to import.
+filename = ../mesh/threehex8.exo
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Nondimensionalization
+normalizer.length_scale = 1.0*m
+normalizer.shear_modulus = 4.0e5*Pa
+normalizer.relaxation_time = 1.0*s
+
+[pylithapp.timedependent.formulation.time_step]
+# Define the total time for the simulation and the default time step size.
+total_time = 0.0*s ; total time of simulation
+
+# ----------------------------------------------------------------------
+# PETSc and solution settings.
+# ----------------------------------------------------------------------
+# Set the solver options.
+
+[pylithapp.petsc]
+
+# Preconditioner settings.
+# ksp_type = gmres
+pc_type = asm
+# comment this out for fieldSplit
+sub_pc_factor_shift_type = nonzero
+
+# Convergence parameters.
+ksp_rtol = 1.0e-10
+ksp_atol = 1.0e-22
+ksp_gmres_restart = 150
+ksp_max_it = 1500
+
+# Linear solver monitoring options.
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+ksp_monitor_singular_value = true
+log_summary = true
+
+# Field split
+# [pylithapp.timedependent.formulation]
+# split_fields = True
+# use_custom_constraint_pc = True
+# matrix_type = aij
+
+[pylithapp.petsc]
+# fs_pc_type = fieldsplit
+# fs_pc_fieldsplit_real_diagonal =
+# fs_pc_fieldsplit_type = multiplicative
+# fs_fieldsplit_0_pc_type = ml
+# fs_fieldsplit_1_pc_type = ml
+# fs_fieldsplit_2_pc_type = ml
+# fs_fieldsplit_3_pc_type = ml
+# fs_fieldsplit_0_ksp_type = preonly
+# fs_fieldsplit_1_ksp_type = preonly
+# fs_fieldsplit_2_ksp_type = preonly
+# fs_fieldsplit_3_ksp_type = preonly
+
+# Nonlinear solver monitoring options.
+snes_rtol = 1.0e-12
+snes_atol = 1.0e-22
+snes_max_it = 25
+snes_monitor = true
+snes_view = true
+snes_converged_reason = true
+
+# Uncomment to launch gdb when starting PyLith.
+# start_in_debugger = true
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/visc_e10_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/visc_e10_ng_ni_v03.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/config/visc_e10_ng_ni_v03.cfg 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,128 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This file is for a simulation with the following characteristics:
+# 1. The material is Maxwell viscoelastic with a viscosity of
+# 1.3e10 Pa-s.
+# 2. There is no gravity.
+# 3. There are no initial stresses.
+# 4. Velocity BC are applied at a rate of 0.03 mm/yr for 10 years.
+# 5. The linear solver is used.
+#
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to this problem.
+# The general parameters are specificed in the pylithapp.cfg
+# file which PyLith reads by default.
+#
+# To run the simulation:
+# pylith visc_e10_ng_ni_v03.cfg
+#
+# Output will be directed to the directory ../results.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
+# 'z_pos', and 'z_neg'.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+
+[pylithapp.timedependent.formulation]
+# Set the output to an array of 1 output manager.
+# We will output the solution over the domain.
+output = [domain]
+
+# Use adaptive time stepping.
+[pylithapp.timedependent.formulation]
+time_step = pylith.problems.TimeStepAdapt
+
+# Change the total simulation time to 10 years, and set the maximum time
+# step size to 0.1 years.
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+max_dt = 0.1*year
+stability_factor = 0.1
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [viscous_material]
+materials.viscous_material = pylith.materials.MaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.viscous_material]
+label = Viscous material
+id = 1
+db_properties.iohandler.filename = ../spatialdb/mat_viscous_e10.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+output.cell_info_fields = [density,mu,lambda,maxwell_time]
+output.cell_data_fields = [total_strain,stress,viscous_strain]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# +x face
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0]
+label = face_xpos
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +x
+db_rate.values = [displacement-rate-x,rate-start-time]
+db_rate.data = [-0.03*mm/year,0.0*year]
+
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = face_xneg
+db_initial.label = Dirichlet BC on -x
+
+# -z face
+[pylithapp.timedependent.bc.z_neg]
+bc_dof = [2]
+label = face_zneg
+db_initial.label = Dirichlet BC on -z
+
+# +z face
+[pylithapp.timedependent.bc.z_pos]
+bc_dof = [2]
+label = face_zpos
+db_initial.label = Dirichlet BC on +z
+
+# -y face
+[pylithapp.timedependent.bc.y_neg]
+bc_dof = [0, 1]
+label = face_yneg
+db_rate = spatialdata.spatialdb.SimpleDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.iohandler.filename = ../spatialdb/basal_rate_bc_03.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain]
+skip = 0
+writer.filename = ../results/visc_e10_ng_ni_v03.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.viscous_material.output]
+# Average values over quadrature points.
+cell_filter = pylith.meshio.CellFilterAvgMesh
+skip = 0
+writer.filename = ../results/visc_e10_ng_ni_v03-viscous_material.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,24 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA = \
+ geometry.jou \
+ mesh_threehex8.jou \
+ threehex8.exo
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/geometry.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/geometry.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/geometry.jou 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,26 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 3 m x 1 m x 1 m
+# -1.5 m <= x <= 1.5 m
+# -0.5 m <= y <= 0.5 m
+# -0.5 m <= z <= 0.5 m
+reset
+brick x 3 y 1 z 1
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/mesh_threehex8.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/mesh_threehex8.jou (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/mesh_threehex8.jou 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,81 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 1
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Create block for material
+# ----------------------------------------------------------------------
+block 1 volume 1
+block 1 name "plastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 6
+# We remove the nodes in face_yneg so there is no conflict in BC.
+group "face_xpos" remove node in surface 3
+nodeset 11 group face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 4
+# We remove the nodes in face_yneg so there is no conflict in BC.
+group "face_xneg" remove node in surface 3
+nodeset 12 group face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 5
+nodeset 13 group face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 3
+nodeset 14 group face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 2
+nodeset 15 group face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 1
+nodeset 16 group face_zpos
+nodeset 16 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "threehex8.exo" dimension 3 overwrite
+
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/threehex8.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/mesh/threehex8.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/results/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/results/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/results/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,22 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA =
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/Makefile.am
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/Makefile.am (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/Makefile.am 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,43 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA = \
+ basal_bc_3.spatialdb \
+ basal_rate_bc_0015.spatialdb \
+ basal_rate_bc_012.spatialdb \
+ basal_rate_bc_03.spatialdb \
+ initial_stress.spatialdb \
+ mat_elastic.spatialdb \
+ mat_plastic_dp1.spatialdb \
+ mat_plastic_dp2.spatialdb \
+ mat_plastic_dp3.spatialdb \
+ mat_plastic_vm1.spatialdb \
+ mat_viscous_e01.spatialdb \
+ mat_viscous_e02.spatialdb \
+ mat_viscous_e03.spatialdb \
+ mat_viscous_e04.spatialdb \
+ mat_viscous_e05.spatialdb \
+ mat_viscous_e06.spatialdb \
+ mat_viscous_e07.spatialdb \
+ mat_viscous_e08.spatialdb \
+ mat_viscous_e09.spatialdb \
+ mat_viscous_e10.spatialdb \
+ mat_viscous_e11.spatialdb \
+ mat_viscous_e12.spatialdb
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_bc_3.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_bc_3.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_bc_3.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = displacement-x displacement-y displacement-z
+ value-units = mm mm mm
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Ux (mm)
+// (5) Uy (mm)
+// (6) Uz (mm)
+-1.5000 0.0 0.0 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.3 0.0 0.0
+ 1.5000 0.0 0.0 -0.3 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_0015.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_0015.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_0015.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = displacement-rate-x displacement-rate-y displacement-rate-z rate-start-time
+ value-units = mm/year mm/year mm/year year
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Vx (mm/year)
+// (5) Vy (mm/year)
+// (6) Vz (mm/year)
+// (7) t0 (year)
+-1.5000 0.0 0.0 0.0000 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.0000 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.0015 0.0 0.0 0.0
+ 1.5000 0.0 0.0 -0.0015 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_012.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_012.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_012.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = displacement-rate-x displacement-rate-y displacement-rate-z rate-start-time
+ value-units = mm/year mm/year mm/year year
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Vx (mm/year)
+// (5) Vy (mm/year)
+// (6) Vz (mm/year)
+// (7) t0 (year)
+-1.5000 0.0 0.0 0.000 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.000 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.012 0.0 0.0 0.0
+ 1.5000 0.0 0.0 -0.012 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_03.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_03.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/basal_rate_bc_03.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of the
+// displacement field for Dirichlet boundary conditions.
+//
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = displacement-rate-x displacement-rate-y displacement-rate-z rate-start-time
+ value-units = mm/year mm/year mm/year year
+ num-locs = 4
+ data-dim = 1 // locations form a line
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) Vx (mm/year)
+// (5) Vy (mm/year)
+// (6) Vz (mm/year)
+// (7) t0 (year)
+-1.5000 0.0 0.0 0.000 0.0 0.0 0.0
+-0.0100 0.0 0.0 0.000 0.0 0.0 0.0
+ 0.0100 0.0 0.0 -0.030 0.0 0.0 0.0
+ 1.5000 0.0 0.0 -0.030 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/initial_stress.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/initial_stress.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/initial_stress.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,36 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the initial stresses for the model
+// The stresses increase linearly with depth.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of stress components
+
+ // names of the stress components
+ value-names = stress-xx stress-yy stress-zz stress-xy stress-yz stress-xz
+
+ value-units = Pa Pa Pa Pa Pa Pa // units
+
+ num-locs = 2 // number of locations: 1 value each at the upper and lower
+ // surfaces. Linear interpolation should be used.
+
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) stress-xx (Pa)
+// (5) stress-yy (Pa)
+// (6) stress-zz (Pa)
+// (7) stress-xy (Pa)
+// (8) stress-yz (Pa)
+// (9) stress-xz (Pa)
+0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+0.0 -1.0 0.0 -100.0 -100.0 -100.0 0.0 0.0 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_elastic.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_elastic.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,27 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for an elastic material. In this case, the material
+// properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3 // number of material property values
+ value-names = density vs vp // names of the material property values
+ value-units = kg/m**3 m/s m/s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp1.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp1.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp1.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp2.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp2.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp2.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 5.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp3.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp3.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_dp3.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 10.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_vm1.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_vm1.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_plastic_vm1.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of material property values
+ value-names = density vs vp friction-angle cohesion dilatation-angle // names of the material property values
+ value-units = kg/m**3 m/s m/s degree Pa degree // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) friction-angle (degree)
+// (8) cohesion (Pa)
+// (8) dilatation-angle (degree)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 0.0 5.0e1 0.0
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e01.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e01.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e01.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e1
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e02.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e02.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e02.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e2
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e03.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e03.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e03.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e3
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e04.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e04.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e04.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e4
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e05.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e05.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e05.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e5
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e06.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e06.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e06.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e6
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e07.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e07.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e07.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e7
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e08.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e08.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e08.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e8
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e09.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e09.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e09.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e9
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e10.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e10.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e10.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e10
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e11.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e11.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e11.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e11
Added: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e12.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e12.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex8/spatialdb/mat_viscous_e12.spatialdb 2010-08-27 04:35:47 UTC (rev 17129)
@@ -0,0 +1,28 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a power-law viscoelastic material. In this case,
+// the material properties are uniform.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m**3 m/s m/s Pa*s // units
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa*s)
+0.0 0.0 0.0 10.0 200.0 346.4101615137754 1.3e12
More information about the CIG-COMMITS
mailing list