[cig-commits] r18345 - short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide

[brad at geodynamics.org]

Author: brad
Date: 2011-05-11 16:45:17 -0700 (Wed, 11 May 2011)
New Revision: 18345

Modified:
   short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg
Log:
Adjustments to friction sliding test.

Modified: short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg	2011-05-11 23:41:56 UTC (rev 18344)
+++ short/3D/PyLith/branches/v1.5-stable/tests/2d/frictionslide/pylithapp.cfg	2011-05-11 23:45:17 UTC (rev 18345)
@@ -7,10 +7,10 @@
 # Turn on some journals to show progress.
 [pylithapp.journal.info]
 timedependent = 1
-explicit = 1
+implicit = 1
 petsc = 1
 meshiocubit = 1
-explicitelasticity = 1
+implicitelasticity = 1
 quadrature2d = 1
 fiatlagrange = 1
 faultcohesivedyn = 1
@@ -18,7 +18,7 @@
 [pylithapp.journal.debug]
 pylithapp = 1
 problem = 1
-explicit = 1
+implicit = 1
 
 # ----------------------------------------------------------------------
 # mesh_generator
@@ -50,8 +50,8 @@
 solver = pylith.problems.SolverNonlinear
 
 [pylithapp.timedependent.implicit.time_step]
-total_time = 10.0*s
-dt = 0.2*s
+total_time = 2.0*s
+dt = 0.05*s
 
 # ----------------------------------------------------------------------
 # materials
@@ -77,7 +77,7 @@
 [pylithapp.timedependent.bc]
 xpos = pylith.bc.Neumann
 xneg = pylith.bc.Neumann
-ypos = pylith.bc.Neumann
+ypos = pylith.bc.DirichletBC
 yneg = pylith.bc.DirichletBC
 
 # Dirichlet BC on -y
@@ -86,19 +86,16 @@
 label = yneg
 
 
-# Neumann BC on +y
+# Dirichlet BC on +y
 [pylithapp.timedependent.bc.ypos]
+bc_dof = [0, 1]
 label = ypos
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
 
-db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Neumann BC on +y
-db_initial.values = [traction-shear, traction-normal]
-db_initial.data = [-6.0*MPa, 0.0*MPa]
+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]
 
-
 # Neumann BC on +x
 [pylithapp.timedependent.bc.xpos]
 label = xpos
@@ -109,7 +106,7 @@
 db_initial = spatialdata.spatialdb.UniformDB
 db_initial.label = Neumann BC on +x
 db_initial.values = [traction-shear, traction-normal]
-db_initial.data = [6.0*MPa, 0.0*MPa]
+db_initial.data = [0.0*MPa, 0.0*MPa]
 
 
 # Neumann BC on -x
@@ -122,7 +119,7 @@
 db_initial = spatialdata.spatialdb.UniformDB
 db_initial.label = Neumann BC on -x
 db_initial.values = [traction-shear, traction-normal]
-db_initial.data = [6.0*MPa, 0.0*MPa]
+db_initial.data = [0.0*MPa, 0.0*MPa]
 
 
 # ----------------------------------------------------------------------
@@ -160,7 +157,7 @@
 ksp_gmres_restart = 100
 
 ksp_monitor = true
-ksp_view = true
+#ksp_view = true
 ksp_converged_reason = true
 
 # SNES
@@ -169,7 +166,7 @@
 snes_max_it = 500
 
 snes_monitor = true
-snes_view = true
+#snes_view = true
 snes_converged_reason = true
 
 # Friction
@@ -196,6 +193,6 @@
 
 [pylithapp.timedependent.materials.elastic.output]
 cell_info_fields = []
-cell_data_fields = []
-cell_filter = pylith.meshio.CellFilterAvgMesh
+#cell_data_fields = []
+#cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.time_format = %05.2f