Using NASTRAN, for a plate structure, the shear force, axial force, and moment are given by
SF = integral of XY_membrane_force
AFy = integral of Y_membrane_force
M = integral of (neutral axis-z)*Y_membrane_force
The neutral axis location can be estimated from the location of maximum xy shear
I used FEMAP to output the centroid of each element, and the corresponding element membrane forces/shears. I used MATLAB to locate the elements at a station, and intergrated using trapazoid command to get the desired forces/moments. Plotted the results verses station to be similar to a beam diagram.
I think the results can also be obtained from integrating the corresponding stresses multiplied by plate thickness.