[cig-commits] r16837 - in short/3D/PyLith/trunk/examples/bar_shearwave: hex8 quad4 tet4 tri3
brad at geodynamics.org
brad at geodynamics.org
Sun May 30 13:58:09 PDT 2010
Author: brad
Date: 2010-05-30 13:58:09 -0700 (Sun, 30 May 2010)
New Revision: 16837
Added:
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_ratestateageing.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_slipweakening.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_staticfriction.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematic.cfg
Removed:
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg
Modified:
short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg
Log:
Removed static quad4 examples (superceded by 3d/hx8 examples). Adjusted settings for consistency and discussion in manual.
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,11 +1,28 @@
# -*- Python -*-
+[pylithapp]
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
-[pylithapp]
+# Impose 1.0 m of left-lateral slip at time 0.1 s with slip rate
+# following Brune's far-field time function. The rise time (95% of the
+# final slip) is 2.0 s.
+#
+# The fault slip causes a shear wave to propagate down the bar. The
+# shear wave is absorbed at the lateral ends of the bar using the
+# absorbing dampers boundary condition.
# ----------------------------------------------------------------------
+# RUNNING THE SIMULATION
+# ----------------------------------------------------------------------
+
+# This is a self-contained simulation configuration file.
+#
+# To run the simulation:
+# pylith
+
+# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
[pylithapp.journal.info]
@@ -199,19 +216,20 @@
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-skip = 3
+vertex_data_fields = [displacement,velocity]
+skip = 1
writer.filename = output/shearwave.vtk
writer.time_format = %05.2f
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
+skip = 1
writer.filename = output/shearwave-fault.vtk
writer.time_format = %05.2f
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-skip = 33
+skip = 19
writer.filename = output/shearwave-statevars.vtk
writer.time_format = %05.2f
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am 2010-05-30 20:58:09 UTC (rev 16837)
@@ -20,15 +20,13 @@
matprops.spatialdb \
mesh_quad4_200m.jou \
pylithapp.cfg \
- shearwave_kinematicslip.cfg \
+ kinematic.cfg \
kinematic_slip.spatialdb \
kinematic_risetime.spatialdb \
kinematic_sliptime.spatialdb \
- dynamicfault.cfg \
- static_stick.cfg \
- static_sliding.cfg \
- shearwave_staticfriction.cfg \
- shearwave_slipweakening.cfg \
- shearwave_ratestateageing.cfg
+ dynamic.cfg \
+ dynamic_staticfriction.cfg \
+ dynamic_slipweakening.cfg \
+ dynamic_ratestateageing.cfg
# End of file
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README 2010-05-30 20:58:09 UTC (rev 16837)
@@ -8,17 +8,13 @@
* Kinematic fault slip (same as problem in tri3, tet4, and hex8
subdirectories)
- * Static problem with fault friction without slip
+ * Dynamic fault slip with static friction fault constitutive model
- * Static problem with fault friction with slip
+ * Dynamic fault slip with linear slip-weakening fault constitutive model
- * Dynamic problem with static friction fault constitutive model
+ * Dynamic fault slip with rate- and state-friction fault constitutive model
- * Dynamic problem with linear slip-weakening fault constitutive model
- * Dynamic problem with rate- and state-friction fault constitutive model
-
-
Files common to all example problems
- geometry.jou: CUBIT journal file defining the geometry of the mesh
@@ -30,30 +26,22 @@
Files associated with the example problems:
+ Kinematic fault slip
- - shearwave_kinematicslip.cfg: Parameters for simulation
+ - kinematic.cfg: Parameters for simulation
- kinematic_risetime.spatialdb: Spatial database for slip rise time
- kinematic_slip.spatialdb: Spatial database for final slip
- kinematic_sliptime.spatialdb: Spatial database for slip initiation time
- + Static problem with fault friction without slip
+ + Dynamic fault slip with static friction fault constitutive model
- - static_stick.cfg: Parameters for simulation
+ - dynamic.cfg: Parameters common to dynamic friction problems
+ - dynamic_staticfriction.cfg: Parameters for simulation
- + Static problem with fault friction with slip
+ + Dynamic fault slip with linear slip-weakening fault constitutive model
- - static_sliding.cfg: Parameters for simulation
+ - dynamic.cfg: Parameters common to dynamic friction problems
+ - dynamic_slipweakening.cfg: Parameters for simulation
- + Dynamic problem with static friction fault constitutive model
+ + Dynamic fault slip with rate- and state-friction fault constitutive model
- - dynamicfault.cfg: Parameters common to dynamic friction problems
- - shearwave_staticfriction.cfg: Parameters for simulation
-
- + Dynamic problem with linear slip-weakening fault constitutive model
-
- - dynamicfault.cfg: Parameters common to dynamic friction problems
- - shearwave_slipweakening.cfg: Parameters for simulation
-
- + Dynamic problem with rate- and state-friction fault constitutive model
-
- - dynamicfault.cfg: Parameters common to dynamic friction problems
- - shearwave_ratestateageing.cfg: Parameters for simulation
+ - dynamic.cfg: Parameters common to dynamic friction problems
+ - dynamic_ratestateageing.cfg: Parameters for simulation
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic.cfg (from rev 16833, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -0,0 +1,49 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies the general parameters common to the dynamic
+# fault (friction interface) simulations in this directory.
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on progress for the dynamic fault.
+[pylithapp.journal.info]
+faultcohesivedyn = 1
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+# Overwrite the use of explicit time stepping with a full Jacobian
+# matrix set in pylithapp.cfg with a lumped Jacobian matrix. The
+# lumped Jacobian reduces memory use and uncouples the degrees of
+# freedom in the solve, permitting use of a very fast, simple solver.
+formulation = pylith.problems.ExplicitLumped
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+# Change fault to dynamic fault interface.
+fault = pylith.faults.FaultCohesiveDyn
+
+[pylithapp.timedependent.interfaces.fault]
+# Specify the initial tractions on the fault using a uniform DB.
+# shear: 6.1 MPa (right-lateral)
+# normal 10 MPa (compressive)
+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 = [-6.1*MPa, -10.0*MPa]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault.output]
+# Add slip rate to the default output fiels of slip and traction.
+vertex_data_fields=[slip,slip_rate,traction]
+skip = 1
+writer.time_format = %05.2f
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_ratestateageing.cfg (from rev 16833, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_ratestateageing.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_ratestateageing.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -0,0 +1,66 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+# Use constant initial tractions on the fault to drive fault slip
+# subject to rate- and state-friction on the fault. Deformation
+# propagates to the edges of the bar where it is absorbed using
+# absorbing dampers boundary conditions. The steady state solution is
+# a constant slip rate on the fault with uniform strain in the bar.
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATION
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for an explicit
+# time-stepping problem using a dynamic fault interface condition with
+# static friction.
+#
+# To run the simulation:
+# pylith dynamicfault.cfg staticfriction.cfg
+#
+# Output will be directed to the directory output.
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+# Use rate- and state-friction with the ageing law for evolution of
+# the state variable.
+friction = pylith.friction.RateStateAgeing
+
+# Set the friction model parameters.
+# reference coefficient of friction: 0.6
+# reference slip rate: 1.0e-06 m/s
+# slip weakening parameter: 0.037 m
+# a: 0.0125
+# b: 0.0172
+# cohesion: 0 Pa
+friction.db_properties = spatialdata.spatialdb.UniformDB
+friction.db_properties.label = Rate State Ageing
+friction.db_properties.values = [reference-friction-coefficient,reference-slip-rate,characteristic-slip-distance,constitutive-parameter-a,constitutive-parameter-b,cohesion]
+friction.db_properties.data = [0.6,1.0e-6*m/s,0.0370*m,0.0125,0.0172,0.0*Pa]
+
+# Set spatial database for the initial value of the state variable.
+friction.db_initial_state = spatialdata.spatialdb.UniformDB
+friction.db_initial_state.label = Rate State Ageing State
+friction.db_initial_state.values = [state-variable]
+friction.db_initial_state.data = [92.7*s]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Set filenames for output.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/shearwave-ratestateageing.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/shearwave-ratestateageing-fault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/shearwave-ratestateageing-statevars.vtk
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_slipweakening.cfg (from rev 16833, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_slipweakening.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_slipweakening.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -0,0 +1,65 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+# Use constant initial tractions on the fault to drive fault slip
+# subject to linear slip-weakening friction on the fault. Deformation
+# propagates to the edges of the bar where it is absorbed using
+# absorbing dampers boundary conditions. The steady state solution is
+# a constant slip rate on the fault with uniform strain in the bar.
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATION
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for an explicit
+# time-stepping problem using a dynamic fault interface condition with
+# linear slip-weakening friction.
+#
+# To run the simulation:
+# pylith dynamicfault.cfg slipweakening.cfg
+#
+# Output will be directed to the directory output.
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+# Use slip-weakening friction
+friction = pylith.friction.SlipWeakening
+
+# Set the slip-weakening friction model parameters.
+# static coefficient of friction: 0.6
+# dynamic coefficient of friction: 0.5
+# slip-weakening parameter: 0.2 m
+# cohesion: 0 Pa
+friction.db_properties = spatialdata.spatialdb.UniformDB
+friction.db_properties.label = Slip weakening
+friction.db_properties.values = [static-coefficient,dynamic-coefficient,slip-weakening-parameter,cohesion]
+friction.db_properties.data = [0.6,0.5,0.2*m,0.0*Pa]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/slipweakening.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/slipweakening-fault.vtk
+
+# Add friction model parameters to fault info file in addition to
+# default values.
+vertex_info_fields = [strike_dir,normal_dir,initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
+
+# Add output of state variables cumulative slip and previous slip to default
+# values.
+vertex_data_fields = [slip,traction,cumulative_slip,previous_slip]
+
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/slipweakening-statevars.vtk
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_staticfriction.cfg (from rev 16833, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_staticfriction.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamic_staticfriction.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -0,0 +1,53 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+# Use constant initial tractions on the fault to drive fault slip
+# subject to static friction on the fault. Deformation propagates to
+# the edges of the bar where it is absorbed using absorbing dampers
+# boundary conditions. The steady state solution is a constant slip
+# rate on the fault with uniform strain in the bar.
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATION
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for an explicit
+# time-stepping problem using a dynamic fault interface condition with
+# static friction.
+#
+# To run the simulation:
+# pylith dynamicfault.cfg staticfriction.cfg
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+# Use static friction
+friction = pylith.friction.StaticFriction
+
+# Set the static friction model parameters.
+# static coefficient of friction: 0.6
+# cohesion: 0 Pa
+friction.db_properties = spatialdata.spatialdb.UniformDB
+friction.db_properties.label = Static friction
+friction.db_properties.values = [friction-coefficient,cohesion]
+friction.db_properties.data = [0.6,0.0*Pa]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Set filenames for output.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/staticfriction.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/staticfriction-fault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/staticfriction-statevars.vtk
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,78 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the general parameters common to the dynamic
-# fault (friction interface) simulations in this directory.
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-# Turn on progress for the dynamic fault.
-[pylithapp.journal.info]
-faultcohesivedyn = 1
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Overwrite the use of explicit time stepping with a full Jacobian
-# matrix set in pylithapp.cfg with a lumped Jacobian matrix. The
-# lumped Jacobian reduces memory use and uncouples the degrees of
-# freedom in the solve, permitting use of a very fast, simple solver.
-formulation = pylith.problems.ExplicitLumped
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc]
-# Change the default Dirichlet BC to the AbsorbingDampers BC
-x_pos = pylith.bc.AbsorbingDampers
-x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing BC on +x face of bar
-
-# Name of nodeset in CUBIT
-label = face_xpos
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing BC on -x face of bar
-
-# Name of nodeset in CUBIT
-label = face_xneg
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces]
-# Change fault to dynamic fault interface.
-fault = pylith.faults.FaultCohesiveDyn
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault.output]
-# Add slip rate to the default output fiels of slip and traction.
-vertex_data_fields=[slip,slip_rate,traction]
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematic.cfg (from rev 16833, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematic.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematic.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -0,0 +1,83 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+# Impose 1.0 m of left-lateral slip at time 0.1 s with slip rate
+# following Brune's far-field time function. The rise time (95% of the
+# final slip) is 2.0 s.
+#
+# The fault slip causes a shear wave to propagate down the bar. The
+# shear wave is absorbed at the lateral ends of the bar using the
+# absorbing dampers boundary condition.
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATION
+# ----------------------------------------------------------------------
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for an explicit
+# time-stepping problem using a fault with kinematically-specified
+# slip.
+#
+# To run the simulation:
+# pylith shearwave_kinematicslip.cfg
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on progress for the kinematic fault.
+[pylithapp.journal.info]
+faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+# Uncomment the line below to overwrite the use of explicit time
+# stepping with a full Jacobian matrix set in pylithapp.cfg with a
+# lumped Jacobian matrix. The lumped Jacobian reduces memory use and
+# uncouples the degrees of freedom in the solve, permitting use of a
+# very fast, simple solver.
+#formulation = pylith.problems.ExplicitLumped
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault]
+# Switch to Brune slip time function
+eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+
+# Database specifying the final slip.
+slip.iohandler.filename = kinematic_slip.spatialdb
+
+# Database specifying rise_time.
+rise_time.iohandler.filename = kinematic_risetime.spatialdb
+
+# Database specifying time at which slip begins at each point.
+slip_time.iohandler.filename = kinematic_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Set the filename for each type of output.
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/kinematic.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/kinematic-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/kinematic-statevars.vtk
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -95,6 +95,43 @@
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc]
+# Change the default Dirichlet BC to the AbsorbingDampers BC
+x_pos = pylith.bc.AbsorbingDampers
+x_neg = pylith.bc.AbsorbingDampers
+
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing BC on +x face of bar
+
+# Name of nodeset in CUBIT
+label = face_xpos
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 1-D Lagrange cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+# Absorbing BC on -x face of bar
+
+# Name of nodeset in CUBIT
+label = face_xneg
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 1-D Lagrange cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
[pylithapp.timedependent.bc.all_nofault]
# Dirichlet boundary condition on all vertices except fault vertices
@@ -107,6 +144,7 @@
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.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.
@@ -118,8 +156,6 @@
# 2. If you use the default fault ID, make sure that none of your
# material ID's are equal to 100.
-[pylithapp.timedependent.interfaces.fault]
-
# Material id
id = 100
@@ -169,17 +205,12 @@
# writer.time_format - Format of time stamp (decimal point is removed)
# We will set the output filename in the problem specific .cfg file.
[pylithapp.problem.formulation.output.output]
-skip = 3
+vertex_data_fields = [displacement,velocity]
+skip = 1
writer.time_format = %05.2f
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
-writer.time_format = %05.2f
-
-# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-skip = 33
+skip = 19
writer.time_format = %05.2f
# Averate state variables over quadrature points.
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,121 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# PROBLEM DESCRIPTION
-# ----------------------------------------------------------------------
-
-# Impose 1.0 m of left-lateral slip at time 0.1 s with slip rate
-# following Brune's far-field time function. The rise time (95% of the
-# final slip) is 2.0 s.
-#
-# The fault slip causes a shear wave to propagate down the bar. The
-# shear wave is absorbed at the lateral ends of the bar using the
-# absorbing dampers boundary condition.
-
-# ----------------------------------------------------------------------
-# RUNNING THE SIMULATION
-# ----------------------------------------------------------------------
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for an explicit
-# time-stepping problem using a fault with kinematically-specified
-# slip.
-#
-# To run the simulation:
-# pylith shearwave_kinematicslip.cfg
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-# Turn on progress for the kinematic fault.
-[pylithapp.journal.info]
-faultcohesivekin = 1
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Uncomment the line below to overwrite the use of explicit time
-# stepping with a full Jacobian matrix set in pylithapp.cfg with a
-# lumped Jacobian matrix. The lumped Jacobian reduces memory use and
-# uncouples the degrees of freedom in the solve, permitting use of a
-# very fast, simple solver.
-#formulation = pylith.problems.ExplicitLumped
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc]
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-x_pos = pylith.bc.AbsorbingDampers
-x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing boundary condition on +x face of bar
-
-# Name of nodeset from CUBIT
-label = face_xpos
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing boundary condition on -x face of bar
-
-# Name of nodeset from CUBIT
-label = face_xneg
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-
-# Switch to Brune slip time function
-eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
-
-[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
-
-# Database specifying the final slip.
-slip.iohandler.filename = kinematic_slip.spatialdb
-
-# Database specifying rise_time.
-rise_time.iohandler.filename = kinematic_risetime.spatialdb
-
-# Database specifying time at which slip begins at each point.
-slip_time.iohandler.filename = kinematic_sliptime.spatialdb
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Set the filename for each type of output.
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/kinematicslip.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/kinematicslip-fault.vtk
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/kinematicslip-statevars.vtk
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,74 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# PROBLEM DESCRIPTION
-# ----------------------------------------------------------------------
-
-# Use constant initial tractions on the fault to drive fault slip
-# subject to rate- and state-friction on the fault. Deformation
-# propagates to the edges of the bar where it is absorbed using
-# absorbing dampers boundary conditions. The steady state solution is
-# a constant slip rate on the fault with uniform strain in the bar.
-
-# ----------------------------------------------------------------------
-# RUNNING THE SIMULATION
-# ----------------------------------------------------------------------
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for an explicit
-# time-stepping problem using a dynamic fault interface condition with
-# static friction.
-#
-# To run the simulation:
-# pylith dynamicfault.cfg shearwave_staticfriction.cfg
-#
-# Output will be directed to the directory output.
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-
-# Specify the initial tractions on the fault using a uniform DB.
-# shear: 6.1 MPa (right-lateral)
-# normal 10 MPa (compressive)
-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 = [-6.1*MPa, -10.0*MPa]
-
-# Use rate- and state-friction with the ageing law for evolution of
-# the state variable.
-friction = pylith.friction.RateStateAgeing
-
-# Set the friction model parameters.
-# reference coefficient of friction: 0.6
-# reference slip rate: 1.0e-06 m/s
-# slip weakening parameter: 0.037 m
-# a: 0.0125
-# b: 0.0172
-# cohesion: 0 Pa
-friction.db_properties = spatialdata.spatialdb.UniformDB
-friction.db_properties.label = Rate State Ageing
-friction.db_properties.values = [reference-friction-coefficient,reference-slip-rate,characteristic-slip-distance,constitutive-parameter-a,constitutive-parameter-b,cohesion]
-friction.db_properties.data = [0.6,1.0e-6*m/s,0.0370*m,0.0125,0.0172,0.0*Pa]
-
-# Set spatial database for the initial value of the state variable.
-friction.db_initial_state = spatialdata.spatialdb.UniformDB
-friction.db_initial_state.label = Rate State Ageing State
-friction.db_initial_state.values = [state-variable]
-friction.db_initial_state.data = [92.7*s]
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Set filenames for output.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/shearwave-ratestateageing.vtk
-
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/shearwave-ratestateageing-fault.vtk
-
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/shearwave-ratestateageing-statevars.vtk
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,73 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# PROBLEM DESCRIPTION
-# ----------------------------------------------------------------------
-
-# Use constant initial tractions on the fault to drive fault slip
-# subject to linear slip-weakening friction on the fault. Deformation
-# propagates to the edges of the bar where it is absorbed using
-# absorbing dampers boundary conditions. The steady state solution is
-# a constant slip rate on the fault with uniform strain in the bar.
-
-# ----------------------------------------------------------------------
-# RUNNING THE SIMULATION
-# ----------------------------------------------------------------------
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for an explicit
-# time-stepping problem using a dynamic fault interface condition with
-# linear slip-weakening friction.
-#
-# To run the simulation:
-# pylith dynamicfault.cfg shearwave_slipweakening.cfg
-#
-# Output will be directed to the directory output.
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-
-# Specify the initial tractions on the fault using a uniform DB.
-# shear: 6.1 MPa (right-lateral)
-# normal 10 MPa (compressive)
-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 = [-6.1*MPa, -10.0*MPa]
-
-# Use slip-weakening friction
-friction = pylith.friction.SlipWeakening
-
-# Set the slip-weakening friction model parameters.
-# static coefficient of friction: 0.6
-# dynamic coefficient of friction: 0.5
-# slip-weakening parameter: 0.2 m
-# cohesion: 0 Pa
-friction.db_properties = spatialdata.spatialdb.UniformDB
-friction.db_properties.label = Slip weakening
-friction.db_properties.values = [static-coefficient,dynamic-coefficient,slip-weakening-parameter,cohesion]
-friction.db_properties.data = [0.6,0.5,0.2*m,0.0*Pa]
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/shearwave-slipweakening.vtk
-
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/shearwave-slipweakening-fault.vtk
-
-# Add friction model parameters to fault info file in addition to
-# default values.
-vertex_info_fields = [strike_dir,normal_dir,initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
-
-# Add output of state variables cumulative slip and previous slip to default
-# values.
-vertex_data_fields = [slip,traction,cumulative_slip,previous_slip]
-
-
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/shearwave-slipweakening-statevars.vtk
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,61 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# PROBLEM DESCRIPTION
-# ----------------------------------------------------------------------
-
-# Use constant initial tractions on the fault to drive fault slip
-# subject to static friction on the fault. Deformation propagates to
-# the edges of the bar where it is absorbed using absorbing dampers
-# boundary conditions. The steady state solution is a constant slip
-# rate on the fault with uniform strain in the bar.
-
-# ----------------------------------------------------------------------
-# RUNNING THE SIMULATION
-# ----------------------------------------------------------------------
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for an explicit
-# time-stepping problem using a dynamic fault interface condition with
-# static friction.
-#
-# To run the simulation:
-# pylith dynamicfault.cfg shearwave_staticfriction.cfg
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-
-# Specify the initial tractions on the fault using a uniform DB.
-# shear: 6.1 MPa (right-lateral)
-# normal 10 MPa (compressive)
-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 = [-6.1*MPa, -10.0*MPa]
-
-# Use static friction
-friction = pylith.friction.StaticFriction
-
-# Set the static friction model parameters.
-# static coefficient of friction: 0.6
-# cohesion: 0 Pa
-friction.db_properties = spatialdata.spatialdb.UniformDB
-friction.db_properties.label = Static friction
-friction.db_properties.values = [friction-coefficient,cohesion]
-friction.db_properties.data = [0.6,0.0*Pa]
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Set filenames for output.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/dynamic-staticfriction.vtk
-
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/dynamic-staticfriction-fault.vtk
-
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/dynamic-staticfriction-statevars.vtk
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,147 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# PROBLEM DESCRIPTION
-# ----------------------------------------------------------------------
-
-# Impose Dirichlet BC on the lateral ends of the bar to create uniform
-# shear. Friction on the fault controlled by a static friction fault
-# constitutive model is insufficient to keep the falt locked, so slip
-# occurs.
-
-# This problem differs significantly from most real applications
-# involving fault slip, because the deformation is uniformly
-# distributed. As a result, when computing the fault slip associated
-# with the fault constitutive model imposed limits on fault shear
-# tractions, a few hundred iterations are needed in the nonlinear
-# solver. In most real applications only about 10 iterations are
-# needed for convergence.
-
-# ----------------------------------------------------------------------
-# RUNNING THE SIMULATION
-# ----------------------------------------------------------------------
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for a static problem
-# with fault slip controlled by a fault constitutive model.
-#
-# To run the simulation:
-# pylith static_sliding.cfg
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-# Turn on progress for the dynamic fault.
-[pylithapp.journal.info]
-faultcohesivedyn = 1
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Change to an implicit time-stepping formulation and the nonlinear
-# solver.
-formulation = pylith.problems.Implicit
-formulation.solver = pylith.problems.SolverNonlinear
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and
-# 'all_nofault'.
-bc = [x_pos,x_neg,all_nofault]
-
-[pylithapp.timedependent.formulation.time_step]
-# Switch to a static simulation.
-total_time = 0.0*s
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc]
-# Switch from absorbing boundary condition to Dirichlet BC.
-x_pos = pylith.bc.DirichletBC
-x_neg = pylith.bc.DirichletBC
-
-[pylithapp.timedependent.bc.x_pos]
-# Dirichlet BC on +x face of bar
-
-# Name of nodeset in CUBIT.
-label = face_xpos
-
-# Pin the y-dof.
-bc_dof = [1]
-
-# Use a uniform DB to impose the Dirichlet BC. Impose 320 m of
-# displacement in the y direction.
-db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Dirichlet BC +x edge
-db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*m,320.0*m]
-
-[pylithapp.timedependent.bc.x_neg]
-# Dirichlet BC on -x face of bar. Use the default (zero displacement) values.
-
-# Name of nodeset in CUBIT
-label = face_xneg
-
-# Pin the y-dof.
-bc_dof = [1]
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------x
-[pylithapp.timedependent.interfaces]
-# Change fault to dynamic fault interface.
-fault = pylith.faults.FaultCohesiveDyn
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Impose initial compressive normal tractions of -100 MPa. The initial
-# shear traction is zero.
-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*MPa, -100.0*MPa]
-
-# Use the static friction model.
-friction = pylith.friction.StaticFriction
-
-# Set static friction model parameters using a uniform DB. Set the
-# static coefficient of friction to 0.6 and cohesion to 0.0 Pa.
-friction.db_properties = spatialdata.spatialdb.UniformDB
-friction.db_properties.label = Static friction
-friction.db_properties.values = [friction-coefficient,cohesion]
-friction.db_properties.data = [0.6,0.0*Pa]
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-
-# Friction sensitivity solve used to compute the increment in slip
-# associated with changes in the Lagrange multiplier imposed by the
-# fault constitutive model.
-friction_pc_type = asm
-friction_sub_pc_factor_shift_type = nonzero
-friction_ksp_max_it = 25
-friction_ksp_gmres_restart = 30
-# Uncomment to view details of friction sensitivity solve.
-#friction_ksp_monitor = true
-#friction_ksp_view = true
-friction_ksp_converged_reason = true
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Set filenames for output.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/static_sliding.vtk
-
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/static_sliding-fault.vtk
-
-# Add slip rate to the default output fiels of slip and traction.
-vertex_data_fields=[slip,slip_rate,traction]
-
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/static_sliding-statevars.vtk
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,141 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# PROBLEM DESCRIPTION
-# ----------------------------------------------------------------------
-
-# Impose Dirichlet BC on the lateral ends of the bar to create uniform
-# shear. Friction on the fault controlled by a static friction fault
-# constitutive model is sufficient to keep the falt locked (no slip).
-
-# ----------------------------------------------------------------------
-# RUNNING THE SIMULATION
-# ----------------------------------------------------------------------
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for a static problem
-# with a locked fault controlled by a fault constitutive model.
-#
-# To run the simulation:
-# pylith static_stick.cfg
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-# Turn on progress for the dynamic fault.
-[pylithapp.journal.info]
-faultcohesivedyn = 1
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Change to an implicit time-stepping formulation and the nonlinear
-# solver.
-formulation = pylith.problems.Implicit
-formulation.solver = pylith.problems.SolverNonlinear
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and
-# 'all_nofault'.
-bc = [x_pos,x_neg,all_nofault]
-
-[pylithapp.timedependent.formulation.time_step]
-# Switch to a static simulation.
-total_time = 0.0*s
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc]
-# Switch from absorbing boundary condition to Dirichlet BC.
-x_pos = pylith.bc.DirichletBC
-x_neg = pylith.bc.DirichletBC
-
-[pylithapp.timedependent.bc.x_pos]
-# Dirichlet BC on +x face of bar
-
-# Name of nodeset in CUBIT.
-label = face_xpos
-
-# Pin the y-dof.
-bc_dof = [1]
-
-# Use a uniform DB to impose the Dirichlet BC. Impose 100 m of
-# displacement in the y direction.
-db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Dirichlet BC +x edge
-db_initial.values = [displacement-y]
-db_initial.data = [100.0*m]
-
-[pylithapp.timedependent.bc.x_neg]
-# Dirichlet BC on -x face of bar. Use the default (zero displacement) values.
-
-# Name of nodeset in CUBIT.
-label = face_xneg
-
-# Pin the y-dof.
-bc_dof = [1]
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces]
-# Change fault to dynamic fault interface.
-fault = pylith.faults.FaultCohesiveDyn
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Impose initial compressive normal tractions of -100 MPa. The initial
-# shear traction is zero.
-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*MPa, -100.0*MPa]
-
-# Use the static friction model.
-friction = pylith.friction.StaticFriction
-
-# Set static friction model parameters using a uniform DB. Set the
-# static coefficient of friction to 0.6 and cohesion to 0.0 Pa.
-friction.db_properties = spatialdata.spatialdb.UniformDB
-friction.db_properties.label = Static friction
-friction.db_properties.values = [friction-coefficient,cohesion]
-friction.db_properties.data = [0.6,0.0*Pa]
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-
-# Friction sensitivity solve used to compute the increment in slip
-# associated with changes in the Lagrange multiplier imposed by the
-# fault constitutive model. In this problem the fault remains locked,
-# so the fault constitutive model does not impose any changes in the
-# Lagrange multipliers, so the friction sensitivity solve will
-# converge in the first iteration.
-friction_pc_type = asm
-friction_sub_pc_factor_shift_type = nonzero
-friction_ksp_max_it = 25
-friction_ksp_gmres_restart = 30
-# Uncomment to view details of friction sensitivity solve.
-#friction_ksp_monitor = true
-#friction_ksp_view = true
-friction_ksp_converged_reason = true
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Set filenames for output.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/static_stick.vtk
-
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/static_stick-fault.vtk
-
-# Add slip rate to the default output fiels of slip and traction.
-vertex_data_fields=[slip,slip_rate,traction]
-
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/static_stick-statevars.vtk
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,11 +1,28 @@
# -*- Python -*-
+[pylithapp]
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
-[pylithapp]
+# Impose 1.0 m of left-lateral slip at time 0.1 s with slip rate
+# following Brune's far-field time function. The rise time (95% of the
+# final slip) is 2.0 s.
+#
+# The fault slip causes a shear wave to propagate down the bar. The
+# shear wave is absorbed at the lateral ends of the bar using the
+# absorbing dampers boundary condition.
# ----------------------------------------------------------------------
+# RUNNING THE SIMULATION
+# ----------------------------------------------------------------------
+
+# This is a self-contained simulation configuration file.
+#
+# To run the simulation:
+# pylith
+
+# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
[pylithapp.journal.info]
@@ -221,19 +238,20 @@
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-skip = 3
+vertex_data_fields = [displacement,velocity]
+skip = 2
writer.filename = output/shearwave.vtk
writer.time_format = %05.2f
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
+skip = 2
writer.filename = output/shearwave-fault.vtk
writer.time_format = %05.2f
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.pos.output]
-skip = 33
+skip = 29
writer.filename = output/shearwave-statevars-pos.vtk
writer.time_format = %05.2f
@@ -242,7 +260,7 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.neg.output]
-skip = 33
+skip = 29
writer.filename = output/shearwave-statevars-neg.vtk
writer.time_format = %05.2f
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg 2010-05-30 20:56:34 UTC (rev 16836)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg 2010-05-30 20:58:09 UTC (rev 16837)
@@ -1,11 +1,28 @@
# -*- Python -*-
+[pylithapp]
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
-[pylithapp]
+# Impose 1.0 m of left-lateral slip at time 0.1 s with slip rate
+# following Brune's far-field time function. The rise time (95% of the
+# final slip) is 2.0 s.
+#
+# The fault slip causes a shear wave to propagate down the bar. The
+# shear wave is absorbed at the lateral ends of the bar using the
+# absorbing dampers boundary condition.
# ----------------------------------------------------------------------
+# RUNNING THE SIMULATION
+# ----------------------------------------------------------------------
+
+# This is a self-contained simulation configuration file.
+#
+# To run the simulation:
+# pylith
+
+# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
[pylithapp.journal.info]
@@ -186,7 +203,7 @@
ksp_rtol = 1.0e-8
ksp_atol = 1.0e-10
ksp_max_it = 50
-ksp_gmres_restart = 10
+ksp_gmres_restart = 25
#ksp_monitor = true
#ksp_view = true
@@ -203,19 +220,20 @@
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-skip = 3
+vertex_data_fields = [displacement,velocity]
+skip = 2
writer.filename = output/shearwave.vtk
writer.time_format = %05.2f
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
+skip = 2
writer.filename = output/shearwave-fault.vtk
writer.time_format = %05.2f
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-skip = 33
+skip = 29
writer.filename = output/shearwave-statevars.vtk
writer.time_format = %05.2f
More information about the CIG-COMMITS
mailing list