Non-local strain regularisation techniques have been developed to mitigate mesh-dependency issues that develop in conventional finite element analyses involving strain softening. However, implementation of the non-local method generally requires high additional computational resources. In view of this, this paper introduces a more efficient procedure to implement the non-local method in Abaqus Coupled Eulerian-Lagrangian (CEL) analyses. The non-local method is examined for two boundary value problems: pipe penetration and vane shear tests in strain softening clays under undrained conditions. Results show that high computational efficiency is achieved by taking advantage of the stationary Eulerian element and the built-in mapping algorithm in CEL, even in the simulation of large-scale boundary value problems. The non-local method is effective in the pipe penetration problem where the strain localisation freely evolves within the soil domain, whereas this method fails to regularise strain localisation in the vane shear problem as the geometry governs the location of the strain localisation.