[cig-commits] r11453 - in short/3D/PyLith/trunk/doc/userguide: boundaryconditions boundaryconditions/figs runpylith
brad at geodynamics.org
brad at geodynamics.org
Sun Mar 16 15:03:38 PDT 2008
Author: brad
Date: 2008-03-16 15:03:37 -0700 (Sun, 16 Mar 2008)
New Revision: 11453
Added:
short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.eps
short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.fig
Modified:
short/3D/PyLith/trunk/doc/userguide/boundaryconditions/boundaryconditions.lyx
short/3D/PyLith/trunk/doc/userguide/runpylith/runpylith.lyx
Log:
Worked on boundary conditions chapter.
Modified: short/3D/PyLith/trunk/doc/userguide/boundaryconditions/boundaryconditions.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/boundaryconditions/boundaryconditions.lyx 2008-03-15 18:37:54 UTC (rev 11452)
+++ short/3D/PyLith/trunk/doc/userguide/boundaryconditions/boundaryconditions.lyx 2008-03-16 22:03:37 UTC (rev 11453)
@@ -1,4 +1,4 @@
-#LyX 1.5.1 created this file. For more info see http://www.lyx.org/
+#LyX 1.5.2 created this file. For more info see http://www.lyx.org/
\lyxformat 276
\begin_document
\begin_header
@@ -116,47 +116,90 @@
\end_layout
\begin_layout Subsection
-Boundary Condition Groups
+Arrays of Boundary Condition Components
\end_layout
\begin_layout Standard
-A boundary condition group associates a type of boundary condition with
- a set of vertices via a string label.
- In previous versions of PyLith it was necessary to specify containers that
- defined the number of groups and associated label for each group.
- This was necessary because previous versions of Pyre did not support dynamic
- arrays of components, and it was necessary to predefine these arrays.
- The current version of Pythia (0.8.1.5) does support this, however, and it
- is now possible to define boundary condition groups using a
-\family typewriter
-.cfg
-\family default
- file, a
-\family typewriter
-.pml
-\family default
- file, or on the command-line.
+A dynamic array of boundary condition components associates a name (string)
+ with each boundary condition.
+ This dynamic array of boundary conditions replaces the boundary condition
+ container in PyLith 1.0.
User-defined containers are no longer necessary, and the predefined containers
are no longer available (or necessary).
- The default boundary condition for each group is the Dirichlet boundary
- condition.
- It is possible to specify as many boundary condition groups as desired.
- If no parameters are specified for a group, it is not used.
+ The default boundary condition for each component in the array is the Dirichlet
+Points object.
+ Other boundary conditions can be bound to the names itesm in the array
+ via a .cfg file, .pml file, or the command line.
+ The parameters for the boundary condition are set using the name of the
+ boundary condition.
+ An example of setting the array of boundary condition components and changing
+ the types of boundary conditions in a .cfg file:
\end_layout
+\begin_layout LyX-Code
+[pylithapp.problem]
+\end_layout
+
+\begin_layout LyX-Code
+bc = [x_neg,x_pos,y_pos,z_neg] ; Array of boundary conditions
+\end_layout
+
+\begin_layout LyX-Code
+# Default boundary condition is DirichletPoints
+\end_layout
+
+\begin_layout LyX-Code
+# Keep default value for bc.x_neg
+\end_layout
+
+\begin_layout LyX-Code
+bc.x_pos = pylith.bc.DirichletBoundary ; change BC type to DirichletBoundary
+\end_layout
+
+\begin_layout LyX-Code
+bc.y_pos = pylith.bc.AbsorbingDampers ; change BC type to AbsorbingDampers
+\end_layout
+
+\begin_layout LyX-Code
+bc.z_neg = pylith.bc.Neumann ; change BC type to Neumann (traction)
+\end_layout
+
\begin_layout Section
Dirichlet Boundary Conditions
\end_layout
\begin_layout Standard
-Dirichlet boundary conditions prescribe a fixed value of displacement on
- a subset of the vertices of the finite-element mesh.
+Dirichlet boundary conditions in PyLith prescribe the displacement of a
+ subset of the vertices of the finite-element mesh,
+\begin_inset Formula \[
+u(t)=u_{o}+(t-t_{ref})\dot{u}_{o},\]
+
+\end_inset
+
+where
+\begin_inset Formula $u_{o}$
+\end_inset
+
+is the prescribed displacement at time
+\begin_inset Formula $t_{ref}$
+\end_inset
+
+ and
+\begin_inset Formula $\dot{u}_{o}$
+\end_inset
+
+ is the rate of change of displacement (velocity).
While Dirichlet boundary conditions can be applied to any vertex, usually
they are applied to vertices on the lateral and bottom boundaries of the
domain.
- Currently, time dependent displacements are not supported.
- Allowing time variation of the prescribed displacements on the boundaries
- will likely be available in a future release.
+ Time dependent displacements with a uniform velocity are now supported.
+ The default velocity is zero (fixed displacement).
+ There are two types of Dirichlet boundary conditions, DirichletPoints and
+ DirichletBoundary.
+ Both provide identical constraints on the solution, but DirichletBoundary
+ is limited to vertices of a simply-connected surface, which allows diagnostic
+ output of the prescribed displacements.
+ DirichletPoints can be applied to a set of unconnected vertices.
\end_layout
\begin_layout Subsection
@@ -164,23 +207,22 @@
\end_layout
\begin_layout Standard
-The principal parameters for the Dirichlet boundary condition are:
+The properties and components common to both the DirichletPoints and DirichletBo
+undary boundary conditions are:
\end_layout
\begin_layout Description
-id This is an integer identifier for the boundary condition.
- It is not used in the current implementation and may be omitted.
+label Label of the group of vertices associated with the boundary condition.
\end_layout
\begin_layout Description
-name Name or label for the boundary condition.
- This associates sets of vertices with the boundary condition in addition
- to being used in error and diagnostic reports.
+db Spatial database specifying the spatial variation in the displacement
+ field (default is FixedDOFDB).
\end_layout
\begin_layout Description
-db Spatial database specifying the spatial variation in the displacement
- field.
+rate_db Spatial database specifying rate of change in the displacement field
+ (default is FixedDOFDB).
\end_layout
\begin_layout Description
@@ -188,9 +230,49 @@
is 0).
\end_layout
+\begin_layout Description
+reference_t Reference time for rate of change of values (
+\begin_inset Formula $t_{ref},$
+\end_inset
+
+ default value is 0.0 s).
+\end_layout
+
\begin_layout Standard
-An example of setting these parameters in a
+DirichletBoundary contains an additional component:
+\end_layout
+
+\begin_layout Description
+output Manager for output of displacements on boundary with specified displaceme
+nts.
+\end_layout
+
+\begin_layout Standard
+By default the output manager does not output any information.
+ The specified displacements and velocities can be output by including
+\begin_inset Quotes eld
+\end_inset
+
+displacements
+\begin_inset Quotes erd
+\end_inset
+
+ and
+\begin_inset Quotes eld
+\end_inset
+
+velocities
+\begin_inset Quotes erd
+\end_inset
+
+ in the output managers
\family typewriter
+vertex_info_fields
+\family default
+ array parameter.
+ An example of setting the Dirichlet boundary condition parameters in a
+
+\family typewriter
.cfg
\family default
file is:
@@ -201,7 +283,7 @@
\end_layout
\begin_layout LyX-Code
-bc = pylith.bc.BCTwoSides
+bc = [mybc]
\end_layout
\begin_layout LyX-Code
@@ -209,44 +291,66 @@
\end_layout
\begin_layout LyX-Code
-[pylithapp.problem.bc.pos]
+[pylithapp.problem.bc.mybc]
\end_layout
\begin_layout LyX-Code
-label = face A
+label = group A
\end_layout
\begin_layout LyX-Code
-db.iohandler.filename = disp_A.spatialdb
+fixed_dof = [2] ; fixed displacement in z direction
\end_layout
\begin_layout LyX-Code
-db.query_type = nearest
+reference_t = 3.0*yr
\end_layout
\begin_layout LyX-Code
-fixed_dof = [2]
+db = spatialdata.spatialdb.SimpleDB ; change db to SimpleDB
\end_layout
-\begin_layout Standard
-where we have specified that we want to use the
-\family typewriter
-BCTwoSides
-\family default
- boundary condition container, the group of vertices has the label
+\begin_layout LyX-Code
+db.iohandler.filename = disp_A.spatialdb
+\end_layout
+
+\begin_layout LyX-Code
+db.query_type = nearest ; change query type to nearest point algorithm
+\end_layout
+
+\begin_layout LyX-Code
+rate_db.values = [
\begin_inset Quotes eld
\end_inset
-face A,'' the spatial database query should use the nearest point algorithm,
- and motion in the z-direction is constrained.
- See Appendix
+dof-2
+\begin_inset Quotes erd
+\end_inset
+
+]
+\end_layout
+
+\begin_layout LyX-Code
+rate_db.data = [1.0e-06] ; velocity is 1.0e-06 m/s
+\end_layout
+
+\begin_layout Standard
+We have created an array with one boundary condition, mybc.
+ The group of vertices associated with the boundary condition is group A.
+ For the database associated with the displacement at the reference time
+ (3.0 yr), we change the spatial database from the default FixedDOFDB to
+ a SimpleDB.
+ We set the filename and query type for the database.
+ For the rate of change of values, we use the FixedDOFDB (which is a special
+ case of the UniformDB) and specify the velocity in the z-direction to be
+ 1.0e-06 m/s.
+ See section
\begin_inset LatexCommand ref
-reference "sec:Spatialdata:SimpleIOAscii"
+reference "sec:spatial:databases"
\end_inset
- for a complete list of the spatial database settings.
-
+ for a discussion of the different types of spatial databases available.
\end_layout
\begin_layout Subsection
@@ -334,7 +438,7 @@
\begin_layout Standard
\family typewriter
-dof-1
+dof-0
\end_layout
\end_inset
@@ -356,7 +460,7 @@
\begin_layout Standard
\family typewriter
-dof-2
+dof-1
\end_layout
\end_inset
@@ -378,7 +482,7 @@
\begin_layout Standard
\family typewriter
-dof-3
+dof-2
\end_layout
\end_inset
@@ -395,8 +499,7 @@
\begin_inset Caption
\begin_layout Standard
-Values in spatial database used as parameters in the elastic material constituti
-ve models.
+Values in spatial database used for Dirichlet boundary conditions.
\end_layout
\end_inset
@@ -415,25 +518,22 @@
\begin_layout Standard
Neumann boundary conditions are surface tractions applied over a subset
- of the mesh vertices.
- As for Dirichlet boundary conditions, each boundary condition group must
- be a simply-connected surface.
+ of the mesh.
+ As with the DirichletBoundary condition, each Neumann boundary condition
+ can only be applied to a simply-connected surface.
The surface over which the tractions are applied always has a spatial dimension
- that is one less than the dimension of the problem domain.
+ that is one less than the dimension of the finite-element mesh.
Traction values are computed at the integration points of each cell on
- the surface, using values from the spatial database.
- These are then integrated over the surface to obtain the equivalent nodal
- forces.
- An example problem using Neumann boundary conditions may be found in
-\family typewriter
-examples/twocells/twoquad4
-\family default
-.
- The Neumann example may be run by typing:
-\end_layout
+ the surface, using values from a spatial database.
+ The tractions are integrated over each cell and assembled to obtain the
+ forces applied at the vertices.
+ See section
+\begin_inset LatexCommand ref
+reference "sub:Tutorial-twoquad4-traction"
-\begin_layout LyX-Code
-pylith axialtract.cfg
+\end_inset
+
+ for a tutorial that uses Neumann boundary conditions.
\end_layout
\begin_layout Subsection
@@ -441,11 +541,12 @@
\end_layout
\begin_layout Standard
-The principal Neumann boundary condition properties are:
+The Neumann boundary condition properties and components are:
\end_layout
\begin_layout Description
-label This is a string used to identify the boundary in the mesh.
+label Name of the group of vertices defining the mesh boundary for the Neumann
+ boundary condition.
\end_layout
\begin_layout Description
@@ -458,17 +559,13 @@
upward direction differs from the default.
\end_layout
-\begin_layout Standard
-The principal Neumann boundary condition components are:
-\end_layout
-
\begin_layout Description
db The spatial database specifying the spatial variation in the applied
tractions.
\end_layout
\begin_layout Description
-output The output manager associated with diagnostic output.
+output The output manager associated with diagnostic output (traction vector).
\end_layout
\begin_layout Description
@@ -476,14 +573,29 @@
Since we are integrating over a surface that is one dimension lower than
the problem domain, this would typically be set to something like
\family typewriter
-quadrature2din3d
+Quadrature2Din3D
\family default
(for a three-dimensional problem).
\end_layout
\begin_layout Standard
-An example of setting these parameters in a
+By default the output manager does not output any information.
+ The specified tractions can be output in global coordinates by including
+
+\begin_inset Quotes eld
+\end_inset
+
+tractions
+\begin_inset Quotes erd
+\end_inset
+
+ in the output managers
\family typewriter
+cell_info_fields
+\family default
+ array parameter.
+ An example of setting these parameters in a
+\family typewriter
.cfg
\family default
file is:
@@ -494,15 +606,11 @@
\end_layout
\begin_layout LyX-Code
-formulation = pylith.problems.Implicit
-\end_layout
-
-\begin_layout LyX-Code
bc = [x_neg,x_pos,y_neg]
\end_layout
\begin_layout LyX-Code
-bc.x_pos = pylith.bc.Neumann
+bc.x_pos = pylith.bc.Neumann ; Change BC type to Neumann
\end_layout
\begin_layout LyX-Code
@@ -514,11 +622,11 @@
\end_layout
\begin_layout LyX-Code
-label = x_pos
+label = x_pos ; Name of group of vertices for +x boundary
\end_layout
\begin_layout LyX-Code
-db = spatialdata.spatialdb.SimpleDB
+db = spatialdata.spatialdb.SimpleDB ; default
\end_layout
\begin_layout LyX-Code
@@ -530,7 +638,7 @@
\end_layout
\begin_layout LyX-Code
-db.query_type = nearest
+db.query_type = nearest ; default
\end_layout
\begin_layout LyX-Code
@@ -550,62 +658,25 @@
\end_layout
\begin_layout Standard
-These settings correspond to the two-dimensional example problem in
-\family typewriter
-examples/twocells/twoquad4
-\family default
-, as described in
+These settings correspond to the example problem described in section
\begin_inset LatexCommand ref
reference "sub:Tutorial-twoquad4-traction"
\end_inset
.
- The group of vertices to which Neumann boundary conditions are to be applied
- has the label
-\begin_inset Quotes eld
-\end_inset
-
-
+ It is necessary to set the boundary condition type to
\family typewriter
-x_pos
-\family default
-'' and the spatial database query should use the nearest point algorithm.
- See Appendix
-\begin_inset LatexCommand ref
-reference "sec:Spatialdata:SimpleIOAscii"
-
-\end_inset
-
- for a complete list of the spatial database settings.
- Note that it is necessary to explicitly set the boundary condition type
- for this group to
-\family typewriter
pylith.bc.Neumann
\family default
, since the default value is
-\begin_inset Quotes eld
-\end_inset
-
-
\family typewriter
-dirichletpoints
+DirichletPoints
\family default
-
-\begin_inset Quotes erd
-\end_inset
-
.
- It is also necessary to define the quadrature information.
- Since this is a two-dimensional problem, we must integrate over edges in
- two dimensions.
- The quadrature is therefore
-\family typewriter
-Quadrature1Din2D
-\family default
-, and the cell dimension is one.
- The cell quadrature order has been set to two.
- Since constant tractions are used for this particular problem, a quadrature
+ Because the problem is two-dimensional, we must supply a quadrature object
+ to integrate over edges in two dimensions.
+ Constant tractions are used for this particular problem, so a quadrature
order of one would have been sufficient; however, for problems involving
more complex variations (e.g., a linear variation), a quadrature order of
two will provide more accurate results.
@@ -620,10 +691,10 @@
\begin_layout Standard
The spatial database files for the Neumann boundary condition specify the
applied tractions.
- The number of traction components is equal to the number of degrees of
- freedom for the problem.
- The tractions are specified in a boundary-local coordinate system.
- The traction value names for the different types of problems are:
+ The number of traction components is equal to the spatial dimension for
+ the problem.
+ The tractions are specified in a local coordinate system for the boundary.
+ The names of the components of the traction vector are:
\end_layout
\begin_layout Description
@@ -657,21 +728,12 @@
\end_layout
\begin_layout Standard
-To resolve any ambiguities in specifying the shear tractions in three-dimensiona
-l problems, users may make use of the
+Ambiguities in specifying the shear tractions in 3D problems are resolved
+ using the
\family typewriter
up_dir
\family default
- parameter, which gives the components of the vector generally thought to
- represent
-\begin_inset Quotes eld
-\end_inset
-
-up
-\begin_inset Quotes erd
-\end_inset
-
-.
+ parameter.
In the case of a horizontal surface, users will need to pick an alternative
vector, as the default
\family typewriter
@@ -703,9 +765,43 @@
waves.
If the waves arriving at the absorbing boundary are relatively small in
amplitude compared to the amplitudes of primary interest, this boundary
- condition gives reasonable results.
+ condition gives reasonable results.Neumann Boundary Condition Parameters
\end_layout
+\begin_layout Standard
+The Neumann boundary condition properties and components are:
+\end_layout
+
+\begin_layout Description
+label Name of the group of vertices defining the mesh boundary for the absorbing
+ boundary condition.
+\end_layout
+
+\begin_layout Description
+up_dir This is a 3-vector that provides a hint for the direction perpendicular
+ to the horizontal tangent direction that is not collinear with the direction
+ normal to the surface.
+ The default value is (0,0,1), which assumes that the z-axis is positive
+ upward.
+ This vector is only needed for three-dimensional problems where the positive
+ upward direction differs from the default.
+\end_layout
+
+\begin_layout Description
+db The spatial database specifying the material properties for the seismic
+ velocities.
+\end_layout
+
+\begin_layout Description
+quadrature The quadrature object to be used for numerical integration.
+ Since we are integrating over a surface that is one dimension lower than
+ the problem domain, this would typically be set to something like
+\family typewriter
+Quadrature2Din3D
+\family default
+ (for a three-dimensional problem).
+\end_layout
+
\begin_layout Subsection
Finite Element Implementation of Absorbing Boundary
\end_layout
@@ -942,9 +1038,9 @@
For fault interfaces, dislocations in 1D correspond to fault-opening (and
closing), in 2D lateral-slip and fault opening, and in 3D lateral-slip,
reverse-slip, and fault opening.
- The current release of PyLith supports kinematic (prescribed) slip.
+ The current release of PyLith only supports kinematic (prescribed) slip.
Support for fault constitutive models for frictional interface conditions
- will be supported in a future release (Winter 2007).
+ will be supported in a future release (Winter 2008).
\end_layout
\begin_layout Subsection
@@ -1074,143 +1170,60 @@
\end_layout
\begin_layout Subsection
-Fault Surface Container
+Fault Implementation
\end_layout
\begin_layout Standard
-The fault surface container associates a fault interface condition with
- a set of vertices via a string label.
- Because the number of fault surfaces depends on the problem, and Pyre does
- not support dynamic arrays of components, a generic fault surface container
- is not yet available.
- We supply a container for a single fault surface,
-\family typewriter
-pylith.faults.SingleFault
-\family default
-.
- For problems with multiple fault surfaces, the user must supply an appropriate
- fault surface container (see
-\family typewriter
-examples/templates/containers/faults.odb
-\family default
- for a template fault surface container).
- This
-\family typewriter
-.odb
-\family default
- file should be placed in the directory with the other simulation input
- files.
- The default fault surface container is an empty container, so you must
- specify a fault surface container to include fault interfaces.
- Each simply connected fault surface should be placed in a separate container.
- Additionally, fault surfaces may not share vertices.
- At fault intersections, one fault must end one cell before the intersection.
- Future releases of PyLith may relax this requirement.
+In order to create relative motion across the fault surface in the finite-elemen
+t mesh, additional degrees of freedom are added along with adjustment of
+ the topology of the mesh.
+ These additional degreees of freedom are associated with cohesive cells.
+ These zero-volume cells allow control of the relative motion between vertices
+ on the two sides of the fault.
+ PyLith automatically adds cohesive cells for each fault surface.
+ Figure
+\begin_inset LatexCommand ref
+reference "fig:fault:cohesive:quad4"
+
+\end_inset
+
+ illustrates the results of inserting a cohesive cell in a mesh consisting
+ of two quadrilateral cells.
\end_layout
\begin_layout Standard
\noindent
\align center
-\begin_inset Float table
+\begin_inset Float figure
placement H
wide false
sideways false
status open
\begin_layout Standard
-\noindent
\align center
-\begin_inset Tabular
-<lyxtabular version="3" rows="2" columns="3">
-<features>
-<column alignment="left" valignment="top" leftline="true" width="0.88in">
-<column alignment="center" valignment="top" leftline="true" width="1.3in">
-<column alignment="center" valignment="top" leftline="true" rightline="true" width="2in">
-<row topline="true" bottomline="true">
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
+\begin_inset Graphics
+ filename /Users/brad/src/cig/pylith/doc/userguide/boundaryconditions/figs/quad4cohesive.eps
-\begin_layout Standard
-
-\series bold
-Container
-\end_layout
-
\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-\begin_layout Standard
-\series bold
-Number of Fault Surfaces
\end_layout
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
\begin_layout Standard
-
-\series bold
-Container Facilities
-\end_layout
-
-\end_inset
-</cell>
-</row>
-<row topline="true" bottomline="true">
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Standard
-
-\family typewriter
-SingleFault
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Standard
-1
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Standard
-f
-\family typewriter
-ault
-\end_layout
-
-\end_inset
-</cell>
-</row>
-</lyxtabular>
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
\begin_inset Caption
\begin_layout Standard
+Example of cohesive cell inserted into a mesh of two quadrilateral cells.
+ The zero thickness cohesive cell (shown with dashed lines) controls slip
+ on the fault via the relative motion between vertices 3 and 7 and 2 and
+ 6.
\begin_inset LatexCommand label
-name "tab:fault:containers"
+name "fig:fault:cohesive:quad4"
\end_inset
-Predefined fault surface containers.
+
\end_layout
\end_inset
@@ -1223,29 +1236,7 @@
\end_layout
-\begin_layout Subsection
-Fault Implementation
-\end_layout
-
\begin_layout Standard
-In order to create relative motion across the fault surface in the finite-elemen
-t mesh, additional degrees of freedom must be added along with adjustment
- of the topology of the mesh.
- These additional degreees of freedom are associated with cohesive cells.
- These zero-volume cells allow control of the relative motion between vertices
- on the two sides of the fault.
- PyLith automatically adds cohesive cells for each fault surface.
- Figure
-\begin_inset LatexCommand ref
-reference "fig:fault:cohesive:quad4"
-
-\end_inset
-
- illustrates the results of inserting a cohesive cell in a mesh consisting
- of two quadrilateral cells.
-\end_layout
-
-\begin_layout Standard
\noindent
\align center
\begin_inset Float figure
@@ -1257,7 +1248,7 @@
\begin_layout Standard
\align center
\begin_inset Graphics
- filename figs/quad4cohesive.eps
+ filename /Users/brad/src/cig/pylith/doc/userguide/boundaryconditions/figs/quad4cohesivekin.eps
\end_inset
@@ -1268,12 +1259,18 @@
\begin_inset Caption
\begin_layout Standard
-Example of cohesive cell inserted into a mesh of two quadrilateral cells.
+Example of cohesive cell using Lagrange multipliers inserted into a mesh
+ of two quadrilateral cells.
The zero thickness cohesive cell (shown with dashed lines) controls slip
- on the fault via the relative motion between vertices 3 and 7 and 2 and
- 6.
+ on the fault.
+ The Lagrange mutipliers for a kinematic earthquake source are associated
+ with the degrees of freedom for vertices 8 and 9.
+ The constraint associated with vertex 8 controls the relative motion between
+ vertices 2 and 6, and the constraint associated with vertex 9 controls
+ the relative motion between vertices 3 and 7.
+
\begin_inset LatexCommand label
-name "fig:fault:cohesive:quad4"
+name "fig:fault:cohesive:quad4:kin"
\end_inset
@@ -1341,8 +1338,31 @@
the conditioning of the system).
\end_layout
+\begin_layout Description
+output Manager for output of diagnostic and data fields for the fault.
+\end_layout
+
\begin_layout Standard
-An example of setting these parameters in a
+By default the output manager outputs both diagnostic information (fault
+ normal direction, final slip, and slip initiation time) and the slip at
+ each time step.
+ Table
+\begin_inset LatexCommand ref
+reference "tab:fault:output"
+
+\end_inset
+
+ lists the fields available for output.
+ The fault coordinate system is shown in Figure
+\begin_inset LatexCommand ref
+reference "fig:fault:slip:motions"
+
+\end_inset
+
+.
+ The vectors in the fault coordinate system can be transformed to the global
+ coordinate system using the direction vectors in the diagnostic output.
+ An example of setting these parameters in a
\family typewriter
.cfg
\family default
@@ -1354,7 +1374,7 @@
\end_layout
\begin_layout LyX-Code
-interfaces = pylith.faults.SingleFault
+interfaces = [fault]
\end_layout
\begin_layout LyX-Code
@@ -1366,7 +1386,7 @@
\end_layout
\begin_layout LyX-Code
-fault = pylith.faults.FaultCohesiveKin
+fault = pylith.faults.FaultCohesiveKin ; default
\end_layout
\begin_layout LyX-Code
@@ -1374,19 +1394,20 @@
\end_layout
\begin_layout LyX-Code
-label = fault A
+label = fault A ; Group of vertices defining the fault surface
\end_layout
\begin_layout LyX-Code
-id = 100 // This is the default value.
+id = 100 ; Value for material identifier associated with fault's cohesive
+ cells
\end_layout
\begin_layout LyX-Code
-up_dir = [0, 0, 1] // This is the default value.
+up_dir = [0, 0, 1] ; default
\end_layout
\begin_layout LyX-Code
-normal_dir = [1, 0, 0] // This is the default value.
+normal_dir = [1, 0, 0] ; default
\end_layout
\begin_layout LyX-Code
@@ -1406,8 +1427,7 @@
\end_layout
\begin_layout Standard
-We have specified that we want to use the SingleFault container.
- The group of vertices has the label ``fault A.'' We replicate the default
+The group of vertices has the label ``fault A.'' We replicate the default
values for the fault ``up'' direction, and the desired general direction
for the fault normal.
These settings apply to a 2D fault surface embedded within a 3D mesh, so
@@ -1415,9 +1435,310 @@
reference cells.
The spatial database for elastic properties is used to determine the approximat
e shear modulus and condition the equations for faster convergence rates.
-
\end_layout
+\begin_layout Standard
+\begin_inset Float table
+placement H
+wide false
+sideways false
+status open
+
+\begin_layout Standard
+\noindent
+\align center
+\begin_inset Tabular
+<lyxtabular version="3" rows="8" columns="3">
+<features>
+<column alignment="center" valignment="top" leftline="true" width="0">
+<column alignment="center" valignment="top" leftline="true" rightline="true" width="0">
+<column alignment="center" valignment="top" rightline="true" width="0">
+<row topline="true" bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\series bold
+Field Type
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\series bold
+Field
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\series bold
+Description
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+vertex_info_fields
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+normal_dir
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Direction of fault normal in global coordinate system
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+strike_dir
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Direction of fault strike in global coordinate system
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+dip_dir
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Up-dip direction on hanging wall in global coordinate system
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+final_slip
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Vector of final slip (in fault coordinate system)
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true" bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+slip_time
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Time at which slip begins
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+vertex_data_fields
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+slip
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Slip vector at time step (in fault coordinate system)
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\family typewriter
+tractions_change
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Change in fault tractions (in fault coordinate system)
+\end_layout
+
+\end_inset
+</cell>
+</row>
+</lyxtabular>
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Caption
+
+\begin_layout Standard
+Fields available in output of fault information.
+\begin_inset LatexCommand label
+name "tab:fault:output"
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
\begin_layout Subsection
Kinematic Earthquake Source
\end_layout
@@ -1445,7 +1766,14 @@
multipliers results in a complex implementation.
However, the use of Lagrange multiplier constraints with cohesive cells
makes this very easy; we simply add the additional degrees of freedom associate
-d with the Lagrange multipliers to the cohesive cells.
+d with the Lagrange multipliers to the cohesive cells as shown in Figure
+
+\begin_inset LatexCommand ref
+reference "fig:fault:cohesive:quad4:kin"
+
+\end_inset
+
+.
As a result, the fault implementation is completely confined to the cohesive
cell.
Furthermore, the Lagrange multiplier constraints correspond to forces required
Added: short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.eps
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.eps (rev 0)
+++ short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.eps 2008-03-16 22:03:37 UTC (rev 11453)
@@ -0,0 +1,399 @@
+%!PS-Adobe-2.0 EPSF-2.0
+%%Title: quad4cohesivekin.fig
+%%Creator: fig2dev Version 3.2 Patchlevel 4
+%%CreationDate: Sun Mar 16 14:32:56 2008
+%%For: brad at eagle.local (Brad Aagaard)
+%%BoundingBox: 0 0 456 265
+%%Magnification: 1.0000
+%%EndComments
+/$F2psDict 200 dict def
+$F2psDict begin
+$F2psDict /mtrx matrix put
+/col-1 {0 setgray} bind def
+/col0 {0.000 0.000 0.000 srgb} bind def
+/col1 {0.000 0.000 1.000 srgb} bind def
+/col2 {0.000 1.000 0.000 srgb} bind def
+/col3 {0.000 1.000 1.000 srgb} bind def
+/col4 {1.000 0.000 0.000 srgb} bind def
+/col5 {1.000 0.000 1.000 srgb} bind def
+/col6 {1.000 1.000 0.000 srgb} bind def
+/col7 {1.000 1.000 1.000 srgb} bind def
+/col8 {0.000 0.000 0.560 srgb} bind def
+/col9 {0.000 0.000 0.690 srgb} bind def
+/col10 {0.000 0.000 0.820 srgb} bind def
+/col11 {0.530 0.810 1.000 srgb} bind def
+/col12 {0.000 0.560 0.000 srgb} bind def
+/col13 {0.000 0.690 0.000 srgb} bind def
+/col14 {0.000 0.820 0.000 srgb} bind def
+/col15 {0.000 0.560 0.560 srgb} bind def
+/col16 {0.000 0.690 0.690 srgb} bind def
+/col17 {0.000 0.820 0.820 srgb} bind def
+/col18 {0.560 0.000 0.000 srgb} bind def
+/col19 {0.690 0.000 0.000 srgb} bind def
+/col20 {0.820 0.000 0.000 srgb} bind def
+/col21 {0.560 0.000 0.560 srgb} bind def
+/col22 {0.690 0.000 0.690 srgb} bind def
+/col23 {0.820 0.000 0.820 srgb} bind def
+/col24 {0.500 0.190 0.000 srgb} bind def
+/col25 {0.630 0.250 0.000 srgb} bind def
+/col26 {0.750 0.380 0.000 srgb} bind def
+/col27 {1.000 0.500 0.500 srgb} bind def
+/col28 {1.000 0.630 0.630 srgb} bind def
+/col29 {1.000 0.750 0.750 srgb} bind def
+/col30 {1.000 0.880 0.880 srgb} bind def
+/col31 {1.000 0.840 0.000 srgb} bind def
+/col32 {0.004 0.004 0.004 srgb} bind def
+/col33 {0.973 0.973 0.973 srgb} bind def
+/col34 {0.184 0.207 0.281 srgb} bind def
+/col35 {0.445 0.496 0.656 srgb} bind def
+/col36 {0.977 0.250 0.250 srgb} bind def
+/col37 {0.766 0.000 0.012 srgb} bind def
+/col38 {0.977 0.730 0.410 srgb} bind def
+/col39 {0.941 0.496 0.000 srgb} bind def
+/col40 {0.977 0.977 0.449 srgb} bind def
+/col41 {0.879 0.879 0.000 srgb} bind def
+/col42 {0.375 0.777 0.055 srgb} bind def
+/col43 {0.230 0.488 0.031 srgb} bind def
+/col44 {0.199 0.730 0.977 srgb} bind def
+/col45 {0.117 0.430 0.566 srgb} bind def
+/col46 {0.785 0.547 0.977 srgb} bind def
+/col47 {0.379 0.004 0.684 srgb} bind def
+
+end
+save
+newpath 0 265 moveto 0 0 lineto 456 0 lineto 456 265 lineto closepath clip newpath
+-9.0 291.2 translate
+1 -1 scale
+
+/cp {closepath} bind def
+/ef {eofill} bind def
+/gr {grestore} bind def
+/gs {gsave} bind def
+/sa {save} bind def
+/rs {restore} bind def
+/l {lineto} bind def
+/m {moveto} bind def
+/rm {rmoveto} bind def
+/n {newpath} bind def
+/s {stroke} bind def
+/sh {show} bind def
+/slc {setlinecap} bind def
+/slj {setlinejoin} bind def
+/slw {setlinewidth} bind def
+/srgb {setrgbcolor} bind def
+/rot {rotate} bind def
+/sc {scale} bind def
+/sd {setdash} bind def
+/ff {findfont} bind def
+/sf {setfont} bind def
+/scf {scalefont} bind def
+/sw {stringwidth} bind def
+/tr {translate} bind def
+/tnt {dup dup currentrgbcolor
+ 4 -2 roll dup 1 exch sub 3 -1 roll mul add
+ 4 -2 roll dup 1 exch sub 3 -1 roll mul add
+ 4 -2 roll dup 1 exch sub 3 -1 roll mul add srgb}
+ bind def
+/shd {dup dup currentrgbcolor 4 -2 roll mul 4 -2 roll mul
+ 4 -2 roll mul srgb} bind def
+ /DrawEllipse {
+ /endangle exch def
+ /startangle exch def
+ /yrad exch def
+ /xrad exch def
+ /y exch def
+ /x exch def
+ /savematrix mtrx currentmatrix def
+ x y tr xrad yrad sc 0 0 1 startangle endangle arc
+ closepath
+ savematrix setmatrix
+ } def
+
+/$F2psBegin {$F2psDict begin /$F2psEnteredState save def} def
+/$F2psEnd {$F2psEnteredState restore end} def
+
+$F2psBegin
+10 setmiterlimit
+0 slj 0 slc
+ 0.06000 0.06000 sc
+%
+% Fig objects follow
+%
+%
+% here starts figure with depth 50
+% Polyline
+7.500 slw
+n 450 3300 m 1350 3300 l 1350 4200 l 450 4200 l
+ cp gs col45 s gr
+% Polyline
+n 4050 3300 m 4950 3300 l 4950 4200 l 4050 4200 l
+ cp gs col45 s gr
+% Polyline
+n 5850 3300 m 6750 3300 l 6750 4200 l 5850 4200 l
+ cp gs col45 s gr
+% Polyline
+ [60] 0 sd
+n 5250 3300 m 5550 3300 l 5550 4200 l 5250 4200 l
+ cp gs col45 s gr [] 0 sd
+% Polyline
+2 slj
+n 300 1125 m 300 2325 l 2700 2325 l 2700 1125 l
+ cp gs col45 s gr
+% Polyline
+n 1500 1125 m
+ 1500 2325 l gs col45 s gr
+% Polyline
+n 4200 1125 m 4200 2325 l 6600 2325 l 6600 1125 l
+ cp gs col45 s gr
+% Polyline
+n 5400 1125 m
+ 5400 2325 l gs col45 s gr
+% Polyline
+0 slj
+n 1650 3300 m 2550 3300 l 2550 4200 l 1650 4200 l
+ cp gs col45 s gr
+% here ends figure;
+%
+% here starts figure with depth 40
+% Ellipse
+7.500 slw
+n 6600 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 1650 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 1350 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 1350 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 2550 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 2550 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 1650 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 450 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 450 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 4050 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 4050 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 4950 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 4950 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5850 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5850 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 6750 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 6750 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5400 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5400 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5362 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5438 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5362 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5438 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5250 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5550 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5250 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5550 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5400 3300 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 5400 4200 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+/Helvetica-Bold ff 180.00 scf sf
+1500 600 m
+gs 1 -1 sc (Original Mesh) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+4050 2550 m
+gs 1 -1 sc (0) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+6675 2550 m
+gs 1 -1 sc (4) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+6675 1050 m
+gs 1 -1 sc (5) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+4050 1050 m
+gs 1 -1 sc (1) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+2775 1050 m
+gs 1 -1 sc (5) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+150 1050 m
+gs 1 -1 sc (1) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+150 2550 m
+gs 1 -1 sc (0) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+2775 2550 m
+gs 1 -1 sc (4) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+300 3225 m
+gs 1 -1 sc (1) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+1350 4425 m
+gs 1 -1 sc (2) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+1650 4425 m
+gs 1 -1 sc (2) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+300 4425 m
+gs 1 -1 sc (0) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+1350 3225 m
+gs 1 -1 sc (3) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+1650 3225 m
+gs 1 -1 sc (3) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+2625 3225 m
+gs 1 -1 sc (5) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+2625 4425 m
+gs 1 -1 sc (4) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+6825 3225 m
+gs 1 -1 sc (5) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+3975 4425 m
+gs 1 -1 sc (0) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+3900 3225 m
+gs 1 -1 sc (1) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+6825 4425 m
+gs 1 -1 sc (4) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+1500 2550 m
+gs 1 -1 sc (2) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+1500 1050 m
+gs 1 -1 sc (3) col0 sh gr
+/Helvetica-Bold ff 180.00 scf sf
+3375 4800 m
+gs 1 -1 sc (Exploded view of meshes) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 180.00 scf sf
+5400 600 m
+gs 1 -1 sc (Mesh with Cohesive Cell using Lagrange Multipliers) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5250 1275 m
+gs 1 -1 sc (3) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5475 1275 m
+gs 1 -1 sc (7) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5250 2250 m
+gs 1 -1 sc (2) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5475 2250 m
+gs 1 -1 sc (6) col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5400 4425 m
+gs 1 -1 sc (8) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5400 3225 m
+gs 1 -1 sc (9) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5625 3225 m
+gs 1 -1 sc (7) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5175 3225 m
+gs 1 -1 sc (3) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5925 3225 m
+gs 1 -1 sc (7) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+4875 3225 m
+gs 1 -1 sc (3) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5925 4425 m
+gs 1 -1 sc (6) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5625 4425 m
+gs 1 -1 sc (6) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5175 4425 m
+gs 1 -1 sc (2) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+4875 4425 m
+gs 1 -1 sc (2) dup sw pop 2 div neg 0 rm col0 sh gr
+/Helvetica-Bold ff 150.00 scf sf
+5400 2550 m
+gs 1 -1 sc (8) dup sw pop 2 div neg 0 rm col0 sh gr
+% Ellipse
+n 300 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+/Helvetica-Bold ff 150.00 scf sf
+5400 1050 m
+gs 1 -1 sc (9) dup sw pop 2 div neg 0 rm col0 sh gr
+% Ellipse
+n 300 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 1500 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 1500 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 2700 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 2700 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 4200 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 4200 2325 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% Ellipse
+n 6600 1125 37 37 0 360 DrawEllipse gs col39 1.00 shd ef gr gs col39 s gr
+
+% here ends figure;
+$F2psEnd
+rs
+showpage
Added: short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.fig
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.fig (rev 0)
+++ short/3D/PyLith/trunk/doc/userguide/boundaryconditions/figs/quad4cohesivekin.fig 2008-03-16 22:03:37 UTC (rev 11453)
@@ -0,0 +1,122 @@
+#FIG 3.2
+Landscape
+Center
+Inches
+Letter
+100.00
+Single
+-2
+1200 2
+0 32 #010101
+0 33 #f9f9f9
+0 34 #2f3548
+0 35 #727fa8
+0 36 #fa4040
+0 37 #c40003
+0 38 #fabb69
+0 39 #f17f00
+0 40 #fafa73
+0 41 #e1e100
+0 42 #60c70e
+0 43 #3b7d08
+0 44 #33bbfa
+0 45 #1e6e91
+0 46 #c98cfa
+0 47 #6101af
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 300 1125 37 37 300 1125 337 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 300 2325 37 37 300 2325 337 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1500 2325 37 37 1500 2325 1537 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1500 1125 37 37 1500 1125 1537 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2700 1125 37 37 2700 1125 2737 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2700 2325 37 37 2700 2325 2737 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4200 1125 37 37 4200 1125 4237 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4200 2325 37 37 4200 2325 4237 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6600 1125 37 37 6600 1125 6637 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6600 2325 37 37 6600 2325 6637 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1650 3300 37 37 1650 3300 1687 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1350 3300 37 37 1350 3300 1387 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1350 4200 37 37 1350 4200 1387 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2550 3300 37 37 2550 3300 2587 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 2550 4200 37 37 2550 4200 2587 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 1650 4200 37 37 1650 4200 1687 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 450 3300 37 37 450 3300 487 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 450 4200 37 37 450 4200 487 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4050 3300 37 37 4050 3300 4087 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4050 4200 37 37 4050 4200 4087 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4950 4200 37 37 4950 4200 4987 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 4950 3300 37 37 4950 3300 4987 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5850 3300 37 37 5850 3300 5887 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5850 4200 37 37 5850 4200 5887 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6750 4200 37 37 6750 4200 6787 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 6750 3300 37 37 6750 3300 6787 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 1125 37 37 5400 1125 5437 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 2325 37 37 5400 2325 5437 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5362 2325 37 37 5362 2325 5399 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5438 2325 37 37 5438 2325 5475 2325
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5362 1125 37 37 5362 1125 5399 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5438 1125 37 37 5438 1125 5475 1125
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5250 3300 37 37 5250 3300 5287 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5550 3300 37 37 5550 3300 5587 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5250 4200 37 37 5250 4200 5287 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5550 4200 37 37 5550 4200 5587 4200
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 3300 37 37 5400 3300 5437 3300
+1 3 0 1 39 39 40 -1 20 0.000 1 0.0000 5400 4200 37 37 5400 4200 5437 4200
+2 3 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 5
+ 300 1125 300 2325 2700 2325 2700 1125 300 1125
+2 1 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 2
+ 1500 1125 1500 2325
+2 3 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 5
+ 4200 1125 4200 2325 6600 2325 6600 1125 4200 1125
+2 1 0 1 45 7 50 -1 -1 0.000 2 0 -1 0 0 2
+ 5400 1125 5400 2325
+2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
+ 450 3300 1350 3300 1350 4200 450 4200 450 3300
+2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
+ 1650 3300 2550 3300 2550 4200 1650 4200 1650 3300
+2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
+ 4050 3300 4950 3300 4950 4200 4050 4200 4050 3300
+2 2 0 1 45 7 50 -1 -1 0.000 0 0 -1 0 0 5
+ 5850 3300 6750 3300 6750 4200 5850 4200 5850 3300
+2 2 1 1 45 7 50 -1 -1 4.000 0 0 -1 0 0 5
+ 5250 3300 5550 3300 5550 4200 5250 4200 5250 3300
+4 1 0 40 -1 18 12 0.0000 4 195 1245 1500 600 Original Mesh\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 4050 2550 0\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 6675 2550 4\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 6675 1050 5\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 4050 1050 1\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 2775 1050 5\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 150 1050 1\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 150 2550 0\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 2775 2550 4\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 300 3225 1\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 1350 4425 2\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 1650 4425 2\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 300 4425 0\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 1350 3225 3\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 1650 3225 3\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 2625 3225 5\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 2625 4425 4\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 6825 3225 5\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 3975 4425 0\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 3900 3225 1\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 6825 4425 4\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 1500 2550 2\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 1500 1050 3\001
+4 1 0 40 -1 18 12 0.0000 4 180 2160 3375 4800 Exploded view of meshes\001
+4 1 0 40 -1 18 12 0.0000 4 195 4560 5400 600 Mesh with Cohesive Cell using Lagrange Multipliers\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 5250 1275 3\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 5475 1275 7\001
+4 0 0 40 -1 18 10 0.0000 4 105 90 5250 2250 2\001
+4 0 0 40 -1 18 10 0.0000 4 105 75 5475 2250 6\001
+4 1 0 40 -1 18 10 0.0000 4 105 90 5400 4425 8\001
+4 1 0 40 -1 18 10 0.0000 4 105 75 5400 3225 9\001
+4 1 0 40 -1 18 10 0.0000 4 105 75 5625 3225 7\001
+4 1 0 40 -1 18 10 0.0000 4 105 90 5175 3225 3\001
+4 1 0 40 -1 18 10 0.0000 4 105 75 5925 3225 7\001
+4 1 0 40 -1 18 10 0.0000 4 105 90 4875 3225 3\001
+4 1 0 40 -1 18 10 0.0000 4 105 75 5925 4425 6\001
+4 1 0 40 -1 18 10 0.0000 4 105 75 5625 4425 6\001
+4 1 0 40 -1 18 10 0.0000 4 105 90 5175 4425 2\001
+4 1 0 40 -1 18 10 0.0000 4 105 90 4875 4425 2\001
+4 1 0 40 -1 18 10 0.0000 4 105 90 5400 2550 8\001
+4 1 0 40 -1 18 10 0.0000 4 105 75 5400 1050 9\001
Modified: short/3D/PyLith/trunk/doc/userguide/runpylith/runpylith.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/runpylith/runpylith.lyx 2008-03-15 18:37:54 UTC (rev 11452)
+++ short/3D/PyLith/trunk/doc/userguide/runpylith/runpylith.lyx 2008-03-16 22:03:37 UTC (rev 11453)
@@ -1,4 +1,4 @@
-#LyX 1.5.1 created this file. For more info see http://www.lyx.org/
+#LyX 1.5.2 created this file. For more info see http://www.lyx.org/
\lyxformat 276
\begin_document
\begin_header
@@ -1282,6 +1282,11 @@
\end_layout
\begin_layout Section
+\begin_inset LatexCommand label
+name "sec:spatial:databases"
+
+\end_inset
+
Databases for Boundaries, Interfaces, and Material Properties
\end_layout
More information about the cig-commits
mailing list