[cig-commits] r14178 - long/3D/SNAC/trunk/Snac/plugins/plastic

[cstark at geodynamics.org]

Author: cstark
Date: 2009-03-01 15:01:19 -0800 (Sun, 01 Mar 2009)
New Revision: 14178

Modified:
   long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c
Log:
Major/minor bug fix in Constitutive.c

Corrected recent (14173) change to handling of cohesion etc estimation from piecewise linear function of plastic strain.
Terminal values for extreme plastic strain were set for penultimate piecewise node, not the final node.
Hardening should also have been set to zero.  


Modified: long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c
===================================================================
--- long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c	2009-03-01 18:52:08 UTC (rev 14177)
+++ long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c	2009-03-01 23:01:19 UTC (rev 14178)
@@ -134,7 +134,10 @@
 #endif
 					}
 			}
-			/*ccccc*/
+			/*
+			*  Calculate material props (friction angle, cohesion, etc) from piecewise-linear
+			*    functions defined in terms of plastic strain.
+			*/
                         for( i = 0; i < material->nsegments; i++ ) {
 				const double pl1 = material->plstrain[i];
 				const double pl2 = material->plstrain[i+1];
@@ -154,10 +157,10 @@
 				     *    has pl strain larger than piece  -  then set physical params to maxima for this piece
 				     *  Purpose:  to prevent crazy behavior for extreme plastic strains
 				     */
-				    frictionAngle = material->frictionAngle[i];
-				    dilationAngle = material->dilationAngle[i];
-				    cohesion = material->cohesion[i];
-				    hardening = (material->cohesion[i+1] - material->cohesion[i]) / (pl2 - pl1);
+				    frictionAngle = material->frictionAngle[material->nsegments];
+				    dilationAngle = material->dilationAngle[material->nsegments];
+				    cohesion = material->cohesion[material->nsegments];
+				    hardening = 0.0;
 				}
                         }