[cig-commits] r19816 - short/3D/PyLith/trunk/examples/3d/hex8

[brad at geodynamics.org]

Author: brad
Date: 2012-03-19 12:46:24 -0700 (Mon, 19 Mar 2012)
New Revision: 19816

Modified:
   short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg
Log:
Allow tensile yield in Drucker-Prager example.

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg	2012-03-19 19:45:52 UTC (rev 19815)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg	2012-03-19 19:46:24 UTC (rev 19816)
@@ -57,9 +57,9 @@
 # time stepping.
 time_step = pylith.problems.TimeStepAdapt
 
+[pylithapp.timedependent.implicit.time_step]
 # Change the total simulation time to 700 years, and set the maximum time
 # step size to 5 years.
-[pylithapp.timedependent.implicit.time_step]
 total_time = 700.0*year
 max_dt = 5.0*year
 stability_factor = 1.0 ; use time step equal to stable value from materials
@@ -83,13 +83,19 @@
 db_properties.db_A = spatialdata.spatialdb.SimpleDB
 db_properties.db_B = spatialdata.spatialdb.SimpleDB
 
+# Extend yield surface into tensile region (cone becomes line) to
+# allow tensile yielding. In most real applications, gravitational
+# body forces will keep the material in compression so the default
+# value (False) may be used.
+allow_tensile_yield = True
+
 # Provide the values to be obtained from each database and the database
 # name.
 [pylithapp.timedependent.materials.lower_crust.db_properties]
 values_A = [density,vs,vp]   ; Elastic properties.
 db_A.label = Elastic properties
 db_A.iohandler.filename = spatialdb/mat_elastic.spatialdb
-values_B = [friction-angle,cohesion,dilatation-angle]   ; Drucker-Prager properties.
+values_B = [friction-angle,cohesion,dilatation-angle]
 db_B.label = Drucker-Prager properties
 db_B.iohandler.filename = spatialdb/mat_druckerprager.spatialdb