Fix trace in DeviatoricStress equation.

pypr · Feb 18, 2020 · 381f4e6 · 381f4e6
1 parent e3d5e10
commit 381f4e6
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/pysph/sph/solid_mech/basic.py b/pysph/sph/solid_mech/basic.py
@@ -472,7 +472,7 @@ def loop(self, d_idx, d_s00, d_s01, d_s02, d_s11, d_s12, d_s22, d_v00,
         omega22 = 0.0
 
         tmp = 2.0 * d_G[0]
-        trace = 1.0 / 3.0 * (eps00 + eps11)
+        trace = 1.0 / 3.0 * (eps00 + eps11 + eps22)
 
         # S_00
         d_as00[d_idx] = tmp*( eps00 - trace ) + \