Hi everyone

I have written a MATLAB code to calculate the plane strain ans stress in a planar FEM problem.

My element is quadrilateral.

After resulting the displacement vector, I use that set of data to compute the strain and stress vectors using the following equation:

epsilon=[strain-displacement matrix]*{displacement vector}

I conduct this equation for every single Gauss point, therefore I end up 4 sets of vectors for every single element.

My question is that as I compared my results with the software, I realized that my results showing the strain and stress magnitudes on the Gauss points.

How should I map my results to the node to get the nodal strain and stress values?

I would appreciate your valuable comments.

Thanks.