Why are the stress, strain, moment and curvature discontinuous across element boundaries?
This is true for slffea-1.0, but for slffea-1.1, nodal averaging is used. Even though I feel it is more honest to assign color based upon Gauss points for each individual element, the output data files were getting much too big, and I had to switch to nodal averaging. (As an example, I had a brick output file which was 10 MB, and I reduced it to 2.6 MB by using nodal averaging) It should be noted that two elements which do not have nodal averaging are the beam and truss elements because stress, strain, moment and curvature are local to the coordinate system of each individual element. This is also true for the shell, but I explain why nodal averaging is justified in the README for this element.