kindred/create
1
1
Fork 0
Code Issues 37 Pull Requests Actions Packages Projects 9 Releases 2 Wiki Activity

FEM: results, move vonMises from result reader to result tools and use new result stress float lists

Browse Source
This commit is contained in:
Bernd Hahnebach
2019-03-15 21:06:30 +01:00
committed by Yorik van Havre
parent 3f58cebca4
commit d95e781d8f
3 changed files with 35 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
src/Mod/Fem/feminout/importCcxFrdResults.py
View File

@@ -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 = (

20
src/Mod/Fem/feminout/importToolsFem.py
View File

@@ -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]],

35
src/Mod/Fem/femresult/resulttools.py
View File

@@ -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: