[cig-commits] r16539 - in short/3D/PyLith/trunk: libsrc/friction playpen/friction/bar_shearwave/quad4 tests/2d/quad4

Mon Apr 12 19:18:08 PDT 2010

Author: surendra
Date: 2010-04-12 19:18:08 -0700 (Mon, 12 Apr 2010)
New Revision: 16539

Added:
   short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave-ratestateageing.cfg
Modified:
   short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc
   short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/pylithapp.cfg
   short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou
Log:
Added full-scale test for Rate and State friction with Ageing Law

Modified: short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc
===================================================================

--- short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc	2010-04-12 23:32:41 UTC (rev 16538)
+++ short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc	2010-04-13 02:18:08 UTC (rev 16539)
@@ -349,7 +349,7 @@
   stateVars[s_state] = thetaTpdtVertex;
 
   std::cout << "STATEVAR before: " << thetaTVertex << ", after: " << thetaTpdtVertex << std::endl;
-
+  std::cout << "TIME-STEP: " << dt << std::endl;
   PetscLogFlops(6);
 } // _updateStateVars
 

Modified: short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/pylithapp.cfg	2010-04-12 23:32:41 UTC (rev 16538)
+++ short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/pylithapp.cfg	2010-04-13 02:18:08 UTC (rev 16539)
@@ -7,18 +7,18 @@
 # ----------------------------------------------------------------------
 [pylithapp.journal.info]
 timedependent = 1
-explicit = 1
-petsc = 1
-meshiocubit = 1
-explicitelasticity = 1
+#explicit = 1
+#petsc = 1
+#meshiocubit = 1
+#explicitelasticity = 1
 #quadrature2d = 1
 #fiatlagrange = 1
-faultcohesivedyn = 1
+#faultcohesivedyn = 1
 
 [pylithapp.journal.debug]
-pylithapp = 1
-problem = 1
-explicit = 1
+#pylithapp = 1
+#problem = 1
+#explicit = 1
 
 # ----------------------------------------------------------------------
 # mesh_generator
@@ -85,14 +85,14 @@
 ksp_gmres_restart = 250
 snes_max_it = 300
 
-ksp_monitor = true
-ksp_view = true
+#ksp_monitor = true
+#ksp_view = true
 ksp_converged_reason = true
 
 snes_rtol = 1.0e-8
 snes_atol = 1.0e-15
 snes_monitor = true
-snes_view = true
+#snes_view = true
 snes_converged_reason = true
 
 #log_summary = true

Added: short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave-ratestateageing.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave-ratestateageing.cfg	                        (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave-ratestateageing.cfg	2010-04-13 02:18:08 UTC (rev 16539)
@@ -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 = [-7500.0*MPa, -10000.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

Modified: short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou	2010-04-12 23:32:41 UTC (rev 16538)
+++ short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou	2010-04-13 02:18:08 UTC (rev 16539)
@@ -75,3 +75,4 @@
 # Export exodus file
 # ----------------------------------------------------------------------
 export mesh "mesh.exo" dimension 2 overwrite
+

