[cig-commits] r16095 - in short/3D/PyLith/trunk/playpen/friction: twohex8 twoquad4
surendra at geodynamics.org
surendra at geodynamics.org
Thu Dec 10 23:57:57 PST 2009
Author: surendra
Date: 2009-12-10 23:57:57 -0800 (Thu, 10 Dec 2009)
New Revision: 16095
Modified:
short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg
short/3D/PyLith/trunk/playpen/friction/twohex8/twohex8.mesh
short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg
short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg
short/3D/PyLith/trunk/playpen/friction/twoquad4/twoquad4.mesh
Log:
Modified hex8 shear-sliding friction example to verify with analytical solution
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg 2009-12-10 01:52:12 UTC (rev 16094)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg 2009-12-11 07:57:57 UTC (rev 16095)
@@ -35,9 +35,13 @@
formulation = pylith.problems.Implicit
formulation.solver = pylith.problems.SolverNonlinear
-# Set bc to an array with 3 boundary conditions: 'x_neg' and 'x_pos'.
-bc = [x_neg,x_pos]
+# Set bc to an array with 2 boundary conditions: 'x_neg' and 'x_pos'.
+bc = [x_neg,x_pos,y_pos,y_neg,z_pos,z_neg]
bc.x_pos = pylith.bc.DirichletBoundary
+bc.y_pos = pylith.bc.Neumann
+bc.y_neg = pylith.bc.Neumann
+bc.z_pos = pylith.bc.Neumann
+bc.z_neg = pylith.bc.Neumann
# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
@@ -88,7 +92,7 @@
# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
-# We are prescribing the 2 (z) and 1 (y) and fixing 0 (x) degrees of freedom.
+# We are fixing 0 (x) and prescribing 1 (y) and 2 (z)degrees of freedom.
bc_dof = [0,1,2]
# The nodes associated with this boundary condition have the name
@@ -98,8 +102,97 @@
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Dirichlet BC +x edge
db_initial.values = [displacement-x,displacement-y,displacement-z]
-db_initial.data = [0.0*m,1.0*m,0.5*m]
+#db_initial.data = [0.0*m,1.0*m,1.0*m]
+db_initial.data = [0.0*m,0.6*m,0.8*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-horiz,traction-shear-vert,traction-normal]
+#db_initial.data = [-6000.0*MPa,0.0*MPa,0.0*MPa]
+db_initial.data = [-3600.0*MPa,0.0*MPa,0.0*MPa]
+
+# Set cell type to quadrilateral (2-d Lagrange).
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+
+
+# 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-horiz,traction-shear-vert,traction-normal]
+#db_initial.data = [-6000.0*MPa,0.0*MPa,0.0*MPa]
+db_initial.data = [-3600.0*MPa,0.0*MPa,0.0*MPa]
+
+# Set cell type to quadrilateral (2-d Lagrange).
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+
+
+# Boundary conditions to be applied to the positive z-side of the mesh.
+[pylithapp.timedependent.bc.z_pos]
+
+# The nodes associated with this boundary condition have the name
+# 'z_pos' in the mesh file.
+label = z_pos
+
+up_dir = [0, 1, 0]
+
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Neumann BC +z edge
+db_initial.values = [traction-shear-horiz,traction-shear-vert,traction-normal]
+#db_initial.data = [0.0*MPa,-6000.0*MPa,0.0*MPa]
+db_initial.data = [-4800.0*MPa,0.0*MPa,0.0*MPa]
+
+# Set cell type to quadrilateral (2-d Lagrange).
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+
+
+# Boundary conditions to be applied to the negative z-side of the mesh.
+[pylithapp.timedependent.bc.z_neg]
+
+# The nodes associated with this boundary condition have the name
+# 'z_neg' in the mesh file.
+label = z_neg
+
+up_dir = [0, 1, 0]
+
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Neumann BC -z edge
+db_initial.values = [traction-shear-horiz,traction-shear-vert,traction-normal]
+#db_initial.data = [0.0*MPa,-6000.0*MPa,0.0*MPa]
+db_initial.data = [-4800.0*MPa,0.0*MPa,0.0*MPa]
+
+# Set cell type to quadrilateral (2-d Lagrange).
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+
+
+
+
+
+
+
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/twohex8.mesh
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/twohex8.mesh 2009-12-10 01:52:12 UTC (rev 16094)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/twohex8.mesh 2009-12-11 07:57:57 UTC (rev 16095)
@@ -103,10 +103,19 @@
type = vertices
count = 6
indices = {
- 0 2 4 6 8 10
+ 0 2 4 6 8 10
}
}
+ group = {
+ name = y_pos
+ type = vertices
+ count = 6
+ indices = {
+ 1 3 5 7 9 11
+ }
+ }
+
// y_neg with fault vertices removed
group = {
name = y_neg_nofault
@@ -126,6 +135,15 @@
}
}
+ group = {
+ name = z_pos
+ type = vertices
+ count = 6
+ indices = {
+ 6 7 8 9 10 11
+ }
+ }
+
// z_neg with fault vertices removed
group = {
name = z_neg_nofault
Modified: short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg 2009-12-10 01:52:12 UTC (rev 16094)
+++ short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg 2009-12-11 07:57:57 UTC (rev 16095)
@@ -61,10 +61,10 @@
formulation.solver = pylith.problems.SolverNonlinear
# Set bc to an array with 2 boundary conditions: 'x_neg', 'x_pos'
-bc = [x_neg,x_pos,y_pos,y_neg]
+bc = [x_neg,x_pos,y_pos_fault,y_neg_fault]
bc.x_pos = pylith.bc.DirichletBoundary
-bc.y_pos = pylith.bc.Neumann
-bc.y_neg = pylith.bc.Neumann
+bc.y_pos_fault = pylith.bc.Neumann
+bc.y_neg_fault = pylith.bc.Neumann
# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
@@ -133,27 +133,27 @@
# Boundary conditions to be applied to the positive y-side of the mesh.
-[pylithapp.timedependent.bc.y_pos]
+[pylithapp.timedependent.bc.y_pos_fault]
# The nodes associated with this boundary condition have the name
-# 'y_pos' in the mesh file.
-label = y_pos
+# 'y_pos_fault' in the mesh file.
+label = y_pos_fault
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]
+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]
+[pylithapp.timedependent.bc.y_neg_fault]
# The nodes associated with this boundary condition have the name
-# 'y_neg' in the mesh file.
-label = y_neg
+# 'y_neg_fault' in the mesh file.
+label = y_neg_fault
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Neumann BC -y edge
Modified: short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg 2009-12-10 01:52:12 UTC (rev 16094)
+++ short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg 2009-12-11 07:57:57 UTC (rev 16095)
@@ -61,8 +61,10 @@
formulation.solver = pylith.problems.SolverNonlinear
# 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_neg,x_pos_disp,x_pos_tract,y_pos_fault,y_neg_fault]
bc.x_pos_tract = pylith.bc.Neumann
+bc.y_pos_fault = pylith.bc.Neumann
+bc.y_neg_fault = pylith.bc.Neumann
# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
@@ -133,12 +135,49 @@
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.0*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 positive y-side of the mesh.
+[pylithapp.timedependent.bc.y_pos_fault]
+
+# The nodes associated with this boundary condition have the name
+# 'y_pos_fault' in the mesh file.
+label = y_pos_fault
+
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Neumann BC +y edge
+db_initial.values = [traction-shear,traction-normal]
+db_initial.data = [-1000.0*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_fault]
+
+# The nodes associated with this boundary condition have the name
+# 'y_neg_fault' in the mesh file.
+label = y_neg_fault
+
+db_initial = spatialdata.spatialdb.UniformDB
+db_initial.label = Neumann BC -y edge
+db_initial.values = [traction-shear,traction-normal]
+db_initial.data = [-1000.0*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/trunk/playpen/friction/twoquad4/twoquad4.mesh
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twoquad4/twoquad4.mesh 2009-12-10 01:52:12 UTC (rev 16094)
+++ short/3D/PyLith/trunk/playpen/friction/twoquad4/twoquad4.mesh 2009-12-11 07:57:57 UTC (rev 16095)
@@ -84,7 +84,7 @@
}
// This group of vertices may be used to specify boundary conditions.
- // There are 3 vertices corresponding to indices 0, 2, and 4.
+ // There are 2 vertices corresponding to indices 0 and 4.
group = {
name = y_neg
type = vertices
@@ -96,6 +96,19 @@
}
// This group of vertices may be used to specify boundary conditions.
+ // There are 3 vertices corresponding to indices 0, 2, and 4.
+ group = {
+ name = y_neg_fault
+ type = vertices
+ count = 3
+ indices = {
+ 0
+ 2
+ 4
+ }
+ }
+
+ // This group of vertices may be used to specify boundary conditions.
// There are 2 vertices corresponding to indices 0 and 1.
group = {
name = x_neg
@@ -120,7 +133,7 @@
}
// This group of vertices may be used to specify boundary conditions.
- // There is 1 vertex corresponding to index 3.
+ // There is 2 vertex corresponding to index 3 and 5.
group = {
name = y_pos
type = vertices
@@ -130,4 +143,17 @@
5
}
}
+
+ // This group of vertices may be used to specify boundary conditions.
+ // There is 3 vertex corresponding to index 1, 3 and 5.
+ group = {
+ name = y_pos_fault
+ type = vertices
+ count = 3
+ indices = {
+ 1
+ 3
+ 5
+ }
+ }
}
More information about the CIG-COMMITS
mailing list