[cig-commits] r6338 - in
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet:
. vistest
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed Mar 21 14:38:27 PDT 2007
Author: willic3
Date: 2007-03-21 14:38:26 -0700 (Wed, 21 Mar 2007)
New Revision: 6338
Added:
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/README
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/genmax.cfg
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/linmaxESF.cfg
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermax.cfg
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermaxZT.cfg
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d-debug.cfg
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d.cfg
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.bc
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.connect
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.coord
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.fuldat
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.hist
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.prop
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.statevar
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.time
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestGenmax.prop
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestLinmaxESF.prop
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestNL.time
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlaw.prop
short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlawZT.prop
Log:
Directory containing simple tests for all of the viscoelastic models
presently available. The default model is a linear Maxwell model
based on the Zienkiewiecz & Taylor formulation. Other models are
obtained simply by putting the appropriate *.cfg file on the
command-line.
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/README
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/README 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/README 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,10 @@
+This is meant to be a very simple test of all the viscoelastic models
+presently available. The default model is a linear Maxwell viscoelastic
+model, which is obtained by simply typing 'pylith3dapp.py' in this
+directory. Additional models can be used by typing:
+'pylith3dapp.py xxx.cfg'
+where xxx is the name of the desired model (just look at the names of
+the xxx.cfg files in this directory).
+
+Notice that the nonlinear models (power-law) make use of a different
+time stepping description file designed for nonlinear material models.
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/genmax.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/genmax.cfg 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/genmax.cfg 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,8 @@
+# Configuration file for the generalized Maxwell model. To use,
+# simply name this file on the command line:
+#
+# pylith3dapp.py genmax.cfg
+
+[pylith3d.scanner]
+materialPropertiesInputFile = vistestGenmax.prop
+outputFileRoot = vistestGenmax
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/linmaxESF.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/linmaxESF.cfg 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/linmaxESF.cfg 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,8 @@
+# Configuration file for the linear Maxwell model with an ESF
+# formulation. To use, simply name this file on the command line:
+#
+# pylith3dapp.py linmaxESF.cfg
+
+[pylith3d.scanner]
+materialPropertiesInputFile = vistestLinmaxESF.prop
+outputFileRoot = vistestLinmaxESF
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermax.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermax.cfg 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermax.cfg 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,9 @@
+# Configuration file for the power-law Maxwell model using an ESF
+# formulation. To use, simply name this file on the command line:
+#
+# pylith3dapp.py powermax.cfg
+
+[pylith3d.scanner]
+materialPropertiesInputFile = vistestPowerlaw.prop
+outputFileRoot = vistestPowerlaw
+timeStepInputFile = vistestNL.time
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermaxZT.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermaxZT.cfg 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/powermaxZT.cfg 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,10 @@
+# Configuration file for the power-law Maxwell model using a
+# formulation based on Zienkiewicz & Taylor. To use, simply name this
+# file on the command line:
+#
+# pylith3dapp.py powermaxZT.cfg
+
+[pylith3d.scanner]
+materialPropertiesInputFile = vistestPowerlawZT.prop
+outputFileRoot = vistestPowerlawZT
+timeStepInputFile = vistestNL.time
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d-debug.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d-debug.cfg 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d-debug.cfg 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,14 @@
+[pylith3d]
+# scanner (pylith3d) options
+scanner.fileRoot = vistest
+scanner.asciiOutput = full
+scanner.ucdOutput = ascii
+# PETSc options
+log_summary = true
+pc_type = bjacobi
+sub_pc_type = ilu
+ksp_monitor = true
+ksp_view = true
+ksp_rtol = 1.0e-9
+start_in_debugger = true
+debugger_pause = 100
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d.cfg 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/pylith3d.cfg 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,113 @@
+[pylith3d]
+# PETSc options
+log_summary = true
+pc_type = bjacobi
+sub_pc_type = ilu
+ksp_monitor = true
+ksp_view = true
+ksp_rtol = 1.0e-9
+# start_in_debugger = true
+# debugger_pause = 100
+
+[pylith3d.scanner]
+fileRoot = vistest
+asciiOutput = full
+ucdOutput = ascii
+
+# Example of a keyword=value file to be used with TECTON.
+# In this example, the default values are listed, but commented out.
+# To change a default value, uncomment the appropriate entry and
+# edit the value.
+# With the present implementation, there should not be any spaces
+# before the keyword.
+#
+# Non-default parameters to be used for SCEC BM 2.
+#
+# Scaling factors applied to Winkler forces. These factors may be
+# used as a quick and easy way of changing the density or gravitational
+# acceleration when Winkler forces are used to simulate gravity.
+#
+#winklerScaleX = 1.0
+#winklerScaleY = 1.0
+#winklerScaleZ = 1.0
+#
+#
+# Parameters controlling stress integration and numerical computation
+# of the tangent material matrix. These default values should be
+# reasonable for most cases.
+#
+#stressTolerance = 1.0e-12*Pa
+#minimumStrainPerturbation = 1.0e-7
+#initialStrainPerturbation = 1.0e-1
+#
+#
+# Parameters controlling the solution of the linear problem at each
+# iteration. This is now all controlled by PETSc command-line
+# arguments. The only option now specified in this file is whether to
+# use the solution from the previous time step as the starting guess
+# for the current time step.
+#
+# usePreviousDisplacementFlag = 0
+#
+#
+# Quadrature order for the problem. The options are:
+# Full: Quadrature order that should give the exact
+# element matrices when the elements are
+# geometrically undistorted.
+# Reduced: Quadrature order that is one order less than
+# full quadrature. This option should be used
+# with caution.
+# Selective: Uses Hughes' b-bar formulation to perform
+# reduced quadrature on the dilatational parts of
+# the strain-displacement matrix. This can be
+# useful in nearly-incompressible problems.
+#
+#quadratureOrder = Full
+#
+#
+# Gravitational acceleration in each direction.
+#
+#gravityX = 0.0*m/(s*s)
+#gravityY = 0.0*m/(s*s)
+#gravityZ = 0.0*m/(s*s)
+#
+#
+# Factors controlling computation of prestresses. When gravity is being
+# used, an automatic computation option may be used, with the option of
+# using alternative values for Poisson's ratio and Young's modulus.
+# Each prestress component may also be scaled. This option is only
+# useful if the prestresses are read from a file (and not automatically
+# computed).
+#
+#prestressAutoCompute = False
+#prestressAutoChangeElasticProps = False
+#prestressAutoComputePoisson = 0.49
+#prestressAutoComputeYoungs = 1.0e30*Pa
+#
+#prestressScaleXx = 1.0
+#prestressScaleYy = 1.0
+#prestressScaleZz = 1.0
+#prestressScaleXy = 1.0
+#prestressScaleXz = 1.0
+#prestressScaleYz = 1.0
+#
+#
+# Scaling factors applied to differential Winkler forces. Differential
+# Winkler forces are those applied across a slippery node interface, and
+# are generally used to keep the fault locked at certain times. These
+# factors control the magnitudes and provide a simple way of scaling the
+# forces so the fault remains sufficiently 'locked'.
+#
+#winklerSlipScaleX = 1.0
+#winklerSlipScaleY = 1.0
+#winklerSlipScaleZ = 1.0
+#
+#
+# Unit numbers used by f77. These defaults should work for most Unix
+# systems, but may be altered if necessary.
+#
+#f77StandardInput = 5
+#f77StandardOutput = 6
+#f77FileInput = 10
+#f77AsciiOutput = 11
+#f77PlotOutput = 12
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.bc
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.bc 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.bc 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,57 @@
+displacement_units = m
+velocity_units = m/s
+force_units = newton
+ 1 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 2 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 6 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 7 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 10 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 14 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 16 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 21 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 22 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 3 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 4 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 5 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 8 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 11 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 13 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 15 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 19 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 24 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 1 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 2 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 4 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 8 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 10 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 18 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 19 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 20 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 23 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 3 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 5 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 6 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 7 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 12 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 13 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 16 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 17 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 26 1 1 0 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 1 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 5 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 6 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 8 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 12 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 15 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 18 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 21 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 27 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 2 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 3 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 4 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 7 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 11 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 14 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 17 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 20 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
+ 25 1 1 1 0.00000000E+00 0.00000000E+00 1.00000000E+00
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.connect
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.connect 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.connect 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,96 @@
+ 1 5 1 0 8 15 27 28
+ 2 5 1 0 8 27 18 28
+ 3 5 1 0 6 16 21 29
+ 4 5 1 0 8 28 23 19
+ 5 5 1 0 4 34 11 19
+ 6 5 1 0 7 14 22 32
+ 7 5 1 0 4 25 11 34
+ 8 5 1 0 2 33 14 20
+ 9 5 1 0 8 24 15 28
+ 10 5 1 0 5 13 26 30
+ 11 5 1 0 2 23 10 33
+ 12 5 1 0 6 12 16 29
+ 13 5 1 0 27 12 6 29
+ 14 5 1 0 15 5 12 30
+ 15 5 1 0 28 30 29 9
+ 16 5 1 0 27 28 15 29
+ 17 5 1 0 27 29 15 12
+ 18 5 1 0 15 29 28 30
+ 19 5 1 0 15 30 12 29
+ 20 5 1 0 18 21 1 31
+ 21 5 1 0 27 6 21 29
+ 22 5 1 0 28 29 31 9
+ 23 5 1 0 18 28 27 31
+ 24 5 1 0 18 31 27 21
+ 25 5 1 0 27 31 28 29
+ 26 5 1 0 27 29 21 31
+ 27 5 1 0 21 22 1 31
+ 28 5 1 0 16 7 22 32
+ 29 5 1 0 29 32 31 9
+ 30 5 1 0 21 29 16 31
+ 31 5 1 0 21 31 16 22
+ 32 5 1 0 16 31 29 32
+ 33 5 1 0 16 32 22 31
+ 34 5 1 0 23 33 2 20
+ 35 5 1 0 28 9 33 34
+ 36 5 1 0 19 34 20 4
+ 37 5 1 0 23 19 28 20
+ 38 5 1 0 23 20 28 33
+ 39 5 1 0 28 20 19 34
+ 40 5 1 0 28 34 33 20
+ 41 5 1 0 11 35 3 24
+ 42 5 1 0 34 9 35 28
+ 43 5 1 0 19 28 24 8
+ 44 5 1 0 11 19 34 24
+ 45 5 1 0 11 24 34 35
+ 46 5 1 0 34 24 19 28
+ 47 5 1 0 34 28 35 24
+ 48 5 1 0 22 10 1 31
+ 49 5 1 0 14 2 10 33
+ 50 5 1 0 32 33 31 9
+ 51 5 1 0 22 32 14 31
+ 52 5 1 0 22 31 14 10
+ 53 5 1 0 14 31 32 33
+ 54 5 1 0 14 33 10 31
+ 55 5 1 0 11 17 3 35
+ 56 5 1 0 25 7 17 32
+ 57 5 1 0 34 32 35 9
+ 58 5 1 0 11 34 25 35
+ 59 5 1 0 11 35 25 17
+ 60 5 1 0 25 35 34 32
+ 61 5 1 0 25 32 17 35
+ 62 5 1 0 14 32 7 25
+ 63 5 1 0 33 9 32 34
+ 64 5 1 0 20 34 25 4
+ 65 5 1 0 14 20 33 25
+ 66 5 1 0 14 25 33 32
+ 67 5 1 0 33 25 20 34
+ 68 5 1 0 33 34 32 25
+ 69 5 1 0 15 13 5 30
+ 70 5 1 0 24 3 13 35
+ 71 5 1 0 28 35 30 9
+ 72 5 1 0 15 28 24 30
+ 73 5 1 0 15 30 24 13
+ 74 5 1 0 24 30 28 35
+ 75 5 1 0 24 35 13 30
+ 76 5 1 0 26 17 7 32
+ 77 5 1 0 13 3 17 35
+ 78 5 1 0 30 35 32 9
+ 79 5 1 0 26 30 13 32
+ 80 5 1 0 26 32 13 17
+ 81 5 1 0 13 32 30 35
+ 82 5 1 0 13 35 17 32
+ 83 5 1 0 10 18 1 31
+ 84 5 1 0 23 8 18 28
+ 85 5 1 0 33 28 31 9
+ 86 5 1 0 10 33 23 31
+ 87 5 1 0 10 31 23 18
+ 88 5 1 0 23 31 33 28
+ 89 5 1 0 23 28 18 31
+ 90 5 1 0 16 26 7 32
+ 91 5 1 0 12 5 26 30
+ 92 5 1 0 29 30 32 9
+ 93 5 1 0 16 29 12 32
+ 94 5 1 0 16 32 12 26
+ 95 5 1 0 12 32 29 30
+ 96 5 1 0 12 30 26 32
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.coord
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.coord 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.coord 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,36 @@
+coord_units = m
+ 1 0.00000000E+00 0.00000000E+00 -2.40000000E+04
+ 2 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 3 2.40000000E+04 2.40000000E+04 0.00000000E+00
+ 4 2.40000000E+04 0.00000000E+00 0.00000000E+00
+ 5 2.40000000E+04 2.40000000E+04 -2.40000000E+04
+ 6 0.00000000E+00 2.40000000E+04 -2.40000000E+04
+ 7 0.00000000E+00 2.40000000E+04 0.00000000E+00
+ 8 2.40000000E+04 0.00000000E+00 -2.40000000E+04
+ 9 1.20000000E+04 1.20000000E+04 -1.20000000E+04
+ 10 0.00000000E+00 0.00000000E+00 -1.20000000E+04
+ 11 2.40000000E+04 1.20000000E+04 0.00000000E+00
+ 12 1.20000000E+04 2.40000000E+04 -2.40000000E+04
+ 13 2.40000000E+04 2.40000000E+04 -1.20000000E+04
+ 14 0.00000000E+00 1.20000000E+04 0.00000000E+00
+ 15 2.40000000E+04 1.20000000E+04 -2.40000000E+04
+ 16 0.00000000E+00 2.40000000E+04 -1.20000000E+04
+ 17 1.20000000E+04 2.40000000E+04 0.00000000E+00
+ 18 1.20000000E+04 0.00000000E+00 -2.40000000E+04
+ 19 2.40000000E+04 0.00000000E+00 -1.20000000E+04
+ 20 1.20000000E+04 0.00000000E+00 0.00000000E+00
+ 21 0.00000000E+00 1.20000000E+04 -2.40000000E+04
+ 22 0.00000000E+00 1.20000000E+04 -1.20000000E+04
+ 23 1.20000000E+04 0.00000000E+00 -1.20000000E+04
+ 24 2.40000000E+04 1.20000000E+04 -1.20000000E+04
+ 25 1.20000000E+04 1.20000000E+04 0.00000000E+00
+ 26 1.20000000E+04 2.40000000E+04 -1.20000000E+04
+ 27 1.20000000E+04 1.20000000E+04 -2.40000000E+04
+ 28 1.80000000E+04 6.00000000E+03 -1.80000000E+04
+ 29 6.00000000E+03 1.80000000E+04 -1.80000000E+04
+ 30 1.80000000E+04 1.80000000E+04 -1.80000000E+04
+ 31 6.00000000E+03 6.00000000E+03 -1.80000000E+04
+ 32 6.00000000E+03 1.80000000E+04 -6.00000000E+03
+ 33 6.00000000E+03 6.00000000E+03 -6.00000000E+03
+ 34 1.80000000E+04 6.00000000E+03 -6.00000000E+03
+ 35 1.80000000E+04 1.80000000E+04 -6.00000000E+03
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.fuldat
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.fuldat 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.fuldat 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,21 @@
+# Time steps for which full output is desired.
+#
+ 1
+ 2
+ 3
+ 4
+ 5
+ 6
+ 7
+ 8
+ 9
+ 10
+ 20
+ 30
+ 40
+ 50
+ 60
+ 70
+ 80
+ 90
+ 100
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.hist
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.hist 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.hist 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,22 @@
+# Sample time history definition file for TECTON.
+# Each history consists of two or more lines. The first line defines
+# the number of points defining the history and the default load
+# value for the history. Subsequent lines define time, load pairs.
+#
+# Dummy load history file. No load histories are actually required
+# for SCEC benchmark 1.
+#
+# The definitions are as follows:
+#
+# npoints = The number of points defining the given history.
+# defval = The value assigned to every point by default. The
+# default value is then overridden by load values.
+# time = The time value for a given point.
+# load = The load value for a given point.
+#
+# npoints(i) defval(i)
+# time(j,i) load(j,i), j=1,points(i), i=1,nhistories
+#
+ 2 0.0
+ 0.0 1.0
+ 4.0e8 1.0
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.prop
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.prop 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.prop 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,18 @@
+# Sample material properties definition file.
+# The number of properties for each material model have been
+# previously defined.
+# At present, the materials are assumed to be listed in order, so that
+# material number is determined by position in this file. This may
+# change in the future.
+# Also, the end of each material description is signified by setting
+# the endMaterial flag to 'True'.
+# Also at present, the materialType corresponds to a class, so the
+# first letter is capitalized.
+#
+# Material number 1
+materialType = 'IsotropicLinearMaxwellViscoelastic'
+density = 3000.0*kg/m**3 # Density
+youngsModulus = 7.5e10*Pa # Young's modulus
+poissonsRatio = 0.25 # Poisson's ratio
+viscosity = 1.0e18*Pa*s # Viscosity
+endMaterial = True
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.statevar
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.statevar 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.statevar 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,34 @@
+# Sample file defining which state variables are to be output for
+# the elastic and time dependent solutions.
+#
+# The istatout array specifies output options for each individual
+# state variable. At present there are a maximum of 24 possible
+# state variables, and this number may increase with the addition
+# of new material models. There are three types of state variable
+# output:
+#
+# 1 Total accumulated values for the current time step
+# 2 Incremental values from the previous step to the current
+# 3 Rates computed from the previous step to the current
+#
+# Present state variables occur in groups of 6, corresponding to
+# the number of stress/strain components, although this may change
+# in the future. The present groups are:
+#
+# 1-6: Cauchy stress
+# 7-12: Total strain
+# 13-18: Viscous strain
+# 18-24: Plastic strain
+#
+# Three lines of input are required, corresponding to the three
+# types of state variable output. For each line the user must
+# enter:
+# The number of state variables to output for this type (nstatout).
+# Note that the value of nstatout may be zero, in which case no
+# further output is needed for that line.
+# The state variables to output for this type (nstatout values).
+#
+#nstatout, istatout(i),i=1,nstatout
+ 12 1 2 3 4 5 6 7 8 9 10 11 12
+ 12 1 2 3 4 5 6 7 8 9 10 11 12
+ 0
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.time
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.time 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistest.time 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,64 @@
+# Sample time step group definition file for TECTON. This file
+# contains several integer and real parameters that must be defined
+# for each group. Note that the elastic solution corresponds to
+# group number 0 and must always be defined. Some of the parameters
+# have no meaning for the elastic solution but they are defined
+# anyway. In the fortran code, each of the parameters (other than n)
+# is stored in a separate array of dimension nintg+1, where nintg
+# is the number of time step groups. The definitions are as follows:
+#
+# n = time step group number. The elastic solution
+# corresponds to group number 0, and must always be
+# defined.
+# maxstp = the number of steps in the group (this is always 1 for
+# the elastic solution).
+# delt = time step size. This is a real value with units of
+# time.
+# alfa = amount of implicitness. This is a real dimensionless
+# parameter with values that can range from 0.0 (fully
+# explicit) to 1.0 (fully implicit). The value is
+# generally set to 0.5.
+# maxit = maximum number of equilibrium iterations before stiffness
+# matrix is reformed.
+# ntdinit= number of time steps between initial reformation of
+# stiffness matrix (before any iterations). A negative
+# value indicates that reformation should occur only once
+# for the first step in the group. A value of zero indicates
+# that reformation should never occur.
+# lgdef = large deformation solution flag. This is an integer
+# parameter with 3 possible values:
+# 0: linear strain
+# 1: large strain, but use only the linear contribution
+# to the stiffness matrix. This sometimes gives
+# better convergence.
+# 2: large strain using the nonlinear contribution to the
+# stiffness matrix.
+#
+# Convergence criteria for the nonlinear iterative solution:
+#
+# utol = convergence tolerance for displacements. This is a
+# dimensionless real value.
+# ftol = convergence tolerance for forces. This is a
+# dimensionless real value.
+# etol = convergence tolerance for energy. This is a
+# dimensionless real value.
+# itmax = maximum number of equilibrium iterations. This is an
+# integer quantity.
+#
+# Note that the convergence criteria defined below are appropriate for
+# a linear elastic problem. The values would be much different for a
+# nonlinear or time-dependent problem.
+#
+# The only parameter in this list requiring units is the time step
+# size, with units of time. I generally make sure everything is
+# converted to seconds, although I frequently specify time in years.
+#
+# Time step information for SCEC benchmark 1.
+#
+#
+time_units = year # This should be converted to seconds.
+#
+# n maxstp delt alfa maxit ntdinit lgdef utol ftol etol itmax
+#
+ 0 1 0.0d0 5.d-1 1001 4 0 1.0d+00 1.0d+0 1.0d+00 1
+ 1 100 0.1d0 5.d-1 10 -1 0 9.0d-01 9.0d-01 9.0d-01 2
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestGenmax.prop
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestGenmax.prop 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestGenmax.prop 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,23 @@
+# Sample material properties definition file.
+# The number of properties for each material model have been
+# previously defined.
+# At present, the materials are assumed to be listed in order, so that
+# material number is determined by position in this file. This may
+# change in the future.
+# Also, the end of each material description is signified by setting
+# the endMaterial flag to 'True'.
+# Also at present, the materialType corresponds to a class, so the
+# first letter is capitalized.
+#
+# Material number 1
+materialType = 'IsotropicLinearGenMaxwellViscoelastic'
+density = 3000.0*kg/m**3 # Density
+youngsModulus = 7.5e10*Pa # Young's modulus
+poissonsRatio = 0.25 # Poisson's ratio
+shearRatio1 = 0.3 # Shear Ratio 1
+viscosity1 = 1.0e20 # Viscosity 1
+shearRatio2 = 0.3 # Shear Ratio 2
+viscosity2 = 1.0e19 # Viscosity 2
+shearRatio3 = 0.4 # Shear Ratio 3
+viscosity3 = 1.0e18 # Viscosity 3
+endMaterial = True
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestLinmaxESF.prop
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestLinmaxESF.prop 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestLinmaxESF.prop 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,18 @@
+# Sample material properties definition file.
+# The number of properties for each material model have been
+# previously defined.
+# At present, the materials are assumed to be listed in order, so that
+# material number is determined by position in this file. This may
+# change in the future.
+# Also, the end of each material description is signified by setting
+# the endMaterial flag to 'True'.
+# Also at present, the materialType corresponds to a class, so the
+# first letter is capitalized.
+#
+# Material number 1
+materialType = 'IsotropicLinearMaxwellViscoelasticESF'
+density = 3000.0*kg/m**3 # Density
+youngsModulus = 7.5e10*Pa # Young's modulus
+poissonsRatio = 0.25 # Poisson's ratio
+viscosity = 1.0e18*Pa*s # Viscosity
+endMaterial = True
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestNL.time
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestNL.time 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestNL.time 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,64 @@
+# Sample time step group definition file for TECTON. This file
+# contains several integer and real parameters that must be defined
+# for each group. Note that the elastic solution corresponds to
+# group number 0 and must always be defined. Some of the parameters
+# have no meaning for the elastic solution but they are defined
+# anyway. In the fortran code, each of the parameters (other than n)
+# is stored in a separate array of dimension nintg+1, where nintg
+# is the number of time step groups. The definitions are as follows:
+#
+# n = time step group number. The elastic solution
+# corresponds to group number 0, and must always be
+# defined.
+# maxstp = the number of steps in the group (this is always 1 for
+# the elastic solution).
+# delt = time step size. This is a real value with units of
+# time.
+# alfa = amount of implicitness. This is a real dimensionless
+# parameter with values that can range from 0.0 (fully
+# explicit) to 1.0 (fully implicit). The value is
+# generally set to 0.5.
+# maxit = maximum number of equilibrium iterations before stiffness
+# matrix is reformed.
+# ntdinit= number of time steps between initial reformation of
+# stiffness matrix (before any iterations). A negative
+# value indicates that reformation should occur only once
+# for the first step in the group. A value of zero indicates
+# that reformation should never occur.
+# lgdef = large deformation solution flag. This is an integer
+# parameter with 3 possible values:
+# 0: linear strain
+# 1: large strain, but use only the linear contribution
+# to the stiffness matrix. This sometimes gives
+# better convergence.
+# 2: large strain using the nonlinear contribution to the
+# stiffness matrix.
+#
+# Convergence criteria for the nonlinear iterative solution:
+#
+# utol = convergence tolerance for displacements. This is a
+# dimensionless real value.
+# ftol = convergence tolerance for forces. This is a
+# dimensionless real value.
+# etol = convergence tolerance for energy. This is a
+# dimensionless real value.
+# itmax = maximum number of equilibrium iterations. This is an
+# integer quantity.
+#
+# Note that the convergence criteria defined below are appropriate for
+# a linear elastic problem. The values would be much different for a
+# nonlinear or time-dependent problem.
+#
+# The only parameter in this list requiring units is the time step
+# size, with units of time. I generally make sure everything is
+# converted to seconds, although I frequently specify time in years.
+#
+# Time step information for SCEC benchmark 1.
+#
+#
+time_units = year # This should be converted to seconds.
+#
+# n maxstp delt alfa maxit ntdinit lgdef utol ftol etol itmax
+#
+ 0 1 0.0d0 5.d-1 1001 4 0 1.0d+00 1.0d+0 1.0d+00 1
+ 1 100 0.1d0 5.d-1 1 1 0 9.0d-01 9.0d-01 9.0d-01 2
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlaw.prop
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlaw.prop 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlaw.prop 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,19 @@
+# Sample material properties definition file.
+# The number of properties for each material model have been
+# previously defined.
+# At present, the materials are assumed to be listed in order, so that
+# material number is determined by position in this file. This may
+# change in the future.
+# Also, the end of each material description is signified by setting
+# the endMaterial flag to 'True'.
+# Also at present, the materialType corresponds to a class, so the
+# first letter is capitalized.
+#
+# Material number 1
+materialType = 'IsotropicPowerLawMaxwellViscoelastic'
+density = 3000.0*kg/m**3 # Density
+youngsModulus = 7.5e10*Pa # Young's modulus
+poissonsRatio = 0.25 # Poisson's ratio
+powerLawExponent = 3.50 # Exponent
+viscosityCoefficient = 5.0e9*s*Pa**3.5 # EMHU
+endMaterial = True
Added: short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlawZT.prop
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlawZT.prop 2007-03-21 20:23:32 UTC (rev 6337)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/examples/lintet/vistest/vistestPowerlawZT.prop 2007-03-21 21:38:26 UTC (rev 6338)
@@ -0,0 +1,19 @@
+# Sample material properties definition file.
+# The number of properties for each material model have been
+# previously defined.
+# At present, the materials are assumed to be listed in order, so that
+# material number is determined by position in this file. This may
+# change in the future.
+# Also, the end of each material description is signified by setting
+# the endMaterial flag to 'True'.
+# Also at present, the materialType corresponds to a class, so the
+# first letter is capitalized.
+#
+# Material number 1
+materialType = 'IsotropicPowerLawMaxwellViscoelasticZT'
+density = 3000.0*kg/m**3 # Density
+youngsModulus = 7.5e10*Pa # Young's modulus
+poissonsRatio = 0.25 # Poisson's ratio
+powerLawExponent = 3.50 # Exponent
+viscosityCoefficient = 5.0e9*s*Pa**3.5 # EMHU
+endMaterial = True
More information about the cig-commits
mailing list