[cig-commits] r18299 - in short/3D/PyLith/branches/v1.5-stable/tests/2d: . frictionslide
brad at geodynamics.org
brad at geodynamics.org
Thu Apr 28 16:08:12 PDT 2011
Author: brad
Date: 2011-04-28 16:08:12 -0700 (Thu, 28 Apr 2011)
New Revision: 18299
Added:
short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/
short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg
short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/quad4.mesh
short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/ratestate_stable.cfg
Log:
Started work on friction sliding test.
Added: short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg 2011-04-28 23:08:12 UTC (rev 18299)
@@ -0,0 +1,169 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+petsc = 1
+meshiocubit = 1
+explicitelasticity = 1
+quadrature2d = 1
+fiatlagrange = 1
+faultcohesivedyn = 1
+
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+explicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+reader = pylith.meshio.MeshIOAscii
+
+[pylithapp.mesh_generator.reader]
+filename = quad4.mesh
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+
+formulation = pylith.problems.Implicit
+
+normalizer = spatialdata.units.NondimElasticQuasistatic
+normalizer.length_scale = 1.0*m
+normalizer.relaxation_time = 1.0*s
+
+bc = [top,bottom]
+interfaces = [fault]
+materials = [elastic]
+
+[pylithapp.timedependent.implicit]
+solver = pylith.problems.SolverNonlinear
+
+[pylithapp.timedependent.implicit.time_step]
+total_time = 12.0*s
+dt = 0.05*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 0
+db_properties = spatialdata.spatialdb.UniformDB
+db_properties.label = Elastic properties
+db_properties.values = [density, vp, vs]
+db_properties.data = [2500.0*km/m**3, 5.1962*km/s, 3.0*km/s]
+
+
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.top]
+bc_dof = [0, 1]
+label = top
+
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on +y
+db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
+db_rate.data = [1.0e-6*m/s,0.0*m/s,0.0*s]
+
+[pylithapp.timedependent.bc.bottom]
+bc_dof = [0, 1]
+label = bottom
+
+db_rate = spatialdata.spatialdb.UniformDB
+db_rate.label = Dirichlet rate BC on -y
+db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
+db_rate.data = [-1.0e-6*m/s,0.0*m/s,0.0*s]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveDyn
+
+[pylithapp.timedependent.interfaces.fault]
+id = 100
+label = 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]
+
+
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+
+# Preconditioner settings.
+pc_type = asm
+sub_pc_factor_shift_type = nonzero
+
+# KSP
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-12
+ksp_max_it = 500
+ksp_gmres_restart = 100
+
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+
+# SNES
+snes_rtol = 1.0e-8
+snes_atol = 1.0e-12
+snes_max_it = 500
+
+snes_monitor = true
+snes_view = true
+snes_converged_reason = true
+
+# Friction
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
+#log_summary = true
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+vertex_data_fields = [displacement,velocity]
+writer.time_format = %05.2f
+
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_data_fields=[slip,slip_rate,traction]
+writer.time_format = %05.2f
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_info_fields = []
+cell_data_fields = []
+cell_filter = pylith.meshio.CellFilterAvgMesh
Added: short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/quad4.mesh
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/quad4.mesh (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/quad4.mesh 2011-04-28 23:08:12 UTC (rev 18299)
@@ -0,0 +1,85 @@
+// Original mesh
+//
+// 2 ----- 5 ----- 8 -----11 -----14
+// | | | | |
+// | 4 | 5 | 6 | 7 |
+// | | | | |
+// 1 ----- 4 ----- 7 -----10 -----13
+// | | | | |
+// | 0 | 1 | 2 | 3 |
+// | | | | |
+// 0 ----- 3 ----- 6 ----- 9 -----12
+//
+mesh = {
+ dimension = 2
+ use-index-zero = true
+ vertices = {
+ dimension = 2
+ count = 15
+ coordinates = {
+ 0 -2.0 -1.0
+ 1 -2.0 0.0
+ 2 -2.0 1.0
+ 3 -1.0 -1.0
+ 4 -1.0 0.0
+ 5 -1.0 1.0
+ 6 0.0 -1.0
+ 7 0.0 0.0
+ 8 0.0 1.0
+ 9 +1.0 -1.0
+ 10 +1.0 0.0
+ 11 +1.0 1.0
+ 12 +2.0 -1.0
+ 13 +2.0 0.0
+ 14 +2.0 1.0
+ }
+ }
+ cells = {
+ count = 8
+ num-corners = 4
+ simplices = {
+ 0 0 3 4 1
+ 1 3 6 7 4
+ 2 6 9 10 7
+ 3 9 12 13 10
+ 4 1 4 5 2
+ 5 4 7 8 5
+ 6 7 10 11 8
+ 7 10 13 14 11
+ }
+ material-ids = {
+ 0 0
+ 1 0
+ 2 0
+ 3 0
+ 4 0
+ 5 0
+ 6 0
+ 7 0
+ }
+ }
+ group = {
+ name = fault
+ type = vertices
+ count = 5
+ indices = {
+ 1 4 7 10 13
+ }
+ }
+ group = {
+ name = top
+ type = vertices
+ count = 5
+ indices = {
+ 2 5 8 11 14
+ }
+ }
+ group = {
+ name = bottom
+ type = vertices
+ count = 5
+ indices = {
+ 0 3 6 9 12
+ }
+ }
+}
Added: short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/ratestate_stable.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/ratestate_stable.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/ratestate_stable.cfg 2011-04-28 23:08:12 UTC (rev 18299)
@@ -0,0 +1,45 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+# Rate-state friction with stable sliding.
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+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.01 m
+# a: 0.016
+# b: 0.012
+# 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.01*m,0.016,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]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Set filenames for output.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/ratestate_stable.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/ratestate_stable-fault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/ratestate_stable-statevars.vtk
More information about the CIG-COMMITS
mailing list