NURBS plasticity: Yield surface evolution and implicit stress integration for isotropic hardening

This paper extends the non-uniform rational basis spline (NURBS) plasticity framework of Coombs et al. (2016) to include isotropic hardening of the yield surfaces. The approach allows any smooth isotropic yield envelope to be represented by a NURBS surface. The key extension provided by this paper is that the yield surface can expand or contract through the movement of control points linked to the level of inelastic straining experienced by the material. The model is integrated using a fully implicit backward Euler algorithm that constrains the return path to the yield surface and allows the derivation of the algorithmic consistent tangent to ensure optimum convergence of the global equilibrium equations. This provides a powerful framework for elasto-plastic constitutive models where, unlike the majority of models presented in the literature, the underlying numerical algorithm (and implemented code) remains unchanged for different yield surfaces. The performance of the algorithm is demonstrated, and validated, using both material point and boundary values simulations including plane stress, plane strain and three-dimensional examples for different yield criteria.