[cig-commits] r11647 - short/3D/PyLith/trunk/doc/userguide/tutorials/twoquad4
brad at geodynamics.org
brad at geodynamics.org
Sun Mar 30 18:47:40 PDT 2008
Author: brad
Date: 2008-03-30 18:47:39 -0700 (Sun, 30 Mar 2008)
New Revision: 11647
Modified:
short/3D/PyLith/trunk/doc/userguide/tutorials/twoquad4/twoquad4.lyx
Log:
Updated twoquad4 tutorial to account for addition of fault with zero slip.
Modified: short/3D/PyLith/trunk/doc/userguide/tutorials/twoquad4/twoquad4.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/tutorials/twoquad4/twoquad4.lyx 2008-03-31 00:31:37 UTC (rev 11646)
+++ short/3D/PyLith/trunk/doc/userguide/tutorials/twoquad4/twoquad4.lyx 2008-03-31 01:47:39 UTC (rev 11647)
@@ -752,142 +752,145 @@
\end_layout
\begin_layout Subsection
-\begin_inset LatexCommand label
-name "sub:Tutorial-twoquad4-traction"
-
-\end_inset
-
-Axial Traction Example
+Kinematic Fault Slip Example
\end_layout
\begin_layout Standard
-The third example demonstrates the use of Neumann (traction) boundary conditions.
- Constant tractions are applied to the right edge of the mesh, while displacemen
-ts normal to the boundaries are held fixed along the left and bottom edges
- of the mesh.
+The next example problem is a left-lateral fault slip applied between the
+ two square cells using kinematic cohesive cells.
+ The left and right boundaries are held fixed in the x and y directions.
Parameter settings that override or augment those in
\family typewriter
pylithapp.cfg
\family default
are contained in the file
\family typewriter
-axialtract.cfg
+dislocation.cfg
\family default
.
These settings are:
\end_layout
\begin_layout Description
-pylithapp.timedependent Specifies an implicit formulation for the problem
- and specifies the array of boundary conditions.
- The boundary condition type for
-\family typewriter
-x_pos
-\family default
- is explicitly set to
-\family typewriter
-Neumann
-\family default
-, since the default boundary condition type is
-\family typewriter
-Dirichlet
-\family default
-.
+pylithapp.journal.info Turns on journaling for 1D quadrature (used for 2D
+ faults) and for cohesive kinematic faults.
\end_layout
\begin_layout Description
+pylithapp.timedependent Specifies an implicit formulation for the problem,
+ the array of boundary conditions, and the array of interfaces.
+\end_layout
+
+\begin_layout Description
pylithapp.timedependent.bc.x_neg Specifies the boundary conditions for the
left side of the mesh, defining which degrees of freedom are being constrained
- (x) and giving the label (defined in
+ (x and y), giving the label (
\family typewriter
+x_neg
+\family default
+, defined in
+\family typewriter
twoquad4.mesh
\family default
-) defining the points desired.
- In this case, rather than specifying a spatial database file with values
- for the Dirichlet boundary conditions, we use the default spatial database
- (FixedDOFDB) for the Dirichlet boundary condition, which sets the displacements
- to zero for all time.
+) defining the points desired, and assigning a label to the boundary condition
+ set.
+ Instead of specifying a spatial database file for the values of the Dirichlet
+ boundary condition, we use the default spatial database (FixedDOFDB) for
+ the Dirichlet boundary condition, which sets the displacements to zero
+ for all time.
\end_layout
\begin_layout Description
-pylithapp.timedependent.bc.x_pos Specifies the Neumann boundary conditions
- for the right side of the mesh, giving the label (defined in
+pylithapp.timedependent.bc.x_pos Specifies the boundary conditions for the
+ right side of the mesh, defining which degrees of freedom are being constrained
+ (x and y), giving the label (
\family typewriter
-twoquad4.mesh
+x_neg
\family default
-) defining the points desired, assigning a label to the boundary condition
- set, and giving the name of the spatial database with the traction vectors
- for the Neumann boundary condition (
+, defined in
\family typewriter
-axialtract.spatialdb
+twoquad4.mesh
\family default
-).
+) defining the points desired, and assigning a label to the boundary condition
+ set.
+ We use the FixedDOFDB for this boundary condition as well, which sets the
+ displacements to zero for all time.
\end_layout
\begin_layout Description
-pylithapp.timedependent.bc.y_neg Specifies the boundary conditions for the
- bottom two corners of the mesh, defining which degrees of freedom are being
- constrained (y) and giving the label (defined in
+pylithapp.timedependent.interfaces Gives the label (defined in
\family typewriter
twoquad4.mesh
\family default
-) defining the points desired.
- In this case, we again use the FixedDOFDB, which sets the displacements
- to zero for all time.
+) defining the points on the fault, provides quadrature information, and
+ then gives database names for material properties (needed for conditioning),
+ fault slip, peak fault slip rate, and fault slip time.
\end_layout
\begin_layout Description
pylithapp.problem.formulation.output.output.writer Gives the base filename for
VTK output (
\family typewriter
-axialtract.vtk
+dislocation.vtk
\family default
).
\end_layout
\begin_layout Description
-pylithapp.timedependent.materials.material.output Defines the filter to be used
- when writing cell state variables (average the quadrature points for the
- cell) and gives the base filename for state variable output files (
+pylithapp.timedependent.interfaces.fault.output.writer Gives the base filename
+ for cohesive cell output files (
\family typewriter
-axialtract-statevars.vtk
+dislocation-fault.vtk
\family default
).
\end_layout
\begin_layout Description
-pylithapp.timedependent.bc.x_pos.output Gives the field to be output for the
-
+pylithapp.timedependent.materials.material.output Defines the filter to be used
+ when writing cell state variables (average the quadrature points for the
+ cell) and gives the base filename for state variable output files (
\family typewriter
-x_pos
+dislocation-statevars.vtk
\family default
- boundary (
-\family typewriter
-tractions
-\family default
-), and gives the base filename for
-\family typewriter
-x_pos
-\family default
- boundary output (
-\family typewriter
-axialtract-tractions.vtk
-\family default
).
\end_layout
\begin_layout Standard
-The traction vectors for the Neumann boundary conditions are given in the
+The fault example requires three additional database files that were not
+ needed for the simple displacement examples.
+ The first file (
+\family typewriter
+dislocation_slip.spatialdb
+\family default
+) specifies 0.01 m of left-lateral fault slip for the entire fault.
+ The data dimension is zero since the same data are applied to all points
+ in the set.
+ It is also necessary to specify the peak slip rate, which is done in the
file
\family typewriter
-axialtract.spatialdb
+dislocation_sliprate.spatialdb
\family default
-, as specified in
+.
+ The peak slip rate of 1.0e6 m/s creates a step function for the slip time
+ history.
+ Finally, we must provide the time at which slip begins.
+ The elastic solution begins at
+\begin_inset Formula $t=-dt$
+\end_inset
+
+, so this is the value given (
+\begin_inset Formula $t=-1\sec$
+\end_inset
+
+) in
\family typewriter
-axialtract.cfg
+dislocation_sliptime.spatialdb
\family default
.
- The files containing common information (
+\end_layout
+
+\begin_layout Standard
+The files containing common information (
\family typewriter
twoquad4.mesh
\family default
@@ -901,25 +904,35 @@
\family default
) along with the problem-specific files (
\family typewriter
-axialtract.cfg
+\size small
+dislocation.cfg
\family default
,
\family typewriter
-axialtract.spatialdb
+dislocation_slip.spatialdb
\family default
+,
+\family typewriter
+dislocation_sliprate.spatialdb
+\family default
+,
+\family typewriter
+dislocation_sliptime.spatialdb
+\family default
+\size default
) provide a complete description of the problem, and we can then run this
example by typing
\end_layout
\begin_layout LyX-Code
-pylith axialtract.cfg
+pylith dislocation.cfg
\end_layout
\begin_layout Standard
-Once the problem has run, four files will be produced.
+Once the problem has run, five files are produced.
The first file is named
\family typewriter
-axialtract_t0000000.vtk
+dislocation_t0000000.vtk
\family default
.
The
@@ -932,20 +945,27 @@
mesh vertices.
The second file is named
\family typewriter
-axialtract-statevars_t0000000.vtk
+dislocation-statevars_t0000000.vtk
\family default
.
This file contains the state variables for each cell.
The default fields are the total strain and stress fields.
As specified in
\family typewriter
-axialtract.cfg
+dislocation.cfg
\family default
, these values are averaged over each cell.
- The third file (
+ The third file is named
\family typewriter
-axialtract-statevars_info.vtk
+dislocation-fault_t0000000.vtk
\family default
+.
+ This file gives the specified fault slip for each vertex on the fault,
+ along with the computed traction change for the cohesive cell.
+ The fourth file (
+\family typewriter
+dislocation-statevars_info.vtk
+\family default
) gives the material properties used for the problem.
Since we have not specified which properties to write, the default properties
(
@@ -963,17 +983,10 @@
) are written.
The final file (
\family typewriter
-axialtract-tractions_info.vtk
+dislocation-fault_info.vtk
\family default
-) gives the
-\family typewriter
-x
-\family default
- and
-\family typewriter
-y
-\family default
- components of traction applied at each integration point.
+) provides information such as the normal direction, final slip, and slip
+ time for each vertex on the fault.
All of the
\family typewriter
.vtk
@@ -981,20 +994,12 @@
files may be used with a number of visualization packages.
If the problem ran correctly, you should be able to generate a figure such
as
-\begin_inset LatexCommand vref
-reference "fig:twoquad4-axialtract"
+\begin_inset LatexCommand ref
+reference "fig:twoquad4-disloc"
\end_inset
, which was generated using ParaView.
- The results may be compared against the analytical solution derived in
- Section
-\begin_inset LatexCommand ref
-reference "sub:Analytical-Constant-Traction"
-
-\end_inset
-
-.
\end_layout
\begin_layout Standard
@@ -1006,7 +1011,7 @@
\begin_layout Standard
\align center
\begin_inset Graphics
- filename figs/axialtract.jpg
+ filename figs/dislocation.jpg
lyxscale 50
scale 45
@@ -1019,9 +1024,9 @@
\begin_inset Caption
\begin_layout Standard
-Color contours and vectors of displacement for the axial traction example.
+Color contours and vectors of displacement for the kinematic fault example.
\begin_inset LatexCommand label
-name "fig:twoquad4-axialtract"
+name "fig:twoquad4-disloc"
\end_inset
@@ -1033,155 +1038,148 @@
\end_layout
-\begin_layout Standard
-
-\end_layout
-
\end_inset
\end_layout
\begin_layout Subsection
-Kinematic Fault Slip Example
+\begin_inset LatexCommand label
+name "sub:Tutorial-twoquad4-traction"
+
+\end_inset
+
+Axial Traction Example
\end_layout
\begin_layout Standard
-The next example problem is a left-lateral fault slip applied between the
- two square cells using kinematic cohesive cells.
- The left and right boundaries are held fixed in the x and y directions.
+The third example demonstrates the use of Neumann (traction) boundary conditions.
+ Constant tractions are applied to the right edge of the mesh, while displacemen
+ts normal to the boundaries are held fixed along the left and bottom edges
+ of the mesh.
Parameter settings that override or augment those in
\family typewriter
pylithapp.cfg
\family default
are contained in the file
\family typewriter
-dislocation.cfg
+axialtract.cfg
\family default
.
These settings are:
\end_layout
\begin_layout Description
-pylithapp.journal.info Turns on journaling for 1D quadrature (used for 2D
- faults) and for cohesive kinematic faults.
+pylithapp.timedependent Specifies an implicit formulation for the problem
+ and specifies the array of boundary conditions.
+ The boundary condition type for
+\family typewriter
+x_pos
+\family default
+ is explicitly set to
+\family typewriter
+Neumann
+\family default
+, since the default boundary condition type is
+\family typewriter
+Dirichlet
+\family default
+.
\end_layout
\begin_layout Description
-pylithapp.timedependent Specifies an implicit formulation for the problem,
- the array of boundary conditions, and the array of interfaces.
-\end_layout
-
-\begin_layout Description
pylithapp.timedependent.bc.x_neg Specifies the boundary conditions for the
left side of the mesh, defining which degrees of freedom are being constrained
- (x and y), giving the label (
+ (x) and giving the label (defined in
\family typewriter
-x_neg
-\family default
-, defined in
-\family typewriter
twoquad4.mesh
\family default
-) defining the points desired, and assigning a label to the boundary condition
- set.
- Instead of specifying a spatial database file for the values of the Dirichlet
- boundary condition, we use the default spatial database (FixedDOFDB) for
- the Dirichlet boundary condition, which sets the displacements to zero
- for all time.
+) defining the points desired.
+ In this case, rather than specifying a spatial database file with values
+ for the Dirichlet boundary conditions, we use the default spatial database
+ (FixedDOFDB) for the Dirichlet boundary condition, which sets the displacements
+ to zero for all time.
\end_layout
\begin_layout Description
-pylithapp.timedependent.bc.x_pos Specifies the boundary conditions for the
- right side of the mesh, defining which degrees of freedom are being constrained
- (x and y), giving the label (
+pylithapp.timedependent.bc.x_pos Specifies the Neumann boundary conditions
+ for the right side of the mesh, giving the label (defined in
\family typewriter
-x_neg
+twoquad4.mesh
\family default
-, defined in
+) defining the points desired, assigning a label to the boundary condition
+ set, and giving the name of the spatial database with the traction vectors
+ for the Neumann boundary condition (
\family typewriter
-twoquad4.mesh
+axialtract.spatialdb
\family default
-) defining the points desired, and assigning a label to the boundary condition
- set.
- We use the FixedDOFDB for this boundary condition as well, which sets the
- displacements to zero for all time.
+).
\end_layout
\begin_layout Description
-pylithapp.timedependent.interfaces Gives the label (defined in
+pylithapp.timedependent.bc.y_neg Specifies the boundary conditions for the
+ bottom two corners of the mesh, defining which degrees of freedom are being
+ constrained (y) and giving the label (defined in
\family typewriter
twoquad4.mesh
\family default
-) defining the points on the fault, provides quadrature information, and
- then gives database names for material properties (needed for conditioning),
- fault slip, peak fault slip rate, and fault slip time.
+) defining the points desired.
+ In this case, we again use the FixedDOFDB, which sets the displacements
+ to zero for all time.
\end_layout
\begin_layout Description
pylithapp.problem.formulation.output.output.writer Gives the base filename for
VTK output (
\family typewriter
-dislocation.vtk
+axialtract.vtk
\family default
).
\end_layout
\begin_layout Description
-pylithapp.timedependent.interfaces.fault.output.writer Gives the base filename
- for cohesive cell output files (
+pylithapp.timedependent.materials.material.output Defines the filter to be used
+ when writing cell state variables (average the quadrature points for the
+ cell) and gives the base filename for state variable output files (
\family typewriter
-dislocation-fault.vtk
+axialtract-statevars.vtk
\family default
).
\end_layout
\begin_layout Description
-pylithapp.timedependent.materials.material.output Defines the filter to be used
- when writing cell state variables (average the quadrature points for the
- cell) and gives the base filename for state variable output files (
+pylithapp.timedependent.bc.x_pos.output Gives the field to be output for the
+
\family typewriter
-dislocation-statevars.vtk
+x_pos
\family default
+ boundary (
+\family typewriter
+tractions
+\family default
+), and gives the base filename for
+\family typewriter
+x_pos
+\family default
+ boundary output (
+\family typewriter
+axialtract-tractions.vtk
+\family default
).
\end_layout
\begin_layout Standard
-The fault example requires three additional database files that were not
- needed for the simple displacement examples.
- The first file (
-\family typewriter
-dislocation_slip.spatialdb
-\family default
-) specifies 0.01 m of left-lateral fault slip for the entire fault.
- The data dimension is zero since the same data are applied to all points
- in the set.
- It is also necessary to specify the peak slip rate, which is done in the
+The traction vectors for the Neumann boundary conditions are given in the
file
\family typewriter
-dislocation_sliprate.spatialdb
+axialtract.spatialdb
\family default
-.
- The peak slip rate of 1.0e6 m/s creates a step function for the slip time
- history.
- Finally, we must provide the time at which slip begins.
- The elastic solution begins at
-\begin_inset Formula $t=-dt$
-\end_inset
-
-, so this is the value given (
-\begin_inset Formula $t=-1\sec$
-\end_inset
-
-) in
+, as specified in
\family typewriter
-dislocation_sliptime.spatialdb
+axialtract.cfg
\family default
.
-\end_layout
-
-\begin_layout Standard
-The files containing common information (
+ The files containing common information (
\family typewriter
twoquad4.mesh
\family default
@@ -1195,35 +1193,25 @@
\family default
) along with the problem-specific files (
\family typewriter
-\size small
-dislocation.cfg
+axialtract.cfg
\family default
,
\family typewriter
-dislocation_slip.spatialdb
+axialtract.spatialdb
\family default
-,
-\family typewriter
-dislocation_sliprate.spatialdb
-\family default
-,
-\family typewriter
-dislocation_sliptime.spatialdb
-\family default
-\size default
) provide a complete description of the problem, and we can then run this
example by typing
\end_layout
\begin_layout LyX-Code
-pylith dislocation.cfg
+pylith axialtract.cfg
\end_layout
\begin_layout Standard
-Once the problem has run, five files are produced.
+Once the problem has run, six files will be produced.
The first file is named
\family typewriter
-dislocation_t0000000.vtk
+axialtract_t0000000.vtk
\family default
.
The
@@ -1236,27 +1224,20 @@
mesh vertices.
The second file is named
\family typewriter
-dislocation-statevars_t0000000.vtk
+axialtract-statevars_t0000000.vtk
\family default
.
This file contains the state variables for each cell.
The default fields are the total strain and stress fields.
As specified in
\family typewriter
-dislocation.cfg
+axialtract.cfg
\family default
, these values are averaged over each cell.
- The third file is named
+ The third file (
\family typewriter
-dislocation-fault_t0000000.vtk
+axialtract-statevars_info.vtk
\family default
-.
- This file gives the specified fault slip for each vertex on the fault,
- along with the computed traction change for the cohesive cell.
- The fourth file (
-\family typewriter
-dislocation-statevars_info.vtk
-\family default
) gives the material properties used for the problem.
Since we have not specified which properties to write, the default properties
(
@@ -1272,12 +1253,34 @@
density
\family default
) are written.
- The final file (
+ The fourth file (
\family typewriter
-dislocation-fault_info.vtk
+axialtract-tractions_info.vtk
\family default
-) provides information such as the normal direction, final slip, and slip
- time for each vertex on the fault.
+) gives the
+\family typewriter
+x
+\family default
+ and
+\family typewriter
+y
+\family default
+ components of traction applied at each integration point.
+ The file
+\family typewriter
+axialtract-fault_info.vtk
+\family default
+ gives the fault normal direction,
+ slip (zero) and time of slip (-1 second) for the fault.
+ The final file
+\family typewriter
+axialtract-fault_t0000000.vtk
+\family default
+ gives the slip vector
+ and change in tractions on the fault surface in the fault coordinate
+ system (along-strike and opening).
+ The second component of the change in traction vector corresponds to
+ the applied traction, as expected.
All of the
\family typewriter
.vtk
@@ -1285,12 +1288,20 @@
files may be used with a number of visualization packages.
If the problem ran correctly, you should be able to generate a figure such
as
+\begin_inset LatexCommand vref
+reference "fig:twoquad4-axialtract"
+
+\end_inset
+
+, which was generated using ParaView.
+ The results may be compared against the analytical solution derived in
+ Section
\begin_inset LatexCommand ref
-reference "fig:twoquad4-disloc"
+reference "sub:Analytical-Constant-Traction"
\end_inset
-, which was generated using ParaView.
+.
\end_layout
\begin_layout Standard
@@ -1302,7 +1313,7 @@
\begin_layout Standard
\align center
\begin_inset Graphics
- filename figs/dislocation.jpg
+ filename figs/axialtract.jpg
lyxscale 50
scale 45
@@ -1315,9 +1326,9 @@
\begin_inset Caption
\begin_layout Standard
-Color contours and vectors of displacement for the kinematic fault example.
+Color contours and vectors of displacement for the axial traction example.
\begin_inset LatexCommand label
-name "fig:twoquad4-disloc"
+name "fig:twoquad4-axialtract"
\end_inset
@@ -1329,6 +1340,10 @@
\end_layout
+\begin_layout Standard
+
+\end_layout
+
\end_inset
More information about the cig-commits
mailing list