[cig-commits] r16690 - short/3D/PyLith/trunk/examples/3d/hex8

[willic3 at geodynamics.org]

Author: willic3
Date: 2010-05-11 18:58:06 -0700 (Tue, 11 May 2010)
New Revision: 16690

Modified:
   short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step02.cfg
Log:
Made some minor changes to where things are defined, and added more
description.



Modified: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg	2010-05-12 00:59:42 UTC (rev 16689)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg	2010-05-12 01:58:06 UTC (rev 16690)
@@ -43,24 +43,20 @@
 [pylithapp.timedependent]
 
 # Set materials to an array of 2 materials:
-#   'elastic' and 'inelastic'
-materials = [elastic,inelastic]
+#   'upper_crust' and 'lower_crust'
+materials = [upper_crust,lower_crust]
 
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.inelastic = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
+[pylithapp.timedependent.materials.upper_crust]
+label = Upper crust material
 id = 1
 db_properties.iohandler.filename = spatialdb/mat_elastic.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
 
-[pylithapp.timedependent.materials.inelastic]
-label = Inelastic material
+[pylithapp.timedependent.materials.lower_crust]
+label = Lower crust material
 id = 2
-db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
+db_properties.iohandler.filename = spatialdb/mat_elastic.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
 

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg	2010-05-12 00:59:42 UTC (rev 16689)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg	2010-05-12 01:58:06 UTC (rev 16690)
@@ -6,7 +6,11 @@
 # ----------------------------------------------------------------------
 
 #
-# [ADD PHYSICAL DESCRIPTION OF THE PROBLEM BEING SOLVED]
+# This is a purely elastic static problem using only Dirichlet (displacement)
+# boundary conditions. The lower (minimum z) boundary is held fixed in the
+# z-direction. On the positive and negative x-boundaries, displacements are
+# applied corresponding to 1 m of shear displacement in the positive and
+# negative y-directions, and 1 m of axial (x) compression on each side.
 
 # ----------------------------------------------------------------------
 # RUNNING THE SIMULATON
@@ -91,18 +95,14 @@
 label = 17
 writer.filename = step01/step01-groundsurf.vtk
 
-# Give basename for VTK output of elastic state variables.
-[pylithapp.timedependent.materials.elastic.output]
+# Give basename for VTK output of upper_crust state variables.
+[pylithapp.timedependent.materials.upper_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = step01/step01-elastic.vtk
+writer.filename = step01/step01-upper_crust.vtk
 
-# Give basename for VTK output of inelastic state variables.
-[pylithapp.timedependent.materials.inelastic.output]
-# Default (elastic) info and data fields do not include viscoelastic
-# information, so we include it here.
-cell_info_fields = [density,mu,lambda,maxwell_time]
-cell_data_fields = [total_strain,stress,viscous_strain]
+# Give basename for VTK output of lower_crust state variables.
+[pylithapp.timedependent.materials.lower_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = step01/step01-inelastic.vtk
+writer.filename = step01/step01-lower_crust.vtk

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step02.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step02.cfg	2010-05-12 00:59:42 UTC (rev 16689)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step02.cfg	2010-05-12 01:58:06 UTC (rev 16690)
@@ -1,6 +1,23 @@
 # -*- Python -*-
 [pylithapp]
 
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# This is a purely elastic static problem using both Dirichlet (displacement)
+# boundary conditions and Neumann (traction) boundary conditions. The lower
+# (minimum z) boundary is held fixed in the z-direction. On the negative
+# x-boundary, displacements are held fixed in the x and y-directions.
+# On the positive x-boundary, tractions corresponding to 1 MPa of shear in
+# the +y direction and 1 MPa of normal stress in the -x direction are
+# applied.
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
 # This is not a self-contained simulation configuration file. This
 # file only specifies parameters specific to tutorial step02.
 # The general parameters are specificed in the pylithapp.cfg
@@ -27,49 +44,51 @@
 output.subdomain = pylith.meshio.OutputSolnSubset
 
 # ----------------------------------------------------------------------
-# materials
+# boundary conditions
 # ----------------------------------------------------------------------
+# Set the parameters for Dirichlet boundary conditions applied on the
+# -x and -z faces of the box, and the Neumann boundary conditions
+# applied on the +x face of the box.
+#
+# We fix the x and y degrees of freedom on the -x face, and
+# fix the z degree of freedom on the bottom (-z) face.
+#
+# We apply axial and shear tractions on the +x face.
+#
+# For x_pos, we change the database from ZeroDispDB (which
+# specifies a zero value) to SimpleDB (which uses a data file to
+# specify a spatial variation), and we change the boundary condition type
+# to Neumann (traction).
+#
+# For x_neg and z_neg, we keep the default ZeroDispDB, since we just want zero
+# displacements on those boundaries.
+#
 
-# We specifically define the database info for the inelastic material
-# for each model, since the type of spatialdb can be different.
+# The label corresponds to the nodeset ID in CUBIT.
 
-[pylithapp.timedependent.materials.inelastic]
-db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-# For x_neg, we use the default Dirichlet boundary condition type along
-# with the default ZeroDispDB spatial database type. This will provide
-# zero displacements in the x and y directions for that boundary.
-# For x_pos, we change the boundary condition type to Neumann (traction),
-# and we specify shear and normal tractions using a SimpleDB.
-# For z_neg, we keep the default ZeroDispDB, since we just want zero
-# z-displacements on that boundary.
-
-# Change the BC type to Neumann for x_pos.
+# +x face -- first change bc type to Neumann
 [pylithapp.timedependent.bc]
 x_pos = pylith.bc.Neumann
 
-# Specify zero displacement (Dirichlet) BC for x_neg.
-# The label corresponds to the nodeset ID in CUBIT.
-[pylithapp.timedependent.bc.x_neg]
-bc_dof = [0, 1]
-label = 12
-db_initial.label = Dirichlet BC on -x
-
-# Specify a spatial database (SimpleDB) to describe the Neumann BC on x_pos.
-# The label corresponds to the nodeset ID in CUBIT.
+# +x face
 [pylithapp.timedependent.bc.x_pos]
 label = 11
 db_initial = spatialdata.spatialdb.SimpleDB
 db_initial.label = Neumann BC on +x
 db_initial.iohandler.filename = spatialdb/tractions_axial_shear.spatialdb
+
 # We must specify quadrature information for the cell faces.
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 quadrature.cell.quad_order = 2
 
+# -x face
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0, 1]
+label = 12
+db_initial.label = Dirichlet BC on -x
+
+# -z face
 [pylithapp.timedependent.bc.z_neg]
 bc_dof = [2]
 label = 15
@@ -87,16 +106,12 @@
 label = 17 ; nodeset for subdomain
 writer.filename = step02/step02-groundsurf.vtk
 
-# Give basename for VTK output of elastic state variables.
-[pylithapp.timedependent.materials.elastic.output]
+# Give basename for VTK output of upper_crust state variables.
+[pylithapp.timedependent.materials.upper_crust.output]
 cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = step02/step02-elastic.vtk
+writer.filename = step02/step02-upper_crust.vtk
 
-# Give basename for VTK output of inelastic state variables.
-[pylithapp.timedependent.materials.inelastic.output]
-# Default (elastic) info and data fields do not include viscoelastic
-# information, so we include it here.
-cell_info_fields = [density,mu,lambda,maxwell_time]
-cell_data_fields = [total_strain,stress,viscous_strain]
+# Give basename for VTK output of lower_crust state variables.
+[pylithapp.timedependent.materials.lower_crust.output]
 cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = step02/step02-inelastic.vtk
+writer.filename = step02/step02-lower_crust.vtk