Mesh bias and shear band inclination in standard and non-standard continua

A severe, spurious dependence of numerical simulations on the mesh size and orientation can be observed in elasto-plastic models with a non-associated flow rule. This is due to the loss of ellipticity and may also cause a divergence in the incremental-iterative solution procedure. This paper first analyses the dependence of the shear band inclination in a biaxial test on the mesh size as well as on the mesh orientation. Next, a Cosserat continuum model, which has been employed successfully for strain-softening plasticity, is proposed to prevent loss of ellipticity. Now, numerical solutions result for shear band formation which are independent of the size and the orientation of the discretisation.