Principal Stress Plots
Moderator: bernd
Forum rules
and Helpful information for the FEM forum
and Helpful information for the FEM forum
Principal Stress Plots
I saw a few posts (can't locate them anymore) where people struggled to plot principal stress directions and magnitude (i.e. line vectors) in ParaView. This is a regression, because it worked back in 2019.
I managed to fix the problem (irrespective of material), as shown below:
My biggest challenge is how to generate a pull request to get this into FC0.20
I managed to fix the problem (irrespective of material), as shown below:
My biggest challenge is how to generate a pull request to get this into FC0.20
Re: Principal Stress Plots
I do have an account, but merging etc. before the pull request would cost me all day, so herewith the two affected files:
Re: Principal Stress Plots
Re: Principal Stress Plots
thanks!
before I can merge, I need more info:
- why did you delete 3 stresses in FemVTKTools.cpp: https://github.com/FreeCAD/FreeCAD/pull ... 128f50d576
(Why do they harm and are not necessary?)
- in resulttools.py you changed
calculate_principal_stress_reinforced
to
calculate_principal_stress
(https://github.com/FreeCAD/FreeCAD/pull ... s/935b0c1a)
Are you sure this has no negative side effect on simulations for reinforced concrete etc?
Re: Principal Stress Plots
1) I moved the three stresses up to the vector section, because that’s what they are and that’s what you need to display them in ParaView.
2) I wrote the concrete routines and am convinced that this works. Actually, the old concrete routine is now used for all materials and I just renamed it by dropping “concrete” in the name.
Thanks. Harry
2) I wrote the concrete routines and am convinced that this works. Actually, the old concrete routine is now used for all materials and I just renamed it by dropping “concrete” in the name.
Thanks. Harry
Re: Principal Stress Plots
That's not right then. See my above file and below:
Code: Select all
std::map<std::string, std::string> _getFreeCADMechResultVectorProperties() {
// see src/Mod/Fem/femobjects/_FemResultMechanical
// App::PropertyVectorList will be a list of vectors in vtk
std::map<std::string, std::string> resFCVecProp;
resFCVecProp["DisplacementVectors"] = "Displacement";
// the following three are filled only if there is a reinforced mat object
// https://forum.freecadweb.org/viewtopic.php?f=18&t=33106&start=70#p296317
// https://forum.freecadweb.org/viewtopic.php?f=18&t=33106&p=416006#p412800
resFCVecProp["PS1Vector"] = "Major Principal Stress";
resFCVecProp["PS2Vector"] = "Intermediate Principal Stress";
resFCVecProp["PS3Vector"] = "Minor Principal Stress";