Improved algorithms are proposed for a gradient plasticity theory in which the Laplacian of an invariant plastic strain measure enters the yield function. Particular attention is given to the type of finite elements that can be used within the format of gradient-dependent plasticity. Assuming a weak satisfaction of the yield function, mixed finite elements are developed, in which the invariant plastic strain measure and the displacements are discretized.
Two families of finite elements are developed: one in which the invariant plastic strain measure is interpolated using C1-continuous polynomials, and one in which penalty-enhanced C0-continuous interpolants are used. The performance of both families of finite elements is assessed numerically in one-dimensional and two-dimensional boundary value problems.
The regularizing effect of the used gradient enhancement in computations of elasto-plastic solids is demonstrated, both for mesh refinement and for the directional bias of the grid lines.