[cig-commits] r7670 - in long/3D/Gale/trunk: . input
walter at geodynamics.org
walter at geodynamics.org
Sat Jul 14 01:50:57 PDT 2007
Author: walter
Date: 2007-07-14 01:50:56 -0700 (Sat, 14 Jul 2007)
New Revision: 7670
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/extension.xml
long/3D/Gale/trunk/input/shortening.xml
Log:
r1876 at earth: boo | 2007-07-14 01:46:53 -0700
Make extension and shortening work with the new friction. Shortening is a bit broken, getting hung up on the second time step.
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1871
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1876
Modified: long/3D/Gale/trunk/input/extension.xml
===================================================================
--- long/3D/Gale/trunk/input/extension.xml 2007-07-14 00:54:17 UTC (rev 7669)
+++ long/3D/Gale/trunk/input/extension.xml 2007-07-14 08:50:56 UTC (rev 7670)
@@ -122,7 +122,6 @@
<param name="DofLayout">pressureDofLayout</param>
<param name="LinkedDofInfo">pressureLinkedDofs</param>
</struct>
-
<struct name="StressField">
<param name="Type">StressField</param>
<param name="StrainRateField">StrainRateField</param>
@@ -132,7 +131,6 @@
<param name="Mesh">mesh-linear</param>
<param name="IC">stressICs</param>
</struct>
-
<struct name="cellLayout">
<param name="Type">SingleCellLayout</param>
</struct>
@@ -366,7 +364,7 @@
<struct name="crust">
<param name="Type">RheologyMaterial</param>
<param name="Shape">crustShape</param>
- <param name="density">100.0</param>
+ <param name="density">1.0</param>
<list name="Rheology">
<param>crustViscosity</param>
<param>yielding</param>
@@ -377,7 +375,7 @@
<struct name="pdms">
<param name="Type">RheologyMaterial</param>
<param name="Shape">pdmsShape</param>
- <param name="density">100.0</param>
+ <param name="density">1.0</param>
<list name="Rheology">
<param>pdmsViscosity</param>
<param>storeViscosity</param>
@@ -553,6 +551,9 @@
<param name="type">FrictionVC</param>
<param name="wall">bottom</param>
<param name="StaticFriction">0.1</param>
+ <param name="includeUpperX">False</param>
+ <param name="StressField">StressField</param>
+ <param name="PressureField">PressureField</param>
<list name="variables">
<struct>
<param name="name">vx</param>
Modified: long/3D/Gale/trunk/input/shortening.xml
===================================================================
--- long/3D/Gale/trunk/input/shortening.xml 2007-07-14 00:54:17 UTC (rev 7669)
+++ long/3D/Gale/trunk/input/shortening.xml 2007-07-14 08:50:56 UTC (rev 7670)
@@ -122,6 +122,15 @@
<param name="DofLayout">pressureDofLayout</param>
<param name="LinkedDofInfo">pressureLinkedDofs</param>
</struct>
+ <struct name="StressField">
+ <param name="Type">StressField</param>
+ <param name="StrainRateField">StrainRateField</param>
+ <param name="Context">context</param>
+ <param name="ConstitutiveMatrix">constitutiveMatrix</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="Mesh">mesh-linear</param>
+ <param name="IC">stressICs</param>
+ </struct>
<struct name="cellLayout">
<param name="Type">SingleCellLayout</param>
</struct>
@@ -287,6 +296,12 @@
<param name="nonLinearTolerance">nonLinearTolerance</param>
<param name="makeConvergenceFile">false</param>
</struct>
+ <struct name="buoyancyForceTerm">
+ <param name="Type">BuoyancyForceTerm</param>
+ <param name="ForceVector">mom_force</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="gravity">9.8</param>
+ </struct>
<struct name="backgroundShape">
<param name="Type">Everywhere</param>
</struct>
@@ -299,43 +314,68 @@
<param name="startZ">minZ</param>
<param name="endZ">maxZ</param>
</struct>
- <struct name="crustShape">
+ <struct name="sandShape">
<param name="Type">Intersection</param>
<list name="shapes">
<param>backgroundShape</param>
<param>!beadsShape</param>
</list>
</struct>
- <struct name="crustViscosity">
+ <struct name="sandStrainWeakening">
+ <param name="Type">StrainWeakening</param>
+ <param name="TimeIntegrator">timeIntegrator</param>
+ <param name="MaterialPointsSwarm">materialSwarm</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>
+ <param name="healingRate">0.0</param>
+ </struct>
+ <struct name="sandYielding">
+ <param name="Type">MohrCoulomb</param>
+ <param name="PressureField">PressureField</param>
+ <param name="VelocityGradientsField">VelocityGradientsField</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ <param name="Context">context</param>
+ <param name="StrainWeakening">sandStrainWeakening</param>
+ <param name="cohesion">10.0</param>
+ <param name="cohesionAfterSoftening">1.0</param>
+ <param name="frictionCoefficient">0.726542528005</param>
+ <param name="frictionCoefficientAfterSoftening">0.600860619027</param>
+ <param name="minimumYieldStress">0.00001</param>
+ </struct>
+ <struct name="sandViscosity">
<param name="Type">MaterialViscosity</param>
- <param name="eta0">10.0</param>
+ <param name="eta0">1.0e2</param>
</struct>
- <struct name="strainWeakening">
+ <struct name="beadsStrainWeakening">
<param name="Type">StrainWeakening</param>
<param name="TimeIntegrator">timeIntegrator</param>
<param name="MaterialPointsSwarm">materialSwarm</param>
- <param name="softeningStrain">0.1</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>
<param name="healingRate">0.0</param>
</struct>
- <struct name="yielding">
+ <struct name="beadsYielding">
<param name="Type">MohrCoulomb</param>
<param name="PressureField">PressureField</param>
<param name="VelocityGradientsField">VelocityGradientsField</param>
<param name="MaterialPointsSwarm">materialSwarm</param>
<param name="Context">context</param>
- <param name="StrainWeakening">strainWeakening</param>
+ <param name="StrainWeakening">beadsStrainWeakening</param>
<param name="cohesion">10.0</param>
- <param name="cohesionAfterSoftening">0.0001</param>
- <param name="frictionCoefficient">0.0</param>
- <param name="frictionCoefficientAfterSoftening">0.0</param>
+ <param name="cohesionAfterSoftening">1.0</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">1.0</param>
+ <param name="eta0">1.0e2</param>
</struct>
<struct name="storeViscosity">
<param name="Type">StoreVisc</param>
@@ -345,12 +385,13 @@
<param name="Type">StoreStress</param>
<param name="MaterialPointsSwarm">materialSwarm</param>
</struct>
- <struct name="crust">
+ <struct name="sand">
<param name="Type">RheologyMaterial</param>
- <param name="Shape">crustShape</param>
+ <param name="Shape">sandShape</param>
+ <param name="density">1560</param>
<list name="Rheology">
- <param>crustViscosity</param>
- <param>yielding</param>
+ <param>sandViscosity</param>
+ <param>sandYielding</param>
<param>storeViscosity</param>
<param>storeStress</param>
</list>
@@ -358,8 +399,10 @@
<struct name="beads">
<param name="Type">RheologyMaterial</param>
<param name="Shape">beadsShape</param>
+ <param name="density">1480</param>
<list name="Rheology">
<param>beadsViscosity</param>
+ <param>beadsYielding</param>
<param>storeViscosity</param>
<param>storeStress</param>
</list>
@@ -372,7 +415,7 @@
<list name="origin">
<param>0.0</param>
</list>
- <param name="offset">1.5</param>
+ <param name="offset">.3</param>
<param name="gradient">0.1763269807</param>
</struct>
<struct name="escapedRoutine">
@@ -405,7 +448,7 @@
<param>Underworld_DumpSwarm</param>
<param>Underworld_VTKOutput</param>
</list>
- <param name="maxTimeSteps">30</param>
+ <param name="maxTimeSteps">500</param>
<param name="outputEvery">1</param>
<param name="dumpEvery">1</param>
<param name="outputPath">./output.shortening</param>
@@ -414,10 +457,10 @@
<param name="minX">0.0f</param>
<param name="minY">0.0f</param>
<param name="minZ">0.0f</param>
- <param name="maxX">2.0f</param>
- <param name="maxY">0.35f</param>
+ <param name="maxX">0.40f</param>
+ <param name="maxY">0.035f</param>
<param name="maxZ">1.0f</param>
- <param name="elementResI">32</param>
+ <param name="elementResI">64</param>
<param name="elementResJ">16</param>
<param name="elementResK">4</param>
<param name="allowUnbalancing">True</param>
@@ -514,7 +557,39 @@
</struct>
</list>
</struct>
+ <struct>
+ <param name="type">FrictionVC</param>
+ <param name="wall">bottom</param>
+ <param name="StaticFriction">0.34432761329</param>
+ <param name="includeUpperX">False</param>
+ <param name="StressField">StressField</param>
+ <param name="PressureField">PressureField</param>
+ <list name="variables">
+ <struct>
+ <param name="name">vx</param>
+ <param name="type">double</param>
+ <param name="value">0</param>
+ </struct>
+ </list>
+ </struct>
</list>
</struct>
+
+ <struct name="stressICs">
+ <param name="type">CompositeVC</param>
+ <list name="vcList">
+ <struct>
+ <param name="type">AllNodesVC</param>
+ <list name="variables">
+ <struct>
+ <param name="name">stress</param>
+ <param name="type">double</param>
+ <param name="value">0.0</param>
+ </struct>
+ </list>
+ </struct>
+ </list>
+ </struct>
+
<param name="checkpointEvery">1</param>
</StGermainData>
More information about the cig-commits
mailing list