[CIG-SHORT] Question about traction application (2)

Li, Teng tengli2 at illinois.edu
Fri Mar 31 22:30:45 PDT 2017 

Hi,
Here I attach the code of this problem for your convenience. And also, please see the attached mesh picture.

[pylithapp]
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
# Turn on some journals to show progress.
[pylithapp.journal.info]
pylithapp = 1
timedependent = 1
explicit = 1
pets = 1
meshiocubit = 1
explicitelasticity = 1
quadrature2d = 1
fiatlagrange = 1
faultcohesivedyn = 1
[pylithapp.journal.debug]
pylithapp = 1
problem = 1
explicit = 1

# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
# uncomment to get very verbose mesh information
#debug = 1

# Change the default mesh reader to the CUBIT reader.
reader = pylith.meshio.MeshIOCubit
#refiner = pylith.topology.RefineUniform
[pylithapp.mesh_generator.reader]
# Set filename of mesh to read and dimension for coordinate system.
filename =inter_long.exo
coordsys.space_dim = 2

# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
# Set the spatial dimension of the problem.
dimension = 2

# Change to an explicit time stepping formulation with no elastic prestep.
formulation = pylith.problems.Explicit
elastic_prestep = False

# Nondimensionalize problem using wave propagation parameters.
normalizer = spatialdata.units.NondimElasticDynamic
normalizer.shear_wave_speed = 1.0*km/s

# Set bc to an array with 4 boundary conditions: 'x_neg', 'x_pos', 'y_neg', and 'y_pos'
bc = [x_pos,x_neg,y_neg,y_pos]

# Set interfaces to an array with 3 faults: 'fault' 'fault1' 'fault2'.
interfaces = [fault,fault1]

# Set materials to an array with 1 material 'elastic'.
materials = [elastic]

[pylithapp.timedependent.formulation.time_step]
# Set the total time of the simulation and the time step.
total_time = 0.2*s
dt = 0.00098*s
# This is not a self-contained simulation configuration file. This
# file only specifies the general parameters common to the dynamic
# fault (friction interface) simulations in this directory.
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
[pylithapp.timedependent.materials]
# Use elastic plane strain bulk constitutive model.
elastic = pylith.materials.ElasticPlaneStrain

[pylithapp.timedependent.materials.elastic]

# Label for material
label = Elastic material

# Nodeset id from CUBIT
id = 1

# Spatial database with physical properties for elastic material
db_properties.label = Elastic properties
db_properties.iohandler.filename = matprops.spatialdb

# Set the basis functions and quadrature:
# 2-D Lagrange cell with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
quadrature.cell.quad_order = 2

# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc]
# Change the default Dirichlet BC to the AbsorbingDampers BC
x_pos = pylith.bc.AbsorbingDampers
x_neg = pylith.bc.AbsorbingDampers
y_neg = pylith.bc.AbsorbingDampers
y_pos = pylith.bc.AbsorbingDampers


#-----------------------------------------------------------------------
# x positive
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.x_pos]
# Absorbing BC on +x face of bar

# Name of nodeset in CUBIT
label = face_xpos

# Specify label and filename for spatial database with physical properties
db.label = Absorbing BC +x
db.iohandler.filename = matprops.spatialdb

# Specify the basis functions and quadrature:
# 1-D Lagrange cell in 2-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2


#-----------------------------------------------------------------------
# x negative
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.x_neg]
# Absorbing BC on -x face of bar

# Name of nodeset in CUBIT
label = face_xneg

# Specify label and filename for spatial database with physical properties
db.label = Absorbing BC -x
db.iohandler.filename = matprops.spatialdb

# Specify the basis functions and quadrature:
# 1-D Lagrange cell in 2-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2


#-----------------------------------------------------------------------
# y negative
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.y_neg]
# Absorbing BC on -y face of bar

# Name of nodeset in CUBIT
label = face_yneg

# Specify label and filename for spatial database with physical properties
db.label = Absorbing BC -y
db.iohandler.filename = matprops.spatialdb

# Specify the basis functions and quadrature:
# 1-D Lagrange cell in 2-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2



#-----------------------------------------------------------------------
# y positive
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.y_pos]
# Absorbing BC on +y face of bar

# Name of nodeset in CUBIT
label = face_ypos

# Specify label and filename for spatial database with physical properties
db.label = Absorbing BC +y
db.iohandler.filename = matprops.spatialdb

# Specify the basis functions and quadrature:
# 1-D Lagrange cell in 2-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2


# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------

# ----------------------------------------------------------------------
# fault
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces]
# Change fault to dynamic fault interface.
fault = pylith.faults.FaultCohesiveDyn
fault1 = pylith.faults.FaultCohesiveDyn

[pylithapp.timedependent.interfaces.fault]

# Material id
id = 100

# Name of nodeset from CUBIT
label = fault
# edge = inter

# Specify the basis functions and quadrature:
# 1-D Lagrange cell in 2-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2

# slip-weakening friction
friction = pylith.friction.SlipWeakening
friction.label = Slip weakening


