[cig-commits] r16572 - short/3D/PyLith/trunk/tests/2d/quad4
surendra at geodynamics.org
surendra at geodynamics.org
Wed Apr 21 13:08:52 PDT 2010
Author: surendra
Date: 2010-04-21 13:08:52 -0700 (Wed, 21 Apr 2010)
New Revision: 16572
Modified:
short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionCompression.py
short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionOpening.py
short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py
short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py
short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningCompression.py
short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningOpening.py
short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearSliding.py
short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearStick.py
short/3D/PyLith/trunk/tests/2d/quad4/friction_compression.cfg
short/3D/PyLith/trunk/tests/2d/quad4/friction_opening.cfg
short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg
short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg
short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_compression.cfg
short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_opening.cfg
short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_sliding.cfg
short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_stick.cfg
Log:
Finished tests for .vtk files of slip-weakening friction. Added checks for friction parameters which are hooked to fault output. Updated corresponding static-friction tests also to check friction parameters
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionCompression.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionCompression.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionCompression.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","friction_coefficient","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -128,6 +128,7 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
+ frictionCoefficient = 0.6
nvertices = self.faultMesh['nvertices']
@@ -146,6 +147,13 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "friction_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = frictionCoefficient
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionOpening.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionOpening.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionOpening.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","friction_coefficient","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -129,6 +129,7 @@
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
slip = (0.0, 1.0)
+ frictionCoefficient = 0.6
nvertices = self.faultMesh['nvertices']
@@ -147,6 +148,13 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "friction_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = frictionCoefficient
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
field[:,0] = slip[0]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","friction_coefficient","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -128,6 +128,7 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
+ frictionCoefficient = 0.6
uy_l = 1.0
len = 8000.0
@@ -155,6 +156,13 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "friction_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = frictionCoefficient
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
field[:,0] = slip[0]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","friction_coefficient","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -128,6 +128,7 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+7)
+ frictionCoefficient = 0.6
nvertices = self.faultMesh['nvertices']
@@ -146,6 +147,13 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "friction_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = frictionCoefficient
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningCompression.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningCompression.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningCompression.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","static_coefficient","dynamic_coefficient","slip_weakening_parameter","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -91,7 +91,7 @@
return
filename = "%s-fault_t0000000.vtk" % self.outputRoot
- fields = ["slip", "traction"]
+ fields = ["slip", "traction","cumulative_slip","previous_slip"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -128,6 +128,9 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
+ staticCoefficient = 0.6
+ dynamicCoefficient = 0.59
+ slipWeakeningParameter = 0.2
nvertices = self.faultMesh['nvertices']
@@ -146,6 +149,21 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "static_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = staticCoefficient
+
+ elif name == "dynamic_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = dynamicCoefficient
+
+ elif name == "slip_weakening_parameter":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = slipWeakeningParameter
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
@@ -154,6 +172,12 @@
field[:,0] = 0.0
field[:,1] = -2.2e+6
+ elif name == "cumulative_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
+ elif name == "previous_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
else:
raise ValueError("Unknown fault field '%s'." % name)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningOpening.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningOpening.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningOpening.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","static_coefficient","dynamic_coefficient","slip_weakening_parameter","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -91,7 +91,7 @@
return
filename = "%s-fault_t0000000.vtk" % self.outputRoot
- fields = ["slip", "traction"]
+ fields = ["slip", "traction","cumulative_slip","previous_slip"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -129,6 +129,9 @@
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
slip = (0.0, 1.0)
+ staticCoefficient = 0.6
+ dynamicCoefficient = 0.59
+ slipWeakeningParameter = 0.2
nvertices = self.faultMesh['nvertices']
@@ -147,6 +150,21 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "static_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = staticCoefficient
+
+ elif name == "dynamic_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = dynamicCoefficient
+
+ elif name == "slip_weakening_parameter":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = slipWeakeningParameter
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
field[:,0] = slip[0]
@@ -155,6 +173,12 @@
elif name == "traction":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
+ elif name == "cumulative_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
+ elif name == "previous_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
else:
raise ValueError("Unknown fault field '%s'." % name)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearSliding.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearSliding.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearSliding.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","static_coefficient","dynamic_coefficient","slip_weakening_parameter","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -91,7 +91,7 @@
return
filename = "%s-fault_t0000000.vtk" % self.outputRoot
- fields = ["slip", "traction"]
+ fields = ["slip", "traction","cumulative_slip","previous_slip"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -128,6 +128,9 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
+ staticCoefficient = 0.6
+ dynamicCoefficient = 0.59
+ slipWeakeningParameter = 0.2
uy_l = 1.0
len = 8000.0
@@ -137,6 +140,7 @@
p_mu = p_density*p_vs**2
D = uy_l - len * sigma_f/p_mu
slip = (D, 0.0)
+ cumulativeSlip = D
nvertices = self.faultMesh['nvertices']
@@ -155,6 +159,21 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "static_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = staticCoefficient
+
+ elif name == "dynamic_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = dynamicCoefficient
+
+ elif name == "slip_weakening_parameter":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = slipWeakeningParameter
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
field[:,0] = slip[0]
@@ -165,6 +184,13 @@
field[:,0] = 0.6e+6
field[:,1] = -1.0e+6
+ elif name == "cumulative_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = cumulativeSlip
+
+ elif name == "previous_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
else:
raise ValueError("Unknown fault field '%s'." % name)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearStick.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearStick.py 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestSlipWeakeningShearStick.py 2010-04-21 20:08:52 UTC (rev 16572)
@@ -77,7 +77,7 @@
return
filename = "%s-fault_info.vtk" % self.outputRoot
- fields = ["strike_dir", "normal_dir", "initial_traction"]
+ fields = ["strike_dir", "normal_dir", "initial_traction","static_coefficient","dynamic_coefficient","slip_weakening_parameter","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -91,7 +91,7 @@
return
filename = "%s-fault_t0000000.vtk" % self.outputRoot
- fields = ["slip", "traction"]
+ fields = ["slip", "traction","cumulative_slip","previous_slip"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
@@ -128,6 +128,9 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+7)
+ staticCoefficient = 0.6
+ dynamicCoefficient = 0.59
+ slipWeakeningParameter = 0.2
nvertices = self.faultMesh['nvertices']
@@ -146,6 +149,21 @@
field[:,0] = initialTraction[0]
field[:,1] = initialTraction[1]
+ elif name == "static_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = staticCoefficient
+
+ elif name == "dynamic_coefficient":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = dynamicCoefficient
+
+ elif name == "slip_weakening_parameter":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+ field[:] = slipWeakeningParameter
+
+ elif name == "cohesion":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
@@ -154,6 +172,12 @@
field[:,0] = 1.0e+6
field[:,1] = -1.0e+7
+ elif name == "cumulative_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
+ elif name == "previous_slip":
+ field = numpy.zeros( (nvertices, 1), dtype=numpy.float64)
+
else:
raise ValueError("Unknown fault field '%s'." % name)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_compression.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_compression.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_compression.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -168,7 +168,7 @@
# Give basename for VTK fault output.
[friction_compression.timedependent.interfaces.fault.output]
writer.filename = friction_compression-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction,friction_coefficient,cohesion]
# Give basename for VTK output of state variables.
[friction_compression.timedependent.materials.material.output]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_opening.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_opening.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_opening.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -195,7 +195,7 @@
# Give basename for VTK fault output.
[friction_opening.timedependent.interfaces.fault.output]
writer.filename = friction_opening-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction,friction_coefficient,cohesion]
# Give basename for VTK output of state variables.
[friction_opening.timedependent.materials.material.output]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -231,7 +231,7 @@
# Give basename for VTK fault output.
[friction_shear_sliding.timedependent.interfaces.fault.output]
writer.filename = friction_shear_sliding-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction,friction_coefficient,cohesion]
# Give basename for VTK output of state variables.
[friction_shear_sliding.timedependent.materials.material.output]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -245,7 +245,7 @@
# Give basename for VTK fault output.
[friction_shear_stick.timedependent.interfaces.fault.output]
writer.filename = friction_shear_stick-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction,friction_coefficient,cohesion]
# Give basename for VTK output of state variables.
[friction_shear_stick.timedependent.materials.material.output]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_compression.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_compression.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_compression.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -132,7 +132,7 @@
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.61,0.6,0.2*m,0.0*Pa]
+friction.db_properties.data = [0.6,0.59,0.2*m,0.0*Pa]
# ----------------------------------------------------------------------
# PETSc
@@ -169,15 +169,11 @@
# Give basename for VTK fault output.
[slipweakening_compression.timedependent.interfaces.fault.output]
writer.filename = slipweakening_compression-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+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]
# Give basename for VTK output of state variables.
[slipweakening_compression.timedependent.materials.material.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = slipweakening_compression-elastic.vtk
-# Give basename for VTK friction output.
-[slipweakening_compression.timedependent.interfaces.fault.friction.output]
-writer.filename = slipweakening_compression-statevars.vtk
-vertex_info_fields = [static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
-
Modified: short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_opening.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_opening.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_opening.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -155,10 +155,12 @@
db_initial_tractions.values = [traction-shear,traction-normal]
db_initial_tractions.data = [0.0*Pa, -1.0*MPa]
+friction = pylith.friction.SlipWeakening
+
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]
+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.59,0.2*m,0.0*Pa]
# ----------------------------------------------------------------------
# PETSc
@@ -195,7 +197,8 @@
# Give basename for VTK fault output.
[slipweakening_opening.timedependent.interfaces.fault.output]
writer.filename = slipweakening_opening-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+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]
# Give basename for VTK output of state variables.
[slipweakening_opening.timedependent.materials.material.output]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_sliding.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_sliding.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -191,10 +191,12 @@
db_initial_tractions.values = [traction-shear,traction-normal]
db_initial_tractions.data = [0.0*Pa, -1.0*MPa]
+friction = pylith.friction.SlipWeakening
+
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]
+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.59,0.2*m,0.0*Pa]
# ----------------------------------------------------------------------
# PETSc
@@ -231,7 +233,8 @@
# Give basename for VTK fault output.
[slipweakening_shear_sliding.timedependent.interfaces.fault.output]
writer.filename = slipweakening_shear_sliding-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+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]
# Give basename for VTK output of state variables.
[slipweakening_shear_sliding.timedependent.materials.material.output]
Modified: short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_stick.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_stick.cfg 2010-04-21 19:16:50 UTC (rev 16571)
+++ short/3D/PyLith/trunk/tests/2d/quad4/slipweakening_shear_stick.cfg 2010-04-21 20:08:52 UTC (rev 16572)
@@ -204,10 +204,12 @@
db_initial_tractions.values = [traction-shear,traction-normal]
db_initial_tractions.data = [0.0*Pa, -10.0*MPa]
+friction = pylith.friction.SlipWeakening
+
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]
+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.59,0.2*m,0.0*Pa]
# ----------------------------------------------------------------------
# PETSc
@@ -245,7 +247,8 @@
# Give basename for VTK fault output.
[slipweakening_shear_stick.timedependent.interfaces.fault.output]
writer.filename = slipweakening_shear_stick-fault.vtk
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+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]
# Give basename for VTK output of state variables.
[slipweakening_shear_stick.timedependent.materials.material.output]
More information about the CIG-COMMITS
mailing list