[cig-commits] r15237 - in short/3D/PyLith/trunk/doc/userguide: fileformats materials
brad at geodynamics.org
brad at geodynamics.org
Sun Jun 14 08:56:57 PDT 2009
Author: brad
Date: 2009-06-14 08:56:57 -0700 (Sun, 14 Jun 2009)
New Revision: 15237
Modified:
short/3D/PyLith/trunk/doc/userguide/fileformats/fileformats.lyx
short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx
Log:
Added time history format. Updated materials. Charles still needs to update the viscoelastic stuff in materials.
Modified: short/3D/PyLith/trunk/doc/userguide/fileformats/fileformats.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/fileformats/fileformats.lyx 2009-06-14 15:03:32 UTC (rev 15236)
+++ short/3D/PyLith/trunk/doc/userguide/fileformats/fileformats.lyx 2009-06-14 15:56:57 UTC (rev 15237)
@@ -1333,6 +1333,158 @@
\begin_layout Section
\begin_inset CommandInset label
LatexCommand label
+name "sec:Spatialdata:TimeHistoryIO"
+
+\end_inset
+
+Time History Database Files
+\end_layout
+
+\begin_layout Standard
+Time history database files contain a header describing the numer of points
+ in the time history and the units for the time stamps followed by a list
+ with pairs of time stamps and amplitude values.
+ The amplitude at an arbitraty point in time is computed via interpolation
+ of the values in the database.
+ This means that the time history database must span the range of time values
+ of interest.
+ The points in the time history must also be ordered in time.
+\end_layout
+
+\begin_layout LyX-Code
+// This time history database specifies temporal variation in
+\end_layout
+
+\begin_layout LyX-Code
+// amplitude.
+ In this case we prescribe a triangular slip time
+\end_layout
+
+\begin_layout LyX-Code
+// history.
+
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+// Comments can appear almost anywhere in these files and are
+\end_layout
+
+\begin_layout LyX-Code
+// delimited with two slashes (//) just like in C++.
+ All text and
+\end_layout
+
+\begin_layout LyX-Code
+// whitespace after the delimiter on a given line is ignored.
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+// The next line is the magic header for spatial database files
+\end_layout
+
+\begin_layout LyX-Code
+// in ASCII format.
+\end_layout
+
+\begin_layout LyX-Code
+#TIME HISTORY ascii
+\end_layout
+
+\begin_layout LyX-Code
+TimeHistory { // start specifying the database parameters
+\end_layout
+
+\begin_layout LyX-Code
+ num-points = 5 // number of points in time history
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+ // Specify the units used in the time stamps.
+\end_layout
+
+\begin_layout LyX-Code
+ time-unit = year
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+} // end of TimeHistory header
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+// The time history values are listed after the parameters.
+\end_layout
+
+\begin_layout LyX-Code
+// Columns time and amplitude where the amplitude values are unitless.
+\end_layout
+
+\begin_layout LyX-Code
+ 0.0 0.00
+\end_layout
+
+\begin_layout LyX-Code
+ 2.0 1.00
+\end_layout
+
+\begin_layout LyX-Code
+ 6.0 4.00
+\end_layout
+
+\begin_layout LyX-Code
+10.0 2.00
+\end_layout
+
+\begin_layout LyX-Code
+11.0 0.00
+\end_layout
+
+\begin_layout Standard
+\begin_inset Float figure
+placement H
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\begin_inset Caption
+
+\begin_layout Plain Layout
+Format of time history database files.
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Section
+\begin_inset CommandInset label
+LatexCommand label
name "sec:FileFormat:TimeStepUser"
\end_inset
Modified: short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx 2009-06-14 15:03:32 UTC (rev 15236)
+++ short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx 2009-06-14 15:56:57 UTC (rev 15237)
@@ -559,6 +559,8 @@
\begin_inset Text
\begin_layout Plain Layout
+
+\family typewriter
mu, lambda, density
\end_layout
@@ -568,6 +570,8 @@
\begin_inset Text
\begin_layout Plain Layout
+
+\family typewriter
total_strain, stress
\end_layout
@@ -588,6 +592,8 @@
\begin_inset Text
\begin_layout Plain Layout
+
+\family typewriter
mu, lambda, density, maxwell_time
\end_layout
@@ -597,9 +603,17 @@
\begin_inset Text
\begin_layout Plain Layout
-total_strain, viscous_strain, stress
+
+\family typewriter
+total_strain,
\end_layout
+\begin_layout Plain Layout
+
+\family typewriter
+stress, viscous_strain
+\end_layout
+
\end_inset
</cell>
</row>
@@ -617,6 +631,8 @@
\begin_inset Text
\begin_layout Plain Layout
+
+\family typewriter
mu, lambda, density, shear_ratio, maxwell_time
\end_layout
@@ -626,9 +642,17 @@
\begin_inset Text
\begin_layout Plain Layout
-total_strain, viscous_strain, stress
+
+\family typewriter
+total_strain,
\end_layout
+\begin_layout Plain Layout
+
+\family typewriter
+stress, viscous_strain
+\end_layout
+
\end_inset
</cell>
</row>
@@ -647,8 +671,46 @@
\begin_layout Plain Layout
+\family typewriter
+mu, lambda, density,
\end_layout
+\begin_layout Plain Layout
+
+\family typewriter
+shear_ratio_1,
+\end_layout
+
+\begin_layout Plain Layout
+
+\family typewriter
+shear_ratio_2,
+\end_layout
+
+\begin_layout Plain Layout
+
+\family typewriter
+shear_ratio_3,
+\end_layout
+
+\begin_layout Plain Layout
+
+\family typewriter
+maxwell_time_1,
+\end_layout
+
+\begin_layout Plain Layout
+
+\family typewriter
+maxwell_time_2,
+\end_layout
+
+\begin_layout Plain Layout
+
+\family typewriter
+maxwell_time_3
+\end_layout
+
\end_inset
</cell>
<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
@@ -656,8 +718,16 @@
\begin_layout Plain Layout
+\family typewriter
+total_strain, stress, viscous_strain_1,viscous_strain_2,
\end_layout
+\begin_layout Plain Layout
+
+\family typewriter
+viscous_strain_3
+\end_layout
+
\end_inset
</cell>
</row>
@@ -1398,9 +1468,8 @@
\begin_inset Tabular
-<lyxtabular version="3" rows="4" columns="3">
+<lyxtabular version="3" rows="4" columns="2">
<features>
-<column alignment="center" valignment="top" width="0.68in">
<column alignment="left" valignment="top" width="2.85in">
<column alignment="center" valignment="top" width="2.47in">
<row>
@@ -1410,17 +1479,6 @@
\begin_layout Plain Layout
\series bold
-Available?
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-
-\series bold
Model Name
\end_layout
@@ -1443,15 +1501,6 @@
\begin_inset Text
\begin_layout Plain Layout
-Yes
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
MaxwellIsotropic3D
\end_layout
@@ -1472,15 +1521,6 @@
\begin_inset Text
\begin_layout Plain Layout
-Yes
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
GenMaxwellIsotropic3D
\end_layout
@@ -1501,15 +1541,6 @@
\begin_inset Text
\begin_layout Plain Layout
-Yes
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
PowerLaw3D
\end_layout
@@ -2218,8 +2249,8 @@
may be employed for both a linear Maxwell model and a power-law Maxwell
model.
- Note that this formulation is not presently employed for linear
- viscoelastic models, but it may be employed in future versions.
+ Note that this formulation is not presently employed for linear viscoelastic
+ models, but it may be employed in future versions.
For the viscoelastic materials considered here, the viscous volumetric
strains are zero (incompressible flow), and it is convenient to separate
the general stress-strain relationship at time
@@ -2872,6 +2903,444 @@
of a power-law Maxwell viscoelastic material.
\end_layout
+\begin_layout Standard
+\noindent
+\align center
+\begin_inset Note Note
+status open
+
+\begin_layout Plain Layout
+CHARLES FIX THESE TABLES
+\end_layout
+
+\end_inset
+
+
+\begin_inset Float table
+placement H
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\noindent
+\align center
+\begin_inset Caption
+
+\begin_layout Plain Layout
+Values in spatial database used as parameters in the linear Maxwell viscoelastic
+ material constitutive model.
+
+\family typewriter
+\series bold
+FIX THIS TABLE
+\end_layout
+
+\end_inset
+
+
+\begin_inset Tabular
+<lyxtabular version="3" rows="4" columns="2">
+<features>
+<column alignment="center" valignment="middle" width="0.85in">
+<column alignment="center" valignment="middle" width="2.47in">
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\series bold
+Parameter
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\align center
+
+\series bold
+Name in Spatial Database
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $v_{p}$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+Vp
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $v_{s}$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+Vs
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $\rho$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+density
+\end_layout
+
+\end_inset
+</cell>
+</row>
+</lyxtabular>
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\noindent
+\align center
+\begin_inset Float table
+placement H
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\noindent
+\align center
+\begin_inset Caption
+
+\begin_layout Plain Layout
+Values in spatial database used as parameters in the generalized linear
+ Maxwell viscoelastic material constitutive model.
+
+\family typewriter
+\series bold
+FIX THIS TABLE
+\end_layout
+
+\end_inset
+
+
+\begin_inset Tabular
+<lyxtabular version="3" rows="4" columns="2">
+<features>
+<column alignment="center" valignment="middle" width="0.85in">
+<column alignment="center" valignment="middle" width="2.47in">
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\series bold
+Parameter
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\align center
+
+\series bold
+Name in Spatial Database
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $v_{p}$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+Vp
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $v_{s}$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+Vs
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $\rho$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+density
+\end_layout
+
+\end_inset
+</cell>
+</row>
+</lyxtabular>
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\noindent
+\align center
+\begin_inset Float table
+placement H
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\noindent
+\align center
+\begin_inset Caption
+
+\begin_layout Plain Layout
+Values in spatial database used as parameters in the nonlinear power-law
+ viscoelastic material constitutive model.
+\family typewriter
+\series bold
+FIX THIS TABLE
+\end_layout
+
+\end_inset
+
+
+\begin_inset Tabular
+<lyxtabular version="3" rows="4" columns="2">
+<features>
+<column alignment="center" valignment="middle" width="0.85in">
+<column alignment="center" valignment="middle" width="2.47in">
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\series bold
+Parameter
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\align center
+
+\series bold
+Name in Spatial Database
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $v_{p}$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+Vp
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $v_{s}$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+Vs
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\begin_inset Formula $\rho$
+\end_inset
+
+
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+density
+\end_layout
+
+\end_inset
+</cell>
+</row>
+</lyxtabular>
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
\begin_layout Section
Initial State Variables
\end_layout
@@ -2889,33 +3358,30 @@
Instead, gravity must be 'turned on' at the beginning of the simulation.
Unfortunately, this results in unrealistic amounts of deformation at the
beginning of a simulation.
- An alternative is to provide an initial state for the region under consideratio
-n.
+ An alternative is to provide an initial state variables for the region
+ under consideration.
This allows the specification of a set of state variables that is consistent
with the prior application of gravitational body forces.
In a more general sense, initial values for state variables may be used
to provide values that are consistent with any set of conditions that occurred
prior to the beginning of a simulation.
- The current release of PyLith only allows the specification of initial
- stresses.
- In future releases, however, we will allow a more general specification
- of initial state values.
+ The current release of PyLith allows the specification of initial stresses
+ and strains for elastic materials, but only initial stresses for anelastic
+ materials.
+ The next release will support initial stresses, strains, and state variables
+ for all materials.
\end_layout
\begin_layout Subsection
-Specification of Initial State Values
+Specification of Initial State Variables
\end_layout
\begin_layout Standard
Since state variables are specific to a given material, initial values for
state variables are specified as part of the material description.
By default, initial state values are not used.
- To override this behavior, the
-\family typewriter
-use_initial_state
-\family default
- flag must be set to true and a spatial database with values the initial
- stress tensor provided as in this example from the tutorial in Section
+ To override this behavior specify a spatial database for the initial stress,
+ strain, and/or state variables as in the example from the tutorial in Section
\begin_inset CommandInset ref
LatexCommand ref
@@ -2931,8 +3397,120 @@
\end_layout
\begin_layout LyX-Code
-initial_state_db.iohandler.filename = initial_state.spatialdb
+db_initial_stress = spatialdata.spatialdb.SimpleDB
\end_layout
+\begin_layout LyX-Code
+db_initial_stress.iohandler.filename = initial_stress.spatialdb
+\end_layout
+
+\begin_layout Standard
+\noindent
+\align center
+\begin_inset Float table
+placement H
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\noindent
+\align center
+\begin_inset Caption
+
+\begin_layout Plain Layout
+Values in spatial database for initial state variables for 3-D problems.
+ 2-D problems use only the relevant values.
+\end_layout
+
+\end_inset
+
+
+\begin_inset Tabular
+<lyxtabular version="3" rows="3" columns="2">
+<features>
+<column alignment="center" valignment="middle" width="0.85in">
+<column alignment="center" valignment="middle" width="2.47in">
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\series bold
+State Variable
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+\align center
+
+\series bold
+Values in Spatial Database
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+initial stress
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+stress-xx, stress-yy, stress-zz, stress-xy, stress-yz, stress-xz
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+initial strain
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Plain Layout
+
+\family typewriter
+strain-xx, strain-yy, strain-zz, strain-xy, strain-yz, strain-xz
+\end_layout
+
+\end_inset
+</cell>
+</row>
+</lyxtabular>
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
\end_body
\end_document
More information about the CIG-COMMITS
mailing list