[cig-commits] r17160 - in short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27: config mesh

[willic3 at geodynamics.org]

Author: willic3
Date: 2010-08-31 20:26:33 -0700 (Tue, 31 Aug 2010)
New Revision: 17160

Modified:
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou
   short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/threehex27.exo
Log:
Updated threehex27 tests to work properly, and added nodeset defining
central layer of nodes added for quadratic cells.



Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/elas_ng_ni_v03.cfg	2010-09-01 03:26:33 UTC (rev 17160)
@@ -33,8 +33,11 @@
 [pylithapp.timedependent]
 # Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
 # 'z_pos', and 'z_neg'.
-bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+# bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
 
+# Option for fixing central layer of vertices.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg,z_mid]
+
 [pylithapp.timedependent.formulation]
 # Set the output to an array of 1 output managers.
 # We will output the solution over the domain.
@@ -58,6 +61,10 @@
 db_properties.iohandler.filename = ../spatialdb/mat_elastic.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
+quadrature.cell.degree = 2
+
+# Option for changing the integration order.
+# quadrature.cell.quad_order = 2
 output.cell_info_fields = [density,mu,lambda]
 output.cell_data_fields = [total_strain,stress]
 
@@ -93,6 +100,12 @@
 label = face_zpos
 db_initial.label = Dirichlet BC on +z
 
+# face_zmid face
+[pylithapp.timedependent.bc.z_mid]
+bc_dof = [2]
+label = face_zmid
+db_initial.label = Dirichlet BC on face_zmid
+
 # -y face
 [pylithapp.timedependent.bc.y_neg]
 bc_dof = [0, 1]

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp1_ng_ni_v03.cfg	2010-09-01 03:26:33 UTC (rev 17160)
@@ -37,8 +37,11 @@
 
 # Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
 # 'z_pos', and 'z_neg'.
-bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+# bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
 
+# Option for fixing central layer of vertices.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg,z_mid]
+
 [pylithapp.timedependent.formulation]
 # Set the output to an array of 1 output manager.
 # We will output the solution over the domain.
@@ -63,6 +66,10 @@
 db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp1.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
+quadrature.cell.degree = 2
+
+# Option for changing the integration order.
+# quadrature.cell.quad_order = 2
 output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
 output.cell_data_fields = [total_strain,stress,plastic_strain]
 
@@ -98,6 +105,12 @@
 label = face_zpos
 db_initial.label = Dirichlet BC on +z
 
+# face_zmid face
+[pylithapp.timedependent.bc.z_mid]
+bc_dof = [2]
+label = face_zmid
+db_initial.label = Dirichlet BC on face_zmid
+
 # -y face
 [pylithapp.timedependent.bc.y_neg]
 bc_dof = [0, 1]

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp2_ng_ni_v03.cfg	2010-09-01 03:26:33 UTC (rev 17160)
@@ -37,8 +37,11 @@
 
 # Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
 # 'z_pos', and 'z_neg'.
-bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+# bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
 
+# Option for fixing central layer of vertices.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg,z_mid]
+
 [pylithapp.timedependent.formulation]
 # Set the output to an array of 1 output manager.
 # We will output the solution over the domain.
@@ -63,6 +66,10 @@
 db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp2.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
+quadrature.cell.degree = 2
+
+# Option for changing the integration order.
+# quadrature.cell.quad_order = 2
 output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
 output.cell_data_fields = [total_strain,stress,plastic_strain]
 
@@ -98,6 +105,12 @@
 label = face_zpos
 db_initial.label = Dirichlet BC on +z
 
+# face_zmid face
+[pylithapp.timedependent.bc.z_mid]
+bc_dof = [2]
+label = face_zmid
+db_initial.label = Dirichlet BC on face_zmid
+
 # -y face
 [pylithapp.timedependent.bc.y_neg]
 bc_dof = [0, 1]

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_dp3_ng_ni_v03.cfg	2010-09-01 03:26:33 UTC (rev 17160)
@@ -37,8 +37,11 @@
 
 # Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
 # 'z_pos', and 'z_neg'.
-bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+# bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
 
+# Option for fixing central layer of vertices.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg,z_mid]
+
 [pylithapp.timedependent.formulation]
 # Set the output to an array of 1 output manager.
 # We will output the solution over the domain.
@@ -63,6 +66,10 @@
 db_properties.iohandler.filename = ../spatialdb/mat_plastic_dp3.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
+quadrature.cell.degree = 2
+
+# Option for changing the integration order.
+# quadrature.cell.quad_order = 2
 output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
 output.cell_data_fields = [total_strain,stress,plastic_strain]
 
