Abstract The evolution of the shape of limit surfaces during plastic excursions is a key component of the inelastic response of structural members. However, stress-resultant-based computational plasticity formulations rarely account for any change in shape of the loading or ultimate strength surfaces. This research documents the evolution of the shape, size, and position of the loading and ultimate strength surfaces in stress-resultant space for wide-flange steel shapes subjected to cyclic proportional and nonproportional loading having a combination of axial force and either strong- or weak-axis flexure. The evolution is based on numerical integration of the cross-section stresses assuming a cyclic uniaxial stress–strain curve with a trilinear backbone that accounts for the Bauschinger effect. Surface evolution is investigated both with and without residual stresses. The results indicate that the evolution in shape of the loading surface, while relatively path-independent, is substantial, particularly at high levels of plastic excursion, where the loading surface essentially vanishes. The evolution in shape of the ultimate strength surface, in contrast, is mild. The research also documents the relative magnitudes of the components of cyclic plastic deformation. Recommendations are provided for the estimation of the direction of cyclic plastic flow through the calculation of gradients to the limit surfaces within the assumption of using an associated flow rule in stress-resultant space.