[cig-commits] r16711 - short/3D/PyLith/trunk/examples/bar_shearwave/quad4
brad at geodynamics.org
brad at geodynamics.org
Fri May 14 21:30:23 PDT 2010
Author: brad
Date: 2010-05-14 21:30:22 -0700 (Fri, 14 May 2010)
New Revision: 16711
Modified:
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-sliding.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-stick.cfg
Log:
Worked on bar shear wave examples.
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-15 04:30:04 UTC (rev 16710)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-15 04:30:22 UTC (rev 16711)
@@ -1,105 +1,28 @@
# -*- Python -*-
[pylithapp]
+# 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.
+
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
[pylithapp.journal.info]
-timedependent = 1
-#explicit = 1
-#petsc = 1
-#meshiocubit = 1
-#explicitelasticity = 1
-#quadrature2d = 1
-#fiatlagrange = 1
-#faultcohesivedyn = 1
+faultcohesivedyn = 1
-[pylithapp.journal.debug]
-#pylithapp = 1
-#problem = 1
-#explicit = 1
-
# ----------------------------------------------------------------------
-# mesh_generator
+# faults
# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-reader = pylith.meshio.MeshIOCubit
+[pylithapp.timedependent.interfaces]
+# Change fault to dynamic fault interface.
+fault = pylith.faults.FaultCohesiveDyn
-[pylithapp.mesh_generator.reader]
-filename = bar_quad4_200m.exo
-coordsys.space_dim = 2
-
# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-dimension = 2
-
-# Set interfaces to an array with 1 fault: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 1 material 'elastic'.
-materials = [elastic]
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-normalizer.shear_wave_speed = 1.0*km/s
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 0.0*s
-dt = 0.05*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials]
-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.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
-
-# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
[pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-# Change the preconditioner settings.
-sub_pc_factor_shift_type = nonzero
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-15
-ksp_max_it = 200
-ksp_gmres_restart = 250
-snes_max_it = 300
-
-#ksp_monitor = true
-#ksp_view = true
-ksp_converged_reason = true
-
-snes_rtol = 1.0e-8
-snes_atol = 1.0e-15
-#snes_monitor = true
-#snes_view = true
-snes_converged_reason = true
-
-#log_summary = true
-
# Friction sensitivity solve
friction_pc_type = asm
friction_sub_pc_factor_shift_type = nonzero
@@ -109,31 +32,9 @@
#friction_ksp_view = true
friction_ksp_converged_reason = true
-
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-# We set three output manager parameters for each type of output.
-# skip - Number of time steps to skip between output
-# writer.filename - Root name for VTK file
-# writer.time_format - Format of time stamp (decimal point is removed)
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-vertex_data_fields=[displacement,velocity]
-skip = 3
-writer.time_format = %05.2f
-
-# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
+# Add slip rate to the default output fiels of slip and traction.
vertex_data_fields=[slip,slip_rate,traction]
-skip = 3
-writer.time_format = %05.2f
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-skip = 33
-writer.time_format = %05.2f
-
-# Averate state variables over cell
-cell_filter = pylith.meshio.CellFilterAvgMesh
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg 2010-05-15 04:30:04 UTC (rev 16710)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/kinematicfault.cfg 2010-05-15 04:30:22 UTC (rev 16711)
@@ -1,95 +1,34 @@
# -*- Python -*-
-
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
-
[pylithapp]
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for an explicit
+# time-stepping problem using a fault with kinematically-specified
+# slip.
+
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
[pylithapp.journal.info]
-timedependent = 1
-#explicit = 1
-#petsc = 1
-meshiocubit = 1
-#explicitelasticity = 1
-#quadrature2d = 1
-#fiatlagrange = 1
-#faultcohesivekin = 1
+faultcohesivekin = 1
-[pylithapp.journal.debug]
-pylithapp = 1
-problem = 1
-explicit = 1
-
# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1
-reader = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.reader]
-filename = bar_quad4_200m.exo
-coordsys.space_dim = 2
-
-# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
dimension = 2
-# Change to an explicit time stepping formulation
-#formulation = pylith.problems.Explicit
-formulation = pylith.problems.ExplicitLumped
+# Uncomment the line below to switch to explicit time-stepping with a
+# lumped Jacobian matrix.
+#formulation = pylith.problems.ExplicitLumped
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'all_nofault'.
-bc = [x_pos,x_neg,all_nofault]
-
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
# Change the DirichletPoints BC to the AbsorbingDampers BC
bc.x_pos = pylith.bc.AbsorbingDampers
bc.x_neg = pylith.bc.AbsorbingDampers
-# Set interfaces to an array with 1 fault: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 1 material 'elastic'.
-materials = [elastic]
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.05*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials]
-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.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.x_pos]
# Absorbing boundary condition on +x face of bar
@@ -122,102 +61,38 @@
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2
-[pylithapp.timedependent.bc.all_nofault]
-# Dirichlet boundary condition on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
[pylithapp.timedependent.interfaces.fault]
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# 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
-
# Switch to Brune slip time function
eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
# Database specifying the final slip.
-slip.iohandler.filename = shearwave_slip.spatialdb
+slip.iohandler.filename = kinematic_slip.spatialdb
# Database specifying rise_time.
-rise_time.iohandler.filename = shearwave_risetime.spatialdb
+rise_time.iohandler.filename = kinematic_risetime.spatialdb
# Database specifying time at which slip begins at each point.
-slip_time.iohandler.filename = shearwave_sliptime.spatialdb
+slip_time.iohandler.filename = kinematic_sliptime.spatialdb
# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-
-# Change the preconditioner settings.
-sub_pc_factor_shift_type = nonzero
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-12
-ksp_max_it = 50
-ksp_gmres_restart = 75
-
-ksp_monitor = true
-ksp_view = true
-#log_summary = true
-
-
-# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-# We set three output manager parameters for each type of output.
-# skip - Number of time steps to skip between output
-# writer.filename - Root name for VTK file
-# writer.time_format - Format of time stamp (decimal point is removed)
+# Set the filename for each type of output.
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-skip = 3
-writer.filename = output/shearwave.vtk
-writer.time_format = %05.2f
+writer.filename = output/kinematicfault.vtk
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
-writer.filename = output/shearwave-fault.vtk
-writer.time_format = %05.2f
+writer.filename = output/kinematicfault-fault.vtk
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-skip = 33
-writer.filename = output/shearwave-statevars.vtk
-writer.time_format = %05.2f
-
-# Averate state variables over cell
-cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = output/kinematicfault-statevars.vtk
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2010-05-15 04:30:04 UTC (rev 16710)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2010-05-15 04:30:22 UTC (rev 16711)
@@ -1,22 +1,22 @@
# -*- Python -*-
-
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
-
[pylithapp]
+# This is not a self-contained simulation configuration file. This
+# file only specifies the general parameters common to the simulations
+# in this directory.
+
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
+# Turn on some journals to show progress.
[pylithapp.journal.info]
timedependent = 1
-#explicit = 1
-#petsc = 1
+explicit = 1
+petsc = 1
meshiocubit = 1
-#explicitelasticity = 1
-#quadrature2d = 1
-#fiatlagrange = 1
-#faultcohesivekin = 1
+explicitelasticity = 1
+quadrature2d = 1
+fiatlagrange = 1
[pylithapp.journal.debug]
pylithapp = 1
@@ -27,10 +27,14 @@
# 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
[pylithapp.mesh_generator.reader]
+# Set filename of mesh to read and dimension for coordinate system.
filename = bar_quad4_200m.exo
coordsys.space_dim = 2
@@ -40,20 +44,18 @@
[pylithapp.timedependent]
dimension = 2
-# Change to an explicit time stepping formulation
-#formulation = pylith.problems.Explicit
-formulation = pylith.problems.ExplicitLumped
+# Change to an explicit time stepping formulation. We use explicit
+# time-stepping with a consistent (full) Jacobian matrix.
+formulation = pylith.problems.Explicit
# Nondimensionalize problem using wave propagation parameters.
normalizer = spatialdata.units.NondimElasticDynamic
+normalizer.shear_wave_speed = 1.0*km/s
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
-bc = [x_pos,x_neg,y_pos]
+# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and
+# 'all_nofault'.
+bc = [x_pos,x_neg,all_nofault]
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
@@ -90,39 +92,7 @@
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing boundary condition on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# 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
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing boundary condition on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-# 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
-
-[pylithapp.timedependent.bc.y_pos]
+[pylithapp.timedependent.bc.all_nofault]
# Dirichlet boundary condition on all vertices except fault vertices
# Pin x and z degrees of freedom
@@ -159,65 +129,55 @@
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2
-# Switch to Brune slip time function
-eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
-
-[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
-
-# Database specifying the final slip.
-slip.iohandler.filename = shearwave_slip.spatialdb
-
-# Database specifying rise_time.
-rise_time.iohandler.filename = shearwave_risetime.spatialdb
-
-# Database specifying time at which slip begins at each point.
-slip_time.iohandler.filename = shearwave_sliptime.spatialdb
-
# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
[pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-# Change the preconditioner settings.
+# Preconditioner settings.
+pc_type = asm
sub_pc_factor_shift_type = nonzero
+# Convergence parameters.
ksp_rtol = 1.0e-8
ksp_atol = 1.0e-12
-ksp_max_it = 50
-ksp_gmres_restart = 75
+ksp_max_it = 500
+ksp_gmres_restart = 100
+snes_max_it = 500
+# Linear solver monitoring options.
ksp_monitor = true
ksp_view = true
-#log_summary = true
+ksp_converged_reason = true
+# Nonlinear solver monitoring options.
+snes_monitor = true
+snes_view = true
+snes_converged_reason = true
+# PETSc summary -- useful for performance information.
+log_summary = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-# We set three output manager parameters for each type of output.
+# We set two output manager parameters for each type of output.
# skip - Number of time steps to skip between output
-# writer.filename - Root name for VTK file
# writer.time_format - Format of time stamp (decimal point is removed)
-
-# Give basename for VTK domain output of solution over domain.
+# We will set the output filename in the problem specific .cfg file.
[pylithapp.problem.formulation.output.output]
skip = 3
-writer.filename = output/shearwave.vtk
writer.time_format = %05.2f
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
skip = 3
-writer.filename = output/shearwave-fault.vtk
writer.time_format = %05.2f
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
skip = 33
-writer.filename = output/shearwave-statevars.vtk
writer.time_format = %05.2f
-# Averate state variables over cell
+# Averate state variables over quadrature points.
cell_filter = pylith.meshio.CellFilterAvgMesh
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-sliding.cfg 2010-05-15 04:30:04 UTC (rev 16710)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-sliding.cfg 2010-05-15 04:30:22 UTC (rev 16711)
@@ -1,22 +1,34 @@
# -*- Python -*-
-
[pylithapp]
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for a static problem
+# using a fault with a fault constitutive model.
+
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-# Change to an explicit time stepping formulation
+# Change to an implicit time-stepping formulation and the nonlinear
+# solver.
formulation = pylith.problems.Implicit
formulation.solver = pylith.problems.SolverNonlinear
# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
bc = [x_pos,x_neg,all_nofault]
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc]
+# Switch from absorbing boundary condition to Dirichlet BC.
+x_pos = pylith.bc.DirichletBC
+x_neg = pylith.bc.DirichletBC
+
[pylithapp.timedependent.bc.x_pos]
# Dirichlet BC on +x face of bar
@@ -29,7 +41,6 @@
db_initial.values = [displacement-x,displacement-y]
db_initial.data = [0.0*m,32.0*km]
-
[pylithapp.timedependent.bc.x_neg]
# Dirichlet BC on -x face of bar
@@ -38,38 +49,12 @@
bc_dof = [1]
-[pylithapp.timedependent.bc.all_nofault]
-# Dirichlet BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
# ----------------------------------------------------------------------
# faults
-# ----------------------------------------------------------------------
-
+# ----------------------------------------------------------------------x
# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-
-fault = pylith.faults.FaultCohesiveDyn
-
[pylithapp.timedependent.interfaces.fault]
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# 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
-
db_initial_tractions = spatialdata.spatialdb.UniformDB
db_initial_tractions.label = Initial fault tractions
db_initial_tractions.values = [traction-shear,traction-normal]
@@ -88,12 +73,12 @@
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-writer.filename = output/shear-sliding.vtk
+writer.filename = output/static-sliding.vtk
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/shear-sliding-fault.vtk
+writer.filename = output/static-sliding-fault.vtk
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/shear-sliding-statevars.vtk
+writer.filename = output/static-sliding-statevars.vtk
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-stick.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-stick.cfg 2010-05-15 04:30:04 UTC (rev 16710)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static-stick.cfg 2010-05-15 04:30:22 UTC (rev 16711)
@@ -1,22 +1,34 @@
# -*- Python -*-
-
[pylithapp]
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for a static problem
+# using a fault with a fault constitutive model.
+
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-# Change to an explicit time stepping formulation
+# Change to an implicit time-stepping formulation and the nonlinear
+# solver.
formulation = pylith.problems.Implicit
formulation.solver = pylith.problems.SolverNonlinear
# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
bc = [x_pos,x_neg,all_nofault]
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc]
+# Switch from absorbing boundary condition to Dirichlet BC.
+x_pos = pylith.bc.DirichletBC
+x_neg = pylith.bc.DirichletBC
+
[pylithapp.timedependent.bc.x_pos]
# Dirichlet BC on +x face of bar
@@ -29,7 +41,6 @@
db_initial.values = [displacement-x,displacement-y]
db_initial.data = [0.0*m,10.0*km]
-
[pylithapp.timedependent.bc.x_neg]
# Dirichlet BC on -x face of bar
@@ -38,38 +49,11 @@
bc_dof = [1]
-[pylithapp.timedependent.bc.all_nofault]
-# Dirichlet BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
-
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-
-fault = pylith.faults.FaultCohesiveDyn
-
[pylithapp.timedependent.interfaces.fault]
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# 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
-
db_initial_tractions = spatialdata.spatialdb.UniformDB
db_initial_tractions.label = Initial fault tractions
db_initial_tractions.values = [traction-shear,traction-normal]
@@ -88,12 +72,12 @@
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-writer.filename = output/shear-stick.vtk
+writer.filename = output/static-stick.vtk
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/shear-stick-fault.vtk
+writer.filename = output/static-stick-fault.vtk
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/shear-stick-statevars.vtk
+writer.filename = output/static-stick-statevars.vtk
More information about the CIG-COMMITS
mailing list