[cig-commits] r11594 - in long/3D/Gale/trunk: . input/benchmarks

[walter at geodynamics.org]

Author: walter
Date: 2008-03-26 13:57:46 -0700 (Wed, 26 Mar 2008)
New Revision: 11594

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/input/benchmarks/shortening.xml
Log:
 r2047 at earth:  boo | 2008-03-20 12:21:34 -0700
 Add a viscous corner to shortening.xml.  Still does not work



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2046
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2047

Modified: long/3D/Gale/trunk/input/benchmarks/shortening.xml
===================================================================
--- long/3D/Gale/trunk/input/benchmarks/shortening.xml	2008-03-26 20:57:41 UTC (rev 11593)
+++ long/3D/Gale/trunk/input/benchmarks/shortening.xml	2008-03-26 20:57:46 UTC (rev 11594)
@@ -131,6 +131,13 @@
       <param name="Mesh">mesh-linear</param>
       <param name="IC">stressICs</param>
     </struct>
+    <struct name="ViscosityField">
+      <param name="Type">ViscosityField</param>
+      <param name="Context">context</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="Mesh">mesh-linear</param>
+      <param name="ConstitutiveMatrix">constitutiveMatrix</param>
+    </struct> 
     <struct name="cellLayout">
       <param name="Type">SingleCellLayout</param>
     </struct>
@@ -314,19 +321,33 @@
       <param name="startZ">minZ</param>
       <param name="endZ">maxZ</param>
     </struct>
+
+    <struct name="cornerShape">
+      <param name="Type">Box</param>
+      <param name="startX">0.39</param>
+      <param name="endX">maxX</param>
+      <param name="startY">0.0</param>
+      <param name="endY">0.01</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+
     <struct name="sandShape">
       <param name="Type">Intersection</param>
       <list name="shapes">
         <param>backgroundShape</param>
         <param>!beadsShape</param>
+        <param>!cornerShape</param>
       </list>
     </struct>
+
+
     <struct name="sandStrainWeakening">
       <param name="Type">StrainWeakening</param>
       <param name="TimeIntegrator">timeIntegrator</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="initialSofteningStrain">0.0</param>
-      <param name="finalSofteningStrain">0.01</param>
+      <param name="initialSofteningStrain">0.5</param>
+      <param name="finalSofteningStrain">1.0</param>
       <param name="initialDamageFraction">0.0</param>
       <param name="initialDamageWavenumber">0.5</param>
       <param name="initialDamageFactor">0.5</param>
@@ -341,10 +362,10 @@
       <param name="StrainWeakening">sandStrainWeakening</param>
       <param name="cohesion">10.0</param>
       <param name="cohesionAfterSoftening">0.1</param>
-      <param name="frictionCoefficient">0.726542528005</param>
-      <param name="frictionCoefficientAfterSoftening">0.600860619027</param>
-      <param name="minimumYieldStress">0.00001</param>
+      <param name="frictionCoefficient">0.404026225835</param>
+      <param name="frictionCoefficientAfterSoftening">0.363970234266</param>
     </struct>
+
     <struct name="sandViscosity">
       <param name="Type">MaterialViscosity</param>
       <param name="eta0">100000.0</param>
@@ -371,12 +392,15 @@
       <param name="cohesionAfterSoftening">0.1</param>
       <param name="frictionCoefficient">0.404026225835</param>
       <param name="frictionCoefficientAfterSoftening">0.363970234266</param>
-      <param name="minimumYieldStress">0.00001</param>
     </struct>
     <struct name="beadsViscosity">
       <param name="Type">MaterialViscosity</param>
       <param name="eta0">100000.0</param>
     </struct>
+    <struct name="cornerViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">1.0</param>
+    </struct>
     <struct name="storeViscosity">
       <param name="Type">StoreVisc</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
@@ -407,6 +431,16 @@
         <param>storeStress</param>
       </list>
     </struct>
+    <struct name="corner">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">cornerShape</param>
+      <param name="density">1560</param>
+      <list name="Rheology">
+        <param>cornerViscosity</param>
+        <param>storeViscosity</param>
+        <param>storeStress</param>
+      </list>
+    </struct>
     <struct name="surfaceAdaptor">
       <param name="Type">SurfaceAdaptor</param>
       <param name="mesh">mesh-linear</param>
@@ -581,7 +615,7 @@
           <struct>
             <param name="name">vx</param>
             <param name="type">double</param>
-            <param name="value">0</param>
+            <param name="value">0.0</param>
           </struct>
         </list>
       </struct>
@@ -604,10 +638,16 @@
     </list>
   </struct>
 
+  <param name="StepFunctionLowerOffset">0.0</param>
+  <param name="StepFunctionUpperOffset">0.035</param>
+  <param name="StepFunctionValue">-1.0</param>
+  <param name="StepFunctionDim">0</param>
+  <param name="StepFunctionLessThan">False</param>
+
   <param name="checkpointEvery">1</param>
   <param name="dtFactor">1.0</param>
-  <param name="journal.info">True</param>
-  <param name="journal.debug">True</param>
-  <param name="journal-level.info">2</param>
-  <param name="journal-level.debug">2</param>
+<!--   <param name="journal.info">True</param> -->
+<!--   <param name="journal.debug">True</param> -->
+<!--   <param name="journal-level.info">2</param> -->
+<!--   <param name="journal-level.debug">2</param> -->
 </StGermainData>