[cig-commits] r16690 - short/3D/PyLith/trunk/examples/3d/hex8
willic3 at geodynamics.org
willic3 at geodynamics.org
Tue May 11 18:58:06 PDT 2010
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
More information about the CIG-COMMITS
mailing list