[cig-commits] r15806 - short/3D/PyLith/branches/pylith-friction/playpen/friction

[brad at geodynamics.org]

Author: brad
Date: 2009-10-13 13:04:01 -0700 (Tue, 13 Oct 2009)
New Revision: 15806

Modified:
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.cfg
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg
Log:
Fixed parameters for shear example. Use shear traction on +x face.

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.cfg	2009-10-13 20:02:56 UTC (rev 15805)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-axial.cfg	2009-10-13 20:04:01 UTC (rev 15806)
@@ -21,7 +21,7 @@
 timedependent = 1
 implicit = 1
 petsc = 1
-solverlinear = 1
+solvernonlinear = 1
 meshioascii = 1
 homogeneous = 1
 elasticityimplicit = 1

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg	2009-10-13 20:02:56 UTC (rev 15805)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4-shear.cfg	2009-10-13 20:04:01 UTC (rev 15806)
@@ -21,7 +21,7 @@
 timedependent = 1
 implicit = 1
 petsc = 1
-solverlinear = 1
+solvernonlinear = 1
 meshioascii = 1
 homogeneous = 1
 elasticityimplicit = 1
@@ -58,9 +58,11 @@
 dimension = 2
 normalizer.length_scale = 1.0*m
 formulation = pylith.problems.Implicit
+formulation.solver = pylith.problems.SolverNonlinear
 
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'y_neg' and 'x_pos'.
-bc = [x_neg,y_neg,x_pos]
+# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos_disp', 'x_pos_tract'
+bc = [x_neg,x_pos_disp,x_pos_tract]
+bc.x_pos_tract = pylith.bc.Neumann
 
 # Set interfaces to an array with 1 fault: 'fault'.
 interfaces = [fault]
@@ -105,37 +107,38 @@
 [pylithapp.timedependent.bc.x_neg]
 
 # We are fixing the 0 (x) degree of freedom.
-bc_dof = [0]
+bc_dof = [0,1]
 
 # The nodes associated with this boundary condition have the name
 # 'x_neg' in the mesh file.
 label = x_neg
 
-# Boundary conditions to be applied to the negative y-side of the mesh.
-[pylithapp.timedependent.bc.y_neg]
+# Boundary conditions to be applied to the posiditve x-side of the mesh.
+[pylithapp.timedependent.bc.x_pos_disp]
 
-# We are fixing the 1 (y) degree of freedom.
-bc_dof = [1]
+# We are fixing the 0 (x) degree of freedom.
+bc_dof = [0]
 
 # The nodes associated with this boundary condition have the name
 # 'x_pos' in the mesh file.
-label = y_neg
+label = x_pos
 
 # Boundary conditions to be applied to the positive x-side of the mesh.
-[pylithapp.timedependent.bc.x_pos]
+[pylithapp.timedependent.bc.x_pos_tract]
 
-# We are fixing the 0 (x) degrees of freedom.
-bc_dof = [0]
-
 # The nodes associated with this boundary condition have the name
 # 'x_pos' in the mesh file.
 label = x_pos
 
 db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Dirichlet BC +x edge
-db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*m,-0.001*m]
+db_initial.label = Neumann BC +x edge
+db_initial.values = [traction-shear,traction-normal]
+db_initial.data = [1000*MPa,0.0*MPa,]
 
+# Set cell type to quadrilateral (1-d Lagrange).
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+
 # ----------------------------------------------------------------------
 # faults
 # ----------------------------------------------------------------------
@@ -184,13 +187,18 @@
 sub_pc_factor_shift_positive_definite = 0
 sub_pc_factor_shift_nonzero = 
 
-#ksp_monitor = true
-#ksp_view = true
+ksp_monitor = true
+ksp_view = true
 #log_summary = true
 ksp_max_it = 100
 ksp_gmres_restart = 50
 #start_in_debugger = true
 
+snes_monitor = true
+snes_view = true
+ksp_converged_reason = true
+snes_converged_reason = true
+
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------