[cig-commits] r8028 - short/3D/PyLith/trunk/examples/3d/hex8
brad at geodynamics.org
brad at geodynamics.org
Tue Sep 25 12:32:05 PDT 2007
Author: brad
Date: 2007-09-25 12:32:05 -0700 (Tue, 25 Sep 2007)
New Revision: 8028
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/box_hex8_1000m.exo
short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou
short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou
short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
Log:
Changed dislocation BC to x & y DOF fixed on +x and -x faces and z DOF fixed on -z face. Updated CUBIT journal script to use groups.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/box_hex8_1000m.exo
===================================================================
(Binary files differ)
Modified: short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2007-09-25 18:29:09 UTC (rev 8027)
+++ short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2007-09-25 19:32:05 UTC (rev 8028)
@@ -40,17 +40,23 @@
# The label corresponds to the nodeset ID in CUBIT.
[pylithapp.timedependent.bc.x_pos]
-fixed_dof = [0, 1, 2]
+fixed_dof = [0, 1]
label = 11
db.label = Dirichlet BC on +x
db.iohandler.filename = fixeddisp_zero.spatialdb
[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0, 1, 2]
+fixed_dof = [0, 1]
label = 12
db.label = Dirichlet BC on -x
db.iohandler.filename = fixeddisp_zero.spatialdb
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+db.iohandler.filename = fixeddisp_zero.spatialdb
+
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou 2007-09-25 18:29:09 UTC (rev 8027)
+++ short/3D/PyLith/trunk/examples/3d/hex8/geometry.jou 2007-09-25 19:32:05 UTC (rev 8028)
@@ -43,5 +43,5 @@
# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
-imprint volume all with volume all
-merge volume all
+imprint all with volume all
+merge all
Modified: short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou 2007-09-25 18:29:09 UTC (rev 8027)
+++ short/3D/PyLith/trunk/examples/3d/hex8/mesh_hex8_1000m.jou 2007-09-25 19:32:05 UTC (rev 8028)
@@ -15,7 +15,7 @@
# ----------------------------------------------------------------------
# Set discretization size
# ----------------------------------------------------------------------
-volume all size 2000
+volume all size 1000
# ----------------------------------------------------------------------
# Generate the mesh
@@ -37,40 +37,64 @@
# ----------------------------------------------------------------------
# Create nodeset for fault
# ----------------------------------------------------------------------
-nodeset 10 surface fault_surface fault_surface at A
-nodeset 10 name "fault surface"
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+nodeset 10 group fault
+nodeset 10 name "fault"
# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
-nodeset 11 surface 20 28
+group "face_xpos" add node in surface 20
+group "face_xpos" add node in surface 28
+nodeset 11 group face_xpos
nodeset 11 name "face xpos"
# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
-nodeset 12 surface 30 38
+group "face_xneg" add node in surface 30
+group "face_xneg" add node in surface 38
+nodeset 12 group face_xneg
nodeset 12 name "face xneg"
# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
-nodeset 13 surface 21 27 33 35
+group "face_ypos" add node in surface 21
+group "face_ypos" add node in surface 27
+group "face_ypos" add node in surface 33
+group "face_ypos" add node in surface 35
+nodeset 13 group face_ypos
nodeset 13 name "face ypos"
# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
-nodeset 14 surface 23 25 31 37
+group "face_yneg" add node in surface 23
+group "face_yneg" add node in surface 25
+group "face_yneg" add node in surface 31
+group "face_yneg" add node in surface 37
+nodeset 14 group face_yneg
nodeset 14 name "face yneg"
# ----------------------------------------------------------------------
# Create nodeset for -z face
# ----------------------------------------------------------------------
-nodeset 15 surface 12 16
-nodeset 15 name "face -z"
+group "face_zneg" add node in surface 12
+group "face_zneg" add node in surface 16
+nodeset 15 group face_zneg
+nodeset 15 name "face zneg"
# ----------------------------------------------------------------------
+# Create nodeset for -z face w/o fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 group face_zneg_nofault
+nodeset 16 name "face zneg nofault"
+
+# ----------------------------------------------------------------------
# Export exodus file
# ----------------------------------------------------------------------
export mesh "box_hex8_1000m.exo" dimension 3 overwrite
Modified: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2007-09-25 18:29:09 UTC (rev 8027)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2007-09-25 19:32:05 UTC (rev 8028)
@@ -11,6 +11,7 @@
petsc = 1
solverlinear = 1
meshiocubit = 1
+meshimporter = 1
implicitelasticity = 1
faultcohesivekin = 1
meshimporter = 1
@@ -84,6 +85,7 @@
# Set the solver options.
[pylithapp.petsc]
+ksp_rtol = 1.0e-8
pc_type = bjacobi
ksp_monitor = true
ksp_view = true
More information about the cig-commits
mailing list