From e0c0e8de005127f61cd377fa1fc88505987f3a24 Mon Sep 17 00:00:00 2001 From: Bernd Hahnebach Date: Mon, 2 Jul 2018 09:02:22 +0200 Subject: [PATCH] FEM: principal stress calculation, fix for direction of principal stress --- src/Mod/Fem/feminout/importToolsFem.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) 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()