[cig-commits] r9214 - short/3D/PyLith/trunk/tests/3d/hex8
willic3 at geodynamics.org
willic3 at geodynamics.org
Sun Feb 3 08:23:04 PST 2008
Author: willic3
Date: 2008-02-03 08:23:04 -0800 (Sun, 03 Feb 2008)
New Revision: 9214
Added:
short/3D/PyLith/trunk/tests/3d/hex8/axial-small-tractions.cfg
short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1-tract-8cells.cfg
short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1a-tract-8cells.cfg
short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax2-tract-8cells.cfg
short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax3-tract-8cells.cfg
short/3D/PyLith/trunk/tests/3d/hex8/axialtractions-small.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/axialxyz-small.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1a.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell2.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell3.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/pylithapp.cfg
Modified:
short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am
short/3D/PyLith/trunk/tests/3d/hex8/eightcells.mesh
Log:
Added files to test generalized Maxwell model.
This directory is still a mess, and needs to be reorganized.
Things still to add:
Results from FEAP.
Plots showing fit between PyLith and FEAP results.
Analytical solution for regular Maxwell model uniaxial traction test.
Modified: short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am 2008-02-01 21:06:52 UTC (rev 9213)
+++ short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am 2008-02-03 16:23:04 UTC (rev 9214)
@@ -21,11 +21,23 @@
box.exo \
matprops.spatialdb \
matprops-maxwell.spatialdb \
+ axial-small-tractions.cfg \
axialelasticisotropic.cfg \
+ axialgenmax1-tract-8cells.cfg \
+ axialgenmax1a-tract-8cells.cfg \
+ axialgenmax2-tract-8cells.cfg \
+ axialgenmax3-tract-8cells.cfg \
axialmaxwellisotropic.cfg \
+ axialtractions-small.spatialdb \
+ axialxyz.spatialdb \
axialx.spatialdb \
axialy.spatialdb \
axialz.spatialdb \
+ genmaxwell1.spatialdb \
+ genmaxwell1a.spatialdb \
+ genmaxwell2.spatialdb \
+ genmaxwell3.spatialdb \
+ pylithapp.cfg \
shearelasticisotropic.cfg \
shearplanestrain.cfg \
shearmaxwellisotropic.cfg \
@@ -46,6 +58,50 @@
noinst_TMP = \
axialelasticisotropic_t0.vtk \
axialelasticisotropic_t1.vtk \
+ axialgenmax1-tract-8cells_t0000000.vtk \
+ axialgenmax1-tract-8cells_t11045160000000.vtk \
+ axialgenmax1-tract-8cells_t12623040000000.vtk \
+ axialgenmax1-tract-8cells_t14200920000000.vtk \
+ axialgenmax1-tract-8cells_t1577880000000.vtk \
+ axialgenmax1-tract-8cells_t15778800000000.vtk \
+ axialgenmax1-tract-8cells_t3155760000000.vtk \
+ axialgenmax1-tract-8cells_t4733640000000.vtk \
+ axialgenmax1-tract-8cells_t6311520000000.vtk \
+ axialgenmax1-tract-8cells_t7889400000000.vtk \
+ axialgenmax1-tract-8cells_t9467280000000.vtk \
+ axialgenmax1a-tract-8cells_t0000000.vtk \
+ axialgenmax1a-tract-8cells_t11045160000000.vtk \
+ axialgenmax1a-tract-8cells_t12623040000000.vtk \
+ axialgenmax1a-tract-8cells_t14200920000000.vtk \
+ axialgenmax1a-tract-8cells_t1577880000000.vtk \
+ axialgenmax1a-tract-8cells_t15778800000000.vtk \
+ axialgenmax1a-tract-8cells_t3155760000000.vtk \
+ axialgenmax1a-tract-8cells_t4733640000000.vtk \
+ axialgenmax1a-tract-8cells_t6311520000000.vtk \
+ axialgenmax1a-tract-8cells_t7889400000000.vtk \
+ axialgenmax1a-tract-8cells_t9467280000000.vtk \
+ axialgenmax2-tract-8cells_t0000000.vtk \
+ axialgenmax2-tract-8cells_t11045160000000.vtk \
+ axialgenmax2-tract-8cells_t12623040000000.vtk \
+ axialgenmax2-tract-8cells_t14200920000000.vtk \
+ axialgenmax2-tract-8cells_t1577880000000.vtk \
+ axialgenmax2-tract-8cells_t15778800000000.vtk \
+ axialgenmax2-tract-8cells_t3155760000000.vtk \
+ axialgenmax2-tract-8cells_t4733640000000.vtk \
+ axialgenmax2-tract-8cells_t6311520000000.vtk \
+ axialgenmax2-tract-8cells_t7889400000000.vtk \
+ axialgenmax2-tract-8cells_t9467280000000.vtk \
+ axialgenmax3-tract-8cells_t0000000.vtk \
+ axialgenmax3-tract-8cells_t11045160000000.vtk \
+ axialgenmax3-tract-8cells_t12623040000000.vtk \
+ axialgenmax3-tract-8cells_t14200920000000.vtk \
+ axialgenmax3-tract-8cells_t1577880000000.vtk \
+ axialgenmax3-tract-8cells_t15778800000000.vtk \
+ axialgenmax3-tract-8cells_t3155760000000.vtk \
+ axialgenmax3-tract-8cells_t4733640000000.vtk \
+ axialgenmax3-tract-8cells_t6311520000000.vtk \
+ axialgenmax3-tract-8cells_t7889400000000.vtk \
+ axialgenmax3-tract-8cells_t9467280000000.vtk \
axialmaxwellisotropic_t0.vtk \
axialmaxwellisotropic_t1.vtk \
shearelasticisotropic_t0.vtk \
Added: short/3D/PyLith/trunk/tests/3d/hex8/axial-small-tractions.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/axial-small-tractions.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/axial-small-tractions.cfg 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,70 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+debug = 1
+importer = pylith.meshio.MeshIOAscii
+
+[pylithapp.mesh_generator.importer]
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc.z_pos = pylith.bc.Neumann
+
+[pylithapp.timedependent.bc.z_pos]
+label = z_pos
+up_dir=[0.0, 1.0, 0.0]
+db.label = Neumann BC +z edge
+db.iohandler.filename = axialtractions-small.spatialdb
+db.query_type = nearest
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+id = 2
+label = z_neg
+db = spatialdata.spatialdb.SimpleDB
+db.label = Dirichlet BC -z edge
+db.iohandler.filename = axialxyz-small.spatialdb
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0]
+id = 3
+label = x_neg
+db = spatialdata.spatialdb.SimpleDB
+db.label = Dirichlet BC -x edge
+db.iohandler.filename = axialxyz-small.spatialdb
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0]
+id = 4
+label = x_pos
+db = spatialdata.spatialdb.SimpleDB
+db.label = Dirichlet BC +x edge
+db.iohandler.filename = axialxyz-small.spatialdb
+
+[pylithapp.timedependent.bc.y_neg]
+fixed_dof = [1]
+id = 5
+label = y_neg
+db = spatialdata.spatialdb.SimpleDB
+db.label = Dirichlet BC -y edge
+db.iohandler.filename = axialxyz-small.spatialdb
+
+[pylithapp.timedependent.bc.y_pos]
+fixed_dof = [1]
+id = 6
+label = y_pos
+db = spatialdata.spatialdb.SimpleDB
+db.label = Dirichlet BC +y edge
+db.iohandler.filename = axialxyz-small.spatialdb
+
Added: short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1-tract-8cells.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1-tract-8cells.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1-tract-8cells.cfg 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,38 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.importer]
+filename = eightcells.mesh
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.5*year
+default_dt = 0.05*year
+dimension = 3
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCSixSides
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+material = pylith.materials.GenMaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.material]
+label = Generalized Maxwell viscoelastic material
+id = 0
+db.iohandler.filename = genmaxwell1.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output.writer]
+filename = axialgenmax1-tract-8cells.vtk
Added: short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1a-tract-8cells.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1a-tract-8cells.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax1a-tract-8cells.cfg 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,38 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.importer]
+filename = eightcells.mesh
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.5*year
+default_dt = 0.05*year
+dimension = 3
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCSixSides
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+material = pylith.materials.GenMaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.material]
+label = Generalized Maxwell viscoelastic material
+id = 0
+db.iohandler.filename = genmaxwell1a.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output.writer]
+filename = axialgenmax1a-tract-8cells.vtk
Added: short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax2-tract-8cells.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax2-tract-8cells.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax2-tract-8cells.cfg 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,38 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.importer]
+filename = eightcells.mesh
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.5*year
+default_dt = 0.05*year
+dimension = 3
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCSixSides
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+material = pylith.materials.GenMaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.material]
+label = Generalized Maxwell viscoelastic material
+id = 0
+db.iohandler.filename = genmaxwell2.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output.writer]
+filename = axialgenmax2-tract-8cells.vtk
Added: short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax3-tract-8cells.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax3-tract-8cells.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/axialgenmax3-tract-8cells.cfg 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,38 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.importer]
+filename = eightcells.mesh
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.5*year
+default_dt = 0.05*year
+dimension = 3
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCSixSides
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+material = pylith.materials.GenMaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.material]
+label = Generalized Maxwell viscoelastic material
+id = 0
+db.iohandler.filename = genmaxwell3.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output.writer]
+filename = axialgenmax3-tract-8cells.vtk
Added: short/3D/PyLith/trunk/tests/3d/hex8/axialtractions-small.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/axialtractions-small.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/axialtractions-small.spatialdb 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = horiz-shear-traction vert-shear-traction normal-traction
+ value-units = Pa Pa Pa
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 +1.0 0.0 0.0 1.0e3
Added: short/3D/PyLith/trunk/tests/3d/hex8/axialxyz-small.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/axialxyz-small.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/axialxyz-small.spatialdb 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,15 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = dof-0 dof-1 dof-2
+ value-units = m m m
+ num-locs = 2
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 +1.0 0.0 0.0 0.1
+0.0 0.0 -1.0 0.0 0.0 0.0
Modified: short/3D/PyLith/trunk/tests/3d/hex8/eightcells.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/eightcells.mesh 2008-02-01 21:06:52 UTC (rev 9213)
+++ short/3D/PyLith/trunk/tests/3d/hex8/eightcells.mesh 2008-02-03 16:23:04 UTC (rev 9214)
@@ -74,4 +74,36 @@
6 7 8 15 16 17 24 25 26
}
}
+ group = {
+ name = y_neg
+ type = vertices
+ count = 9
+ indices = {
+ 0 3 6 9 12 15 18 21 24
+ }
+ }
+ group = {
+ name = y_pos
+ type = vertices
+ count = 9
+ indices = {
+ 2 5 8 11 14 17 20 23 26
+ }
+ }
+ group = {
+ name = z_neg
+ type = vertices
+ count = 9
+ indices = {
+ 0 1 2 3 4 5 6 7 8
+ }
+ }
+ group = {
+ name = z_pos
+ type = vertices
+ count = 9
+ indices = {
+ 18 19 20 21 22 23 24 25 26
+ }
+ }
}
Added: short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1.spatialdb 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 9
+ value-names = density vs vp shear_ratio_1 shear_ratio_2 shear_ratio_3 viscosity_1 viscosity_2 viscosity_3
+ value-units = kg/m^3 m/s m/s None None None Pa*s Pa*s Pa*s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258 1.0 0.0 0.0 1.0e18 1.0 1.0
Added: short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1a.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1a.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell1a.spatialdb 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 9
+ value-names = density vs vp shear_ratio_1 shear_ratio_2 shear_ratio_3 viscosity_1 viscosity_2 viscosity_3
+ value-units = kg/m^3 m/s m/s None None None Pa*s Pa*s Pa*s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258 0.5 0.0 0.0 1.0e18 1.0 1.0
Added: short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell2.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell2.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell2.spatialdb 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 9
+ value-names = density vs vp shear_ratio_1 shear_ratio_2 shear_ratio_3 viscosity_1 viscosity_2 viscosity_3
+ value-units = kg/m^3 m/s m/s None None None Pa*s Pa*s Pa*s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258 0.5 0.5 0.0 1.0e18 1.0e19 1.0
Added: short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell3.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell3.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/genmaxwell3.spatialdb 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 9
+ value-names = density vs vp shear_ratio_1 shear_ratio_2 shear_ratio_3 viscosity_1 viscosity_2 viscosity_3
+ value-units = kg/m^3 m/s m/s None None None Pa*s Pa*s Pa*s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258 0.3 0.3 0.3 1.0e18 1.0e17 1.0e19
Added: short/3D/PyLith/trunk/tests/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/3d/hex8/pylithapp.cfg 2008-02-03 16:23:04 UTC (rev 9214)
@@ -0,0 +1,34 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshioascii = 1
+meshiocubit = 1
+meshimporter = 1
+homogeneous = 1
+implicitelasticity = 1
+quadrature3d = 1
+fiatsimplex = 1
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+# Set the solver options.
+[pylithapp.petsc]
+ksp_rtol = 1.0e-12
+ksp_max_it = 100
+ksp_gmres_restart = 50
+pc_type = bjacobi
+ksp_monitor = true
+ksp_view = true
+log_summary = true
+# start_in_debugger = true
More information about the cig-commits
mailing list