[cig-commits] r18924 - short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes

Fri Sep 16 20:44:01 PDT 2011

Author: willic3
Date: 2011-09-16 20:44:00 -0700 (Fri, 16 Sep 2011)
New Revision: 18924

Added:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou
Removed:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.jou
Modified:
   short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou
Log:
Moved refined tet mesh journal file and made some changes to geometry file.



Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou
===================================================================

--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou	2011-09-16 19:13:45 UTC (rev 18923)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/geometry.jou	2011-09-17 03:44:00 UTC (rev 18924)
@@ -1,20 +1,13 @@
-## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
-## Cubit Version 10.2
-## Cubit Build 8
-## Revised 10/17/2006 15:13:41 MST
-## Running 06/16/2008 02:47:47 PM
-## Command Options:
-## -warning = On
-## -information = On
-
 # Description:
 # This journal creates defines the geometry for a mesh with dimensions of
-# 2000 km (x), 1000 km (y), and 400 km (z).  The volume is divided into 3
+# 2000 km (x), 1000 km (y), and 400 km (z).  The volume is divided into 2
 # horizontal layers:
-# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
-# a middle layer from 20 to 40 km, marking the creeping section;
+# a top layer from 0 to 40 km, marking the surface elastic layer;
 # a bottom layer from 40 to 400 km, which approximates the viscoelastic
 # half-space.
+# The fault plane is also cut horizontally at 20 km depth, representing
+# the division between the upper locked section and the lower
+# creeping section.
 # The volume is cut vertically by the fault surface, which lies along the
 # x-axis, and by a symmetry plane, which lies along the y-axis.
 # The fault is restricted to the upper two (elastic) layers.
@@ -43,7 +36,7 @@
 create planar surface with plane yplane offset 0
 surface 9 name "symmetry_plane"
 create planar surface with plane zplane offset -20000
-surface 10 name "fault_bottom"
+surface 10 name "locking_transition"
 
 # ----------------------------------------------------------------------
 # Cut volume with interfaces.
@@ -52,17 +45,17 @@
 webcut volume all with plane fault_surface
 webcut volume all with plane material_interface
 webcut volume all with plane symmetry_plane
-webcut volume all with plane fault_bottom
+webcut volume 2 with plane locking_transition
+webcut volume 16 with plane locking_transition
 
 # ----------------------------------------------------------------------
-# Imprint all volumes, merging surfaces.
+# Remove surface bodies.
 # ----------------------------------------------------------------------
-imprint all with volume all
-merge all
+delete body 3 4 5 7 8 9 11 12 14 17 18 20 21 23
 
 # ----------------------------------------------------------------------
-# Remove surface bodies.
+# Imprint all volumes, merging surfaces.
 # ----------------------------------------------------------------------
-delete body 2 3 4 5 7 8 9 11 12 14 16 17 18 20 21 23 26 27 29 31
 imprint all with volume all
 merge all
