[cig-commits] r12684 - in short/3D/PyLith/trunk/examples/3d: hex8 tet4
brad at geodynamics.org
brad at geodynamics.org
Tue Aug 19 16:37:19 PDT 2008
Author: brad
Date: 2008-08-19 16:37:18 -0700 (Tue, 19 Aug 2008)
New Revision: 12684
Added:
short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/3d/tet4/dislocation.cfg
short/3D/PyLith/trunk/examples/3d/tet4/mesh_tet4_1000m.lagrit
short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg
short/3D/PyLith/trunk/examples/3d/tet4/shearxy.cfg
short/3D/PyLith/trunk/examples/3d/tet4/tet4_1000m_ascii.gmv
short/3D/PyLith/trunk/examples/3d/tet4/tet4_1000m_binary.gmv
Log:
Merged materials on each side of fault and added viscoelastic material for bottom half in 3d/tet4. Adjusted solver and debugging settings to be the same for 3d/tet4 and 3d/hex8.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2008-08-19 22:51:24 UTC (rev 12683)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2008-08-19 23:37:18 UTC (rev 12684)
@@ -11,12 +11,10 @@
petsc = 1
solverlinear = 1
meshiocubit = 1
-meshimporter = 1
implicitelasticity = 1
faultcohesivekin = 1
-meshimporter = 1
+fiatlagrange = 1
#quadrature3d = 1
-#fiatlagrange = 1
# ----------------------------------------------------------------------
# mesh_generator
Modified: short/3D/PyLith/trunk/examples/3d/tet4/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/dislocation.cfg 2008-08-19 22:51:24 UTC (rev 12683)
+++ short/3D/PyLith/trunk/examples/3d/tet4/dislocation.cfg 2008-08-19 23:37:18 UTC (rev 12684)
@@ -93,18 +93,10 @@
filename = dislocation-fault.vtk
# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic_neg1.output]
+[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = dislocation-statevars-neg1.vtk
+writer.filename = dislocation-statevars-elastic.vtk
-[pylithapp.timedependent.materials.elastic_pos1.output]
+[pylithapp.timedependent.materials.viscoelastic.output]
cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = dislocation-statevars-pos1.vtk
-
-[pylithapp.timedependent.materials.elastic_neg2.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = dislocation-statevars-pos2.vtk
-
-[pylithapp.timedependent.materials.elastic_pos2.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = dislocation-statevars-neg2.vtk
+writer.filename = dislocation-statevars-viscoelastic.vtk
Added: short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb 2008-08-19 23:37:18 UTC (rev 12684)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the distribution of material
+// properties for a Maxwell viscoelastic material. In this case,
+// the material properties are uniform.
+//
+// The viscosity corresponds to a Maxwell time of 100 years.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4 // number of material property values
+ value-names = density vs vp viscosity // names of the material property values
+ value-units = kg/m^3 m/s m/s Pa-s // units (conversions not supported)
+ num-locs = 1 // number of locations
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+// (7) viscosity (Pa-s)
+0.0 0.0 0.0 2500.0 3000.0 5291.502622129181 7.10046e+19
Modified: short/3D/PyLith/trunk/examples/3d/tet4/mesh_tet4_1000m.lagrit
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/mesh_tet4_1000m.lagrit 2008-08-19 22:51:24 UTC (rev 12683)
+++ short/3D/PyLith/trunk/examples/3d/tet4/mesh_tet4_1000m.lagrit 2008-08-19 23:37:18 UTC (rev 12684)
@@ -193,6 +193,26 @@
pset / boundary_zp / surface / surf_zp /
*
* ----------------------------------------------------------------------
+* MERGE MATERIALS
+* ----------------------------------------------------------------------
+eltset / elastic1 / itetclr / eq / 1
+eltset / elastic2 / itetclr / eq / 2
+eltset / elastic / union / elastic1 elastic2
+cmo / setatt / box / itetclr / eltset,get,elastic / 1
+eltset / elastic1 / delete
+eltset / elastic2 / delete
+eltset / elastic / delete
+
+eltset / visco1 / itetclr / eq / 3
+eltset / visco2 / itetclr / eq / 4
+eltset / visco / union / visco1 visco2
+cmo / setatt / box / itetclr / eltset,get,visco / 2
+eltset / visco1 / delete
+eltset / visco2 / delete
+eltset / visco / delete
+cmo / printatt / box / itetclr / minmax
+
+* ----------------------------------------------------------------------
* WRITE MESH
* ----------------------------------------------------------------------
*
Modified: short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg 2008-08-19 22:51:24 UTC (rev 12683)
+++ short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg 2008-08-19 23:37:18 UTC (rev 12684)
@@ -15,6 +15,7 @@
implicitelasticity = 1
faultcohesivekin = 1
fiatsimplex = 1
+#quadrature3d = 1
# ----------------------------------------------------------------------
# mesh_generator
@@ -52,57 +53,42 @@
# materials
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-# Define the total time for the simulation and the default time step size.
-total_time = 0.0*s ; total time of simulation
-default_dt = 1.0*s ; time step
-# Set materials to an array of 4 elastic isotropic 3-D (default) materials.
-materials = elastic_neg1,elastic_neg2,elastic_pos1,elastic_pos2
+# Set materials to an array of 2 materials:
+# 'elastic' and 'viscoelastic'
+materials = [elastic,viscoelastic]
-# The lines below set the parameters for each material in the
-# materials array.
+# Change from default material (elastic, isotropic, 3-D material) to
+# Maxwell viscoelastic, isotropic 3-D material.
+materials.viscoelastic = pylith.materials.MaxwellIsotropic3D
-[pylithapp.timedependent.materials.elastic_pos1]
-label = Elastic material +x (upper)
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
id = 1
db.iohandler.filename = mat_elastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
-[pylithapp.timedependent.materials.elastic_neg1]
-label = Elastic material -x (upper)
+[pylithapp.timedependent.materials.viscoelastic]
+label = Viscoelastic material
id = 2
-db.iohandler.filename = mat_elastic.spatialdb
+db.iohandler.filename = mat_viscoelastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature3D
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
-[pylithapp.timedependent.materials.elastic_pos2]
-label = Elastic material +x (lower)
-id = 3
-db.iohandler.filename = mat_elastic.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = tetrahedron
-
-[pylithapp.timedependent.materials.elastic_neg2]
-label = Elastic material -x (lower)
-id = 4
-db.iohandler.filename = mat_elastic.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = tetrahedron
-
# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
# Set the solver options.
[pylithapp.petsc]
-pc_type = bjacobi
+ksp_rtol = 1.0e-8
+pc_type = asm
ksp_monitor = true
ksp_view = true
-#log_summary = true
+log_summary = true
ksp_max_it = 100
ksp_gmres_restart = 50
+# start_in_debugger = true
Modified: short/3D/PyLith/trunk/examples/3d/tet4/shearxy.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/shearxy.cfg 2008-08-19 22:51:24 UTC (rev 12683)
+++ short/3D/PyLith/trunk/examples/3d/tet4/shearxy.cfg 2008-08-19 23:37:18 UTC (rev 12684)
@@ -74,18 +74,10 @@
writer.filename = shearxy-groundsurf.vtk
# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic_neg1.output]
+[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = shearxy-statevars-neg1.vtk
+writer.filename = shearxy-statevars-elastic.vtk
-[pylithapp.timedependent.materials.elastic_pos1.output]
+[pylithapp.timedependent.materials.viscoelastic.output]
cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = shearxy-statevars-pos1.vtk
-
-[pylithapp.timedependent.materials.elastic_neg2.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = shearxy-statevars-pos2.vtk
-
-[pylithapp.timedependent.materials.elastic_pos2.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = shearxy-statevars-neg2.vtk
+writer.filename = shearxy-statevars-viscoelastic.vtk
Modified: short/3D/PyLith/trunk/examples/3d/tet4/tet4_1000m_ascii.gmv
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/tet4_1000m_ascii.gmv 2008-08-19 22:51:24 UTC (rev 12683)
+++ short/3D/PyLith/trunk/examples/3d/tet4/tet4_1000m_ascii.gmv 2008-08-19 23:37:18 UTC (rev 12684)
@@ -1,6 +1,6 @@
gmvinput ascii
codename LaGriT
-simdate 06/22/07
+simdate /19/08
nodes 245
0.00000E+000 0.00000E+000 0.00000E+000 1.00000E+003 1.00000E+003 1.00000E+003 2.00000E+003 2.00000E+003 2.00000E+003 3.00000E+003
3.00000E+003 3.00000E+003 0.00000E+000 0.00000E+000 0.00000E+000 1.00000E+003 1.00000E+003 1.00000E+003 2.00000E+003 2.00000E+003
@@ -1057,55 +1057,53 @@
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4
4 4 4 4 4
endflag
-material 5 0
+material 3 0
mr_elast
mr_elast
-mr_visco
-mr_visco
ERROR
- 3 4 1 2 1 4 2 2 2 4 4 3 1 4 1 2 2 1 1 4
- 3 3 3 2 4 3 2 2 2 2 2 2 4 3 2 3 2 1 4 4
- 3 4 3 1 2 3 2 3 3 4 4 3 2 3 2 2 2 2 2 2
- 3 4 1 4 1 4 4 1 4 3 2 3 2 3 1 1 2 4 2 1
- 2 2 2 3 2 2 3 3 1 2 2 1 2 1 2 1 1 4 2 4
- 1 1 2 2 2 2 4 3 2 1 4 4 3 1 4 4 4 1 2 2
- 2 4 1 4 2 4 2 4 3 3 3 3 4 4 4 1 3 4 4 3
- 4 4 4 4 4 4 3 3 3 4 3 1 3 3 1 3 2 2 2 1
- 1 3 3 1 4 1 4 3 3 1 3 1 3 1 3 1 2 2 2 4
- 4 2 2 2 2 2 4 2 2 2 4 4 2 2 4 4 4 4 4 1
- 3 3 3 4 3 4 3 3 1 3 3 1 1 4 1 4 1 4 2 4
- 4 3 4 4 1 3 1 4 2 4 1 1 1 1 3 1 1 4 2 1
- 1 1 1 2 1 1 1 1 2 1 1 4 4 2 2 2 4 4 4 3
- 4 2 4 2 4 4 3 1 2 3 3 3 3 3 3 3 3 2 3 4
- 1 4 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 2 4 4
- 1 1 3 3 2 2 1 1 2 1 1 2 2 1 1 4 3 3 3 1
- 1 2 2 2 2 2 2 1 4 2 3 2 2 4 1 2 3 4 2 2
- 2 4 1 2 2 1 4 1 3 2 4 3 3 3 3 1 3 3 3 3
- 3 3 1 3 3 3 3 1 1 1 4 1 4 4 4 1 3 2 4 4
- 4 2 4 2 4 2 2 2 2 1 1 3 3 3 1 3 1 2 1 1
- 1 1 1 1 1 1 1 1 2 4 1 4 1 2 2 2 4 2 2 1
- 2 2 2 2 2 1 4 2 1 1 4 3 3 3 3 3 3 1 4 3
- 1 3 3 2 3 3 4 4 4 4 4 4 4 3 4 3 3 4 4 3
- 4 3 4 3 3 4 3 1 1 1 1 1 1 4 1 1 1 3 2 1
- 3 4 4 4 3 2 2 2 1 2 2 2 2 2 2 4 2 4 2 2
- 2 2 2 2 2 2 2 1 3 1 1 3 1 1 4 4 3 3 3 3
- 1 3 1 3 4 4 3 3 4 4 4 2 1 4 4 4 4 1 2 1
- 3 3 2 1 2 1 1 1 1 1 3 1 1 1 1 4 4 2 3 3
- 4 4 4 3 1 1 1 4 1 2 1 4 1 4 3 4 2 2 3 4
- 2 2 2 3 4 4 4 2 4 4 3 2 3 4 4 3 3 4 4 3
- 4 4 3 3 2 1 3 2 4 4 4 4 4 4 4 4 4 2 4 2
- 2 2 3 3 2 2 2 1 1 2 1 1 3 3 1 3 3 3 1 1
- 1 1 1 1 1 1 2 2 2 1 2 2 4 2 2 4 2 2 1 1
- 2 2 3 1 3 3 4 3 3 3 3 3 1 3 3 3 1 3 3 1
- 1 3 1 4 4 4 4 4 4 4 4 4 4 4 1 1 4 4 4 3
- 3 1 1 1 3 3 1 1 4 4 4 4 3 1 4 3 3 2 2 2
- 4 2 2 1 4 2 2 2 2 4 4 2 2 2 4 2 3 3 3 1
- 3 3 3 3 2 2 3 2 2 1 3 3 3 3 3 3 2 2 2 4
- 4 4 2 4 4 3 2 2 4 4 1 1 1 1 1 1 1 1 1 1
- 2 1 2 1 2 2 1 3 3 1 1 3 1 1 1 1 1 2 2 1
- 2 1 2 2 4 2 2 2 2 2 2 1 2 1 2 3 3 3 3 3
- 3 3 3 3 3 3 1 1 3 3 1 1 1 1 1 4 4 4 4 4
- 1 4 4 4 1 1 1 1 1 1 1 2
+ 2 2 1 1 1 2 1 1 1 2 2 2 1 2 1 1 1 1 1 2
+ 2 2 2 1 2 2 1 1 1 1 1 1 2 2 1 2 1 1 2 2
+ 2 2 2 1 1 2 1 2 2 2 2 2 1 2 1 1 1 1 1 1
+ 2 2 1 2 1 2 2 1 2 2 1 2 1 2 1 1 1 2 1 1
+ 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 1 1 2 1 2
+ 1 1 1 1 1 1 2 2 1 1 2 2 2 1 2 2 2 1 1 1
+ 1 2 1 2 1 2 1 2 2 2 2 2 2 2 2 1 2 2 2 2
+ 2 2 2 2 2 2 2 2 2 2 2 1 2 2 1 2 1 1 1 1
+ 1 2 2 1 2 1 2 2 2 1 2 1 2 1 2 1 1 1 1 2
+ 2 1 1 1 1 1 2 1 1 1 2 2 1 1 2 2 2 2 2 1
+ 2 2 2 2 2 2 2 2 1 2 2 1 1 2 1 2 1 2 1 2
+ 2 2 2 2 1 2 1 2 1 2 1 1 1 1 2 1 1 2 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 2 2 2 2
+ 2 1 2 1 2 2 2 1 1 2 2 2 2 2 2 2 2 1 2 2
+ 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 2 2
+ 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 1
+ 1 1 1 1 1 1 1 1 2 1 2 1 1 2 1 1 2 2 1 1
+ 1 2 1 1 1 1 2 1 2 1 2 2 2 2 2 1 2 2 2 2
+ 2 2 1 2 2 2 2 1 1 1 2 1 2 2 2 1 2 1 2 2
+ 2 1 2 1 2 1 1 1 1 1 1 2 2 2 1 2 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 2 1 1 1
+ 1 1 1 1 1 1 2 1 1 1 2 2 2 2 2 2 2 1 2 2
+ 1 2 2 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
+ 2 2 2 2 2 2 2 1 1 1 1 1 1 2 1 1 1 2 1 1
+ 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1
+ 1 1 1 1 1 1 1 1 2 1 1 2 1 1 2 2 2 2 2 2
+ 1 2 1 2 2 2 2 2 2 2 2 1 1 2 2 2 2 1 1 1
+ 2 2 1 1 1 1 1 1 1 1 2 1 1 1 1 2 2 1 2 2
+ 2 2 2 2 1 1 1 2 1 1 1 2 1 2 2 2 1 1 2 2
+ 1 1 1 2 2 2 2 1 2 2 2 1 2 2 2 2 2 2 2 2
+ 2 2 2 2 1 1 2 1 2 2 2 2 2 2 2 2 2 1 2 1
+ 1 1 2 2 1 1 1 1 1 1 1 1 2 2 1 2 2 2 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1
+ 1 1 2 1 2 2 2 2 2 2 2 2 1 2 2 2 1 2 2 1
+ 1 2 1 2 2 2 2 2 2 2 2 2 2 2 1 1 2 2 2 2
+ 2 1 1 1 2 2 1 1 2 2 2 2 2 1 2 2 2 1 1 1
+ 2 1 1 1 2 1 1 1 1 2 2 1 1 1 2 1 2 2 2 1
+ 2 2 2 2 1 1 2 1 1 1 2 2 2 2 2 2 1 1 1 2
+ 2 2 1 2 2 2 1 1 2 2 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 2 2 1 1 2 1 1 1 1 1 1 1 1
+ 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2
+ 2 2 2 2 2 2 1 1 2 2 1 1 1 1 1 2 2 2 2 2
+ 1 2 2 2 1 1 1 1 1 1 1 1
cycleno 0
probtime 0.00000E+000
endgmv
Modified: short/3D/PyLith/trunk/examples/3d/tet4/tet4_1000m_binary.gmv
===================================================================
(Binary files differ)
More information about the cig-commits
mailing list