[cig-commits] r7397 -
short/3D/PyLith/trunk/examples/twocells/twoquad4
willic3 at geodynamics.org
willic3 at geodynamics.org
Fri Jun 22 16:41:06 PDT 2007
Author: willic3
Date: 2007-06-22 16:41:05 -0700 (Fri, 22 Jun 2007)
New Revision: 7397
Added:
short/3D/PyLith/trunk/examples/twocells/twoquad4/README
short/3D/PyLith/trunk/examples/twocells/twoquad4/axialdisp.jpg
short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation.jpg
short/3D/PyLith/trunk/examples/twocells/twoquad4/sheardispl.jpg
Modified:
short/3D/PyLith/trunk/examples/twocells/twoquad4/axialdisp.cfg
short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation.cfg
short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation_disp.spatialdb
Log:
Two-quad example.
Everything seems to work correctly.
Added: short/3D/PyLith/trunk/examples/twocells/twoquad4/README
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twoquad4/README 2007-06-22 22:30:42 UTC (rev 7396)
+++ short/3D/PyLith/trunk/examples/twocells/twoquad4/README 2007-06-22 23:41:05 UTC (rev 7397)
@@ -0,0 +1,53 @@
+This directory contains a set of 3 elastic plane strain example
+problems. All problems share the same mesh and material properties,
+with the only differences being in the boundary conditions. The mesh
+consists of two quadrilaterals forming a rectangular mesh. The
+information describing the examples may be divided into four groups:
+
+1. Global mesh information: described in file twoquad4.mesh
+
+2. Parameter settings common to all examples: described in file
+pylithapp.cfg.
+
+3. Material properties common to all examples: described in the
+database file matprops.cfg.
+
+4. Example-specific parameters and boundary condition specifications.
+
+The examples and associated files are:
+
+1. axialdisp: Axial displacements are applied to the left and right
+edges of the mesh.
+ a. axialdisp.cfg: Parameter settings for this problem.
+ b. axialdisp.spatialdb: Boundary condition database for this
+ problem.
+
+2. sheardisp: The left and right edges of the mesh are pinned in the
+x-direction. The left edge has a negative y-displacement, while the
+right edge has a positive y-displacement.
+ a. sheardisp.cfg: Parameter settings for this problem.
+ b. sheardisp.spatialdb: Boundary condition database for this
+ problem.
+
+3. dislocation: Fault displacements are applied between the two
+elements, while the left and right edges are fixed in x and y.
+ a. dislocation.cfg: Parameter settings for this problem.
+ b. dislocation.spatialdb: Boundary condition database for this
+ problem.
+ c. dislocation_disp.spatialdb: Fault slip database for this problem.
+ d. dislocation_sliprate.spatialdb: Fault slip rate database for this
+ problem.
+ e. dislocation_sliptime.spatialdb: Fault slip time database for this
+ problem.
+
+All of the files are documented. The 'pylithapp.cfg' file is read
+automatically when running PyLith. To add to or override these
+settings, you may place an additional .cfg file on the command line.
+To run the examples, you can type:
+pylith axialdisp.cfg
+pylith sheardisp.cfg
+pylith dislocation.cfg
+
+Each of these commands will produce a corresponding VTK file, which may be viewed with visualization packages such as ParaView or Mayavi. The
+jpeg files in this directory show results for the different examples
+visualized using ParaView.
Modified: short/3D/PyLith/trunk/examples/twocells/twoquad4/axialdisp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twoquad4/axialdisp.cfg 2007-06-22 22:30:42 UTC (rev 7396)
+++ short/3D/PyLith/trunk/examples/twocells/twoquad4/axialdisp.cfg 2007-06-22 23:41:05 UTC (rev 7397)
@@ -46,7 +46,7 @@
# The name of the file containing the spatial database for the BC
# specification.
-db.iohandler.filename = axialx.spatialdb
+db.iohandler.filename = axialdisp.spatialdb
db.query_type = linear
@@ -65,7 +65,7 @@
# The name of the file containing the spatial database for the BC
# specification.
-db.iohandler.filename = axialx.spatialdb
+db.iohandler.filename = axialdisp.spatialdb
db.query_type = linear
@@ -84,7 +84,7 @@
# The name of the file containing the spatial database for the BC
# specification.
-db.iohandler.filename = axialy.spatialdb
+db.iohandler.filename = axialdisp.spatialdb
db.query_type = linear
# ----------------------------------------------------------------------
Added: short/3D/PyLith/trunk/examples/twocells/twoquad4/axialdisp.jpg
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/examples/twocells/twoquad4/axialdisp.jpg
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Modified: short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation.cfg 2007-06-22 22:30:42 UTC (rev 7396)
+++ short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation.cfg 2007-06-22 23:41:05 UTC (rev 7397)
@@ -31,7 +31,7 @@
formulation = pylith.problems.Implicit
# This is a container for a boundary condition specified at a point.
-bc = pylith.bc.BCSingle
+bc = pylith.bc.BCFourSides
# This is a container for a single fault.
interfaces = pylith.faults.SingleFault
@@ -40,25 +40,43 @@
# boundary conditions
# Provide information on the boundary conditions.
# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.bc]
+# Boundary conditions to be applied to the negative x-side of the mesh.
+[pylithapp.timedependent.bc.x_neg]
+
# We are fixing the 0 (x) and 1 (y) degrees of freedom.
fixed_dof = [0, 1]
-# Give an ID of 10 to this boundary condition.
-id = 10
+# The nodes associated with this boundary condition have the name
+# 'x_neg' in the mesh file.
+label = x_neg
+# We are assigning the label 'Dirichlet BC' to the database.
+db.label = Dirichlet BC -x edge
+
+# The name of the file containing the spatial database for the BC
+# specification.
+db.iohandler.filename = dislocation_disp.spatialdb
+db.query_type = linear
+
+
+# Boundary conditions to be applied to the positive x-side of the mesh.
+[pylithapp.timedependent.bc.x_pos]
+
+# We are fixing the 0 (x) and 1 (y) degrees of freedom.
+fixed_dof = [0, 1]
+
# The nodes associated with this boundary condition have the name
-# 'end points' in the mesh file.
-label = end points
+# 'x_pos' in the mesh file.
+label = x_pos
# We are assigning the label 'Dirichlet BC' to the database.
-db.label = Dirichlet BC
+db.label = Dirichlet BC +x edge
# The name of the file containing the spatial database for the BC
# specification.
db.iohandler.filename = dislocation_disp.spatialdb
-
+db.query_type = linear
# ----------------------------------------------------------------------
# faults
# Provide information on the fault (interface).
@@ -71,12 +89,13 @@
# Define fault properties.
[pylithapp.timedependent.interfaces.fault]
-# Give an ID of 20 to the group of nodes defining the fault.
-id = 20
-
# The nodes associated with this fault have the name 'fault' in the mesh file.
label = fault
+# We need the material database for the fault to properly condition the
+# system of equations.
+mat_db.iohandler.filename = matprops.spatialdb
+
# The quadrature for a 2D fault is 1D with a linear shape.
quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
quadrature.cell.shape = line
Added: short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation.jpg
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation.jpg
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Modified: short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation_disp.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation_disp.spatialdb 2007-06-22 22:30:42 UTC (rev 7396)
+++ short/3D/PyLith/trunk/examples/twocells/twoquad4/dislocation_disp.spatialdb 2007-06-22 23:41:05 UTC (rev 7397)
@@ -14,12 +14,12 @@
// The constraint values (displacements) have units of meters.
value-units = m m
- // The values are specified at one spatial location.
- num-locs = 1
+ // The values are specified at two spatial locations.
+ num-locs = 2
- // The dimension of the spatial distribution is 0, since the same data
- // is specified for all points in the set.
- data-dim = 0
+ // The dimension of the spatial distribution is 1, since the data
+ // is being specified at points.
+ data-dim = 1
// The spatial dimension of the database is 2.
space-dim = 2
@@ -37,7 +37,8 @@
// This is where the data is specified.
// As described in dislocation.cfg, this database will be used to specify
-// data at a point (pylith.bc.BCSingle).
+// data along an edge (pylith.bc.BCFourSides).
// The entries are:
// X-coord, Y-coord, dof-0 (x constraint), dof-1 (y constraint).
-0.0 0.0 0.0 0.0
+-2.0 0.0 0.0 0.0
+ 2.0 0.0 0.0 0.0
Added: short/3D/PyLith/trunk/examples/twocells/twoquad4/sheardispl.jpg
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/examples/twocells/twoquad4/sheardispl.jpg
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
More information about the cig-commits
mailing list