[cig-commits] r4403 -
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive
baagaard at geodynamics.org
baagaard at geodynamics.org
Wed Aug 23 12:39:30 PDT 2006
Author: baagaard
Date: 2006-08-23 12:39:30 -0700 (Wed, 23 Aug 2006)
New Revision: 4403
Removed:
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.keyval
Modified:
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.bc
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.coord
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.fuldat
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.par
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.prop
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.statevar
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.time
short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/tutor.py
Log:
Updated to split cube problem. Fixed units in .par file.
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.bc
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.bc 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.bc 2006-08-23 19:39:30 UTC (rev 4403)
@@ -1,4 +1,4 @@
-# Boundary conditions for split block problem.
+# Boundary conditions for split cube problem.
#
displacement_units = m
velocity_units = m/s
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.coord
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.coord 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.coord 2006-08-23 19:39:30 UTC (rev 4403)
@@ -1,4 +1,4 @@
-coord_units = km
+coord_units = m
1 5.00000000E+02 0.00000000E+00 0.00000000E+00
2 1.00000000E+03 0.00000000E+00 0.00000000E+00
3 0.00000000E+00 1.00000000E+03 1.00000000E+03
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.fuldat
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.fuldat 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.fuldat 2006-08-23 19:39:30 UTC (rev 4403)
@@ -1,5 +1,9 @@
-# Time steps for which full output is desired for SCEC benchmark 5.
+# File containing time steps at which full output is desired.
#
- 10
- 50
- 100
+# Comment lines begin with '#'
+#
+# Note: Time step 0 (elastic solution) is always included in the
+# output.
+#
+# List the time steps, one per line.
+#
Deleted: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.keyval
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.keyval 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.keyval 2006-08-23 19:39:30 UTC (rev 4403)
@@ -1,98 +0,0 @@
-# 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
-#prestressAutoChangeElasticProperties = 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
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.par
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.par 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.par 2006-08-23 19:39:30 UTC (rev 4403)
@@ -15,7 +15,7 @@
#
# cscale ix iy iz
#-----------------------------------------------------------------------
- 1000.0d0 1 2 3
+ 1.0d0 1 2 3
#-----------------------------------------------------------------------
#
# Coordinate units:
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.prop
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.prop 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.prop 2006-08-23 19:39:30 UTC (rev 4403)
@@ -1,14 +1,27 @@
-# 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.
+# File containing material properties.
#
+# Comment lines begin with '#'
+#
+# The material type and material property values are specified using a
+# "keyword = value" syntax. The keywords for the different material
+# types are given below. The units for values with dimensions should
+# follow the value as shown in the examples below.
+#
+# Materials and keywords:
+# Isotropic linear elastic
+# materialType ='IsotropicLinearElastic'
+# density
+# youngsModulus
+# poissonsRatio
+# endMaterial ='True' (flag indicating end of material)
+# Isotropic linear maxwell viscoelastic
+# materialType ='IsotropicLinearMaxwellViscoelastic'
+# density
+# youngsModulus
+# poissonsRatio
+# viscosity
+# endMaterial ='True' (flag indicating end of material)
+#
# Material number 1
materialType = 'IsotropicLinearElastic'
density = 3000.0*kg/m**3 # Density
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.statevar
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.statevar 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.statevar 2006-08-23 19:39:30 UTC (rev 4403)
@@ -1,34 +1,23 @@
-# Sample file defining which state variables are to be output for
-# the elastic and time dependent solutions.
+# File specifying which state variables to output.
#
-# 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:
+# Comment lines begin with '#'
#
-# 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
+# State variables occur in groups of 6, corresponding to the number of
+# stress/strain components. The present groups are:
+# 1-6: Cauchy stress
+# 7-12: Total strain
+# 13-18: Viscous strain
+# 18-24: Plastic strain
#
-# 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:
+# Lines:
+# (1) Total accumulated values for the current time step
+# (2) Incremental values (previous to current)
+# (3) Rate values (previous to current)
#
-# 1-6: Cauchy stress
-# 7-12: Total strain
-# 13-18: Viscous strain
-# 18-24: Plastic strain
+# Columns (per line):
+# (1) Number of state variables to output (0 ≤ value ≤ 24)
+# (2)+ State variable number to output (1 ≤ value ≤ 24)
#
-# 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
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.time
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.time 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/splitcube.time 2006-08-23 19:39:30 UTC (rev 4403)
@@ -1,63 +1,40 @@
-# 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:
+# File containing time stepping parameters.
#
-# 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.
+# Comment lines begin with '#'
#
-# Convergence criteria for the nonlinear iterative solution:
+# First, specify units used in values with dimensions of time.
#
-# 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.
+time_units = year
#
-# 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 stepping parameters are given in groups. The elastic solution
+# corresponds to group 0 and must always be defined. Although some of
+# the parameters do not have any meaning for the elastic solution,
+# they must be present anyway.
#
-# Time step information for SCEC benchmark 1.
+# Columns:
+# (1) Time step group number (=0 for elastic solution).
+# (2) The number of time steps in the group (=1 for elastic solution).
+# (3) Time step size (given in units of time_units).
+# (4) Amount of implicitness. Real dimensionless parameter that
+# ranges from 0.0 (fully explicit) to 1.0 (fully implicit). The
+# value is generally set to 0.5.
+# (5) Maximum number of equilibrium iterations before stiffness
+# matrix is reformed.
+# (6) Number of time steps between initial reformation of stiffness
+# matrix.
+# <0 Indicates that reformation should occur only for the first
+# step in each time step group.
+# =0 Indicates that reformation should never occur.
+# (7) Large deformation solution flag
+# =0 Linear strain
+# =1 Large strain but use only linear contribution to the
+# stiffness matrix (sometimes results in better convergence)
+# =2 Large strain and use nonlinear contribution to the
+# stiffness matrix
+# (8) Convergece tolerance for displacements (dimensionless value)
+# (9) Convergece tolerance for forces (dimensionless value)
+# (10) Convergece tolerance for energy (dimensionless value)
+# (11) Maximum number of equilibrium iterations
#
-#
-time_units = year # This should be converted to seconds.
-#
-# n maxstp delt alfa maxit ntdinit lgdef utol ftol etol itmax
-#
0 1 0.e+0 5.e-1 1001 4 0 1.0e+00 1.0e+0 1.0e+00 1
Modified: short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/tutor.py
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/tutor.py 2006-08-23 03:34:03 UTC (rev 4402)
+++ short/3D/PyLith/branches/pylith-0.8/tutorials/splitcube/archive/tutor.py 2006-08-23 19:39:30 UTC (rev 4403)
@@ -67,17 +67,17 @@
fragments = [ [root1], [".0"],
[".gmesh", ".mesh", ".mesh.split"],
["", \
- ".time.00000", ".time.00010", ".time.00050", ".time.00100"],
+ ".time.00000"],
[".inp"] ]
files['output'] += assembleStrings(fragments)
elif "viz1" == step:
root1 = root + "_1"
fragments = [ [root1], [".0"],
- [".mesh", ".mesh.time.00010"], [".inp"] ]
+ [".mesh", ".mesh.time.00000"], [".inp"] ]
files['input'] += assembleStrings(fragments)
fragments = [ [root1], [".0"],
- ["t00010"], [".inp"] ]
+ ["t00000"], [".inp"] ]
files['output'] += assembleStrings(fragments)
elif "run2" == step:
fragments = [ [root], [".coord", ".connect", ".split", ".bc",
@@ -98,17 +98,17 @@
fragments = [ [root1], [".0", ".1"],
[".gmesh", ".mesh", ".mesh.split"],
["", \
- ".time.00000", ".time.00010", ".time.00050", ".time.00100"],
+ ".time.00000"],
[".inp"] ]
files['output'] += assembleStrings(fragments)
elif "viz2" == step:
root1 = root + "_2"
fragments = [ [root1], [".0", ".1"],
- [".mesh", ".mesh.time.00010"], [".inp"] ]
+ [".mesh", ".mesh.time.00000"], [".inp"] ]
files['input'] += assembleStrings(fragments)
fragments = [ [root1], [".0", ".1"],
- [".t00010"], [".inp"] ]
+ [".t00000"], [".inp"] ]
files['output'] += assembleStrings(fragments)
print files
else:
More information about the cig-commits
mailing list