The coupled Eulerian–Lagrangian (CEL) method implemented in Abaqus is an established tool for modelling large deformations in numerical geomechanics. As shown in previous work, it can be extended to a hydro‐mechanically coupled scheme by exploiting the similarity of the heat balance equation to the mass balance equation of fluids. However, the distinction between effective and total contact stresses has not been possible in this approach, which hinders its application to problems with interface friction as frictional stresses are calculated based on total normal contact stress. An approach to circumvent this problem is presented in this work. Its relevance is demonstrated by the simulation of vibratory pile driving model tests in water‐saturated sand. The implementation by means of user‐subroutines is available and works for any total‐stress analysis.