[cig-commits] r16423 - short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4
brad at geodynamics.org
brad at geodynamics.org
Mon Mar 15 16:37:21 PDT 2010
Author: brad
Date: 2010-03-15 16:37:21 -0700 (Mon, 15 Mar 2010)
New Revision: 16423
Modified:
short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg
Log:
Fixed errors in settings for lumped, dynamic example.
Modified: short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg 2010-03-15 22:05:22 UTC (rev 16422)
+++ short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg 2010-03-15 23:37:21 UTC (rev 16423)
@@ -5,12 +5,17 @@
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
-formulation = pylith.problems.Explicit
-#formulation = pylith.problems.ExplicitLumped
-formulation.solver = pylith.problems.SolverNonlinear
+[pylithapp.timedependent]
+#formulation = pylith.problems.Explicit
+#formulation.solver = pylith.problems.SolverNonlinear
+formulation = pylith.problems.ExplicitLumped
+
# Nondimensionalize problem using wave propagation parameters.
normalizer = spatialdata.units.NondimElasticDynamic
+normalizer.mass_density = 1000.0*kg/m**3
+normalizer.shear_wave_speed = 1.0*m/s
+normalizer.wave_period = 1.0*s
bc = [x_pos,x_neg,all_nofault]
@@ -57,7 +62,7 @@
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2
-[pylithapp.timedependent.bc.y_pos]
+[pylithapp.timedependent.bc.all_nofault]
# Dirichlet BC on all vertices except fault vertices
# Pin x and z degrees of freedom
@@ -92,7 +97,7 @@
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]
+db_initial_tractions.data = [-6001.0*MPa, -10000.0*MPa]
friction = pylith.friction.StaticFriction
More information about the CIG-COMMITS
mailing list