[cig-commits] r22288 - short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4
brad at geodynamics.org
brad at geodynamics.org
Fri Jun 14 15:48:10 PDT 2013
Author: brad
Date: 2013-06-14 15:48:10 -0700 (Fri, 14 Jun 2013)
New Revision: 22288
Added:
short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/traction_change.timedb
Modified:
short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/pylithapp.cfg
short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg
short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg
short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg
Log:
Updated for more complex spontanoues rupture nucleation.
Modified: short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2013-06-14 22:45:52 UTC (rev 22287)
+++ short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2013-06-14 22:48:10 UTC (rev 22288)
@@ -44,9 +44,9 @@
# Set the spatial dimension of the problem.
dimension = 2
-# Change to an explicit time stepping formulation. We use explicit
-# time-stepping with a consistent (full) Jacobian matrix.
+# Change to an explicit time stepping formulation with no elastic prestep.
formulation = pylith.problems.Explicit
+elastic_prestep = False
# Nondimensionalize problem using wave propagation parameters.
normalizer = spatialdata.units.NondimElasticDynamic
Modified: short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg 2013-06-14 22:45:52 UTC (rev 22287)
+++ short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup.cfg 2013-06-14 22:48:10 UTC (rev 22288)
@@ -19,17 +19,31 @@
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)
+# Specify the initial tractions and a temporal perturbation on the
+# fault using a uniform DB.
+#
+# Initial
+# shear: 75.0 MPa (right-lateral)
+# normal 120 MPa (compressive)
+# Perturbation
+# shear: 25.0 MPa (right-lateral)
traction_perturbation = pylith.faults.TractPerturbation
[pylithapp.timedependent.interfaces.fault.traction_perturbation]
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Initial fault tractions
db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [-6.1*MPa, -10.0*MPa]
+db_initial.data = [-75.0*MPa, -120.0*MPa]
+db_change = spatialdata.spatialdb.UniformDB
+db_change.label = Traction perturbation
+db_change.values = [traction-shear,traction-normal,change-start-time]
+db_change.data = [-25.0*MPa, 0.0*MPa, 0.0*s]
+
+th_change = spatialdata.spatialdb.TimeHistory
+th_change.label = Traction perturbation time history
+th_change.filename = traction_change.timedb
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg 2013-06-14 22:45:52 UTC (rev 22287)
+++ short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_ratestateageing.cfg 2013-06-14 22:48:10 UTC (rev 22288)
@@ -44,13 +44,13 @@
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]
+friction.db_properties.data = [0.6,1.0e-6*m/s,0.02*m,0.008,0.012,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]
+friction.db_initial_state.data = [1.60623899921e+9*s]
# ----------------------------------------------------------------------
# output
@@ -60,7 +60,8 @@
writer.filename = output/ratestateageing.vtk
[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/ratestateageing-fault.vtk
+vertex_data_fields = [slip, slip_rate, traction, state_variable]
+writer.filename = output/ratestateageing-faut.vtk
[pylithapp.timedependent.materials.elastic.output]
writer.filename = output/ratestateageing-statevars.vtk
Modified: short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg 2013-06-14 22:45:52 UTC (rev 22287)
+++ short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/spontaneousrup_slipweakening.cfg 2013-06-14 22:48:10 UTC (rev 22288)
@@ -56,9 +56,7 @@
# default values.
vertex_info_fields = [strike_dir,normal_dir,traction_initial_value,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]
+vertex_data_fields = [slip,slip_rate,traction]
[pylithapp.timedependent.materials.elastic.output]
Added: short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/traction_change.timedb
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/traction_change.timedb (rev 0)
+++ short/3D/PyLith/branches/v1.7-trunk/examples/bar_shearwave/quad4/traction_change.timedb 2013-06-14 22:48:10 UTC (rev 22288)
@@ -0,0 +1,107 @@
+#TIME HISTORY ascii
+TimeHistory {
+ num-points = 102
+ time-units = s
+}
+ 0.000000e+00 0.000000e+00
+ 1.000000e-02 4.078037e-22
+ 2.000000e-02 2.932663e-11
+ 3.000000e-02 1.218502e-07
+ 4.000000e-02 7.849176e-06
+ 5.000000e-02 9.549752e-05
+ 6.000000e-02 5.049429e-04
+ 7.000000e-02 1.658363e-03
+ 8.000000e-02 4.044422e-03
+ 9.000000e-02 8.088308e-03
+ 1.000000e-01 1.407778e-02
+ 1.100000e-01 2.214807e-02
+ 1.200000e-01 3.230207e-02
+ 1.300000e-01 4.444363e-02
+ 1.400000e-01 5.841131e-02
+ 1.500000e-01 7.400641e-02
+ 1.600000e-01 9.101426e-02
+ 1.700000e-01 1.092190e-01
+ 1.800000e-01 1.284131e-01
+ 1.900000e-01 1.484033e-01
+ 2.000000e-01 1.690133e-01
+ 2.100000e-01 1.900854e-01
+ 2.200000e-01 2.114797e-01
+ 2.300000e-01 2.330741e-01
+ 2.400000e-01 2.547623e-01
+ 2.500000e-01 2.764530e-01
+ 2.600000e-01 2.980681e-01
+ 2.700000e-01 3.195412e-01
+ 2.800000e-01 3.408162e-01
+ 2.900000e-01 3.618464e-01
+ 3.000000e-01 3.825927e-01
+ 3.100000e-01 4.030232e-01
+ 3.200000e-01 4.231117e-01
+ 3.300000e-01 4.428374e-01
+ 3.400000e-01 4.621836e-01
+ 3.500000e-01 4.811377e-01
+ 3.600000e-01 4.996902e-01
+ 3.700000e-01 5.178344e-01
+ 3.800000e-01 5.355658e-01
+ 3.900000e-01 5.528822e-01
+ 4.000000e-01 5.697828e-01
+ 4.100000e-01 5.862684e-01
+ 4.200000e-01 6.023408e-01
+ 4.300000e-01 6.180030e-01
+ 4.400000e-01 6.332586e-01
+ 4.500000e-01 6.481120e-01
+ 4.600000e-01 6.625682e-01
+ 4.700000e-01 6.766324e-01
+ 4.800000e-01 6.903103e-01
+ 4.900000e-01 7.036079e-01
+ 5.000000e-01 7.165313e-01
+ 5.100000e-01 7.290867e-01
+ 5.200000e-01 7.412804e-01
+ 5.300000e-01 7.531188e-01
+ 5.400000e-01 7.646082e-01
+ 5.500000e-01 7.757550e-01
+ 5.600000e-01 7.865654e-01
+ 5.700000e-01 7.970456e-01
+ 5.800000e-01 8.072018e-01
+ 5.900000e-01 8.170397e-01
+ 6.000000e-01 8.265654e-01
+ 6.100000e-01 8.357846e-01
+ 6.200000e-01 8.447028e-01
+ 6.300000e-01 8.533254e-01
+ 6.400000e-01 8.616578e-01
+ 6.500000e-01 8.697051e-01
+ 6.600000e-01 8.774721e-01
+ 6.700000e-01 8.849638e-01
+ 6.800000e-01 8.921848e-01
+ 6.900000e-01 8.991395e-01
+ 7.000000e-01 9.058323e-01
+ 7.100000e-01 9.122673e-01
+ 7.200000e-01 9.184485e-01
+ 7.300000e-01 9.243798e-01
+ 7.400000e-01 9.300647e-01
+ 7.500000e-01 9.355070e-01
+ 7.600000e-01 9.407098e-01
+ 7.700000e-01 9.456765e-01
+ 7.800000e-01 9.504101e-01
+ 7.900000e-01 9.549135e-01
+ 8.000000e-01 9.591895e-01
+ 8.100000e-01 9.632406e-01
+ 8.200000e-01 9.670695e-01
+ 8.300000e-01 9.706784e-01
+ 8.400000e-01 9.740695e-01
+ 8.500000e-01 9.772450e-01
+ 8.600000e-01 9.802067e-01
+ 8.700000e-01 9.829564e-01
+ 8.800000e-01 9.854958e-01
+ 8.900000e-01 9.878265e-01
+ 9.000000e-01 9.899498e-01
+ 9.100000e-01 9.918671e-01
+ 9.200000e-01 9.935795e-01
+ 9.300000e-01 9.950880e-01
+ 9.400000e-01 9.963935e-01
+ 9.500000e-01 9.974969e-01
+ 9.600000e-01 9.983987e-01
+ 9.700000e-01 9.990996e-01
+ 9.800000e-01 9.995999e-01
+ 9.900000e-01 9.999000e-01
+ 1.000000e+00 1.000000e+00
+ 1.500000e+01 1.000000e+00
More information about the CIG-COMMITS
mailing list