[cig-commits] r7210 - short/3D/PyLith/trunk/tests/3d/hex8
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed Jun 13 12:13:33 PDT 2007
Author: willic3
Date: 2007-06-13 12:13:33 -0700 (Wed, 13 Jun 2007)
New Revision: 7210
Added:
short/3D/PyLith/trunk/tests/3d/hex8/matprops-twocells.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/shearelastic-twocells.cfg
short/3D/PyLith/trunk/tests/3d/hex8/shearxyz-twocells.spatialdb
short/3D/PyLith/trunk/tests/3d/hex8/twocells.mesh
Modified:
short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am
short/3D/PyLith/trunk/tests/3d/hex8/matprops-maxwell.spatialdb
Log:
Added two cells example for debugging.
Fixed missing viscosity for Maxwell model.
Modified: short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am 2007-06-13 19:03:16 UTC (rev 7209)
+++ short/3D/PyLith/trunk/tests/3d/hex8/Makefile.am 2007-06-13 19:13:33 UTC (rev 7210)
@@ -26,9 +26,13 @@
axialz.spatialdb \
matprops.spatialdb \
matprops-maxwell.spatialdb \
+ matprops-twocells.spatialdb \
shearelasticisotropic.cfg \
+ shearelastic-twocells.cfg \
shearmaxwellisotropic.cfg \
- shearxyz.spatialdb
+ shearxyz.spatialdb \
+ shearxyz-twocells.spatialdb \
+ twocells.mesh
noinst_TMP = \
axialelasticisotropic.vtk \
Modified: short/3D/PyLith/trunk/tests/3d/hex8/matprops-maxwell.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/matprops-maxwell.spatialdb 2007-06-13 19:03:16 UTC (rev 7209)
+++ short/3D/PyLith/trunk/tests/3d/hex8/matprops-maxwell.spatialdb 2007-06-13 19:13:33 UTC (rev 7210)
@@ -11,4 +11,4 @@
space-dim = 3
}
}
-0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258
+0.0 0.0 0.0 2700.0 3333.333333333333 5773.502691896258 1.0e18
Added: short/3D/PyLith/trunk/tests/3d/hex8/matprops-twocells.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/matprops-twocells.spatialdb 2007-06-13 19:03:16 UTC (rev 7209)
+++ short/3D/PyLith/trunk/tests/3d/hex8/matprops-twocells.spatialdb 2007-06-13 19:13:33 UTC (rev 7210)
@@ -0,0 +1,14 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = density vs vp
+ value-units = kg/m^3 m/s m/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
Added: short/3D/PyLith/trunk/tests/3d/hex8/shearelastic-twocells.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/shearelastic-twocells.cfg 2007-06-13 19:03:16 UTC (rev 7209)
+++ short/3D/PyLith/trunk/tests/3d/hex8/shearelastic-twocells.cfg 2007-06-13 19:13:33 UTC (rev 7210)
@@ -0,0 +1,81 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshioascii = 1
+homogeneous = 1
+implicitelasticity = 1
+quadrature3d = 1
+fiatsimplex = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+debug = 1
+importer = pylith.meshio.MeshIOAscii
+
+[pylithapp.mesh_generator.importer]
+filename = twocells.mesh
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.0*s
+default_dt = 1.0*s
+dimension = 3
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCSixSides
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+material = pylith.materials.ElasticIsotropic3D
+
+[pylithapp.timedependent.materials.material]
+label = elastic material
+id = 0
+db.iohandler.filename = matprops-twocells.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1, 2]
+id = 3
+label = x_neg
+db.label = Dirichlet BC -x edge
+db.iohandler.filename = shearxyz-twocells.spatialdb
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1, 2]
+id = 4
+label = x_pos
+db.label = Dirichlet BC +x edge
+db.iohandler.filename = shearxyz-twocells.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+pc_type = jacobi
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+filename = shearelastic-twocells.vtk
Added: short/3D/PyLith/trunk/tests/3d/hex8/shearxyz-twocells.spatialdb
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/shearxyz-twocells.spatialdb 2007-06-13 19:03:16 UTC (rev 7209)
+++ short/3D/PyLith/trunk/tests/3d/hex8/shearxyz-twocells.spatialdb 2007-06-13 19:13:33 UTC (rev 7210)
@@ -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
+ }
+}
++1.0 0.0 0.0 0.0 +0.01 0.0
+-1.0 0.0 0.0 0.0 -0.01 0.0
Added: short/3D/PyLith/trunk/tests/3d/hex8/twocells.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/3d/hex8/twocells.mesh 2007-06-13 19:03:16 UTC (rev 7209)
+++ short/3D/PyLith/trunk/tests/3d/hex8/twocells.mesh 2007-06-13 19:13:33 UTC (rev 7210)
@@ -0,0 +1,50 @@
+mesh = {
+ dimension = 3
+ use-index-zero = true
+ vertices = {
+ dimension = 3
+ count = 12
+ coordinates = {
+ 0 -1.0 -0.5 -0.5
+ 1 -1.0 0.5 -0.5
+ 2 0.0 -0.5 -0.5
+ 3 0.0 0.5 -0.5
+ 4 1.0 -0.5 -0.5
+ 5 1.0 0.5 -0.5
+ 6 -1.0 -0.5 0.5
+ 7 -1.0 0.5 0.5
+ 8 0.0 -0.5 0.5
+ 9 0.0 0.5 0.5
+ 10 1.0 -0.5 0.5
+ 11 1.0 0.5 0.5
+ }
+ }
+ cells = {
+ count = 2
+ num-corners = 8
+ simplices = {
+ 0 0 2 3 1 6 8 9 7
+ 1 2 4 5 3 8 10 11 9
+ }
+ material-ids = {
+ 0 0
+ 1 0
+ }
+ }
+ group = {
+ name = x_neg
+ type = vertices
+ count = 4
+ indices = {
+ 0 1 6 7
+ }
+ }
+ group = {
+ name = x_pos
+ type = vertices
+ count = 4
+ indices = {
+ 4 5 10 11
+ }
+ }
+}
More information about the cig-commits
mailing list