The figure data is divided between two folders corresponding to the two numerical experiments presented in the manuscript. Visualisations were produced in Paraview, and the line graphs were generated using MATLAB.

4.1 Plane strain compression (biaxial) test

A .vtu is provided for each of the four levels of axial strain. Figures 3 and 6a were produced by viewing the 'internal hardening' field, while Figure 5a depicts the 'y' component of the 'rotation' vector converted into degrees. Figure 5b was produced by generating a 2D Delaunay triangulation on the final .vtu, on which two spatial gradients were calculated: the gradient of the internal hardening, and the gradient of the rotation angle. The former gradient was normalised (defining the vector n), and the plotted curvature k followed from taking its dot product with the latter gradient.

The band centrelines were manually sampled (along an isocontour of n), yielding the coordinates provided in the four .csv files (x,y,z coordinates in metres). These were passed into the plotSpline.m MATLAB script to apply a Savitsky-Golay filter and then generate the centreline orientation plot given as Figure 4. A window size of 5 was used in all cases.

An additional .vtk for the classical result contains the data used to generate Figure 6b. Here the equivalent plastic strain is the scalar field named 'alpha'.

4.2 Plane strain column collapse

The force-displacement curve (Figure 8) was generated from the data given in the .csv (displacement in metres and specific weight in kilonewtons per cubic metre).

Five .vtk files are provided for each of the snapshots. The accumulated plastic multiplier shown in Figure 9 is the scalar field 'alpha', and the curvature shown in Figure 10 follows from the same steps followed for previous example (swapping internal hardening for 'alpha', and rotation for 'phi_y').