[cig-commits] r16712 - short/3D/PyLith/trunk/examples/bar_shearwave/quad4
brad at geodynamics.org
brad at geodynamics.org
Fri May 14 22:02:34 PDT 2010
Author: brad
Date: 2010-05-14 22:02:33 -0700 (Fri, 14 May 2010)
New Revision: 16712
Added:
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
Removed:
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault_shearwave.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.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/static-sliding.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-stick.cfg
Modified:
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
Log:
More work on quad4 shear wave examples.
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am 2010-05-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am 2010-05-15 05:02:33 UTC (rev 16712)
@@ -20,15 +20,15 @@
matprops.spatialdb \
mesh_quad4_200m.jou \
pylithapp.cfg \
- kinematicfault.cfg \
+ shearwave_kinematicslip.cfg \
kinematic_slip.spatialdb \
kinematic_risetime.spatialdb \
kinematic_sliptime.spatialdb \
dynamicfault.cfg \
- static-stick.cfg \
- static-sliding.cfg \
- dynamicfault_shearwave.cfg \
- shearwave-slipweakening.cfg \
- shearwave-ratestateageing.cfg
+ static_stick.cfg \
+ static_sliding.cfg \
+ shearwave_staticfriction.cfg \
+ shearwave_slipweakening.cfg \
+ shearwave_ratestateageing.cfg
# End of file
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -28,6 +28,7 @@
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
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault_shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault_shearwave.cfg 2010-05-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault_shearwave.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -1,124 +0,0 @@
-# -*- Python -*-
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-formulation = pylith.problems.Explicit
-formulation.solver = pylith.problems.SolverNonlinear
-#formulation = pylith.problems.ExplicitLumped
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-#normalizer.mass_density = 1000.0*kg/m**3
-#normalizer.shear_wave_speed = 1.0*m/s
-#normalizer.wave_period = 1.0*s
-
-bc = [x_pos,x_neg,all_nofault]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.05*s
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing BC on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# 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
-
-# Nodeset id from CUBIT
-label = 12
-
-# 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 BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-
-fault = pylith.faults.FaultCohesiveDyn
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# 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
-
-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]
-
-friction = pylith.friction.StaticFriction
-
-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
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-vertex_data_fields = [displacement,velocity]
-writer.filename = output/shearwave.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-vertex_data_fields = [slip,slip_rate,traction]
-writer.filename = output/shearwave-fault.vtk
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/shearwave-statevars.vtk
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg 2010-05-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -1,98 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# 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.
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-faultcohesivekin = 1
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-dimension = 2
-
-# Uncomment the line below to switch to explicit time-stepping with a
-# lumped Jacobian matrix.
-#formulation = pylith.problems.ExplicitLumped
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing boundary condition on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# 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
-
-# Nodeset id from CUBIT
-label = 12
-
-# 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/kinematicfault.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/kinematicfault-fault.vtk
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/kinematicfault-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-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave-ratestateageing.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -1,127 +0,0 @@
-# -*- Python -*-
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-formulation = pylith.problems.Explicit
-formulation.solver = pylith.problems.SolverNonlinear
-#formulation = pylith.problems.ExplicitLumped
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-#normalizer.mass_density = 1000.0*kg/m**3
-#normalizer.shear_wave_speed = 1.0*m/s
-#normalizer.wave_period = 1.0*s
-
-bc = [x_pos,x_neg,all_nofault]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.05*s
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing BC on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# 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
-
-# Nodeset id from CUBIT
-label = 12
-
-# 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 BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-
-fault = pylith.faults.FaultCohesiveDyn
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# 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
-
-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]
-
-friction = pylith.friction.RateStateAgeing
-
-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,0.000001*m/s,0.0370*m,0.0125,0.0172,0.0*Pa]
-
-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
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/shearwave-ratestateageing.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/shearwave-ratestateageing-fault.vtk
-
-# Give basename for VTK output of state variables.
-[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-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave-slipweakening.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -1,124 +0,0 @@
-# -*- Python -*-
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-formulation = pylith.problems.Explicit
-formulation.solver = pylith.problems.SolverNonlinear
-#formulation = pylith.problems.ExplicitLumped
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-#normalizer.mass_density = 1000.0*kg/m**3
-#normalizer.shear_wave_speed = 1.0*m/s
-#normalizer.wave_period = 1.0*s
-
-bc = [x_pos,x_neg,all_nofault]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 1.0*s
-dt = 0.05*s
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing BC on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# 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
-
-# Nodeset id from CUBIT
-label = 12
-
-# 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 BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-
-fault = pylith.faults.FaultCohesiveDyn
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# 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
-
-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]
-
-friction = pylith.friction.SlipWeakening
-
-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
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/shearwave-slipweakening.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-vertex_info_fields = [strike_dir,normal_dir,initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
-vertex_data_fields = [slip,traction,cumulative_slip,previous_slip]
-writer.filename = output/shearwave-slipweakening-fault.vtk
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/shearwave-slipweakening-statevars.vtk
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg (from rev 16711, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_kinematicslip.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -0,0 +1,108 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# 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
+#
+# Output will be directed to the directory output.
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+
+# Uncomment the line below to switch to explicit time-stepping with a
+# lumped Jacobian matrix.
+#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
+
+# Nodeset id from CUBIT
+label = 11
+
+# 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
+
+# Nodeset id from CUBIT
+label = 12
+
+# 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
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg (from rev 16708, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave-ratestateageing.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -0,0 +1,127 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+formulation = pylith.problems.Explicit
+formulation.solver = pylith.problems.SolverNonlinear
+#formulation = pylith.problems.ExplicitLumped
+
+# Nondimensionalize problem using wave propagation parameters.
+normalizer = spatialdata.units.NondimElasticDynamic
+#normalizer.mass_density = 1000.0*kg/m**3
+#normalizer.shear_wave_speed = 1.0*m/s
+#normalizer.wave_period = 1.0*s
+
+bc = [x_pos,x_neg,all_nofault]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 12.0*s
+dt = 0.05*s
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing BC on +x face of bar
+
+# Nodeset id from CUBIT
+label = 11
+
+# 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
+
+# Nodeset id from CUBIT
+label = 12
+
+# 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 BC on all vertices except fault vertices
+
+# Pin x and z degrees of freedom
+bc_dof = [0]
+
+# Nodeset id from CUBIT
+label = 13
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+
+# Provide information on the fault (interface).
+[pylithapp.timedependent.interfaces]
+
+fault = pylith.faults.FaultCohesiveDyn
+
+[pylithapp.timedependent.interfaces.fault]
+
+# Material id
+id = 100
+
+# Nodeset id from CUBIT
+label = 10
+
+# 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
+
+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]
+
+friction = pylith.friction.RateStateAgeing
+
+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,0.000001*m/s,0.0370*m,0.0125,0.0172,0.0*Pa]
+
+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
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/shearwave-ratestateageing.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/shearwave-ratestateageing-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/shearwave-ratestateageing-statevars.vtk
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg (from rev 16708, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave-slipweakening.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -0,0 +1,124 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+formulation = pylith.problems.Explicit
+formulation.solver = pylith.problems.SolverNonlinear
+#formulation = pylith.problems.ExplicitLumped
+
+# Nondimensionalize problem using wave propagation parameters.
+normalizer = spatialdata.units.NondimElasticDynamic
+#normalizer.mass_density = 1000.0*kg/m**3
+#normalizer.shear_wave_speed = 1.0*m/s
+#normalizer.wave_period = 1.0*s
+
+bc = [x_pos,x_neg,all_nofault]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 1.0*s
+dt = 0.05*s
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing BC on +x face of bar
+
+# Nodeset id from CUBIT
+label = 11
+
+# 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
+
+# Nodeset id from CUBIT
+label = 12
+
+# 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 BC on all vertices except fault vertices
+
+# Pin x and z degrees of freedom
+bc_dof = [0]
+
+# Nodeset id from CUBIT
+label = 13
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+
+# Provide information on the fault (interface).
+[pylithapp.timedependent.interfaces]
+
+fault = pylith.faults.FaultCohesiveDyn
+
+[pylithapp.timedependent.interfaces.fault]
+
+# Material id
+id = 100
+
+# Nodeset id from CUBIT
+label = 10
+
+# 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
+
+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]
+
+friction = pylith.friction.SlipWeakening
+
+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
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/shearwave-slipweakening.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
+vertex_data_fields = [slip,traction,cumulative_slip,previous_slip]
+writer.filename = output/shearwave-slipweakening-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/shearwave-slipweakening-statevars.vtk
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg (from rev 16708, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault_shearwave.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -0,0 +1,103 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# 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.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+# Switch to a lumped Jacobian matrix.
+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 BC on +x face of bar
+
+# Nodeset id from CUBIT
+label = 11
+
+# 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
+
+# Nodeset id from CUBIT
+label = 12
+
+# 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 BC on all vertices except fault vertices
+
+# Pin x and z degrees of freedom
+bc_dof = [0]
+
+# Nodeset id from CUBIT
+label = 13
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+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]
+
+friction = pylith.friction.StaticFriction
+
+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-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-sliding.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -1,84 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for a static problem
-# using a fault with a fault constitutive model.
-
-# ----------------------------------------------------------------------
-# 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 'y_pos'.
-bc = [x_pos,x_neg,all_nofault]
-
-[pylithapp.timedependent.formulation.time_step]
-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
-
-# Nodeset id from CUBIT
-label = 11
-
-bc_dof = [1]
-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,32.0*km]
-
-[pylithapp.timedependent.bc.x_neg]
-# Dirichlet BC on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-bc_dof = [1]
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------x
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces.fault]
-
-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, -10000.0*MPa]
-
-friction = pylith.friction.StaticFriction
-
-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*m,0.0*Pa]
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/static-sliding.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/static-sliding-fault.vtk
-
-# Give basename for VTK output of state variables.
-[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-15 04:30:22 UTC (rev 16711)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-stick.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -1,83 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies the specific parameters for a static problem
-# using a fault with a fault constitutive model.
-
-# ----------------------------------------------------------------------
-# 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 'y_pos'.
-bc = [x_pos,x_neg,all_nofault]
-
-[pylithapp.timedependent.formulation.time_step]
-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
-
-# Nodeset id from CUBIT
-label = 11
-
-bc_dof = [1]
-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,10.0*km]
-
-[pylithapp.timedependent.bc.x_neg]
-# Dirichlet BC on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-bc_dof = [1]
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-
-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, -10000.0*MPa]
-
-friction = pylith.friction.StaticFriction
-
-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*m,0.0*Pa]
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/static-stick.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/static-stick-fault.vtk
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/static-stick-statevars.vtk
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg (from rev 16711, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-sliding.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -0,0 +1,93 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# 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 dynamicfault.cfg static_sliding.cfg
+#
+# Output will be directed to the directory output.
+
+# ----------------------------------------------------------------------
+# 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 'y_pos'.
+bc = [x_pos,x_neg,all_nofault]
+
+[pylithapp.timedependent.formulation.time_step]
+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
+
+# Nodeset id from CUBIT
+label = 11
+
+bc_dof = [1]
+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
+
+# Nodeset id from CUBIT
+label = 12
+
+bc_dof = [1]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------x
+# Provide information on the fault (interface).
+[pylithapp.timedependent.interfaces.fault]
+
+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]
+
+friction = pylith.friction.StaticFriction
+
+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*m,0.0*Pa]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/static_sliding.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/static_sliding-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/static_sliding-statevars.vtk
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg (from rev 16711, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-stick.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg 2010-05-15 05:02:33 UTC (rev 16712)
@@ -0,0 +1,92 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
+# 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 dynamicfault.cfg static-stick.cfg
+#
+# Output will be directed to the directory output.
+
+# ----------------------------------------------------------------------
+# 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 'y_pos'.
+bc = [x_pos,x_neg,all_nofault]
+
+[pylithapp.timedependent.formulation.time_step]
+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
+
+# Nodeset id from CUBIT
+label = 11
+
+bc_dof = [1]
+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,100.0*m]
+
+[pylithapp.timedependent.bc.x_neg]
+# Dirichlet BC on -x face of bar
+
+# Nodeset id from CUBIT
+label = 12
+
+bc_dof = [1]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+
+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]
+
+friction = pylith.friction.StaticFriction
+
+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*m,0.0*Pa]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/static_stick.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/static_stick-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/static_stick-statevars.vtk
More information about the CIG-COMMITS
mailing list