FEM: principal stress calculation, fix for direction of principal stress

FreeCAD · Jul 3, 2018 · 0f8566e · 0f8566e
1 parent 1bd2463
commit 0f8566e
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/Mod/Fem/feminout/importToolsFem.py b/src/Mod/Fem/feminout/importToolsFem.py
@@ -426,9 +426,9 @@ def calculate_von_mises(i):
 
 
 def calculate_principal_stress(i):
-    sigma = np.array([[i[0], i[3], i[4]],
-                      [i[3], i[1], i[5]],
-                      [i[4], i[5], i[2]]])
+    sigma = np.array([[i[0], i[3], i[5]],
+                      [i[3], i[1], i[4]],
+                      [i[5], i[4], i[2]]])  # https://forum.freecadweb.org/viewtopic.php?f=18&t=24637&start=10#p240408
     # compute principal stresses
     eigvals = list(np.linalg.eigvalsh(sigma))
     eigvals.sort()