[cig-commits] r7201 - in short/3D/PyLith/branches/pylith-0.8: .
tests tests/sheartest1 tests/sheartest1/hex8
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed Jun 13 09:40:02 PDT 2007
Author: willic3
Date: 2007-06-13 09:40:02 -0700 (Wed, 13 Jun 2007)
New Revision: 7201
Added:
short/3D/PyLith/branches/pylith-0.8/tests/
short/3D/PyLith/branches/pylith-0.8/tests/axialtest1/
short/3D/PyLith/branches/pylith-0.8/tests/axialtest2/
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d-debug.cfg
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d.cfg
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.bc
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.connect
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.coord
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.fuldat
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.prop
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.statevar
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.time
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/quad4/
short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/tet4/
short/3D/PyLith/branches/pylith-0.8/tests/sheartest2/
Log:
Started creating tests for PyLith 0.8 that can be compared directly
against those for 1.0.
Added: short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d-debug.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d-debug.cfg 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d-debug.cfg 2007-06-13 16:40:02 UTC (rev 7201)
@@ -0,0 +1,17 @@
+[pylith3d.petsc]
+# 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 options
+[pylith3d]
+fileRoot = sheartest1
+asciiOutput = full
+ucdOutput = ascii
+# materialPropertiesInputFile = vistestLinmaxESF.prop
Added: short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d.cfg
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d.cfg 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/pylith3d.cfg 2007-06-13 16:40:02 UTC (rev 7201)
@@ -0,0 +1,115 @@
+[pylith3d.petsc]
+# PETSc options
+log_summary = true
+pc_type = bjacobi
+sub_pc_type = ilu
+ksp_monitor = true
+ksp_view = true
+ksp_rtol = 1.0e-11
+# start_in_debugger = true
+# debugger_pause = 100
+
+[pylith3d]
+fileRoot = sheartest1
+asciiOutput = full
+ucdOutput = ascii
+# materialPropertiesInputFile = vistestLinmaxESF.prop
+
+# 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 5.
+#
+# 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
+#f77UcdOutput = 13
Added: short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.bc
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.bc 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.bc 2007-06-13 16:40:02 UTC (rev 7201)
@@ -0,0 +1,56 @@
+# Boundary conditions for traction test problem.
+#
+displacement_units = m
+velocity_units = m/s
+force_units = newton
+#
+ 1 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 2 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 3 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 4 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
+ 7 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
+ 9 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 10 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 11 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
+ 13 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 14 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
+ 16 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 17 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
+ 19 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 20 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
+ 22 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 23 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 24 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 25 1 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 101 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 102 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 103 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 104 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 105 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 106 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 107 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 108 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 109 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 110 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 111 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 112 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 113 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 114 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 115 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 116 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 117 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 118 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 119 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 120 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 121 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 122 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 123 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 124 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
+ 125 1 1 1 0.00000000E+00 1.00000000E+00 0.00000000E+00
Added: short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.connect
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.connect 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.connect 2007-06-13 16:40:02 UTC (rev 7201)
@@ -0,0 +1,95 @@
+# Sample connectivity input file for 3D problem. The element number is
+# not strictly necessary, but makes it easier to examine the file.
+# This file defines the element type, material type, and infinite
+# element info as well as the connectivity.
+# The entries are as follows:
+# iel = element number.
+# ietype = element type from the following list:
+# 1 = linear hexahedron (8 nodes)
+# 2 = "wrick" -- linear hex with one set of collapsed nodes
+# (7 nodes)
+# 3 = wedge -- linear hex with two sets of collapsed nodes
+# (6 nodes)
+# 4 = pyramid -- linear hex with 4 nodes collapsed to a point
+# (5 nodes)
+# 5 = linear tetrahedron (4 nodes)
+# 6 = quadratic hexahedron (20 nodes)
+# 7 = quadratic "wrick" -- quadratic hex with 3 nodes
+# along one edge collapsed to a point (18 nodes)
+# 8 = quadratic wedge -- quadratic hex with 3 sets of
+# collapsed nodes (15 nodes)
+# 9 = quadratic pyramid -- quadratic hex with 9 nodes
+# collapsed to a point (13 nodes)
+# 10 = quadratic tetrahedron (10 nodes)
+# mat = material type for element.
+# infin = infinite domain element code.
+# ia(i) = node numbers for element.
+#
+# At present, comments are not allowed within the connectivity
+# information.
+#
+# N IETYPE MAT INFIN N1 N2 N3 N4 N5 N6 N7 N8
+ 1 1 1 0 1 26 27 2 6 31 32 7
+ 2 1 1 0 2 27 28 3 7 32 33 8
+ 3 1 1 0 3 28 29 4 8 33 34 9
+ 4 1 1 0 4 29 30 5 9 34 35 10
+ 5 1 1 0 6 31 32 7 11 36 37 12
+ 6 1 1 0 7 32 33 8 12 37 38 13
+ 7 1 1 0 8 33 34 9 13 38 39 14
+ 8 1 1 0 9 34 35 10 14 39 40 15
+ 9 1 1 0 11 36 37 12 16 41 42 17
+ 10 1 1 0 12 37 38 13 17 42 43 18
+ 11 1 1 0 13 38 39 14 18 43 44 19
+ 12 1 1 0 14 39 40 15 19 44 45 20
+ 13 1 1 0 16 41 42 17 21 46 47 22
+ 14 1 1 0 17 42 43 18 22 47 48 23
+ 15 1 1 0 18 43 44 19 23 48 49 24
+ 16 1 1 0 19 44 45 20 24 49 50 25
+ 17 1 1 0 26 51 52 27 31 56 57 32
+ 18 1 1 0 27 52 53 28 32 57 58 33
+ 19 1 1 0 28 53 54 29 33 58 59 34
+ 20 1 1 0 29 54 55 30 34 59 60 35
+ 21 1 1 0 31 56 57 32 36 61 62 37
+ 22 1 1 0 32 57 58 33 37 62 63 38
+ 23 1 1 0 33 58 59 34 38 63 64 39
+ 24 1 1 0 34 59 60 35 39 64 65 40
+ 25 1 1 0 36 61 62 37 41 66 67 42
+ 26 1 1 0 37 62 63 38 42 67 68 43
+ 27 1 1 0 38 63 64 39 43 68 69 44
+ 28 1 1 0 39 64 65 40 44 69 70 45
+ 29 1 1 0 41 66 67 42 46 71 72 47
+ 30 1 1 0 42 67 68 43 47 72 73 48
+ 31 1 1 0 43 68 69 44 48 73 74 49
+ 32 1 1 0 44 69 70 45 49 74 75 50
+ 33 1 1 0 51 76 77 52 56 81 82 57
+ 34 1 1 0 52 77 78 53 57 82 83 58
+ 35 1 1 0 53 78 79 54 58 83 84 59
+ 36 1 1 0 54 79 80 55 59 84 85 60
+ 37 1 1 0 56 81 82 57 61 86 87 62
+ 38 1 1 0 57 82 83 58 62 87 88 63
+ 39 1 1 0 58 83 84 59 63 88 89 64
+ 40 1 1 0 59 84 85 60 64 89 90 65
+ 41 1 1 0 61 86 87 62 66 91 92 67
+ 42 1 1 0 62 87 88 63 67 92 93 68
+ 43 1 1 0 63 88 89 64 68 93 94 69
+ 44 1 1 0 64 89 90 65 69 94 95 70
+ 45 1 1 0 66 91 92 67 71 96 97 72
+ 46 1 1 0 67 92 93 68 72 97 98 73
+ 47 1 1 0 68 93 94 69 73 98 99 74
+ 48 1 1 0 69 94 95 70 74 99 100 75
+ 49 1 1 0 76 101 102 77 81 106 107 82
+ 50 1 1 0 77 102 103 78 82 107 108 83
+ 51 1 1 0 78 103 104 79 83 108 109 84
+ 52 1 1 0 79 104 105 80 84 109 110 85
+ 53 1 1 0 81 106 107 82 86 111 112 87
+ 54 1 1 0 82 107 108 83 87 112 113 88
+ 55 1 1 0 83 108 109 84 88 113 114 89
+ 56 1 1 0 84 109 110 85 89 114 115 90
+ 57 1 1 0 86 111 112 87 91 116 117 92
+ 58 1 1 0 87 112 113 88 92 117 118 93
+ 59 1 1 0 88 113 114 89 93 118 119 94
+ 60 1 1 0 89 114 115 90 94 119 120 95
+ 61 1 1 0 91 116 117 92 96 121 122 97
+ 62 1 1 0 92 117 118 93 97 122 123 98
+ 63 1 1 0 93 118 119 94 98 123 124 99
+ 64 1 1 0 94 119 120 95 99 124 125 100
Added: short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.coord
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.coord 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.coord 2007-06-13 16:40:02 UTC (rev 7201)
@@ -0,0 +1,142 @@
+# Sample coordinate file for 3D problem. The tentative format is to
+# describe the units (distance) at the top of the file, followed by
+# the coordinate information. The node number is not strictly
+# necessary (and will not be stored), but it makes it easier to see
+# the node-location correspondence when looking at this file.
+# In fortran, the nodal coordinate array is dimensioned as:
+# x(nsd,numnp), where nsd is the number of spatial dimensions (3 for
+# this 3D example), and numnp is the number of nodes.
+#
+coord_units = m
+#
+# List the coordinates below.
+# Note: The current implementation does not allow comments within
+# the coordinate list.
+#
+# Node X-coord Y-coord Z-coord
+#
+ 1 -2.00000000E+02 -2.00000000E+02 -2.00000000E+02
+ 2 -2.00000000E+02 -2.00000000E+02 -1.00000000E+02
+ 3 -2.00000000E+02 -2.00000000E+02 0.00000000E+00
+ 4 -2.00000000E+02 -2.00000000E+02 1.00000000E+02
+ 5 -2.00000000E+02 -2.00000000E+02 2.00000000E+02
+ 6 -2.00000000E+02 -1.00000000E+02 -2.00000000E+02
+ 7 -2.00000000E+02 -1.00000000E+02 -1.00000000E+02
+ 8 -2.00000000E+02 -1.00000000E+02 0.00000000E+00
+ 9 -2.00000000E+02 -1.00000000E+02 1.00000000E+02
+ 10 -2.00000000E+02 -1.00000000E+02 2.00000000E+02
+ 11 -2.00000000E+02 0.00000000E+00 -2.00000000E+02
+ 12 -2.00000000E+02 0.00000000E+00 -1.00000000E+02
+ 13 -2.00000000E+02 0.00000000E+00 0.00000000E+00
+ 14 -2.00000000E+02 0.00000000E+00 1.00000000E+02
+ 15 -2.00000000E+02 0.00000000E+00 2.00000000E+02
+ 16 -2.00000000E+02 1.00000000E+02 -2.00000000E+02
+ 17 -2.00000000E+02 1.00000000E+02 -1.00000000E+02
+ 18 -2.00000000E+02 1.00000000E+02 0.00000000E+00
+ 19 -2.00000000E+02 1.00000000E+02 1.00000000E+02
+ 20 -2.00000000E+02 1.00000000E+02 2.00000000E+02
+ 21 -2.00000000E+02 2.00000000E+02 -2.00000000E+02
+ 22 -2.00000000E+02 2.00000000E+02 -1.00000000E+02
+ 23 -2.00000000E+02 2.00000000E+02 0.00000000E+00
+ 24 -2.00000000E+02 2.00000000E+02 1.00000000E+02
+ 25 -2.00000000E+02 2.00000000E+02 2.00000000E+02
+ 26 -1.00000000E+02 -2.00000000E+02 -2.00000000E+02
+ 27 -1.00000000E+02 -2.00000000E+02 -1.00000000E+02
+ 28 -1.00000000E+02 -2.00000000E+02 0.00000000E+00
+ 29 -1.00000000E+02 -2.00000000E+02 1.00000000E+02
+ 30 -1.00000000E+02 -2.00000000E+02 2.00000000E+02
+ 31 -1.00000000E+02 -1.00000000E+02 -2.00000000E+02
+ 32 -1.00000000E+02 -1.00000000E+02 -1.00000000E+02
+ 33 -1.00000000E+02 -1.00000000E+02 0.00000000E+00
+ 34 -1.00000000E+02 -1.00000000E+02 1.00000000E+02
+ 35 -1.00000000E+02 -1.00000000E+02 2.00000000E+02
+ 36 -1.00000000E+02 0.00000000E+00 -2.00000000E+02
+ 37 -1.00000000E+02 0.00000000E+00 -1.00000000E+02
+ 38 -1.00000000E+02 0.00000000E+00 0.00000000E+00
+ 39 -1.00000000E+02 0.00000000E+00 1.00000000E+02
+ 40 -1.00000000E+02 0.00000000E+00 2.00000000E+02
+ 41 -1.00000000E+02 1.00000000E+02 -2.00000000E+02
+ 42 -1.00000000E+02 1.00000000E+02 -1.00000000E+02
+ 43 -1.00000000E+02 1.00000000E+02 0.00000000E+00
+ 44 -1.00000000E+02 1.00000000E+02 1.00000000E+02
+ 45 -1.00000000E+02 1.00000000E+02 2.00000000E+02
+ 46 -1.00000000E+02 2.00000000E+02 -2.00000000E+02
+ 47 -1.00000000E+02 2.00000000E+02 -1.00000000E+02
+ 48 -1.00000000E+02 2.00000000E+02 0.00000000E+00
+ 49 -1.00000000E+02 2.00000000E+02 1.00000000E+02
+ 50 -1.00000000E+02 2.00000000E+02 2.00000000E+02
+ 51 0.00000000E+00 -2.00000000E+02 -2.00000000E+02
+ 52 0.00000000E+00 -2.00000000E+02 -1.00000000E+02
+ 53 0.00000000E+00 -2.00000000E+02 0.00000000E+00
+ 54 0.00000000E+00 -2.00000000E+02 1.00000000E+02
+ 55 0.00000000E+00 -2.00000000E+02 2.00000000E+02
+ 56 0.00000000E+00 -1.00000000E+02 -2.00000000E+02
+ 57 0.00000000E+00 -1.00000000E+02 -1.00000000E+02
+ 58 0.00000000E+00 -1.00000000E+02 0.00000000E+00
+ 59 0.00000000E+00 -1.00000000E+02 1.00000000E+02
+ 60 0.00000000E+00 -1.00000000E+02 2.00000000E+02
+ 61 0.00000000E+00 0.00000000E+00 -2.00000000E+02
+ 62 0.00000000E+00 0.00000000E+00 -1.00000000E+02
+ 63 0.00000000E+00 0.00000000E+00 0.00000000E+00
+ 64 0.00000000E+00 0.00000000E+00 1.00000000E+02
+ 65 0.00000000E+00 0.00000000E+00 2.00000000E+02
+ 66 0.00000000E+00 1.00000000E+02 -2.00000000E+02
+ 67 0.00000000E+00 1.00000000E+02 -1.00000000E+02
+ 68 0.00000000E+00 1.00000000E+02 0.00000000E+00
+ 69 0.00000000E+00 1.00000000E+02 1.00000000E+02
+ 70 0.00000000E+00 1.00000000E+02 2.00000000E+02
+ 71 0.00000000E+00 2.00000000E+02 -2.00000000E+02
+ 72 0.00000000E+00 2.00000000E+02 -1.00000000E+02
+ 73 0.00000000E+00 2.00000000E+02 0.00000000E+00
+ 74 0.00000000E+00 2.00000000E+02 1.00000000E+02
+ 75 0.00000000E+00 2.00000000E+02 2.00000000E+02
+ 76 1.00000000E+02 -2.00000000E+02 -2.00000000E+02
+ 77 1.00000000E+02 -2.00000000E+02 -1.00000000E+02
+ 78 1.00000000E+02 -2.00000000E+02 0.00000000E+00
+ 79 1.00000000E+02 -2.00000000E+02 1.00000000E+02
+ 80 1.00000000E+02 -2.00000000E+02 2.00000000E+02
+ 81 1.00000000E+02 -1.00000000E+02 -2.00000000E+02
+ 82 1.00000000E+02 -1.00000000E+02 -1.00000000E+02
+ 83 1.00000000E+02 -1.00000000E+02 0.00000000E+00
+ 84 1.00000000E+02 -1.00000000E+02 1.00000000E+02
+ 85 1.00000000E+02 -1.00000000E+02 2.00000000E+02
+ 86 1.00000000E+02 0.00000000E+00 -2.00000000E+02
+ 87 1.00000000E+02 0.00000000E+00 -1.00000000E+02
+ 88 1.00000000E+02 0.00000000E+00 0.00000000E+00
+ 89 1.00000000E+02 0.00000000E+00 1.00000000E+02
+ 90 1.00000000E+02 0.00000000E+00 2.00000000E+02
+ 91 1.00000000E+02 1.00000000E+02 -2.00000000E+02
+ 92 1.00000000E+02 1.00000000E+02 -1.00000000E+02
+ 93 1.00000000E+02 1.00000000E+02 0.00000000E+00
+ 94 1.00000000E+02 1.00000000E+02 1.00000000E+02
+ 95 1.00000000E+02 1.00000000E+02 2.00000000E+02
+ 96 1.00000000E+02 2.00000000E+02 -2.00000000E+02
+ 97 1.00000000E+02 2.00000000E+02 -1.00000000E+02
+ 98 1.00000000E+02 2.00000000E+02 0.00000000E+00
+ 99 1.00000000E+02 2.00000000E+02 1.00000000E+02
+ 100 1.00000000E+02 2.00000000E+02 2.00000000E+02
+ 101 2.00000000E+02 -2.00000000E+02 -2.00000000E+02
+ 102 2.00000000E+02 -2.00000000E+02 -1.00000000E+02
+ 103 2.00000000E+02 -2.00000000E+02 0.00000000E+00
+ 104 2.00000000E+02 -2.00000000E+02 1.00000000E+02
+ 105 2.00000000E+02 -2.00000000E+02 2.00000000E+02
+ 106 2.00000000E+02 -1.00000000E+02 -2.00000000E+02
+ 107 2.00000000E+02 -1.00000000E+02 -1.00000000E+02
+ 108 2.00000000E+02 -1.00000000E+02 0.00000000E+00
+ 109 2.00000000E+02 -1.00000000E+02 1.00000000E+02
+ 110 2.00000000E+02 -1.00000000E+02 2.00000000E+02
+ 111 2.00000000E+02 0.00000000E+00 -2.00000000E+02
+ 112 2.00000000E+02 0.00000000E+00 -1.00000000E+02
+ 113 2.00000000E+02 0.00000000E+00 0.00000000E+00
+ 114 2.00000000E+02 0.00000000E+00 1.00000000E+02
+ 115 2.00000000E+02 0.00000000E+00 2.00000000E+02
+ 116 2.00000000E+02 1.00000000E+02 -2.00000000E+02
+ 117 2.00000000E+02 1.00000000E+02 -1.00000000E+02
+ 118 2.00000000E+02 1.00000000E+02 0.00000000E+00
+ 119 2.00000000E+02 1.00000000E+02 1.00000000E+02
+ 120 2.00000000E+02 1.00000000E+02 2.00000000E+02
+ 121 2.00000000E+02 2.00000000E+02 -2.00000000E+02
+ 122 2.00000000E+02 2.00000000E+02 -1.00000000E+02
+ 123 2.00000000E+02 2.00000000E+02 0.00000000E+00
+ 124 2.00000000E+02 2.00000000E+02 1.00000000E+02
+ 125 2.00000000E+02 2.00000000E+02 2.00000000E+02
Added: short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.fuldat
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.fuldat 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.fuldat 2007-06-13 16:40:02 UTC (rev 7201)
@@ -0,0 +1,12 @@
+# Time steps for which full output is desired for simple shear test.
+#
+ 1
+ 2
+ 3
+ 4
+ 5
+ 6
+ 7
+ 8
+ 9
+ 10
Added: short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.prop
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.prop 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.prop 2007-06-13 16:40:02 UTC (rev 7201)
@@ -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/tests/sheartest1/hex8/sheartest1.statevar
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.statevar 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.statevar 2007-06-13 16:40:02 UTC (rev 7201)
@@ -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/tests/sheartest1/hex8/sheartest1.time
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.time 2007-06-13 16:19:02 UTC (rev 7200)
+++ short/3D/PyLith/branches/pylith-0.8/tests/sheartest1/hex8/sheartest1.time 2007-06-13 16:40:02 UTC (rev 7201)
@@ -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.0d+0 5.d-1 1001 4 0 1.0d+00 1.0d+0 1.0d+00 1
+ 1 10 0.1d+0 5.d-1 1 1 0 1.0d-03 1.0d-3 1.0d-06 4
More information about the cig-commits
mailing list