[cig-commits] r16030 - short/3D/PyLith/branches/pylith-friction/playpen/friction
brad at geodynamics.org
brad at geodynamics.org
Tue Nov 24 10:17:40 PST 2009
Author: brad
Date: 2009-11-24 10:17:39 -0800 (Tue, 24 Nov 2009)
New Revision: 16030
Added:
short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic2D.spatialdb
short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic3D.spatialdb
Removed:
short/3D/PyLith/branches/pylith-friction/playpen/friction/matprops.spatialdb
short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-matprops.spatialdb
short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-neumannBC.cfg
Modified:
short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-axial.cfg
short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-opening.cfg
short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8.mesh
short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.cfg
short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-dirichletBC.cfg
short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg
Log:
Minor cleanup of files.
Copied: short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic2D.spatialdb (from rev 16029, short/3D/PyLith/branches/pylith-friction/playpen/friction/matprops.spatialdb)
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic2D.spatialdb (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic2D.spatialdb 2009-11-24 18:17:39 UTC (rev 16030)
@@ -0,0 +1,46 @@
+#SPATIAL.ascii 1
+
+// This database is used to specify the material properties for all of the
+// examples in this directory.
+
+// This follows the format for a Simple DB (the only type presently available).
+SimpleDB {
+
+ // There are 3 values specified in the database, corresponding to density,
+ // S-velocity, and P-velocity (values for shear modulus and Lame's constant
+ // are computed from these values.
+ num-values = 3
+ value-names = density vs vp
+
+ // These are the units used to specify density, vs, and vp.
+ value-units = kg/m**3 m/s m/s
+
+ // Values are only specified at a single point since they are constant
+ // throughout the mesh.
+ num-locs = 1
+
+ // The dimension of the spatial distribution is 0, since it is constant
+ // throughout the mesh.
+ data-dim = 0
+
+ // The spatial dimension of the database is 2.
+ space-dim = 2
+
+ // We are specifying the data in a Cartesian coordinate system.
+ cs-data = cartesian {
+
+ // Our units are already in meters, so we can just multiply by one.
+ to-meters = 1.0
+
+ // We are using a 2D Cartesian coordinate system.
+ space-dim = 2
+ }
+}
+// This is where the data is specified.
+// We only need to specify a single point, since the properties are uniform.
+// The values given here will give a shear modulus and Lame's constant both
+// equal to 30 GPa (Poisson's ratio = 0.25).
+// The entries are:
+// X-coord, Y-coord, density, Vs, Vp.
+
+0.0 0.0 2700.0 3333.333333333333 5773.502691896258
Copied: short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic3D.spatialdb (from rev 16029, short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-matprops.spatialdb)
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic3D.spatialdb (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/matelastic3D.spatialdb 2009-11-24 18:17:39 UTC (rev 16030)
@@ -0,0 +1,17 @@
+#SPATIAL.ascii 1
+
+SimpleDB {
+
+ num-values = 3
+ value-names = density vs vp
+ value-units = kg/m**3 m/s m/s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258
Deleted: short/3D/PyLith/branches/pylith-friction/playpen/friction/matprops.spatialdb
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/matprops.spatialdb 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/matprops.spatialdb 2009-11-24 18:17:39 UTC (rev 16030)
@@ -1,46 +0,0 @@
-#SPATIAL.ascii 1
-
-// This database is used to specify the material properties for all of the
-// examples in this directory.
-
-// This follows the format for a Simple DB (the only type presently available).
-SimpleDB {
-
- // There are 3 values specified in the database, corresponding to density,
- // S-velocity, and P-velocity (values for shear modulus and Lame's constant
- // are computed from these values.
- num-values = 3
- value-names = density vs vp
-
- // These are the units used to specify density, vs, and vp.
- value-units = kg/m**3 m/s m/s
-
- // Values are only specified at a single point since they are constant
- // throughout the mesh.
- num-locs = 1
-
- // The dimension of the spatial distribution is 0, since it is constant
- // throughout the mesh.
- data-dim = 0
-
- // The spatial dimension of the database is 2.
- space-dim = 2
-
- // We are specifying the data in a Cartesian coordinate system.
- cs-data = cartesian {
-
- // Our units are already in meters, so we can just multiply by one.
- to-meters = 1.0
-
- // We are using a 2D Cartesian coordinate system.
- space-dim = 2
- }
-}
-// This is where the data is specified.
-// We only need to specify a single point, since the properties are uniform.
-// The values given here will give a shear modulus and Lame's constant both
-// equal to 30 GPa (Poisson's ratio = 0.25).
-// The entries are:
-// X-coord, Y-coord, density, Vs, Vp.
-
-0.0 0.0 2700.0 3333.333333333333 5773.502691896258
Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-axial.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-axial.cfg 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-axial.cfg 2009-11-24 18:17:39 UTC (rev 16030)
@@ -36,7 +36,7 @@
formulation.solver = pylith.problems.SolverNonlinear
# Set bc to an array with 3 boundary conditions: 'x_neg', 'z_neg' and 'x_pos'.
-bc = [x_neg,z_neg,x_pos]
+bc = [x_neg,z_neg,y_neg,x_pos]
bc.x_pos = pylith.bc.Neumann
# Set interfaces to an array with 1 fault: 'fault'.
@@ -61,7 +61,7 @@
# in the mesh file.
id = 1
-db_properties.iohandler.filename = twohex8-matprops.spatialdb
+db_properties.iohandler.filename = matelastic3D.spatialdb
# We are doing 3D quadrature for a hex8.
quadrature.cell = pylith.feassemble.FIATLagrange
@@ -78,8 +78,8 @@
# Boundary conditions to be applied to the negative x-side of the mesh.
[pylithapp.timedependent.bc.x_neg]
-# We are fixing the 0 (x) and 1 (y) degree of freedom.
-bc_dof = [0,1]
+# We are fixing the 0 (x) degree of freedom.
+bc_dof = [0]
# The nodes associated with this boundary condition have the name
# 'x_neg' in the mesh file.
@@ -95,6 +95,16 @@
# 'z_pos' in the mesh file.
label = z_neg_nofault
+# Boundary conditions to be applied to the negative z-side of the mesh.
+[pylithapp.timedependent.bc.y_neg]
+
+# We are fixing the 1 (y) degree of freedom.
+bc_dof = [1]
+
+# The nodes associated with this boundary condition have the name
+# 'z_pos' in the mesh file.
+label = y_neg_nofault
+
# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
Deleted: short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-matprops.spatialdb
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-matprops.spatialdb 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-matprops.spatialdb 2009-11-24 18:17:39 UTC (rev 16030)
@@ -1,17 +0,0 @@
-#SPATIAL.ascii 1
-
-SimpleDB {
-
- num-values = 3
- value-names = density vs vp
- value-units = kg/m**3 m/s m/s
- num-locs = 1
- data-dim = 0
- space-dim = 3
-
- cs-data = cartesian {
- to-meters = 1.0
- space-dim = 3
- }
-}
-0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258
Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-opening.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-opening.cfg 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8-opening.cfg 2009-11-24 18:17:39 UTC (rev 16030)
@@ -37,7 +37,6 @@
# Set bc to an array with 2 boundary conditions: 'x_neg' and 'x_pos'.
bc = [x_neg,x_pos]
-bc.x_pos = pylith.bc.DirichletBoundary
# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
@@ -61,7 +60,7 @@
# in the mesh file.
id = 1
-db_properties.iohandler.filename = twohex8-matprops.spatialdb
+db_properties.iohandler.filename = matelastic3D.spatialdb
# We are doing 3D quadrature for a hex8.
quadrature.cell = pylith.feassemble.FIATLagrange
Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8.mesh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8.mesh 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twohex8.mesh 2009-11-24 18:17:39 UTC (rev 16030)
@@ -80,8 +80,6 @@
}
}
- // This group of vertices may be used to specify boundary conditions.
- // There are 4 vertices corresponding to indices 0, 1, 6, and 7.
group = {
name = x_neg
type = vertices
@@ -91,8 +89,6 @@
}
}
- // This group of vertices may be used to specify boundary conditions.
- // There are 4 vertices corresponding to indices 4, 5, 10, and 11.
group = {
name = x_pos
type = vertices
@@ -102,8 +98,6 @@
}
}
- // This group of vertices may be used to specify boundary conditions.
- // There are 6 vertices corresponding to indices 0, 2, 4, 6, 8, and 10.
group = {
name = y_neg
type = vertices
@@ -113,9 +107,17 @@
}
}
- // This group of vertices may be used to specify boundary conditions.
- // There are 6 vertices corresponding to indices 0, 1, 2, 3, 4, and 5.
+ // y_neg with fault vertices removed
group = {
+ name = y_neg_nofault
+ type = vertices
+ count = 4
+ indices = {
+ 0 4 6 10
+ }
+ }
+
+ group = {
name = z_neg
type = vertices
count = 6
@@ -124,8 +126,7 @@
}
}
- // This group of vertices may be used to specify boundary conditions.
- // There are 4 vertices corresponding to indices 0, 1, 4, and 5.
+ // z_neg with fault vertices removed
group = {
name = z_neg_nofault
type = vertices
Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.cfg 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-opening.cfg 2009-11-24 18:17:39 UTC (rev 16030)
@@ -91,8 +91,8 @@
id = 0
# The properties for this material are given in the spatial database file
-# 'matprops.spatialdb'.
-db_properties.iohandler.filename = matprops.spatialdb
+# 'matelastic2D.spatialdb'.
+db_properties.iohandler.filename = matelastic2D.spatialdb
# Set cell type to quadrilateral (2-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-dirichletBC.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-dirichletBC.cfg 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-dirichletBC.cfg 2009-11-24 18:17:39 UTC (rev 16030)
@@ -91,8 +91,8 @@
id = 0
# The properties for this material are given in the spatial database file
-# 'matprops.spatialdb'.
-db_properties.iohandler.filename = matprops.spatialdb
+# 'matelastic2D.spatialdb'.
+db_properties.iohandler.filename = matelastic2D.spatialdb
# Set cell type to quadrilateral (2-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
Deleted: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-neumannBC.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-neumannBC.cfg 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear-neumannBC.cfg 2009-11-24 18:17:39 UTC (rev 16030)
@@ -1,218 +0,0 @@
-# -*- Python -*-
-
-# The settings in this file (pylithapp.cfg) will be read automatically
-# by pylith, as long as the file is placed in the run directory.
-
-# The settings in this file will override any settings in:
-# PREFIX/etc/pylithapp.cfg
-# $HOME/.pyre/pylithapp/pylithapp.cfg
-
-# The settings in this file will be overridden by any .cfg file given
-# on the command line or by any command line settings.
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-# The settings below turn on journal info for the specified components.
-# If you want less output to stdout, you can turn these off.
-[pylithapp.journal.info]
-timedependent = 1
-implicit = 1
-petsc = 1
-solvernonlinear = 1
-meshioascii = 1
-homogeneous = 1
-elasticityimplicit = 1
-fiatlagrange = 1
-quadrature1d = 1
-faultcohesivedyn = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-# The settings below control the mesh generation (importing mesh info).
-# Turn on debugging output for mesh generation.
-[pylithapp.mesh_generator]
-debug = 0
-
-# This component specification means we are using PyLith ASCII format,
-# and we then specify the filename and number of space dimensions for
-# the mesh.
-[pylithapp.mesh_generator.reader]
-filename = twoquad4-shear.mesh
-coordsys.space_dim = 2
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-# Specify the problem settings.
-# This is a time-dependent problem, so we select this as our problem type.
-# We select a total time of 0 sec, and a time step size of 1 sec, so we
-# are performing a single time step.
-# The spatial dimension for this problem is 2.
-# For an implicit formulation (using implicit.cfg), we will perform 1
-# implicit time step from t = -1.0 to t = 0.0 (elastic solution step).
-[pylithapp.timedependent]
-dimension = 2
-normalizer.length_scale = 1.0*m
-formulation = pylith.problems.Implicit
-formulation.solver = pylith.problems.SolverNonlinear
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos_disp', 'x_pos_tract'
-bc = [x_neg,x_pos_disp,x_pos_tract]
-bc.x_pos_tract = pylith.bc.Neumann
-
-# Set interfaces to an array with 1 fault: 'fault'.
-interfaces = [fault]
-
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 0.0*s
-dt = 1.0*s
-
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-# Specify the material information for the problem.
-# The material type is isotropic elastic formulated for plane strain.
-[pylithapp.timedependent.materials]
-material = pylith.materials.ElasticPlaneStrain
-
-[pylithapp.timedependent.materials.material]
-
-# We give a label of 'elastic material' to this material.
-label = elastic material
-
-# The cells associated with this material are given a material ID of 0
-# in the mesh file.
-id = 0
-
-# The properties for this material are given in the spatial database file
-# 'matprops.spatialdb'.
-db_properties.iohandler.filename = matprops.spatialdb
-
-# Set cell type to quadrilateral (2-d Lagrange).
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-# Provide information on the boundary conditions.
-
-# Boundary conditions to be applied to the negative x-side of the mesh.
-[pylithapp.timedependent.bc.x_neg]
-
-# We are fixing the 0 (x) and 1 (y) degrees of freedom.
-bc_dof = [0,1]
-
-# The nodes associated with this boundary condition have the name
-# 'x_neg' in the mesh file.
-label = x_neg
-
-# Boundary conditions to be applied to the positive x-side of the mesh.
-[pylithapp.timedependent.bc.x_pos_disp]
-
-# We are fixing the 0 (x) degree of freedom.
-bc_dof = [0,1]
-
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
-label = x_pos
-
-# Boundary conditions to be applied to the positive x-side of the mesh.
-[pylithapp.timedependent.bc.x_pos_tract]
-
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
-label = x_pos
-
-db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Neumann BC +x edge
-db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [12000*MPa,0.0*MPa,]
-
-# Set cell type to quadrilateral (1-d Lagrange).
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-
-fault = pylith.faults.FaultCohesiveDynL
-
-# Define fault properties.
-[pylithapp.timedependent.interfaces.fault]
-
-# The nodes associated with this fault have the name 'fault' in the mesh file.
-label = fault
-
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-# The quadrature for a 2D fault is 1D with a linear shape.
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-
-db_initial_tractions = spatialdata.spatialdb.UniformDB
-db_initial_tractions.label = "Initial fault tractions"
-db_initial_tractions.values = [traction-shear,traction-normal]
-db_initial_tractions.data = [0.0*Pa, -10000.0*MPa]
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-# We are using all of the default settings for PETSc except for
-# specifying the block Jacobi preconditioner. Additional PETSc
-# command-line arguments may be found in the PETSc documentation.
-[pylithapp.petsc]
-ksp_rtol = 1.0e-8
-pc_type = asm
-# Change the preconditioner settings (must turn off
-# shift_positive_definite and turn on shift_nonzero).
-sub_pc_factor_shift_positive_definite = 0
-sub_pc_factor_shift_nonzero =
-
-ksp_monitor = true
-ksp_view = true
-#log_summary = true
-ksp_max_it = 100
-ksp_gmres_restart = 50
-#start_in_debugger = true
-
-snes_monitor = true
-snes_view = true
-ksp_converged_reason = true
-snes_converged_reason = true
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output.writer]
-filename = twoquad4-shear-neumannBC.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = twoquad4-shear-neumannBC-fault.vtk
-vertex_info_fields = []
-vertex_data_fields = []
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.material.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = twoquad4-shear-neumannBC-statevars.vtk
Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg 2009-11-24 17:00:53 UTC (rev 16029)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg 2009-11-24 18:17:39 UTC (rev 16030)
@@ -91,8 +91,8 @@
id = 0
# The properties for this material are given in the spatial database file
-# 'matprops.spatialdb'.
-db_properties.iohandler.filename = matprops.spatialdb
+# 'matelastic2D.spatialdb'.
+db_properties.iohandler.filename = matelastic2D.spatialdb
# Set cell type to quadrilateral (2-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
More information about the CIG-COMMITS
mailing list