Changeset 20966

Timestamp:

07/20/16 16:31:44 (9 years ago)

Author:

felicity

Message:

BUG: effective strain rate was not calculated correctly as it was ignoring eps_zz = -eps_xx

File:

• issm/trunk-jpl/src/c/classes/Materials/Matice.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/classes/Materials/Matice.cpp

@@ -644 +644 @@
         }
         else{
-                /* eps_eff^2 = 1/2 (exx^2 + 2*exy^2 )*/
+                /* eps_eff^2 = 1/2 (2*exx^2 + 2*exy^2 ) (since eps_zz = - eps_xx)*/
                 element->StrainRateHO2dvertical(&epsilon2d[0],xyz_list,gauss,vx_input,vy_input);
-                eps_eff = 1./sqrt(2.)*sqrt(epsilon2d[0]*epsilon2d[0] + 2.*epsilon2d[1]*epsilon2d[1]);
+                eps_eff = 1./sqrt(2.)*sqrt(2*epsilon2d[0]*epsilon2d[0] + 2*epsilon2d[1]*epsilon2d[1]);
         }
 
@@ -732 +732 @@
         }
         else{
-                /* eps_eff^2 = 1/2 exx^2*/
+                /* eps_eff^2 = exx^2*/
                 element->StrainRateSSA1d(&epsilon1d,xyz_list,gauss,vx_input);
-                eps_eff = sqrt(epsilon1d*epsilon1d/2.);
+                eps_eff = fabs(epsilon1d);
         }