A virtual surface concept for nested yield surface plasticity

A return mapping algorithm based on the concept of virtual surfaces is presented for stress integration in the numerical solution of boundary value problems using an elastoplastic multiple yield surface model. Both kinematic hardening and isotropic softening yield surfaces are addressed. With the implementation of virtual surfaces, the active yield surface is located before a stress update is implemented. The return mapping algorithm consists of only two steps. The first step consists of performing the intermediate update of the trial stress corresponding to the translation of all yield surfaces inside the active yield surface. Little calculation is involved in this step. The second step consists of performing the conventional explicit one‐step radial return mapping of the updated trial stress on the corresponding active yield surface. The algorithm efficiency resides in the virtual surface based return mapping which preserves the elastic predictor representation of the updated trial stress, and limits iterations to at most two steps.