[cig-commits] r15251 - short/3D/PyLith/trunk/examples/3d/hex8
brad at geodynamics.org
brad at geodynamics.org
Sun Jun 14 21:36:41 PDT 2009
Author: brad
Date: 2009-06-14 21:36:40 -0700 (Sun, 14 Jun 2009)
New Revision: 15251
Added:
short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb
Removed:
short/3D/PyLith/trunk/examples/3d/hex8/initial_state.spatialdb
short/3D/PyLith/trunk/examples/3d/hex8/planestrain.cfg
short/3D/PyLith/trunk/examples/3d/hex8/savageprescott-powerlaw.cfg
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg
short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg
short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg
Log:
Cleaned up parameter settings.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/dislocation.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -28,33 +28,6 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Set materials to an array of 2 materials:
-# 'elastic' and 'maxwell'
-materials = [elastic,maxwell]
-
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.maxwell = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-[pylithapp.timedependent.materials.maxwell]
-label = Viscoelastic material
-id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
@@ -123,10 +96,10 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = dislocation-statevars-elastic.vtk
+writer.filename = dislocation-elastic.vtk
-[pylithapp.timedependent.materials.maxwell.output]
+[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = dislocation-statevars-maxwell.vtk
+writer.filename = dislocation-viscoelastic.vtk
Modified: short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/gravity.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -30,33 +30,6 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Set materials to an array of 2 materials:
-# 'elastic' and 'maxwell'
-materials = [elastic,maxwell]
-
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.maxwell = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-[pylithapp.timedependent.materials.maxwell]
-label = Viscoelastic material
-id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
# Use the default FixedDOFDB, which has a uniform spatial distribution.
@@ -111,9 +84,9 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = gravity-statevars-elastic.vtk
+writer.filename = gravity-elastic.vtk
-[pylithapp.timedependent.materials.maxwell.output]
+[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = gravity-statevars-maxwell.vtk
+writer.filename = gravity-viscoelastic.vtk
Modified: short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/gravity_istress.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -18,9 +18,6 @@
# 'y_neg', and 'z_neg'.
bc = [x_pos,x_neg,y_pos,y_neg,z_neg]
-# Use gravity for this problem.
-use_gravity = true
-
[pylithapp.timedependent.implicit]
# Set the output to an array of 2 output managers.
# We will output the solution over the domain and the ground surface.
@@ -30,44 +27,13 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Set materials to an array of 2 materials:
-# 'elastic' and 'maxwell'
-materials = [elastic,maxwell]
-
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.maxwell = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-[pylithapp.timedependent.materials.maxwell]
-label = Viscoelastic material
-id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-# ----------------------------------------------------------------------
# initial stresses
# ----------------------------------------------------------------------
[pylithapp.timedependent.materials.elastic]
-use_initial_state = True
-initial_state_db.iohandler.filename = initial_state.spatialdb
-# initial_state_db.query_type = linear
+db_initial_stress.iohandler.filename = initial_stress.spatialdb
-[pylithapp.timedependent.materials.maxwell]
-use_initial_state = True
-initial_state_db.iohandler.filename = initial_state.spatialdb
-# initial_state_db.query_type = linear
+[pylithapp.timedependent.materials.viscoelastic]
+db_initial_stress.iohandler.filename = initial_stress.spatialdb
# ----------------------------------------------------------------------
# boundary conditions
@@ -124,12 +90,12 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = gravity_istress-statevars-elastic.vtk
+writer.filename = gravity_istress-elastic.vtk
-[pylithapp.timedependent.materials.maxwell.output]
+[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = gravity_istress-statevars-maxwell.vtk
+writer.filename = gravity_istress-viscoelastic.vtk
[pylithapp.petsc]
ksp_atol = 1.0e-12
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/initial_state.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/initial_state.spatialdb 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/initial_state.spatialdb 2009-06-15 04:36:40 UTC (rev 15251)
@@ -1,30 +0,0 @@
-// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
-//
-// This spatial database specifies the initial stresses for the model
-// The stresses increase linearly with depth.
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 6 // number of stress components
- value-names = stress_xx stress_yy stress_zz stress_xy stress_yz stress_xz // names of the stress components
- value-units = Pa Pa Pa Pa Pa Pa // units
- num-locs = 4 // number of locations: 1 value at the mid-depth for each
- // layer of elements.
- data-dim = 1
- space-dim = 3
- cs-data = cartesian {
- to-meters = 1.0
- space-dim = 3
- }
-}
-// Columns are
-// (1) x coordinate (m)
-// (2) y coordinate (m)
-// (3) z coordinate (m)
-// (4) density (kg/m^3)
-// (5) vs (m/s)
-// (6) vp (m/s)
-0.0 0.0 -3500.0 -3.064578125e+07 -3.064578125e+07 -8.58081875e+07 0.0 0.0 0.0
-0.0 0.0 -2500.0 -2.188984375e+07 -2.188984375e+07 -6.12915625e+07 0.0 0.0 0.0
-0.0 0.0 -1500.0 -1.313390625e+07 -1.313390625e+07 -3.67749375e+07 0.0 0.0 0.0
-0.0 0.0 -500.0 -4.37796875e+06 -4.37796875e+06 -1.22583125e+07 0.0 0.0 0.0
Copied: short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb (from rev 15250, short/3D/PyLith/trunk/examples/3d/hex8/initial_state.spatialdb)
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/hex8/initial_stress.spatialdb 2009-06-15 04:36:40 UTC (rev 15251)
@@ -0,0 +1,30 @@
+// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
+//
+// This spatial database specifies the initial stresses for the model
+// The stresses increase linearly with depth.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6 // number of stress components
+ value-names = stress_xx stress_yy stress_zz stress_xy stress_yz stress_xz // names of the stress components
+ value-units = Pa Pa Pa Pa Pa Pa // units
+ num-locs = 4 // number of locations: 1 value at the mid-depth for each
+ // layer of elements.
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are
+// (1) x coordinate (m)
+// (2) y coordinate (m)
+// (3) z coordinate (m)
+// (4) density (kg/m^3)
+// (5) vs (m/s)
+// (6) vp (m/s)
+0.0 0.0 -3500.0 -3.064578125e+07 -3.064578125e+07 -8.58081875e+07 0.0 0.0 0.0
+0.0 0.0 -2500.0 -2.188984375e+07 -2.188984375e+07 -6.12915625e+07 0.0 0.0 0.0
+0.0 0.0 -1500.0 -1.313390625e+07 -1.313390625e+07 -3.67749375e+07 0.0 0.0 0.0
+0.0 0.0 -500.0 -4.37796875e+06 -4.37796875e+06 -1.22583125e+07 0.0 0.0 0.0
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/planestrain.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/planestrain.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/planestrain.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -1,157 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the boundary
-# conditions. The general parameters are specificed in the pylithapp.cfg
-# file which PyLith reads by default.
-#
-# To run the simulation:
-# pylith hex8_1000m.cfg shearxy.cfg
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Change initial guess for solution from vector of zeros to previous
-# solution.
-formulation.solver.initial_guess_zero = False
-
-# Set bc to an array of 5 boundary conditions.
-bc = [x_pos,x_neg,y_pos,y_neg,z_neg]
-
-# Change time stepping algorithm from uniform time step, to adaptive
-# time stepping.
-implicit.time_step = pylith.problems.TimeStepAdapt
-
-[pylithapp.timedependent.implicit.time_step]
-# Define the total time for the simulation and the default time step size.
-# Overwrite values in pylithapp.cfg
-total_time = 500.0*year ; total time of simulation
-max_dt = 50.0*year ; time step
-stability_factor = 1.0 ; use time step equal to stable value from materials
-
-[pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
-# We will output the solution over the domain and the ground surface.
-output = [domain,subdomain]
-
-# Set subdomain component to OutputSolnSubset (subset of domain).
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Set materials to an array of 2 materials:
-# 'maxwell_top' and 'maxwell_bottom'
-materials = [maxwell_top,maxwell_bottom]
-
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.maxwell_top = pylith.materials.MaxwellIsotropic3D
-materials.maxwell_bottom = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.maxwell_top]
-label = Upper Maxwell material
-id = 1
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-[pylithapp.timedependent.materials.maxwell_bottom]
-label = Lower Maxwell material
-id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-# We change spatial database for initial value from FixedDOFDB (which has
-# a uniform spatial distribution) to SimpleDB (which uses a data file
-# to specify a spatial variation).
-
-# Set the parameters for Dirichlet boundary conditions applied on the
-# +x and -x faces of the box.
-
-# We fix the x degrees of freedom on the +x and -x faces, the y degrees of
-# freedom on the +y and -y faces, and fix the z degree of freedom on the
-# bottom face.
-
-# The label corresponds to the nodeset ID in CUBIT.
-
-[pylithapp.timedependent.bc.x_pos]
-fixed_dof = [0]
-label = 11
-db = spatialdata.spatialdb.SimpleDB
-db.label = Dirichlet BC on +x
-db.iohandler.filename = fixeddisp_axial.spatialdb
-
-[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0]
-label = 12
-db = spatialdata.spatialdb.SimpleDB
-db.label = Dirichlet BC on -x
-db.iohandler.filename = fixeddisp_axial.spatialdb
-
-[pylithapp.timedependent.bc.y_pos]
-fixed_dof = [1]
-label = 13
-db = spatialdata.spatialdb.SimpleDB
-db.label = Dirichlet BC on +y
-db.iohandler.filename = fixeddisp_axial.spatialdb
-
-[pylithapp.timedependent.bc.y_neg]
-fixed_dof = [1]
-label = 14
-db = spatialdata.spatialdb.SimpleDB
-db.label = Dirichlet BC on +y
-db.iohandler.filename = fixeddisp_axial.spatialdb
-
-[pylithapp.timedependent.bc.z_neg]
-fixed_dof = [2]
-label = 15
-db = spatialdata.spatialdb.SimpleDB
-db.label = Dirichlet BC on -z
-db.iohandler.filename = fixeddisp_axial.spatialdb
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain]
-skip = 0
-writer.filename = planestrain.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for VTK domain output of solution over ground surface.
-[pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
-skip = 0
-writer.filename = planestrain-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.maxwell_top.output]
-skip = 0
-cell_info_fields = [density,mu,lambda,maxwell_time]
-cell_data_fields = [total_strain,viscous_strain,stress]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = planestrain-statevars-maxwell_top.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.maxwell_bottom.output]
-skip = 0
-cell_info_fields = [density,mu,lambda,maxwell_time]
-cell_data_fields = [total_strain,viscous_strain,stress]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = planestrain-statevars-maxwell_bottom.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
Modified: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -40,6 +40,33 @@
total_time = 0.0*s ; total time of simulation
# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+# Set materials to an array of 2 materials:
+# 'elastic' and 'viscoelastic'
+materials = [elastic,viscoelastic]
+
+# Change from default material (elastic, isotropic, 3-D material) to
+# Maxwell viscoelastic, isotropic 3-D material.
+materials.viscoelastic = pylith.materials.MaxwellIsotropic3D
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = mat_elastic.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+[pylithapp.timedependent.materials.viscoelastic]
+label = Viscoelastic material
+id = 2
+db_properties.iohandler.filename = mat_maxwell.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
# Set the solver options.
Deleted: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott-powerlaw.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/savageprescott-powerlaw.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/savageprescott-powerlaw.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -1,194 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the boundary and
-# interface conditions. The general parameters are specificed in the
-# pylithapp.cfg file which PyLith reads by default.
-#
-# To run the simulation:
-# pylith savageprescott-powerlaw.cfg
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Set this up as a nonlinear problem.
-implicit.solver = pylith.problems.SolverNonlinear
-
-# Change initial guess for solution from vector of zeros to previous
-# solution.
-formulation.solver.initial_guess_zero = False
-
-# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
-bc = [x_pos,x_neg,z_neg]
-
-# Set interfaces to an array of 1 fault: 'fault'.
-interfaces = [fault]
-
-# Change time stepping algorithm from uniform time step, to adaptive
-# time stepping.
-implicit.time_step = pylith.problems.TimeStepAdapt
-
-[pylithapp.timedependent.implicit.time_step]
-# Define the total time for the simulation and the default time step size.
-# Overwrite values in pylithapp.cfg
-total_time = 700.0*year ; total time of simulation
-max_dt = 50.0*year ; time step
-stability_factor = 1.0 ; use time step equal to stable value from materials
-
-[pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
-# We will output the solution over the domain and the ground surface.
-output = [domain,subdomain]
-
-# Set subdomain component to OutputSolnSubset (subset of domain).
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Set materials to an array of 2 materials:
-# 'elastic' and 'powerlaw'
-materials = [elastic,powerlaw]
-
-# Change from default material (elastic, isotropic, 3-D material) to
-# power-law viscoelastic, isotropic 3-D material.
-materials.powerlaw = pylith.materials.PowerLaw3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-[pylithapp.timedependent.materials.powerlaw]
-label = Power-law material
-id = 2
-db_properties.iohandler.filename = mat_powerlaw.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-
-# We change spatial database for the rate of change from FixedDOFDB
-# (which has zero velocity) to UniformDB (which allows nonzero
-# velocities).
-
-# Set the parameters for the desired boundary conditions.
-#
-# The label corresponds to the nodeset ID in CUBIT.
-
-[pylithapp.timedependent.bc.x_pos]
-fixed_dof = [0, 1]
-label = 11
-db_initial.label = Dirichlet BC on +x
-db_rate = spatialdata.spatialdb.UniformDB
-db_rate.label = Dirichlet rate BC on +x
-db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
-db_rate.data = [0.0*cm/year,-1.0*cm/year,0.0*cm/year]
-
-[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0, 1]
-label = 12
-db_initial.label = Dirichlet BC on -x
-db_rate = spatialdata.spatialdb.UniformDB
-db_rate.label = Dirichlet rate BC on -x
-db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
-db_rate.data = [0.0*cm/year,+1.0*cm/year,0.0*cm/year]
-
-[pylithapp.timedependent.bc.z_neg]
-fixed_dof = [2]
-label = 16
-db_initial.label = Dirichlet BC on -z
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces]
-
-# Set the parameters for the fault interface conditions.
-
-[pylithapp.timedependent.interfaces.fault]
-label = 10
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-output.vertex_info_fields = [normal_dir,strike_dir,dip_dir,final_slip_creep,final_slip_one,final_slip_two,final_slip_three,slip_time_creep,slip_time_one,slip_time_two,slip_time_three]
-
-# Set earthquake ruptures
-eq_srcs = [creep,one,two,three]
-eq_srcs.creep.origin_time = 00.0*year
-eq_srcs.one.origin_time = 200.0*year
-eq_srcs.two.origin_time = 400.0*year
-eq_srcs.three.origin_time = 600.0*year
-
-# one
-[pylithapp.timedependent.interfaces.fault.eq_srcs.one.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# two
-[pylithapp.timedependent.interfaces.fault.eq_srcs.two.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# three
-[pylithapp.timedependent.interfaces.fault.eq_srcs.three.slip_function]
-slip.iohandler.filename = finalslip_rupture.spatialdb
-slip_time.iohandler.filename = sliptime.spatialdb
-
-# creep
-[pylithapp.timedependent.interfaces.fault.eq_srcs.creep]
-slip_function = pylith.faults.ConstRateSlipFn
-origin_time = 0.0*s
-slip_function.slip_rate.iohandler.filename = sliprate_creep.spatialdb
-slip_function.slip_time.iohandler.filename = sliptime.spatialdb
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain]
-skip = 9
-writer.filename = savageprescott-powerlaw.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for VTK domain output of solution over ground surface.
-[pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
-skip = 0
-writer.filename = savageprescott-powerlaw-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for vtk fault rupture output.
-[pylithapp.timedependent.interfaces.fault.output]
-skip = 0
-writer.filename = savageprescott-powerlaw-fault.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-skip = 9
-cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = savageprescott-powerlaw-statevars-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.powerlaw.output]
-cell_info_fields = [density,mu,lambda,viscosity_coeff,power_law_exponent]
-# cell_data_fields = [total_strain,stress,viscous_strain]
-skip = 9
-cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = savageprescott-powerlaw-statevars-powerlaw.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
Modified: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -13,10 +13,6 @@
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-# Change initial guess for solution from vector of zeros to previous
-# solution.
-formulation.solver.initial_guess_zero = False
-
# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
bc = [x_pos,x_neg,z_neg]
@@ -42,33 +38,25 @@
# Set subdomain component to OutputSolnSubset (subset of domain).
output.subdomain = pylith.meshio.OutputSolnSubset
+[pylithapp.timedependent.materials.viscoelastic]
+output.cell_info_fields = [density,mu,lambda,maxwell_time] ; Maxwell model
+
# ----------------------------------------------------------------------
-# materials
+# Uncomment the lines below to switch to the power-law rheology and
+# nonlinear solver.
# ----------------------------------------------------------------------
+
[pylithapp.timedependent]
+# Uncomment next line to use the nonlinear solver.
+#implicit.solver = pylith.problems.SolverNonlinear
-# Set materials to an array of 2 materials:
-# 'elastic' and 'maxwell'
-materials = [elastic,maxwell]
+# Change from the Maxwell viscoelastic, isotropic 3-D material to power-law.
+#materials.viscoelastic = pylith.materials.PowerLaw3D
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.maxwell = pylith.materials.MaxwellIsotropic3D
+[pylithapp.timedependent.materials.viscoelastic]
+#db_properties.iohandler.filename = mat_powerlaw.spatialdb
+#output.cell_info_fields = [density,mu,lambda,viscosity_coeff,power_law_exponent]
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-[pylithapp.timedependent.materials.maxwell]
-label = Maxwell viscoelastic material
-id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
@@ -176,14 +164,13 @@
[pylithapp.timedependent.materials.elastic.output]
skip = 9
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = savageprescott-statevars-elastic.vtk
+writer.filename = savageprescott-elastic.vtk
writer.time_format = %04.0f
writer.time_constant = 1.0*year
-[pylithapp.timedependent.materials.maxwell.output]
-cell_info_fields = [density,mu,lambda,maxwell_time]
+[pylithapp.timedependent.materials.viscoelastic.output]
skip = 9
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = savageprescott-statevars-maxwell.vtk
+writer.filename = savageprescott-viscoelastic.vtk
writer.time_format = %04.0f
writer.time_constant = 1.0*year
Modified: short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg 2009-06-15 02:52:09 UTC (rev 15250)
+++ short/3D/PyLith/trunk/examples/3d/hex8/shearxy.cfg 2009-06-15 04:36:40 UTC (rev 15251)
@@ -26,33 +26,6 @@
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Set materials to an array of 2 materials:
-# 'elastic' and 'maxwell'
-materials = [elastic,maxwell]
-
-# Change from default material (elastic, isotropic, 3-D material) to
-# Maxwell viscoelastic, isotropic 3-D material.
-materials.maxwell = pylith.materials.MaxwellIsotropic3D
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = mat_elastic.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-[pylithapp.timedependent.materials.maxwell]
-label = Viscoelastic material
-id = 2
-db_properties.iohandler.filename = mat_maxwell.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-
-# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
# We change spatial database for initial value from FixedDOFDB (which has
@@ -103,9 +76,9 @@
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = shearxy-statevars-elastic.vtk
+writer.filename = shearxy-elastic.vtk
-[pylithapp.timedependent.materials.maxwell.output]
+[pylithapp.timedependent.materials.viscoelastic.output]
cell_info_fields = [density,mu,lambda,maxwell_time]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = shearxy-statevars-maxwell.vtk
+writer.filename = shearxy-viscoelastic.vtk
More information about the CIG-COMMITS
mailing list