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

[surendra at geodynamics.org]

Author: surendra
Date: 2009-12-05 20:52:07 -0800 (Sat, 05 Dec 2009)
New Revision: 16077

Modified:
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-sliding.cfg
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-stick.cfg
   short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/twoquad4.mesh
Log:
Modified twoquad4-shear-sliding.cfg frition example to check with analytical solution

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-sliding.cfg	2009-12-06 00:30:29 UTC (rev 16076)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-sliding.cfg	2009-12-06 04:52:07 UTC (rev 16077)
@@ -60,9 +60,11 @@
 formulation = pylith.problems.Implicit
 formulation.solver = pylith.problems.SolverNonlinear
 
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos'
-bc = [x_neg,x_pos]
+# Set bc to an array with 2 boundary conditions: 'x_neg', 'x_pos'
+bc = [x_neg,x_pos,y_pos,y_neg]
 bc.x_pos = pylith.bc.DirichletBoundary
+bc.y_pos = pylith.bc.Neumann
+bc.y_neg = pylith.bc.Neumann
 
 # Set interfaces to an array with 1 fault: 'fault'.
 interfaces = [fault]
@@ -126,8 +128,42 @@
 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,3.0*m]
+db_initial.data = [0.0*m,1.0*m]
 
+
+
+# Boundary conditions to be applied to the positive y-side of the mesh.
+[pylithapp.timedependent.bc.y_pos]
+
+# The nodes associated with this boundary condition have the name
+# 'y_pos' in the mesh file.
+label = y_pos
+
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Neumann BC +y edge
+db_initial.values = [traction-shear,traction-normal]
+db_initial.data = [6000*MPa, 0.0*MPa]
+
+# Set cell type to quadrilateral (1-d Lagrange).
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+
+# Boundary conditions to be applied to the negative y-side of the mesh.
+[pylithapp.timedependent.bc.y_neg]
+
+# The nodes associated with this boundary condition have the name
+# 'y_neg' in the mesh file.
+label = y_neg
+
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Neumann BC -y edge
+db_initial.values = [traction-shear,traction-normal]
+db_initial.data = [-6000*MPa, 0.0*MPa]
+
+# Set cell type to quadrilateral (1-d Lagrange).
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+
 # ----------------------------------------------------------------------
 # faults
 # ----------------------------------------------------------------------

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-stick.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-stick.cfg	2009-12-06 00:30:29 UTC (rev 16076)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/shear-stick.cfg	2009-12-06 04:52:07 UTC (rev 16077)
@@ -133,7 +133,7 @@
 db_initial = spatialdata.spatialdb.UniformDB
 db_initial.label = Neumann BC +x edge
 db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [1000*MPa,0.0*MPa,]
+db_initial.data = [1000*MPa,0.0*MPa]
 
 # Set cell type to quadrilateral (1-d Lagrange).
 quadrature.cell = pylith.feassemble.FIATLagrange

Modified: short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/twoquad4.mesh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/twoquad4.mesh	2009-12-06 00:30:29 UTC (rev 16076)
+++ short/3D/PyLith/branches/pylith-friction/playpen/friction/twoquad4/twoquad4.mesh	2009-12-06 04:52:07 UTC (rev 16077)
@@ -124,9 +124,10 @@
   group = {
     name = y_pos
     type = vertices
-    count = 1
+    count = 2
     indices = {
-      3
+      1
+      5
     }
   }
 }