[cig-commits] r16399 - in short/3D/PyLith/trunk/playpen/friction: . RateStateAgeing RateStateAgeing/quad4
surendra at geodynamics.org
surendra at geodynamics.org
Tue Mar 9 17:25:12 PST 2010
Author: surendra
Date: 2010-03-09 17:25:11 -0800 (Tue, 09 Mar 2010)
New Revision: 16399
Added:
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/README
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/bar_quad4_200m.exo
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/geometry.jou
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/matprops.spatialdb
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/mesh_quad4_200m.jou
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/output/
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylith_parameters.txt
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylithapp.cfg
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_risetime.spatialdb
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_slip.spatialdb
short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_sliptime.spatialdb
Log:
Added bar shearwave example for static friction
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/README
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/README (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/README 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,18 @@
+The bar is discretized with quadrilateral cells.
+
+Files associated with the example problem:
+
+ - geometry.jou: CUBIT journal file defining the geometry of the mesh
+ - mesh_quad4_200m.jou: CUBIT journal file for meshing the problem
+ - pylithapp.cfg: Parameter settings for the shear wave dislocation problem
+ - matprops.spatialdb: Material database of physical properties
+ - bar_quad4_200m.exo: Exodus file with mesh generated by CUBIT
+ - shearwave_slip.spatialdb: Spatial database with values for final slip
+ - shearwave_sliprate.spatialdb: Spatial database with values for the
+ peak slip rate
+ - shearwave_sliptime.spatialdb: Spatial database with values for the
+ slip initiation time
+
+To run the problem, enter
+
+ pylith
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/bar_quad4_200m.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/bar_quad4_200m.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/geometry.jou
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/geometry.jou (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/geometry.jou 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,36 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Bar is 400m x 8000m
+# -4000 m <= x <= 4000 m
+# -200 m <= y <= 200 m
+#
+# We create the 2-D bar by extracting the mid-surface from a 3-D bar.
+reset
+brick x 8000 y 400 z 400
+create midsurface volume 1 surface 1 2
+delete volume 1
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+
+# Create fault (yz plane) at x = 0.0
+split surface 8 across location position 0 -200 0 location position 0 200 0
+curve 21 name "fault_surface"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all
+merge all
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/matprops.spatialdb
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/matprops.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/matprops.spatialdb 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ // number of physical properties
+ num-values = 3
+
+ // Names and units of physical properties
+ value-names = density vs vp
+ value-units = kg/m**3 m/s m/s
+
+ // Number of locations where physical properties are given
+ // 1 location -> data is uniform -> data dimension is 0
+ num-locs = 1
+ data-dim = 0
+
+ // Problem is in 2-D
+ space-dim = 2
+
+ // Coordinate system (2-D with coordinates in meters)
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+// Columns are x, y, density, vs, vp
+0.0 0.0 2500.0 1000.0 1732.0
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/mesh_quad4_200m.jou
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/mesh_quad4_200m.jou (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/mesh_quad4_200m.jou 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,63 @@
+## /tools/common/cubit-10.2/bin/clarox
+## Cubit Version 10.2
+## Cubit Build 24
+## Revised 12/15/2006 16:09:40 MST
+## Running 06/18/2007 10:26:50 AM
+## Command Options:
+## -warning = On
+## -information = On
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface all size 200
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 surface 9 10
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in curve 18
+nodeset 11 group face_xpos
+nodeset 11 name "face xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in curve 20
+nodeset 12 group face_xneg
+nodeset 12 name "face xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for all vertices but fault
+# ----------------------------------------------------------------------
+group "vertices_nofault" add node all
+group "vertices_nofault" remove node in fault
+nodeset 13 group vertices_nofault
+nodeset 13 name "vertices_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "bar_quad4_200m.exo" dimension 2 overwrite
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylith_parameters.txt
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylith_parameters.txt (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylith_parameters.txt 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,238 @@
+Application: pylithapp <pylith.apps.PyLithApp.InfoApp object at 0x48ee50>
+ nodes (int) = 1
+ job = job (<pyre.schedulers.Job.Job object at 0x3b56fb0>)
+ arguments (list) = []
+ comments (list) = []
+ environment (list) = []
+ executable (str) =
+ mail (bool) = False
+ name (str) =
+ queue (str) =
+ stderr (str) = stderr.txt
+ stdin (str) = /dev/null
+ stdout (str) = stdout.txt
+ walltime (dimensional) = 0*s
+ launcher = mpich (<mpi.LauncherMPICH.LauncherMPICH object at 0x3b5c4d0>)
+ arguments (list) = []
+ command (str) = mpirun -np ${nodes}
+ dry (bool) = False
+ executable (str) =
+ machinefile (str) = mpirun.nodes
+ nodegen (str) =
+ nodelist (slice) = []
+ nodes (int) = 1
+ mesh_generator = meshimporter (<pylith.topology.MeshImporter.MeshImporter object at 0x3b3cc30>)
+ debug (bool) = False
+ interpolate (bool) = False
+ distributor = partitioner (<pylith.topology.Distributor.Distributor; proxy of <Swig Object of type 'pylith::topology::Distributor *' at 0x3b2f9c0> >)
+ debug (bool) = False
+ partitioner (str) = chaco
+ data_writer = datawritervtkmesh (<pylith.meshio.DataWriterVTKMesh.DataWriterVTKMesh; proxy of <Swig Object of type 'pylith::meshio::DataWriterVTK< pylith::topology::Mesh,pylith::topology::Field< pylith::topology::Mesh > > *' at 0x3b2ff40> >)
+ filename (str) = output.vtk
+ float_precision (int) = 6
+ time_constant (dimensional) = 1*s
+ time_format (str) = %f
+ reader = meshiocubit (<pylith.meshio.MeshIOCubit.MeshIOCubit; proxy of <Swig Object of type 'pylith::meshio::MeshIOCubit *' at 0x3b2ff00> >)
+ filename (str) = bar_quad4_200m.exo
+ coordsys = cscart (<spatialdata.geocoords.CSCart.CSCart; proxy of <Swig Object of type 'spatialdata::geocoords::CSCart *' at 0x3b2ff20> >)
+ space_dim (int) = 2
+ units (str) = m
+ refiner = refiner (<pylith.topology.MeshRefiner.MeshRefiner object at 0x3b6d950>)
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3b56750>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ petsc = petsc (<pylith.utils.PetscManager.PetscManager object at 0x3b56df0>)
+ problem = timedependent (<pylith.problems.TimeDependent.TimeDependent object at 0x3b56850>)
+ dimension (int) = 2
+ bc = bc (<pyre.inventory.FacilityArrayFacility.FacilityArray object at 0x3b9c770>)
+ x_neg = absorbingdampers (<pylith.bc.AbsorbingDampers.AbsorbingDampers; proxy of <Swig Object of type 'pylith::bc::AbsorbingDampers *' at 0x3b2f9c0> >)
+ label (str) = 12
+ up_dir (list) = [0, 0, 1]
+ db = simpledb (<spatialdata.spatialdb.SimpleDB.SimpleDB; proxy of <Swig Object of type 'spatialdata::spatialdb::SimpleDB *' at 0x3b2fbe0> >)
+ label (str) = Absorbing BC -x
+ query_type (str) = nearest
+ iohandler = simpleioascii (<spatialdata.spatialdb.SimpleIOAscii.SimpleIOAscii; proxy of <Swig Object of type 'spatialdata::spatialdb::SimpleIOAscii *' at 0x3b2fc80> >)
+ filename (str) = matprops.spatialdb
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3c315d0>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ quadrature = submeshquadrature (<pylith.feassemble.Quadrature.SubMeshQuadrature; proxy of <Swig Object of type 'pylith::feassemble::Quadrature< pylith::topology::SubMesh > *' at 0x3b2fc80> >)
+ check_conditioning (bool) = False
+ min_jacobian (float) = 1e-06
+ cell = fiatlagrange (<pylith.feassemble.FIATLagrange.FIATLagrange object at 0x3cac090>)
+ degree (int) = 1
+ dimension (int) = 1
+ quad_order (int) = 2
+ x_pos = absorbingdampers (<pylith.bc.AbsorbingDampers.AbsorbingDampers; proxy of <Swig Object of type 'pylith::bc::AbsorbingDampers *' at 0x3b2fc80> >)
+ label (str) = 11
+ up_dir (list) = [0, 0, 1]
+ db = simpledb (<spatialdata.spatialdb.SimpleDB.SimpleDB; proxy of <Swig Object of type 'spatialdata::spatialdb::SimpleDB *' at 0x3b2fbe0> >)
+ label (str) = Absorbing BC +x
+ query_type (str) = nearest
+ iohandler = simpleioascii (<spatialdata.spatialdb.SimpleIOAscii.SimpleIOAscii; proxy of <Swig Object of type 'spatialdata::spatialdb::SimpleIOAscii *' at 0x3b2ff20> >)
+ filename (str) = matprops.spatialdb
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3cc4110>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ quadrature = submeshquadrature (<pylith.feassemble.Quadrature.SubMeshQuadrature; proxy of <Swig Object of type 'pylith::feassemble::Quadrature< pylith::topology::SubMesh > *' at 0x3b2ff20> >)
+ check_conditioning (bool) = False
+ min_jacobian (float) = 1e-06
+ cell = fiatlagrange (<pylith.feassemble.FIATLagrange.FIATLagrange object at 0x3cc4650>)
+ degree (int) = 1
+ dimension (int) = 1
+ quad_order (int) = 2
+ y_pos = dirichletbc (<pylith.bc.DirichletBC.DirichletBC; proxy of <Swig Object of type 'pylith::bc::DirichletBC *' at 0x3b2ff20> >)
+ bc_dof (list) = [0]
+ label (str) = 13
+ up_dir (list) = [0, 0, 1]
+ db_change = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3cd0670>)
+ db_initial = zerodispdb (<pylith.bc.ZeroDispDB.ZeroDispDB; proxy of <Swig Object of type 'spatialdata::spatialdb::UniformDB *' at 0x3b2fbe0> >)
+ data (list) = []
+ label (str) = spatial database
+ values (list) = []
+ db_rate = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3cd0570>)
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3cd07b0>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ th_change = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3cd0710>)
+ checkpoint = checkpointtimer (<pylith.utils.CheckpointTimer.CheckpointTimer object at 0x3b93490>)
+ dt (dimensional) = 9.9e+99*s
+ formulation = explicit (<pylith.problems.Explicit.Explicit; proxy of <Swig Object of type 'pylith::problems::Explicit *' at 0x3b2fce0> >)
+ matrix_type (str) = unknown
+ split_fields (bool) = False
+ view_jacobian (bool) = False
+ jacobian_viewer = formulation (<pylith.topology.JacobianViewer.JacobianViewer object at 0x3d560b0>)
+ filename (str) = jacobian.mat
+ time_constant (dimensional) = 1*s
+ time_format (str) = %f
+ output = singleoutput (<pylith.meshio.SingleOutput.SingleOutput object at 0x3d562d0>)
+ output = outputsoln (<pylith.meshio.OutputSoln.OutputSoln; proxy of <Swig Object of type 'pylith::meshio::OutputManager< pylith::topology::Mesh,pylith::topology::Field< pylith::topology::Mesh > > *' at 0x3b2fc80> >)
+ cell_info_fields (list) = []
+ output_freq (str) = skip
+ skip (int) = 3
+ time_step (dimensional) = 1*s
+ vertex_data_fields (list) = ['displacement']
+ cell_filter = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3d63690>)
+ coordsys = cscart (<spatialdata.geocoords.CSCart.CSCart; proxy of <Swig Object of type 'spatialdata::geocoords::CSCart *' at 0x3b2f3c0> >)
+ space_dim (int) = 3
+ units (str) = m
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3d635d0>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ vertex_filter = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3d637f0>)
+ writer = datawritervtkmesh (<pylith.meshio.DataWriterVTKMesh.DataWriterVTKMesh; proxy of <Swig Object of type 'pylith::meshio::DataWriterVTK< pylith::topology::Mesh,pylith::topology::Field< pylith::topology::Mesh > > *' at 0x3b2fdc0> >)
+ filename (str) = output/shearwave.vtk
+ float_precision (int) = 6
+ time_constant (dimensional) = 1*s
+ time_format (str) = %05.2f
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3d4dff0>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ solver = solvernonlinear (<pylith.problems.SolverNonlinear.SolverNonlinear; proxy of <Swig Object of type 'pylith::problems::SolverNonlinear *' at 0x3b2fc40> >)
+ time_step = timestepuniform (<pylith.problems.TimeStepUniform.TimeStepUniform object at 0x3d4d890>)
+ dt (dimensional) = 0.05*s
+ total_time (dimensional) = 2*s
+ gravity_field = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3b93c50>)
+ interfaces = interfaces (<pyre.inventory.FacilityArrayFacility.FacilityArray object at 0x3bbd1f0>)
+ fault = faultcohesivedynl (<pylith.faults.FaultCohesiveDyn.FaultCohesiveDyn; proxy of <Swig Object of type 'pylith::faults::FaultCohesiveDyn *' at 0x3b2f660> >)
+ fault_mesh_filename (str) = fault.inp
+ id (int) = 100
+ label (str) = 10
+ normal_dir (list) = [1.0, 0.0, 0.0]
+ up_dir (list) = [0.0, 0.0, 1.0]
+ use_fault_mesh (bool) = False
+ db_initial_tractions = uniformdb (<spatialdata.spatialdb.UniformDB.UniformDB; proxy of <Swig Object of type 'spatialdata::spatialdb::UniformDB *' at 0x3b2f340> >)
+ data (list) = ['7000.0*MPa', ' -10000.0*MPa']
+ label (str) = Initial fault tractions
+ values (list) = ['traction-shear', 'traction-normal']
+ friction = staticfriction (<pylith.friction.StaticFriction.StaticFriction; proxy of <Swig Object of type 'pylith::friction::StaticFriction *' at 0x3b2fca0> >)
+ label (str) =
+ db_initial_state = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3d0f130>)
+ db_properties = uniformdb (<spatialdata.spatialdb.UniformDB.UniformDB; proxy of <Swig Object of type 'spatialdata::spatialdb::UniformDB *' at 0x3b2f660> >)
+ data (list) = ['0.6*m']
+ label (str) = Static friction
+ values (list) = ['friction-coefficient']
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3d0f450>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ output = outputfaultdyn (<pylith.meshio.OutputFaultDyn.OutputFaultDyn; proxy of <Swig Object of type 'pylith::meshio::OutputManager< pylith::topology::SubMesh,pylith::topology::Field< pylith::topology::SubMesh > > *' at 0x3b2f340> >)
+ output_freq (str) = skip
+ skip (int) = 3
+ time_step (dimensional) = 1*s
+ vertex_data_fields (list) = ['slip', 'traction']
+ vertex_info_fields (list) = ['normal_dir']
+ cell_filter = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3cf7950>)
+ coordsys = cscart (<spatialdata.geocoords.CSCart.CSCart; proxy of <Swig Object of type 'spatialdata::geocoords::CSCart *' at 0x3b2fbc0> >)
+ space_dim (int) = 3
+ units (str) = m
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3cf7890>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ vertex_filter = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3cf7c30>)
+ writer = datawritervtksubmesh (<pylith.meshio.DataWriterVTKSubSubMesh.DataWriterVTKSubSubMesh; proxy of <Swig Object of type 'pylith::meshio::DataWriterVTK< pylith::topology::SubMesh,pylith::topology::Field< pylith::topology::SubMesh > > *' at 0x3b2f100> >)
+ filename (str) = output/shearwave-fault.vtk
+ float_precision (int) = 6
+ time_constant (dimensional) = 1*s
+ time_format (str) = %05.2f
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3cf7350>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ quadrature = submeshquadrature (<pylith.feassemble.Quadrature.SubMeshQuadrature; proxy of <Swig Object of type 'pylith::feassemble::Quadrature< pylith::topology::SubMesh > *' at 0x3b2f960> >)
+ check_conditioning (bool) = False
+ min_jacobian (float) = 1e-06
+ cell = fiatlagrange (<pylith.feassemble.FIATLagrange.FIATLagrange object at 0x3d03cd0>)
+ degree (int) = 1
+ dimension (int) = 1
+ quad_order (int) = 2
+ materials = materials (<pyre.inventory.FacilityArrayFacility.FacilityArray object at 0x3c22370>)
+ elastic = elasticplanestrain (<pylith.materials.ElasticPlaneStrain.ElasticPlaneStrain; proxy of <Swig Object of type 'pylith::materials::ElasticPlaneStrain *' at 0x3b2f960> >)
+ id (int) = 1
+ label (str) = Elastic material
+ db_initial_state = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3d22dd0>)
+ db_initial_strain = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3d22f10>)
+ db_initial_stress = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3d22d30>)
+ db_properties = simpledb (<spatialdata.spatialdb.SimpleDB.SimpleDB; proxy of <Swig Object of type 'spatialdata::spatialdb::SimpleDB *' at 0x3b2f660> >)
+ label (str) = spatial database
+ query_type (str) = nearest
+ iohandler = simpleioascii (<spatialdata.spatialdb.SimpleIOAscii.SimpleIOAscii; proxy of <Swig Object of type 'spatialdata::spatialdb::SimpleIOAscii *' at 0x3b2f3a0> >)
+ filename (str) = matprops.spatialdb
+ output = outputmatelastic (<pylith.meshio.OutputMatElastic.OutputMatElastic; proxy of <Swig Object of type 'pylith::meshio::OutputManager< pylith::topology::Mesh,pylith::topology::Field< pylith::topology::Mesh > > *' at 0x3b2f3a0> >)
+ cell_data_fields (list) = ['total_strain', 'stress']
+ cell_info_fields (list) = ['mu', 'lambda', 'density']
+ output_freq (str) = skip
+ skip (int) = 33
+ time_step (dimensional) = 1*s
+ cell_filter = cellfilteravgmesh (<pylith.meshio.CellFilterAvgMesh.CellFilterAvgMesh; proxy of <Swig Object of type 'pylith::meshio::CellFilterAvg< pylith::topology::Mesh,pylith::topology::Field< pylith::topology::Mesh > > *' at 0x3b2f980> >)
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3d42810>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ coordsys = cscart (<spatialdata.geocoords.CSCart.CSCart; proxy of <Swig Object of type 'spatialdata::geocoords::CSCart *' at 0x3b2f3a0> >)
+ space_dim (int) = 3
+ units (str) = m
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3d3c390>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ vertex_filter = nullcomponent (<pylith.utils.NullComponent.NullComponent object at 0x3d3cdb0>)
+ writer = datawritervtkmesh (<pylith.meshio.DataWriterVTKMesh.DataWriterVTKMesh; proxy of <Swig Object of type 'pylith::meshio::DataWriterVTK< pylith::topology::Mesh,pylith::topology::Field< pylith::topology::Mesh > > *' at 0x3b2f980> >)
+ filename (str) = output/shearwave-statevars.vtk
+ float_precision (int) = 6
+ time_constant (dimensional) = 1*s
+ time_format (str) = %05.2f
+ perf_logger = memory_logger (<pylith.perf.MemoryLogger.MemoryLogger object at 0x3d22e70>)
+ include_dealloc (bool) = True
+ verbose (bool) = False
+ quadrature = meshquadrature (<pylith.feassemble.Quadrature.MeshQuadrature; proxy of <Swig Object of type 'pylith::feassemble::Quadrature< pylith::topology::Mesh > *' at 0x3b2f640> >)
+ check_conditioning (bool) = False
+ min_jacobian (float) = 1e-06
+ cell = fiatlagrange (<pylith.feassemble.FIATLagrange.FIATLagrange object at 0x3d29610>)
+ degree (int) = 1
+ dimension (int) = 2
+ quad_order (int) = 2
+ normalizer = nondimelasticdynamic (<spatialdata.units.NondimElasticDynamic.NondimElasticDynamic; proxy of <Swig Object of type 'spatialdata::units::Nondimensional *' at 0x3b2f960> >)
+ mass_density (dimensional) = 3000*m**-3*kg
+ shear_wave_speed (dimensional) = 3000*m*s**-1
+ wave_period (dimensional) = 1*s
+ scheduler = scheduler-none (<pyre.schedulers.SchedulerNone.SchedulerNone object at 0x3b5cab0>)
+ dry (bool) = False
+ shell (str) = /bin/sh
+ wait (bool) = True
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/pylithapp.cfg 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,242 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+petsc = 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]
+#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.solver = pylith.problems.SolverNonlinear
+#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 'y_pos'.
+bc = [x_pos,x_neg,y_pos]
+
+# 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
+
+# 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]
+# 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.
+
+# 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]
+db_initial_tractions.data = [7000.0*MPa, -10000.0*MPa]
+
+friction = pylith.friction.StaticFriction
+
+friction.db_properties = spatialdata.spatialdb.UniformDB
+friction.db_properties.label = Static friction
+friction.db_properties.values = [friction-coefficient]
+friction.db_properties.data = [0.6*m] ; units required (use SI so not scaled)
+
+# 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 (must turn off
+# shift_positive_definite and turn on shift_nonzero).
+sub_pc_factor_shift_positive_definite = 0
+sub_pc_factor_shift_nonzero =
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-10
+ksp_max_it = 50
+ksp_gmres_restart = 10
+
+ksp_monitor = true
+ksp_view = true
+log_summary = true
+
+snes_monitor = true
+snes_view = true
+snes_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]
+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
+cell_filter = pylith.meshio.CellFilterAvgMesh
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_risetime.spatialdb
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_risetime.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_risetime.spatialdb 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ // number of values for rise time
+ num-values = 1
+
+ // Name and units for rise time
+ value-names = rise-time
+ value-units = s
+
+ // Number of locations where physical properties are given
+ // 1 location -> data is uniform -> data dimension is 0
+ num-locs = 1
+ data-dim = 0
+
+ // Problem is in 2-D
+ space-dim = 2
+
+ // Coordinate system (2-D with coordinates in meters)
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+// Columns are x, y, rise-time
+0.0 0.0 2.0
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_slip.spatialdb
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_slip.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_slip.spatialdb 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ // number of values for final slip
+ num-values = 2
+
+ // Names and units of final slip values
+ value-names = left-lateral-slip fault-opening
+ value-units = m m
+
+ // Number of locations where physical properties are given
+ // 1 location -> data is uniform -> data dimension is 0
+ num-locs = 1
+ data-dim = 0
+
+ // Problem is in 2-D
+ space-dim = 2
+
+ // Coordinate system (2-D with coordinates in meters)
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+// Columns are x, y, left-lateral-slip, fault-opening
+0.0 0.0 1.0 0.0
Added: short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_sliptime.spatialdb
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_sliptime.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/RateStateAgeing/quad4/shearwave_sliptime.spatialdb 2010-03-10 01:25:11 UTC (rev 16399)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ // number of values for slip initiation time
+ num-values = 1
+
+ // Names and units of slip initiation time
+ value-names = slip-time
+ value-units = s
+
+ // Number of locations where physical properties are given
+ // 1 location -> data is uniform -> data dimension is 0
+ num-locs = 1
+ data-dim = 0
+
+ // Problem is in 2-D
+ space-dim = 2
+
+ // Coordinate system (2-D with coordinates in meters)
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+// Columns are x, y, slip-time
+0.0 0.0 0.1
More information about the CIG-COMMITS
mailing list