[cig-commits] r16523 - short/3D/PyLith/trunk/tests/2d/quad4
surendra at geodynamics.org
surendra at geodynamics.org
Fri Apr 9 14:12:27 PDT 2010
Author: surendra
Date: 2010-04-09 14:12:26 -0700 (Fri, 09 Apr 2010)
New Revision: 16523
Modified:
short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py
short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py
short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg
short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding_soln.py
short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg
short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick_soln.py
short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou
short/3D/PyLith/trunk/tests/2d/quad4/testpylith.py
Log:
Finished implementing tests for .vtk of friction full-scale examples
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearSliding.py 2010-04-09 21:12:26 UTC (rev 16523)
@@ -127,8 +127,17 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
- initialTraction = (0.0, -1.0e+10)
+ initialTraction = (0.0, -1.0e+6)
+ uy_l = 1.0
+ len = 8000.0
+ sigma_f = 0.6e+6
+ p_density = 2500.0
+ p_vs = 3000.0
+ p_mu = p_density*p_vs**2
+ D = uy_l - len * sigma_f/p_mu
+ slip = (D, 0.0)
+
nvertices = self.faultMesh['nvertices']
if name == "strike_dir":
@@ -148,11 +157,13 @@
elif name == "slip":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
+ field[:,0] = slip[0]
+ field[:,1] = slip[1]
elif name == "traction":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
- field[:,0] = 0.0
- field[:,1] = -2.2e+10
+ field[:,0] = 0.6e+6
+ field[:,1] = -1.0e+6
else:
raise ValueError("Unknown fault field '%s'." % name)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/TestFrictionShearStick.py 2010-04-09 21:12:26 UTC (rev 16523)
@@ -127,7 +127,7 @@
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
- initialTraction = (0.0, -1.0e+10)
+ initialTraction = (0.0, -1.0e+7)
nvertices = self.faultMesh['nvertices']
@@ -151,8 +151,8 @@
elif name == "traction":
field = numpy.zeros( (nvertices, 3), dtype=numpy.float64)
- field[:,0] = 0.0
- field[:,1] = -2.2e+10
+ field[:,0] = 1.0e+6
+ field[:,1] = -1.0e+7
else:
raise ValueError("Unknown fault field '%s'." % name)
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding.cfg 2010-04-09 21:12:26 UTC (rev 16523)
@@ -121,7 +121,7 @@
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Dirichlet BC +x edge
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*m,4000.0*m]
+db_initial.data = [0.0*m,1.0*m]
@@ -130,12 +130,12 @@
# The nodes associated with this boundary condition have the name
# 'y_pos_fault' in the mesh file.
-label = 25
+label = 22
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Neumann BC +y edge
db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [-6000*MPa, 0.0*MPa]
+db_initial.data = [-0.6*MPa, 0.0*MPa]
# Set cell type to quadrilateral (1-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
@@ -146,12 +146,12 @@
# The nodes associated with this boundary condition have the name
# 'y_neg_fault' in the mesh file.
-label = 24
+label = 23
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Neumann BC -y edge
db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [-6000*MPa, 0.0*MPa]
+db_initial.data = [-0.6*MPa, 0.0*MPa]
# Set cell type to quadrilateral (1-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
@@ -189,7 +189,7 @@
db_initial_tractions = spatialdata.spatialdb.UniformDB
db_initial_tractions.label = "Initial fault tractions"
db_initial_tractions.values = [traction-shear,traction-normal]
-db_initial_tractions.data = [0.0*Pa, -10000.0*MPa]
+db_initial_tractions.data = [0.0*Pa, -1.0*MPa]
friction.db_properties = spatialdata.spatialdb.UniformDB
friction.db_properties.label = Static friction
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding_soln.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding_soln.py 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_sliding_soln.py 2010-04-09 21:12:26 UTC (rev 16523)
@@ -25,8 +25,8 @@
p_lambda = p_density*p_vp**2 - 2*p_mu
# Uniform stress field (plane strain)
-sxx = -1.2e+10
-sxy = 0.0
+sxx = 0.0
+sxy = 0.6e+6
syy = 0.0
szz = p_lambda/(2*p_lambda+2*p_mu)*(sxx+syy)
@@ -37,9 +37,15 @@
exy = 1.0/(2*p_mu) * (sxy)
-print exx,eyy,exy,ezz,szz
-print -exx*p_lambda/(p_lambda+2*p_mu)
+# Calculate slip
+uy_l = 1.0
+len = 8000.0
+sigma_f = 0.6e+6
+D = uy_l - len * sigma_f/p_mu
+#print exx,eyy,exy,ezz,szz
+#print -exx*p_lambda/(p_lambda+2*p_mu)
+
# ----------------------------------------------------------------------
class AnalyticalSoln(object):
"""
@@ -57,8 +63,8 @@
(nlocs, dim) = locs.shape
disp = numpy.zeros( (nlocs, 3), dtype=numpy.float64)
- disp[:,0] = exx*(locs[:,0]+max(abs(locs[:,0])))
- disp[:,1] = eyy*(locs[:,1]+max(abs(locs[:,1])))
+ disp[0:nlocs/2,1] = 2 * exy * (locs[0:nlocs/2,0] + len/2) + D
+ disp[nlocs/2:nlocs,1] = 2 * exy * (locs[nlocs/2:nlocs,0] + len/2)
return disp
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick.cfg 2010-04-09 21:12:26 UTC (rev 16523)
@@ -127,7 +127,7 @@
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Neumann BC +x edge
db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [1000.0*MPa,0.0*MPa]
+db_initial.data = [1.0*MPa,0.0*MPa]
# Set cell type to quadrilateral (1-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
@@ -141,12 +141,12 @@
# The nodes associated with this boundary condition have the name
# 'y_pos_fault' in the mesh file.
-label = 25
+label = 22
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Neumann BC +y edge
db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [-1000.0*MPa, 0.0*MPa]
+db_initial.data = [-1.0*MPa, 0.0*MPa]
# Set cell type to quadrilateral (1-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
@@ -157,12 +157,12 @@
# The nodes associated with this boundary condition have the name
# 'y_neg_fault' in the mesh file.
-label = 24
+label = 23
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Neumann BC -y edge
db_initial.values = [traction-shear,traction-normal]
-db_initial.data = [-1000.0*MPa, 0.0*MPa]
+db_initial.data = [-1.0*MPa, 0.0*MPa]
# Set cell type to quadrilateral (1-d Lagrange).
quadrature.cell = pylith.feassemble.FIATLagrange
@@ -202,7 +202,7 @@
db_initial_tractions = spatialdata.spatialdb.UniformDB
db_initial_tractions.label = "Initial fault tractions"
db_initial_tractions.values = [traction-shear,traction-normal]
-db_initial_tractions.data = [0.0*Pa, -10000.0*MPa]
+db_initial_tractions.data = [0.0*Pa, -10.0*MPa]
friction.db_properties = spatialdata.spatialdb.UniformDB
friction.db_properties.label = Static friction
Modified: short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick_soln.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick_soln.py 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/friction_shear_stick_soln.py 2010-04-09 21:12:26 UTC (rev 16523)
@@ -25,8 +25,8 @@
p_lambda = p_density*p_vp**2 - 2*p_mu
# Uniform stress field (plane strain)
-sxx = -1.2e+10
-sxy = 0.0
+sxx = 0.0
+sxy = 1.0e+6
syy = 0.0
szz = p_lambda/(2*p_lambda+2*p_mu)*(sxx+syy)
@@ -37,8 +37,8 @@
exy = 1.0/(2*p_mu) * (sxy)
-print exx,eyy,exy,ezz,szz
-print -exx*p_lambda/(p_lambda+2*p_mu)
+#print exx,eyy,exy,ezz,szz
+#print -exx*p_lambda/(p_lambda+2*p_mu)
# ----------------------------------------------------------------------
class AnalyticalSoln(object):
@@ -57,8 +57,7 @@
(nlocs, dim) = locs.shape
disp = numpy.zeros( (nlocs, 3), dtype=numpy.float64)
- disp[:,0] = exx*(locs[:,0]+max(abs(locs[:,0])))
- disp[:,1] = eyy*(locs[:,1]+max(abs(locs[:,1])))
+ disp[:,1] = 2*exy*(locs[:,0]+max(abs(locs[:,0])))
return disp
Modified: short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/mesh.jou 2010-04-09 21:12:26 UTC (rev 16523)
@@ -72,14 +72,6 @@
nodeset 24 name "edge yneg nofault"
# ----------------------------------------------------------------------
-# Create nodeset for +y edge w/o overlap with fault
-# ----------------------------------------------------------------------
-group "edge_ypos_nofault" add node in edge_ypos
-group "edge_ypos_nofault" remove node in fault_x
-nodeset 25 group edge_ypos_nofault
-nodeset 25 name "edge ypos nofault"
-
-# ----------------------------------------------------------------------
# Export exodus file
# ----------------------------------------------------------------------
export mesh "mesh.exo" dimension 2 overwrite
Modified: short/3D/PyLith/trunk/tests/2d/quad4/testpylith.py
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/testpylith.py 2010-04-09 21:01:40 UTC (rev 16522)
+++ short/3D/PyLith/trunk/tests/2d/quad4/testpylith.py 2010-04-09 21:12:26 UTC (rev 16523)
@@ -20,32 +20,32 @@
"""
suite = unittest.TestSuite()
- #from TestAxialDisp import TestAxialDisp
- #suite.addTest(unittest.makeSuite(TestAxialDisp))
+ from TestAxialDisp import TestAxialDisp
+ suite.addTest(unittest.makeSuite(TestAxialDisp))
- #from TestShearDisp import TestShearDisp
- #suite.addTest(unittest.makeSuite(TestShearDisp))
+ from TestShearDisp import TestShearDisp
+ suite.addTest(unittest.makeSuite(TestShearDisp))
- #from TestDislocation import TestDislocation
- #suite.addTest(unittest.makeSuite(TestDislocation))
+ from TestDislocation import TestDislocation
+ suite.addTest(unittest.makeSuite(TestDislocation))
#from TestDislocation2 import TestDislocation2
#suite.addTest(unittest.makeSuite(TestDislocation2))
- #from TestLgDeformRigidBody import TestRigidBody
- #3suite.addTest(unittest.makeSuite(TestRigidBody))
+ from TestLgDeformRigidBody import TestRigidBody
+ suite.addTest(unittest.makeSuite(TestRigidBody))
- #from TestLgDeformTraction import TestTraction
- #suite.addTest(unittest.makeSuite(TestTraction))
+ from TestLgDeformTraction import TestTraction
+ suite.addTest(unittest.makeSuite(TestTraction))
- #from TestFrictionCompression import TestFrictionCompression
- #suite.addTest(unittest.makeSuite(TestFrictionCompression))
+ from TestFrictionCompression import TestFrictionCompression
+ suite.addTest(unittest.makeSuite(TestFrictionCompression))
- #from TestFrictionShearStick import TestFrictionShearStick
- #suite.addTest(unittest.makeSuite(TestFrictionShearStick))
+ from TestFrictionShearStick import TestFrictionShearStick
+ suite.addTest(unittest.makeSuite(TestFrictionShearStick))
- #from TestFrictionShearSliding import TestFrictionShearSliding
- #suite.addTest(unittest.makeSuite(TestFrictionShearSliding))
+ from TestFrictionShearSliding import TestFrictionShearSliding
+ suite.addTest(unittest.makeSuite(TestFrictionShearSliding))
from TestFrictionOpening import TestFrictionOpening
suite.addTest(unittest.makeSuite(TestFrictionOpening))
More information about the CIG-COMMITS
mailing list