@@ -98,6 +105,12 @@
 label = face_zpos
 db_initial.label = Dirichlet BC on +z
 
+# face_zmid face
+[pylithapp.timedependent.bc.z_mid]
+bc_dof = [2]
+label = face_zmid
+db_initial.label = Dirichlet BC on face_zmid
+
 # -y face
 [pylithapp.timedependent.bc.y_neg]
 bc_dof = [0, 1]

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/plas_vm1_ng_ni_v03.cfg	2010-09-01 03:26:33 UTC (rev 17160)
@@ -36,8 +36,11 @@
 
 # Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
 # 'z_pos', and 'z_neg'.
-bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+# bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
 
+# Option for fixing central layer of vertices.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg,z_mid]
+
 [pylithapp.timedependent.formulation]
 # Set the output to an array of 1 output manager.
 # We will output the solution over the domain.
@@ -62,6 +65,10 @@
 db_properties.iohandler.filename = ../spatialdb/mat_plastic_vm1.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
+quadrature.cell.degree = 2
+
+# Option for changing the integration order.
+# quadrature.cell.quad_order = 2
 output.cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
 output.cell_data_fields = [total_strain,stress,plastic_strain]
 
@@ -97,6 +104,12 @@
 label = face_zpos
 db_initial.label = Dirichlet BC on +z
 
+# face_zmid face
+[pylithapp.timedependent.bc.z_mid]
+bc_dof = [2]
+label = face_zmid
+db_initial.label = Dirichlet BC on face_zmid
+
 # -y face
 [pylithapp.timedependent.bc.y_neg]
 bc_dof = [0, 1]

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/pylithapp.cfg	2010-09-01 03:26:33 UTC (rev 17160)
@@ -32,7 +32,7 @@
 
 [pylithapp.mesh_generator.reader]
 # Set filename of mesh to import.
-filename = ../mesh/threehex8.exo
+filename = ../mesh/threehex27.exo
 coordsys.space_dim = 3
 
 # ----------------------------------------------------------------------

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/config/visc_e10_ng_ni_v03.cfg	2010-09-01 03:26:33 UTC (rev 17160)
@@ -34,8 +34,11 @@
 [pylithapp.timedependent]
 # Set bc to an array of 5 boundary conditions: 'x_pos','x_neg', 'y_neg',
 # 'z_pos', and 'z_neg'.
-bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
+# bc = [x_pos,x_neg,y_neg,z_pos,z_neg]
 
+# Option for fixing central layer of vertices.
+bc = [x_pos,x_neg,y_neg,z_pos,z_neg,z_mid]
+
 [pylithapp.timedependent.formulation]
 # Set the output to an array of 1 output manager.
 # We will output the solution over the domain.
@@ -65,6 +68,11 @@
 db_properties.iohandler.filename = ../spatialdb/mat_viscous_e10.spatialdb
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 3
+quadrature.cell.degree = 2
+
+# Option for changing the integration order.
+# quadrature.cell.quad_order = 2
+
 output.cell_info_fields = [density,mu,lambda,maxwell_time]
 output.cell_data_fields = [total_strain,stress,viscous_strain]
 
@@ -100,6 +108,12 @@
 label = face_zpos
 db_initial.label = Dirichlet BC on +z
 
+# face_zmid face
+[pylithapp.timedependent.bc.z_mid]
+bc_dof = [2]
+label = face_zmid
+db_initial.label = Dirichlet BC on face_zmid
+
 # -y face
 [pylithapp.timedependent.bc.y_neg]
 bc_dof = [0, 1]

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou	2010-08-31 18:49:31 UTC (rev 17159)
+++ short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/mesh_threehex27.jou	2010-09-01 03:26:33 UTC (rev 17160)
@@ -76,6 +76,13 @@
 nodeset 16 name "face_zpos"
 
 # ----------------------------------------------------------------------
+# Create nodeset for middle-z face
+# ----------------------------------------------------------------------
+group "face_zmid" add node with z_coord > -0.1 and with z_coord < 0.1
+nodeset 17 group face_zmid
+nodeset 17 name "face_zmid"
+
+# ----------------------------------------------------------------------
 # Export exodus file
 # ----------------------------------------------------------------------
 export mesh "threehex27.exo" dimension 3 overwrite

Modified: short/3D/PyLith/branches/v1.5-stable/tests/3d/plasticity/threehex27/mesh/threehex27.exo
===================================================================
(Binary files differ)