+delete volume 2 16 25 26

Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.jou	2011-09-16 19:13:45 UTC (rev 18923)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.jou	2011-09-17 03:44:00 UTC (rev 18924)
@@ -1,156 +0,0 @@
-## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
-## Cubit Version 10.2
-## Cubit Build 8
-## Revised 10/17/2006 15:13:41 MST
-## Running 06/16/2008 02:47:47 PM
-## Command Options:
-## -warning = On
-## -information = On
-
-# Description:
-# This journal uses the geometry from file geometry.jou to create a mesh
-# with dimensions of 2000 km (x), 1000 km (y), and 400 km (z).  The mesh is
-# divided into 3 horizontal layers:
-# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
-# a middle layer from 20 to 40 km, marking the creeping section;
-# a bottom layer from 40 to 400 km, which approximates the viscoelastic
-# half-space.
-# The mesh is cut vertically by the fault surface, which lies along the
-# x-axis, and by a symmetry plane, which lies along the y-axis.
-# The fault is restricted to the upper two (elastic) layers.
-#
-# The volume is meshed at an initial resolution of 20 km.
-# An inner volume with dimensions of 480 km (x), 240 km (y), and 100 km (z)
-# is then put through 1 level of refinement.
-# The final mesh contains 125760 vertices and 705181 cells.
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry.jou'
-
-# ----------------------------------------------------------------------
-# Set initial interval size to 20 km and set volume meshing scheme to
-# tetmesh
-# ----------------------------------------------------------------------
-volume all size 20000.0
-volume all scheme tetmesh
-
-# ----------------------------------------------------------------------
-# Set surface meshing scheme to qtri with a map base scheme.
-# This gives us regularly-spaced vertices that should have the same
-# locations as those for the hex meshes.
-# ----------------------------------------------------------------------
-surface all scheme qtri base scheme map
-mesh surface all
-
-# ----------------------------------------------------------------------
-# Mesh the volume.
-# ----------------------------------------------------------------------
-mesh volume all
-
-# ----------------------------------------------------------------------
-# Refine region inside inner box.
-# ----------------------------------------------------------------------
-group "inner_box" add tet with x_coord > -240000.0 and with x_coord < 240000.0 and with y_coord > -120000.0 and with y_coord < 120000.0 and with z_coord > -100000.0
-refine tet in inner_box depth 0
-
-# ----------------------------------------------------------------------
-# Smooth the mesh.
-# ----------------------------------------------------------------------
-set debug 91 on
-volume all smooth scheme condition number
-smooth volume all
-
-# ----------------------------------------------------------------------
-# Create blocks for materials
-# ----------------------------------------------------------------------
-block 1 volume 1 6 15 19 25 28 30 32
-block 1 name "elastic"
-block 2 volume 10 13 22 24
-block 2 name "viscoelastic"
-
-# ----------------------------------------------------------------------
-# Create fault nodeset
-# ----------------------------------------------------------------------
-group "fault" add node in fault_surface
-group "fault" add node in fault_surface at A
-group "fault" add node in fault_surface at D
-group "fault" add node in fault_surface at E
-nodeset 10 fault
-nodeset 10 name "fault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +x face
-# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 133 109 139 111 89 87
-nodeset 11 face_xpos
-nodeset 11 name "face_xpos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -x face
-# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 121 147 127 149 97 104
-nodeset 12 face_xneg
-nodeset 12 name "face_xneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face
-# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 106 120 113 128 35 47
-nodeset 13 face_ypos
-nodeset 13 name "face_ypos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face
-# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 144 132 151 140 49 33
-nodeset 14 face_yneg
-nodeset 14 name "face_yneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -z face
-# ----------------------------------------------------------------------
-group "face_zneg" add node in surface 103 92 96 84
-nodeset 15 face_zneg
-nodeset 15 name "face_zneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -z face without fault
-# ----------------------------------------------------------------------
-group "face_zneg_nofault" add node in face_zneg
-group "face_zneg_nofault" remove node in fault
-nodeset 16 face_zneg_nofault
-nodeset 16 name "face_zneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +z face
-# ----------------------------------------------------------------------
-group "face_zpos" add node in surface 78 62 70 54
-nodeset 17 face_zpos
-nodeset 17 name "face_zpos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_yneg_nofault" add node in face_yneg
-group "face_yneg_nofault" remove node in fault
-group "face_yneg_nofault" remove node in face_xpos
-group "face_yneg_nofault" remove node in face_xneg
-nodeset 18 face_yneg_nofault
-nodeset 18 name "face_yneg_nofault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault or x-boundaries
-# ----------------------------------------------------------------------
-group "face_ypos_nofault" add node in face_ypos
-group "face_ypos_nofault" remove node in fault
-group "face_ypos_nofault" remove node in face_xpos
-group "face_ypos_nofault" remove node in face_xneg
-nodeset 19 face_ypos_nofault
-nodeset 19 name "face_ypos_nofault"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "spbm_tet4_graded_20km.exo" dimension 3 overwrite
-export Abaqus "spbm_tet4_graded_20km.inp"

Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou (from rev 18923, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_graded_20km.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/spbm_tet4_refine1_20km.jou	2011-09-17 03:44:00 UTC (rev 18924)
@@ -0,0 +1,156 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/16/2008 02:47:47 PM
+## Command Options:
+## -warning = On
+## -information = On
+
+# Description:
+# This journal uses the geometry from file geometry.jou to create a mesh
+# with dimensions of 2000 km (x), 1000 km (y), and 400 km (z).  The mesh is
+# divided into 3 horizontal layers:
+# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
+# a middle layer from 20 to 40 km, marking the creeping section;
+# a bottom layer from 40 to 400 km, which approximates the viscoelastic
+# half-space.
+# The mesh is cut vertically by the fault surface, which lies along the
+# x-axis, and by a symmetry plane, which lies along the y-axis.
+# The fault is restricted to the upper two (elastic) layers.
+#
+# The volume is meshed at an initial resolution of 20 km.
+# An inner volume with dimensions of 480 km (x), 240 km (y), and 100 km (z)
+# is then put through 1 level of refinement.
+# The final mesh contains 125760 vertices and 705181 cells.
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set initial interval size to 20 km and set volume meshing scheme to
+# tetmesh
+# ----------------------------------------------------------------------
+volume all size 20000.0
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Set surface meshing scheme to qtri with a map base scheme.
+# This gives us regularly-spaced vertices that should have the same
+# locations as those for the hex meshes.
+# ----------------------------------------------------------------------
+surface all scheme qtri base scheme map
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Mesh the volume.
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine region inside inner box.
+# ----------------------------------------------------------------------
+group "inner_box" add tet with x_coord > -240000.0 and with x_coord < 240000.0 and with y_coord > -120000.0 and with y_coord < 120000.0 and with z_coord > -100000.0
+refine tet in inner_box depth 0
+
+# ----------------------------------------------------------------------
+# Smooth the mesh.
+# ----------------------------------------------------------------------
+set debug 91 on
+volume all smooth scheme condition number
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 6 15 19 25 28 30 32
+block 1 name "elastic"
+block 2 volume 10 13 22 24
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface at D
+group "fault" add node in fault_surface at E
+nodeset 10 fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 133 109 139 111 89 87
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 121 147 127 149 97 104
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 106 120 113 128 35 47
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 144 132 151 140 49 33
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 103 92 96 84
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 78 62 70 54
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_yneg_nofault" add node in face_yneg
+group "face_yneg_nofault" remove node in fault
+group "face_yneg_nofault" remove node in face_xpos
+group "face_yneg_nofault" remove node in face_xneg
+nodeset 18 face_yneg_nofault
+nodeset 18 name "face_yneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_ypos_nofault" add node in face_ypos
+group "face_ypos_nofault" remove node in fault
+group "face_ypos_nofault" remove node in face_xpos
+group "face_ypos_nofault" remove node in face_xneg
+nodeset 19 face_ypos_nofault
+nodeset 19 name "face_ypos_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_tet4_graded_20km.exo" dimension 3 overwrite
+export Abaqus "spbm_tet4_graded_20km.inp"

