You can now obtain axial and shear force for 2D quadrilateral elements in the latest version of SAGE-CRISP. These forces are obtained from the stress integrals of these 2D elements. The accuracy of these forces is dependent on the accuracy of the axial and shear stresses in the 2D elements, which are in turn, dependant on element size. The same restrictions apply as with bending moments obtained from these elements.
The following is a typical plot for axial force for a tunnel lining:
The following shows a typical shear force plot for a retaining wall: