diff --git a/src/Mod/Fem/feminout/importCcxFrdResults.py b/src/Mod/Fem/feminout/importCcxFrdResults.py
index b62fa1fb11..a6572232c2 100644
--- a/src/Mod/Fem/feminout/importCcxFrdResults.py
+++ b/src/Mod/Fem/feminout/importCcxFrdResults.py
@@ -104,6 +104,7 @@ def importFrd(filename, analysis=None, result_name_prefix=None):
                     # compact result object, workaround for bug 2873, https://www.freecadweb.org/tracker/view.php?id=2873
                     res_obj = restools.compact_result(res_obj)
                 res_obj = restools.add_disp_apps(res_obj)  # fill DisplacementLengths
+                res_obj = restools.add_von_mises(res_obj)  # fill StressValues
                 res_obj = restools.fill_femresult_stats(res_obj)  # fill Stats
         else:
             error_message = (
diff --git a/src/Mod/Fem/feminout/importToolsFem.py b/src/Mod/Fem/feminout/importToolsFem.py
index fab9cc3f3a..6bb6591bbd 100644
--- a/src/Mod/Fem/feminout/importToolsFem.py
+++ b/src/Mod/Fem/feminout/importToolsFem.py
@@ -314,27 +314,23 @@ def fill_femresult_mechanical(res_obj, result_set):
         if 'stress' in result_set:
             stress = result_set['stress']
             if len(stress) > 0:
-                mstress = []
                 prinstress1 = []
                 prinstress2 = []
                 prinstress3 = []
                 shearstress = []
                 for i in stress.values():
-                    mstress.append(calculate_von_mises(i))
                     prin1, prin2, prin3, shear = calculate_principal_stress(i)
                     prinstress1.append(prin1)
                     prinstress2.append(prin2)
                     prinstress3.append(prin3)
                     shearstress.append(shear)
                 if eigenmode_number > 0:
-                    res_obj.StressValues = list(map((lambda x: x * scale), mstress))
                     res_obj.PrincipalMax = list(map((lambda x: x * scale), prinstress1))
                     res_obj.PrincipalMed = list(map((lambda x: x * scale), prinstress2))
                     res_obj.PrincipalMin = list(map((lambda x: x * scale), prinstress3))
                     res_obj.MaxShear = list(map((lambda x: x * scale), shearstress))
                     res_obj.Eigenmode = eigenmode_number
                 else:
-                    res_obj.StressValues = mstress
                     res_obj.PrincipalMax = prinstress1
                     res_obj.PrincipalMed = prinstress2
                     res_obj.PrincipalMin = prinstress3
@@ -391,22 +387,6 @@ def fill_femresult_mechanical(res_obj, result_set):
 
 
 # helper
-def calculate_von_mises(i):
-    # Von mises stress (http://en.wikipedia.org/wiki/Von_Mises_yield_criterion)
-    s11 = i[0]
-    s22 = i[1]
-    s33 = i[2]
-    s12 = i[3]
-    s31 = i[4]
-    s23 = i[5]
-    s11s22 = pow(s11 - s22, 2)
-    s22s33 = pow(s22 - s33, 2)
-    s33s11 = pow(s33 - s11, 2)
-    s12s23s31 = 6 * (pow(s12, 2) + pow(s23, 2) + pow(s31, 2))
-    vm_stress = sqrt(0.5 * (s11s22 + s22s33 + s33s11 + s12s23s31))
-    return vm_stress
-
-
 def calculate_principal_stress(i):
     sigma = np.array([[i[0], i[3], i[4]],
                       [i[3], i[1], i[5]],
diff --git a/src/Mod/Fem/femresult/resulttools.py b/src/Mod/Fem/femresult/resulttools.py
index 427c68bc76..784ed27a7e 100644
--- a/src/Mod/Fem/femresult/resulttools.py
+++ b/src/Mod/Fem/femresult/resulttools.py
@@ -29,7 +29,7 @@ __url__ = "http://www.freecadweb.org"
 
 import FreeCAD
 import femtools.femutils as femutils
-from math import sqrt
+from math import pow, sqrt
 
 
 ## Removes all result objects and result meshes from an analysis group
@@ -262,6 +262,23 @@ def add_disp_apps(res_obj):
     return res_obj
 
 
+def add_von_mises(res_obj):
+    mstress = []
+    iterator = zip(
+        res_obj.NodeStressXX,
+        res_obj.NodeStressYY,
+        res_obj.NodeStressZZ,
+        res_obj.NodeStressXY,
+        res_obj.NodeStressXZ,
+        res_obj.NodeStressYZ
+    )
+    for Sxx, Syy, Szz, Sxy, Sxz, Syz in iterator:
+        mstress.append(calculate_von_mises((Sxx, Syy, Szz, Sxy, Sxz, Syz)))
+    res_obj.StressValues = mstress
+    FreeCAD.Console.PrintMessage('Added StressValues (von Mises).\n')
+    return res_obj
+
+
 def compact_result(res_obj):
     '''
     compacts result.Mesh and appropriate result.NodeNumbers
@@ -287,6 +304,22 @@ def compact_result(res_obj):
     return res_obj
 
 
+def calculate_von_mises(i):
+    # Von mises stress (http://en.wikipedia.org/wiki/Von_Mises_yield_criterion)
+    s11 = i[0]
+    s22 = i[1]
+    s33 = i[2]
+    s12 = i[3]
+    s31 = i[4]
+    s23 = i[5]
+    s11s22 = pow(s11 - s22, 2)
+    s22s33 = pow(s22 - s33, 2)
+    s33s11 = pow(s33 - s11, 2)
+    s12s23s31 = 6 * (pow(s12, 2) + pow(s23, 2) + pow(s31, 2))
+    vm_stress = sqrt(0.5 * (s11s22 + s22s33 + s33s11 + s12s23s31))
+    return vm_stress
+
+
 def calculate_disp_abs(displacements):
     disp_abs = []
     for d in displacements: