This paper details a new algorithm for performing computational experiments of subglacial granular deformation. The numerical approach allows detailed studies of internal sediment and pore-water dynamics under shear. Feedbacks between sediment grains and pore water can cause rate-dependent strengthening, which additionally contributes to the plastic shear strength of the granular material. Hardening can stabilise patches of the subglacial beds with implications for landform development.