I am dealing with 2D shell modelling: I have some difficulties checking results in correspondence of a small radius bending, since the expanded 3D visualization “overlaps” some elements in the area I am interested into (and where, obviously, maximum PE values are located). From images you can understand what I am meaning:
you will get the results on non-expanded mesh, but with no bending stresses (averaging removes them).
This shell looks rather thick. Perhaps you could use regular solid elements at least for the filet (mixed meshes are possible thanks to tied contact and tie constraints: https://www.youtube.com/watch?v=Kc9Ln8zn3PE). It can be done even in a separate analysis via submodeling. Or you could define composite shell section to have more layers of elements there (then they can be also hidden individually in the results).
CalculiX always reports averaged nodal stresses in field output (there’s no option to obtain integration point stresses this way). You could request history output and get the stresses from all the integration points of the selected elements in the form of a large table (the same form is obtained with the “From History Output by Equation” tool you’ve mentioned). Then there’s this script that could help here (or the workaround in ParaView mentioned in the original post of that thread): Stresses at the integration points in ParaView - #2 by mkraska
