I am writing a FEA solver (using MATLAB) for computing stresses and stresses in a 2D rectangular plate with an elliptic hole (rotated to an angle). I computed the global stiffness matrix using Gauss quadrature (3 gauss points) and have also calculated stress and strain in each element.
How do I calculate the nodal stresses from the element stress values. I think the element stress values that I obtained are calculated at the centroid of each element. How do I extrapolate/interpolate these values to get nodal stress and strains?
Please advise with some resources or tutorials or solved examples.