diff --git a/src/Mod/Fem/feminout/importToolsFem.py b/src/Mod/Fem/feminout/importToolsFem.py
index 7efe97ffd9..7a62bab029 100644
--- 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()