# Set the slip-weakening friction model parameters.
#  static coefficient of friction: 0.6
#  dynamic coefficient of friction: 0.3
#  slip-weakening parameter: 0.2 m
#  cohesion: 0 Pa
friction.db_properties = spatialdata.spatialdb.UniformDB
friction.db_properties.label = Slip weakening
friction.db_properties.values = [static-coefficient,dynamic-coefficient,slip-weakening-parameter,cohesion]
friction.db_properties.data = [0.6,0.3,0.2*m,0.0*Pa]
# Specify the initial tractions and a temporal perturbation on the
# fault using a uniform DB.
#
# Initial
#   shear: 75.0 MPa (right-lateral)
#   normal 120 MPa (compressive)
# Perturbation
#   shear: 25.0 MPa (right-lateral)
traction_perturbation = pylith.faults.TractPerturbation

[pylithapp.timedependent.interfaces.fault.traction_perturbation]
#db_initial = spatialdata.spatialdb.UniformDB
#db_initial.label ZZ= Initial fault tractions
#db_initial.values = [traction-shear,traction-normal]
#db_initial.data = [-75.0*MPa, -120.0*MPa]
db_initial = spatialdata.spatialdb.SimpleDB
db_initial.label = Initial fault tractions

# Spatial database with physical properties for elastic material
db_initial.iohandler.filename = traction.spatialdb




# ----------------------------------------------------------------------
# fault1
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault1]

# Material id
id = 101

# Name of nodeset from CUBIT
label = fault1
edge = inter

# Specify the basis functions and quadrature:
# 1-D Lagrange cell in 2-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2

# slip-weakening friction
friction = pylith.friction.SlipWeakening
friction.label = Slip weakening

# Set the slip-weakening friction model parameters.
#  static coefficient of friction: 0.6
#  dynamic coefficient of friction: 0.3
#  slip-weakening parameter: 0.2 m
#  cohesion: 0 Pa
friction.db_properties = spatialdata.spatialdb.UniformDB
friction.db_properties.label = Slip weakening
friction.db_properties.values = [static-coefficient,dynamic-coefficient,slip-weakening-parameter,cohesion]
friction.db_properties.data = [0.6,0.3,0.2*m,0.0*Pa]
# Specify the initial tractions and a temporal perturbation on the
# fault using a uniform DB.
#
# Initial
#   shear: 75.0 MPa (right-lateral)
#   normal 120 MPa (compressive)
# Perturbation
#   shear: 25.0 MPa (right-lateral)
traction_perturbation = pylith.faults.TractPerturbation

[pylithapp.timedependent.interfaces.fault1.traction_perturbation]
#db_initial = spatialdata.spatialdb.UniformDB
#db_initial.label = Initial fault1 tractions
#db_initial.values = [traction-shear,traction-normal]
#db_initial.data = [-75.0*MPa, -120.0*MPa]
db_initial = spatialdata.spatialdb.SimpleDB
db_initial.label = Initial fault1 tractions

# Spatial database with physical properties for elastic material
db_initial.iohandler.filename = traction1.spatialdb





# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
[pylithapp.petsc]
# PETSc summary -- useful for performance information.
log_summary = true


# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
[pylithapp.problem.formulation]
output = [domain]
#----Domain velocity
[pylithapp.problem.formulation.output.domain]
vertex_data_fields = [displacement,velocity]
skip= 10
#writer = pylith.meshio.DataWriterHDF5
writer.filename = output/domain.vtk

#----Fault
[pylithapp.problem.interfaces.fault.output]
vertex_info_fields = []
vertex_data_fields = [slip,slip_rate,traction]
writer = pylith.meshio.DataWriterHDF5
skip = 0
writer.filename = output/fault.h5

#----Fault1
[pylithapp.problem.interfaces.fault1.output]
vertex_info_fields = []
vertex_data_fields = [slip,slip_rate,traction]
writer = pylith.meshio.DataWriterHDF5
skip = 0
writer.filename = output/fault1.h5

#----Fault2
# [pylithapp.problem.interfaces.fault2.output]
# vertex_info_fields = []
# vertex_data_fields = [slip,slip_rate,traction]
# writer = pylith.meshio.DataWriterHDF5
# skip = 0
# writer.filename = output/fault2.h5

#----Stress
[pylithapp.problem.materials.elastic.output]
cell_info_fields = []
cell_data_fields = [stress]
cell_filter = pylith.meshio.CellFilterAvg
writer.filename = ouptut/material.vtk
skip=150




Teng Li

Master Candidate in Structures

Department of Civil and Environmental Engineering

University of Illinois at Urbana-Champaign

205 North Mathews Ave, Urbana, IL. 61801

Phone:(217)8196210, Email: tengli2 at illinois.edu


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170401/1b8e6c11/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Mesh.jpg
Type: image/jpeg
Size: 33336 bytes
Desc: Mesh.jpg
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170401/1b8e6c11/attachment-0002.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Fault.jpg
Type: image/jpeg
Size: 47207 bytes
Desc: Fault.jpg
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170401/1b8e6c11/attachment-0003.jpg>

More information about the CIG-SHORT mailing list