def setUp(self):
"""
Setup for test.
"""
TestQuad4.setUp(self)
self.mesh['nvertices'] = 81 + 9
self.nverticesO = 81
self.faultMesh = {
'nvertices': 9,
'spaceDim': 2,
'ncells': 8,
'ncorners': 2
}
run_pylith(SWShearSlidingApp)
self.outputRoot = "slipweakening_shear_sliding"
self.soln = AnalyticalSoln()
return
def setUp(self):
"""
Setup for test.
"""
TestQuad4.setUp(self)
self.mesh['nvertices'] = 81+9
self.nverticesO = 81
self.faultMesh = {'nvertices': 9,
'spaceDim': 2,
'ncells': 8,
'ncorners': 2}
run_pylith()
self.outputRoot = "slipweakening_shear_sliding"
self.soln = AnalyticalSoln()
return
class TestSlipWeakeningShearSliding(TestQuad4):
"""
Test suite for testing pylith with 2-D shear sliding with slipweakening.
"""
def setUp(self):
"""
Setup for test.
"""
TestQuad4.setUp(self)
self.mesh['nvertices'] = 81+9
self.nverticesO = 81
self.faultMesh = {'nvertices': 9,
'spaceDim': 2,
'ncells': 8,
'ncorners': 2}
run_pylith()
self.outputRoot = "slipweakening_shear_sliding"
self.soln = AnalyticalSoln()
return
def test_fault_info(self):
"""
Check fault information.
"""
if not self.checkResults:
return
filename = "%s-fault_info.h5" % self.outputRoot
fields = ["strike_dir", "normal_dir", "traction_initial","static_coefficient","dynamic_coefficient","slip_weakening_parameter","cohesion"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
def test_fault_data(self):
"""
Check fault information.
"""
if not self.checkResults:
return
filename = "%s-fault.h5" % self.outputRoot
fields = ["slip", "traction","cumulative_slip","previous_slip"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
def calcDisplacements(self, vertices):
"""
Calculate displacement field given coordinates of vertices.
"""
return self.soln.displacement(vertices, self.nverticesO)
def calcStateVar(self, name, vertices, cells):
"""
Calculate state variable.
"""
ncells = self.mesh['ncells']
pts = numpy.zeros( (ncells, 3), dtype=numpy.float64)
if name == "total_strain":
stateVar = self.soln.strain(pts)
elif name == "stress":
stateVar = self.soln.stress(pts)
else:
raise ValueError("Unknown state variable '%s'." % name)
return stateVar
def calcFaultField(self, name, vertices):
"""
Calculate fault info.
"""
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
staticCoefficient = 0.61
dynamicCoefficient = 0.60
slipWeakeningParameter = 0.2
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)
cumulativeSlip = D
nvertices = self.faultMesh['nvertices']
if name == "strike_dir":
field = numpy.zeros( (1, nvertices, 2), dtype=numpy.float64)
field[0,:,0] = strikeDir[0]
field[0,:,1] = strikeDir[1]
elif name == "normal_dir":
field = numpy.zeros( (1, nvertices, 2), dtype=numpy.float64)
field[0,:,0] = normalDir[0]
field[0,:,1] = normalDir[1]
elif name == "traction_initial":
field = numpy.zeros( (1, nvertices, 2), dtype=numpy.float64)
field[0,:,0] = initialTraction[0]
field[0,:,1] = initialTraction[1]
elif name == "static_coefficient":
field = numpy.zeros( (1, nvertices, 1), dtype=numpy.float64)
field[0,:,0] = staticCoefficient
elif name == "dynamic_coefficient":
field = numpy.zeros( (1, nvertices, 1), dtype=numpy.float64)
field[0,:,0] = dynamicCoefficient
elif name == "slip_weakening_parameter":
field = numpy.zeros( (1, nvertices, 1), dtype=numpy.float64)
field[0,:,0] = slipWeakeningParameter
elif name == "cohesion":
field = numpy.zeros( (1, nvertices, 1), dtype=numpy.float64)
elif name == "slip":
field = numpy.zeros( (1, nvertices, 2), dtype=numpy.float64)
field[0,:,0] = slip[0]
field[0,:,1] = slip[1]
elif name == "traction":
field = numpy.zeros( (1, nvertices, 2), dtype=numpy.float64)
field[0,:,0] = 0.6e+6
field[0,:,1] = -1.0e+6
elif name == "cumulative_slip":
field = numpy.zeros( (1, nvertices, 1), dtype=numpy.float64)
field[0,:,0] = cumulativeSlip
elif name == "previous_slip":
field = numpy.zeros( (1, nvertices, 1), dtype=numpy.float64)
field[0,:,0] = cumulativeSlip
else:
raise ValueError("Unknown fault field '%s'." % name)
return field
class TestSlipWeakeningShearSliding(TestQuad4):
"""
Test suite for testing pylith with 2-D shear sliding with slipweakening.
"""
def setUp(self):
"""
Setup for test.
"""
TestQuad4.setUp(self)
self.mesh['nvertices'] = 81 + 9
self.nverticesO = 81
self.faultMesh = {
'nvertices': 9,
'spaceDim': 2,
'ncells': 8,
'ncorners': 2
}
run_pylith(SWShearSlidingApp)
self.outputRoot = "slipweakening_shear_sliding"
self.soln = AnalyticalSoln()
return
def test_fault_info(self):
"""
Check fault information.
"""
if not self.checkResults:
return
filename = "%s-fault_info.h5" % self.outputRoot
fields = [
"strike_dir", "normal_dir", "traction_initial",
"static_coefficient", "dynamic_coefficient",
"slip_weakening_parameter", "cohesion"
]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
def test_fault_data(self):
"""
Check fault information.
"""
if not self.checkResults:
return
filename = "%s-fault.h5" % self.outputRoot
fields = ["slip", "traction", "cumulative_slip", "previous_slip"]
check_vertex_fields(self, filename, self.faultMesh, fields)
return
def calcDisplacements(self, vertices):
"""
Calculate displacement field given coordinates of vertices.
"""
return self.soln.displacement(vertices, self.nverticesO)
def calcStateVar(self, name, vertices, cells):
"""
Calculate state variable.
"""
ncells = self.mesh['ncells']
pts = numpy.zeros((ncells, 3), dtype=numpy.float64)
if name == "total_strain":
stateVar = self.soln.strain(pts)
elif name == "stress":
stateVar = self.soln.stress(pts)
else:
raise ValueError("Unknown state variable '%s'." % name)
return stateVar
def calcFaultField(self, name, vertices):
"""
Calculate fault info.
"""
strikeDir = (0.0, -1.0)
normalDir = (-1.0, 0.0)
initialTraction = (0.0, -1.0e+6)
staticCoefficient = 0.61
dynamicCoefficient = 0.60
slipWeakeningParameter = 0.2
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)
cumulativeSlip = D
nvertices = self.faultMesh['nvertices']
if name == "strike_dir":
field = numpy.zeros((1, nvertices, 2), dtype=numpy.float64)
field[0, :, 0] = strikeDir[0]
field[0, :, 1] = strikeDir[1]
elif name == "normal_dir":
field = numpy.zeros((1, nvertices, 2), dtype=numpy.float64)
field[0, :, 0] = normalDir[0]
field[0, :, 1] = normalDir[1]
elif name == "traction_initial":
field = numpy.zeros((1, nvertices, 2), dtype=numpy.float64)
field[0, :, 0] = initialTraction[0]
field[0, :, 1] = initialTraction[1]
elif name == "static_coefficient":
field = numpy.zeros((1, nvertices, 1), dtype=numpy.float64)
field[0, :, 0] = staticCoefficient
elif name == "dynamic_coefficient":
field = numpy.zeros((1, nvertices, 1), dtype=numpy.float64)
field[0, :, 0] = dynamicCoefficient
elif name == "slip_weakening_parameter":
field = numpy.zeros((1, nvertices, 1), dtype=numpy.float64)
field[0, :, 0] = slipWeakeningParameter
elif name == "cohesion":
field = numpy.zeros((1, nvertices, 1), dtype=numpy.float64)
elif name == "slip":
field = numpy.zeros((1, nvertices, 2), dtype=numpy.float64)
field[0, :, 0] = slip[0]
field[0, :, 1] = slip[1]
elif name == "traction":
field = numpy.zeros((1, nvertices, 2), dtype=numpy.float64)
field[0, :, 0] = 0.6e+6
field[0, :, 1] = -1.0e+6
elif name == "cumulative_slip":
field = numpy.zeros((1, nvertices, 1), dtype=numpy.float64)
field[0, :, 0] = cumulativeSlip
elif name == "previous_slip":
field = numpy.zeros((1, nvertices, 1), dtype=numpy.float64)
field[0, :, 0] = cumulativeSlip
else:
raise ValueError("Unknown fault field '%s'." % name)
return field